Interpretation of Zc(4025)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(4025)$$\end{document} as the hidden charm tetraquark states via QCD Sum Rules

By using QCD Sum Rules, we found that the charged hidden charm tetraquark [cu][c¯d¯]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[c u][\bar{c} \bar{d}]$$\end{document} states with JP=1-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ J^P = 1^-$$\end{document} and 2+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ 2^+$$\end{document}, which are possible quantum numbers of the newly observed charmonium-like resonance Zc(4025)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(4025)$$\end{document}, have masses of m1-c=(4.54±0.20)GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{1^-}^c = (4.54 \pm 0.20) \, \text {GeV}$$\end{document} and m2+c=(4.04±0.19)GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{2^+}^c = (4.04 \pm 0.19) \, \text {GeV}$$\end{document}. The contributions up to dimension eight in the operator product expansion were taken into account in the calculation. The tetraquark mass of JP=2+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{P} \,{=}\, 2^{+}$$\end{document} state was consistent with the experimental data of Zc(4025)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(4025)$$\end{document}, suggesting the Zc(4025)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(4025)$$\end{document} state to possess the quantum number of JP=2+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P \,{=}\, 2^+$$\end{document}. Extending to the b-quark sector, the corresponding tetraquark masses m1-b=(10.97±0.25)GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{1^-}^b = (10.97 \pm 0.25) \, \text {GeV}$$\end{document} and m2+b=(10.35±0.25)GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{2^+}^b = (10.35 \pm 0.25) \, \text {GeV}$$\end{document} were obtained, which values are testable in future B-factories.


Introduction
A charged charmonium-like state Z c (4025) has been reported by BESIII in the process e + e − → (D * D * ) ± π ∓ [1]. Its mass and width are, respectively, (4,026.3 ± 2.6 ± 3.7) MeV and (24.8 ± 5.7 ± 7.7) MeV. The mass and decay modes imply that it is a charged hidden charm state, which is similar to Z c (3900) [2][3][4]. In addition, the BESIII Collaboration recently announced another structure called Z c (3885) [5] in the invariant mass spectrum of (DD * ) ± . The existing analyses [6,7] favor this state with the quantum number of J P = 1 + . Were the Z c (3885) of the same origin as the Z c (3900), the quantum numbers of the Z c (3900) would be 1 + . However, the quantum number of Z c (4025) so far is not well determined.
This paper utilize standard techniques of the QCD Sum Rules [8][9][10][11][12] to investigate the masses of charged hidden charm tetraquark states with two quantum numbers, i.e. a e-mail: qiaocf@ucas.ac.cn b e-mail: tangl@ucas.ac.cn J P = 1 − and 2 + . The hidden charm tetraquark state for Z c (4025) is investigated through examination of experimental data. Utilizing the QCD Sum Rules, the hidden charm tetraquark states with various quantum numbers have been investigated in Refs. [13][14][15][16][17], yielding significant conclusions. For J P = 1 − , an unstable mass sum rule was obtained [15], where the interpolating current was consistent with the 1 − current. However, stable results were extracted in Refs. [18][19][20]. This paper reanalyzes this case by adding several new ingredients [13] and invoking moderate criteria [14], to determine the available threshold parameter √ s 0 and the Borel window M 2 B . It should be noted that, very recently, the BESIII Collaboration has observed a charged charmonium-like resonance in the processes e + e − → (h c π ± )π ∓ , named Z c (4020) [21]. So far it is still too early to tell whether the Z c (4025) and the Z c (4020) are the same origin or not [22], though many theoretical investigations have already done [7,[23][24][25][26][27][28]. Among them, Braaten et al. [7] interpreted the exotic states as Born-Oppenheimer tetraquarks which are [cc][qq] (color octetoctet) states; Guo et al. [23] suggested the Z c (4025) as a D * D * virtual state; Khemchandani et al. [24] discussed the possibilities of the Z c (4025) being a 1 + or 2 + D * D * bound state in the framework of QCD Sum Rules; and Aceti et al. [25] also argued that the Z c (4025) is a 2 + state, but with D * D * as a bound state.
In Sect. 2, various essential formulas are presented. Numerical analysis and mass extraction are shown in Sect. 3, with conclusions given in Sect. 4.

Formalism
The QCD Sum Rules begin with the two-point correlation functions: The interpolating currents of 1 − and 2 + hidden charm tetraquark states are, respectively, constructed as where a, b, c, . . . are color indices, and C represents the charge conjugation matrix. For j 1 − μ (x), the correlation function has the following Lorentz covariance form: where the subscripts 1 and 0, respectively, denote the quantum numbers of the spin 1 and 0 mesons. The two-point function of the current j 2 + μν (x) has the following Lorentz form [29]: Here 2 (q 2 ) is the part of μν, αβ (q) which exclusively projects onto the 2 + state, and T μν, αβ is the unique Lorentz tensor of the fourth rank constructed from g μν and q μ : which satisfies the following desired properties: where g t μν (q) = (g μν − q μ q ν /q 2 ). On the phenomenological side, after separating out the ground state contribution from the pole term of the i (q 2 ), where i = 1 or 2, the correlation function is expressed as a dispersion integral over a physical regime, where m c i , λ c i and ρ h i (s), respectively, represent the mass, decay constant, and spectral density of the tetraquark state. Here s 0 is the threshold of higher excited states and continuum states, and λ c i is defined in Refs. [13,29]. It is worth to note that the tetraquark state defined here couples to the four-quark interpolating field, which has a special structure differing from the D * D * configuration. Hence the D * D * threshold effect is neglected here.
On the OPE side of the i (q 2 ), the correlation function is expressed as a dispersion relation: /π , and is expressed as where the "· · · " stands for other higher dimension condensates omitted in our work.
To evaluate the spectral density of the OPE side, the "full" propagators S q i j (x) and S Q i j ( p) of a light quark (q = u, d or s) and a heavy quark (Q = c or b) are, respectively, written with the vacuum condensates clearly displayed [10], thus: where the Lorentz indices κ and λ correspond to the indices of an outer gluon field from another propagator, and G represents the outer gluon field [30].
Using the techniques of Refs. [13,14] the spectral density ρ OPE i (s) was calculated up to dimension eight at the leading order in α s . The 1 − tetraquark state spectral density of the OPE side are given by where M B is the Borel parameter introduced by the Borel Matching the OPE side and the phenomenological side of the correlation function (q 2 ), i.e. using quark-hadron duality, and performing the Borel transformation, the sum rule for the mass of the hidden charm tetraquark state is determined to be with where O 6 represents g 3 s G 3 or qq 2 for i = 1 or i = 2.

Numerical analysis
In the numerical calculation, the values of the condensates and the quark masses are used as [14,31,32] In the QCD Sum Rules, to select the appropriate threshold s 0 and the Borel parameter M 2 B , there are two criteria [8][9][10]12]. As the convergence of the OPE must be retained in order to determine their convergence, it is essential to compare the relative contributions of each term to the total contributions on the OPE side.
The second criterion to constrain the M 2 B is that the pole contribution (PC) must be larger than the continuum contribution. Thus, for various values of the M 2 B , it is necessary to analyze the relative pole contribution, defined as the pole contribution divided by the total contribution, i.e. pole plus continuum. To safely eliminate the contributions of the higher excited and continuum states, the PC is generally greater than 50 % [12,14], which is slightly different from the constraint used in [15].
To find a proper value for √ s 0 , we carry out a similar analysis as in Refs. [12,33]. Since the continuum threshold is connected to the mass of the studied state by the relation √ s 0 ∼ m c i +0.5 GeV, various √ s 0 satisfying this constraint are taken into account. Among these values, one needs then to find the proper one, the one which has an optimal window for Borel parameter M 2 B . That is, within this window, the physical quantity, here the tetraquark mass m c i , is independent of the Borel parameter M 2 B as much as possible. Through the above procedure one obtains the central value of √ s 0 . However, in practice, in the QCD Sum Rules calculation, it is normally acceptable to vary the √ s 0 by 0.1 GeV [33], which gives the lower and upper bounds and hence the uncertainties of √ s 0 .

1 − Hidden Charm Tetraquark State
The OPE convergence of the 1 − hidden charm tetraquark state is shown in Fig. 1, which reflects a strong OPE convergence for M 2 B ≥ 2.1 GeV 2 , making it possible to determine the lower limit constraint of the M 2 B .  where M 2 B was 3.4 GeV 2 , and the errors stemmed from the uncertainties of the charm quark mass, the condensates and the threshold parameter √ s 0 .

2 + Hidden Charm Tetraquark State
For the 2 + sector, the OPE convergence is shown in Fig. 4, which reflects a strong OPE convergence for M 2 B ≥ 2.3 GeV 2 , enabling the determination of the lower limit constraint of the M 2 B . The result of the PC is shown in Fig. 5  The dependence of m c 2 + is drawn on the parameter τ in Fig. 6, where τ = 1/M 2 B , and the continuum thresholds √ s 0 are, respectively, taken as 3.9, 4.2, 4.5, 4.8, and 5.1 GeV, from down to up. In Fig. 6, the optimal mass curve of the 2 + hidden charm tetraquark state is shown with where both aforementioned criteria are satisfied. The mass of the 2 + hidden charm tetraquark was determined to be where M 2 B was 3.0 GeV 2 , and the errors stemmed from the uncertainties of the charm quark mass, the condensates and the threshold parameter √ s 0 .

Conclusions
In this paper, we estimated the masses of the hidden charm tetraquark states with J P = 1 − and 2 + , which are pos-sible quantum numbers possessed by the charmonium-like resonance Z c (4025). In our calculations contributions up to dimension eight in the OPE were taken into account. Noticeably, with the 1 − current, the mass was found to be m c 1 − = (4.54 ± 0.20) GeV, thus we deduced that such a tetraquark structure was not the candidate for Z c (4025). As in the discussions in Ref. [15], it may correspond to the charged partner of the charmonium-like state Y (4360) or Y (4660), within the uncertainties. However, in the case of 2 + , we found that m c 2 + = (4.04 ± 0.19) GeV, which is consistent within the errors with the experimental data of the Z + c (4025) resonance. As was mentioned in the introduction, the existing analyses favor Z c (3900) having quantum number of J P = 1 + , and these calculations cannot discriminate Z c (3900) with the Z c (3885). Possibly they may have the same origin. In this work, we calculate the masses of J P = 2 + tetraquark states in the framework of QCD Sum Rules. Our result suggests that the mass of the hidden charm tetraquark state is a bit more than 4 GeV, which in a certain degree agrees with the recent observations of Z c (4025) or Z c (4020) by BESIII Collaboration. For the b-quark sector, by virtue of a similar numerical analysis, the masses of the tetraquark state Addendum: during the finalization of this work, two reports about Z c (4025) appeared, wherein the molecular picture [34] and the initial-single-pion-emission mechanism [35] were employed.