cpcChinese Physics C1674-1137Chinese 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 Ltdcpc_43_2_02410110.1088/1674-1137/43/2/02410143/2/024101PaperNuclear physicsDifferent production sources of light nuclei inultra-relativistic heavy-ion collisions *
Supported by National Natural Science Foundation of China (11505104, 11575100, 11675091)
We systematically study different production sources of light nuclei in ultra-relativistic heavy-ion collisions with a new method, an exclusive quark combination model + an inclusive hadron recombination model. We take deuterons and ^{3}He produced in Pb-Pb collisions at
TeV as examples to show the contribution of different production sources by studying their rapidity densities
, yield ratios and transverse momentum (
) spectra just after hadronization and at the final kinetic freeze-out. We find that about a half of
and a fourth of ^{3}He created just after hadronization can survive after the hadronic evolution process. Nucleons from
resonance decays make a much larger contribution to the regeneration of light nuclei at the hadronic phase stage, and this contribution is about 77% and 90% for
and ^{3}He, respectively, calculated at the final kinetic freeze-out. In addition, we give an explanation for the constant behaviors of yield ratios
and ^{3}He
as a function of the averaged charged multiplicity in Pb-Pb collisions and also provide a possible explanation for the observation that
in Pb-Pb collisions is larger by a factor of about two than in pp collisions at LHC energies.
The production of light nuclei in ultra-relativistic heavy-ion collisions is of importance for many topics in nuclear and particle physics. It does not only help to understand the mechanism of cluster formation in the interior of the fireball in a heavy-ion collision, but can also serve as an effective probe of the fireball freeze-out properties [1-8]. Experimental measurements of light nuclei have been extensively performed at the Relativistic Heavy Ion Collider (RHIC) [9-13], at relatively low-energy collisions by the NA49 Collaboration at the Super Proton Synchrotron (SPS) [14-18], and more recently in very high energy reactions at the Large Hadron Collider (LHC) [19-22]. An interesting phenomenon observed at the LHC is that the yield ratio
in Pb-Pb collisions is larger by a factor of about two than in pp collisions, while the ratios of hadrons such as
and
, etc., do not show significant differences between Pb-Pb and pp collisions [19, 23]. Up to now, there is no satisfactory explanations for such phenomena.
There are two popular production mechanisms used to describe light nucleus formation. One is the thermal model [24-26], and the other is the recombination/coalescence model [27-35]. The mechanism of the recombination/coalescence production for light nuclei was studied already in 1960s and is now rather well understood [27-35]. In the recombination/coalescence models, the light nuclei can be formed by coalescence of nucleons produced just after hadronization and/or those from resonance decays. Since the binding energies of light nuclei are very small (~a few MeV), final-state coalescence, i.e., nucleons that recombining into light nuclei at the final stage of the hadronic phase evolution (at the final kinetic freeze-out), is commonly adopted in different recombination/coalescence models [32-35]. In fact, light nuclei can be formed during the entire dynamical process of hadronic phase evolution. In the beginning of hadronic evolution, when the temperature is high, light nuclei may be formed and then destroyed immediately due to disruptive collisions in the environment. As the temperature decreases, they may be created again during final-state coalescence. An important question is whether light nuclei are mainly produced just after hadronization, i.e., in a relatively high temperature environment, or at a later stage of hadronic evolution via final-state coalescence [36].
In this paper, we make an estimate of how many light nuclei can be formed just after hadronization, and with what probability the formed light nuclei can coming from until the final kinetic freeze-out. Also, we determine the proportion of the finally observed light nuclei that survive the early formation, and those that are regenerated in the hadronic evolution stage. We employ a hadron recombination model to make a systematic study of the production of light nuclei just after hadronization and at the final kinetic freeze-out, based on an exclusive description of the directly-produced and final-state hadrons using the quark combination model SDQCM developed by the ShanDong group [37-39].
The paper is organized as follows. In Section 2, we introduce the hadron recombination model which describes the light nucleus formation. In Section 3, we calculate the rapidity densities
, yield ratios and
spectra of
and ^{3}He just after hadronization and at the final kinetic freeze-out in Pb-Pb collisions at
2.76 TeV. We present the results and discussions in this section. Section 4 summarizes our findings.
The hadron recombination model
After hadronization of the partonic system produced early in high energy heavy-ion collisions, various kinds of hadrons are produced, and light nuclei can then be formed through the recombination of protons and neutrons. The hadron recombination model deals with how nucleons recombine into different light nuclei. In this model, the momentum distributions of
and ^{3}He can be expressed as follows
is the two-nucleon joint coordinate-momentum distribution for the
cluster, and
is the three-nucleon joint coordinate-momentum distribution for the
cluster. Kernel functions
and
denote the probability densities for
and
, with momenta
and
at
and
, to recombine into a
of momentum
, and for
,
and
, with momenta
,
and
at
,
and
, to recombine into a ^{3}He of momentum
, respectively. Eqs. (1) and (2) are the starting point for describing the production of light nuclei in high energy reactions based on the basic recombination/coalescence idea. All coordinate and momentum variables in Eqs. (1) and (2) are three-dimensional and the integration refers to the whole coordinate and momentum space.
Joint distributions
and
are the number densities that satisfy
and
Here,
is the number of all possible
clusters and
is the number of all the possible
clusters in the bulk hadronic system that we consider.
and
are the number of protons and of neutrons in the corresponding hadronic system. We rewrite
and
, so that the coordinate-momentum distributions are normalized to unity, which we denote by the superscript (
). In terms of these normalized joint distributions, we have
Kernel functions
and
carry the kinematical and dynamical information of the nucleons recombining into light nuclei, but their precise forms are ambiguous due to their complicated non-perturbative properties. Despite this, we know that the kernel functions should have the following three properties. The first is that they must satisfy momentum conservation, so they should contain the term
or
. The energy conservation in the classical recombination process of two/three on-shell nucleons is broken. However, considering the nucleon interaction dynamics, the nucleons entering light nuclei are intermediate scattering states which are not on the mass-shell, so that the exact energy-momentum conservation can also be satisfied. In practice, since the momentum interval between two/three nucleons in recombination is usually very small
0.2 GeV/c, the relative magnitude to which energy conservation is broken
is only a few percent even for on-shell nucleon recombination. The second property is that the kernel functions should depend on the intrinsic quantum numbers of the formed light nuclei, which are denoted by
and
. Considering the simple spin counting, we have
and
. The last and the most important property is the coordinate and momentum dependence of the kernel functions. Even though they can not be solved explicitly from first principles, we know that they should increase as the relative space positions and relative momenta of the recombined nucleons decrease. In general, the coordinate and momentum dependence of the kernel functions may be coupled to each other. In this paper, we do not study such coupling effects and in the following, we consider only the simple case where the coordinate and momentum dependence of the kernel functions are decoupled from each other, i.e. they are factorized. Based on the above three properties, we can write the kernel functions as follows.
The coordinate and momentum parts are given from the Wigner function method as Gaussians
,
,
, and
. The superscript ‘
’ in the coordinate and momentum variables denotes the coordinates and momenta of the nucleons in the rest frame of the forming light nuclei. The widths of the Gaussian distributions are:
fm,
GeV,
fm, and
GeV [34].
and
are the root-mean-square radii for
and ^{3}He, and are
fm and
fm [40], respectively.
We further assume that the normalized joint distributions of the nucleons are coordinate and momentum factorized as follows
Substituting Eqs. (5-8) into Eqs. (3) and (4), we have
We denote the coordinate integrals in Eqs. (9) and (10) as
and
to get
where
From Eqs. (13) and (14), one can easily see that the
coefficients originate from the overlap integral between the nuclear spatial wave function and the spatial distribution of nucleons which form the light nuclei.
We change the two integral variables in Eq. (13) to
and
and those in Eq. (14) to
,
, and
. Then we can write
and the normalization constraints
To evaluate
and
, we assume that the coordinate variables are decoupled in the coordinate distributions as
The normalized distributions of the relative space positions of the nucleons are
,
and
[41].
is the Gaussian width of the relative space positions of the nucleons, and is about 6.1 fm just after hadronization, and 7.0 fm at the final kinetic freeze-out [41, 42]. Note that
used here is the effective radius of the fireball as it evolves till the freeze-out of light nuclei in the static picture, and is not equivalent to the geometrical radius of the fireball [43].
is linked to the effective correlation length of nucleons in the fireball in the static picture. That is why we adopt the value of this radius as obtained from experiential measurements of the femtoscopic radii [42], rather than that extracted from hydrodynamic or transport models at the final kinetic freeze-out [43, 44]. Substituting Eqs. (19) and (20) into Eqs. (15) and (16), we have
Assuming the instantaneous recombination in the rest frame of the formed
or
cluster, we can easily get the numerical results for
and
.
The production of light nuclei just after hadronization and at the final kinetic freeze-out can be computed using Eqs. (11) and (12), as the joint-momentum distributions
and
are known. Ignoring the nucleon correlations, we have
and
, where the single nucleon momentum distributions
and
just after hadronization can be easily computed by SDQCM. When resonance decays and hadronic scattering effects are included, we get
and
at the final kinetic freeze-out.
As a brief summary of this section, we would like to emphasize that with the exclusive quark combination model SDQCM, used to obtain momentum distributions of different light, strange and heavy flavor hadrons, we can apply the hadron recombination model to investigate the production of light nuclei, hypernuclei and other molecular states. There are no extra free parameters in this method apart from the assumptions and/or approximations explicitly stated above, such as the coordinate and momentum factorization of the kernel and joint distribution functions, and decoupling of different coordinate variables in the coordinate distributions. This newly proposed method, an exclusive quark combination model + an inclusive hadron recombination model, is an attempt to mimic the quantum mechanical process of recombination/coalescence. It is very convenient for studying different production sources of light nuclei, hypernuclei and other molecular states.
Results and discussion
In this section, we present a systematic study of the production of light nuclei at midrapidity in central Pb-Pb collisions at
TeV. First, we employ SDQCM to compute rapidity density
and
distribution of various hadrons. SDQCM is unique for its ability to provide yields, rapidity distributions and
spectra of different hadrons, and has been successfully applied in studies of
annihilation, p-p reactions and relativistic heavy-ion collisions [23, 37-39, 45-47]. The starting point of SDQCM is a system consisting of different flavors of constituent quarks and antiquarks, whose phase space distributions are given as input. The quarks and antiquarks combine into different species of directly-produced hadrons according to a sophisticated combination rule [37-39]. After hadronization, directly-produced hadrons undergo a hadronic evolution stage during which short lifetime resonances decay and scattering among different hadrons occurs. Hadronic evolution involves many effects and is very complicated. We do not investigate all effects here, but only focus on those that influence hadron yields and shapes of hadron
spectra. For hadron yields, we only consider resonance decay contributions and neglect scattering effects, which are comparatively small [48, 49]. We directly use the UrQMD results at the LHC energies for hadrons before and after hadronic evolution [50, 51] to correct the shape of hadron
distributions. We calculate the production of relatively stable hadrons, such as
,
,
mesons and
,
,
,
baryons etc., so that the above approximation for handling the effects of the hadronic evolution is sufficient.
To compute momentum distributions of directly-produced hadrons with SDQCM, we need the momentum distributions of constituent quarks just before hadronization. As we focus on the midrapidity region at the LHC energy, the influence of net quarks from incident nuclei is ignored. Also, the isospin symmetry for
and
quarks is assumed. The input values for the rapidity densities of
,
and
quarks are
and
, respectively. The longitudinal rapidity distributions for
,
and
quarks are chosen to be uniform since we focus on hadron production in the rapidity window
. For the transverse momentum distributions of light and strange quarks, we use an exponential parametrization
as we focus on thermal hadron production. The slope parameter
for
,
quarks is 0.27 GeV, and 0.33 GeV for the
quark. These values are slightly smaller than those obtained in our previous work [52], where the production of various light, strange and charm hadrons was studied and the corresponding constituent quark distributions included in the collective radial flows generated during the partonic evolution stage and the hadronic evolution stage. In this study, constituent quark distributions carry only the radial flows created during partonic evolution. The radial flows generated during the hadronic evolution stage are added to hadrons produced after hadronization.
With the above inputs, the rapidity densities
and
distributions for directly-produced hadrons, i.e., those hadrons produced just after hadronization, can be calculated with SDQCM. Including very short lifetime
,
,
and
decays and considering hadronic scattering, we obtain
and
distributions of hadrons at the final kinetic freeze-out. We want to point out that some processes, for example strong decays of some hadrons such as
, electromagnetic decays and weak decays, are not included because their lifetimes are relatively long and they decay after the kinetic freeze-out. After considering all strong and electromagnetic decays, we have
and
distributions for the total final state hadrons, which can be compared to the LHC experimental data. For more information on how SDQCM models hadronization, see Refs. [38, 39, 52, 53].
Rapidity densities
of hadrons in central Pb-Pb collisions at
TeV are listed in Table 1. The experimental data in the second column are from Refs. [54-57]. The third column presents the results for directly-produced hadrons. The fourth column gives the results for hadrons at the final kinetic freeze-out. The last column are the results for the total final state hadrons, which agree with the experimental data, and are comparable to our previous results calculated by analytical formulas [52] and the results given in [58]. The results for directly produced
and
, those created at the final kinetic freeze-out and the total final state ones are the same, since there are no other decay contributions. The directly-produced
ratio is about 0.4, which reveals the strangeness in heavy-ion collisions and is comparable to our previous studies [39]. Directly-produced protons are about half of the total final state protons, while the other half come from
decays. The lifetime of
is so short that they decay before the final kinetic freeze-out occurs. As a result, protons at kinetic freeze-out are identical to the total final state protons. This is the reason why some coalescence models use the momentum distribution of protons measured experimentally as input when studying the production of light nuclei [34].
baryons are directly produced or come from
and
decays, so that there are more of them in the total final state than after the kinetic freeze-out, and both are more numerous than directly-produced ones.
can decay to
, and the decay takes place after kinetic freeze-out. The directly-produced
and those after the kinetic freeze-out are the same, and both are less in number than the total final state ones. For more information of hadron decays, see Ref. [59].
Rapidity densities dN/dy of identified hadrons in central Pb-Pb collisions at
TeV. The experimental data are from Refs. [54-57]. The last three columns are the results calculated with SDQCM for DIrectly-produced hadrons, hadrons at KInetic freeze-out and Total Final State hadrons.
hadron
experimental data
SDQCM
DI
KI
TFS
55
501
737
55
501
737
22
88
111
22
88
111
0
0
106
13.2
13.2
13.2
16
33
33
16
33
33
7.8
17
25
−
7.8
17
25
2.60
2.60
3.90
2.60
2.60
3.90
0.52
0.52
0.52
0.52
0.52
0.52
distributions for hadrons are also computed with SDQCM and the results are shown in Fig. 1. The filled symbols are the experimental data from Refs. [54-57]. Dashed lines are the results for directly-produced hadrons just after hadronization. Solid lines and dot-dashed lines are for the total final state hadrons, including strong and electromagnetic decays in cases when hadronic interactions are turned on and off, respectively. The three lines for
baryons are almost the same as there are no resonance decay contributions to
and hadronic interactions of multistrange hadrons are very small. The same holds for
mesons. Distributions for directly-produced
,
,
,
and
(dashed lines) are harder than those for the final state with decay contributions only (dot-dashed lines), especially for
and
mesons. Hadronic interactions after hadronization push
spectra of hadrons, especially for
and
, to higher values since more collective radial flows are created at this stage. The effects of hadronic interactions on
distribution are very complicated since many mesons and baryons decay to
, and we do not plot it in Fig.1. Solid lines for total final state hadrons
,
,
,
,
,
and
agree very well with the experimental data. A good reproduction of the measurements for the total final state, kinetic freeze-out and directly-produced protons is necessary for the study of light nuclei.
(color online) p_{T} distributions of hadrons at midrapidity in central Pb-Pb collisions at
TeV. The filled symbols are the experimental data from Refs. [54-57]. Dashed lines are the results for directly-produced hadrons just after hadronization. Solid lines and dot-dashed lines are for the total final state hadrons including strong and electromagnetic decays in cases when the hadronic interactions are turned on and off.
In view of the good performance of SDQCM for hadrons^{①)}
The detailed calculations of hadrons produced in central Pb-Pb collisions at TeV with SDQCM can be found in our previous work [52].
, we turn to the study of light nuclei at midrapidity in central Pb-Pb collisions at
TeV. We perform the calculations using the hadron recombination model introduced in Sec. 2, for the stages just after hadronization and at the final kinetic freeze-out. We take into account the destructive and regenerative collisions during the production of light nuclei in the intervening time based on an effective method, rather than using the conventional transport and/or kinetic models such as in Refs. [60, 61]. We compute the production of light nuclei just after hadronization and before hadronic evolution using the distributions of directly-produced nucleons. These light nuclei are formed immediately and are denoted by JAH (Just After Hadronization). The JAH light nuclei undergo hadronic scatterings with hadrons. The influence of hadronic scattering on JAH light nuclei are taken into account by modifying the coordinates and momenta of the constituent-nucleons affected by the scattering. After hadronic evolution, i.e., at the final kinetic freeze-out, we check to see if the constituent-nucleons still meet the requirements of binding together into light nuclei. If that is the case, the JAH light nuclei are kept, otherwise they are discarded and the corresponding constituent-nucleons are promoted as real nucleons.
Recalling Eqs. (11) and (12), we calculate rapidity densities
and
spectra of
and ^{3}He just after hadronization as well as at the final kinetic freeze-out as we know
and
. We count the initially produced light nuclei just after hadronization by evaluating the probability for recombination of every possible directly-produced
pair or
cluster to recombine into a
or ^{3}He. If an initially produced light nucleus is destroyed during hadronic evolution, the constituent nucleons are separated in phase space. Otherwise, they are still close enough. Taking into account this feature, we re-calculate the number of light nuclei at the final kinetic freeze-out using the evolved distributions of these directly-produced nucleons (after accounting for the volume expansion of the hadronic system and the flow generation during the hadronic evolution stage). We treat these light nuclei as a part of the initially produced light nuclei. We also calculate the number of light nuclei at the final kinetic freeze-out formed from nucleons from
resonance decays and denote these light nuclei as regenerated during the hadronic evolution stage. The input nucleon momentum distributions just after hadronization are
and
, where the superscript
denotes directly-produced nucleons. The input momentum distributions at the final kinetic freeze-out are as follows.
The superscript
denotes nucleons at the final kinetic freeze-out, which are composed of two parts. One is from the directly-produced nucleons denoted by the superscript
, and the other is from
resonance decays denoted by the superscript
, i.e.,
Note that
and
are momentum distributions for nucleons produced just after hadronization that have survived till kinetic freeze-out. They carry the radial flows generated during the hadronic evolution stage besides those from the partonic stage. The first terms in Eqs. (25) and (26) contribute to the surviving light nuclei, while all the other terms contribute to regenerated light nuclei.
With the momentum distributions of directly-produced protons and those at final kinetic freeze-out (i.e., total final state ones) shown in Fig. 1, we can compute rapidity densities
and
spectra of
and ^{3}He. Momentum distributions of neutrons are the same as those of protons due to isospin symmetry. The third and fourth columns in Table 2 are the
results for just after hadronization and at the final kinetic freeze-out, respectively. The experimental data in the second column are from Ref. [19]. From Table 2, one can see that the results computed for just after hadronization are much lower than the experimental data. The result for
calculated at the final kinetic freeze-out agree well with the data, while that of ^{3}He is smaller than the lower limit of the data. The reason for underestimating of ^{3}He is that we only include the dominant production contribution for ^{3}He, while the other contributions such as
He and decay contributions from
H are not included. The number of surviving
and ^{3}He formed just after hadronization and including the regenerated ones at the final kinetic freeze-out are shown in the last two columns inTable 2. About half of
produced just after hadronization survive hadronic evolution, but only about one fourth of ^{3}He do so. This is because we treat
as a cluster of two quasi-free constituent nucleons and ^{3}He as a cluster of three quasi-free constituent nucleons when considering hadronic scattering effects. The dissociation probability for
is proportional to 2 times the nucleon-hadron rescattering cross section, while that for ^{3}He is proportional to three times the same cross section. It is therefore much easier for ^{3}He to be destroyed compared to
. The fraction of disassociated JAH deuterons is about 1/2, and that of JAH ^{3}He is about 3/4, instead of 1. This is because the rate of inter-particle interactions/scatterings in the hadronic phase decreases quickly as the hadronic system expands. From this point, it ismeaningful to discuss JAH light nuclei production. For the total light nuclei at the final kinetic freeze-out, about 77% of
and 90% of ^{3}He are regenerated. These numerical results explicitly illustrate the different production sources of
and ^{3}He.
Rapidity densities dN/dy of d and ^{3}He in central Pb-Pb collisions at
TeV. The experimental data are from Ref. [19]. The third and fourth columns are the results for just after hadronization and at the final kinetic freeze-out. The surviving light nuclei and the regenerated nuclei at final kinetic freeze-out are given in the last two columns.
experimental data
just after hadronization
kinetic
surviving at kinetic
regenerated at kinetic
d
(9.82±0.04±1.58)×10^{−2}
4.18×10^{−2}
9.05×10^{−2}
2.11×10^{−2}
6.94×10^{−2}
^{3}He
(2.76±0.09±0.62)×10^{−4}
0.68×10^{−4}
1.62×10^{−4}
0.17×10^{−4}
1.45×10^{−4}
A series of interesting observables are yield ratios for light nuclei, such as
, ^{3}He
and ^{3}He
. Recalling Eqs. (11) and (12), we have
To analyze the behavior of the yield ratios, analytical expressions for
and
are required. We consider the Lorentz contraction effect, and replace
and
with
and
, respectively, where
is the effective Lorentz contraction factor related to the coordinate transformation from the laboratory frame to the center-of mass frame. In this case, we have
One can easily see from Eqs. (30) and (31) that the
factors for
and ^{3}He decrease with increasing fireball radius at freeze-out. This means that it becomes more difficult for a certain
or
cluster to form a light nucleus when the freeze-out fireball becomes larger. Considering that
is about 6~7 fm in heavy-ion collisions at the LHC energies [41, 42], and that
fm and
fm, both
and
are much larger than 1. As
is slightly larger than 1, we get that
and
. Momentum integrals in Eqs. (27-29) are very weakly dependent on the collision centrality as the normalized nucleon momentum distributions do not show significant differences in different centrality collisions.
and
depend only on the intrinsic quantum numbers of light nuclei and are independent of the collision energy and centrality. Therefore, one can see that the ratios in Eqs. (27-29) are approximately proportional to nucleon densities instead of nucleon numbers, and this can well explain the constant behavior of
and ^{3}He
as function of the average charge multiplicity observed by the ALICE Collaboration [19]. The calculated values of
, ^{3}He
and ^{3}He
for just after hadronization and at the final kinetic freeze-out are given in the third and fourth columns in Table 3. The experimental data in the second column in Table 3 are from Ref. [19].
Yield ratios d/p, ^{3}He/p and ^{3}He/d in central Pb-Pb collisions at
TeV. The experimental data in the second column are from Ref. [19]. The third and fourth columns are the results for just after hadronization and at the final kinetic freeze-out. The fifth column is the result at the final kinetic freeze-out when hadronic interactions are turned off. The last column are the experimental data for pp collisions at
TeV.
ratio
data
just after hadronization
kinetic
kinetic (noHI)
data for pp 7 TeV
d/p
(3.2188±0.0142±0.6067)×10^{−3}
2.61×10^{−3}
2.74×10^{−3}
1.31×10^{−3}
(1.63±0.20)×10^{−3}
^{3}He/p
(1.0611±0.0359±0.2570)×10^{−5}
0.43×10^{−5}
0.49×10^{−5}
0.22×10^{−5}
−
^{3}He/d
(2.81±0.86)×10^{−3}
1.63×10^{−3}
1.79×10^{−3}
1.63×10^{−3}
−
In order to better understand the effects of hadronic interactions during the hadronic evolution stage on the production of light nuclei, we also calculate these yield ratios at the final kinetic freeze-out in the case when hadronic interactions are turned off after hadronization. In that case, there is no formation nor destruction of light nuclei during the hadronic system evolution, and the number of light nuclei is just that after hadronization. The number of protons includes those from resonance decays, apart the directly-produced ones, i.e.,
The numerical results for this case are given in the fifth column in Table 3, which shows that d/p,^{3}He/p and ^{3}He/d ratios decrease by about 48%, 45% and 91% , respectively, when hadronic interactions are turned off compared to those including hadronic interactions. This decrease can be also seen in pp collisions, where hadronic interactions after hadronization disappear. In column six of Table 3, the experimental result for d/p in pp collisions at
TeV [19] is given, and is seen to be comparable with the numerical result for d/pin Pb-Pb collisions without hadronic interaction effects. Additional measurements of ^{3}He/p and ^{3}He/d in pp collisions could help to confirm the effect of negligable hadronic interactions in the hadronic system evolution.
distributions for
and ^{3}He are shown in Fig. 2. The open circles and diamonds are the experimental data for
and ^{3}He, respectively, measured by the ALICE Collaboration in Ref. [19]. In Fig. 2(a), the calculated results for just after hadronization are shown by a dashed line, and those at the final kinetic freeze-out by a solid line. For both
and ^{3}He,
spectra calculated just after hadronization are softer than at the final kinetic freeze-out. This is because the
distribution of directly-produced protons is softer than that at the kinetic freeze-out, as it acquires more radial flows during hadronic evolution. In Fig. 2(b),
distributions of the surviving and regenerated
and ^{3}He at the final kinetic freeze-out are represented with a dotted line and dot-dashed line, respectively. The solid lines denoting the total
and ^{3}He at the kinetic freeze-out are shown for guidence.
(color online) p_{T} distributions of d and ^{3}He at midrapidity in central Pb-Pb collisions at
TeV. The open circles and diamonds are the experimental data for d and ^{3}He from [19]. (a) The solid lines are the results at the final kinetic freeze-out, and the dashed lines are for just after hadronization. (b) The solid, dotted and dash-dotted lines are the results for the total, surviving and regenerated light nuclei at the final kinetic freeze-out.
Summary
We have proposed in this paper a new method, an exclusive quark combination model + an inclusive hadron recombination model, for a systematic study of different production sources of light nuclei in ultra-relativistic heavy-ion collisions. Using the exclusive quark combination model SDQCM, we obtained momentum distributions for directly-produced hadrons and those at the final kinetic freeze-out, including the hadronic scattering and resonance decays. These momentum distributions were taken as input for the inclusive hadron recombination model to calculate the production of light nuclei. In the hadron recombination model, we analytically deduced the formulas for formation of d and ^{3}He, and explicitly stated the assumptions used, such as coordinate and momentum factorization for kernel functions and nucleon joint distributions, and the decoupling of different coordinate variables in the coordinate distributions. We computed
and
spectra of
and ^{3}He as well as their ratios just after hadronization and at the final kinetic freeze-out in Pb-Pb collisions at
TeV. Light nuclei produced by the directly-produced nucleon combination and those combined from nucleons from
resonance decays (or those obtained by recombination of nucleons from decays with directly-produced ones) are clearly distinguished in our method.
The results of our calculations showed that about a half of d and a fourth of ^{3}He created just after hadronization cansurvive after the hadronic evolution process. Nucleons from
resonance decays contribute considerably to the regeneration of light nuclei, and their proportion at the final kinetic freeze-out is about 77% and 90% for
and ^{3}He, respectively. Ratios such as
, ^{3}He
and ^{3}He
were found to be effective observables to probe different production sources of light nuclei. Values of
, ^{3}He
and ^{3}He
just after hadronization are all smaller than those at the final kinetic freeze-out.
, ^{3}He
and ^{3}He
were also found to be effective probes of hadronic interaction effects on light nucleus formation after hadronization. Values of
^{3}He
and ^{3}He
when the hadronic interactions are turned off reduce to about 48%, 45% and 91% of those with hadronic interactions. The constant behavior of
and ^{3}He
as function of the average charge multiplicity in Pb-Pb collisions was explained. A possible explanation of the observed difference in
data in Pb-Pb collisions and in pp collisions at the LHC was also provided. All these results open new insights for understanding of the production of light nuclei in ultra-relativistic heavy-ion collisions.
The authors thank Qiang Zhao, Lie-Wen Chen and Kai-Jia Sun for helpful discussions.
ReferencesGossetJ.GutbrodH. H.MeyerG.PoskanzerA. M.SandovalA.StockR.Westfa11G. D.MekjianA. Z.KapustaJoseph I.NagleJ. L.KumarB. S.KusnezovD.SorgeH.MattielloR.ChenL. W.KoC. M.LiB. A.DanielewiczP.SchuckP.ScheiblR.HeinzU.SchmidtM.RöpkeG.SchulzH.AdlerC.et al(STAR Collaboration)AfanasievS.et al(PHENIX Collaboration)AbelevB. I.et al(STAR Collaboration)arXiv: 0909.0566 [nucl-ex]AdamczykL.et al(STAR Collaboration)YuNing(for the STAR Collaboration)arXiv: 1704.04335 [nuclex]BeardenI. G.et alAfanasevS. V.et al(NA49 Collaboration)AnticicT.et al(NA49 Collaboration)AnticicT.et al(NA49 Collaboration)AnticicT.et al(NA49 Collaboration)AdamJ.et al(ALICE Collaboration)AnielskiJ.(for the ALICE Collaboration)LeaRamona(for the ALICE Collaboration)TrogoloStefano(for the ALICE Collaboration)ShaoF. L.WangG. J.WangR. Q.LiH. H.SongJ.CsernaiL. P.KapustaJ. I.AndronicA.Braun-MunzingerP.StachelJ.StoeckerH.ChatterjeeS.MohantyB.ButlerS. T.PearsonC. A.SchwarzschildA.ZupancicC.SatoH.YazakiK.MrowczynskiS.DoverC. B.HeinzU.SchnedermannE.XueL.MaY. G.ChenJ. H.ZhangS.ChenG.ChenH.WuJ.LiD. S.WangM. J.SunK. J.ChenL. W.YinX. J.KoC. M.SunY. F.ZhuL. L.BufalinoS.(for the ALICE Collaboration)XieQ. B.LiuX. M.ShaoF. L.XieQ. B.WangQ.ShaoC. E.SongJ.ShaoF. L.XieQ. B.RopkeG.MrowczynskiS.AdamJ.et al(ALICE Collaboration)LinZ. W.KoC. M.PalSubrataKolbPeter F.RappRalfGouX. R.ShaoF. L.WangR. Q.LiH. H.SongJ.WangR. Q.ShaoF. L.SongJ.XieQ. B.LiangZ. T.SongJ.LiH. H.WangR. Q.ShaoF. L.BassS. A.DumitruA.NonakaC.BassS. A.SteinheimerJ.BleicherM.KnospeA. G.MarkertC.WernerK.SteinheimerJ.BleicherM.WangR. Q.SongJ.ShaoF. LWangR. Q.SongJ.SunK. J.ChenL. W.LiG.ShaoF. LAbelevB.et al(ALICE Collaboration)AbelevB.et al(ALICE Collaboration)AbelevB.et al(ALICE Collaboration)AbelevB.et al(ALICE Collaboration)arXiv: 1404.0495v2 [nucl-ex]AndronicABraun-MunzingerPRedlichKStachelJPatrignaniC.et alParticle Data GroupOhY.LinZ. W.KoC. M.ChoS.SongT.LeeS. H.