No Arabic abstract
Many models of heavy ion collisions employ relativistic hydrodynamics to describe the system evolution at high densities. The Cooper-Frye formula is applied in most of these models to turn the hydrodynamical fields into particles. However, the number of particles obtained from the Cooper-Frye formula is not always positive-definite. Physically negative contributions of the Cooper-Frye formula are particles that stream backwards into the hydrodynamical region. We quantify the Cooper-Frye negative contributions in a coarse-grained transport approach, which allows to compare them to the actual number of underlying particles crossing the transition hypersurface. It is found that the number of underlying inward crossings is much smaller than the one the Cooper-Frye formula gives under the assumption of equilibrium distribution functions. The magnitude of Cooper-Frye negative contributions is also investigated as a function of hadron mass, collision energy in the range $E_{rm lab} = 5-160A$ GeV, and collision centrality. The largest negative contributions we find are around 13% for the pion yield at midrapidity at $E_{rm lab} = 20A$ GeV collisions.
In most heavy ion collision simulations involving relativistic hydrodynamics, the Cooper-Frye formula is applied to transform the hydrodynamical fields to particles. In this article the so-called negative contributions in the Cooper-Frye formula are studied using a coarse-grained transport approach. The magnitude of negative contributions is investigated as a function of hadron mass, collision energy in the range of $E_{rm lab} = 5$--$160A$ GeV, collision centrality and the energy density transition criterion defining the hypersurface. The microscopic results are compared to negative contributions expected from hydrodynamical treatment assuming local thermal equilibrium. The main conclusion is that the number of actual microscopic particles flying inward is smaller than the negative contribution one would expect in an equilibrated scenario. The largest impact of negative contributions is found to be on the pion rapidity distribution at midrapidity in central collisions. For this case negative contributions in equilibrium constitute 8--13% of positive contributions depending on collision energy, but only 0.5--4% in cascade calculation. The dependence on the collision energy itself is found to be non-monotonous with a maximum at 10-20$A$ GeV.
The problem of spectra formation in hydrodynamic approach to A+A collisions is considered within the Boltzmann equations. It is shown analytically and illustrated by numerical calculations that the particle momentum spectra can be presented in the Cooper-Frye form despite freeze-out is not sharp and has the finite temporal width. The latter is equal to the inverse of the particle collision rate at points $(t_{sigma}({bf r},p),{bf r})$ of the maximal emission at a fixed momentum $p$. The set of these points forms the hypersurfaces $t_sigma({bf r},p)$ which strongly depend on the values of $p$ and typically do not enclose completely the initially dense matter. This is an important difference from the standard Cooper-Frye prescription (CFp), with a common freeze-out hypersurface for all $p$, that affects significantly the predicted spectra. Also, the well known problem of CFp as for negative contributions to the spectra from non-space-like parts of the freeze-out hypersurface is naturally eliminated in this improved prescription.
Microscopic transport approaches are the tool to describe the non-equilibrium evolution in low energy collisions as well as in the late dilute stages of high-energy collisions. Here, a newly developed hadronic transport approach, SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) is introduced. The overall bulk dynamics in low energy heavy ion collisions is shown including the excitation function of elliptic flow employing several equations of state. The implications of this new approach for dilepton production are discussed and preliminary results for afterburner calculations at the highest RHIC energy are presented and compared to previous UrQMD results. A detailed understanding of a hadron gas with vacuum properties is required to establish the baseline for the exploration of the transition to the quark-gluon plasma in heavy ion collisions at high net baryon densities.
We develop a new dynamical model for high energy heavy-ion collisions in the beam energy region of the highest net-baryon densities on the basis of non-equilibrium microscopic transport model JAM and macroscopic 3+1D hydrodynamics by utilizing a dynamical initialization method. In this model,dynamical fluidization of a system is controlled by the source terms of the hydrodynamic fields. In addition, time dependent core-corona separation of hot regions is implemented. We show that our new model describes multiplicities and mean transverse mass in heavy-ion collisions within a beam energy region of $3<sqrt{s_{NN}}<30$ GeV. Good agreement of the beam energy dependence of the $K^+/pi^+$ ratio is obtained, which is explained by the fact that a part of the system is not thermalized in our core-corona approach.
To explore the structure of the QCD phase diagram in high baryon density domain, several high-energy nuclear collision experiments in a wide range of beam energies are currently performed or planned using many accelerator facilities. In these experiments search for a first-order phase transition and the QCD critical point is one of the most important topics. To find the signature of the phase transition, experimental data should be compared to appropriate dynamical models which quantitatively describe the process of the collisions. In this study we develop a new dynamical model on the basis of the non-equilibrium hadronic transport model JAM and 3+1D hydrodynamics. We show that the new model reproduce well the experimental beam-energy dependence of hadron yields and particle ratio by the partial thermalization of the system in our core-corona approach.