Studying the properties of baryon-rich quark-gluon plasma (QGP) is the main focus of the beam energy scan (BES) experiments [1-3] at the Relativistic Heavy Ion Collider (RHIC) as well as at the future Facility for Antiproton and Ion Research (FAIR) and the Nuclotron-based Ion Collider Facility (NICA). There have been many theoretical studies on the QCD phase diagram and its structure, based on, for example, the lattice QCD [4], holographic models [5], and effective field theories [6]. Also, the QCD phase diagram has been explored in the hydrodynamic framework by the BEST Collaboration [7] and used in machine learning methods to study its effects in heavy ion collisions [8]. In particular, using the hydrodynamic [9-12] or transport [13] model with an equation of state that has a first-order partonic to hadronic matter phase transition, it has been shown that the produced matter in heavy ion collisions could be mechanically unstable due to the spinodal instability associated with the first-order phase transition. This would then lead to an amplification of the density inhomogeniety, and the resulting density fluctuations could enhance the production of composite particles, such as hadrons and nuclei, which can then be used as a signal of the first-order phase transition. For example, using a partonic transport model based on a three-flavor Nambu-Jona-Lasinio (NJL) model [14, 15] for Au+Au collisions with the initial quark temperature and density in the spinodal region of the corresponding strong-interaction QCD phase diagram, it was shown in Ref. [13] that the quark density fluctuation after the evolution would be larger for the case with a first-order chiral phase transition than for the case without a first-order chiral phase transition in the equation of state. These results are illustrated in panels (b) and (a) of Fig. 1, respectively.

Fig. 1Full-screen View

(Color online) Density distribution of strongly interacting matter in a heavy ion collision after its expansion for the cases of crossover transition (panel (a)) and first-order chiral phase transition (panel (b)). Also shown for illustration of the latter case are deuterons and tritons produced from the density fluctuating hadronic matter and their yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}=N_{\text{t}}{N}_{\text{p}}/N_{\text{d}}^2$, which depends on the magnitude of neutron density distribution as discussed in the text.

Also shown in Fig. 1 are the light nuclei produced from the the hadronic matter resulting from the clumping quark matter after hadronization and further evolution as well as their yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}=N_{\text{t}}{N}_{\text{p}}/N_{\text{d}}^2$, where Np, N_d, and N_t denote, respectively, the proton, deuteron, and triton numbers. According to Ref. [16], where the coalescence model [17-20] is used for light nuclei production, this yield ratio is related to the neutron relative density fluctuation $\text{Δ}n=\langle {(\delta n)}^2\rangle /{\langle n\rangle }^2$, with $\langle n\rangle$ and $\langle {(\delta n)}^2\rangle$ being the average neutron density and its variance, respectively, by ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}\approx \frac{1}{2\sqrt{3}}(1+\text{Δ}n)$. One thus expects an enhanced ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ if the neutron density fluctuation remains present at the kinetic freeze-out of produced hadronic matter. A similar density fluctuation effect on the yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ also holds if deuteron and triton are produced thermally or statistically at the kinetic freeze out of the hadronic matter [21, 22]. It was further argued in Ref. [23] that the critical fluctuations in the vicinity of the critical end point (CEP) of a first-order phase transition could also lead to an enhancement of the yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ as shown in Ref. [24]. An enhanced ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ has also been interpreted as a signal for the CEP according to Refs. [25, 26], which showed that light nuclei production would be enhanced as a result of their increased binding energies when the attractive part of nucleon-nucleon potential due to the σ meson exchange becomes deeper because of its decreasing mass as the system is close to the CEP. Reviews on this very interesting topic of the relation of light nuclei and the related topic of net baryon number fluctuations to the QCD phase transition can be found, respectively, in Ref. [27] and Refs. [28-30]. For a recent review on experimental study of the QCD phase diagram in relativistic heavy-ion collisions, it can be found in Ref. [31].

Recently, the STAR Collaboration has published in Physical Review Letters [34] the yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ from the data measured during the first phase of BES experiments. This result is based on the analysis led by Daniel Cebra, Matthew Harasty, Hui Liu, Xiaofeng Luo, Nu Xu, Ning Yu, and Dingwei Zhang. As shown in Fig. 2 for Au+Au collisions in the center-of-mass energy range of $\sqrt{s_{\text{NN}}}=7.7-200$ GeV, this yield ratio is enhanced at 19.6 and 27 GeV in the most central collisions, although it shows a monotonic collision energy dependence. This result has been compared with theoretical predictions of the coalescence model for deuteron and triton production based on kinetically freeze-out protons and neutrons from microscopic models for relativistic heavy ion collisions. The latter include the hybrid model based on the (3+1)D viscous hydrodynamic model MUSIC for the quark-gluon plasma and the UrQMD transport model for the hadronic matter [32, 35] as well as a multiphase transport (AMPT) model that includes both the partonic and hadronic phases [33]. With a smooth crossover transition between the quark-gluon plasma and hadronic matter, an essentially collision energy independent ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ is predicted by these models for both peripheral and central Au+Au collisions. There have been attempts to extend the study of Ref. [13] based on the NJL model for the quark matter to include its hadronization and the evolution of the resulting hadronic matter via the AMPT model [21, 36]. These studies have shown that the quark density fluctuation can largely survive hadronization but the hadron density fluctuation is somewhat washed out by their scatterings. The resulting yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ from the nucleon coalescence model based on kinetically freeze-out nucleons can, however, still be enhanced if the system goes through the spinodal region of the QCD phase diagram from the NJL model. Because of the low critical temperature and high baryon chemical potential in the NJL model, the enhanced ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ only happens in central Au+Au collisions at center-of-mass energies lower than where the peak is seen in the STAR data. More realistic models for the QCD equation of state are thus needed in this transport model study. These equations of state can also be used in the MUSIC+UrQMD hybrid model to study the collision energy dependence of ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$ and to extract the critical temperature and baryon chemical potential from the STAR data.

Fig. 2Full-screen View

(Color online) Collision energy, p_T dependence of the yield ratio $N_{\text{t}}{N}_{\text{p}}/N_{\text{d}}^2$ in Au+Au collisions at RHIC for 0%-10% central (left panel) and 40%-80% peripheral (right panel) collisions. Dashed lines are the coalescence baselines obtained from the coalescence-inspired fit. Shaded areas denote the calculations from the MUSIC+UrQMD hybrid [32] and the AMPT [33] model. Taken from Ref. [34]

Besides the yield ratio ${\mathcal{O}}_{\text{p}-\text{d}-\text{t}}$, there are other yield ratios of light nuclei that are also sensitive to nucleon density fluctuations and can be used to probe the QCD phase diagram. For example, the yield ratio $\frac{N_{\alpha }{N}_{\text{p}}}{N_{{}^3\text{He}}{N}_{\text{d}}}\approx \frac{2\sqrt{2}}{9\sqrt{3}}(1+\text{Δ}p)$ is sensitive to the proton density fluctuation Δ p [24] and also to the closeness to the critical point [37]. To avoid the smearing effect of hadronic scatterings on the nucleon density fluctuations and to directly probe the large density fluctuation during the quark to hadronic matter first-order phase transition, one can consider the yield ratio of hadrons, such as $\frac{N_p{N}_{{\overline{K}}^0}}{N_{{\pi }^{+}}{N}_{\Lambda }}$ and $\frac{N_{K^{+}}{N}_{{\Xi }^{-}}}{N_{\varphi }{N}_{\Lambda }}$, which can be shown to be sensitive to the up quark and strange quark density fluctuations, respectively, if they are produced through quark coalescence [38, 39] and their numbers do not change during the hadronic evolution as in the statistical hadronization model [40, 41]. It is worthwhile to mention that in intermediate-energy nuclear collisions, density fluctuations of produced warm nucleonic matter can be used to study the nuclear liquid-gas phase transition [42] and light nuclei produced in these collisions can be employed to probe the density dependence of nuclear symmetry energy [43, 44], which has importance implications in both nuclear structure and nuclear astrophysics.