Y ( 4260 ) → γ + X ( 3872 ) in the diquarkonium picture

The observed Y ( 4260 ) → γ + X ( 3872 ) decay is a natural consequence of the diquark–antidiquark description of Y and X resonances. In this note we attempt an estimate of the transition rate through a non-relativistic calculation of the electric dipole term of a diquarkonium bound state. Combin-ing with BESIII data, upper bounds to B ( Y → J /(cid:3) + π + π) and to (cid:5)( Y → μ + μ − ) are obtained. We expect to confront these results with forthcoming data from electron–positron and hadron colliders.


Introduction
Exotic, hidden charm, mesons known as X, Y, Z resonances have been interpreted in [1,2] as tetraquarks, namely states made by two diquark pairs [cq][cq ] with q, q light quarks. Each pair is in color 3 or3 configuration, spin s,s = 1, 0 and relative orbital momentum L = 0, 1. The scheme has met with some degree of success at explaining the rich phenomenology which has emerged from electron-positron and hadron collider experiments. More information is expected in the future data from LHCb, BES III, and Belle II.
The long-standing conviction, based on consideration of large-N QCD, that tetraquark states could only materialize in the form of hadronic resonances too broad to be experimentally resolved, has been recently proven incorrect in [3,4]. Tetraquarks in large-N QCD have been further studied in [5][6][7]. The recent discovery of two pentaquark states of opposite parity [8] has reinforced the case for a new spectroscopic series of hadrons, in which diquarks (antidiquarks) replace antiquarks (quarks) in the classical scheme [9,10].
Recently, a new paradigm for the spin-spin interactions in hidden-charm tetraquarks has been proposed, which assumes the dominance of spin-spin couplings inside the diquark or the antidiquark [11]. This simple ansatz reproduces the mass ordering of the three, well identified, spin 1 + states, a e-mail: luciano.maiani@cern.ch X (3872), Z (3900) and Z (4020) and the pattern of their observed decays. In addition in Ref. [11] the diquark spin assignments of L = 1 states is discussed, pointing out that Y (4260) has the same spin distribution as X (3872) namely States are in the basis |s,s; L where s (s) is the diquark (antidiquark) spin and L the relative orbital angular momentum.
A similar scheme has been extended to exotic, hidden beauty mesons [12,13], and shown to give a consistent picture of the decays of ϒ(10890) into ϒ(nS)π + π − or h b (n P)π + π − , which occur via the intermediate Z b , Z b states [14].
The suppression of spin-spin interactions between a quark and an antiquark in different diquarks, underlined in [11], suggests that the overlap of the two constituents is very small, as if diquark and antidiquark were well separated entities inside the hadron. In the present paper we pursue this idea to the extreme consequences by considering the approximation where a diquark and an antidiquark can be described as pointlike. X, Y, Z would be, in this case, bound particleantiparticle systems, that we call diquarkonia for brief. We shall see that this extremely simplified picture leads to a reasonable approximation to the mass spectrum of S and P wave tetraquarks.
The diquarkonium picture has been introduced by Ali et al. [15] to study the production and decay of the Y (10890) considered as the b-tetraquark analog to the Y (4008). The annihilation of a diquarkonium made by s,s = 0 has been treated in [15] as the annihilation of a pair of spinless, pointlike particles. The extension to Y (4260) → μ + μ − deserves further consideration, given that the diquark and the antidiquark in the Y have not the same spin and the coupling to the photon is not simply determined by the charges.
We study the diquarkonia mass spectrum in the nonrelativistic approximation, using the Cornell potential previously applied to charmonia [16][17][18][19] and then, equipped with the corresponding wave functions, we compute the predicted rate of the ED1 allowed transition with L = 1, s = 0, which arises naturally from (1). Using the masses of the identified X, Y, Z states, we find parameters of the potential rather similar to the Cornell parameters and confirm the identification of the Z (4430) as the first radial excitation of Z (3900).
We compute the rates of the radiative transition for isospin I = 0, 1 of X (3872) and Y (4260). Assuming X (3872) to be an isospin singlet, we find and we compare this result to the available experimental information [20]. The rate of the radiative decay (2) has been computed in Ref. [21] in the molecular scheme, describing Y (4260) and X (3872) as D D 1 and D D bound states, respectively. The resulting rate turns out to be considerably smaller than the values indicated in (3) or (4).

Diquark masses
For S-wave states diquarkonia one writes the rest frame Hamiltonian where s (s) denotes the quark (antiquark) spin and M cq is the effective diquark mass. In the |s,s J basis, S-wave tetraquarks with J P = 1 + are described [11] by M cq can be estimated from the X (3872) and Z (4020) masses, subtracting the spin-spin contributions As a first approximation, we shall use M cq as input mass in the Schrödinger equation that gives the diquarkonia wave functions and masses.
In the case of charmonium, the input charm quark mass in the Schrödinger equation is obtained from the leptonic width (J/ → e + e − ); see [16][17][18]. In our case, the leptonic width of Y (4260) is not available yet and we shall be content to use the value (8) as input. We have verified that the various quantities are little sensitive (only to a few percent) to variations of the input diquark mass around this value.

Bound state masses
The simplest description of diquarkonia is in terms of a nonrelativistic potential, V (r ). For this first exploration we take the Cornell potential [16][17][18][19] with one chromo-Coulombic and one confining term For charmonia, one finds [19] A = 0.47, ν = 0.19 GeV 2 (charmonium spectrum). (10) For diquarkonia, we leave the parameters as free variables to be determined by comparison of diquarkonia eigenvalues 1S, 2S and 2P, to the mass differences of the J = 1 states, X (3872) or Z (3900), Y (4260) and Z (4430), subtracted of spin dependent terms. The subtraction is straightforward for the S-wave states, but for P-waves it requires the determination of not well known spin-orbit couplings [11], which introduces a non-negligible uncertainty.
Let us assume, as in [1,11], that we can write where in the r.h.s. we have introduced the eigenvalues of the Schrödinger equation, M 0 (1S), etc. Explicitly, spin interaction terms are obtained from a parametrization of the constituent quark Hamiltonian, which generalizes Eq. (5) to include orbital angular momentum excitation [11] 1 Signs are chosen so that, for B c , a, κ qc positive, energy increases for increasing L 2 and S 2 . As remarked in [11], this Hamiltonian is not the most general one as it does not include tensor terms which are known to be important in charmonium. The Hamiltonian describes well the J = 1 states but it could not be reliable for states with higher J .
Obvious manipulations lead to and we read M 00 is the analog of M 00 for the first radial excitation. Therefore and We use the mass values summarized in Table 1 [22] and take the value a = 73 MeV from the fit to the masses of Y states in [11] 2 to which we attribute a theoretical error estimated to be not less than 50 %. We find We solve numerically the Schrödinger equation [24] using the diquark mass in (8).
Results for the mass differences are reported in Fig. 1, in the plane of the eigenvalue differences 2S − 1S and 2P − 1S. The result for the Cornell potential with charmonium parameters is given by the round dot, whereas the squared box with errors corresponds to the eigenvalue differences estimated in (18). Lines indicate the results computed with fixed A while varying ν. Approximate agreement with the mass formula point is obtained for The difference 2S − 1S is well reproduced for both sets of parameters, Eqs. (10) and (19), reinforcing the case for Z (4430) to be the first radial excitation of Z (3900) [2,25]. The difference between the parameters in (10) and (19)   The squared box with errors corresponds to the eigenvalue differences estimated in (18). Lines indicate the results computed with fixed A and varying ν. The round dot represents the result for the Cornell potential with charmonium parameters given in Eq. (10) the diquark is not as pointlike as the c quark and, therefore, less sensitive to the short distance effects embodied by the Coulomb term.

The ED1 transition
We consider the process as the ED1 transition from a P-wave to a S-wave tetraquark with the same spin structure. Diquarks are taken as pointlike objects of electric charge Q The Hamiltonian (radiation gauge) is where A is the vector potential, x the coordinate and v the relative velocity of the particles in the center of mass system, with the diquark reduced mass and M 2q given by (8). In the dipole approximation where we set A(x) ≈ A(0), the matrix element for the decay is where and q are the polarization vector and momentum of the photon, ω = E f − E i its energy and m and k label the spin states of X and Y , respectively. The total rate is obtained by (25) = e 2 d 3 q (2π) 3 2ω where we used with n i = q i /ω.

Diquarkonium wave functions and transition radius
Consider first diquarkonia with a given flavor composition, In the non-relativistic approximation, state vectors corresponding to Y (P-wave) or X (S-wave) are written as where d and d c (or d m and d m c ) are the destruction operators of diquark and antidiquark with spin S = 0 (S = 1) and R(r ) the radial wave functions. We have made explicit the color index a = 1, 2, 3. The normalization factors are obtained from (non-relativistic) identities of the form to wit where (27) has been used and the number of colors is N = 3. The transition radius is then computed between normalized states to be and we have introduced the reduced radial wave functions of the 1S and 2P wave functions y(r ) = r R(r ) computed numerically [24]. Finally, we consider the general isospin structure of Y (4260) and X (3872), defining and obtain

Diquarkonium rate
With (26) and (36), we obtain Note that 0 ≤ Q 2 eff ≤ (4/3) 2 , with zero attained when Y = Y u and X = X d or vice versa, and the maximum when Y and X have only u-flavor.
As indicated by data, we take X (3872) close to a pure I = 0 state. For the two sets of parameters of the potential, Eqs. (10) and (19), we summarize in Table 2 (i) the numerical Table 2 Transition radius and corresponding decay widths for Y → γ X . In parentheses the branching ratio, assuming Y (4260) = 120 MeV [22] Charm. potential, Eq. (10 values of the transition radius and (ii) the rate for Y (4260) with I = 0, 1.
With the indicated numerical value of the radius, we are at the border of the dipole approximation, since ω r 0.8, not so much smaller than one. The situation, however, is not so different from the radiative transition χ c2 → J/ γ , which has ω r = 0.86, with estimated ∼10 % corrections; see [26].
From the value of (Y → X γ ) we can also estimate (Y → e − e + ). We use the well-known formula for the peak cross section with the experimental determination [27] σ (e − e + → Y (4260) → X γ ) = 0.33 pb B(X → π + π − J/ )

Conclusions
We estimated the transition rates (Y (4260) → γ + X (3872)) under both the assumptions that Y is an isospin singlet or a triplet bound state confined by a Cornell-like potential, a diquarkonium. We observe that the mass formula of the constituent quark model gives results for the mass differences between the radial excitations being closer to the results computed through an inter-diquark potential linearly rising with the distance and no chromo-Coulombic term. The results obtained, together with an upper bound estimate of the Y electronic width, can be confronted with future data, from electron-positron and hadron colliders.