Search

We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions. BHAC-QGP is based on the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes through the solution of the equations of general-relativistic magnetohydrodynamics. Like the mother code, BHAC-QGP uses Adaptive Mesh Refinement (AMR), which allows for a dynamic adjustment of the resolution in regions of the computational domain where a particularly high accuracy is needed. We here discuss a number of applications of BHAC-QGP to Au-Au collisions at Relativistic Heavy-Ion Collider (RHIC) energies and show that the code is able to reproduce results of other simulations of these scenarios, but with much higher accuracy., Comment: 17 pages, 18 figures

We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions in the presence of electromagnetic fields. It is derived from the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes in a general-relativistic magnetohydrodynamical description. As the original Black Hole Accretion Code, BHAC-QGP benefits from the use of Adaptive Mesh Refinement (AMR), which allows us to dynamically adjust the resolution where necessary, and makes use of time-dependent Milne coordinates and the ultrarelativistic equation of state, $P = e/3$. We demonstrate that BHAC-QGP accurately passes a number of systematic and rigorous tests., Comment: 18 pages, 7 figures

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS., Comment: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

The polarization of the $\Lambda$ hyperon has become an important probe of the Quark-Gluon Plasma produced in relativistic heavy-ion collisions. Recently, it has been found that polarization receives a substantial contribution from a local equilibrium term proportional to the symmetric derivative of the four-temperature vector, the thermal shear tensor. We show that, at very high energies, this term can restore the agreement between the experimental measurements and the predictions of the hydrodynamic model, provided that the hadronization hypersurface is isothermal. We review the theoretical derivation of this new term, discuss numerical computations at RHIC and LHC energies, and compare them with the experimental data. We also present the effect of feed-down corrections., Comment: 6 pages, 4 figures, presented at the 20th international conference on Strangeness in Quark Matter (SQM 2022) v2: Added ancillary files for the figures. Figure 2 has been corrected

To assess the degree of equilibration of the matter created in heavy-ion reactions at low to intermediate beam energies, a hadronic transport approach (SMASH) is employed. By using a coarse-graining method, we compute the energy momentum tensor of the system at fixed time steps and evaluate the degree of isotropy of the diagonal terms and the relative magnitude of the off-diagonal terms. This study focuses mostly on Au+Au collisions in the energy range $\sqrt{s_{\rm NN}}$ = 2.4-7.7 GeV, but central collisions of lighter ions like C+C, Ar+KCl and Ag+Ag are considered as well. We find that the conditions concerning local equilibration for a hydrodynamic description are reasonably satisfied in a large portion of the system for a significant amount of time (several fm/c) when considering the average evolution of many events, yet they are rarely fulfilled on an event by event basis. This is relevant for the application of hybrid approaches at low beam energies as they are or will be reached by the HADES experiment at GSI, the future CBM experiment at FAIR as well as the beam energy scan program at RHIC., Comment: Ancillary files with plots data and workflow provided. Version 2 corresponding to version accepted for publication by EPJC

The freeze-out process in heavy ion collisions is driven by the competition between the scattering rate and the expansion rate of the matter. We analyse the expansion rate $\Theta$ (often called Hubble flow) in relativistic heavy ion collisions in the FAIR and RHIC-BES energy regimes and compare it to the scattering rate $\Gamma$ using the UrQMD transport model. We observe that the time evolution of the system is clearly separated into a compression phase and an expansion phase with time dependent $\Theta_\parallel$ and $\Theta_\perp$. The calculated values of the Hubble expansion at kinetic decoupling are in line with previous simple estimates by statistical hadronization models with a Siemens-Rasmussen type emission source. However, the actual shape of the expanding matter is, as expected, found to be between a spherically symmetric and a purely longitudinal expansion. We confirm for the first time in a microscopic simulation that the decoupling hypersurface is indeed determined by the competition of the expansion rate $\Theta$ and the scattering rate $\Gamma$ as suggested previously. This suggests that, in the range of collision energies explored in this study, simple iso-thermal/iso-energy density criteria that are often used in hybrid models to couple hydrodynamic and transport simulations may not capture the true decoupling hyper-surface.

The HADES experiment at GSI has recently provided data on the flow coefficients $v_1,...,v_4$ for protons in Au+Au reactions at $E_{\rm lab} = 1.23$~$A$GeV (or $\sqrt{s_\mathrm{NN}}=2.4$ GeV). This data allows to estimate the shear viscosity over entropy ratio, $\eta/s$ at low energies via a coarse graining analysis of the UrQMD transport simulations of the flow harmonics in comparison to the experimental data. By this we can provide for the first time an estimate of $\eta/s\approx0.65\pm0.15$ (or $(8\pm2)\,(4\pi)^{-1}$) at such low energies., Comment: 6 pages, 5 figures

We introduce a novel approach based on elastic and inelastic scattering rates to extract the hyper-surface of the chemical freeze-out from a hadronic transport model in the energy range from E$_\mathrm{lab}=1.23$ AGeV to $\sqrt{s_\mathrm{NN}}=62.4$ GeV. For this study, the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model combined with a coarse-graining method is employed. The chemical freeze-out distribution is reconstructed from the pions through several decay and re-formation chains involving resonances and taking into account inelastic, pseudo-elastic and string excitation reactions. The extracted average temperature and baryon chemical potential are then compared to statistical model analysis. Finally we investigate various freeze-out criteria suggested in the literature. We confirm within this microscopic dynamical simulation, that the chemical freeze-out at all energies coincides with $\langle E\rangle/\langle N\rangle\approx1$ GeV, while other criteria, like $s/T^3=7$ and $n_\mathrm{B}+n_\mathrm{\bar{B}}\approx0.12$ fm$^{-3}$ are limited to higher collision energies., Comment: 9 pages, 10 figures

Using the UrQMD/coarse graining approach we explore the kinetic freeze-out stage in central Au + Au collisions at various energies. These studies allow us to obtain detailed information on the thermodynamic properties (e.g. temperature and chemical potential) of the system during the kinetic decoupling stage. We explore five relevant collision energies in detail, ranging from $\sqrt{s_{NN}}=2.4\,\mathrm{GeV}$ (GSI-SIS) to $\sqrt{s_{NN}}=200\,\mathrm{GeV}$ (RHIC). By adopting a standard Hadron Resonance Gas equation of state, we determine the average temperature $\langle T \rangle$ and the average baryon chemical potential $\langle\mu_{\mathrm{B}}\rangle$ on the space-time hyper-surface of last interaction. The results highlight the nature of the kinetic freeze-out as a continuous process. This differential decoupling is an important aspect often missed when summarizing data as single points in the phase diagram as e.g. done in Blast-Wave fits. We compare the key properties of the system derived by using our approach with other models and we briefly review similarities and differences., Comment: Version accepted by JPG

At the earliest times after a heavy-ion collision, the magnetic field created by the spectator nucleons will generate an extremely strong, albeit rapidly decreasing in time, magnetic field. The impact of this magnetic field may have detectable consequences, and is believed to drive anomalous transport effects like the Chiral Magnetic Effect (CME). We detail an exploratory study on the effects of a dynamical magnetic field on the hydrodynamic medium created in the collisions of two ultrarelativistic heavy-ions, using the framework of numerical ideal MagnetoHydroDynamics (MHD) with the ECHO-QGP code. In this study, we consider a magnetic field captured in a conducting medium, where the conductivity can receive contributions from the electromagnetic conductivity $\sigma$ and the chiral magnetic conductivity $\sigma_{\chi}$. We first study the elliptic flow of pions, which we show is relatively unchanged by the introduction of a magnetic field. However, by increasing the magnitude of the magnetic field, we find evidence for an enhancement of the elliptic flow in peripheral collisions. Next, we explore the impact of the chiral magnetic conductivity on electric charges produced at the edges of the fireball. This initial $\sigma_\chi$ can be understood as a long-wavelength effective description of chiral fermion production. We then demonstrate that this chiral charge, when transported by the MHD medium, produces a charge dipole perpendicular to the reaction plane which extends a few units in rapidity. Assuming charge conservation at the freeze-out surface, we show that the produced charge imbalance can have measurable effects on some experimental observables, like $v_1$ or $\langle \sin \phi \rangle$. This demonstrates the ability of a MHD fluid to transport the signature of the initial chiral magnetic fields to late times., Comment: Version accepted by EPJC with minor modifications. Computation of charge dependent pion spectra partially changed. Added a few additional explanations. Other minor changes

We present a study of the elliptic flow and $R_{AA}$ of D and anti-D mesons in Au+Au collisions at FAIR energies. We propagate the charm quarks and the D mesons following a previously applied Langevin dynamics. The evolution of the background medium is modeled in two different ways: (I) we use the UrQMD hydrodynamics + Boltzmann transport hybrid approach including a phase transition to QGP and (II) with the coarse-graining approach employing also an equation of state with QGP. The latter approach has previously been used to describe di-lepton data at various energies very successfully. This comparison allows us to explore the effects of partial thermalization and viscous effects on the charm propagation. We explore the centrality dependencies of the collisions, the variation of the decoupling temperature and various hadronization parameters. We find that the initial partonic phase is responsible for the creation of most of the D mesons elliptic flow and that the subsequent hadronic interactions seem to play only a minor role. This indicates that D mesons elliptic flow is a smoking gun for a partonic phase at FAIR energies. However, the results suggest that the magnitude and the details of the elliptic flow strongly depend on the dynamics of the medium and on the hadronization procedure, which is related to the medium properties as well. Therefore, even at FAIR energies the charm quark might constitute a very useful tool to probe the Quark-Gluon Plasma and investigate its physics.

We present an improved version of the ECHO-QGP numerical code, which self-consistently includes for the first time the effects of electromagnetic fields within the framework of relativistic magnetohydrodynamics (RMHD). We discuss results of its application in relativistic heavy-ion collisions in the limit of infinite electrical conductivity of the plasma. After reviewing the relevant covariant $3\!+\!1$ formalisms, we illustrate the implementation of the evolution equations in the code and show the results of several tests aimed at assessing the accuracy and robustness of the implementation. After providing some estimates of the magnetic fields arising in non-central high-energy nuclear collisions, we perform full RMHD simulations of the evolution of the Quark-Gluon Plasma in the presence of electromagnetic fields and discuss the results. In our ideal RMHD setup we find that the magnetic field developing in non-central collisions does not significantly modify the elliptic-flow of the final hadrons. However, since there are uncertainties in the description of the pre-equilibrium phase and also in the properties of the medium, a more extensive survey of the possible initial conditions as well as the inclusion of dissipative effects are indeed necessary to validate this preliminary result., Comment: 19 pages, 14 figures, vers. 3: final version, with minor changes in the text (introduction and conclusions)

Initial fluctuations in hydrodynamic fields such as energy density or flow velocity give access to understanding initial state and equilibration physics as well as thermodynamic and transport properties. We provide evidence that the fluid dynamic propagation of fluctuations of realistic size can be based on a background-fluctuation splitting and a systematic perturbative expansion in the fluctuating fields. Initial conditions are characterized by a Bessel-Fourier expansion for single events, event-by-event correlations and probability distributions. The evolution equations can be solved order-by-order in the expansion which allows to study the fluid dynamical propagation of single modes, the study of interaction effects between modes, the determination of the associated particle spectra and the generalization of the whole program to event-by-event correlations and distributions., Comment: poceedings of the XXIV Quark Matter conference (2014)

We present a numerical code modeling the evolution of the medium formed in relativistic heavy ion collisions, ECHO-QGP. The code solves relativistic hydrodynamics in $(3+1)-$D, with dissipative terms included within the framework of Israel-Stewart theory; it can work both in Minkowskian and in Bjorken coordinates. Initial conditions are provided through an implementation of the Glauber model (both Optical and Monte Carlo), while freezeout and particle generation are based on the Cooper-Frye prescription. The code is validated against several test problems and shows remarkable stability and accuracy with the combination of a conservative (shock-capturing) approach and the high-order methods employed. In particular it beautifully agrees with the semi-analytic solution known as Gubser flow, both in the ideal and in the viscous Israel-Stewart case, up to very large times and without any ad hoc tuning of the algorithm., Comment: 4 pages, 6 figures, XXIV Quark Matter conference proceedings (2014)

High angular resolution spectra obtained with the Hubble Space Telescope Imaging Spectrograph (HST/STIS) provide rich morphological and kinematical information about the stellar jet phenomenon, which allows us to test theoretical models efficiently. In this work, numerical simulations of stellar jets in the propagation region are executed with the PLUTO code, by adopting inflow conditions that arise from former numerical simulations of magnetized outflows, accelerated by the disk-wind mechanism in the launching region. By matching the two regions, information about the magneto-centrifugal accelerating mechanism underlying a given astrophysical object can be extrapolated by comparing synthetic and observed position-velocity diagrams (PVDs). We show that quite different jets, like those from the young T Tauri stars DG-Tau and RW-Aur, may originate from the same disk-wind model for different configurations of the magnetic field at the disk surface. This result supports the idea that all the observed jets may be generated by the same mechanism., Comment: 15 pages, 18 figures, accepted for publication by A&A

We provide evidence from full numerical solutions that the hydrodynamical evolution of initial density fluctuations in heavy ion collisions can be understood order-by-order in a perturbative series in deviations from a smooth and azimuthally symmetric background solution. To leading linear order, modes with different azimuthal wave numbers do not mix. Quadratic and higher order corrections are small and can be understood as overtones with corresponding wave numbers., Comment: 8 pages, 4 figures