We briefly overview motivations, some recent results and challenges in studying light nuclei production in relativistic heavy ion collisions.

We present a sampling method for the transition from relativistic hydrodynamics to particle transport, commonly referred to as particlization, which preserves the local conservation of energy, momentum, baryon number, strangeness, and electric charge microcanonically, i.e., in every sample. The proposed method is essential for studying fluctuations and correlations by means of stochastic hydrodynamics. It is also useful for studying small systems. The method is based on Metropolis sampling applied to particles within distinct patches of the switching space-time surface, where hydrodynamic and kinetic evolutions are matched.

We present a sampling method for the transition from relativistic hydrodynamics to particle transport, commonly referred to as particlization, which preserves the local event-by-event conservation of energy, momentum, baryon number, strangeness, and electric charge. The proposed method is essential for studying fluctuations and correlations by means of stochastic hydrodynamics. It is also useful for studying small systems. The method is based on Metropolis sampling applied to particles within distinct patches of the hypersurface, where hydrodynamics and kinetic evolutions are matched.

Many simulations of relativistic heavy-ion collisions involve the switching from relativistic hydrodynamics to kinetic particle transport. This switching entails the sampling of particles from the distribution of energy, momentum, and conserved currents provided by hydrodynamics. Usually, this sampling ensures the conservation of these quantities only on the average, i.e., the conserved quantities may actually fluctuate among the sampled particle configurations and only their averages over many such configurations agree with their values from hydrodynamics. Here we apply a recently invented method [D. Oliinychenko and V. Koch, Phys. Rev. Lett. 123, 182302 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.182302] to ensure conservation laws for each sampled configuration in spatially compact regions (patches) and study their effects: From the well-known (micro-)canonical suppression of means and variances to little studied (micro-)canonical correlations and higher-order fluctuations. Most of these effects are sensitive to the patch size. Many of them do not disappear even in the thermodynamic limit, when the patch size goes to infinity. The developed method is essential for particlization of stochastic hydrodynamics. It is useful for studying the chiral magnetic effect, small systems, and in general for fluctuation and correlation observables.

We demonstrate that proton and pion flow measurements in heavy-ion collisions at incident energies ranging from 1 to 20 GeV per nucleon in the fixed target frame can be used for an accurate determination of the symmetric nuclear matter equation of state at baryon densities equal 2-4 times nuclear saturation density n0. We simulate Au+Au collisions at these energies using a hadronic transport model with an adjustable vector mean-field potential dependent on baryon density nB. We show that the mean field can be parametrized to reproduce a given density dependence of the speed of sound at zero temperature cs2(nB,T=0), which we vary independently in multiple density intervals to probe the differential sensitivity of heavy-ion observables to the equation of state at these specific densities. Recent flow data from the STAR experiment at the center-of-mass energies sNN={3.0,4.5}GeV can be described by our model, and a Bayesian analysis of these data indicates a hard equation of state at nB∈(2,3)n0 and a possible phase transition at nB∈(3,4)n0. More data at sNN=2-5GeV, as well as a more thorough analysis of the model systematic uncertainties will be necessary for a more precise conclusion.

We present a method that may allow an estimate of the value of the speed of sound as well as its logarithmic derivative with respect to the baryon number density in matter created in heavy-ion collisions. To this end, we utilize well-known observables: cumulants of the baryon number distribution. In analyses aimed at uncovering the phase diagram of strongly interacting matter, cumulants gather considerable attention as their qualitative behaviour along the explored range of collision energies is expected to aid in detecting the QCD critical point. We show that the cumulants may also reveal the behavior of the speed of sound in the temperature and baryon chemical potential plane. We demonstrate the applicability of such estimates within two models of nuclear matter, and explore what might be understood from known experimental data.

Cumulants of baryon number are given considerable attention in analyses of heavy-ion collision experiments as possible signatures of the QCD critical point. In this work, we show that the values of the lowest three cumulants can also be utilized to recover information about the isothermal speed of sound and its logarithmic derivative with respect to the baryon number density. This result provides a new method for obtaining information about fundamental properties of nuclear matter studied in heavy-ion collisions, with consequences for both the search for the QCD critical point and neutron star studies. While the approximations and the model comparison we considered apply to experiments at low energies, the approach itself can be used at any collision energy provided that measurements of cumulants of baryon number distribution as well as their temperature dependence are available.

We present a method that may allow an estimate of the value of the speed of sound as well as its logarithmic derivative with respect to the baryon number density in matter created in heavy-ion collisions. To this end, we use well-known observables: cumulants of the baryon number distribution. In analyses aimed at uncovering the phase diagram of strongly interacting matter, cumulants gather considerable attention as their qualitative behavior along the explored range of collision energies is expected to aid in detecting the QCD critical point. We show that the cumulants may also reveal the behavior of the speed of sound in the temperature and baryon chemical potential plane. We demonstrate the applicability of such estimates within two models of nuclear matter and explore what might be understood from known experimental data.

The deuteron yield in Pb+Pb collisions at sNN=2.76TeV is consistent with thermal production at a freeze out temperature of T=155MeV. The existence of deuterons with binding energy of 2.2 MeV at this temperature was described as "snowballs in hell" [P. Braun-Münzinger, B. Dönigus, and N. Löher, CERN Courier, August 2015]. We provide a microscopic explanation of this phenomenon, utilizing relativistic hydrodynamics and switching to a hadronic afterburner at the above-mentioned temperature of T=155MeV. The measured deuteron pT spectra and coalescence parameter B2(pT) are reproduced without free parameters, only by implementing experimentally known cross sections of deuteron reactions with hadrons, most importantly πd↔πnp.