We study the properties of the strongly-coupled quark-gluon plasma with a multistage model of heavy ion collisions that combines the TRENTo initial condition ansatz, free-streaming, viscous relativistic hydrodynamics, and a relativistic hadronic transport. A model-to-data comparison with Bayesian inference is performed, revisiting assumptions made in previous studies. The role of parameter priors is studied in light of their importance towards the interpretation of results. We emphasize the use of closure tests to perform extensive validation of the analysis workflow before comparison with observations. Our study combines measurements from the Large Hadron Collider and the Relativistic Heavy Ion Collider, achieving a good simultaneous description of a wide range of hadronic observables from both colliders. The selected experimental data provide reasonable constraints on the shear and the bulk viscosities of the quark-gluon plasma at T ∼ 150-250 MeV, but their constraining power degrades at higher temperatures T 250 MeV. Furthermore, these viscosity constraints are found to depend significantly on how viscous corrections are handled in the transition from hydrodynamics to the hadronic transport. Several other model parameters, including the free-streaming time, show similar model sensitivity, while the initial condition parameters associated with the TRENTo ansatz are quite robust against variations of the particlization prescription. We also report on the sensitivity of individual observables to the various model parameters. Finally, Bayesian model selection is used to quantitatively compare the agreement with measurements for different sets of model assumptions, including different particlization models and different choices for which parameters are allowed to vary between RHIC and LHC energies. CONTENTS Pratt-Torrieri-Bernhard 10 D. Hadronic transport 11 IV. Specifying prior knowledge 11 V. Bayesian Parameter Estimation with a Statistical Emulator 13 A. Overview of Bayesian Parameter Estimation 13 B. Physical model emulator 14 C. Treatment of uncertainties 16 D. Sampling of the posterior 17 E. Maximizing the posterior 17 VI. Closure Tests 17 A. Validating Bayesian inference with closure tests 18 B. Guiding analyses with closure tests 18 37 A. Full posterior of model parameters 37 B. Posterior for LHC and RHIC independently 37 C. Validation of principal component analysis 37 D. Experimental covariance matrix 38 E. Reducing experimental uncertainty 39 F. Bulk relaxation time 39 G. Comparison to previous studies 40 1. Physics models 41 2. Prior distributions 42 3. Experimental data 42 H. Multistage model validation 42 1. Validation of second-order viscous hydrodynamics implementation 42 a. Validation against cylindrically symmetric external solution 43 2. SMASH 43 3. Comparison of JETSCAPE with hic-eventgen 45 4. The σ meson 46 5. Sampling particles on mass-shell 47 6. QCD equations of state with different hadron resonance gases 47 References 48

The deuteron yield in Pb+Pb collisions at √ s NN = 2.76 TeV is consistent with thermal production at a freeze out temperature of T = 155 MeV. The existence of deuterons with binding energy of 2.2 MeV at this temperature was described as "snowballs in hell" [P. Braun-Münzinger, B. Dönigus, and N. Löher, CERN Courier, August 2015]. We provide a microscopic explanation of this phenomenon, utilizing relativistic hydrodynamics and switching to a hadronic afterburner at the above-mentioned temperature of T = 155 MeV. The measured deuteron p T spectra and coalescence parameter B 2 (p T) are reproduced without free parameters, only by implementing experimentally known cross sections of deuteron reactions with hadrons, most importantly π d ↔ π np.

No abstract

Background: Simulations by transport codes are indispensable for extracting valuable physical information from heavy-ion collisions. Pion observables such as the π − /π + yield ratio are expected to be sensitive to the symmetry energy at high densities.Purpose: To evaluate, understand and reduce the uncertainties in transport-code results originating from different approximations in handling the production of ∆ resonances and pions. Methods:We compare ten transport codes under controlled conditions for a system confined in a box, with periodic boundary conditions, and initialized with nucleons at saturation density and at 60 MeV temperature. The reactions N N ↔ N∆ and ∆ ↔ N π are implemented, but the Pauli blocking and the mean-field potential are deactivated in the present comparison. Thus these are cascade calculations including pions and ∆ resonances. Results are compared to those from the two reference cases of a chemically equilibrated ideal gas mixture and of the rate equation. Results:For the numbers of ∆ and π, deviations from the reference values are observed in many codes, and they depend significantly on the size of the time step. These deviations are tied to different ways in ordering the sequence of reactions, such as collisions and decays, that take place in the same time step. Better agreements with the reference values are seen in the reaction rates and the number ratios among the isospin species of ∆ and π. Both the reaction rates and the number ratios are, however, affected by the correlations between particle positions, which are absent in the Boltzmann equation, but are induced by the way particle scatterings are treated in many of the transport calculations. The uncertainty in the transport-code predictions of the π − /π + ratio, after letting the existing ∆ resonances decay, is found to be within a few percent for the system initialized at n/p = 1.5. Conclusions:The uncertainty in the final π − /π + ratio in this simplified case of particles in a box is sufficiently small so that it does not strongly impact constraining the high-density symmetry energy from heavy-ion collisions. Most of the sources of

No abstract

The stopping of baryons in heavy ion collisions at beam momenta of p lab = 20 − 160A GeV is lacking a quantitative description within theoretical calculations. Heavy ion reactions at these energies are experimentally explored at the Super Proton Synchrotron (SPS) and the Relativistic Heavy Ion Collider (RHIC) and will be studied at future facilities such as FAIR and NICA. Since the net baryon density is determined by the amount of stopping, this is the pre-requisiste for any investigation of other observables related to structures in the QCD phase diagram such as a first-order phase transition or a critical endpoint. In this work we employ a string model for treating hadron-hadron interactions within a hadronic transport approach (SMASH, Simulating Many Accelerated Strongly-interacting Hadrons). Free parameters of the string excitation and decay are tuned to match experimental measurements in elementary proton-proton collisions. Afterwards, the model is applied to heavy ion collisions, where the experimentally observed change of the shape of the proton rapidity spectrum from a single peak structure to a double peak structure with increasing beam energy is reproduced. Heavy ion collisions provide the opportunity to study the formation process of string fragments in terms of formation times and reduced interaction cross-sections for pre-formed hadrons. A good agreement with the measured rapidity spectra of protons and pions is achieved while insights on the fragmentation process are obtained. In the future, the presented approach can be used to create event-by-event initial conditions for hybrid calculations.

Motivated by the theory of relativistic hydrodynamic fluctuations we make use of the Green-Kubo formula to compute the electrical conductivity and the (second-order) relaxation time of the electric current of an interacting hadron gas. We use the recently developed transport code SMASH to numerically solve the coupled set of Boltzmann equations implementing realistic hadronic interactions. In particular, we explore the role of the resonance lifetimes in the determination of the electrical relaxation time. As opposed to a previous calculation of the shear viscosity we observe that the presence of resonances with lifetimes of the order of the mean-free time does not appreciably affect the relaxation of the electric current fluctuations. We compare our results to other approaches describing similar systems, and provide the value of the electrical conductivity and the relaxation time for a hadron gas at temperatures between T = 60 MeV and T = 150 MeV.

No abstract