Refine
Year of publication
Document Type
- Article (18)
- Conference Proceeding (4)
- Doctoral Thesis (1)
- Preprint (1)
Language
- English (24)
Has Fulltext
- yes (24)
Is part of the Bibliography
- no (24)
Keywords
- transport theory (3)
- Relativistic heavy-ion collisions (2)
- correlations (2)
- hot spots (2)
- initial state (2)
- relativistic hydrodynamics (2)
- Boltzmann equation (1)
- Boltzmann-Gleichung (1)
- CLVisc (1)
- Dynamic transport (1)
Institute
A mechanism for locally density-dependent dynamic parton rearrangement and fusion has been implemented into the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) approach. The same mechanism has been previously built in the Quark Gluon String Model (QGSM). This rearrangement and fusion approach based on parton coalescence ideas enables the description of multi-particle interactions, namely 3 -> 3 and 3 -> 2, between (pre)hadronic states in addition to standard binary interactions. The UrQMD model (v2.3) extended by these additional processes allows to investigate implications of multi-particle interactions on the reaction dynamics of ultrarelativistic heavy ion collisions. The mechanism, its implementation and first results of this investigation are presented and discussed.
The effects of the onset of deconfinement on longitudinal and transverse flow are studied. First, we analyze longitudinal pion spectra from Elab = 2A GeV to √sNN = 200 GeV within Landau’s hydrodynamical model and the UrQMD transport approach. From the measured data on the widths of the pion rapidity spectra, we extract the sound velocity c2s in the early 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 Ebeam = 30A GeV. This softening of the EoS is compatible with the assumption of the formation of a mixed phase at the onset of deconfinement. Furthermore, the energy excitation function of elliptic flow (v2) from Ebeam = 90A MeV to √sNN = 200 GeV is explored within the UrQMD framework and discussed in the context of the available data. The transverse flow should also be sensitive to changes in the equation of state. Therefore, the underestimation of elliptic flow by the UrQMD model calculation above Elab = 30A GeV might also be explained by assuming a phase transition from a hadron gas to the quark gluon plasma around this energy. This would be consistent with the model calculations, indicating a transition from hadronic matter to “string matter” in this energy range.
In this thesis the first fully integrated Boltzmann+hydrodynamics approach to relativistic heavy ion reactions has been developed. After a short introduction that motivates the study of heavy ion reactions as the tool to get insights about the QCD phase diagram, the most important theoretical approaches to describe the system are reviewed. To model the dynamical evolution of the collective system assuming local thermal equilibrium ideal hydrodynamics seems to be a good tool. Nowadays, the development of either viscous hydrodynamic codes or hybrid approaches is favoured. For the microscopic description of the hadronic as well as the partonic stage of the evolution transport approaches have beeen successfully applied, since they generate the full phse-space dynamics of all the particles. The hadron-string transport approach that this work is based on is the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) approach. It constitutes an effective solution of the relativistic Boltzmann equation and is restricted to binary collisions of the propagated hadrons. Therefore, the Boltzmann equation and the basic assumptions of this model are introduced. Furthermore, predictions for the charged particle multiplicities at LHC energies are made. The next step is the development of a new framework to calculate the baryon number density in a transport approach. Time evolutions of the net baryon number and the quark density have been calculated at AGS, SPS and RHIC energies and the new approach leads to reasonable results over the whole energy range. Studies of phase diagram trajectories using hydrodynamics are performed as a first move into the direction of the development of the hybrid approach. The hybrid approach that has been developed as the main part of this thesis is based on the UrQMD transport approach with an intermediate hydrodynamical evolution for the hot and dense stage of the collision. The initial energy and baryon number density distributions are not smooth and not symmetric in any direction and the initial velocity profiles are non-trivial since they are generated by the non-equilibrium transport approach. The fulll (3+1) dimensional ideal relativistic one fluid dynamics evolution is solved using the SHASTA algorithm. For the present work, three different equations of state have been used, namely a hadron gas equation of state without a QGP phase transition, a chiral EoS and a bag model EoS including a strong first order phase transition. For the freeze-out transition from hydrodynamics to the cascade calculation two different set-ups are employed. Either an in the computational frame isochronous freeze-out or an gradual freeze-out that mimics an iso-eigentime criterion. The particle vectors are generated by Monte Carlo methods according to the Cooper-Frye formula and UrQMD takes care of the final decoupling procedure of the particles. The parameter dependences of the model are investigated and the time evolution of different quantities is explored. The final pion and proton multiplicities are lower in the hybrid model calculation due to the isentropic hydrodynamic expansion while the yields for strange particles are enhanced due to the local equilibrium in the hydrodynamic evolution. The elliptic flow values at SPS energies are shown to be in line with an ideal hydrodynamic evolution if a proper initial state is used and the final freeze-out proceeds gradually. The hybrid model calculation is able to reproduce the experimentally measured integrated as well as transverse momentum dependent $v_2$ values for charged particles. The multiplicity and mean transverse mass excitation function is calculated for pions, protons and kaons in the energy range from $E_{\rm lab}=2-160A~$GeV. It is observed that the different freeze-out procedures have almost as much influence on the mean transverse mass excitation function as the equation of state. The experimentally observed step-like behaviour of the mean transverse mass excitation function is only reproduced, if a first order phase transition with a large latent heat is applied or the EoS is effectively softened due to non-equilibrium effects in the hadronic transport calculation. The HBT correlation of the negatively charged pion source created in central Pb+Pb collisions at SPS energies are investigated with the hybrid model. It has been found that the latent heat influences the emission of particles visibly and hence the HBT radii of the pion source. The final hadronic interactions after the hydrodynamic freeze-out are very important for the HBT correlation since a large amount of collisions and decays still takes place during this period.
There is little doubt that Quantumchromodynamics (QCD) is the theory which describes strong interaction physics. Lattice gauge simulations of QCD predict that in the m,T plane there is a line where a transition from confined hadronic matter to deconfined quarks takes place. The transition is either a cross over (at low m) or of first order (at high m). It is the goal of the present and future heavy ion experiment at RHIC and FAIR to study this phase transition at different locations in the m,T plane and to explore the properties of the deconfined phase. It is the purpose of this contribution to discuss some of the observables which are considered as useful for this purpose.
The goal of heavy ion reactions at low beam energies is to explore the QCD phase diagram at high net baryon chemical potential. To relate experimental observations with a first order phase transition or a critical endpoint, dynamical approaches for the theoretical description have to be developed. In this summary of the corresponding plenary talk, the status of the dynamical modeling including the most recent advances is presented. The remaining challenges are highlighted and promising experimental measurements are pointed out.
While the existence of a strongly interacting state of matter, known as “quark-gluon plasma” (QGP), has been established in heavy ion collision experiments in the past decade, the task remains to map out the transition from the hadronic matter to the QGP. This is done by measuring the dependence of key observables (such as particle suppression and elliptic flow) on the collision energy of the heavy ions. This procedure, known as "beam energy scan", has been most recently performed at the Relativistic Heavy Ion Collider (RHIC).
Utilizing a Boltzmann+hydrodynamics hybrid model, we study the collision energy dependence of initial state eccentricities and the final state elliptic and triangular flow. This approach is well suited to investigate the relative importance of hydrodynamics and hadron transport at different collision energies.
Simulating Many Accelerated Strongly-interacting Hadrons (SMASH) is a new hadronic transport approach designed to describe the non-equilibrium evolution of heavy-ion collisions. The production of strange particles in such systems is enhanced compared to elementary reactions (Blume and Markert 2011), providing an interesting signal to study. Two different strangeness production mechanisms are discussed: one based on resonances and another using forced canonical thermalization. Comparisons to experimental data from elementary collisions are shown.
A primordial state of matter consisting of free quarks and gluons that existed in the early universe a few microseconds after the Big Bang is also expected to form in high-energy heavy-ion collisions. Determining the equation of state (EoS) of such a primordial matter is the ultimate goal of high-energy heavy-ion experiments. Here we use supervised learning with a deep convolutional neural network to identify the EoS employed in the relativistic hydrodynamic simulations of heavy ion collisions. High-level correlations of particle spectra in transverse momentum and azimuthal angle learned by the network act as an effective EoS-meter in deciphering the nature of the phase transition in quantum chromodynamics. Such EoS-meter is model-independent and insensitive to other simulation inputs including the initial conditions for hydrodynamic simulations.
We present a systematic study of the normalized symmetric cumulants, NSC(n,m), at the eccentricity level in proton-proton interactions at within a wounded hot spot approach. We focus our attention on the influence of spatial correlations between the proton constituents, in our case gluonic hot spots, on this observable. We notice that the presence of short-range repulsive correlations between the hot spots systematically decreases the values of and in mid- to ultra-central collisions while increases them in peripheral interactions. In the case of we find that, as suggested by data, an anti-correlation of and in ultra-central collisions, i.e. , is possible within the correlated scenario while it never occurs without correlations when the number of gluonic hot spots is set to three. We attribute this fact to the decisive role of correlations on enlarging the probability of interaction topologies that reduce the value of and, eventually, make it negative. Further, we explore the dependence of our conclusions on the number of hot spots, the values of the hot spot radius and the repulsive core distance. Our results add evidence to the idea that considering spatial correlations between the subnucleonic degrees of freedom of the proton may have a strong impact on the initial state properties of proton-proton interactions [1].
In this proceeding we review our recent work using supervised learning with a deep convolutional neural network (CNN) to identify the QCD equation of state (EoS) employed in hydrodynamic modeling of heavy-ion collisions given only final-state particle spectra ρ(pT, Ф). We showed that there is a traceable encoder of the dynamical information from phase structure (EoS) that survives the evolution and exists in the final snapshot, which enables the trained CNN to act as an effective “EoS-meter” in detecting the nature of the QCD transition.
We investigate the long-standing question of the effect of proton-antiproton annihilation on the (anti-)proton yield, while respecting detailed balance for the five-body back-reaction for the first time in a full microscopic description of the late stages of heavy-ion collisions. This is achieved by employing a stochastic collision criterion in a hadronic transport approach (SMASH), which is used to account for the regeneration of (anti-)protons via 5π→p¯p. We investigate Au+Au and Pb+Pb collisions from √sNN=17.3GeV−5.02 TeV in a viscous hybrid approach. Our results show that back-reactions happen for a fraction of 15%–20% of all annihilations, independent of the beam energy or centrality of the system. The inclusion of the back-reaction results in the regeneration of half of the (anti-)proton yield lost to annihilations at midrapidity. We also find that, concerning the multiplicities, treating the back-reaction as a chain of two-body reactions is equivalent to a single 5-to-2 reaction.
Motivated by a recent finding of an exact solution of the relativistic Boltzmann equation in a Friedmann–Robertson–Walker spacetime, we implement this metric into the newly developed transport approach Simulating Many Accelerated Strongly-interacting Hadrons (SMASH). We study the numerical solution of the transport equation and compare it to this exact solution for massless particles. We also compare a different initial condition, for which the transport equation can be independently solved numerically. Very nice agreement is observed in both cases. Having passed these checks for the SMASH code, we study a gas of massive particles within the same spacetime, where the particle decoupling is forced by the Hubble expansion. In this simple scenario we present an analysis of the freeze-out times, as function of the masses and cross sections of the particles. The results might be of interest for their potential application to relativistic heavy-ion collisions, for the characterization of the freeze-out process in terms of hadron properties.
We describe two independent frameworks which provide unambiguous determinations of the deconfinement and the decoupling conditions of a relativistic gas at finite temperature. First, we use the Polyakov-Nambu-Jona–Lasinio model to compute meson and baryon masses at finite temperature and determine their melting temperature as a function of their strangeness content. Second, we analyze a simple expanding gas within a Friedmann-Robertson-Walker metric, which admits a well-defined decoupling mechanism. We examine the decoupling time as a function of the particle mass and cross section. We find evidences of an inherent dependence of the hadronization and freeze-out conditions on flavor, and on mass and cross section, respectively.
We present results on Hanbury Brown-Twiss (HBT) radii extracted from the Ultra-relativistic Molecular Dynamics (UrQMD) approach to relativistic heavy ion collisions. The present investigation provides a comparison of results from pure hadronic transport calculations to a Boltzmann + Hydrodynamic hybrid approach with an intermediate hydrodynamic phase. For the hydrodynamic phase different Equations of State (EoS) have been employed, i.e. bag model, hadron resonance gas and a chiral EoS. The influence of various freeze-out scenarios has been investigated and shown to be negligible if hadronic rescatterings after the hydrodynamic evolution are included. Furthermore, first results of the source tilt from azimuthal sensitive HBT and the direct extraction from the transport model are presented and exhibit a very good agreement with E895 data at AGS.
Previous calculations of the shear viscosity to entropy density ratio in the hadron gas have failed to reach a consensus, with η/s predictions differing by almost an order of magnitude. This work addresses and solves this discrepancy by providing an independent extraction of η/s using the newly-developed SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) transport code and the Green-Kubo formalism. We compare the results from SMASH with numerical solutions of the Boltzmann equation for various systems using the Chapman-Enskog expansion as well as previous results in the literature. Substantial deviations of the coefficient are found between transport approaches mainly based on resonance propagation with finite lifetime (such as SMASH) and other (semi-analytical) approaches with energy-dependent cross-sections, where interactions do not introduce a timescale other than the inverse scattering rate. Our conclusion is that long- lived resonances strongly affect the transport properties of the system, resulting in significant differences in η/s with respect to other approaches where binary collisions dominate. We argue that the relaxation time of the system —which characterizes the shear viscosity— is determined by the interplay between the mean- free time and the lifetime of resonances. We finally show how an artificial shortening of the resonance lifetimes or the addition of a background elastic cross section nicely interpolate between the two discrepant results.
The state-of-the-art pattern recognition method in machine learning (deep convolution neural network) is used to identify the equation of state (EoS) employed in the relativistic hydrodynamic simulations of heavy ion collisions. High-level correlations of particle spectra in transverse momentum and azimuthal angle learned by the network act as an effective EoS-meter in deciphering the nature of the phase transition in QCD. The EoS-meter is model independent and insensitive to other simulation inputs including the initial conditions and shear viscosity for hydrodynamic simulations. Through this study we demonstrate that there is a traceable encoder of the dynamical information from the phase structure that survives the evolution and exists in the final snapshot of heavy ion collisions and one can exclusively and effectively decode these information from the highly complex final output with machine learning when traditional methods fail. Besides the deep neural network, the performance of traditional machine learning classifiers are also provided.
Effects of a phase transition on HBT correlations in an integrated Boltzmann+hydrodynamics approach
(2009)
A systematic study of HBT radii of pions, produced in heavy ion collisions in the intermediate energy regime (SPS), from an integrated (3+1)d Boltzmann+hydrodynamics approach is presented. The calculations in this hybrid approach, incorporating an hydrodynamic stage into the Ultra-relativistic Quantum Molecular Dynamics transport model, allow for a comparison of different equations of state retaining the same initial conditions and final freeze-out. The results are also compared to the pure cascade transport model calculations in the context of the available data. Furthermore, the effect of different treatments of the hydrodynamic freeze-out procedure on the HBT radii are investigated. It is found that the HBT radii are essentially insensitive to the details of the freeze-out prescription as long as the final hadronic interactions in the cascade are taken into account. The HBT radii RL and RO and the RO/RS ratio are sensitive to the EoS that is employed during the hydrodynamic evolution. We conclude that the increased lifetime in case of a phase transition to a QGP (via a Bag Model equation of state) is not supported by the available data.
We apply a coupled transport-hydrodynamics model to discuss the production of multi-strange meta-stable objects in Pb + Pb reactions at the FAIR facility. In addition to making predictions for yields of these particles we are able to calculate particle dependent rapidity and momentum distributions. We argue that the FAIR energy regime is the optimal place to search for multi-strange baryonic object (due to the high baryon density, favoring a distillation of strangeness). Additionally, we show results for strangeness and baryon density fluctuations. Using the UrQMD model we calculate the strangeness separation in phase space which might lead to an enhanced production of MEMOs compared to models that assume global thermalization.
ϕ-meson production in In–In collisions at Elab=158A GeV: Evidence for relics of a thermal phase
(2010)
Yields and transverse mass distributions of the ϕ-mesons reconstructed in the ϕ→μ+μ− channel in In+In collisions at Elab=158A GeV are calculated within an integrated Boltzmann+hydrodynamics hybrid approach based on the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) transport model with an intermediate hydrodynamic stage. The analysis is performed for various centralities and a comparison with the corresponding NA60 data in the muon channel is presented. We find that the hybrid model, that embeds an intermediate locally equilibrated phase subsequently mapped into the transport dynamics according to thermal phase-space distributions, gives a good description of the experimental data, both in yield and slope. On the contrary, the pure transport model calculations tend to fail in catching the general properties of the ϕ meson production: not only the yield, but also the slope of the mT spectra, compare poorly with the experimental observations at top SPS energies.
In this work, largely based on [J. L. Albacete, A. Soto-Ontoso, Hot spots and the hollowness of proton-proton interactions at high energies, arXiv:1605.09176; J. L. Albacete, H. Petersen, A. Soto-Ontoso, Correlated wounded hot spots in proton-proton interactions, arXiv:1612.06274], we present a novel initial state geometry for proton-proton interactions. We rely on gluonic hot spots as effective degrees of freedom whose transverse positions inside the proton are correlated. We explore the impact of these non-trivial spatial correlations on the eccentricity and triangularity of the system following a Monte Carlo Glauber approach.
We compute the fermion spin distribution in the vortical fluid created in off-central high energy heavy-ion collisions. We employ the event-by-event (3+1)D viscous hydrodynamic model. The spin polarization density is proportional to the local fluid vorticity in quantum kinetic theory. As a result of strong collectivity, the spatial distribution of the local vorticity on the freeze-out hyper-surface strongly correlates to the rapidity and azimuthal angle distribution of fermion spins. We investigate the sensitivity of the local polarization to the initial fluid velocity in the hydrodynamic model and compute the global polarization of Λ hyperons by the AMPT model. The energy dependence of the global polarization agrees with the STAR data.
We present a systematic study on the influence of spatial correlations between the proton constituents, in our case gluonic hot spots, their size and their number on the symmetric cumulant SC(2,3), at the eccentricity level, within a Monte Carlo Glauber framework [J.L. Albacete, H. Petersen, A. Soto-Ontoso, Symmetric cumulants as a probe of the proton substructure at LHC energies, Phys. Lett. B778 (2018) 128–136. arXiv:1707.05592, doi:10.1016/j.physletb.2018.01.011]. When modeling the proton as composed by 3 gluonic hot spots, the most common assumption in the literature, we find that the inclusion of spatial correlations is indispensable to reproduce the negative sign of SC(2,3) in the highest centrality bins as dictated by data. Further, the subtle interplay between the different scales of the problem is discussed. To conclude, the possibility of feeding a 2+1D viscous hydrodynamic simulation with our entropy profiles is exposed.
Microscopic transport approaches are the tool to describe the non-equilibrium evolution in low energy collisions as well as in the late dilute stages of high-energy collisions. Here, a newly developed hadronic transport approach, SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) is introduced. The overall bulk dynamics in low energy heavy ion collisions is shown including the excitation function of elliptic flow employing several equations of state. The implications of this new approach for dilepton production are discussed and preliminary results for afterburner calculations at the highest RHIC energy are presented and compared to previous UrQMD results. A detailed understanding of a hadron gas with vacuum properties is required to establish the baseline for the exploration of the transition to the quark-gluon plasma in heavy ion collisions at high net baryon densities.
The intense photon fluxes from relativistic nuclei provide an opportunity to study photonuclear interactions in ultraperipheral collisions. The measurement of coherently photoproduced π+π−π+π− final states in ultraperipheral Pb-Pb collisions at sNN−−−√=5.02 TeV is presented for the first time. The cross section, dσ/dy, times the branching ratio (ρ→π+π+π−π−) is found to be 47.8±2.3 (stat.)±7.7 (syst.) mb in the rapidity interval |y|<0.5. The invariant mass distribution is not well described with a single Breit-Wigner resonance. The production of two interfering resonances, ρ(1450) and ρ(1700), provides a good description of the data. The values of the masses (m) and widths (Γ) of the resonances extracted from the fit are m1=1385±14 (stat.)±3 (syst.) MeV/c2, Γ1=431±36 (stat.)±82 (syst.) MeV/c2, m2=1663±13 (stat.)±22 (syst.) MeV/c2 and Γ2=357±31 (stat.)±49 (syst.) MeV/c2, respectively. The measured cross sections times the branching ratios are compared to recent theoretical predictions.