Collective Excitation in High-Energy Nuclear Collisions -- In Memory of Professor Lianshou Liu

We celebrate the legacies of our friend and mentor Professor Lianshou Liu who was one of the pioneers for the phenomenology of multi-particle interactions and initiated the physics of relativistic heavy-ion collisions in China. In this article, we discuss some of the recent exciting experimental observations on the collective phenomena including collectivity, chirality, criticality, strangeness production, and thermal equilibrium in high-energy nuclear collisions. Future directions, especially the physics at high baryon density, will be discussed with a focus on the first-order phase boundary and hyperon-nucleon interactions.


Collectivity: Azimuthal Angular Anisotropy in High-Energy Nuclear Collisions
Collective flow defined by the coefficiencies of the Fourier expansion of final particle distribution in momentum space is sensitive to the early stage of nuclear collisions. Specifically, the first three coefficiencies are called directed flow (v 1 ), elliptic flow (v 2 ), and triangular flow (v 3 ), respectively. Directed flow is sensitive to the Equation of State (EoS) of the medium; elliptic flow is sensitive to the degree of freedom, partonic or hadronic level, and degree of equilibrium of the medium; triangular flow is sensitive to the initial geometry fluctuations. A comprehensive set of measurements have been achieved in RHIC-STAR experiment of nuclear collisions [1][2][3][4][5][6][7][8][9]. The Number of Constituent Quark (NCQ) scaling of v n observed at high energy collisions (>20 GeV) indicates that the partonic collectivity has been built-up [1][2][3]8,10]. In particular, the D meson also follows the NCQ scaling [2,10,11], suggesting that the charm quark collectivity is at the same level as that of u, d, and s quarks; therefore, the created medium reaches (nearly) equilibrium.
The main motivation of Beam Energy Scan (BES) program is to explore the QCD phase diagram and search for the possible phase boundary and critical point. The first phase of BES program (BES-I) at STAR experiment covers collision energy √ s NN = 7.7-62.4 GeV.
Lots of interesting phenomena have been observed; here, we focus on the collective flow v n measurements. Figure 1 summarizes the directed, elliptic, and triangular flow relevant observations from STAR BES-I. The v 1 slope of net-baryon near mid-rapidity as a function of collision energy is regarded as a possible signal of first-order phase transition. The nonmonotonic energy dependence of v 1 slope is associated with the phase transition and the minimum of v 1 slope is called the "softest point collapse" [12]. In the experiment, as neutrons are hard to be measured, we use net-proton as a proxy of net-baryon. The left panel of Figure 1 shows the v 1 slope of net-proton, net-Kaon, and net-Λ as a function of collision energy. The non-monotonic behavior is observed for net-proton and net-Λ, and the minimum occurs near √ s NN of 10-20 GeV. On the other hand, the net-Kaon slope shows a monotonic increase from low to high energy. Above √ s NN = 20 GeV, the net-Kaon slope is overlapping with net-proton and net-Λ, while the divergence happens below 20 GeV. Further investigation is needed to understand the physical mechanism of this divergence at the low energy region. In addition, the centrality dependence measurements of net-particle v 1 slope are crucial to understand the effects from mechanism which is not related to phase transition [13][14][15]. The v 2 of particles and corresponding anti-particles is also found to be significantly different below 20 GeV as shown in the middle panel of Figure 1. Even though the observed NCQ scaling is not valid for particles and anti-particles, roughly scaling still works for particle group and anti-particle group separately. One explanation of the v 2 difference is the transported effect. The final particles are a mixture of transported and produced particles, while the dominant part produces particles at high energy collisions, such as 200 GeV at RHIC and a few TeV at LHC. With the decrease of collision energy, there are more final particles from transportation. Since the transported particles undergo the whole collision evolution, the collectivity should be different with the produced particles. Several model studies try to explain the v 2 difference of particles and antiparticles and succeed partly [16][17][18][19][20][21][22]. More detailed measurements as a function of centrality and for (multi-)strange hadrons help us pick out the right mechanism. A similar difference is also observed for v 3 , as shown in the right panel of Figure 1. The collision energy and particle type dependence are the same as v 2 . Recently, the v 1 and v 2 results of identified particles from √ s NN = 3 GeV suggest that the hadronic interactions are dominant at such energy [23][24][25].
Comprehensive measurements of collective flow from collision energy <20 GeV, such as STAR experiment BES-II (3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20) and experiments at SPS and SIS18 [26][27][28], help us further constrain the phase transition boundary.  The energy dependent experimental data on the collective flow can be used to extract important properties of the medium. The temperature dependence of the event-averaged shear viscosity-to-entropy ratio 4πη/s [30] is shown in Figure 2. In the left panel, chemical freeze-out temperature from each energy [31] is used and normalized to that from 200 GeV Au+Au collisions. As one can see, in the high energy limit, √ s NN = 39-200 GeV, the ratio reaches unity, the quantum limit, implying that the medium created in such collisions is dominated by partonic interactions with a minimum value of 4πη/s. At lower collision energies, on the other hand, hadronic interactions are dominant, and the medium shows a rapid increase of the viscosity-to-entropy ratio. The right panel is taken from Ref. [30], where the temperature evolution of the shear ratio is shown as a function of the scaled temperature T/T c . Here, T c represents the critical temperature in the calculation [32,33]. The observed V-shaped feature is quite similar to what is described in Ref. [34] and can be taken as experimental evidence of the expected crossover transition in QCD. The effective values of shear viscosity-to-entropy ratio, 4πη/s, shown as a function of the scaled temperature. The horizontal dashed line indicates the quantum lower limit. Left panel: the extracted 4πη/s from the energy dependence of the measured v 2 [6] and v 3 [35], shown as the scaled chemical freeze-out temperature T ch /T ch (200 GeV). Right panel: temperature evolution of 4πη/s, extracted from Bayesian analyses [32,33].

Chirality: Chiral Magnetic Effect in High-Energy Nuclear Collisions
The QCD vacuum can exist in states of varying topological charge or Chern-Simons winding number [36][37][38]. Sphaleron transitions between those states, strongly suppressed under normal conditions, may have appreciable probability at high temperatures, such as those achieved in high-energy nuclear collisions. Local metastable domains of gluon fields with non-zero topological charges may thus form in those collisions [39]. Interactions of quarks with those gluon fields can cause an imbalance in chirality, which may manifest as an electric current along a strong magnetic field, a phenomenon called the Chiral Magnetic Effect (CME) [40]. Because those local domains break the parity (P) and charge-parity (CP) symmetries and because those domains arise ultimately from vacuum fluctuations, out of which our universe is believed to come into being, an observation of the CME in high-energy nuclear collisions may unravel the mystery of the matter-antimatter asymmetry of our universe [41].
The CME has been extensively sought at the RHIC at BNL and the LHC at CERN [42,43]. A signature of the CME is back-to-back charge separation along the magnetic field, produced mainly by spectator protons in relativistic heavy-ion collisions [44]. Because the magnetic field is on average perpendicular to the reaction plane (RP), a commonly used observable is the three-point azimuthal correlator [45], where φ α and φ β are the azimuthal angles of particles α and β of either the opposite sin (OS) or the same sign (SS). The ψ RP is the azimuthal angle of the RP and is usually reconstructed from final-state particles, whose inaccuracy is corrected by a resolution factor [46]. Equivalently, γ can also be calculated by the three-particle correlator where v 2,c is the elliptic flow parameter of the third particle c. Several other observables have been proposed and were found to be similar to ∆γ [47].
Charge separation measurements revealed strong ∆γ signals [48][49][50][51], approximately independent of collision energy except at lower RHIC energies where the signal dies off [52]. However, major backgrounds exist that arise from two-particle correlations coupled with the finite elliptic flow of those background correlation sources [53][54][55], Here, N 2p and v 2,2p refer, respectively, to the number and elliptic flow of those two-particle correlation sources, such as resonances and jets, and φ 2p refers to their azimuths; N α,β are the α, β particle multiplicities. Experimental indication of major backgrounds comes from small-system measurements [56,57]. Clear experimental evidence of those backgrounds is observed by the STAR measurement of ∆γ as a function of pair invariant mass [58,59]. It became clear because backgrounds are dominant that data-driven methods must be invoked to reliably extract the possible CME signal. An initial attempt was completed by STAR [60] by analyzing data as a function of the event-by-event ellipticity of the particles of interest (i.e., α and β). A similar event-shape engineering (ESE) analysis [61] was carried out by ALICE [62] and CMS [63] by analyzing data as a function of v 2 of the particles of interest at midrapidity in events selected according to the ellipticity in the forward-rapidity region. Data are consistent with vanishing CME signals and upper limits, on the order of 20% of the inclusive ∆γ measurements, have been extracted.
To eliminate large backgrounds, experiments often rely on comparative or relative measurements. One such relative measurement is the isobar 44 96 Ru+ 44 96 Ru and 40 96 Zr+ 40 96 Zr collisions conducted in 2018 at RHIC [64]. Because of the equal mass numbers of the isobar nuclei, the physics backgrounds are expected to be the same; the 10% difference in the atomic numbers is expected to result in a difference in the magnetic fields and thus in the CME signal strengths between the two systems [65]. A blind analysis was conducted [66], and the results [67] are shown in Figure 3 with an unprecedented precision of 0.4%. The Ru+Ru/Zr+Zr ratios of ∆γ/v 2 (motivated by Equation (2)) are all below unity, which has been understood to be due mainly to the multiplicity difference between the two isobar systems. This difference was predicted by energy density functional theory calculations [68,69] to arise from the smaller size of the 44 96 Ru nucleus compared to that of the 40 96 Zr nucleus [70]. If the ∆γ is inversely proportional to multiplicity, then the baseline for the double ratio would be the bottom dashed line in Figure 3. The measured double ratios are all above this line, seemingly suggesting finite CME signals, as also pointed out in Ref. [71]. However, the inverse multiplicity scaling is only approximate because the background in ∆γ scales with the number of correlation sources (see Equation (2)), which may not be strictly proportional to multiplicity. A more realistic baseline may be the Ru+Ru/Zr+Zr ratio of the pair excess r = (N OS − N SS )/N OS [67] indicated by the middle dashed line. If so, then there is no evidence of a CME signal in the isobar data. To complicate the matter further, background contamination beyond that of Equation (2) exists because of genuine three-particle correlations and nonflow contamination in v 2 [72]. Preliminary estimates [73] of those nonflow effects are shown by the shaded bands in Figure 3, indicating the final baselines for the double ratios measured by the full event and subevent methods. The isobar data are consistent with these baselines within approximately one standard deviation, suggesting that the CME signal in isobar collisions cannot be larger than a few percent-the 0.4% data precision translates into a f CME = 2-3% CME fraction of the inclusive ∆γ measurement. Another comparative method consists of measurements with respect to the spectator plane (SP) and the participant plane (PP) [76,77]. Because the SP is better aligned with the perpendicular orientation of the magnetic field and the PP is determined by the elliptic flow harmonic direction, ∆γ measurements w.r.t. SP and PP contain different contributions from the CME signal and flow background [76]. One can obtain the CME fraction by f CME ≡ where a is the ratio of v 2 measurement w.r.t. SP to that w.r.t. PP, which quantifies the angular spread between SP and PP, and A is the ratio of ∆γ measurement w.r.t. SP to that w.r.t. PP. Figure 4 shows the f CME and the absolute ∆γ CME signal strength in Au+Au collisions at 200 GeV by STAR [78]. While the peripheral 50-80% data are consistent with vanishing CME, the midcentral 20-50% results indicate a finite CME signal of the order of 10%, with a significance ranging from 1 to 3 standard deviations. While the flow-induced background is removed, the results are still contaminated by nonflow contributions. Model estimates [72] of those nonflow contributions indicate a consistent-with-zero contribution to the full-event result and a negative contribution to the subevent results. Estimation utilizing real data is ongoing. If nonflow contamination is small or even negative, then the positive f CME observation by STAR may indeed indicate a finite CME signal. Note this does not contradict the null signal from the isobar data because the signal to background ratio in isobar collisions is expected to be significantly smaller, perhaps by a factor of three [79], than in Au+Au collisions, due to the smaller multiplicity (thus larger background) and the smaller magnetic field strength (thus smaller CME strength) in isobar collisions [79]; the f CME ∼ 10% signal in Au+Au collisions and the less than a few percent signal in isobar collisions are, in fact, consistent with each other. In summary, the CME-sensitive ∆γ measurements in heavy-ion collisions are dominated by flow-induced backgrounds arising from two-particle correlations coupled with elliptic flow v 2 . Further background contamination arises from three-particle correlations and nonflow effects in v 2 . Methods have been developed to eliminate those backgrounds, including eventshape engineering, isobar collisions, and measurements w.r.t. spectator plane SP and participant plane PP. While the former two have yielded null CME signal with the present statistics, the comparative measurements w.r.t. SP and PP indicate a hint of f CME ∼ 10% CME signal in Au+Au collisions at 200 GeV, with a significance of 1-3 standard deviations. In the future, an order of magnitude increase in Au+Au statistics is anticipated from 2023 and 2025 by STAR. This would present a powerful data set to either identify the CME or put a stringent upper limit on it. Such an outcome would significantly advance our understanding of QCD and our universe.

High Order Moments and Search for the QCD Critical Point
Exploring the phase structure of hot and dense QCD matter is one of the main goals of heavy-ion collisions. In the large baryon density region, the QCD critical point is the endpoint of the first-order phase transition boundary and a landmark in the QCD phase diagram. Mapping out the first-order phase boundary and pinning down the location of the QCD critical point will enhance our understanding of universe evolution and the structure of visible matters. During the last two decades, many efforts have been made both experimentally [80][81][82][83][84][85][86][87][88][89][90][91] and theoretically [92][93][94][95][96][97][98][99][100][101] to determine the location of the QCD critical point. In this review, we will focus on the measurement of high-order moments with the data taken in the first phase of the beam energy scan program at RHIC. For a quantum system, the thermodynamic pressure [102] can be expressed via temperature T and chemical potentials µ i , where (i = B, Q, S) stands for conserved quantum numbers of baryon (B), electric charge (Q), and strangeness (S) in heavy-ion collisions: Here, k, g i , and ± signs stand for the particle momentum, the degeneracy, and the nature of Fermion (+) or Boson (−). The sum goes over all particles in the equilibrium with masses up to 2.5 GeV from the PDG.
For a conserved quantity, the cumulants in thermal equilibrium can be expressed through the generalized thermodynamic susceptibilities: and the n th order generalized susceptibilities are the derivatives of the pressure: Fluctuations and correlations among conserved charges, i.e., baryon number (B), electric charge (Q) and strangeness (S), are sensitive observables to probe the QCD phase structure. Experimental proxies, such as net-protons, net-kaons, are used for measurements of mean (M), variance (σ 2 ), skewness (s), and kurtosis (κ) of conversed charges [82,83,85]. As an example, the connections between moments and cumulant (C n ) ratios of conserved charges are listed as the following: where i = B, Q and S are the conserved charges, C i n is the n th order cumulant, and χ i n (T, µ i ) is the generalized susceptibility defined as the n th order derivative of pressure with respect to chemical potential µ i .
One of the most important advantages of Equations (6)-(8) is that they connect the measurements, the experimentally measured moments on the left sides of the equations, to the ratios of the thermodynamic susceptibilities (and ratios of cumulants) from theoretical calculations [98,[103][104][105][106][107]. These constructions offer not only a sensitive probe for studying the QCD phase structure in high-energy collisions but also to test the limits of thermalization in such collisions. The latter is a necessary step for us to understand the emerging macroscopic features from the violent microscopic scatterings.
We start first from experimental observations. As a summary, Figure 5 shows recent results on the fourth-order net-proton (filled red circles) and proton (open squares) high moments in central Au+Au collisions from the RHIC BES program and the HADES experiment [86,89,108] together with model comparisons. The thin solid red and dot-dashed blue lines depict a qualitative prediction for the behavior of the fourth-order net-proton cumulant ratio C 4 /C 2 (κσ 2 ) due to an evolution in the vicinity of a critical region [94]; the locations of the peak of the dot-dashed blue curve and the dip of the solid red curve are chosen only for a qualitative comparison. Intriguing non-monotonic behavior is observed in the data, and while large error bars prevent one from making a more decisive conclusion, the behavior of the fourth-order netproton cumulant ratio in the energy range √ s NN = 7.7-27 GeV seems to significantly deviate from the non-critical baseline provided by models. In particular, both the hadronic transport model UrQMD (cascade mode) [109,110] (gold band) and a thermal model in the canonical ensemble [111] (dot-dashed red line) predict, for decreasing collision energy, a monotonic suppression in the fourth-order cumulant ratio C 4 /C 2 due to baryon number conservation, in contrast to the behavior tentatively seen in data. In addition, the experimental data at √ s NN ≤ 20 GeV indicate an excess of two-proton correlations as compared to a non-critical baseline including effects due to, e.g., baryon number conservation [112]. At 3 GeV, the agreement between data and transport model calculation, see the cross in the figure, implies that hadronic interaction dominates the properties of the medium [89,113]. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV. The green band, in Figure 5, covers collision energy √ s NN = 7.7-19.6 GeV and is the estimated statistical uncertainty from the RHIC beam energy scan phase II program. As one can see, in order to complete the beam energy scan and determine if the QCD critical point exists, the CBM experiment at FAIR is necessary to fill the energy gap between √ s NN = 3-8 GeV [114]; see the hatched region in Figure 5. Finally, it is worth noting that, in low-energy collisions, net-proton cumulants have also been recently proposed as a means for extracting the speed of sound and its logarithmic derivative [115].   [111] and transport model UrQMD [109,110] are shown as the red dot dashed line and gold band, respectively. The thin solid red and dot-dashed blue lines depict a qualitative prediction for the behavior of the fourth-order net-proton cumulant due to an evolution in the vicinity of a critical region [94].

Limits of Thermalization in High-Energy Nuclear Collisions
In the previous subsection, the energy dependence of the high order moments of netprotons for the search of the QCD critical point in the beam energy scan program at RHIC is discussed. These high order moments, on the other hand, can also be used to test the limits of thermalization in such high-energy collisions [116]. In the past, yields of hadrons from high-energy nuclear collisions have been used to fit to the results of thermal model calculations in order to extract the freeze-out parameters; see, for example, discussions in Ref. [117]. The mean value of hadron yield is the first order of its multiplicity distributions. Certainly, the thermal analysis with the first moment M i is necessary, but not sufficient, for understanding the dynamics that lead to the macroscopic thermal behavior in high-energy nuclear collisions. Experimentals of high order moments have to be used in the analysis. Below, we report recent progress [118].
In Equations (6)-(8), the terms on the left and in the middle are experimentally measured quantities while the generalized thermodynamic susceptibilities, shown on the right side of these equations, can be extracted from thermodynamic calculations with a given ensemble. By comparing the associated cumulant ratios with experimental data, one then can determine freeze-out temperature and chemical potential. In addition, the test of the limit of thermalization can be performed at any order where data are available. This approach has been used in the tests of Lattice QCD calculations [103][104][105][106][107]119] and in a comparison with experimental results [118].
Here, we present a study using STAR recent published results of net-proton, net-Kaon and net-charge [82,84,85] and compared to a thermal calculation [116] with grand canonical ensemble (GCE). The resulting differences in freeze-out temperature and chemical potential are displayed in Figure 6. Left-and middle-panels show the experimental data of central (0-5%) and peripheral (70-80%) Au+Au collisions. Right-panels display the results of the thermal comparison with central collisions generated by the hadronic transport model UrQMD [109,110].   Figure 6. Comparison of the freeze-out parameters obtained by fitting to different orders of high moments of net-proton, net-Kaon, and net-charge from RHIC data [82,84,85]. The differences in freeze-out temperature (∆T) and baryon chemical potential (∆µ B ) are shown in the top and bottom row, respectively. Triangles, dots, and squares represent the difference of net-proton cumulants between 2nd and 1st, 3rd and 2nd, and 3rd and 1st. Left and middle-columns are results from experimental data. The results from hadronic transport model UrQMD [109,110] simulations are shown in the right column. Yellow bands indicate 1σ mark in differences.
The energy dependence of the net-proton fourth-order moment, κσ 2 from central Au+Au collisions (low panel) are compared with the fitting results (top panel) in Figure 7. Similarly, at high energy, both fitting results shown in the top panel and net-proton kurtosis are consistent with the Grand Canonical limits, and clear deviations appear at energy lower than 40 GeV. In the net-proton kurtosis [86], the expected criticality pattern [94], as a function of energy, seems evident; see the dot-dashed line in the lower panel of Figure 7. It is necessary to point out that criticality is a non-equilibrium phenomenon. The observed deviation from Grand Canonical limit in low energy collisions could be caused by the system passing through the critical point, but it could also be that the system is out of equilibrium. The observation of departures from thermodynamic equilibrium in the final state opens up new directions in the study of heavy-ion collisions. Further analyzing the new experimental precision data [120] with Canonical Ensemble, where the effect of baryon number conservation is included, is necessary in order to understand the underlying physics in the finite baryon density region.  Figure 7. The top panel shows the energy dependence of χ 2 /ndf from the thermal fits with either 13 data points (green triangles) or 11 data points (gray squares). The bottom panel shows the energy dependence of the net-proton fourth-order cumulant ratios C 4 /C 2 = κσ 2 [82]. The colour filling in both top and bottom panels tentatively divides the range of collision energies into regions that are clearly in agreement with the predictions of an ideal resonance gas, and therefore cannot contain the QCD critical point. Dot-dashed line is a fit result and indicates possible oscillation in the fourth-order net-proton cumulant ratios [94].

Strange Quark Probes of Parton Dynamics and QCD Interactions
Strange quark is a unique probe of the Quark-Gluon Plasma produced in high energy nucleus-nucleus collisions. Rafelski and Muller pointed out [121,122] that, in a hot QGP with a temperature above 160 MeV, the strange quark with a mass of 100 MeV/c 2 can be abundantly produced via gluon-gluon fusion to a strange quark pair. Strangeness enhancement, an increased production of strange hadrons, multi-strange hyperons in particular, has been considered a signature of equilibrated Quark-Gluon Plasma (QGP). On the other hand, if a hadronic gas system can interact and maintain thermal equilibrium hypothetically for a sufficiently long time, strangeness production can also be significantly enhanced to reach an equilibrium. Therefore, it is important to investigate the underlying dynamics for strangeness production.
Strange hadrons can be identified by their decay topologies and such particle identification capability can be extended to momentum much higher than the limit from traditional time-of-flight detectors. Measurements of elliptic flow and nuclear modification factor of identified particles up to the intermediate transverse momentum region of 5 GeV/c have been essential in the discovery of the quark number scaling [123,124]. The quark number scaling is a manifestation of the quark coalescence dynamics for particle formation at the chemical freeze-out stage. The effective degrees of freedom before the hadronization must be dominated by parton dynamics, and the partonic hydrodynamics are responsible for the development of the azimuthal angular anisotropy in the flow measurements. The change of paradigm from a QGP of free quarks and gluons in a QCD bag to a QGP of almost perfect fluid of strongly interacting quarks and gluons is one of the most important achievements of the heavy ion collision physics in the past two decades. The coalescence dynamics showed that the partonic degrees of freedom dominate the evolution dynamics of the QCD matter created in high energy nucleus-nucleus collisions.
The φ meson and Ω hyperon have played a unique role in probing the partonic dynamics of the QGP. These strange particles do not have significant hadronic interaction cross sections with other hadrons during the hadronic evolution after the chemical freeze-out. Therefore, they can carry imprints from dynamics of the partonic phase of the QGP evolution. Using coalescence dynamics, these particles have been used to probe strange quark properties in the QGP at the chemical freeze-out [125,126]. Such measurements have also been conducted with the STAR BES data to probe variations of the strange quark properties as the collision energy decreases [127]. We note that the ratio of Ω toΩ increases as the colliding energy decreases from the STAR BES measurement, and the ratio is significantly above one at low energy. There is a net baryon number in Ω hyperons although strangeness conservation dictates that strange and anti-strange quarks must be produced in pairs. The dynamics of baryon number transport to Ω hyperons are a subject of interest. Measurements of correlations between Ω and other particles such as Kaon and anti-hyperons could shed light on the roles of strangeness and baryon number conservation in baryon transport dynamics [128].
The large number of strange hyperons produced in one single heavy ion collision also provided a unique opportunity to investigate hyperon-nucleon, hyperon-hyperon interactions through correlation measurements. The STAR experiment carried out the first meaningful measurement of a Λ-Λ correlation function and extracted the interaction parameters between Λs [129]. Motivated by calculations of the HAL QCD Collaboration [130], STAR measured the p-Ω correlation function from Au+Au collisions [131]. The STAR result, based on the ratio of correlations from two centrality bins as suggested by the theoretical calculation, slightly favors the existence of a bound p-Ω system. ALICE also measured correlations of p-Ξ and p-Ω from p + p collisions at 13 TeV energy. The ALICE data indicated a strong p-Ω attractive interaction though no signal was observed for the existence of a bound state [132]. Existing measurements of hyperon correlations are not precise enough to provide definitive conclusions on the interaction strength between hyperons and the binding energy of possible composite hyperon states. The order of magnitude increase in the collision data samples in the coming years will enable us to utilize the heavy ion collisions as hyperon factories to experimentally address an important topic of hyperon interactions currently incomplete in QCD descriptions of strong interactions.
Because we can identify hyperons at high transverse momentum (p T ) region using their decay topologies, the particle-type (baryon versus meson) dependence of the nuclear modification factors can yield important insight on the parton energy loss in the QGP as well. Figure 8 shows the nuclear modification factors, R CP ratio of normalized yields of central to peripheral collisions, as a function of p T for charged hadrons, Kaons and Λ hyperons [123]. The R CP values for mesons and hyperons seem to approach each other around p T of 6 GeV/c. The disappearance of the particle-type dependence could be an important feature of the R CP : it would indicate a minimum p T above which the jet quenching dynamics may dominate, a landmark p T cut separating jet probe from hadrodynamical and coalescence regions. Improved measurements of the nuclear modification factors above this important p T cut of 6 GeV/c for identified particles in the next phase of RHIC scientific program from 2023-2025 could bring new insight on parton energy loss dynamics at RHIC.  Au+Au collisions. Solid and dot-dashed lines stand for the expected scaling of binary collisions and number of participants, respectively. There is a clear particle-type (baryon versus meson) dependence below intermediate p T presumably due to hydrodynamic flow and coalescence dynamics. The particletype dependence seems to disappear above p T of 6 GeV/c where partonic jet quenching dynamics may start to dominate.
Nucleus-nucleus collisions also opened an important venue to study hyper-nuclei dynamics, especially from the STAR BES data to a lower energy regime at future facilities. Driven by baryon density and the strangeness production threshold, both the thermal model [133,134] and transport model [135] have predicted the production of light nuclei and hyper-nuclei to peak around 3 ≤ √ s NN ≤ 10 GeV in high-energy nuclear collisions. To remove trivial factors including chemical potential and canonical effects, double ratios are often used; for instance, the ratio involving Hypertriton, Helium-3, Λ, and proton yields: The experimental results of S 3 are shown in Figure 9 along with those from model calculations [135][136][137][138]. New precision data from STAR experiment show a gradual increase as a function of the collision energy, and the value seems to approach the equilibrium limit in collisions at the LHC [139,140]. It is also interesting to note that the limiting value is about 2/3, which is commonly used in calculations of Λ-N interactions [141]. In the high baryon density region where √ s NN ≤ 10 GeV, the double ratio is further away from the thermal limit, suggesting a clear density effect. More precision measurements are needed in order to understand the Y-N and Y-Y interactions as well as their implication to the inner structure of compact stars. For recent discussions on Y-N and Y-N-N interactions, see Ref. [142] and references therein. The next generation experiment CBM at FAIR with an unprecedentedly interaction rate capability will be very important for such measurements after RHIC.

Outlook: Physics at High Baryon Density
Since the discovery of the new form of matter, the strongly coupled Quark-Gluon Plasma (QGP) [143], created in high-energy nuclear collisions in the early 2000s, scientists have been asking: "What is the structure of the QCD phase diagram in the high baryon density region?" and "Is there a QCD critical point?". Model studies have shown that a first-order phase boundary is expected at the finite baryon chemical potential µ B , while at vanishing µ B , there is a smooth crossover between QGP and hadronic matter. Thermodynamically, a critical point ought to be there at the end of the first-order phase boundary; see Figure 10. More discussions on experimental results and Lattice calculations can be found in Refs. [88,99,101,144].
The results, including the observations on collectivity, chirality, critical point, and strangeness production, from the first campaigning of the beam energy scan at RHIC show that partonic activities in central Au+Au collisions persist from √ s NN = 200 GeV to 39 GeV (corresponding to µ B /T ≤ 3) while in collisions at 3 GeV (corresponding to µ B /T ∼ 7) hadronic interactions dominate [23,86,87,89,124,125]. In addition, charge separation, a measure of the strength of the effect of CME, is vanishing in low energy collisions [52]. The QCD critical point, if existing, should be accessible in collisions at energy between 3-39 GeV. RHIC's second phase energy scan has completed data taking, and more than 10 fold statistics of Au + Au collisions have been collected for energies between 7.7-19.6 GeV. In order to complete the physics program of beam energy scan and study the phase structure of nuclear matter in high baryon density regions, future experiments such as CBM at FAIR are necessary. It will not only aid the search for the QCD critical point, but also extend the research to hyper-nuclear production which is important for studying the basic Y − N interactions. The hyperon-nuclear interaction is one of the vital connections between the physics of high-energy nuclear collisions and inner structure of compact stars.