The New Charm-Strange Resonances in the Channel
Abstract
We evaluate the masses and decay constants of the and open-charm tetraquarks and molecular states from QCD spectral sum rules (QSSR) by using QCD Laplace sum rule (LSR). This method takes into account the stability criteria where the factorised perturbative NLO corrections and the contributions of quark and gluon condensates up to dimension-6 in the OPE are included. We confront our results with the invariant mass recently reported by LHCb from decays. We expect that the resonance near the threshold can be originated from the molecule and/or scattering. The scalar state and the resonance (if ) can emerge from a minimal mixing model, with a tiny mixing angle , between a scalar Tetramole (superposition of nearly degenerated hypothetical molecules and compact tetraquarks states with the same quantum numbers), having a mass MeV, and the first radial excitation of the molecule with mass MeV. In an analogous way, the and the (if ) could be a mixture between the vector Tetramole , with a mass MeV, and its first radial excitation having a mass MeV with an angle . A (non)-confirmation of these statements requires experimental findings of the quantum numbers of the resonances at and MeV.
keywords:
QCD sum rules , Perturbative and non-perturbative QCD , Exotic hadrons , Masses and decay constants.1 Introduction
In this work, based on the paper in Ref. [1], we attempt to estimate, from LSR, the masses and couplings of the and molecules and compact tetraquarks states for interpreting the recent LHCb data from decays [2, 3], where one finds two prominent peaks (units of MeV):
We have studied in Ref. [4] the masses and couplings of the molecule and of the corresponding tetraquark states decaying into but not into and we found the lowest ground state masses:
We have used this result to interpret the nature of the (2317) compiled by PDG [5] where the existence of a pole at this energy has been recently confirmed from lattice calculations of scattering amplitudes [6].
Scalar states () | Vector states () |
Tetraquarks | |
Molecules | |
For the molecular state, we can interchange the and quarks in the interpolating current and deduce from SU(2) symmetry that the molecule mass is degenerated with the one. Compared with the LHCb data, one may invoke that this charged molecule can be responsible of the bump near the DK threshold around 2.4 GeV but is too light to explain the peaks.
2 The Laplace sum rule (LSR)
We shall work with the Finite Energy version of the QCD Inverse Laplace sum rules (LSR) and their ratios [15, 16, 17, 18, 19, 20, 21, 22, 23, 25, 24, 26, 27]:
(1) | |||||
(2) |
where and are the on-shell / pole charm and running strange quark masses, is the LSR variable, is the degree of moments, is the threshold of the “QCD continuum” which parametrizes, from the discontinuity of the Feynman diagrams, the spectral function is evaluated by the calculation of the scalar correlator defined as:
(3) |
where are the interpolating currents for the tetraquarks and molecules states. The superscript refers to the spin of the particles.
2.1 The Interpolating Operators
We shall be concerned with the interpolating given in Table 1. The lowest order (LO) perturbative (PT) QCD expressions - including the quark and gluon condensates contributions up to dimension-six condensates of the corresponding two-point spectral functions - the NLO PT corrections, the QCD input parameters and further details of the QSSR calculations for those interpolating operators are given in the Ref.[1].
3 Tetraquarks and Molecules
The sum rule analysis for the and states present similar features for all currents in Table 1. Then, we show only explicitly the analysis of the tetraquark channel for a better understanding on the extraction of our results.
3.1 - and -stabilities
We show in Fig.1a) the - and - dependence of the mass obtained from ratio of moments . The analysis of the coupling from the moment is shown in Fig. 1b). The results stabilize at GeV-2 (inflexion point for the mass and minimum for the coupling).
From Fig.1a), we extract the mass as the mean value of the one for 12 GeV2 (beginning of the inflexion point) and of the one at beginning of -stability of about 18 GeV2. We use this (physical) mass value in to draw Fig. 1b). We check the range of -values where the above-mentioned stabilities have been obtained by confronting Figs. 1a) and b). Here, one can easily check that this range of -values is the same for the mass and coupling. If the range does not coincide, we take the common range of and redo the extraction of the mass.
One can also see that the range of -stabilities coincide in Fig. 1a) (inflexion points) and in Fig. 1b) (minimas). It is obvious that the value of from the minimum is more precise which we re-use to fix the final value of the mass.
a)

b)

3.2 -stability
In Fig. 2, we show the -dependence of the results for given =18 GeV2 and =0.49 GeV-2. One finds a common stability for GeV, which confirms the result in Ref. [4].
a)

b)

4 The First Radial Excitation
We extend the analysis in Ref. [4] by using a “Two resonances” + “QCD continuum” parametrization of the spectral function. To enhance the contribution of the 1st radial excitation, we shall also work with the ratio of moments in addition to for getting the masses. We use the same criteria involving the stability points in and the optimal results are given in Table 2. We observe that the mass-splittings between the first radial excitation and the lowest ground state are in order of MeV, which is much bigger than MeV typically used for ordinary mesons.
5 Understanding LHCb Experimental Data
Our results indicate that the molecules and tetraquark states leading to the same final states are almost degenerated in masses. Therefore, we expect that the “physical state” is a combination of almost degenerated molecules and tetraquark states with the same quantum numbers which we shall call: Tetramole .
5.1 The and states
Taking literally our results in Table 2, one can see that we have three (almost) degenerate states:
and their couplings to the corresponding currents are almost the same:
We assume that the physical state, hereafter called Tetramole , is a superposition of these nearly degenerated hypothetical states having the same quantum numbers. Taking its mass and coupling as (quadratic) means of the previous numbers, we obtain:
The tetramole is a good candidate for explaining the though its mass is slightly lighter.
One can also see from Table 2 that the radial excitation mass and coupling are :
which is the lightest first radial excitation. Assuming that the bump is a scalar state (), we attempt to use a two-component minimal mixing model between the Tetramole and the radially excited molecule:
We reproduce the data with a tiny mixing angle :
5.2 The and states
As one can see in Table 2, there are four degenerate states with masses around MeV, and couplings around keV. We assume again that the (unmixed) physical state is a combination of these hypothetical states. We evaluate the mass and coupling of this Tetramole as the (geometric) means:
where one may notice that it can contribute to the state but its mass is slightly lower. Looking at to the radial excitations in Table 2, one can see that they are almost degenerated around 4.5 GeV from which one can extract the masses and couplings (geometric mean) of the spin 1 Tetramoles:
Then, we may consider a minimal two-component mixing of the spin 1 Tetramole () with its 1st radial excitation to explain the state and the bump assuming that the latter is a spin 1 state. The data can be fitted with a tiny mixing angle :
A (non)-confirmation of these two minimal mixing models requires an experimental identification of the quantum numbers of the bumps at 3150 and 3350 MeV.
5.3 Final results
Our final results are obtained at the stability points of the set of parameters and they are summarized in Table 2. One can notice that, for some molecule and tetraquark states, the ground state mass values are above GeV which are too far to contribute to the LHCb observations in invariant mass. In such cases, the sum rule results are discarded. One can find the full analysis of different sources of errors, as well as an interesting discussion on the relevance of NLO calculations for sum rules in Ref.[1].
Observables | (MeV) | (keV) | (MeV) | (keV) |
---|---|---|---|---|
States | ||||
Molecule | ||||
2402(42) | 254(48) | 3678(310) | 211(51) | |
2808(41) | 405(33) | 4626(252) | 568(167) | |
5258(113) | 664(57) | |||
6270(160) | 249(18) | |||
Tetraquark | ||||
2736(21) | 345(28) | 4586(268) | 359(81) | |
2675(65) | 498(43) | 4593(289) | 547(95) | |
5704(149) | 713(66) | |||
5917(98) | 538(41) | |||
States | ||||
Molecule | ||||
2676(47) | 191(21) | 4582(414) | 157(71) | |
2744(41) | 216(22) | 4662(269) | 237(63) | |
5377(166) | 351(31) | |||
5358(153) | 255(23) | |||
Tetraquark | ||||
2666(32) | 285(29) | 4571(213) | 258(82) | |
2593(31) | 259(25) | 4541(345) | 243(68) | |
5542(139) | 416(38) | |||
5698(175) | 412(43) |
6 Summary and conclusions
Motivated by the recent LHCb data on the invariant mass from decay, we have systematically calculated the masses and couplings of some possible configurations of the molecules and tetraquarks states using QCD Laplace sum rules (LSR) within stability criteria where we have added to the LO perturbative term, the NLO radiative corrections, and the contributions from quark and gluon condensates up to dimension-6 in the OPE.
The peak around the threshold can be due to scattering amplitude the lowest mass molecule.
The () and (if it is a state) can e.g result from a mixing of the Tetramole () with the 1st radial excitation of the molecule state with a tiny mixing angle .
The () and (if it is a state) can result from a mixing of the Tetramole () with its 1st radial excitation with a tiny mixing angle .
More data on the precise quantum numbers of the and states are needed for testing the previous two minimal mixing models proposal.
References
- [1] R.M. Albuquerque, S. Narison, D. Rabetiarivony and G. Randriamanatrika, Nucl. Phys. A 1007 (2021) 122113.
- [2] [LHCb Collab.] R. Aaij et al., Phys.Rev.Lett. 125, 242001 (2020).
- [3] [LHCb Collab.], R. Aaij et al., Phys.Rev. D102, 112003 (2020).
- [4] R. M. Albuquerque et al., Int. J. Mod. Phys. A 31 (2016) 17, 1650093.
- [5] M. Tanabashi et al. (Particle Data Group), Phys. Rev. D 98 (2018) 030001 and 2019 update.
- [6] G.K.C. Cheung et al., arXiv: 2008.06432v1 [hep-lat] (2020).
- [7] H. G. Dosch, M. Jamin and B. Stech, Z. Phys. C42 (1989) 167.
- [8] M. Jamin and M. Neubert, Phys. Lett. B238 (1990) 387.
- [9] J.-R. Zhang, arXiv: 2008.07295 [hep-ph] (2020).
- [10] H.-X. Chen et al., Chin. Phys. Lett. 37 (2020) 101201.
- [11] Z.-W. Wang, Int. J. Mod. Phys. A35 (2020) 2050187.
- [12] W. Chen et al., Phys. Rev. D95 (2017) 114005.
- [13] S.S. Agaev, K. Azizi and H. Sundu, arXiv: 2008.13027 [hep-ph] (2020).
- [14] H. Mutuk, arXiv: 2009.02492 [hep-ph] (2020).
- [15] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B147 (1979) 385.
- [16] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B147 (1979) 448.
- [17] S. Narison, QCD spectral sum rules, World Sci. Lect. Notes Phys. 26 (1989) 1.
- [18] S. Narison, QCD as a theory of hadrons, Cambridge Monogr. Part. Phys. Nucl. Phys. Cosmol. 17 (2004) 1-778 [hep-ph/0205006].
- [19] S. Narison, Phys. Rept. 84 (1982) 263;
- [20] S. Narison, Acta Phys. Pol. B 26(1995) 687;
- [21] B.L. Ioffe, Prog. Part. Nucl. Phys. 56 (2006) 232.
- [22] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 1.
- [23] E. de Rafael, les Houches summer school, hep-ph/9802448 (1998).
- [24] F.J Yndurain, The Theory of Quark and Gluon Interactions, 3rd edition, Springer (1999).
- [25] R.A. Bertlmann, Acta Phys. Austriaca 53, (1981) 305.
- [26] P. Pascual and R. Tarrach, QCD: renormalization for practitioner, Springer 1984.
- [27] H.G. Dosch, Non-pertubative Methods, ed. Narison, World Scientific (1985).