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

No abstract

In relativistic heavy-ion collisions, the production of heavy quarks at large transverse momenta is strongly suppressed compared to proton-proton collisions. In addition an unexpectedly large azimuthal anisotropy was observed for the emission of charmed hadrons in non-central collisions. Both observations pose challenges to the theoretical understanding of the coupling between heavy quarks and the quark-gluon plasma produced in these collisions. Transport models for the evolution of heavy quarks in a QCD medium offer the opportunity to study these effects -two of the most successful approaches are based on the linearized Boltzmann transport equation and the Langevin equation. In this work, we develop a hybrid transport model that combines the strengths of both of these approaches: heavy quarks scatter with medium partons using matrix-elements calculated in perturbative QCD, while between these discrete hard scatterings they evolve using a Langevin equation with empirical transport coefficients to capture the non-perturbative soft part of the interaction. With the hybrid transport model coupled to a state-of-the-art event-by-event bulk evolution model based on 2+1D relativistic viscous fluid dynamics, we study the azimuthal anisotropy and nuclear modification factor of heavy quarks in Pb+Pb collisions at √ s = 5.02 TeV. The parameters of our model are calibrated using a Bayesian analysis comparing to available D-meson and B-meson data at the LHC. Using the calibrated model, we study the implications on heavy-flavor transport properties and predict novel observables. arXiv:1806.08848v1 [nucl-th]

Several transport models have been employed in recent years to analyze heavy-flavor meson spectra in high-energy heavy-ion collisions. Heavy-quark transport coefficients extracted from these models with their default parameters vary, however, by up to a factor of 5 at high momenta. To investigate the origin of this large theoretical uncertainty, a systematic comparison of heavy-quark transport coefficients is carried out between various transport models. Within a common scheme devised for the nuclear modification factor of charm quarks in a brick medium of a quark-gluon plasma, the systematic uncertainty of the extracted drag coefficient among these models is shown to be reduced to a factor of 2, which can be viewed as the smallest intrinsic systematical error band achievable at present time. This indicates the importance of a realistic hydrodynamic evolution constrained by bulk hadron spectra and of heavy-quark hadronization for understanding the final heavy-flavor hadron spectra and extracting heavy-quark drag coefficient. The transverse transport coefficient is less constrained due to the influence of the underlying mechanism for heavy-quark medium interaction. Additional constraints on transport models such as energy loss fluctuation and transverse-momentum broadening can further reduce theoretical uncertainties in the extracted transport coefficients.

No abstract