Analysis of the mass and width of Y(4274) as axialvector molecule-like state

In this article, we assign Y(4274) to the color octet–octet type axialvector molecule-like state with JPC=1++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{PC}=1^{++}$$\end{document} tentatively, and construct the color octet–octet type axialvector current to study its mass and width with the QCD sum rules in details. The predicted mass favors assigning the Y(4274) to a color octet–octet type molecule-like state, but the predicted width strongly disfavors assigning the Y(4274) to the color octet–octet type molecule-like state. The Y(4274) may be the conventional charmonium state χc1(3P)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _{c1}(\mathrm 3P)$$\end{document}, and it is important to observe the decay Y(4274)→J/ψω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y(4274)\rightarrow J/\psi \omega $$\end{document} to diagnose the nature of Y(4274).

Recently, the LHCb collaboration performed the first full amplitude analysis of the decays B + → J/ψφ K + and confirmed the two old particles Y (4140) and Y (4274) in the J/ψφ mass spectrum with statistical significances 8.4σ and 6.0σ , respectively, and determined the quantum numbers to be J PC = 1 ++ with statistical significances 5.7σ and 5.8σ , a e-mail: zgwang@aliyun.com respectively [7,8] (1) The LHCb collaboration determined the quantum numbers of the Y (4274) to be J PC = 1 ++ , which rules out the 0 −+ molecule assignment, which is consistent with our previous work [5]. There have been several possible assignments, such as the color sextet-sextet type cscs tetraquark state [9,10], the conventional orbitally excited state χ c1 (3P) [11,12], the color triplet-triplet type 1 √ 6 (uū + dd − 2ss)cc tetraquark state [13], etc.
In Ref. [16], we study the axialvector hidden charm and hidden bottom tetraquark states in details with the QCD sum rules and obtain the mass M cucd,J P =1 + = (4.44 ± 0.19) GeV for the diquark-antidiquark type tetraquark state. In Ref. [17], Chen and Zhu study the vector and axialvector charmonium-like tetraquark states with the QCD sum rules in a systematic way and obtain the mass M cucd,J P =1 + = (4.16 ± 0.10) GeV for the diquark-antidiquark type tetraquark state. In Ref. [16], we choose the input parameters m c (1 GeV), qq (1 GeV), qg s σ Gq (1 GeV), while in Ref. [17], Chen and Zhu choose the input parameters m c (m c ), qq (1 GeV), qg s σ Gq (1 GeV). The different predictions for the Cγ 5 ⊗ γ μ C type axialvector tetraquark state in Ref. [16] and Ref. [17] originate from the different choice of the c-quark mass. If we take different choices of the heavy quark masses as a source of uncertainties, the predicted mass is about M cucd,J P =1 + = (4.06-4.63) GeV.
In Ref. [18], we distinguish the charge conjugations of the interpolating currents, study the diquark-antidiquark type axialvector tetraquark states in a systematic way with the QCD sum rules by taking into account the energy scale dependence of the QCD spectral densities for the first time, and obtain the predictions M X (3872) = 3.87 +0.09 −0.09 GeV and M Z c (3900) =3.91 +0.11 −0.09 GeV. In Ref. [19], Albuquerque et al. take into account the next-to-leading order and next-to-next-to-leading order factorizable radiative corrections to the perturbative terms, and obtain the predication M cucd,J P =1 + = (3.888 ± 0.130) GeV, which also depends on special choice of the energy scale μ, in other words, the M S mass m c (μ). The non-factorizable radiative corrections are still needed to make precise predictions. In leading order approximation, M X (3872) = 3.87 +0.09 −0.09 GeV and M Z c (3900) = 3.91 +0. 11 −0.09 GeV [18], after taking into account the next-to-leading order and next-to-next-to-leading order factorizable radiative corrections to the perturbative terms, M cucd,J P =1 + = (3.888 ± 0.130) GeV [19], the predicted masses only change slightly. On the other hand, including the next-to-leading order and next-to-next-to-leading order factorizable radiative corrections to the perturbative terms leads to the value of the pole residue λ cucd,J P =1 + undergoes the replacement λ cucd,J P =1 + → 1.09 λ cucd,J P =1 + . According to Refs. [18,19], the masses of the ground state diquarkantidiquark type axialvector tetraquark states cucd are about 3.9 GeV.
In Ref. [20], Chen et al. assign Z c (4200) to be the ground state axialvector tetraquark state cucd, calculate its decay width with the QCD sum rules, and obtain the value Z c (4200) = 435 ± 180 MeV. In Ref. [9], Chen et al. assign X (4140) to be the ground state axialvector tetraquark state cscs. If Z c (4200) and X (4140) are the color triplet-triplet Cγ 5 ⊗ γ μ C type axialvector tetraquark states, it is more natural in the case that the X (4140) has larger mass than Z c (4200).
In Ref. [21], we assign Z c (4200) to the color octet-octet type axialvector molecule-like stateūλ a ccλ a d, where the λ a are the Gell-Mann matrices, and we construct the color octet-octet type axialvector current to study its mass (width) with the QCD sum rules by calculating the vacuum condensates up to dimension 10 (5)  −132 MeV from the Belle collaboration [14,15], and favor assigning the Z c (4200) to be the color octet-octet type molecule-like state with J PC = 1 +− . Moreover, we study the energy scale dependence of the QCD spectral density of the molecule-like state in details and suggest an empirical energy scale formula to determine the ideal energy scale, in other words, to determine the ideal c-quark mass.
Also in Ref. [21], we discuss the possible assignments of Z c (3900), Z c (4200) and Z (4430) as the ground state color triplet-triplet diquark-antidiquark type tetraquark states with J PC = 1 +− in details. The QCD sum rules support assigning Z c (3900) and Z (4430) to the ground state and the first radial excited state of the diquark-antidiquark type axialvector tetraquark states with J PC = 1 +− , respectively [18,22]. If we assign Z c (4200) and Y (4274) to molecule-like states with J PC = 1 +− and 1 ++ , respectively, the mass difference M Y (4274) − M Z c (4200) ≈ 77 MeV. This is reasonable, as the SU (3) breaking effects are very small for the four-quark systems [16,23,24]. In this article, we assign Y (4274) to the color octet-octet type molecule-like state tentatively, study its mass and decay width with the QCD sum rules in details, where the meson-like states D a s and D a * s have the same quark constituents as the mesons D s and D * s respectively, but they are in the color octet representation, a corresponds to the Gell-Mann matrices.
The article is arranged as follows: we derive the QCD sum rules for the mass and width of the color octet-octet type axialvector molecule-like state Y (4274) in Sect. 2 and in Sect. 3, respectively; Sect. 4 is reserved for our conclusion.

The mass of the color octet-octet type axialvector molecule-like state
In the following, we write down the two-point correlation function μν ( p) in the QCD sum rules, where the λ a are the Gell-Mann matrices in the color space. We construct the color octet-octet type current J μ (x) to study the molecule-like state Y (4274). One can consult Refs. [21,[25][26][27] for more literatures on the color octet-octet type currents. Under charge conjugation transform C, the current J μ (x) has the property, At the phenomenological side, we insert a complete set of intermediate hadronic states with the same quantum numbers as the current operator J μ (x) into the correlation function μν ( p) to obtain the hadronic representation [28][29][30], and isolate the ground state contribution, In the following, we briefly outline the operator product expansion for the correlation function μν ( p). We contract the quark fields s and c in the correlation function μν ( p) with the Wick theorem, and we obtain the result where and t n = λ n 2 , D α = ∂ α − ig s G n α t n [30], then compute the integrals both in the coordinate space and in the momentum space, and we obtain the correlation function μν ( p), therefore the QCD spectral density through dispersion relation. For technical details, one can consult Ref. [18]. Now we take the quark-hadron duality below the continuum threshold s 0 and perform Borel transform with respect to the variable P 2 = −p 2 to obtain the QCD sum rule: where The explicit expressions of the QCD spectral densities ρ 0 (s), ρ 3 (s), ρ 4 (s), ρ 5 (s), ρ 6 (s), ρ 7 (s), ρ 8 (s) and ρ 10 (s) are given in the Appendix. Even in the leading order approximation, the strong coupling constant g 2 s (μ) = 4πα s (μ) appears according to the equation of motion D α G βα a = g s q=u,d,sq γ β t a q, see the terms g 2 s ss 2 in the spectral density ρ 6 (s). So we have to consider the energy scale dependence of the QCD sum rules, the preferred c-quark mass is the M S mass m c (μ).
The on-shell quark propagator has no infrared divergences in perturbation theory, and this provides a perturbative definition of the quark mass [31,32]. But the pole mass cannot be used to arbitrarily high accuracy because of nonperturbative infrared effects in QCD. The pole mass m c and the M S mass m c (m c ) have the relation, The value m c (m c ) = 1.275 ± 0.025 GeV from the Particle Data Group corresponds to m c = 1.67 ± 0.07 GeV [33]. If we take the pole mass, then 2 m c > M J/ψ > M η c , at the phenomenological side of the QCD sum rules for J/ψ and η c , If we want to obtain nonzero values, we have to choose smaller pole mass, 2 m c < M η c < M J/ψ [34,35]. For an observable particle such as the electron, the physical mass appears as the pole mass, irrespective of the leading order, next-to-leading order, next-to-next-to-leading order, . . ., radiative corrections are concerned. In the leading order approximation, m c = m c (m c ), however, m c (m c ) originates from the radiative corrections and renormalization, which are beyond the leading order approximation. So in the leading order approximation, the definition of the pole mass m c is of arbitrary. Again, we can see that the preferred cquark mass is the M S mass m c (μ). Moreover, the full quark propagator has no pole because the quarks are confined. The pole mass corresponds to a non-confined particle, at the QCD side of the QCD sum rules for J/ψ and η c , the heavy quarks c andc are confined particles. We derive Eq. (11) with respect to τ = 1 T 2 , then eliminate the pole residue λ Y (4274) to obtain the QCD sum rule for the mass, Now we choose the input parameters at the QCD side of the QCD sum rules. We take the vacuum condensates to have the standard values qq = −(0.24 ± 0.01 GeV) 3 [33]. Moreover, we take into account the energy-scale dependence of the quark condensate, mixed quark condensate and M S masses from the renormalization group equation [33,37], 4 9 , sg s σ Gs (μ) = sg s σ Gs (Q) α s (Q) α s (μ) 2 27 , 12 25 , where t = log μ 2 2 , = 213 MeV, 296 MeV and 339 MeV for the flavors n f = 5, 4 and 3, respectively [33].
As the quark masses m c (μ), m s (μ), the quark condensate ss (μ), the mixed condensate sg s σ Gs (μ) all depend on the energy scale μ, the QCD spectral density ρ(s) depends on the energy scale μ, we have to determine the energy scales of the QCD sum rules for those molecule-like states in a consistent way.
The hidden charm (or bottom) four-quark systems qq QQ can be described by a double-well potential. In the four-quark system qq QQ, the heavy quark Q serves as a static well potential and combines with the light quark q to form a heavy diquark D i in color antitriplet, q + Q → D i [18,[22][23][24][38][39][40][41][42], or combines with the light antiquarkq to form a heavy meson in color singlet (meson-like state in color octet), q + Q →q Q (q λ a Q) [21,27,43]; the heavy antiquarkQ serves as another static well potential and combines with the light antiquarkq to form a heavy antidiquark D i in color triplet,q +Q → D i [18,[22][23][24][38][39][40][41][42], or combines with the light quark q to form a heavy meson in color singlet (meson-like state in color octet), q +Q →Qq (Qλ a q) [21,27,43], where i is color index, the λ a are the Gell-Mann matrices. Then the two heavy quarks Q andQ stabilize the four-quark systems qq QQ, just as in the case of the (μ − e + )(μ + e − ) molecule in QED [44]. The four-quark systems qq QQ are characterized by the effective heavy quark mass M Q and the virtuality V = [18,[21][22][23][24]27,[38][39][40][41][42][43]. It is natural to take the energy scale μ = V . The M Q is just an empirical parameter to determine the optimal energy scales of the QCD spectral densities, and has no relation to the pole massm Q or the M S mass m Q (μ).
In Refs. [18,[21][22][23][24]27,[38][39][40][41][42][43], we observe that there exist three universal values for the effective heavy quark masses M Q , which correspond to the compact tetraquark states, molecular states, molecule-like states, respectively. The empirical energy scale formula We evolve all the input parameters in the QCD spectral density to the special energy scale determined by the empirical formula, In Ref. [21], we obtain the effective mass M c = 1.98 GeV for the molecule-like states. Then we re-checked the numerical calculations and corrected a small error concerning the mixed condensate, the updated value is M c = 2.01 GeV. From the empirical energy scale formula, we can obtain the energy We take into account all uncertainties of the input parameters, and we obtain the values of the mass and pole residue, which are shown explicitly in Fig. 1, In Fig. 1, we plot the mass and pole residue of the Y (4274) with variation of the Borel parameter T 2 at a larger interval than the Borel window. From the figure, we can see that there appear platforms, criterion 3 is also satisfied. Now the three criteria are all satisfied, it is reliable to extract the ground state mass. The predicted mass M Y (4274) = (4.27 ± 0.09) GeV is consistent with the experimental value 4273.3 ± 8.3 +17.2 −3.6 MeV from the LHCb collaboration [7,8], which supports assigning Y (4274) to the color octet-octet typesλ a ccλ a s molecule-like state.
In Ref. [19], Albuquerque et al. study the hidden-charm and hidden-bottom molecular states and tetraquark states by taking into account the next-to-leading order and next-tonext-to-leading order radiative corrections to the perturbative terms from the factorizable Feynman diagrams (without including the non-factorizable Feynman diagrams). The numerical results indicate that the predicted masses are slightly modified, while the decay constants (which relate to the pole residues) are modified significantly, the largest modification amounts to multiplying the decay constants by a factor 1.8. So we expect that the predication M Y (4274) = (4.27±0.09) GeV survives approximately even if the next-toleading order radiative corrections to the perturbative terms are taken into account. Moreover, at the present time, even the next-to-leading order factorizable contributions are not available for the color octet-octet type molecule-like states, it is a challenging work to calculate both the next-to-leading order factorizable and non-factorizable Feynman diagrams.

The width of the color octet-octet type axialvector molecule-like state
We can study the strong decay Y (4274) → J/ψφ with the three-point correlation function αμν ( p, q), where the currents interpolate the mesons J/ψ and φ(1020) according to the current-hadron couplings, f J/ψ and f φ are the decay constants, the ξ α and ζ μ are polarization vectors of the mesons J/ψ and φ(1020), respectively. At the phenomenological side, we insert a complete set of intermediate hadronic states with the same quantum numbers as the current operators J J/ψ α (x), J φ μ (y), J ν (0) into the threepoint correlation function αμν ( p, q) and isolate the ground state contributions to obtain the result where p = p+q, G Y J/ψφ is the hadronic coupling constant, which is defined by In this article, we choose the tensor structure ε αμνλ p λ to study the coupling constant G Y J/ψφ . The two unknown functions ρ Y φ ( p 2 , t, p 2 ) and ρ Y J/ψ (t, q 2 , p 2 ) parameterize transitions between the ground states and the higher resonances or the continuum states, the net effects can be parameterized by C Y φ and C Y J/ψ , In calculations, we take the C Y φ and C Y J/ψ as free parameters, and vary them to eliminate the contaminations to obtain Borel platforms [45,46].
We carry out the operator product expansion up to the vacuum condensates of dimension 5 and neglect the gluon condensate, which plays a minor important role. We obtain the QCD spectral density through the dispersion relation, take the quark-hadron duality below the continuum thresholds, then set p 2 = p 2 and take double Borel transforms with respect to the variables P 2 = −p 2 and Q 2 = −q 2 respectively to obtain the QCD sum rule where s 0 Y and s 0 φ are the continuum threshold parameters for Y (4274) and φ(1020), respectively.
The hadronic parameters are taken as M φ = 1.019461 GeV, M J/ψ = 3.0969 GeV [33], f J/ψ = 0.418 GeV [47], If the radiative corrections to the perturbative term of the correlation function αμν ( p, q) also amount to multiplying a factor about 1.8, as the color octet-octet type current J μ (x) is also presented, the value of the quantity f φ f J/ψ λ Y G Y J/ψφ at the hadronic side in the QCD sum rules in Eq. (26) changes according to the rule, In this article, we take the values f J/ψ = 0.418 GeV [47] and f φ = 0.253 GeV, which include next-to-leading order radiative corrections. The factors 1.8 come from the radiative corrections to the two-point correction function and threepoint correlation function cancel out with each other, the net modification of the hadronic coupling constant G Y J/ψφ is estimated to be tiny, just like the hadronic coupling constants D * Dπ and B * Bπ , the net effects of the radiative corrections can be neglected [48]. Now it is easy to obtain the decay width, where . It is difficult to assign Y (4274) to the color octet-octet type moleculelike statesλ a ccλ a s. In Ref. [23], we assign Y (4140) to the diquark-antidiquark type tetraquark state cscs with J PC = 1 ++ , and we study the mass and pole residue with the QCD sum rules in detail. The predicted mass disfavors assigning Y (4140) to the J PC = 1 ++ diquark-antidiquark type tetraquark state cscs. Y (4140) and Y (4274) have the same quantum numbers except for the masses and widths, the QCD sum rules also disfavor assigning Y (4274) to the J PC = 1 ++ diquark-antidiquark type tetraquark state cscs.

Conclusion
In this article, we assign Y (4274) to the color octet-octet type axialvector molecule-like state with J PC = 1 ++ tentatively, and construct the color octet-octet type axialvector current to study its mass and width with the QCD sum rules in details. The predicted mass M Y (4274) = (4.27 ± 0.09) GeV is consistent with the experimental value 4273.3 ± 8.3 +17.2 −3.6 MeV from the LHCb collaboration, and favors assigning Y (4274) to the color octetoctet type molecule-like statesλ a ccλ a s. The predicted width (Y (4274) → J/ψφ) = 1.8 GeV is much larger than the experimental value 56 ± 11 +8 −11 MeV from the LHCb collaboration and disfavors assigning the Y (4274) to be the color octet-octet type molecule-like state strongly. The Y (4274) may be the conventional charmonium state χ c1 (3P), and the