-
Quarks and gluons are the fundamental degrees of freedom of quantum chromodynamics (QCD). Because of the nonperturbative feature of QCD, we can only observe confined colorless hadrons. A normal hadron has two modes: a meson is made up of one quark and one antiquark, and a baryon is made up of three (anti)quarks. Multiquark hadrons made up of more than three quarks were proposed at the beginning of the construction of the quark model by Gell-Mann and Zweig [1–3]. However, the existence of tetraquarks and pentaquarks was not proven until the observation of
XYZ states [4], hidden-charmPc states [5–8], doubly-charmT+cc [9, 10] and fully-charm tetraquark states [11], etc. [12–18].Five
J/ψϕ structuresX(4140) ,X(4274) ,X(4500) ,X(4685) , andX(4700) in theB+→J/ψϕK+ decay process were observed by the LHCb Collaboration [19–21], CDF Collaboration [22, 23], CMS Collaboration [24], D0 Collaboration [25], and BARAR Collaboration [26].X(4140) andX(4274) are considered as thecsˉcˉs tetraquark ground states, whereasX(4500) andX(4700) are considered as thecsˉcˉs tetraquark excited states, in various theoretical methods [27–35]. In Ref. [36],X(4685) was also considered as the axial vector 2S radial excitedcsˉcˉs tetraquark state. In Refs. [37–40], the mass spectra of the S-wave and D-wavecsˉcˉs tetraquarks in different excitation structures are calculated using the QCD sum rules method. There have been many theoretical studies on the inner structure of these X's, such as the molecular states [41–59], compact or diquark-antidiquark states [27, 28, 60–68], cusp effects [69–71], dynamically generated resonances [72, 73], conventional charmonium [74], and hybrid charmonium states [42, 43]. However, overall, the inner structure ofX(4140) ,X(4274) ,X(4500) ,X(4685) , andX(4700) remains an open question.In the molecular picture, a
Xcsˉcˉs is formed by a strange-charmed mesonD+s (D−s ) and aD∗−s (D∗+s ), while aX(3872) is formed by a charmed mesonD0 (D∗0 ,D+ ,D− ) and a¯D0 (¯D∗0 ,D∗− ,D∗+ ). In the tetraquark picture, aXcsˉcˉs is formed by a spin triplet diquark[cs]1 (spin singlet diquark[cs]0 ) and a spin singlet antidiquark[ˉcˉs]0 (spin triplet antidiquark[ˉcˉs]1 ), while aX(3872) is formed by a diquark[cq]1 ([cq]0 ) and a[ˉcˉq]0 ([ˉcˉq]1 ), q for u and d quarks. Although light quarks u and d inX(3872) are replaced with s quarks inXcsˉcˉs , their inner structures may or may not be the same. This motivates the present study, in which we examine whether the approach we proposed in Ref. [75] can also be applied to theXcsˉcˉs case and thus find a way to distinguish the two internal structures with heavy ion measurements. In this work, we try to distinguish the two aforementioned possible inner structures ofXcsˉcˉs , i.e., a loose hadronic molecule or a compact tetraquark, by studying its production in heavy ion collisions. The remainder of this paper is organized as follows. In section II, we introduce the generation mechanism ofXcsˉcˉs into the AMPT model corresponding to its two possible inner structures following the production ofX(3872) described in Ref. [75]. In section III, we examine the production ofXcsˉcˉs as a function of centrality, transverse momentum, and rapidity. A volume effect is found, which can be a probe of the inner structure ofXcsˉcˉs . A summary and outlook are presented in section IV. -
In this study, we generate a total of one million minimum bias events for Pb-Pb collisions at
√sNN= 5.02 TeV by using the framework developed in Ref. [75]. We introduce the production mechanism to produceXcsˉcˉs for its two possible configurations, i.e., the hadronic molecular configurations and the tetraquark configurations into the default version (v1.26t9b) of the AMPT transport model [76]. Given thatXcsˉcˉs contains (anti-)charm quarks and (anti-)strange quarks, we need to generate a reasonable number of individual charm and strange quarks in the partonic phase. On top of the default version of AMPT, we modify the factor K ([77]) to enhance the initial c andˉc spectra because of a lack of some channels related to initial heavy quarks. The AMPT calculation gives a reasonable (order-of-magnitude) description of the experimental data [78] for the total yield ofD++D+∗ in the lowpT region (see upper panel of Fig. 1). For the strange quarks, an upper limit on the relative production of strange to non-strange quarks in AMPT is set to 0.6 because of the strangeness enhancement effect (see [79]), and our calculations also give a reasonable (order-of-magnitude) description of the experimental data [80] for the yield ofD+s mesons (see lower panel of Fig. 1). The main purpose of this work is to distinguish two inner structures ofXcsˉcˉs through their significantly different production rates. The difference in D andD+s meson production between our calculation and experimental data should not influence the relative yield between two inner structures and thus cannot change the qualitative results.Figure 1. (color online) Upper panel: total production of
D++D+∗ from the ALICE Collaboration [78]; lower panel: production ofD+s from the ALICE Collaboration [80]. The bands reflect the uncertainty due to the constituent composition, as discussed around Eq. (1), and are obtained by varying the composition fraction by ±10%.We use the same production mechanism developed in Ref. [75] for the hadronic molecule and tetraquark configurations of the
Xcsˉcˉs . For the molecular picture, the charmed-strange mesons are collected after the hadronization process. Then,D+s (D−s ) andD−∗s (D+∗s ) are coalesced (similar to the hadronization process mentioned in [76]) to form the "molecule"Xcsˉcˉs according to the following conditions: the relative distance within the region[5fm,7fm] and invariant mass within the region[2MD+s,2MD+∗s] . For the tetraquark picture, the "tetra"Xcsˉcˉs is formed via two steps. (i) First, diquarks (cs ) and diquarks (ˉcˉs ) are formed by matching a (anti-)charm quark with the nearest (in both position space and momentum space) (anti-)strange quark in the parton. (ii) Then, these (anti)diquarks are coalesced to form theXcsˉcˉs according to the following conditions: the relative distance<1 fm and invariant mass within the region[2M[cs]1,2M[cs]0] (the spin triplet and singlet diquark masses are defined in Refs. [30, 31]). Owing to a lack of spin information in the AMPT model for the formation of the charmed-strange mesons and (anti)diquarks, the relative yield ratios are estimated using the thermal model:R(AB)≡Yield(A)Yield(B)=e−(mA−mB)/Tfreezeout,
(1) where
mA andmB represent the masses of hadrons A and B, respectively. Here,Tfreezeout≃160MeV is the freeze-out temperature. For the hadronic picture, A and B are theD+s andD+∗s mesons, respectively. For the tetraquark picture, A and B are the spin triplet and singlet diquark, respectively. This estimate indicates a composition of30% (70% ) forD+s (D+∗s ) and a composition of35% (65% ) for spin triplet(singlet) diquarks. We also vary the composition between20% (80% ) and40% (60% ) to show the uncertainty bands. -
Within this simulation framework, we use the Monte Carlo method to generate a total of one million minimum bias events for Pb-Pb collisions at
√sNN=5.02TeV . The inclusive yield ofXcsˉcˉs is found to be approximately 42000 in the molecular picture and approximately 200 in the tetraquark picture. As a benchmark for comparison, we also estimate the yield ofX(3872) within the same framework (see the production mechanism in Ref. [75]; the yield should be multiplied a factor of1/4 owing to wavefunction normalization for both the molecular and tetraquark pictures). The inclusive yield ofX(3872) is found to be approximately 171000 in the molecular picture and approximately 600 in the tetraquark picture. The yield ofXcsˉcˉs is approximately1/4 of that ofX(3872) . Compared with the experimental data ofX(3872) measured by the CMS collaboration for PbPb collisions at√sNN=5.02TeV [81], our finding suggests that an observable signal ofXcsˉcˉs can be measured in heavy ion collisions at the LHC energy.One can also find that the production in the molecular picture significantly exceeds that in the tetraquark picture, by a factor of 200 — a 2-order-of-magnitude difference. This result may be understood as follows: the
c−ˉc ands−ˉs quarks must be pair produced in the initial conditions of heavy ion collisions and then expand and cool with the bulk flow; the moleculeXcsˉcˉs needs a large volume to be formed, while the tetraquarkXcsˉcˉs needs a compact volume to be formed; thus, the probability of the formation of hadron molecules is far higher than that for the tetraquark state.We plot the
Xcsˉcˉs production as a function of centrality in Pb-Pb collisions at√sNN=5.02TeV for the hadronic molecular state and tetraquark state in Fig. 2. One can find that the yield of theXcsˉcˉs in the molecular picture is 2 orders of magnitude larger than that in the tetraquark picture. From the central collision region to the peripheral collision region, the production first increases then decreases for both the molecular state and the tetraquark state, and the slope of the decrease is far larger in the molecular state than in the tetraquark state. This results from a competing effect between the volume of the bulk system and the size ofXcsˉcˉs . For central collisions, the number of (anti-)charm and (anti-)strange quarks is large, the bulk volume is large, and its evolution time is long; thus, the (anti-)charm and (anti-)strange quarks separate sufficiently, which benefits the production of a large-size molecular state. For the peripheral collisions, both the number of (anti-)charm and (anti-)strange quarks and the size of the fireball are small; as such, the evolution time of the fireball is short, which benefits the production of small-sized tetraquark states. This size effect could help to explore the internal structure ofXcsˉcˉs through different collision systems, e.g., Pb-Pb, Au-Au, Xe-Xe, Cu-Cu, O-O, and d-A/p-A.Figure 2. (color online) Centrality dependence of the
Xcsˉcˉs in Pb-Pb collisions at√sNN=5.02TeV for the hadronic molecular configuration (red solid boxes) and tetraquark configuration (blue shaded boxes). The bands reflect both statistical uncertainty from our simulations and the uncertainty due to the constituent composition, as discussed around Eq. (1), and are obtained by varying the composition fraction by ±10%.In Fig. 3, we present the rapidity and the transverse momentum distributions of
Xcsˉcˉs . One can find that the distribution for both the hadronic molecular state and the tetraquark state is similar to that of the usual hadrons [82, 83]. We also show the elliptic flow coefficientv2 ofXcsˉcˉs as a function of the transverse momentumpT in Fig. 4. The elliptic flow is sensitive to the geometry of the initial fireball and the generation mechanism ofXcsˉcˉs .Figure 3. (color online) Rapidity y and transverse momentum
pT distribution of theXcsˉcˉs yield in Pb-Pb collisions at√sNN=5.02TeV for the hadronic molecular configuration (red solid boxes) and the tetraquark configuration (blue shaded boxes). The bands are determined as described in Fig. 2.Figure 4. (color online) Elliptic flow coefficient
v2 versus transverse momentumpT for the producedXcsˉcˉs in minimum bias Pb-Pb collisions at√sNN=5.02TeV , predicted from our computation for the hadronic molecule picture. The bands are determined as described in Fig. 2. -
In this work, we studied the yields of
Xcsˉcˉs for Pb-Pb collisions at√sNN=5.02TeV by introducing the production mechanism of two possible configurations, i.e., the hadronic molecular state and tetraquark state, into the AMPT model. We found that the production in the molecular picture exceeds that in the tetraquark picture by two orders of magnitude. The centrality distribution of the yields ofXcsˉcˉs shows a strongly decreasing trend for the hadronic molecular state and a mild change for the tetraquark state. This system size dependence could be a good probe for the inner structure ofXcsˉcˉs . We also showed the rapidity and the transverse momentum distributions ofXcsˉcˉs production, as well as its elliptic flow coefficient, as a function of the transverse momentum, which can be tested in future experimental measurements. In Ref. [80], a strangeness enhancement effect in heavy ion collisions was found by ALICE collaboration, which could be evidence for quark-gluon plasma. We expect a similar effect to be found in the ratio ofXcsˉcˉs toX(3872) , which will be studied in our future work. -
The authors would like to thank Dr. J. Liao, E. Wang, Q. Wang, and H. Xing for the helpful discussion.
Production of Xcsˉcˉs in heavy ion collisions
- Received Date: 2023-01-18
- Available Online: 2023-05-15
Abstract: The yields of