Researcher profile

Jan Steinheimer

Jan Steinheimer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 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&#39;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 &#34;after-burner&#34;. 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.