Search

The chemical freeze-out curve in heavy-ion collisions is investigated in the context of QCD critical point (CP) search at finite baryon densities. Taking the hadron resonance gas picture at face value, chemical freeze-out points at a given baryochemical potential provide a lower bound on the possible temperature of the QCD CP. We first verify that the freeze-out data in heavy-ion collisions are well described by a constant energy per particle curve, $E/N = \rm const$, under strangeness neutrality conditions ($\mu_S \neq 0$, $\mu_Q \neq 0$). We then evaluate the hypothetical freeze-out curve based on this criterion in the absence of strangeness neutrality ($\mu_S = 0$, $\mu_Q = 0$) and confront it with recent predictions on the CP location. We find that recent estimates based on Yang-Lee edge singularities from lattice QCD data on coarse lattices place the CP significantly below the freeze-out curve and are thus disfavored by the heavy-ion data. On the other hand, predictions based on functional methods and holography place the CP slightly above the freeze-out curve, indicating that the QCD CP may be located very close to the chemical freeze-out in A+A collisions at $\sqrt{s_{NN}} = 3.5 \div 5$ GeV., Comment: 7 pages, 3 figures

We analyze particle number fluctuations in the crossover region near the critical endpoint of a first-order phase transition by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our previous study [V.A. Kuznietsov et al., Phys. Rev. C 105, 044903 (2022)] by incorporating longitudinal collective flow. The scaled variance of particle number distribution inside different coordinate and momentum space acceptances is computed through ensemble averaging and found to agree with earlier results obtained using time averaging, validating the ergodic hypothesis for fluctuation observables. Presence of a sizable collective flow is found to be essential for observing large fluctuations from the critical point in momentum space acceptances. We discuss our findings in the context of heavy-ion collisions., Comment: 13 pages, 10 figures

Subensemble Acceptance Method (SAM) [1,2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{s_{NN}}=2.4$ GeV which are found to be consistent with the binomial momentum space acceptance [5]., Comment: 4 pages, 2 figures. Proceedings of the XXXth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2023), Houston, Texas, USA, September 3-9, 2023

It is well known that isospin symmetry is fulfilled to a good approximation in strong interactions, as confirmed in low-energy scattering experiments and in mass spectra of both light and heavy hadrons. In collisions of nuclei with an equal number of protons and neutrons, isospin symmetry imposes that the number of produced charged kaons should equal the number of neutral ones. The NA61/SHINE experiment at CERN recently reported an excess of charged over neutral kaon production in high-energy nucleus-nucleus collisions. Here, we argue that the measured charge-to-neutral kaon ratio of about 1.2 indicates an unexpectedly large violation of isospin symmetry. Using well-established models for hadron production, we demonstrate that known symmetry-breaking effects and the initial nuclei containing more neutrons than protons lead only to a small (few per cent) deviation from unity at high energies. Thus, they cannot explain the measurements. The significance of the isospin symmetry violation beyond the known effects is 5.5$ \cdot \sigma$ when errors quoted by the experiments are used and 8.1$ \cdot \sigma$ for the PDG-like scaled errors. New systematic, high-precision measurements and theoretical efforts are needed to establish the origin of the observed large isospin-symmetry breaking., Comment: 18 pages, 2 figures. Theoretical paper related to arXiv:2312.06572 [nucl-ex]

Molecular dynamics simulations are performed for a finite non-relativistic system of particles with Lennard-Jones potential. We study the effect of liquid-gas mixed phase on particle number fluctuations in coordinate subspace. A metastable region of the mixed phase, the so-called nucleation region, is analyzed in terms of a non-interacting cluster model. Large fluctuations due to spinodal decomposition are observed. They arise due to the interplay between the size of the acceptance region and that of the liquid phase. These effects are studied with a simple geometric model. The model results for the scaled variance of particle number distribution are compared with those obtained from the direct molecular dynamic simulations., Comment: 13 pages, 9 figures

We investigate the development of the directed, $v_1$, and elliptic flow, $v_2$, in heavy ion collisions in mid-central Au+Au reactions at $E_\mathrm{lab}=1.23 A$ GeV. We demonstrate that the elliptic flow of hot and dense matter is initially positive ($v_2>0$) due to the early pressure gradient. This positive $v_2$ transfers its momentum to the spectators, which leads to the creation of the directed flow $v_1$. In turn, the spectator shadowing of the in-plane expansion leads to a preferred decoupling of hadrons in the out-of-plane direction and results in a negative $v_2$ for the observable final state hadrons. We propose a measurement of $v_1-v_2$ flow correlations and of the elliptic flow of dileptons as methods to pin down this evolution pattern. The elliptic flow of the dileptons allows then to determine the early-state EoS more precisely, because it avoids the strong modifications of the momentum distribution due to shadowing seen in the protons. This opens the unique opportunity for the HADES and CBM collaborations to measure the Equation-of-State directly at 2-3 times nuclear saturation density., Comment: 8 pages, 5 figures

Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks., Comment: 5 pages, 1 fugure, contribution to the proceedings of Strangeness in Quark Matter 2022

Experiments motivated by predictions of quantum mechanics indicate non-trivial correlations between spacelike-separated measurements. The phenomenon is referred to as a violation of strong-locality and, after Einstein, called ghostly action at a distance. An intriguing and previously unasked question is how the evolution of an assembly of particles to equilibrium-state relates to strong-locality. More specifically, whether, with this respect, indistinguishable particles differ from distinguishable ones. To address the question, we introduce a Markov-chain based framework over a finite set of microstates. For the first time, we formulate conditions needed to obey the particle transport- and strong-locality for indistinguishable particles. Models which obey transport-locality and lead to equilibrium-state are considered. We show that it is possible to construct models obeying and violating strong-locality both for indistinguishable particles and for distinguishable ones. However, we find that only for distinguishable particles strongly-local evolution to equilibrium is possible without breaking the microstate-symmetry. This is the strongest symmetry one can impose and leads to the shortest equilibration time. We hope that the results presented here may provide a new perspective on a violation of strong-locality, and the developed framework will help in future studies. Specifically they may help to interpret results on high-energy nuclear collisions indicating a fast equilibration of indistinguishable particles.

We study fluctuations of particle number in the presence of critical point by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid in a periodic box. The numerical solution of the $N$-body problem naturally incorporates all correlations, exact conservation laws, and finite size effects, allowing us to study the fluctuation signatures of the critical point in a dynamical setup. We find that large fluctuations associated with the critical point are observed when measurements are performed in coordinate subspace, but, in the absence of collective flow and expansion, are essentially washed out when momentum cuts are imposed instead. We put our findings in the context of event-by-event fluctuations in heavy-ion collisions., Comment: 16 pages, 9 figures, to be published in Physical Review C

General formulas are presented for higher order cumulants of the conserved charge statistical fluctuations inside the mixed phase. As a particular example the van der Waals model in the grand canonical ensemble is used. The higher order measures of the conserved charge fluctuations up to the hyperkurtosis are calculated in a vicinity of the critical point (CP). The analysis includes both the mixed phase region and the pure phases on the phase diagram. It is shown that even-order fluctuation measures, e.g. scaled variance, kurtosis, and hyperkurtosis, have only positive values in the mixed phase, and go to infinity at the CP. For odd-order measures, such as skewness and hyperskewness, the regions of positive and negative values are found near the left and right binodals, respectively. The obtained results are discussed in a context of the event-by-event fluctuation measurements in heavy-ion collisions., Comment: 10 pages, 3 figures

We investigate simultaneous effects of finite system size and global charge conservation on thermal fluctuations in the vicinity of a critical point. For that we consider a finite interacting system which exchanges particles with a finite reservoir (thermostat), comprising a statistical ensemble that is distinct from the common canonical and grand canonical ensembles. As a particular example the van der Waals model is used. The global charge conservation effects strongly influence the cumulants of particle number distribution when the system size is comparable to that of the reservoir. If the system size is large enough to capture all the physics associated with the interactions, the global charge conservation effects can be accurately described and corrected for analytically, within a recently developed subensemble acceptance method. The finite size effects start to play a significant role when the correlation length grows large due to proximity of the critical point or when the system is small enough to be comparable to an eigenvolume of an individual particle. We discuss our results in the context of fluctuation measurements in heavy-ion collisions., Comment: 11 pages, 5 figures

The role of repulsive interactions in statistical systems of Bose particles is investigated. Three different phenomenological frameworks are considered: a mean field model, an excluded volume model, and a model with a medium dependent effective mass. All three models are tuned to yield similar equations of state, with only minor deviations from the ideal Bose gas at small chemical potentials. Our analysis indicates, however, that these models lead to qualitatively different results for the Bose-Einstein condensation phenomenon. We discuss the different aspects of this phenomenon, namely, an onset of the Bose-Einstein condensation, particle number fluctuations, and a behavior of the Bose condensate. The obtained results can be helpful for interpreting the lattice QCD data at small temperature and large isospin chemical potential and the data on multiple pion production in high energy nuclear collisions., Comment: 10 pages, 4 figures

The paper briefly presents history, status, and plans of the search for the critical structures - the onset of fireball, the onset of deconfinement, and the deconfinement critical point - in high energy nucleus-nucleus collisions. First, the basic ideas are introduced, the history of the observation of strongly interacting matter in heavy ion collisions is reviewed, and the path towards the quark-gluon plasma discovery is sketched. Then the status of the search for critical structures is discussed - the discovery of the onset of deconfinement, indications for the onset of fireball, and still inconclusive results concerning the deconfinement critical point. Finally, an attempt to formulate priorities for future measurements - charm quarks vs the onset of deconfinement and detailed study of the onset of fireball - closes the paper., Comment: 38 pages, 17 figures

We derive the relation between cumulants of a conserved charge measured in a subvolume of a thermal system and the corresponding grand-canonical susceptibilities, taking into account exact global conservation of that charge. The derivation is presented for an arbitrary equation of state, with the assumption that the subvolume is sufficiently large to be close to the thermodynamic limit. Our framework -- the subensemble acceptance method (SAM) -- quantifies the effect of global conservation laws and is an important step toward a direct comparison between cumulants of conserved charges measured in central heavy ion collisions and theoretical calculations of grand-canonical susceptibilities, such as lattice QCD. As an example, we apply our formalism to net-baryon fluctuations at vanishing baryon chemical potentials as encountered in collisions at the LHC and RHIC., Comment: 7 pages, 2 figures, published version

The QCD equation of state at finite temperature and densities of conserved charges is considered in the framework of a Hagedorn bag-like model, incorporating both the finite sizes of hadrons as well as their exponential mass spectrum. Augmented with non-zero masses of quarks and gluons in the bag spectrum, the model yields a fair quantitative description of lattice data on thermodynamic functions, the conserved charges susceptibilities, and Fourier coefficients of net-baryon density. Both at zero and finite baryon densities a broad crossover transition between hadronic and quark-gluon matter is observed. The model thus provides a thermodynamically consistent construction of a crossover equation of state for finite baryon number, electric charge and strangeness chemical potentials, which can be used in fluid dynamical simulations of heavy-ion collisions., Comment: 5 pages, 1 figure. Contribution to the Strangeness in Quark Matter 2019 (SQM 2019) conference proceedings

The binomial acceptance correction procedure is studied for particle number distributions detected in high energy reactions in finite regions of the momentum space. We present acceptance correction formulas for scaled variance, skewness, and kurtosis. Our considerations include various specific types of particles - positively or negatively charged, baryons and antibaryons - as well as conserved charges, namely, the net baryon number and electric charge. A simple model with effects of exact charge conservation, namely the Bessel distribution, is studied in some detail where effects of multi-particle correlations are present. The accuracy of the binomial filter is studied with UrQMD model simulations of inelastic proton-proton reactions. Binomial acceptance correction procedure works well when used inside a small region of phase space as well as for certain efficiency corrections, in particular for constructing net proton fluctuation from net baryon ones. Its performance is less accurate when applied to obtain UrQMD fluctuations in a finite rapidity window from fluctuations in the full $4\pi$ space., Comment: 17 pages, 7 figures

The nuclear liquid-gas transition at normal nuclear densities, $n \sim n_0 = 0.16$ fm$^{-3}$, and small temperatures, $T \sim 20$ MeV, has a large influence on analytic properties of the QCD grand-canonical thermodynamic potential. A classical van der Waals equation is used to determine these unexpected features due to dense cold matter qualitatively. The existence of the nuclear matter critical point results in thermodynamic branch points, which are located at complex chemical potential values, for $T > T_c \simeq 20$ MeV, and exhibit a moderate model dependence up to rather large temperatures $T \lesssim 100$ MeV. The behavior at higher temperatures is studied using the van der Waals hadron resonance gas (vdW-HRG) model. The baryon-baryon interactions have a decisive influence on the QCD thermodynamics close to $\mu_B = 0$. In particular, nuclear matter singularities limit the radius of convergence $r_{\mu_B/T}$ of the Taylor expansion in $\mu_B/T$, with $r_{\mu_B/T} \sim 2-3$ values at $T \sim 140-170$ MeV obtained in the vdW-HRG model., Comment: 11 pages, 5 figures

The production of backward nucleons, $N(180^\circ)$, at $180^\circ$ in the nuclear target rest frame in proton-nucleus ($\mathrm{p}+A$) collisions is studied. The backward nucleons appearing outside of the kinematically allowed range of proton-nucleon ($\mathrm{p}+N$) reactions are shown to be due to secondary reactions of heavy baryonic resonances produced inside the nucleus. Baryonic resonances $R$ created in primary $\mathrm{p}+N$ reactions can change their masses and momenta due to successive collisions $R+N\rightarrow R +N $ with other nuclear nucleons. Two distinct mechanisms and kinematic restrictions are studied: the reaction $R+N\rightarrow N(180^\circ)+N$ and the resonance decay $R\rightarrow N(180^\circ)+\pi$. Simulations of $\mathrm{p}+A$ collisions using the Ultra-relativistic Quantum Molecular Dynamics model support these mechanisms and are consistent with available data on proton backward production., Comment: 16 pages, 15 figures

Thermodynamic properties of systems with repulsive interactions, are considered in the grand canonical ensemble. The analytic structure of the excluded-volume model in the complex plane of the system chemical potential (fugacity) is elaborated, based on the fact that the pressure function can be given in terms of the Lambert W-function. Even though the excluded volume model has no phase transitions at real values of the chemical potential, it does exhibit a branch cut singularity in the complex plane, thus limiting the convergence range of the Taylor expansion in the chemical potential. Close similarities to analytic properties of the other models with repulsive interactions, such as a cluster expansion model, the mean-field model, and the ideal Fermi gas model, are pointed out. As an example, repulsive baryonic interactions in a hadron gas, with a focus on the fugacity/virial and Taylor expansion methods used in lattice QCD, are presented. The asymptotic behavior of the Fourier expansion coefficients in these various models suggests that the singular part of net baryonic density can to leading order be universally expressed in terms of polylogarithms., Comment: 11 pages, 4 figures

Thermodynamic functions, the (higher-order) fluctuations and correlations of conserved charges at $\mu_B = 0$, and the Fourier coefficients of net-baryon density at imaginary $\mu_B$, are considered in the framework of a Hagedorn bag-like model with a crossover transition. The qualitative behavior of these observables is found to be compatible with lattice QCD results. Fair quantitative description of the lattice data is obtained when quasiparticle-type quarks and gluons with non-zero masses are introduced into the bag spectrum. The equation of state of the model exhibits a smooth and wide crossover transition., Comment: 42 pages, 13 figures. Appendix added, accepted for publication in Physical Review C

The non-congruent liquid-gas phase transition (LGPT) in asymmetric nuclear matter is studied using the recently developed Quantum van der Waals model in the grand canonical ensemble. Different values of the electric-to-baryon charge ratio, $Q/B$, are considered. This non-congruent LGPT exhibits several features which are not present in the congruent LGPT of symmetric nuclear matter. These include a continuous phase transformation, a change in the location of the critical point, and the separation of the critical point and the endpoints. The effects which are associated with the non-congruent LGPT become negligible for the following cases: when $Q/B$ approaches its limiting values, $0.5$ or $0$, or if quantum statistical effects can be neglected. The latter situation is realized when the particle degeneracy attains large values, $g\gtrsim 10$., Comment: 24 pages, 20 figures

Different scenarios for modeling resonances in a thermal model description of hadron yields measured in heavy-ion collisions are explored: the zero-width approximation, the energy independent Breit-Wigner scheme, and the energy dependent Breit-Wigner (eBW) scheme. Application of the eBW scheme leads to a notable suppression in the proton yields, stemming mainly from a reduced feeddown from $\Delta$ resonances because of the threshold effects. A significantly improved agreement of thermal model with hadron yields measured in Pb-Pb collisions at $\sqrt{s_{_{NN}}} = 2.76$ TeV by the ALICE collaboration is obtained in the eBW scheme at $T \simeq 155$ MeV, indicating a possible resolution of the so-called 'proton anomaly'. The results obtained show that there are significant systematic uncertainties in the thermal model due to the modeling of broad resonances., Comment: 9 pages, 4 figures. Accepted for publication in Physical Review C, source code is available at https://github.com/vlvovch/1807.02079

The Monte Carlo (MC) procedure for sampling the hadron yields within hadron resonance gas (HRG) model is presented. The effects of excluded-volume due to the finite hadron eigenvolumes and of exact charge conservation within the canonical ensemble (CE) formulation are simultaneously taken into account with the help of the importance sampling technique combined with the rejection sampling. The MC procedure allows one to calculate arbitrary moments of the event-by-event hadron yields. Note that the CE formulation for the excluded-volume HRG has not been considered before. The MC simulations coincide with the known analytic results in the thermodynamic limit for the excluded-volume HRG in the grand canonical ensemble and for the CE of non-interacting particles. The MC procedure is applied to study the simultaneous excluded-volume and CE effects. These effects are considered within the full HRG to calculate the particle number fluctuations and to estimate the finite size effects., Comment: 29 pages, 5 figures. To be published in Physical Review C

The quantum van der Waals (QvdW) extension of the ideal hadron resonance gas (HRG) model which includes the attractive and repulsive interactions between baryons -- the QvdW-HRG model -- is applied to study the behavior of the baryon number related susceptibilities in the crossover temperature region. Inclusion of the QvdW interactions leads to a qualitatively different behavior of susceptibilities, in many cases resembling lattice QCD simulations. It is shown that for some observables, in particular for $\chi_{11}^{BQ} / \chi_2^B$, effects of the QvdW interactions essentially cancel out. It is found that the inclusion of the finite resonance widths leads to an improved description of $\chi_2^B$, but it also leads to a worse description of $\chi_{11}^{BQ} / \chi_2^B$, as compared to the lattice data. On the other hand, inclusion of the extra, unconfirmed baryons into the hadron list leads to a simultaneous improvement in the description of both observables., Comment: 4 pages, 1 figure. Proceedings for the Strangeness in Quark Matter (SQM2017) conference, Utrecht, the Netherlands, 10-15 July 2017

Event-by-event multiplicity fluctuations in nucleus-nucleus collisions are studied within the relativistic transport models: EPOS, PHSD, and UrQMD. As measures of particle number fluctuations we consider the scaled variances $\omega[X]$ for positive, negative, and all charged hadrons, and the strongly intensive quantities $\Delta[K^+,\pi^+],$ $\Sigma[K^+,\pi^+]$ for $K^+$ and $\pi^+$ yields. At the SPS energy range the fluctuation measures are calculated for proton-proton, Ar+Sc, and Pb+Pb collisions. Comparison with recent NA61/SHINE and older NA49 measurements of the multiplicity fluctuations is done. A validity of the model of independent sources, a role of the experimental acceptance, and the centrality selection procedure are studied.Event-by-event multiplicity fluctuations in nucleus-nucleus collisions are studied within the relativistic transport models: EPOS, PHSD, and UrQMD. As measures of particle number fluctuations we consider the scaled variances $\omega[X]$ for positive, negative, and all charged hadrons, and the strongly intensive quantities $\Delta[K^+,\pi^+],$ $\Sigma[K^+,\pi^+]$ for $K^+$ and $\pi^+$ yields. At the SPS energy range the fluctuation measures are calculated for proton-proton, Ar+Sc, and Pb+Pb collisions. Comparison with recent NA61/SHINE and older NA49 measurements of the multiplicity fluctuations is done. A validity of the model of independent sources, a role of the experimental acceptance, and the centrality selection procedure are studied., Comment: 9 pages, 10 figures, Accepted for publication in Journal of Physics G

The beam energy dependence of net baryon number susceptibilities is studied in the framework of the hadron resonance gas model with the attractive and repulsive van der Waals interactions between baryons. The collision energy dependences for the skewness $S\sigma$ and kurtosis $\kappa\sigma^2$ deviate significantly from the Poisson baseline and demonstrate the existence of rich structures at moderate collision energies. This behavior may result from the critical end point of the nuclear liquid-gas first order phase transition. In particular, $\kappa\sigma^2$ shows a non-monotonic energy dependence, and, in contrast to the standard scenario for the QCD critical point, it does not decrease at low collision energies. It is also found that the measurable net proton fluctuations differ significantly from the net baryon fluctuations when interactions between baryons cannot be neglected. The results are compared with the experimental net proton number fluctuations measured by the STAR collaboration., Comment: 6 pages, 3 figures. Extended discussion, error bands for calculations added, version accepted for publication in Physical Review C

The quantum mechanical Beth-Uhlenbeck (BU) approach for repulsive hard-core interactions between baryons is applied to the thermodynamics of a hadron gas. The second virial coefficient $a_2$ -- the "excluded volume" parameter -- calculated within the BU approach is found to be temperature dependent, and it differs dramatically from the classical excluded volume (EV) model result. At temperatures $T =100-200$ MeV, the widely used classical EV model underestimates the EV parameter for nucleons at a given value of the nucleon hard-core radius by large factors of 3-4. Previous studies, which employed the hard-core radii of hadrons as an input into the classical EV model, have to be re-evaluated using the appropriately rescaled EV parameters. The BU approach is used to model the repulsive baryonic interactions in the hadron resonance gas (HRG) model. Lattice data for the second and fourth order net baryon susceptibilities are described fairly well when the temperature dependent BU baryonic excluded volume parameter corresponds to nucleon hard-core radii of $r_c = 0.25-0.3$ fm. Role of the attractive baryonic interactions is also considered. It is argued that HRG model with a constant baryon-baryon EV parameter $v_{NN} \simeq 1$ fm$^3$ provides a simple yet efficient description of baryon-baryon interaction in the crossover temperature region., Comment: 22 pages, 4 figures, accepted for publication in Phys. Rev. C

The quantum statistical Clausius-based equation of state is used to describe the system of interacting nucleons. The interaction parameters $a$, $b$, and $c$ of the model are fixed by the empirically known nuclear ground state properties and nuclear incompressibility modulus. The model is generalized to describe the baryon-baryon interactions in the hadron resonance gas (HRG). The predictions of such a Clausius-HRG model are confronted with the lattice QCD data at zero and at small chemical potentials, and are also contrasted with the standard van der Waals approach. It is found that the behavior of the lattice QCD observables in a high-temperature hadron gas is sensitive to the nuclear matter properties. An improved description of the nuclear incompressibility factor correlates with an improved description of the lattice QCD data in the crossover transition region., Comment: 8 pages, 3 figures, submitted to EPJ A Topical Issue on "Frontiers in nuclear, heavy ion and strong field physics" dedicated to Walter Greiner

A generalization of the quantum van der Waals equation of state for a multi-component system in the grand canonical ensemble is proposed. The model includes quantum statistical effects and allows to specify the parameters characterizing repulsive and attractive forces for each pair of particle species. The model can be straightforwardly applied to the description of asymmetric nuclear matter and also for mixtures of interacting nucleons and nuclei. Applications of the model to the equation of state of an interacting hadron resonance gas are discussed., Comment: 20 pages, 5 figures, presentation improved, misprints corrected, version accepted for publication in Physical Review C

An extension of the ideal hadron resonance gas (HRG) model is constructed which includes the attractive and repulsive van der Waals (VDW) interactions between baryons. This VDW-HRG model yields the nuclear liquid-gas transition at low temperatures and high baryon densities. The VDW parameters $a$ and $b$ are fixed by the ground state properties of nuclear matter, and the temperature dependence of various thermodynamic observables at zero chemical potential are calculated within VDW-HRG model. Compared to the ideal HRG model, the inclusion of VDW interactions between baryons leads to a qualitatively different behavior of second and higher moments of fluctuations of conserved charges, in particular in the so-called crossover region $T\sim 140 - 190$ MeV. For many observables this behavior resembles closely the results obtained from lattice QCD simulations. This hadronic model also predicts nontrivial behavior of net-baryon fluctuations in the region of phase diagram probed by heavy-ion collision experiments. These results imply that VDW interactions play a crucial role in thermodynamics of hadron gas. Thus, the commonly performed comparisons of the ideal HRG model with the lattice and heavy-ion data may lead to misconceptions and misleading conclusions., Comment: 7 pages, 3 figures, final version published in Physical Review Letters, minor misprints corrected

The early stage of high multiplicity pp, pA and AA collider is represented by a nearly quarkless, hot, deconfined pure gluon plasma. According to pure Yang-Mills Lattice Gauge Theory, this hot pure glue matter undergoes, at a high temperature, $T_c = 270$ MeV, a first order phase transition into a confined Hagedorn-GlueBall fluid. These new scenario should be characterized by a suppression of high $p_T$ photons and dileptons, baryon suppression and enhanced strange meson production. We propose to observe this newly predicted class of events at LHC and RHIC., Comment: 7 pages, 6 figures

Several theoretical results concerning event-by-event fluctuations are discussed: (1) a role of the global conservation laws and concept of statistical ensembles; (2) strongly intensive measures are introduced; they give a possibility to study e-by-e fluctuations for physical systems with volume fluctuations; (3) identity method for chemical fluctuations in a case of incomplete particle identification; (4) the example of particle number fluctuations in a vicinity of the critical point., Comment: invited talk at the CPOD 2014 conference, 25 pages

A brief history of the observation of the onset of deconfinement - the beginning of the creation of quark gluon plasma in nucleus-nucleus collisions with increasing collision energy - is presented. It starts with the measurement of hadron mass spectrum and the Hagedorn's hypothesis of the limiting temperature of hadronic matter (the Hagedorn temperature). Then the conjecture that the Hagedorn temperature is the phase transition temperature was formulated with the crucial Hagedorn participation. It was confirmed by the observation of the onset of deconfinement in lead-lead collisions at the CERN SPS energies., Comment: arXiv admin note: substantial text overlap with arXiv:1201.0485

Strongly intensive measures $\Delta$ and $\Sigma$ are used to study event-by-event fluctuations of hadron multiplicities in nucleus-nucleus collisions. The effects of resonance decays are investigated within statistical model and relativistic transport model. Two specific examples are considered: resonance decays to two positively charged particles (e.g., $\Delta^{++}\rightarrow p+ \pi^+$) and to $\pi^+\pi^-$-pairs. (e.g., $\rho^0\rightarrow \pi^-+\pi^+$). It is shown that resonance abundances at the chemical freeze-out can be estimated by measuring the fluctuations of the number of stable hadrons. These model results are compared to the full hadron-resonance gas analysis within both the grand canonical and canonical ensemble. The ultra-relativistic quantum molecular dynamics (UrQMD) model of nucleus-nucleus collisions is used to illustrate the role of global charge conservation, centrality selection, and limited experimental acceptance., Comment: expanded discussions; added references; version published in J. Phys. G

The strongly intensive measures $\Delta[P_T,N]$ and $\Sigma[P_T,N]$ are used to study the event-by-event fluctuations of the transverse momentum $P_T$ and particle multiplicity $N$ in nucleus-nucleus collisions. A special normalization for these fluctuation measures ensures that they are dimensionless and yields a common scale required for a quantitative comparison of fluctuations. In this paper basic properties of the $\Delta[P_T, N]$ and $\Sigma[P_T, N]$ measures are tested within different phenomenological models using the Monte Carlo simulations (the so-called fast generators) and analytical solutions. The obtained results are helpful to elucidate the properties of $\Delta[P_T,N]$ and $\Sigma[P_T,N]$ measures., Comment: v2: text changes, Fig. 2 replaced (higher statistics), section with new simulations, Fig. 5 a) and b), added (centrality dependence of fluctuation measures), paper accepted by Phys. Rev. C