Y(4260) -->gamma + X(3872) in the diquarkonium picture

The observed Y(4260)\to \gamma + 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, \Gamma_{\rm rad}, by a non-relativistic calculation of the electric dipole term of a diquarkonium bound state. We compute \Gamma_{\rm rad} for generic composition values of the isospin of X and Y. Specializing to I=0 for X(3872), we find \Gamma_{\rm rad}= 496~keV for Y(4260) with I=0 and \Gamma_{\rm rad}= 179~keV for I=1. Combining with BESIII data, we derive upper bounds to B(Y\to J/\Psi+\pi+\pi) and to \Gamma(Y\to \mu^+ \mu^-). 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 proton-(anti)proton 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]. Tetraquarks in large-N QCD have been further studied in [4]. The recent discovery of two pentaquark states of opposite parity [5] has reinforced the case for a new spectroscopic series of hadrons, in which diquarks (antidiquarks) replace antiquarks (quarks) in the classical scheme [6].
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 [7]. This simple ansatz reproduces the mass ordering of the three, well identified, spin 1 + states, X(3872), Z(3900) and Z(4020) and the pattern of their observed decays. In addition in Ref. [7] 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 [8], and shown to give a consistent picture of the decays of Υ(10890) into Υ(nS)π + π − or h b (nP )π + π − , which occur via the intermediate Z b , Z ′ b states [9].
The suppression of spin-spin interactions between a quark and an antiquark in different diquarks, underlined in [7], 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 particle-antiparticle 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 A. Ali et al. [10] 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 with s,s = 0 has been treated in [10] 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 non-relativistic approximation, using the Cor-nell potential previously applied to charmonia [11,12] 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 = 179 keV (I : 1 → 0) (4) and compare this result to the available experimental information [14].
The rate of the radiative decay (2) has been computed in Ref. [13] in the molecular scheme, describing Y (4260) and X(3872) as DD 1 and DD ⋆ 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 [7] 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 [11]. 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 few percents) to variations of the input diquark mass around this value.

Bound state masses
The simplest description of diquarkonia is in term of a non-relativistic potential, V (r). For this first exploration we take the Cornell potential [11,12] with one Chromo-Coulombic and one confining term For charmonia, one finds [12] A = 0.47, ν = 0.19 GeV 2 (charmonium spectrum) 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 [7], which introduces a non negligible uncertainty.
Let us assume, as in [1,7], that we can write Obvious manipulations lead to and we read (M ′ 00 is the analog of M 00 for the first radial excitation) so that We use the mass values summarized in Tab. 1 [15] and take the value a = 73 MeV from the fit to the masses of Y states in [7] 2 to which we attribute a theoretical error estimated to be not less than 50%. We find M 2S − M 1S = 0.60 ± 0.03 GeV M 2P − M 1S = 0.23 ± 0.07 GeV (18) We solve numerically the Schrödinger equation [16] 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 round dot represents the result for the Cornell potential with charmonium parameters given in Eq. (10) and the squared box with errors corresponds to the eigenvalue differences estimated in (18). Lines indicate the results computed with fixed A and varying ν.
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,4]. The difference between the parameters in (10) and (19) may be due to the inaccuracy of the mass formula or to the fact that the diquark is not as pointlike as the c quark, therefore less sensitive to the short distance effects embodied by the Coulomb term.

The ED1 transition
We consider the process Y (4260) → γ + X(3872) (20) 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 centre 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) with n i = q i /ω.

Diquarkonium wave-functions and transition radius
Consider first diquarkonia with a given flavor composition, e.g.
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 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) = rR(r) computed numerically [16].
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 viceversa 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 Tab. 2 (i) the numerical 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 [19].
The result found in Ref. [14] can be stated as which, assuming [15] Using our result, we predict 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 [18] σ(e − e + → Y (4260) → X γ) = 0.33 pb B(X → π + π − J/Ψ) and the input values in Table 1  for Y isospin equal to 0 or 1, respectively.

Conclusions
We estimated the transition rates Γ(Y (4260) → γ + X(3872)) under the assumption that Y is a diquarkonium confined by a Cornell-like potential, either isospin singlet or triplet. Our results reinforce the case for Z(4430) to be the first radial excitation of Z(3900) [2,4]. Mass differences between states with different orbital excitations computed with a linearly rising potential and no Chromo-Coulombic term approximately agree with the mass formula derived from the constituent quark model in [7]. The results obtained, together with upper bound estimates of B(Y → J/Ψ ππ) and of the Y electronic width, can be confronted with future data, from electron-positron and hadron colliders.