Abstract
In this study, we investigate the , , , , , , , and pentaquark molecular states with and without strangeness via the QCD sum rules in detail, focusing on the light flavor, , breaking effects, and make predictions for new pentaquark molecular states besides assigning , , , , and self-consistently. In the future, we can search for these pentaquark molecular states in the decay of , , and . Furthermore, we discuss high-dimensional vacuum condensates in detail.
Export citation and abstract BibTeX RIS
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Article funded by SCOAP and published under licence by Chinese Physical Society and the Institute of High Energy Physics of the Chinese Academy of Sciences and the Institute of Modern Physics of the Chinese Academy of Sciences and IOP Publishing Ltd
I. INTRODUCTION
In 2015, the LHCb collaboration explored the decay and observed two pentaquark candidates and in the mass spectrum with preferred quantum numbers and , respectively [1]. The Breit-Wigner masses and widths are , MeV, MeV, and MeV, respectively. In 2019, the LHCb collaboration re-investigated a data sample of the decay, which was an order of magnitude larger than that previously analyzed, and observed a narrow pentaquark candidate in the mass spectrum, and confirmed structure , which consisted of two narrow overlapping peaks and [2]. The measured Breit-Wigner masses and widths are
Very recently, the LHCb collaboration reported an evidence of a hidden-charm pentaquark candidate with strangeness in the invariant mass spectrum with a statistical significance of in the decay using the collision data corresponding to a total integrated luminosity of collected in LHCb experiments at centre-of-mass energies of 7, 8, and 13 TeV [3]. Its Breit-Wigner mass and width are
but the spin and the parity have not been determined yet. , , , , and lie slightly below the thresholds of meson-baryon pairs , , , and , respectively. The nearby meson-baryon thresholds are listed clearly in Table 1, because the D and mesons, flavor antitriplet ( , , and ), and and flavor sextets ( , , ) and ( , , ) have been well established [4].
Table 1. The thresholds of the open-charm meson and baryon pairs, where the unit is MeV.
( , ) | (4336, 4337) |
---|---|
( , ) | (4478, 4478) |
( , ) | (4318, 4323) |
( , ) | (4382, 4388) |
( , ) | (4460, 4464) |
( , ) | (4524, 4529) |
( , ) | (4444, 4448) |
( , ) | (4511, 4515) |
( , ) | (4586, 4589) |
( , ) | (4653, 4656) |
As expected, , , , , and have been tentatively assigned to the , , , , and pentaquark molecular states, respectively, based on the contact-range effective field theory [5, 6], one-boson exchange potential model [7-9], (quasipotential) Bethe-Salpeter equation [10-13], effective Lagrangian approach [14, 15], effective-range expansion and resonance compositeness relation [16], Lippmann-Schwinger equation [17, 18], and QCD sum rules [19-23].
At first glance, and , and it is unusual that the molecular state, which involves exchanges of strange mesons, is more tightly bound than the molecular state, which involves exchanges of non-strange mesons. Therefore, we have to introduce coupled-channel effects [7-13, 16-18].
In the QCD sum rules, we usually choose local currents to interpolate tetraquark or pentaquark molecular states having two color-neutral clusters [19-23], which are not necessary to be physical mesons and baryons, and the tetraquark or pentaquark molecular states are not necessary to be loosely bound; they can be compact objects and can lie below or above the corresponding meson-meson or meson-baryon pairs [24].
In contrast, , , , , and can also be tentatively assigned to diquark-diquark-antiquark-type (or diquark-triquark-type) pentaquark states in the diquark-model by exploring their masses [25, 26] ([27]) and decay modes [28-30] ([31]) using an effective Hamiltonian or by investigating their masses [32-35, 37], decay [38], and electromagnetic properties [39] using the QCD sum rules.
In Ref. [40], we suggested the hadronic dressing mechanism to compromise pentaquark and pentaquark molecular interpretations based on the calculations of the QCD sum rules and that pentaquark states may have a diquark-diquark-antiquark-type pentaquark core with the typical size of -type baryon states. Moreover, the strong couplings to meson-baryon pairs lead to some pentaquark molecule Fock components; therefore, pentaquark states may spend a rather long time as molecular states. We can choose either diquark-diquark-antiquark-type currents or color-singlet-color-singlet-type five-quark currents to interpolate pentaquark states.
In the scenario of pentaquark molecular states, , , , and are assigned to the , , , and pentaquark molecular states, which involve charmed baryon states in flavor sextets [5, 7, 10, 14, 16-22], whereas is assigned to the pentaquark molecular state, which involves a charmed baryon state in flavor antitriplet [6, 8, 9, 11-13, 15, 23]. Thus the and belong to different flavor multiplets, and in the present study, we focus on the flavor sextets, .
In this study, we extend our previous research [22] to investigate the masses and pole residues of the , , , , , , , and pentaquark molecular states with the QCD sum rules by accomplishing operator product expansion up to vacuum condensates of dimension 13 in a consistent manner, take into account vacuum condensates , , and neglected in Ref. [22], and revisit the assignments of , , , , and . Furthermore, we pay much attention to the light flavor, , breaking effects.
The remainder of this article is arranged as follows: we describe the QCD sum rules for the masses and pole residues of pentaquark molecular states in Section II. In Section III, we present the numerical results and discussions, and Section IV draws our conclusion.
II. QCD SUM RULES FOR THE PENTAQUARK MOLECULAR STATES
We express the two-point correlation functions, , , and , in the QCD sum rules,
where currents , , , , , , , ,
where i, j, and k are color indices. In the present study, we choose color-singlet-color-singlet-type currents , , and to interpolate the pentaquark molecular states with spin parities , , and , respectively. The currents couple potentially to the pentaquark molecular states having two color-neutral clusters: one has the same quantum numbers as the charmed mesons, and the other has the same quantum numbers as the charmed baryons. They are not physical mesons and baryons, as we choose local five-quark currents, whereas the mesons and the baryons are spatial extended objects and have mean spatial sizes . For example, , , and from the lattice QCD, where subscripts E and M stand for the electric and magnetic radii, respectively[41], , , , , from the self-consistent chiral quark-soliton model [42], , and from the light-front quark model [43]. In the present study, although we refer the color-singlet-color-singlet-type pentaquark states as the pentaquark molecular states, they have average spatial sizes as those of the typical heavy mesons and baryons, and are compact objects. For example, a loosely bound molecular state with physical meson and baryon should have average spatial size , which is extremely large to be interpolated by local currents.
Currents , and couple potentially to , , and , , hidden-charm pentaquark molecular states , , and , , , respectively,
where the , , and are the Dirac and Rarita-Schwinger spinors [22, 32-36].
At the hadron side of correlation functions , , and , we isolate the ground state contributions from the hidden-charm pentaquark molecular states with spin parities , , and , respectively, without contamination based on the current-hadron couplings expressed in Eqs. (8)-(10) and obtain the hadronic representation [22, 32-36],
There are other spinor structures, which are not shown explicitly, and we choose the components corresponding to spinor structures , 1, , , and , in correlation functions , , and , respectively, to investigate the , , and pentaquark molecular states.
Now, we briefly digress to discuss the isospins of the interpolating currents. From Eqs. (4)-(7), we can see clearly that currents , , , and without strangeness have the same isospin structures, and currents , , , and with strangeness also have the same isospin structures. Moreover, they can be transformed into each other with the simple replacement, . It is a good objective to explore the light flavor, , breaking effects. We can rewrite currents and in terms of the isospin eigenstates as
where
the , , , , and are the standard currents for the mesons and baryons, respectively, and superscripts , 1, and stand for the isospins of the interpolating currents. Now let us estimate the isospin breaking effects,
where current . The isospin breaking effects between the vacuum matrix elements, and , at the quark-gluon level are suppressed by a factor , which is of minor importance and can be neglected safely in the large limit. Then, we can estimate (or obtain the conclusion tentatively) that the currents and couple potentially to the spin- pentaquark molecular states with almost degenerated masses but different pole residues; therefore, we will not distinguish isospins and as we are only interested in the molecular masses, just like in previous studies [19-22]. Furthermore, current has the quantum numbers and , and we can add superscript to distinguish the components in the isospin triplets and singlets,
and they couple potentially to the pentaquark molecular states with almost degenerated masses. Again the isospin breaking effects are suppressed by the factor, ; therefore, we will not distinguish isospins and 0.
Now, let us go back to the correlation functions. It is straightforward to obtain the spectral densities at the hadron side through the dispersion relation,
where , , and , and we add subscript H to represent the hadron side; then, we introduce weight functions and to obtain the QCD sum rules at the hadron side,
where is the continuum threshold parameter and is the Borel parameter.
It is also straightforward to accomplish the operator product expansion in the deep Euclidian space-time. For the technical details of performing the operator product expansion for the correlation functions in exploring the multiquark states with hidden charm, one can consult Refs. [32-35, 44-49]. If we contract the quark fields in the correlation functions in Eq. (3) with Wick's theorem, we can observe clearly that there are two heavy quark propagators and three light quark propagators. If each heavy quark line emits a gluon and each light quark line contributes a quark-antiquark pair, we obtain a quark-gluon mixed operator with , d or s, which is of dimension 13, and leads to vacuum condensates and . It is better to take into account the vacuum condensates up to dimension at least. In contrast, the quark-gluon operators can be counted by fine structure constant with orders , where , , 1, , . In the present study, we take truncations and consistently, and the quark-gluon operators of orders with are given full considerations, whereas in the previous study [22], we neglected vacuum condensates , , and owing to their small contributions. Furthermore, we take into account the light flavor, , mass-breaking effects by including the contributions of order consistently.
Now, we briefly digress to discuss the higher dimensional vacuum condensates. In the QED, we deal with a perturbative vacuum, and the vacuum expectation values of normal-ordered electron-photon operators can be set as zero, for example, , , and .
In the QCD, we deal with a non-perturbative vacuum and have to resort to the non-zero vacuum expectation values of normal-ordered quark-gluon operators to describe the hadron properties satisfactorily, for example, , , and , where i, j, m , and n are color indices and α, β, λ , and τ are the Dirac spinor indices. We usually parameterize vacuum matrix elements in terms of
or , etc., where are the Gell-mann matrices. Except for the quark condensates, which indicate spontaneous breaking of the Chiral symmetry through the Gell-Mann-Oakes-Renner relation, [50], vacuum condensates such as , , , are just parameters introduced by hand to describe the non-pertubative vacuum.
We can parameterize the non-perturbative properties in some manner, and then compare them to the experimental data on multiquark states to obtain the optimal values. In the QCD sum rules for multiquark states, plays an important role and influences the convergent behaviors of the operator product expansion and pole contributions remarkably; therefore, remarkably influencing the predictions, and large values of ϱ may destroy the platforms [51]. For example, in the present case, if we take value in the QCD sum rules for the pentaquark molecular state, we obtain uncertainty , which is of the same order of the total uncertainty from other parameters, and is a very bad platform (in other words, no platform at all). In the calculations, we observe that the optimal value is , and vacuum saturation (factorizaiton) is effective in the QCD sum rules for multiquark states [32-37, 44-49, 51].
In the QCD sum rules for the , , and mesons, is typically accompanied with the fine-structure constant, , and plays a small role, and the deviation from vacuum saturation (factorization) , for example, , cannot make much difference in the numerical predictions, although in some cases the values can lead to better QCD sum rules [52, 53].
Once the corresponding analytical spectral densities, and , at the quark-gluon level are obtained, we can take the quark-hadron duality below the continuum threshold, , and introduce weight functions and to obtain the QCD sum rules as follows:
where the explicit expressions of spectral densities and at the quark level are neglected for simplicity.
We differentiate Eq. (21) with respect to and then eliminate pole residues with , , and to obtain the QCD sum rules for the masses of the pentaquark molecular states,
where spectral densities and .
III. NUMERICAL RESULTS AND DISCUSSIONS
We take the standard values of the vacuum condensates, , , , , , and , at the energy scale of [54-56], and take masses and GeV from the Particle Data Group [4]. Furthermore, we take into account the energy-scale dependence of the quark condensates, mixed quark condensates, and masses according to the renormalization group equation [57],
where
, , and for the flavors, , 4, and 3, respectively [4, 57].
In the present study, we investigate the hidden-charm pentaquark molecular states with and without strangeness, and it is good to choose flavor number . We evolve all input parameters to typical or special energy scales μ, which satisfy the energy scale formula or a modified energy scale formula [32-35, 44-49] with the updated value of the effective (or constituent) charmed quark mass, [58]. By comparing with the constituent quark masses based on analysis of the and mass spectrum with the well-known Cornell potential [59], we introduce an uncertainty . Furthermore, we take into account the light flavor, , mass-breaking effects, and prefer the modified energy scale formula, , with s-quark numbers , 1, 2, 3 and effective s-quark mass , which was proved to be effective [60]. Compared to the constituent quark mass, the effective s-quark mass, , seems extremely small, as the effective -quark masses serve as milestones and have been absorbed into energy scale μ, where value embodies the net mass-breaking effects.
We can rewrite the energy scale formula in the form,
where the light flavor mass-breaking effects δ have the value of and the Constants have value , and they are all fitted by the QCD sum rules. μ and embody the light degrees of freedom, whereas embodies the heavy degree of freedom. The hidden-charm tetraquark and pentaquark (molecular) states can be divided into heavy and light degrees of freedom [32-35, 45-49]. The predicted tetraquark and pentaquark (molecular) masses and the pertinent energy scales of the QCD spectral densities have Regge-trajectory-like relations [51].
In Ref. [22], we explored the , , , and pentaquark molecular states with the QCD sum rules in detail, and concluded that the energy scale formula, , can enhance the pole contributions at the hadron side remarkably and improve the convergent behaviors of the operator product expansion notably. In fact, we take the energy scale formula as a constraint on the predicted molecular masses, which should be obeyed in the QCD sum rules. In the present study, we consider the light flavor, , mass-breaking effects and utilize the modified energy scale formula, , to choose the best energy scales of the spectral densities at the quark-gluon level. We search for the best Borel parameters and continuum threshold parameters to satisfy the two fundamental criteria of the QCD sum rules: pole dominance at the hadron side and convergence of the operator product expansion at the QCD side, via trial and error.
Then, we obtain the Borel parameters, continuum threshold parameter , optimal energy scales of the spectral densities at the quark-gluon level, and pole contributions of the ground-state pentaquark molecular states, which are listed in Table 2. From the table, we can clearly see that the contributions from the ground states are about or larger than 40–60%, and the pole dominance criterion is satisfied very well. For the conventional hadrons, QCD spectral densities with and 2 for the mesons and baryons, respectively, and it is easy to satisfy the pole dominance criterion, because the integral,
Table 2. Optimal energy scale μ, Borel parameter , continuum threshold parameter , and pole contribution (pole) for hidden-charm pentaquark molecular states.
pole(%) | |||||
---|---|---|---|---|---|
2.2 | |||||
2.2 | |||||
2.4 | |||||
2.4 | |||||
2.5 | |||||
2.5 | |||||
2.6 | |||||
2.6 |
converges quickly even if we choose a large Borel parameter , where is the threshold, and the uncertainty originating from the continuum threshold parameter, , is small. For the multiquark states, QCD spectral densities with and 5 for the tetraquark and pentaquark (molecular) states, respectively, and it is very difficult to satisfy the pole dominance criterion, because the integral,
converges very slowly even if we choose a rather small Borel parameter . In general, we expect to choose , and the integral (or continuum state) is suppressed by a factor . Thus, for the multiquark states, we have to resort to a much stringent suppression of the continuum states, . One may think that such a small Borel parameter might lead to a bad convergent behavior in the operator product expansion.
In Fig. 1, we plot the absolute values of for the central values of the input parameters listed in Table 2, where the is defined as
the denotes the QCD spectral densities for the vacuum condensates of dimension n, and the total spectral densities, . From the figure, we can see clearly that although the largest contributions do not come from term in some cases, the vacuum condensates and with dimension six serve as a milestone, and the absolute values of the contributions, , with decrease monotonically and quickly with the increase in the dimensions, n. Moreover, the value , and the operator product expansion converges very well. The two basic criteria of the QCD sum rules are satisfied.
In the calculations, we observe that the predicted molecular masses increase monotonically and slowly with the increase in the continuum threshold parameter, , if we fix the Borel parameter, ; in contrast, a larger continuum threshold parameter implies a larger pole contribution. We truncate the continuum threshold parameter, , requiring the same pole contributions, approximately 40–60%, in all QCD sum rules so as to reduce the uncertainties originating from the continuum threshold parameter, .
In a previous study [22], we neglected vacuum condensates , , and , which are of dimensions 7, 10, and 13, respectively, owing to their small contributions. From Fig. 1, we can see clearly that the vacuum condensates of 7, 10, and 13 days play a small role in the Borel windows indeed. We prefer to take into account these contributions because they lead to slightly larger pole contributions, and therefore, more reliable QCD sum rules. In the present study, we intend to explore the -breaking effects, and it is better to take into account these contributions consistently.
Finally, we take into account all uncertainties of the input parameters, and obtain the masses and pole residues of the , , and hidden-charm pentaquark molecular states without and with strangeness, which are summarized explicitly in Table 3 and Fig. 2. In Fig. 2, we plot the predicted masses of the hidden-charm pentaquark molecules without and with strangeness according to variations in the Borel parameter, where the regions between two short vertical lines are the Borel windows. From the figure, we can see clearly that there appear rather flat platforms in the Borel windows and the uncertainties coming from the Borel parameters seem rather small. The uncetainities are compatible with the fact that the Borel parameters are just supplementary parameters, and not physical quantities. Furthermore, in the figure, we also present the experimental values of the masses of , , , , and from the LHCb collaboration [1-3].
Table 3. Predicted masses and pole residues of the hidden-charm pentaquark molecular states with the possible assignments, where double-? denotes that such assignment is not excluded owing to the uncertainties.
Thresholds/MeV | Assignments | ||||
---|---|---|---|---|---|
4318 | ? | ||||
4443 | ? | ||||
4382 | ? | ||||
4510 | ?? | ||||
4460 | ? | ||||
4585 | |||||
4524 | ?? | ||||
4652 |
The pentaquark (molecule) candidates , , , and are observed in the mass spectrum, and their isospins are , whereas pentaquark (molecule) candidate is observed in the mass spectrum, and its isospin is . The present calculations support assigning to the pentaquark molecular state with quantum numbers and , assigning as the pentaquark molecular state with quantum numbers and , assigning as the pentaquark molecular state with quantum numbers and , assigning the as the pentaquark molecular state with quantum numbers and . However, we cannot exclude the possibilities of assigning as the pentaquark molecular state with quantum numbers and and assigning the as the pentaquark molecular state with quantum numbers and owing to the uncertainties, and refer to Table 3 and Fig. 2. For example, it is marginal to assign the as the pentaquark molecular state with the quantum numbers and and as lies at the bottom of the predicted mass of the pentaquark molecular state, see Fig. 2-G.
From Tables 2-3, we can see that the modified energy scale formula, with s-quark numbers , 1, 2, 3 and the effective s-quark mass is satisfied very well [60]. In contrast, the predicted masses for the pentaquark molecular states without and with strangeness have the relation, , which is consistent with the light-flavor, , breaking effects for the heavy baryons in the flavor sextet, , from the Particle Data Group [4].
The present calculations indicate that there may exist ( ), ( ), ( ), and ( ) pentaquark molecular states with the , , and , respectively, which lie near the corresponding thresholds: ( ), ( ), ( ) , and ( ), respectively (see Table 3. The two-body strong decay to the corresponding open-charm meson-baryon pairs, such as ( ), ( ), ( ), and ( ), with the fall-apart mechanism directly, can only take place through the higher tails of the mass distributions and are kinematically suppressed in the phase space, and the widths of these pentaquark molecular states should be narrow. A large width may indicate that may correspond to two or more unresolved structures. More experimental data and theoretical studies are still needed to identify , , , , and , unambiguously.
In the present study, we predict the masses of new pentaquark molecular states, besides reproducing the masses of the existing pentaquark candidates, , , , , and . We can search for the non-strange , , , and pentaquark molecular states with isospin (or ) and with spin parity , , and , respectively in the decay,
and search for the strange , , , and pentaquark molecular states with isospin (or ) and spin parity , , , and , respectively, in the and decays,
Download figure:
Standard imageIV. CONCLUSION
In this study, we investigated the , , , , , , , and pentaquark molecular states with and without strangeness using the QCD sum rules in detail by consistently performing operator product expansion up to vacuum condensates of dimension 13. The modified energy scale formula was used to choose the best energy scales of the spectral densities at the quark-gluon levels and make predictions for the masses of the new pentaquark molecular states, besides reproducing the masses of existing pentaquark candidates , , , , and . The present calculations support assigning as the pentaquark molecular state with the quantum numbers and , assigning the is a pentaquark molecular state with quantum numbers and , as the pentaquark molecular state with quantum numbers and , and as the pentaquark molecular state with quantum numbers and . However, we cannot exclude the possibilities of assigning as the pentaquark molecular state with the quantum numbers and and assigning as the pentaquark molecular state with quantum numbers and due to the uncertainties. In the calculations, we observe that the predicted masses of the pentaquark molecular states without strangeness and with strangeness have a mass gap of approximately , which is consistent with the light-flavor, , breaking effects of the heavy baryons in flavor sextet . We can search for both old and new pentaquark molecular states in the decay of , , and in the future to conduct more robust investigations and shed light on the nature of and states.
Footnotes
- *
Supported by National Natural Science Foundation of China (11775079)