Analysis of the mass and width of the X(4140) as axialvector tetraquark state

In this article, we construct both the [sc]T[s¯c¯]A+[sc]A[s¯c¯]T\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[sc]_T[{\bar{s}}{\bar{c}}]_A+[sc]_A[{\bar{s}}{\bar{c}}]_T$$\end{document} type and [sc]T[s¯c¯]V-[sc]V[s¯c¯]T\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[sc]_T[{\bar{s}}{\bar{c}}]_V-[sc]_V[{\bar{s}}{\bar{c}}]_T$$\end{document} type axialvector currents 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} to study the mass of the X(4140) with the QCD sum rules. The predicted masses support assigning the X(4140) to be the [sc]T[s¯c¯]V-[sc]V[s¯c¯]T\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[sc]_T[{\bar{s}}{\bar{c}}]_V-[sc]_V[{\bar{s}}{\bar{c}}]_T$$\end{document} type axialvector tetraquark state. Then we study the hadronic coupling constant gXJ/ψϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{XJ/\psi \phi }$$\end{document} with the QCD sum rules based on solid quark-hadron duality, and obtain the decay width Γ(X(4140)→J/ψϕ)=86.9±22.6MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma (X(4140)\rightarrow J/\psi \phi )=86.9\pm 22.6\,{\mathrm{MeV}}$$\end{document}, which is in excellent agreement with the experimental data 83±21-14+21MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$83\pm 21^{+21}_{-14} { \text{ MeV }}$$\end{document} from the LHCb collaboration.


Introduction
In 2009, the CDF collaboration observed the X (4140) for the first time in the J/ψφ mass spectrum in the exclusive B + → J/ψ φ K + decays in pp collisions with a statistical significance more than 3.8σ [1]. Then the X (4140) was confirmed by CDF, CMS, D0, LHCb collaborations [1][2][3][4][5][6][7]. The LHCb collaboration performed the first full amplitude analysis of the decays B + → J/ψφ K + and confirmed the two old particles X (4140) and X (4274) in the J/ψφ mass spectrum with statistical significances 8.4σ and 6.0σ , respectively, and determined the spin-parity-changeconjugation to be J PC = 1 ++ with statistical significances 5.7σ and 5.8σ , respectively [6,7]. In Table 1, we present the mass, width, J PC of the X (4140) from the different experiments. Although the width from the LHCb collaboration [6,7] differs from other measurements greatly, the masses from different experiments are consistent with each other. The D * sD * s threshold is 4224.4 MeV from the Particle Data Group [8], which leads to the possible molecule assignment for the X (4140). The X (4140) was observed in the final state J/ψφ, its J PC = 0 ++ , 1 ++ , 2 ++ for the S-wave couplings, and 0 −+ , 1 −+ , 2 −+ , 3 −+ for the P-wave couplings. The a e-mail: zgwang@aliyun.com most popular current to interpolate the D * s meson is J α (x) = s(x)γ α c(x), the most popular current to interpolate the D * sD * s molecular states is J αβ (x) =s(x)γ α c(x)c(x)γ β s(x). We can study the J PC = 0 ++ , 1 −+ , 2 ++ , 1 +− , 1 −− D * sD * s molecular states with the QCD sum rules by using the suitable projectors. The LHCb collaboration determined the quantum numbers of the X (4140) to be J PC = 1 ++ , which rules out the 0 ++ or 2 ++ D * sD * s molecule assignment, but does not rule out the existence of the 0 ++ or 2 ++ D * sD * s molecular states.
In the QCD sum rules for the hidden-charm (or hiddenbottom) tetraquark states and molecular states, the integrals are sensitive to the energy scales μ, where the ρ QC D (s, μ) are the QCD spectral densities, the T 2 are the Borel parameters, the s 0 are the continuum threshold parameters, the predicted masses depend heavily on the energy scales μ. In Refs. [28][29][30][31], we suggest an energy scale formula μ = 2 with the effective Q-quark mass M Q to determine the ideal energy scales of the QCD spectral densities. The formula enhances the pole contributions remarkably, we obtain the pole contributions as large as (40 − 60)% in Refs. [10,[21][22][23], otherwise, the pole contributions are about 40% [24] or 20% [26,27]. The energy scale formula also works well in the QCD sum rules for the hidden-charm pentaquark states [32][33][34].
From  2 and choose the updated value M c = 1.82 GeV [35], we can obtain the optimal energy scales μ = 1.4 GeV and 2.0 GeV for the QCD spectral densities in the QCD sum rules for the Z c (3900) and X (4140), respectively. In Ref. [23], we observe that the mass of the X (4140) can be reproduced at the energy scale μ = 1. T type axialvector currents to study the mass of the X (4140) as the axialvector tetraquark state with the QCD sum rules in details, then study the width of the X (4140) with the QCD sum rules based on the solid quark-hadron duality.
The article is arranged as follows: we derive the QCD sum rules for the mass and width of the X (4140) as axialvector tetraquark state in Sects. 2 and 3 respectively; Sect. 4 is reserved for our conclusion.

The mass of the X (4140) as the axialvector tetraquark state
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. Under charge conjugation (parity) transform C ( P), the currents J μ (x) have the properties, the four vectors T axialvector tetraquark states with J PC = 1 ++ , respectively. The tensor diquark operators have the properties, under parity transform, the tensor diquark operators couple potentially to both the J P = 1 + and 1 − diquark states. We should project out the 1 + or 1 − component by multiplying tensor diquark operators by the axialvector antidiquark operator ε imnsm (x)γ ν Cc T n (x) or vector antidiquark operator ε imnsm (x)γ 5 γ ν Cc T n (x).
At the phenomenological side, we 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 [36][37][38], and isolate the ground state contributions, where the pole residues λ X are defined by 0|J μ (0)|X ( p) = λ X ε μ , the ε μ are the polarization vectors of the axialvector tetraquark states X . Now we briefly outline the operator product expansion for the correlation functions μμ ( p). We contract the quark fields s and c in the correlation functions μμ ( p) with Wick theorem, and obtain the results, where and t n = λ n 2 [38], then compute the integrals both in the coordinate space and in the momentum space, and obtain the correlation functions μμ ( p) (i.e. 1 μμ ( p) and 2 μμ ( p)), therefore the QCD spectral densities through dispersion relation ρ(s) = lim ε→0 . For technical details, one can consult Ref. [39]. Now we 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 QCD sum rules: the subscripts i in the QCD spectral densities ρ i (s) denote the dimensions of the vacuum condensates, the lengthy expressions of the QCD spectral densities are given in Appendix. We derive Eq. (11) with respect to τ = 1 T 2 , then eliminate the pole residues λ X to obtain the QCD sum rules for the masses, At the QCD side, we take the vacuum condensates to be the standard values qq = −(0.24±0.01 GeV) 3 [8]. 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, ss (μ) = ss (Q) α s (Q) α s (μ) 12 33−2n f , where t = log μ 2 2 = 210 MeV, 292 MeV and 332 MeV for the flavors n f = 5, 4 and 3, respectively [8,41], and evolve all the input parameters to the typical energy scales μ satisfying the energy scale formula μ = 2 to extract the masses of the axialvector tetraquark states.
We search for the optimal Borel parameters T 2 and threshold parameters s 0 to satisfy 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, via try and error, and obtain the Borel windows T 2 , threshold parameters s 0 , optimal energy scales of the QCD spectral densities, and pole contributions of the ground states, see Table 3. Now we take a short digression to illustrate how to impose the four criteria to choose the Borel parameters T 2 and continuum threshold parameters s 0 . Firstly, we set M X = 3.9 GeV tentatively, and obtain the energy scale μ = 1.4 GeV according to the energy scale formula. Then we take the continuum threshold parameters to be √ s 0 = (3.9 + 0.5) GeV as the energy gap between the ground state and the first radial excited state is about (0.4 − 0.6) GeV, and obtain the predicted masses M X , pole contributions, and the contributions of the vacuum condensates of dimension 10. We observe that the predicted masses M X are much larger than 3.9 GeV and the pole contributions are much smaller than 50% in the regions where the Borel platforms appear, Table 3 The Borel windows, continuum threshold parameters, ideal energy scales, pole contributions, masses and pole residues for the axialvector tetraquark states furthermore, the contributions of the vacuum condensates of dimension 10 are not small enough. Then we choose the masses M X > 3.9 GeV, say M X = 4.0 GeV, 4.1 GeV, . . . and reiterate the same procedure until obtain the optimal Borel parameters T 2 and continuum threshold parameters s 0 satisfying the four criteria. From Table 3, we can see that the pole dominance criterion is well satisfied. In calculations, we observe that the contributions of the vacuum condensates of dimension 10 are 1% (about 1%) in the QCD sum rules for the current J 1 μ (x)(J 2 μ (x)), the operator product expansion is well convergent. We take into account all uncertainties of the input parameters, and obtain the values of the masses and pole residues of the axialvector tetraquark states, see Table 3  From Table 3, we can see the predicted mass M X = 4.14 ± 0.10 GeV for the [sc] T [sc] V − [sc] V [sc] T axialvector tetraquark state is in excellent agreement with the experimental data 4146.5 ± 4.5 +4.6 −2.8 MeV from the LHCb collaboration [6,7], which supports assigning the X (4140) to be the

The width of the X (4140) as the axialvector tetraquark state
We can study the two-body strong decay X (4140) → J/ψφ with the three-point correlation function αβμ ( p, q), where the currents interpolate the mesons J/ψ, φ(1020) and X (4140) respectively, the f J/ψ and f φ are the decay constants, the ξ α and ζ β are polarization vectors of the mesons J/ψ and φ(1020), respectively. In this section, we will use the notation m X in stead of M X for the special case J μ (0) = J 2 μ (0). At the phenomenological side, we insert a complete set of intermediate hadronic states with the same quantum numbers as the current operators J into the threepoint correlation function αβμ ( p, q) and isolate the ground state contributions to obtain the result, where p = p + q, the g X J/ψφ is the hadronic coupling constant defined by the four functions ρ X φ ( p 2 , p 2 , t), ρ X ψ ( p 2 , t, q 2 ), ρ X J/ψ (t , p 2 , q 2 ) and ρ X φ (t , p 2 , q 2 ) have complex dependence In this article, we choose the tensor structure ε αβμλ p λ to study the g X J/ψφ to avoid the contaminations from the relevant scalar and pseudoscalar mesons according to the non-vanishing coupling constants, where the f χ c0 , f f 0 and λ X 0 are the decay constants of the χ c0 (3414), f 0 (980) and X 0 (J P = 0 − ), respectively. We introduce the parameters C X φ , C X ψ , C X φ and C X J/ψ to parameterize the net effects, Then the correlation function ( p 2 , p 2 , q 2 ) on the phenomenological side can be written as Now we carry out the operator product expansion up to the vacuum condensates of dimension 5 and neglect the tiny contributions of the gluon condensate. The correlation function QC D ( p 2 , p 2 , q 2 ) can be written as through dispersion relation, where the ρ QC D ( p 2 , s, u) is the QCD spectral density, we introduce the subscript QC D to denote the QC D side. We rewrite the correlation function H ( p 2 , p 2 , q 2 ) on the hadron side as through dispersion relation, where the ρ H (s , s, u) is the hadronic spectral density, we introduce the subscript H to denote the hadron side. However, on the QCD side, the QCD spectral density ρ QC D (s , s, u) does not exist, because lim 3 →0 We math the hadron side of the correlation function with the QCD side of the correlation function, and carry out the integral over ds firstly to obtain the solid duality [42], Fig. 3 The hadronic coupling constant g X J/ψφ with variation of the Borel parameter T 2 M X = 4146.5 MeV [6,7], λ X = 4.30 × 10 −2 GeV 5 . At the QCD side, we take the energy scale of the QCD spectral density to be μ = 2 GeV, just like in the two-point QCD sum rules. Then we set the Borel parameters to be T 2 1 = T 2 2 = T 2 for simplicity. The unknown parameters are chosen as C X J/ψ + C X φ = −0.00261 GeV 7 to obtain platform in the Borel window T 2 = (3.6 − 4.6) GeV 2 .
In Fig. 3, we plot the hadronic coupling constant g X J/ψφ with variation of the Borel parameter T 2 . From the figure, we can see that there appears platform in the Borel window indeed. After taking into account the uncertainties of the input parameters, we obtain the hadronic coupling constant g X J/ψφ , g X J/ψφ = −(2.23 ± 0.29).
Now it is easy to obtain the decay width,