Vector tetraquark state candidates: Y(4260 / 4220), Y(4360 / 4320), Y(4390) and Y(4660 / 4630)

In this article, we construct the C⊗γμC\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C \otimes \gamma _\mu C$$\end{document} and Cγ5⊗γ5γμC\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\gamma _5 \otimes \gamma _5\gamma _\mu C$$\end{document} type currents to interpolate the vector tetraquark states, then carry out the operator product expansion up to the vacuum condensates of dimension-10 in a consistent way, and obtain four QCD sum rules. In calculations, we use the formula μ=MY2-(2Mc)2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =\sqrt{M^2_{Y}-(2{\mathbb {M}}_c)^2}$$\end{document} to determine the optimal energy scales of the QCD spectral densities, moreover, we take the experimental values of the masses of the Y(4260 / 4220), Y(4360 / 4320), Y(4390) and Y(4660 / 4630) as input parameters and fit the pole residues to reproduce the correlation functions at the QCD side. The numerical results support assigning the Y(4660 / 4630) to be the C⊗γμC\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C \otimes \gamma _\mu C$$\end{document} type vector tetraquark state cc¯ss¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c\bar{c}s\bar{s}$$\end{document}, assigning the Y(4360 / 4320) to be Cγ5⊗γ5γμC\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\gamma _5 \otimes \gamma _5\gamma _\mu C$$\end{document} type vector tetraquark state cc¯qq¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c\bar{c}q\bar{q}$$\end{document}, and disfavor assigning the Y(4260 / 4220) and Y(4390) to be the pure vector tetraquark states.


Introduction
In 2005, the BaBar collaboration studied the initial-state radiation process e + e − → γ I S R π + π − J/ψ and observed the Y (4260) in the π + π − J/ψ invariant-mass spectrum [1]. Later the Y (4260) was confirmed by the Belle and CLEO collaborations [2,3], the Belle collaboration also observed an evidence for a very broad structure Y (4008) in the π + π − J/ψ mass spectrum. In 2014, the BES collaboration searched for the production of e + e − → ωχ cJ with J = 0, 1, 2, and observed a resonance in the ωχ c0 cross section, the measured mass and width of the resonance are 4230 ± 8 ± 6 MeV and 38 ± 12 ± 2 MeV, respectively [4].
Also in 2016, the BES collaboration precisely measured the cross section of the process e + e − → π + π − J/ψ at center-of-mass energies from 3.77 to 4.60 GeV and observed two resonant structures. The first resonance has a mass of 4222.0±3.1±1. 4 MeV and a width of 44.1±4.3±2.0 MeV, the second one has a mass of 4320.0 ± 10.4 ± 7.0 MeV and a width of 101.4 +25. 3 −19.7 ± 10.2 MeV [6]. The first resonance agrees with the Y (4260) while the second resonance agrees with the Y (4360) according to the uncertainties, the Y (4008) resonance previously observed by the Belle experiment is not confirmed [2].
In 2006, the BaBar collaboration observed an evidence for a broad structure at 4.32 GeV in the π + π − ψ mass spectrum in the process e + e − → π + π − ψ [8]. In 2007, the Belle collaboration studied the initial-state radiation process e + e − → γ I S R π + π − ψ , and observed two structures Y (4360) and Y (4660) in the π + π − ψ invariant-mass spectrum [9,10]. In 2008, the Belle collaboration studied the initial-state radiation process e + e − → γ I S R + c − c and observed a clear peak Y (4630) in the + c − c invariant-mass spectrum [11]. The Y (4360) and Y (4660/4630) were confirmed by the BaBar collaboration [12].
In Ref. [13], Maiani et al. assign the Y (4260) to be the first orbital excitation of a diquark-antidiquark state [cs] S=0 [cs] S=0 based on the effective Hamiltonian with the spin-spin and spin-orbit interactions [42], where the subscript S denotes the diquark spins. In the type-II diquarkantidiquark model [14], where the spin-spin interactions between the quarks and antiquarks are neglected, L. Maiani et al interpret the Y (4360) and Y (4660) as the first radial excitations of the Y (4008) and Y (4260) respectively, and interpret the Y (4008), Y (4260), Y (4290/4220) and Y (4630) as the four ground states with the angular momentum L = 1. One can consult Ref. [15] for detailed reviews of the effective Hamiltonian approach.
In Ref. [16], Ali et al. analyze the hidden-charm P-wave tetraquarks and the newly excited charmed c states in the diquark model using the effective Hamiltonian incorporating the dominant spin-spin, spin-orbit and tensor interactions, and observe that the preferred assignments of the ground state tetraquark states with L = 1 are the Y (4220), Y (4330), In the effective Hamiltonian, an approximately common rest frame for all the components is assumed, while in the dynamical diquark picture, the diquark-antidiquark pair forms promptly at the production point, and rapidly separates due to the kinematics of the production process, then they create a color flux tube or string between them [17].
Another possible assignment of the Y (4260) is the ccg hybrid state [25,26], the lattice calculations indicate that the vector charmonium hybrid has a mass about 4285 MeV, which is quite close to that of the Y (4260) [27]. Furthermore, lattice results strongly indicate that the quarks should form a spin singlet in the low-lying vector hybrid, which in rough agreement with the spin-flip suppression in the annihilation e + e − → Y (4260). An alternative approach is to describe the X , Y , Z mesons as bound states in the Born-Oppenheimer potentials (usually lattice-QCD computed gluon-induced potentials) for a heavy quark and a heavy antiquark [28]. For more literatures on the QQ potentials involving or not involving of the gluonic excitations, one can consult Refs. [43,44].
In Ref. [29], Li and Voloshin suggest that the Y (4260) and Y (4360) are a mixture, with mixing close to maximal, of two states of the hadro-charmonium, one containing a spin-triplet cc pair and the other containing a spin-singlet cc pair. While the Y (4660) can be assigned to be the ψ f 0 (980) hadrocharmonium or molecular state (the two scenarios overlap in this case) [31,32]. For more literatures on the hadrocharmonium states, one can consult Ref. [30].
In the scenario of molecular states, the Y (4260) and Z c (3900) are assigned to be theD D 1 (2420) + DD 1 (2420) andD D * + DD * molecular states respectively [33,34], which are compatible with the processes e + e − → J/ψπ + π − , h c π + π − measured by the BESIII and Belle collaborations. In Ref. [35], Wang et al. confront both the hadronic molecule and the hadro-charmonium interpretations of the Y (4260) with the available experimental data, and conclude that the data support the Y (4260) being dominantly aD D 1 (2420) + DD 1 (2420) hadronic molecule while they challenge the hadro-charmonium interpretation. The assignment of the Y (4260) as the tetraquark state or molecular state means different decay ratio Y (4260) → γ X (3872) [45,46]. Precisely measuring the decay ratio can shed light on the nature of the Y (4260), Z c (3900) and X (3872).
In Table 1, we present the existing predictions of the masses of the vector tetraquark states based on the QCD sum rules. In Refs. [20][21][22]61,62], the vacuum condensates are taken at the energy scale μ = 1 GeV while the M S mass m c (m c ) is taken at the energy scale μ = m c (m c ), the energy scales of the QCD spectral densities are not specified. In Refs. [23,24,60], we take the formula μ = M 2 X/Y /Z − (2M c ) 2 with the effective charmed quark mass M c to determine the energy scales of the QCD spectral densities of the hiddencharm tetraquark states, and evolve the vacuum condensates and the M S mass m c (m c ) to the optimal energy scales μ to extract the tetraquark masses; the hidden-bottom tetraquark states can be studied analogously [63,64].
In this article, we assume that the Y (4260/4220), Y (4360/4320), Y (4390) and Y (4660/4630) are vector tetraquark states, and restudy the C ⊗γ μ C and Cγ 5 ⊗γ 5 γ μ C type vector tetraquark states with the QCD sum rules in details by taking into account the vacuum condensates up to dimension 10 in a consistent way in the operator product expansion, and use the energy scale formula μ = 2 to determine the optimal energy scales of the QCD spectral densities. In the QCD sum rules for the tetraquark states, the terms associate with 1 T 2 , 1 T 4 , 1 T 6 in the QCD spectral densities manifest themselves at small values of the Borel parameter T 2 , we have to choose large values of the T 2 to warrant convergence of the operator product expansion and appearance of the Borel platforms. The higher dimensional vacuum condensates play an important role in determining the Borel windows therefore the ground state masses and pole residues, though they maybe play a less important role in the Borel windows. We should take them into account consistently.
The article is arranged as follows: we derive the QCD sum rules for the masses and pole residues of the vector tetraquark states in Sect. 2; in Sect. 3, we present the numerical results and discussions; Sect. 4 is reserved for our conclusion.

QCD sum rules for the vector tetraquark states
In the following, we write down the two-point correlation functions μν ( p) in the QCD sum rules, where , the i, j, k, m, n are color indexes, the C is the charge conjugation matrix. Under charge conjugation transform C, the currents J μ (x) have the properties, In the non-relativistic diquark-antidiquark model, one often introduces an explicit P-wave between the diquark and antidiquark in the ground state Cγ 5 ⊗γ 5 C type or Cγ 5 ⊗γ μ C type or Cγ μ ⊗ γ ν C type tetraquark state [13][14][15][16]. In this arti-cle, we choose the C ⊗ γ μ C type and Cγ 5 ⊗ γ 5 γ μ C type vector currents to interpolate the vector tetraquark states, the net effects of the relative P-waves are embodied in the underlined γ 5 in the Cγ 5 γ 5 ⊗ γ μ C type and Cγ 5 ⊗ γ 5 γ μ C type currents or in the underlined γ α in the Cγ α γ α ⊗ γ μ C type currents.
The tetraquark states are spatial extended objects, not point-like objects [17], in the QCD sum rules [18][19][20][21][22][23][24] and the effective Hamiltonian approach [13][14][15][16], the finite size effects are neglected, which leads to some uncertainties, while in the potential models, an explicit spatial extended potential between the diquark and antiquark is introduced [65,66]. Now we perform Fierz re-arrangement to the vector currents J 1 μ (x) and J 3 μ (x) both in the color and Dirac-spinor spaces, and obtain the following results, the Fierz re-arrangement of the J 2 μ (x) and J 4 μ (x) can be obtained analogously. The diquark-antidiquark type current with special quantum numbers couples potentially to a special tetraquark state, while the current can be re-arranged to a current as a special superposition of color singlet-singlet type currents, which couple potentially to the meson-meson pairs or molecular states. The diquark-antidiquark type tetraquark state can be taken as a special superposition of a series of meson-meson pairs, and embodies the net effects.
At the hadronic side, we can insert a complete set of intermediate hadronic states with the same quantum numbers as the current operators J μ (x) into the correlation functions μν ( p) to obtain the hadronic representation [67][68][69]. After isolating the ground state contributions of the vector tetraquark states which are supposed to be the Y (4260/4220), Y (4360/4320), Y (4390) and Y (4660/4630), we get the following results, where the pole residues λ Y are defined by In the following, we take the currents J 1 μ (x) and J 3 μ (x) as an example, and briefly outline the operator product expansion for the correlation functions μν ( p) in perturbative QCD. We contract the c and s quark fields in the correlation functions μν ( p) with Wick theorem, obtain the results: 3 (14) where the S i j (x) and C i j (x) are the full s and c quark propagators respectively, and t n = λ n 2 , the λ n is the Gell-Mann matrix [69], then compute the integrals both in the coordinate and momentum spaces, and obtain the correlation functions μν ( p) therefore the spectral densities at the level of quark-gluon degrees of freedom. In Eq. (15), we retain the terms s j σ μν s i originate from the Fierz re-arrangement of the s is j to absorb the gluons emitted from other quark lines to extract the mixed condensate sg s σ Gs [63,64].
Once analytical expressions of the QCD spectral densities are obtained, we can take the quark-hadron duality below the continuum thresholds s 0 and perform Borel transform with respect to the variable P 2 = −p 2 to obtain the following four QCD sum rules: where ρ(s) = ρ 1 (s), ρ 2 (s), ρ 3 (s) and ρ 4 (x), the subscripts i = 0, 3, 4, 5, 6, 7, 8, 10 denote the dimensions of the vacuum condensates, the explicit expressions of the QCD spectral densities ρ i (s) are presented in the Appendix.
In this article, we carry out the operator product expansion to the vacuum condensates up to dimension-10 and discard the perturbative corrections, and assume vacuum saturation for the higher dimensional vacuum condensates. The higher dimensional vacuum condensates are always factorized to lower dimensional vacuum condensates with vacuum saturation in the QCD sum rules, factorization works well in large N c limit. In reality, N c = 3, some (not much) ambiguities maybe come from the vacuum saturation assumption. The higher dimensional vacuum condensates have not been well studied yet.
We derive Eq. (17) with respect to τ = 1 T 2 , then eliminate the pole residues λ Y , and obtain the QCD sum rules for the masses of the vector tetraquark states,

Numerical results and discussions
We take the standard values of the vacuum condensates qq = −(0.24 ± 0.01 GeV) 3 [71], and set m u = m d = 0. Moreover, we take into account the energy-scale dependence of the input parameters on the QCD side, 12 25 , 12 25 , 12 25 , m s (μ) = m s (2GeV) α s (μ) α s (2GeV) where t = log μ 2  [71], and evolve all the input parameters to the optimal energy scales μ to extract the masses of the vector hidden-charm tetraquark states.
In the present QCD sum rules, we search for the ideal Borel parameters T 2 and continuum threshold parameters s 0 to obey the following four criteria: 1. Pole dominance at the phenomenological side; 2. Convergence of the operator product expansion; 3. Appearance of the Borel platforms; 4. Satisfying the energy scale formula, using try and error.
In Refs. [23,60,63,64], we study the energy scale dependence of the QCD sum rules in details and suggest an energy scale formula μ = M 2 X/Y /Z − (2M Q ) 2 with the effective heavy quark mass M Q to determine the energy scales of the QCD spectral densities of the hidden-charm and hiddenbottom tetraquark states, which also works well for the hidden-charm pentaquark states [72,73]. In this article, we take the updated value M c = 1.82 GeV [24].
The resulting Borel parameters or Borel windows T 2 , continuum threshold parameters s 0 , ideal energy scales of the QCD spectral densities, pole contributions of the ground state tetraquark states, and contributions of the vacuum condensates of dimension 7, 8 and 10 in the operator product expansion are shown explicitly in Table 2. From the table, we can see that the first two criteria of the QCD sum rules are satisfied, so we expect to make reasonable predictions.
We take into account all uncertainties of the input parameters, and obtain the values of the masses and pole residues of the vector tetraquark states, which are shown explicitly in Figs. 1 and 2 and Table 3. From Figs. 1 and 2, we can see that there appear platforms in the Borel windows, the criterion 3 is satisfied. From Table 3, we can see that the criterion 4 is also satisfied. Now the four criteria of the QCD sum rules are all satisfied, and we expect to make reliable predictions.
In Fig. 1, we also present the experimental values of the masses of the Y (4260/4220), Y (4360/4320), Y (4390) and Y (4660/4630) [5,71]. From the figure, we can see that the experimental values of the masses M Y (4660/4630) and M Y (4360/4320) can be well reproduced, the present predications support assigning the Y (4660/4630) to be the C ⊗γ μ C type tetraquark state ccss, and assigning the Y (4360/4320) to be the Cγ 5 ⊗ γ 5 γ μ C type tetraquark state ccqq. The mass M Y (4660/4630) lies just below the upper bound of the predicted mass of the C ⊗γ μ C type vector tetraquark state ccqq, which disfavors assigning the Y (4660/4630) to be the C ⊗γ μ C type vector tetraquark state ccqq, however, such an assignment is not excluded.
On the other hand, the Y (4660) can be assigned to be the ψ f 0 (980) hadro-charmonium or molecular state based on fitting the mass distribution of the process e + e − → ψ π + π − [31] or the calculations of the QCD sum rules [32]. More experimental data are still needed to assign the Y (4660) unambiguously.
In this article, we recalculate the QCD sides of the correlation functions for the C ⊗ γ μ C type currents by taking into account the neglected terms due to the approximations involving the higher dimensional vacuum condensates in Ref. [23] and correct a small error in numerical calculations, the present predictions 4.66 ± 0.09/4.59 ± 0.08 GeV are more robust than the values 4.70 +0.14 −0.10 /4.66 +0.17 −0.10 GeV obtained in Ref. [23]. In Fig. 3, we plot the mass and pole residue of the C ⊗ γ μ C type vector tetraquark state ccss with variation of the Borel parameter T 2 for truncations of the operator product expansion, D = 6, 7, 8 and 10. From the figure, we can  Table 3 The masses and pole residues of the vector tetraquark states, where the superscripts 1, 2, 3 and 4 denote the currents  The ground state Cγ 5 ⊗ γ 5 C type and Cγ μ ⊗ γ μ C type hidden-charm tetraquark states ccqq have the masses about 3.85 GeV from the QCD sum rules in which the vacuum condensates up to dimension 10 are taken into account in a consistent way [58][59][60], if an additional P-wave costs about 0.5 GeV, the ground state vector hidden-charm tetraquark states ccqq have the mass about 4.35 GeV, the present prediction 4.34±0.08 GeV is robust. In Ref. [21], Zhang and Huang introduce an explicit P-wave in the currents, and obtain the value 4.32 ± 0.20 GeV by taking into account the vacuum condensates up to dimension 6.
In Ref. [16], Ali et al. study the hidden-charm P-wave tetraquarks and the newly excited charmed c states with the effective Hamiltonian incorporating the dominant spinspin, spin-orbit and tensor interactions, and observe that the Y (4220), Y (4330), Y (4390), Y (4660) can be assigned to be the four ground states with L = 1 by fitting the coefficients in the effective Hamiltonian to the experimental masses. In the effective Hamiltonian approach, the lowest state is the Y (4220), while we cannot obtain such low mass based on the QCD sum rules [18][19][20][21][22][23][24].
The mass M Y (4390) lies just below the upper bound of the predicted mass of the Cγ 5 ⊗ γ 5 γ μ C type vector tetraquark state ccqq, which disfavors assigning the Y (4390) to be the Cγ 5 ⊗ γ 5 γ μ C type vector tetraquark state ccqq, on the other hand, the mass M Y (4390) lies just below the lower bound of the predicted mass of the Cγ 5 ⊗γ 5 γ μ C type vector tetraquark state ccss, which also disfavors assigning the Y (4390) to be the Cγ 5 ⊗ γ 5 γ μ C type vector tetraquark state ccss, however, such assignments are not completely excluded. If we take the energy scale μ = 3.4 GeV for the ccss tetraquark state or μ = 1.8 GeV for the ccqq tetraquark state, the experimental value of the M Y (4390) can be reproduced, however, such energy scales are not consistent with the QCD sum rules for other tetraquark states.
There are no candidates for the Y (4260/4220) in the present calculations. In Refs. [33,34], the Y (4260) and Z c (3900) are assigned to be theD D 1 (2420) + DD 1 (2420) andD D * + DD * molecular states respectively based on the heavy meson (non-relativistic) effective field theory. In Ref. [36], we study the vector molecular states DD 1 (2420) and D * D * 0 (2400) with the QCD sum rules by taking into account the vacuum condensates up to dimension-10 in the operator product expansion in a consistent way, and use the energy scale formula for the molecular states to determine the optimal energy scales of the QCD spectral densities [74,75], and obtain the predications M DD 1 (1 −− ) = 4.36 ± 0.08 GeV and M D * D * 0 (1 −− ) = 4.78±0.07 GeV. The QCD sum rules support assigning the Y (4390) (not the Y (4260/4220)) and Z c (3900) to be the DD 1 and DD * S-wave molecular states, respectively [36,74,75]. While the lattice QCD supports assigning the Y (4260) to be a hybrid state [27]. Furthermore, there have been observed evidences for the X (3872) in the lattice calculations, though its interpretation was not specified [76], there was no evidence for the Z c (3900) in the lattice calculations [77]. There are also other assignments of the Y (4390), for example, the D * (2010)D 1 (2420) molecular state [78].   [5,71], and take the pole residues λ Y as free parameters to fit the following two QCD sum rules, In Table 4, we present the central values of the fitted pole residues. In Fig. 4, we plot the correlation functions with the central values of the fitted pole residues compared to the operator product expansion. From the figure, we can see that the QCD sides of the correlation functions can be well reproduced. From In Refs. [16,79], the Y (4220) is assigned to be the lowest vector tetraquark state in the simple diquark-antidiquark model with the constituent ccqq. In the present work, we can see that the lowest tetraquark states couple potentially to the current J 4 μ (x), not to the currents J 1 μ (x), J 2 μ (x) and J 3 μ (x), now we suppose that the Y (4260/4220) is a pure vector tetraquark state and saturates the QCD sum rules, and study the energy scale dependence of the extracted mass M Y (4260) with the central values of the input parameters in Table 2. In If we take the energy scale formula μ = 2 with the effective mass M c = 1.82 GeV as a constraint, the QCD sum rules only support assigning the Y (4660/4630) and Y (4360/4320) to be the C ⊗ γ μ C type vector tetraquark ccss and Cγ 5 ⊗ γ 5 γ μ C type vector tetraquark state ccqq respectively, and disfavor assigning the Y (4660/4630) (or Y (4390)) to be the C ⊗ γ μ C (or Cγ 5 ⊗ γ 5 γ μ C) type vector tetraquark state ccqq.

Conclusion
In this article, we construct the C ⊗ γ μ C and Cγ 5 ⊗ γ 5 γ μ C type currents to interpolate the vector tetraquark states, then calculate the contributions of the vacuum condensates up to dimension-10 in the operator product expansion in a consistent way, and obtain four QCD sum rules. In calculations, we use the formula μ = M 2 X/Y /Z − (2M c ) 2 to determine the optimal energy scales of the QCD spectral densities, explore the energy scale dependence of the QCD sum rules in details, moreover, we take the experimental values of the masses of the Y (4260/4220), Y (4360/4320), Y (4390) and Y (4660/4630) as input parameters and fit the pole residues to reproduce the correlation functions at the QCD side. The numerical results support assigning the Y (4660/4630) to be the C ⊗ γ μ C type vector tetraquark state ccss, assigning the Y (4360/4320) to be Cγ 5 ⊗ γ 5 γ μ C type vector tetraquark state ccqq, and disfavor assigning the Y (4260/4220) and Y (4390) to be the pure vector tetraquark states.