Radiative and meson decays of $Y(4230)$ in flavor SU(3)

The charmonium-like exotic states $Y(4230)$ and the less known $Y(4320)$, produced in $e^+e^-$ collisions, are sources of positive parity exotic hadrons in association with photons or pseudoscalar mesons. We analyze the radiative and pion decay channels in the compact tetraquark scheme, with a method that proves to work equally well in the most studied $D^*\to \gamma/\pi+D$ decays. The decay of the vector $Y$ into a pion and a $Z_c$ state requires a flip of charge conjugation and isospin that is described appropriately in the formalism used. Rates however are found to depend on the fifth power of pion momentum which would make the final states $\pi Z_c(4020)$ strongly suppressed with respect to $\pi Z_c(3900)$. The agreement with BES III data would be improved considering the $\pi Z_c(4020)$ events to be fed by the tail of the $Y(4320)$ resonance under the $Y(4230)$. These results should renovate the interest in further clarifying the emerging experimental picture in this mass region.


I. INTRODUCTION
The study of final states in e + e − high energy annihilation, with the pioneering contributions by BaBar, Belle and BES collaborations, has opened the way to the new spectroscopy of exotic hadrons.
The so-called Y states, unexpected charmonium-like states created by the initial lepton pair, are efficient sources of positive parity exotic hadrons produced in association with one photon, pion or K meson.
Our results for Z c (3900) and X(3872) are in quantitative agreement with earlier studies of D * decays. The agreement is, of course, welcome but not unexpected and it supports the picture of compact tetraquarks bound by QCD forces.
We find a strong dependence of decay rates from the pion momentum, Γ ∝ q 5 . As a consequence, the decay Y (4230) → πZ c (4020) is strongly suppressed with respect to the decay into πZ c (3900), which does not seem to be supported by the cross sections reported by BES III. One possible explanation could be that the Z c (4020) events come from the second peak, Y (4320). A clarification of the distribution of πD * D * events in the region as well as information on the decay modes of Y (4320) would be very useful.
Production of exotic states in e + e − annihilation goes essentially via Y resonances. It is reasonable to assume that the open strangeness state Z cs (3985) seen in [2] e + e − → K + Z − cs (3985 If the hypothesis is correct, our analysis of K meson transitions, shows that strange members of the two nonets associated to X(3872) and Z c (3900) should both appear in the final states of (5), i.e. the D sD * 0 + c.c. spectrum should include as well the Z cs (4003) recently observed by LHCb in B + decay [15] (the classification of the newly discovered Z cs resonances is considered in [16]). This is a crucial feature that can be tested in higher luminosity experiments.
A J P C = 1 −− resonance, Y (4620), was first observed by BaBar and confirmed by Belle in e + e − annihilation with Initial State Radiation (ISR) [17,18]. BES III has later studied the 4620 structure with higher resolution and shown that it is resolved in two lines, now indicated as Y (4230) and Y (4320) (see [1]).
Y states as P -wave tetraquarks have been described in [11,13]. One expects four states Y 1 , . . . , Y 4 , the two lightest ones with spin composition Valence quark composition is [cq][cq], diquark and antidiquark spin are indicated in parenthesis, L is the orbital angular momentum. It was noted in [13] that the mass difference of Y 1,2 arises from two contrasting contributions: the hyperfine interaction, which pushes Y 1 down, and the spin-orbit interaction, which pushes Y 2 down. We had chosen M 1 < M 2 on the basis of a preliminary indication that the γ + X(3872) decay was associated to Y (4320), since this decay may arise from the Y 2 structure in Eq. (8) and not from Y 1 , Eq. (7).
Later information [19] indicates that the source of the γ + X(3872) decay is instead Y (4230). Consequently, we are led to change the assignment and propose M 2 < M 1 , that is (Γ e is the width to an electron pair). The total width of Y (4230) is estimated in [1] Γ(Y (4230)) = (56.0 ± 3.6 ± 6.9) MeV (11) Data from BES III indicate that Y (4230) is isoscalar [26,27]. Thus, denoting by Y u,d the Y 2 states with uū and dd valence quarks, we take We consider the decays where X, Z and Z ′ are the S-wave tetraquarks The decay (13) as a dipole transition L = 1 → L = 0, ∆S = 0 has been considered in [28]. Here we rederive the result as an introduction to pionic transitions.
We work in the non-relativistic approximation and describe the states with wave functions in spin and coordinate space. In the rest frame of Y 2 N Y is a normalization constant, the spin wave function from (8) and (16) is We indicate with a bar the charge-conjugate quark fields, x, y are diquark and antidiquark coordinates, ξ = x − y the relative coordinate and r = |ξ| the relative radius. The plus sign in S a(+) reminds of the charge conjugation, as defined on the rhs of (20). Considering the decay into X, we take and normalize spin w.f. according to Thus with y(r) = rR(r), and R(r) the radial wave function. Radiative decay. We work in the radiation gauge, A 0 = 0 and ∇ · A = 0. The photon couples to u and to other quarks with the basic lagrangian The elementary transition amplitudes are The minus sign in Mū arises from charge conjugation. In view of large mass denominators, we neglect radiation from the charm quarks. The right-hand sides of these equations contain products of operators acting on the spin and space wave functions of the initial tetraquark multiplied by variables of the electromagnetic field. As usual, we identify Further, we set and the hamiltonian acting on tetraquark wave functions is with E and B the electric and magnetic fields. The first term corresponds to the well known electric dipole transition that changes by one unit the orbital angular momentum, leaving the spin wave function unchanged [29]. One obtains with For isoscalar Y (4230) we use (12). Summing incoherently over the final states X u and X d , see [30], we get π 0 emission. We assume that quarks couple to pions via the isovector, axial vector current 1 : We follow [14] for the definition of the coupling g and The lagrangian contains the time derivative of the pion field. Applying the Legendre transformation, the interaction hamiltonian is The elementary quark transition is determined by The first term corresponds to ∆L = 0, operative in D * → Dπ [14], the second to ∆L = 1, for Y and D 1 pionic decay, dots indicate terms with ∆L > 1. Using charge conjugation symmetry, restricting to the ∆L = 1 term and specializing to the π 0 case, we obtain the hamitonian We note the results of applying the spin operators to the components of the spin wave function we obtain Note that, going from Y to Z or Z ′ , the minus sign between σ u and σū changes the charge conjugation sign of the spin w.f.. Similarly, the minus sign between the u and d term in (38) changes the S u and S d combination from I = 0 (in Y ) to I = 1 (in Z and Z ′ ).
In conclusion, we find and q is the decay momentum, ω Z = M Y − M Z , we have chosen to normalize the radius with ω X , for comparison with Eq. (32) and m u = 308 MeV, from the constituent quark model spectrum of mesons (see e.g. [3,10]).

IV. CHARGE CONJUGATION IN Y AND OTHER TETRAQUARK NONETS
A charge conjugation quantum number can be given to each self conjugate SU(3) f multiplet according to where C denotes the operator of charge conjugation, T the matrix representing the multiplet in SU(3) space andT the transpose matrix. η T is the sign taken by neutral members, but it can be attributed to all members of the multiplet.
In the exact SU(3) f limit, η is conserved in strong and electromagnetic decays. η = −1 is given to the electromagnetic current J µ and to Y µ while η K,π = +1. We extend Y 2 to a full nonet that we write as (omitting the overall normalization for brevity) D α a = (cσ 2 σ α q a );Db = (cσ 2qb ), etc. (46) F (ξ) is the wave function in the relative coordinate, even under ξ → −ξ.
The decay D * + → π 0 D + . From the hamiltonian (36) and Eq. (37), the relevant term in the hamiltonian is with p 0 the decay momentum. Also Γ(D * + → π + D 0 ) = g 2 (p + ) 3 6πf 2 π (decay momentum p + ). We reproduce the results of [14]. We assume D * + decay to be dominated by πD final states and use the D * + total width [32] to estimate the value of g, obtaining D 0 1 → π 0 D * 0 transition. D 1 (2420) is a well identified P -wave, positive parity charmed meson with total spin and angular momentum S = J = 1. We can use its decay into D * π to calibrate the ∆L = 1 hamiltonian (38). In analogy with (19), we write the D 1 wave function as: where the subscript D indicates that the QCD couplings of thecu system are used. The decay is induced by the ξ dependent part of the hamiltonian, restricted to the u term. Proceeding as before, we find and q 1 is the decay momentum, ω 1 = M D1 − M D * and we have assumed that the πD * modes saturate the total width. The transition radius in Eq. (59) is computed in the next Section, see Tab. II. Using the experimental width [32] we find: the error is estimated from the D 0 1 and D ± 1 width errors and variations in the estimated radius.

VI. TRANSITION RADIUS
The transition radius for a diquarkonium was estimated in [28], from the radial wave functions of a diquarkantidiquark system in a confining, QCD potential. We solve numerically the two body, radial Schrödinger equation [33] with potential and diquark mass Couplings are taken from lattice calculation of charmonium spectrum [34] α s = 0.3 k = 0.15 GeV 2 Alternatively, Ref. [28] uses the parameters of the Cornell potential [35] or a pure confinement case: For the D 1 → D * transition, we use the same potentials and M c = 1.7, m u = 0.308 GeV.  Results are reported in Tab. II
The Z c (4020) puzzle. The axial ∆L = 1 transition amplitude has a strong dependence from the pion momentum which reflects in a steep dependence of the rate: Γ ∝ q 5 , see Eq. (44). The pion momentum of Y (4230) → πZ c (4020) implies a suppression factor ∼ 30 with respect to Y (4230) → πZ c (3900), which does not seem to be supported by the cross sections in Tab. I.
Would it be possible that the Z c (4020) events come from the second peak of the structure, Y (4320)? A clarification of the source of D * D * events in the region and of the decay modes of Y (4320) would be very useful.