Analysis of the D D K system with QCD sum rules Zun-Yan Di 1,2, Zhi-Gang Wang 1 1 Department of Physics, North China Electric Power University, Baoding 071003, P. R. China 2 School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, P. R. China arxiv:1901.05196v1 hep-ph] 16 Jan 2019 Abstract In this article, we construct the color singlet-singlet-singlet interpolating current to study the D D K system with QCD sum rules. In calculations, we consider the contributions of the vacuum condensates up to dimension-16 and use the formula µ = M 2 X/Y/Z (2Mc)2 to determine the energy scale of the QCD spectral density. The numerical result M Z = 4.71 +0.19 0.11 GeV indicates that there exists a virtual state or a resonance state Z lying above the D D K threshold to saturate the QCD sum rules. PACS number: 12.39.Mk, 12.38.Lg Key words: Resonance state, QCD sum rules 1 Introduction Since the observation of the X(3872) by the Belle collaboration in 2003 1], more and more exotic hadrons have been observed and confirmed experimentally, such as the charmonium-like X, Y, Z states, hidden-charm pentaquarks, etc 2]. Those exotic hadron states, which cannot be interpreted as the quark-antiquark mesons or three-quark baryons in the naive quark model 3], are good candidates of the multi-quark states. The multi-quark states are color-neutral objects because of the color confinement, and provide an important platform to deepen our understanding of the low energy behaviors of QCD, as no free particles carrying net color charges have ever been experimentally observed. Compared to the conventional hadrons, the dynamics of the multi-quark states is poorly understood and calls for more works. Some exotic hadrons can be understood as hadronic molecular states, which are analogous to the deuteron as a loosely bound state of the proton and neutron. The most impressive example is the originalexotic state, the X(3872), which has been studied as the D D molecularstate by many theoretical groups 4]. Another impressive example is the P c (4380) and P c (4450) pentaquark states observed by the LHCb collaboration in 2015, which are good candidates for the DΣ c, D Σ c, D Σ c molecular states 4]. In additional to the meson-meson type and meson-baryon type molecular state, there maybe also exist meson-meson-meson type molecular states, in other words, there maybe exist three-meson hadronic molecules. In Refs.5, 6], the authors explore the possible existence of three-meson system D D K molecule according to the attractive interactions of the two-body subsystems DK, DK, D K, D K and D D with the Born-Oppenheimer approximation and the fixed center approximation, respectively. In this article, we study the D D K system with QCD sum rules. The QCD sum rules method is a powerful tool in studying the exotic hadrons, and has given many successful descriptions, for example, the mass and width of the Z c (3900) have been successfully reproduced as an axialvector tetraquark state 7, 8]. In QCD sum rules, the operator product expansion is used to expand the time-ordered currents into a series of quark and gluon condensates which parameterize the non-perturbative properties of the QCD vacuum. According to the quark-hadron duality, we can obtain copious information about the hadronic parameters at the phenomenological side 9, 10]. In this article, we construct the color singlet-singlet-singlet interpolating current to study the D D K system. In calculations, we take into account the contributions of the vacuum condensates E-mail: zgwang@aliyun.com. 1
up to dimension-16 in the operator product expansion and use the energy-scale formula µ = M 2 X/Y/Z (2M c) 2 to determine the ideal energy scale of the QCD spectral density. The rest of this article is arranged as follows: in Sect.2, we derive the QCD sum rules for the mass and pole residue of the D D K state; in Sect.3, we present the numerical results and discussions; Sect.4 is reserved for our conclusion. 2 QCD sum rules for the D D K state where In QCD sum rules, we consider the two-point correlation function, Π µν (p) = i d 4 xe ip x 0 T { J µ (x)j ν (0)} 0, (1) J µ (x) = ū m (x)iγ 5 c m (x) c n (x)γ µ d n (x)ū k (x)iγ 5 s k (x), (2) the m, n, k are color indexes. The color singlet-singlet-singlet current operator J µ (x) has the same quantum numbers as the D D K system. On 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 9, 10]. After isolating the ground state contribution Z from the pole terms, we get the following result: Π µν (p) = λ 2 Z M 2 Z p2 ( g µν + p µp ν = Π ( p 2)( g µν + p µp ν p 2 p 2 ) + ) +, (3) where the pole residue λ Z is defined by 0 J µ (0) Z(p) = λ Z ε µ, the ε µ is the polarization vector of the vector hexaquark molecular state Z. At the quark level, we calculate the correlation function Π µν (p) via the operator product expansion method in perturbative QCD. We contract the u, d, s and c quark fields with the Wick theorem, and obtain the result: Π µν (p) = i d 4 xe ip x { ] ] Tr γ µ D nn (x)γ ν C n n ( x) Tr iγ 5 C mm (x)iγ 5 U m m ( x) ] Tr γ µ D nn (x)γ ν C n n ( x) Tr Tr ] iγ 5 S kk (x)iγ 5 U k k ( x) ]}, (4) iγ 5 C mm (x)iγ 5 U m k ( x)iγ 5 S kk (x)iγ 5 U k m ( x) wheretheu ij (x), D ij (x), S ij (x)andc ij (x)arethefullu,d, sandcquarkpropagators,respectively, 2
(for simplicity, the U ij (x) and D ij (x) can be written as P ij (x)), P ij (x) = iδ ijx/ 2π 2 x 4 δ ij qq δ ijx 2 qg s σgq ig sg n αβ tn ij (x/σαβ +σ αβ x/) 12 192 32π 2 x 2 1 8 q jσ αβ q i σ αβ +, (5) S ij (x) = iδ ijx/ 2π 2 x 4 δ ijm s 4π 2 x 2 δ ij ss + iδ ijx/m s ss δ ijx 2 sg s σgs 12 48 192 C ij (x) = + iδ ijx 2 x/m s sg s σgs 1152 ig sg n αβ tn ij (x/σαβ +σ αβ x/) 32π 2 x 2 1 8 s jσ αβ s i σ αβ +, (6) { i (2π) 4 d 4 ke ik x k/+mc k 2 m 2 δ ij g s t n (k/+m c )σ αβ +σ αβ (k/+m c ) ij Gn αβ c 4(k 2 m 2 c) 2 g2 s(t n t m ) ij G n αβ Gn µν(f αβµν +f αµβν +f αµνβ } ) 4(k 2 m 2 c) 5 +, (7) f λαβ = (k/+m c )γ λ (k/+m c )γ α (k/+m c )γ β (k/+m c ), f αβµν = (k/+m c )γ α (k/+m c )γ β (k/+m c )γ µ (k/+m c )γ ν (k/+m c ), (8) and t n = λn 2, the λn is the Gell-Mann matrix 10]. Then we compute the integralsin the coordinate space for the light quark propagators and in the momentum space for the charm quark propagators and obtain the QCD spectral density ρ(s) via taking the imaginary part of the correlation ImΠ(s+iε) function: ρ(s) = lim ε 0 π 7]. In the operator product expansion, we take into account the contributions of vacuum condensates up to dimension-16, and keep the terms which are linear in the strange quark mass m s. We take the truncation k 1 for the operators of the order O(α k s) in a consistent way and discard the perturbative corrections. Furthermore, the condensates qq αsgg π, qq 2 αsgg π and qq 3 αsgg π play a minor important role and are neglected. According to the quark-hadron duality, we match the correlation function Π(p 2 ) obtained at the hadron level and at the quark level 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: ( ) λ 2 Z exp M2 s0 ( Z T 2 = dsρ(s)exp s ) T 2, (9) where the QCD spectral density ρ(s) = ρ 0 (s)+ρ 3 (s)+ρ 4 (s)+ρ 5 (s)+ρ 6 (s)+ρ 8 (s)+ρ 9 (s)+ρ 10 (s)+ρ 11 (s) +ρ 12 (s)+ρ 13 (s)+ρ 14 (s)+ρ 16 (s), (10) the subscripts 0, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14, 16 denote the dimensions of the vacuum condensates, the explicit expressions are neglected for simplicity, the T 2 is the Borel parameter. We derive Eq.(9) with respect to 1 T, and eliminate the pole residue λ 2 Z to extract the QCD sum rule for the mass: M 2 Z = s0 s0 3 Numerical results and discussions d dsd( 1/T 2 ) ρ(s)exp( s dsρ(s)exp ( ) s. (11) T 2 In this section, we perform the numerical analysis. To extract the numerical values of M Z, we take the standard values of the vacuum condensates qq = (0.24±0.01GeV) 3, ss = (0.8 ± T 2 ) 3
0.1) qq, qg s σgq = m 2 0 qq, sg sσgs = m 2 0 ss, m2 0 = (0.8±0.1)GeV2, αsgg π = (0.33GeV) 4 at the energy scale µ = 1GeV 9, 10, 11], choose the MS masses m c (m c ) = (1.275±0.025)GeV, m s (µ = 2GeV) = (0.095 +0.009 0.003 )GeVfromtheParticleDataGroup2], andneglecttheupanddown quark masses, i.e., m u = m d = 0. Moreover, we take into account the energy-scale dependence of the input parameters on the QCD side from the renormalization group equation, ]4 αs 9 (Q) qq (µ) = qq (Q), ]4 αs 9 (Q) ss (µ) = ss (Q), αs (Q) qg s σgq (µ) = qg s σgq (Q) αs (Q) sg s σgs (µ) = sg s σgs (Q) m s (µ) = m s (2GeV) αs (µ) m c (µ) = m c (m c ) α s (m c ) = 1 b 0 t 1 b 1 b 2 0 ] 2 27, ] 2 27, α s (2GeV) logt t ]12 25, + b2 1 ]4 9, ( log 2 t logt 1 ) +b 0 b 2 b 4 0 t2 ], (12) where t = log µ2 Λ, b 2 0 = 33 2n f 12π, b 1 = 153 19n f 24π, b 2 2 = 2857 5033 9 n f+ 325 27 n2 f 128π, Λ = 213MeV, 296MeV 3 and 339MeV for the flavors n f = 5, 4 and 3, respectively 2]. For the hadron mass, it is independent of the energy scale because of its observability. However, in calculations, the perturbative corrections are neglected, the operators of the orders O n (α k s ) with k > 1 or the dimensions n > 16 are discarded, and the higher dimensional vacuum condensates are factorized into lower dimensional ones therefore the energy-scale dependence of the higher dimensional vacuum condensates is modified. We have to consider the energy-scale dependence of the QCD sum rules. In Refs.7, 12, 13, 14], the authors study the energy-scale dependence of the QCD sum rules for the hidden-charm tetraquark states and molecular states in detail, and suggest an energy scale formula µ = MX/Y/Z 2 (2M c) 2 to determine the optimal energy scale. The energy-scale formula enhances the pole contribution remarkably, improves the convergent behaviors in the operator product expansion, and works well for the exotic hadron states. In this article, we study the D D K state Z through constructing the color singlet-singlet-singlet type current based on the color-singlet q q substructure. For the two-meson molecular states, the basic constituent is also the color-singlet q q substructure 14]. Hence, we extend the previous works to study the D D K system by using the energy-scale formula µ = MX/Y/Z 2 (2M c) 2 with the updated value of the effective c-quark mass M c = 1.85GeV to determine the ideal energy scale of the QCD spectral density. At the present time, there is no candidate observed experimentally for the hexaquark state Z with the symbolic quark constituent c cdūsū. While in the scenario of four-quark states, the QCD sum rules indicate that the Z c (3900) and Z(4430) can be tentatively assigned to be the ground state and the first radial excited state of the axialvector four-quark states, respectively 15], while the X(3915) and X(4500) can be tentatively assigned to be the ground state and the first radial excited state of the scalar four-quark states, respectively 16]. The energy gap between the ground state and the first radial excited state of the hidden-charm four-quark states is about 0.6 GeV. 4
PC 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 s 0 =(4.9GeV) 2 s 0 =(5.0GeV) 2 0.2 s 0 =(5.1GeV) 2 s 0 =(5.2GeV) 2 0.1 s 0 =(5.3GeV) 2 s 0 =(5.4GeV) 2 0.0 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 T 2 (GeV 2 ) Z Figure 1: The pole contribution with variation of the Borel parameter T 2. Here, we suppose the energy gap between the ground state and the first radial excited state of the hidden-charm six-quark states is also about 0.6 GeV, and take the continuum threshold parameter to be s 0 = M Z +(0.4 0.6)GeV, which serves as a constraint to obey. In Eq.(11), there are two free parameters: the Borel parameter T 2 and the continuum threshold parameter s 0. The extracted hadron mass is a function of the Borel parameter T 2 and the continuum threshold parameter s 0. To obtain a reliable mass sum rule analysis, we obey two criteria to choose suitable working ranges for the two free parameters. The first criterion is the pole dominance on the phenomenological side, which requires the pole contribution (PC) to be about (40 60)%. The PC is defined as: PC = s0 dsρ(s)exp ( s T 2 ) dsρ(s)exp ( s T 2 ). (13) The second criterion is the convergence of the operator product expansion. To judge the convergence, we calculate the contributions of the vacuum condensates D(n) in the operator product expansion with the formula: s0 dsρ 4m D(n) = 2 n (s)exp ( ) s c T 2 s0 dsρ(s)exp ( ) s, (14) T 2 where the index n denotes the dimension of the vacuum condensates. In Fig. 1, we show the variation of the PC with respect to the Borel parameter T 2 for different values of the continuum threshold parameter s 0 at the energy scale µ = 2.9GeV. From the figure, we can see that the value s 0 5.0GeV is too small to satisfy the pole dominance condition and result in reasonable Borel window for the state Z. To warrant the Borel platform for the mass m Z, we take the value T 2 = (2.8 3.2)GeV 2. In the above Borel window, if we choose the value s0 = (5.1 5.3)GeV, the PC is about (39 63)%. The pole dominance condition is well satisfied. In Fig.2, we plot the absolute contributions of the vacuum condensates D(n) in the operator product expansion at central values of the above input parameters. From the figure, we can 5
1.0 0.8 0.6 D(n) 0.4 0.2 0.0-0.2-2 0 2 4 6 8 10 12 14 16 18 Figure 2: The absolute contributions of the vacuum condensates with dimension n in the operator product expansion. n M(GeV) 7.8 7.2 6.6 6.0 5.4 4.8 4.2 3.6 Z 3.0 2.80 2.85 2.90 2.95 3.00 3.05 3.10 3.15 3.20 T 2 (GeV 2 ) Central value Upper bound Lower bound Figure 3: The mass with variation of the Borel parameter T 2. 6
(10-4 GeV 8 ) 24 21 18 15 12 9 6 3 0 Z -3 2.80 2.85 2.90 2.95 3.00 3.05 3.10 3.15 3.20 T 2 (GeV 2 ) Central value Upper bound Lower bound Figure 4: The pole residue with variation of the Borel parameters T 2. observe that the contribution of the perturbative term D(0) is not the dominant contribution, the contributions of the vacuum condensates of dimensions 3, 6, 8, 9 and 11 are very large. If we take the contribution of the vacuum condensate of dimension 11 as a milestone, the contributions of the vacuum condensates D(n) decrease quickly with the increase of the dimensions n, the operator product expansion is well convergent. Thus, we obtain the values T 2 = (2.8 3.2)GeV 2, s 0 = (5.1 5.3)GeV and µ = 2.9GeV for the state Z. Taking into account all uncertainties of the input parameters, we obtain the values of the mass and pole residue of the state Z: M Z = 4.71 +0.19 0.11 GeV, λ Z = ( 4.60 +1.15 0.69) 10 4 GeV 8, (15) which are shown explicitly in Figs.3 4. The energy-scale formula µ = MX/Y/Z 2 (2M c) 2 and the relation s 0 = M Z +(0.4 0.6)GeV are also well satisfied. The central value M Z = 4.71GeV is about 337MeV above the threshold M K+D+ D = 497.6 + 1865 + 2010 = 4372.6MeV, which indicates that the Z is probably a virtual state, or a resonance state. The virtual state can not be observed experimentally, while the resonance state may be observed experimentally in the future. For some exotic resonances, the authors have combined the effective range expansion, unitarity, analyticity and compositeness coefficient to probe their inner structure in Refs.17, 18]. Their studies indicated that the underlying two-particle component (in the present case, corresponding to three-particle component) play an important or minor role, in other words, there are the other hadronic degrees of freedom inside the corresponding resonance. Hence, a resonance state embodies the net effect. 4 Conclusion In this article, we construct the color singlet-singlet-singlet interpolating current operator to study the D D K system with QCD sum rules by taking into account the contributions of the 7
vacuum condensates up to dimension-16 in the operator product expansion. In numerical calculations, we saturate the hadron side of the QCD sum rule with a hexaquark molecular state, use the energy-scale formula µ = MX/Y/Z 2 (2M c) 2 to determine the energy scale of the QCD spectral density, and search for the ideal Borel parameter T 2 and continuum threshold s 0 to satisfy the pole dominance and convergent of the operator product expansion. Finally, we obtain the mass and pole residue of the corresponding hexaquark molecular state Z. The predicted mass, M Z = 4.71 +0.19 0.11 GeV, which lies above the D D K threshold, indicates that the Z is probably a virtual state, or a resonance state. The virtual state can not be observed experimentally, while the resonance state, which embodies the net effect, may be observed experimentally in the future. Acknowledgements This work is supported by National Natural Science Foundation, Grant Number 11775079. References 1] S. K. Choi, et al., Phys. Rev. Lett. 91 (2003) 262001. 2] C. Patrignani, et al., Chin. Phys. C40 (2016) 100001. 3] S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189. 4] F. K. Guo, C. Hanhart, U. G. Meissner, Q. Wang, Q. Zhao and B. S. Zou, Rev. Mod. Phys. 90 (2018) 015004. 5] L. Ma, Q. Wang and U.-G. Meissner, arxiv:1711.06143. 6] X. L. Ren, B. B. Malabarba, L. S. Geng, K. P. Khemchandani and A. M. Torres, Phys. Lett. B785 (2018) 112. 7] Z. G. Wang and T. Huang, Phys. Rev. D89 (2014) 054019. 8] Z. G. Wang and J. X. Zhang, Eur. Phys. J. C78 (2018) 14. 9] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 385, 448. 10] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 1. 11] P. Colangelo and A. Khodjamirian, At the Frontier of Particle Physics: Handbook of QCD, Vol. 3, ed. M. Shifman (World Scientific, Singapore, 2001), p. 1495. 12] Z. G. Wang, Eur. Phys. J. C74 (2014) 2874. 13] Z. G. Wang and T. Huang, Nucl. Phys. A930 (2014) 63. 14] Z. G. Wang and T. Huang, Eur. Phys. J. C74 (2014) 2891; Z. G. Wang, Eur. Phys. J. C74 (2014) 2963. 15] Z. G. Wang, Commun. Theor. Phys. 63 (2015) 325. 16] Z. G. Wang, Eur. Phys. J. C77 (2017) 78; Z. G. Wang, Eur. Phys. J. A53 (2017) 19. 17] X. W. Kang, Z. H. Guo and J. A. Oller, Phys. Rev. D94 (2016) 014012. 18] R. Gao, Z. H. Guo, X. W. Kang and J. A. Oller, arxiv:1812.07323. 8