Source author record

Jan Steinheimer

Jan Steinheimer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2023arXiv

Spinodal enhancement of fluctuations in nucleus-nucleus collisions

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].

preprint2022arXiv

3-dimensional flow analysis: A novel tool to study the collision geometry and the Equation-of-State

We propose to extend the commonly known flow analysis in the transverse $p_x$-$p_y$ plane to novel flow coefficients based on the angular distribution in the $p_x$-$p_z$ and $p_y$-$p_z$ planes. The new flow coefficients, called $u_n$ and $w_n$ (in addition to $v_n$), turn out to be also highly sensitive to the nuclear Equation-of-State and can be used to explore the EoS in more detail than is possible using only $v_n$. As an example to quantify the effect of the EoS, the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model is used to investigate 20-30\% central Au+Au collisions at E$_\mathrm{lab}=1.23~A$GeV.

preprint2022arXiv

A Chiral Mean-Field Equation-of-State in UrQMD: Effects on the Heavy Ion Compression Stage

It is shown that the initial compression in central heavy ion collisions at beam energies of $E_\mathrm{lab}=1-10A$~GeV depends dominantly on the underlying equation of state and only marginally on the model used for the dynamical description. To do so, a procedure to incorporate any equation of state in the UrQMD transport model is introduced. In particular we compare the baryon density, temperature and pressure evolution as well as produced entropy in a relativistic ideal hydrodynamics approach and the UrQMD transport model, where the same equation of state is used in both approaches. Not only is the compression similar if the same equation of state is used in either dynamical model, but it also strongly depends on the actual equation of state. These results indicate that the equation of state can be studied with observables which are sensitive to the initial compression phase and maximum compression achieved in heavy ion collisions at these beam energies.

preprint2022arXiv

Coalescence, the thermal model and multi-fragmentation: The energy and volume dependence of light nuclei production in heavy ion collisions

We present results of a phase space coalescence approach within the UrQMD transport and -hybrid model for a very wide range of beam energies from SIS to LHC. The coalescence model is able to qualitatively describe the whole range of experimental data with a fixed set of parameters. Some systematic deviations are observed for very low beam energies where the role of feed down from heavier nuclei and multi-fragmentation becomes relevant. The coalescence results are mostly very close to the thermal model fits. However, both the coalescence approach as well as thermal fits are struggling to simultaneously describe the triton multiplicities measured with the STAR and ALICE experiment. The double ratio of $tp/d^2$, in the coalescence approach, is found to be essentially energy and centrality independent for collisions of heavy nuclei at beam energies of $\mathrm{E_{lab}}> 10 A$ GeV. On the other hand the clear scaling of the $d/p^2$ and $t/p^3$ ratios with the systems volume is broken for peripheral collisions, where a canonical treatment and finite size effects become more important.

preprint2022arXiv

Constraining baryon annihilation in the hadronic phase of heavy-ion collisions via event-by-event fluctuations

We point out that the variance of net-baryon distribution normalized by the Skellam distribution baseline, $κ_2[B-\bar{B}]/\langle B+\bar{B}\rangle$, is sensitive to the possible modification of (anti)baryon yields due to $B\bar{B}$ annihilation in the hadronic phase. The corresponding measurements can thus place stringent limits on the magnitude of the $B\bar{B}$ annihilation and its inverse reaction. We perform Monte Carlo simulations of the hadronic phase in Pb-Pb collisions at the LHC via the recently developed subensemble sampler + UrQMD afterburner and show that the effect survives in net-proton fluctuations, which are directly accessible experimentally. The available experimental data of the ALICE Collaboration on net-proton fluctuations disfavors a notable suppression of (anti)baryon yields in $B\bar{B}$ annihilations predicted by the present version of UrQMD if only global baryon conservation is incorporated. On the other hand, the annihilations improve the data description when local baryon conservation is imposed. The two effects can be disentangled by measuring $κ_2[B+\bar{B}]/\langle B+\bar{B}\rangle$, which at the LHC is notably suppressed by annihilations but virtually unaffected by baryon number conservation.

preprint2022arXiv

Correcting the $B_A$ coalescence factor at GSI-HADES and RHIC-BES energies

We investigate the coalescence factors $B_2$ and $B_3$ at low collision energies ($\sqrt{s_\mathrm{NN}}<6$ GeV) with special focus on the HADES and RHIC-BES experiments. It is shown that, in order to properly interpret the coalescence factors $B_A$, two important corrections are necessary: I) $B_2$ has to be calculated using the proton $\times$ neutron yields in the denominator, instead of the square of the proton yield, and II) the primordial proton (neutron) densities have to be used for the normalization and not the final state (free) protons (neutrons). Both effects lead to a drastic reduction of $B_2$ and $B_3$ at low energies. This reduction decreases the discrepancy between the volumes extracted from HBT measurements and the volumes extracted from the coalescence factor ($V\propto 1/B_2$). While at HADES and low RHIC-BES energies these corrections are substantial, they become irrelevant above $\sqrt{s_\mathrm{NN}}>6$ GeV. The proposed correction method is model independent and is only based on the measurement of protons, clusters and charged pions.

preprint2022arXiv

Deuteron Production in Ultra-Relativistic Heavy-Ion Collisions: A Comparison of the Coalescence and the Minimum Spanning Tree Procedure

The formation of deuterons in heavy-ion collisions at relativistic energies is investigated by employing two recently advanced models -- the Minimum Spanning Tree (MST) method and the coalescence model by embedding them in the PHQMD and the UrQMD transport approaches. While the coalescence mechanism combines nucleons into deuterons at the kinetic freeze-out hypersurface, the MST identifies the clusters during the different stages of time evolution. We find that both clustering procedures give very similar results for the deuteron observables in the UrQMD as well as in the PHQMD environment. Moreover, the results agree well with the experimental data on deuteron production in Pb+Pb collisions at $\sqrt{s_{NN}} = 8.8$ GeV (selected for the comparison of the methods and models in this study). A detailed investigation shows that the coordinate space distribution of the produced deuterons differs from that of the free nucleons and other hadrons. Thus, deuterons are not destroyed by additional rescattering.

preprint2022arXiv

Dilepton Signature of a First-Order Phase Transition

The search for a first-order phase transition in strongly interacting matter is one of the major objectives in the exploration of the phase diagram of Quantum Chromodynamics (QCD). In the present work we investigate dilepton radiation from the hot and dense fireballs created in Au-Au collisions at projectile energies of 1-2 $A$GeV for potential signatures of a first-order transition. Toward this end, we employ a hydrodynamic simulation with two different equations of state, with and without a phase transition. The latter is constrained by susceptibilities at vanishing chemical potential from lattice-QCD as well as neutron star properties, while the former is implemented via modification of the mean-fields in the quark phase. We find that the latent heat involved in the first-order transition leads to a substantial increase in the low-mass thermal emission signal, by about a factor of two above the cross-over scenario.

preprint2022arXiv

Photons from relativistic nuclear collisions

Collisions of atomic nuclei at relativistic velocities allow to recreate the conditions encountered in neutron stars or in the early universe micro-seconds after the Big Bang. These reactions are performed in today's largest accelerator facilities, e.g. at CERN in Geneva, at the Relativistic Heavy Ion Collider at Brookhaven, NY or in the planned FAIR facility in Darmstadt Germany. During such a collision the matter is heated up to hundreds of MeV (billions of degrees) and compressed to densities of $3-10$ times the density inside ordinary atomic nuclei (i.e. $10^{17}-10^{18}$ kg/m$^3$). Usually these collisions are studied via the measurement of a multitude of strongly interacting particles, called hadrons, that are emitted at the end of the collision. However, also some photons are created. These photons are of special interest as they allow to look into the early stage of the collisions, because they are only very weakly (namely only electro-magnetically) interacting with the hadrons of the created fireball. This paper elucidates the physics of the photons and what can be learned from them.

preprint2022arXiv

Sunda-arc seismicity: continuing increase of high-magnitude earthquakes since 2004

Spatial and temporal data for earthquakes with magnitude M greater than or equal to 6.5 can provide crucial information about the seismic history and potential for large earthquakes in a region. We analyzed approximately 313,500 events that occurred in the Sunda-arc region during the last 56 years, from 1964 to 2020, reported by the International Seismological Center. We report a persistent increase in the annual number of the events with mb greater than or equal to 6.5. We tested this increase against the null hypothesis and discarded the possibility of the increase being due to random groupings. The trend given by Auto-Regressive Integrated Moving Average suggests continuing increase of such large-magnitude events in the region during the next decade. At the same time, the computed Gutenberg Richter b value shows anomalies that can be related to the occurrence of the mega 2004 Sumatra earthquake, and to possible state of high tectonic stress in the eastern parts of the region.

preprint2022arXiv

The Role of the Hadron-Quark Phase Transition in Core-Collapse Supernovae

The hadron-quark phase transition in quantum chromodyanmics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F\_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\texttt{-}100\,\text{M}_\odot$. We find explosions only for two low-compactness models ($14 \text{M}_\odot$ and $16\,\text{M}_\odot$) with the DD2F\_SF EoS, both with low explosion energies of $\mathord{\sim}10^{50}\,\mathrm{erg}$. These weak explosions are characterised by a neutrino signal with several mini-bursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as $^{90}\mathrm{Zr}$ for the $14\,\text{M}_\odot$ zero-metallicity model and $^{94}\mathrm{Zr}$ for the $16\,\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F\_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the proto-compact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to $4k_\text{B}/\text{baryon}$ and may have a distinctive gravitational wave signature, also for a second-order phase transition.

preprint2021arXiv

Comparison of heavy ion transport simulations: Ag+Ag collisions at Elab = 1.58 AGeV

We compare the microscopic transport models UrQMD, PHSD, PHQMD, and SMASH to make predictions for the upcoming Ag+Ag data at $E_\mathrm{lab}=1.58$~$A$GeV ($\sqrt{s_\mathrm{NN}}=2.55$~GeV) by the HADES collaboration. We study multiplicities, spectra and effective source temperatures of protons, $π^{\pm,0}$, $K^\pm$, the $η$, $Λ+Σ^0$ and the $Ξ^-$ within these models. Despite variations in the detailed implementation of the dynamics in the different models, the employed transport approaches all show consistent multiplicities of the bulk of investigated hadrons. The main differences are in the $Ξ^-$ production, which is treated differently between UrQMD/SMASH on one side employing high mass resonance states with explicit decays to $\mathrm{Resonance}\rightarrow Ξ+K+K$ in contrast to PHSD/PHQMD which account only non-resonant $Ξ$ production channels. A comparison of the spectra, summarized by effective source temperatures, shows that all models provide similar source temperatures around $T_\mathrm{source}=80-95$~MeV, and show substantial radial flow on the order of $\langle v_T\rangle=0.22c-0.3c$ even for such a small system.

preprint2021arXiv

Proton correlations and apparent intermittency in the UrQMD model with hadronic potentials

It is shown that the inclusion of hadronic interactions, and in particular nuclear potentials, in simulations of heavy ion collisions at the SPS energy range can lead to obvious correlations of protons. These correlations contribute significantly to an intermittency analysis as performed at the NA61 experiment. The beam energy and system size dependence is studied by comparing the resulting intermittency index for heavy ion collisions of different nuclei at beam energies of $40A$, $80A$ and $150A$ GeV. The resulting intermittency index from our simulations is similar to the reported values of the NA61 collaboration, if nuclear interactions are included. The observed apparent intermittency signal is the result of the correlated proton pairs with small relative transverse momentum $Δp_{t}$, which would be enhanced by hadronic potentials, and this correlation between the protons is slightly influenced by the coalescence parameters and the relative invariant four-momentum $q_{inv}$ cut.

preprint2020arXiv

Critical point fluctuations: Finite size and global charge conservation effects

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.

preprint2020arXiv

Elliptic flow splitting between protons and antiprotons from hadronic potentials

The difference in elliptic flow $v_{2}$ between protons and antiprotons, produced in $^{197}\text{Au}+^{197}\text{Au}$ collisions at center-of-mass energies $\sqrt{s_{NN}}=5-12~\text{GeV}$, is studied within a modified version of the ultrarelativistic quantum molecular dynamics (UrQMD) model. Two different model scenarios are compared: the cascade mode and the mean field mode which includes potential interactions for both formed and pre-formed hadrons. The model results for the elliptic flow of protons and the relative $v_{2}$ difference between protons and antiprotons obtained from the mean field mode agree with the available experimental data, while the $v_{2}$ difference is near zero for the cascade mode. Our results show that the elliptic flow splitting, observed for particles and antiparticles, can be explained by the inclusion of proper hadronic interactions. In addition, the difference in $v_{2}$ between protons and antiprotons depends on the centrality and the rapidity window. With smaller centrality and/or rapidity acceptance, the observed elliptic flow splitting is more sensitive to the beam energy, indicating a strong net baryon density dependence of the effect. We propose to confirm this splitting at the upcoming experiments from Beam Energy Scan (BES) Phase-\Rmnum{2} at Relativistic Heavy Ion Collider (RHIC), the Compressed Baryonic Matter (CBM) at Facility for Antiproton and Ion Research (FAIR), High Intensity heavy ion Accelerator Facility (HIAF) and Nuclotron-based Ion Collider fAcility (NICA).

preprint2020arXiv

Identifying the nature of the QCD transition in relativistic collision of heavy nuclei with deep learning

Using deep convolutional neural network (CNN), the nature of the QCD transition can be identified from the final-state pion spectra from hybrid model simulations of heavy-ion collisions that combines a viscous hydrodynamic model with a hadronic cascade "after-burner". Two different types of equations of state (EoS) of the medium are used in the hydrodynamic evolution. The resulting spectra in transverse momentum and azimuthal angle are used as the input data to train the neural network to distinguish different EoS. Different scenarios for the input data are studied and compared in a systematic way. A clear hierarchy is observed in the prediction accuracy when using the event-by-event, cascade-coarse-grained and event-fine-averaged spectra as input for the network, which are about 80%, 90% and 99%, respectively. A comparison with the prediction performance by deep neural network (DNN) with only the normalized pion transverse momentum spectra is also made. High-level features of pion spectra captured by a carefully-trained neural network were found to be able to distinguish the nature of the QCD transition even in a simulation scenario which is close to the experiments.

preprint2020arXiv

Influence of centrality definition and detector efficiency on the net-proton kurtosis

We study the influence of the centrality definition and detector efficiency on the net-proton kurtosis for minimum bias Au+Au collisions at a beam energy of $\sqrt{s_{\mathrm{NN}}}= 7.7$ GeV by using the UrQMD model. We find that different ways of defining the centrality lead to different cumulant ratios. Moreover, we demonstrate that the kurtosis is suppressed for central collisions when a wider transverse momentum acceptance is used. Finally, the influence of a detector efficiency on the measured cumulant ratios is estimated.

preprint2020arXiv

QCD equation of state at vanishing and high baryon density: Chiral Mean Field model

The thermodynamic properties of high temperature and high density QCD-matter are studied using the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The CMF model provides a proper description of lattice QCD data, heavy-ions physics, and static neutron stars. The behavior of lines of constant pressure with increase of baryon density is discussed. The rapid change of pressure behavior at $μ_B/T\approx3$ suggests a strong contribution of baryons to thermodynamic properties at this region. The position of this region is very close to the radius of convergence for a Taylor expansion of the QCD pressure. The role of mesons and unstable hadrons in the hydrodynamic expansion of strongly interacting matter is also discussed.

preprint2020arXiv

The QCD phase diagram and statistics friendly distributions

The preliminary STAR data for proton cumulants for central collisions at \sqrt{s}=7.7 GeV are consistent with a two-component proton multiplicity distribution. We show that this two-component distribution is statistics friendly in that factorial cumulants of surprisingly high orders may be extracted with a relatively small number of events. As a consequence the two-component model can be tested and verified right now with the presently available STAR data from the first phase of the RHIC beam energy scan.

preprint2019arXiv

Equation of state for hot QCD and compact stars from a mean field approach

The thermodynamic properties of high temperature and high density QCD-matter are explored within the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The quark sector of the CMF model is tuned to describe the $μ_B=0$ thermodynamics data of lattice QCD. The resulting lines of constant physical variables as well as the baryon number susceptibilities are studied in some detail in the temperature/chemical potential plane. The CMF model predicts three consecutive transitions, the nuclear first-order liquid-vapor phase transition, chiral symmetry restoration, and the cross-over transition to a quark-dominated phase. All three phenomena are cross-over, for most of the $T-μ_B$-plane. The deviations from the free ideal hadron gas baseline at $μ_B=0$ and $T\approx 100-200$ MeV can be attributed to remnants of the liquid-vapor first order phase transition in nuclear matter. The chiral crossing transition determines the baryon fluctuations at much higher $μ_B\approx1.5$ GeV, and at even higher baryon densities $μ_B\approx2.4$ GeV, the behavior of fluctuations is controlled by the deconfinement cross-over. The CMF model also describe well the static properties of high $μ_B$ neutron stars as well as the new neutron star merger observations. The effective EoS presented here describes simultaneously lattice QCD results at $μ_B=0$, as well as observed physical phenomena (nuclear matter and neutron star matter) at $T\cong0$ and high densities, $μ_B>1$ GeV.

preprint2019arXiv

First, second, third and fourth flow harmonics of deuterons and protons in Au+Au reactions at 1.23 A GeV

We explore the directed, elliptic, triangular and quadrangular flow of deuterons in Au+Au reactions at a beam energy of 1.23 AGeV within the UrQMD approach. These investigations are of direct relevance for the HADES experiment at GSI that has recently presented first data on the flow of light clusters in Au+Au collisions at 1.23 AGeV. To address the deuteron flow, UrQMD has been extended to include deuteron formation by coalescence. We find that this ansatz provides a very good description of the measured deuteron flow data, if a hard equation of state is used for the simulation. In addition we show that light cluster formation has a sizable impact on the proton flow and has to be taken into account to obtain reliable results in the forward/backward region. Based on the observed scaling of the flow, which is a natural result of coalescence, we conclude that deuteron production at GSI energies is a final state recombination effect. Finally, we also discuss the scaling relations of the higher order flow components up to $v_4$. We show that $v_3 \sim v_1v_2$ and $v_4 \sim v_2^2$ as function of transverse momentum and that the integrated $v_2^2 \sim v_4$ over the investigated energy range from $E_{lab}$=0.1 AGeV to 40 AGeV.

preprint2016arXiv

Formation of deuterons by coalescence: Consequences on the deuteron number fluctuations

Two scenarios for cluster production have since long been discussed in the literature: i) direct emission of the clusters from a (grand canonical) thermal source or ii) subsequent formation of the clusters by coalescence of single nucleons. While both approaches have been successfully applied in the past it has not yet been clarified which of the two mechanisms dominates the cluster production. We propose to use recently developed event-by-event techniques to study particle multiplicity fluctuations on nuclear clusters and employ this analysis to the deuteron number fluctuations to disentangle the two production mechanisms. We argue that for a grand canonical cluster formation, the cluster fluctuations will follow Poisson distribution, while for the coalescence scenario, the fluctuations will strongly deviate from the Poisson expectation. We estimate the effect to be 10% for the variance and up to a factor of 5 for the kurtosis of the deuteron number multiplicity distribution. Our prediction can be tested in the beam energy scan program at RHIC as well as experiments at the FAIR and NICA facilities.

preprint2016arXiv

Influence of kinematic cuts on the net charge distribution

The higher moments of the net charge distributions, e.g. the skewness and kurtosis, are studied within an infinite hadronic matter calculation in a transport approach. By dividing the box into several parts, the volume dependence of the fluctuations is investigated. After confirming that the initial distributions follow the expectations from a binomial distribution, the influence of quantum number conservation in this case the net charge in the system on the higher moments is evaluated. For this purpose, the composition of the hadron gas is adjusted and only pions and $ρ$ mesons are simulated to investigate the charge conservation effect. In addition, the effect of imposing kinematic cuts in momentum space is analysed. The role of resonance excitations and decays on the higher moments can also be studied within this model. This work is highly relevant to understand the experimental measurements of higher moments obtained in the RHIC beam energy scan and their comparison to lattice results and other theoretical calculations assuming infinite matter.

preprint2015arXiv

Glueballs amass at RHIC and LHC Colliders! - The early quarkless 1st order phase transition at $T=270$ MeV - from pure Yang-Mills glue plasma to GlueBall-Hagedorn states

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.

preprint2015arXiv

Sub-threshold $ϕ$ and $Ξ^-$ production by high mass resonances with UrQMD

We present a possible explanation for the deep sub-threshold, $ϕ$ and $Ξ^-$ production yields measured with the HADES experiment in Ar+KCl reactions at $E_{\mathrm{lab}}=1.76$ A GeV and present predictions for Au+Au reactions at $E_{\mathrm{lab}}=1.23$ A GeV. To explain the surprisingly high yields of $ϕ$ and $Ξ^-$ hadrons we propose new decay channels for high mass baryon resonances. These new decay channels are constrained by elementary $\mathrm{p+p}\rightarrow \mathrm{p+p+}ϕ$ cross sections, and $Ξ^-$ production in p+Nb. Based on the fits to the elementary reactions one obtains a satisfactorily description of $ϕ$ and $Ξ^-$ production in deep sub-threshold Ar+KCl reactions as well as the observed nuclear transparency ratio in proton induced $ϕ$ production in cold nuclear matter. The results implicate that no new medium effects are required to describe the rare strange particle production data in low energy nuclear collisions.

preprint2015arXiv

The beam energy dependence of collective flow in heavy ion collisions

The major goals of heavy ion research are to explore the phase diagram of quantum chromodynamics (QCD) and to investigate the properties of the quark gluon plasma (QGP), a new state of matter created at high temperatures and/or densities. Collective anisotropic flow is one of the most promising observables to gain insights about the properties of the system created in relativistic heavy ion reactions. The current status of the beam energy dependence of the first three Fourier coefficients of the azimuthal distribution of the produced particles $v_1$ to $v_3$ within hybrid transport plus hydrodynamics approaches are summarized.

preprint2014arXiv

Do lattice data constrain the vector interaction strength of QCD?

We show how repulsive interactions of deconfined quarks as well as confined hadrons have an influence on the baryon number susceptibilities and the curvature of the chiral pseudo critical line in effective models of QCD. We discuss implications and constraints for the vector interaction strength from comparisons to lattice QCD and comment on earlier constraints, extracted from the curvature of the transition line of QCD and compact star observables. Our results clearly point to a strong vector repulsion in the hadronic phase and near-zero repulsion in the deconfined phase.

preprint2014arXiv

What the collective flow excitation function can tell about the quark-gluon plasma

Recent STAR data from the RHIC beam energy scan (BES) show that the midrapidity slope dv1/dy of the directed flow v1 of net-protons changes sign twice within the collision energy range 7.7 - 39 GeV. To investigate this phenomenon, we study the collision energy dependence of v1 utilizing a Boltzmann + hydrodynamics hybrid model. Calculations with dynamically evolved initial and final state show no qualitative difference between an equation of state with a cross-over and one with a first-order phase transition, in contrast to earlier pure fluid predictions. Furthermore, our analysis of the elliptic flow v2 shows that pre-equilibrium transport dynamics are partially compensating for the diminished elliptic flow production in the hydrodynamical phase at lower energies, which leads to a qualitative agreement with STAR BES results in midcentral collisions. No compensation from transport is found in our model for integrated v3, which decreases from approx 0.02 at sqrt(sNN) = 27 GeV to approx 0.005 at sqrt(sNN) = 7.7 GeV in midcentral collisions.

preprint2013arXiv

Conserved Charge Fluctuations in a Chiral Hadronic Model including Hadrons and Quarks

In this work the susceptibility coefficients of the strange and non-strange quark number of second and fourth order are presented. The results at zero baryonchemical potential are obtained using a well tested chiral effective model including all known hadron degrees of freedom and additionally incorporating quarks and gluons in a PNJL model approach. Quark number susceptibilities are sensitive to the fundamental degrees of freedom in the model and signal the shift from massive hadrons to light quarks at the deconfinement transition by a sharp rise at the critical temperature. Furthermore, all susceptibilities are found to be largely suppressed by repulsive vector field interactions of the particles. In the hadronic sector vector repulsion of baryon resonances restrains fluctuations to a region determined in lattice QCD. However, in the quark sector above $T_c$ even small vector field interactions of quarks quench all fluctuations almost completely. For this reason, vector field interactions for quarks have to vanish in the deconfinement limit.

preprint2013arXiv

Correlated D-meson decays competing against thermal QGP dilepton radiation

The QGP that might be created in ultrarelativistic heavy-ion collisions is expected to radiate thermal dilepton radiation. However, this thermal dilepton radiation interferes with dileptons originating from hadron decays. In the invariant mass region between the $ϕ$ and $J/Ψ$ peak ($1\,$GeV$\lesssim M_{\ell^+ \ell^-} \lesssim 3 \,$GeV) the most substantial background of hadron decays originates from correlated D$\bar{\mathrm{D}}$-meson decays. We evaluate this background using a Langevin simulation for charm quarks. As background medium we utilize the well-tested UrQMD-hybrid model. The required drag and diffusion coefficients are taken from a resonance approach. The decoupling of the charm quarks from the hot medium is performed at a temperature of $130\,$MeV and as hadronization mechanism a coalescence approach is chosen. This model for charm quark interactions with the medium has already been successfully applied to the study of the medium modification and the elliptic flow at FAIR, RHIC and LHC energies. In this proceeding we present our results for the dilepton radiation from correlated D$\bar{\mathrm{D}}$ decays at RHIC energy in comparison to PHENIX measurements in the invariant mass range between 1 and 3 GeV using different interaction scenarios. These results can be utilized to estimate the thermal QGP radiation.

preprint2013arXiv

Dileptons from correlated D- and $\bar{\text{D}}$-meson decays in the invariant mass range of the QGP thermal radiation using the UrQMD hybrid model

Relativistic heavy-ion collisions produce a hot and dense thermalized medium, that is expected to emit thermal radiation in form of dileptons. These dileptons are not affected by the strong force and are therefore a clean probe for the possible creation of a Quark Gluon Plasma (QGP). However, electroweak decays of open-charm mesons are expected to induce a substantial background in the invariant mass region between the $ϕ$ and $J/Ψ$ peak ($1\,\text{GeV} \lesssim M_{\ell^+ \ell^-} \lesssim 3 \,\text{GeV}$) of the thermal QGP radiation. To evaluate this background radiation we apply a Langevin approach for the transport of charm quarks in the UrQMD (hydrodynamics + Boltzmann) hybrid model. Due to the inclusion of event-by-event fluctuations and a full (3+1)-dimensional hydrodynamic evolution, the UrQMD hybrid approach provides a more realistic model for the evolution of the matter produced in heavy ion collisions as compared to simple homogeneous fireball expansions usually employed before. As drag and diffusion coefficients we use a resonance approach for elastic heavy-quark scattering and assume a decoupling temperature of the charm quarks from the hot medium of $130\,\text{MeV}$. For the hadronization of the charm quarks we employ a coalescence approach at the decoupling temperature from the medium. In this letter we present our calculations of the D-meson correlations and the invariant mass spectra of the dilepton decays in heavy-ion collisions at FAIR, RHIC, and LHC energies using different interaction scenarios.

preprint2013arXiv

Elliptic flow and nuclear modification factors of D-mesons at FAIR in a Hybrid-Langevin approach

The Compressed Baryonic Matter (CBM) experiment at the Facility for Anti-proton and Ion Research (FAIR) will provide new possibilities for charm-quark ($D$-meson) observables in heavy-ion collisions at low collision energies and high baryon densities. To predict the collective flow and nuclear modification factors of charm quarks in this environment, we apply a Langevin approach for the transport of charm quarks in the UrQMD (hydrodynamics + Boltzmann) hybrid model. Due to the inclusion of event-by-event fluctuations and a full (3+1) dimensional hydrodynamical evolution, the UrQMD hybrid approach provides a realistic evolution of the matter produced in heavy-ion collisions. As drag and diffusion coefficients we use a resonance approach for elastic heavy-quark scattering and assume a decoupling temperature of the charm quarks from the hot medium of $130\, \MeV$. Hadronization of the charm quarks to $D$-mesons by coalescence is included. Since the initial charm-quark distribution at FAIR is unknown, we utilize two different initial charm-quark distributions in our approach to estimate the uncertainty of these predictions. We present calculations of the nuclear modification factor, $R_{AA}$, as well as for the elliptic flow, $v_2$, in Pb+Pb collisions at $E_{lab}=25\,\text{AGeV}$. The different medium modifications of $D$-mesons and $\bar{D}$-mesons at high baryon-chemical potential are explored by modified drag- and diffusion-coefficients using the corresponding fugacity factor. Here we find a considerably larger medium modification for $\bar{\mathrm{D}}$- than for $D$-mesons.

preprint2013arXiv

Hadron Formation in Relativistic Nuclear Collisions and the QCD Phase Diagram

We analyze hadrochemical freeze-out in central Pb+Pb collisions at CERN SPS and LHC energies. We determine the effects of baryon and antibaryon annihilation and/or regeneration occurring during the final cascade expansion stage of the collisions by deriving "survival factors" for each hadronic species and all energies considered by employing the UrQMD hybrid model. These survival factors are shown to resemble the pattern of data deviation from the statistical equilibrium calculations with the statistical model. We apply them in the SM data analysis, obtaining new points in the $T-μ_B$ plane which follow closely the parton-hadron phase boundary predicted by lattice QCD.

preprint2013arXiv

Non-thermal p/pi ratio at LHC as a consequence of hadronic final state interactions

Recent LHC data on Pb+Pb reactions at sqrt(s_{NN})=2.7 TeV suggests that the p/pi is incompatible with thermal models. We explore several hadron ratios (K/pi, p/pi, Lambda/pi, Xi/pi) within a hydrodynamic model with hadronic after burner, namely UrQMD 3.3, and show that the deviations can be understood as a final state effect. We propose the p/pi as an observable sensitive on whether final state interactions take place or not. The measured values of the hadron ratios do then allow to gauge the transition energy density from hydrodynamics to the Boltzmann description. We find that the data can be explained with transition energy densities of 840 +- 150 MeV/fm^3.

preprint2013arXiv

Phase Structure of QCD Matter in a Chiral Effective Model with Quarks

Using a unified hadron-quark effective model for the QCD equation of state, this paper studies the phase structure of strongly interacting matter in a wide range of temperature and baryonchemical potential. At small potentials the model yields a smooth cross-over to chirally restored matter with a transition temperature and curvature in line with recent lattice QCD estimates and thermal model fits of freeze-out curves. Trajectories of constant entropy per net baryon number show a clear dependence on the particle composition in the model and on repulsive vector field interactions. Although the model might feature a critical end-point at a rather high baryonchemical potential and low temperature, probing it in heavy-ion collisions might be highly challenging due to a large thermodynamic spread of matter in the collision fireball.

preprint2013arXiv

QCD Equation of State From a Chiral Hadronic Model Including Quark Degrees of Freedom

This work presents an effective model for strongly interacting matter and the QCD equation of state (EoS). The model includes both hadron and quark degrees of freedom and takes into account the transition of chiral symmetry restoration as well as the deconfinement phase transition. At low temperatures $T$ and baryonic densities $ρ_B$ a hadron resonance gas is described using a SU(3)-flavor sigma-omega model and a quark phase is introduced in analogy to PNJL models for higher $T$ and $ρ_B$. In this way, the correct asymptotic degrees of freedom are used in a wide range of $T$ and $ρ_B$. Here, results of this model concerning the chiral and deconfinement phase transitions and thermodynamic model properties are presented. Large hadron resonance multiplicities in the transition region emphasize the importance of heavy-mass resonance states in this region and their impact on the chiral transition behavior. The resulting phase diagram of QCD matter at small chemical potentials is in line with latest lattice QCD and thermal model results.

preprint2013arXiv

Spinodal density enhancements in simulations of relativistic nuclear collisions

We recently introduced a fluid-dynamical model for simulating relativistic nuclear collisions in the presence of a first-order phase transition and made explorative studies of head-on lead-lead collisions. We give here a more detailed account of this novel theoretical tool and carry out more exhaustive studies of the phase-separation dynamics. Extracting the density enhancement caused by the spinodal instabilities, the associated clump size distribution, and the resulting transverse flow velocity, we examine the sensitivity of these quantities to the strength of the gradient term that promotes the phase separation, to the details of the initial density fluctuations that form the seeds for the subsequent amplification, and to the equation of state.

preprint2012arXiv

Baryon Resonances in a Chiral Hadronic Model for the QCD Equation of State

In this paper we study the influence of hadronic resonances on the phase diagram calculated with an effective chiral flavour SU(3) model. We show that varying the couplings of the baryonic resonances to the attractive scalar and the repulsive vector fields has a major impact on the order and location of the chiral phase transition and the possible existence of a critical end point as well as on the thermodynamic properties of the model. Furthermore, we study (strange) quark number fluctuations and show the related susceptibilities both at zero baryochemical potential and when crossing the phase transition line at three different points in the T-mu plane. We obtain the best agreement with current lattice data if we choose a rather strong vector coupling which in our model limits the phase transition to a smooth crossover and implies the non-existence of a critical end point.

preprint2012arXiv

Charm quark transport in Pb+Pb reactions at sqrt{s}_NN=2.76 TeV from a (3+1) dimensional hybrid approach

We implement a Langevin approach for the transport of charm quarks in the UrQMD (hydrodynamics + Boltzmann) hybrid model. Due to the inclusion of event-by-event fluctuations and a full (3+1) dimensional hydrodynamic evolution, this approach provides a more realistic model for the evolution of the matter produced in heavy ion collisions as compared to simple homogeneous fireball expansions usually employed. As drag and diffusion coefficients we use a resonance approach for elastic heavy-quark scattering and assume a decoupling temperature of the charm quarks from the hot medium of 130 MeV. A coalescence approach at the decoupling temperature for the hadronization of the charm quarks to D-mesons is also included. We present calculations of the nuclear modification factor R_AA as well as the elliptic flow v_2 in Pb+Pb collisions at sqrt s_NN=2.76 TeV. The comparison to ALICE measurements shows a very good agreement with our calculations.

preprint2012arXiv

Extraction of the sound velocity from rapidity spectra: Evidence for QGP formation at FAIR/RHIC-BES energies

We analyze longitudinal pion spectra from \sqrt{s_{NN}}= 2 GeV to \sqrt {s_{\rm NN}}=20 GeV within Landau's hydrodynamical model and the UrQMD hybrid approach. From the measured data on the widths of the pion rapidity spectra, we extract the sound velocity c_s^2 in the dense stage of the reactions. It is found that the sound velocity has a local minimum (indicating a softest point in the equation of state, EoS) at \sqrt{s_{NN}}= 4-9 GeV, an energy range accessible at the Facility for Antiproton and Ion Research (FAIR) as well as the RHIC-Beam Energy Scan (RHIC-BES). This softening of the EoS is compatible with the formation of a QGP at the onset of deconfinement. The extracted sound velocities are then used to calculate an excitation function for the mean transverse mass of pions from the hybrid model. We find that, above \sqrt{s_{NN}} \approx 10 GeV, even the lowest c_s^{2} gives a considerably larger <m_T> of pions compared to data.

preprint2012arXiv

Heavy quark transport in heavy ion collisions at RHIC and LHC within the UrQMD transport model

We have implemented a Langevin approach for the transport of heavy quarks in the UrQMD hybrid model. The UrQMD hybrid approach provides a realistic description of the background medium for the evolution of relativistic heavy ion collisions. We have used two different sets of drag and diffusion coefficients, one based on a $T$-Matrix approach and one based on a resonance model for the elastic scattering of heavy quarks within the medium. In case of the resonance model we have investigated the effects of different decoupling temperatures of the heavy quarks from the medium, ranging between $130\,\text{MeV}$ and $180\,\text{MeV}$. We present calculations of the nuclear modification factor $R_{AA}$, as well as of the elliptic flow $v_2$ in Au+Au collisions at $\sqrt{s_{NN}}=200\,\text{GeV}$ and Pb+Pb collisions at $\sqrt{s_{NN}}=2.76\,\text{TeV}$. To make our results comparable to experimental data at RHIC and LHC we have implemented a Peterson fragmentation and a quark coalescence approach followed by the semileptonic decay of the D- and B-mesons to electrons. We find that our results strongly depend on the decoupling temperature and the hadronization mechanism. At a decoupling temperature of $130\,\text{MeV}$ we reach a good agreement with the measurements at both, RHIC and LHC energies, simultaneously for the elliptic flow $v_2$ and the nuclear modification factor $R_{AA}$.

preprint2012arXiv

Spinodal amplification of density fluctuations in fluid-dynamical simulations of relativistic nuclear collisions

Extending a previously developed two-phase equation of state, we simulate head-on relativistic lead-lead collisions with fluid dynamics, augmented with a finite-range term, and study the effects of the phase structure on the evolution of the baryon density. For collision energies that bring the bulk of the system into the mechanically unstable spinodal region of the phase diagram, the density irregularities are being amplified significantly. The resulting density clumping may be exploited as a signal of the phase transition, possibly through an enhanced production of composite particles.

preprint2010arXiv

Conical Emission from Shock Waves in Ne(1-20 AGeV)+U Collisions

The formation and propagation of high-density compression waves, e.g. Mach shock waves, in cold nuclear matter is studied by simulating high-energy nucleus-nucleus collisions of Ne with U in the energy range from E_lab = 0.5 AGeV to 20 AGeV. In an ideal hydrodynamic approach, the high-density shock wave created by the small Ne nucleus passing through the heavy U nucleus is followed by a slower and more dilute Mach shock wave which causes conical emission of particles at the Mach cone angle. The conical emission originates from low-density regions with a small flow velocity comparable to the speed of sound. Moreover, it is shown that the angular distributions of emitted baryons clearly distinguish between a hydrodynamic approach and binary cascade processes used in the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) transport model.

preprint2010arXiv

Strangeness production in hadronic models and recombination models

We present recent results on the production, spectra and elliptic flow of strange particles in dynamic simulations employing hadronic degrees of freedom and from recombination models. The main focus will be on the Ultra-relativistic Molecular Dynamics (UrQMD) Boltzmann approach to relativistic heavy ion collisions and a hybrid approach with intermediate hydrodynamic evolution based on UrQMD (available for download as UrQMD v3.3). Compared to the standard binary collision approach, an enhancement of the strange particle particle yields is found in the hybrid approach due to the assumption of local equilibration. The production origins of the Phi-meson in the hybrid approach are studied in further detail. We also present results on the transverse momentum spectra of baryon to meson ratios of strange particles. Due to the approximate energy independent scaling of this ratio as a function of p_T we argue, that a maximum in these spectra may not be a unique sign for quark coalescence but can be understood in terms of flow and fragmentation.