Open Access

The statistical model with exact conservation of baryon number, electric charge, and strangeness – the Canonical Statistical Model (CSM) – is used to analyze the dependence of yields of light nuclei at midrapidity on charged pion multiplicity at the LHC. The CSM calculations are performed assuming baryon-symmetric matter, using the recently developed Thermal-FIST package. The light nuclei-to-proton yield ratios show a monotonic increase with charged pion multiplicity, with a saturation at the corresponding grand-canonical values in the high-multiplicity limit, in good qualitative agreement with the experimental data measured by the ALICE collaboration in pp and Pb–Pb collisions at different centralities and energies. Comparison with experimental data at low multiplicities shows that exact conservation of charges across more than one unit of rapidity and/or a chemical freeze-out temperature which decreases with the charged pion multiplicity improves agreement with the data.

The thermal fit to preliminary HADES data of Au+Au collisions at sNN=2.4 GeV shows two degenerate solutions at T≈50 MeV and T≈70 MeV. The analysis of the same particle yields in a transport simulation of the UrQMD model yields the same features, i.e. two distinct temperatures for the chemical freeze-out. While both solutions yield the same number of hadrons after resonance decays, the feeddown contribution is very different for both cases. This highlights that two systems with different chemical composition can yield the same multiplicities after resonance decays. The nature of these two minima is further investigated by studying the time-dependent particle yields and extracted thermodynamic properties of the UrQMD model. It is confirmed, that the evolution of the high temperature solution resembles cooling and expansion of a hot and dense fireball. The low temperature solution displays an unphysical evolution: heating and compression of matter with a decrease of entropy. These results imply that the thermal model analysis of systems produced in low energy nuclear collisions is ambiguous but can be interpreted by taking also the time evolution and resonance contributions into account.

The theory of strong interactions — Quantum Chromodynamics (QCD) — is well-defined mathematically. However, direct applications of this theory to experiment are rather limited due to significant technical obstacles. Even some general features of QCD remain unclear to date. Hence, phenomenological input is important and needed for practical applications, e.g. for theoretical analysis of the heavy-ion collision experiments. In this thesis the role of hadronic interactions is studied in the hadron resonance gas (HRG) model — a popular model for the confined phase of QCD. The description of hadronic interactions is based on the famous van der Waals (VDW) equation and its quantum statistical generalization. While this is not the conventional choice for nuclear/hadronic physicspplications, the simplicity of the VDW approach makes it extremely useful. In particular, this framework allows to include the two most basic ingredients of hadron-hadron interaction: the short-range repulsion, modeled by excluded-volume (EV) corrections, and the intermediate range attraction. The first part of the thesis considers just the repulsive EV interactions between hadrons. A hitherto unknown, but surprisingly strong sensitivity of the long known thermal fits to heavy-ion hadron yield data to the choice of hadron eigenvolumes is uncovered. It challenges the robustness of the chemical freeze-out temperature and baryochemical potential determination from the thermal fits. However, at the same time, the extracted value of the entropy per baryon is found to be a robust observable which depends weakly on this systematic uncertainty of the HRG model. A Monte Carlo procedure to treat EV interactions in HRG is also introduced in this thesis. It allows to study simultaneous effects of EV and of exact charge conservation in HRG for the first time. Generalizations of the classical VDW equation are required for its applications in hadronic physics. he grand canonical ensemble (GCE) formulation of the classical VDW equation is presented. Remarkably, this important aspect of the VDW equation was not discovered before. The GCE formulation yields the analytic structure of the critical fluctuations, both in the vicinity of and far off the critical point. These critical fluctuations are presently actively being used as probes for the QCD critical point. Another extension is the hitherto undiscovered generalization of the VDW equation to include quantum Bose-Einstein and Fermi-Dirac statistics. It is performed for both single-component and multi-component fluids. The Fermi-Dirac VDW equation is applied for the first time. It is used to describe nucleons and basic properties of nuclear matter. The quantum statistical generalization of the VDW equation developed in this work is quite general, and can be applied for any fluid. Thus, its applications are not restricted to QCD physics, but may also find themselves in chemistry and/or industry. The quantum statistical VDW equation is used to describe baryonic interactions in full HRG. The VDW parameters $a$ and $b$ are fixed to the nuclear ground state and the predictions of the model are confronted with lattice QCD calculations. The inclusion of baryonic interactions leads to a qualitatively different behavior of the fluctuations of conserved charges in the crossover region. In many cases it resembles the lattice data. These results suggest that hadrons do not melt quickly with increasing temperature, as one could conclude on the basis of the common simple ideal HRG model. Calculations at finite chemical potentials show that the nuclear liquid-gas transition manifests itself by non-trivial fluctuations of the net baryon number in heavy ion collisions. In the final part of the thesis the pure glue initial scenario for high-energy hadron and heavy-ion collisions is explored. This scenario is shown not to spoil the existing agreement of the hadronic and electromagnetic observables description in Pb+Pb collisions at energies available at the CERN Large Hadron Collider. Hydrodynamic calculations suggest that collisions of small-sized nuclei at lower collision energies available at the BNL Relativistic Heavy Ion Collider are promising in the search for the traces of the chemically non-equilibrium gluon-dominated phase transition.

We estimate the feeddown contributions from decays of unstable A=4 and A=5 nuclei to the final yields of protons, deuterons, tritons, 3He, and 4He produced in relativistic heavy-ion collisions at sNN>2.4 GeV, using the statistical model. The feeddown contribution effects do not exceed 5% at LHC and top RHIC energies due to the large penalty factors involved, but are substantial at intermediate collision energies. We observe large feeddown contributions for tritons, 3He, and 4He at sNN≲10 GeV, where they may account for as much as 70% of the final yield at the lower end of the collision energies considered. Sizable (>10%) effects for deuteron yields are observed at sNN≲4 GeV. The results suggest that the excited nuclei feeddown cannot be neglected in the ongoing and future analysis of light nuclei production at intermediate collision energies, including HADES and CBM experiments at FAIR, NICA at JINR, RHIC beam energy scan and fixed-target programmes, and NA61/SHINE at CERN. We further show that the freeze-out curve in the T-μB plane itself is affected significantly by the light nuclei at high baryochemical potential.

The production of light (anti-)(hyper-)nuclei in heavy-ion collisions at the LHC is considered in the framework of the Saha equation, making use of the analogy between the evolution of the early universe after the Big Bang and that of “Little Bangs” created in the lab. Assuming that disintegration and regeneration reactions involving light nuclei proceed in relative chemical equilibrium after the chemical freeze-out of hadrons, their abundances are determined through the famous cosmological Saha equation of primordial nucleosynthesis and show no exponential dependence on the temperature typical for the thermal model. A quantitative analysis, performed using the hadron resonance gas model in partial chemical equilibrium, shows agreement with experimental data of the ALICE collaboration on d, 3He, HΛ3, and 4He yields for a very broad range of temperatures at T≲155 MeV. The presented picture is supported by the observed suppression of resonance yields in central Pb–Pb collisions at the LHC. Keywords: Light (anti-)(hyper-)nuclei production, Saha equation, Partial chemical equilibrium.