We study the dynamics of energy loss and flow of heavy quarks produced in ultra-relativistic heavy-ion collisions within the framework of a Langevin equation coupled to a (2+1)-dimensional viscous hydrodynamic model that simulates the space-time evolution of the produced hot and dense QCD matter. The classical Langevin approach is improved such that, apart from quasi-elastic scatterings, radiative energy loss is incorporated by treating gluon radiation as an additional force term. The hadronization of emitted heavy quarks is simulated via a hybrid fragmentation plus recombination model. Our calculation shows significant contribution from gluon radiation to heavy quark energy loss at high energies, and we find the recombination mechanism is important for heavy flavor meson production at intermediate energies. We present numerical results for the nuclear modification and elliptic flow of D mesons, which are consistent with measurements at both LHC and RHIC; predictions for B mesons are also provided.

A Linearized Boltzmann Transport (LBT) model coupled with hydrodynamical background is established to describe the evolution of jet shower partons and medium excitations in high energy heavy-ion collisions. We extend the LBT model to include both elastic and inelastic processes for light and heavy partons in the quark-gluon plasma. A hybrid model of fragmentation and coalescence is developed for the hadronization of heavy quarks. Within this framework, we investigate how heavy flavor observables depend on various ingredients, such as different energy loss and hadronization mechanisms, the momentum and temperature dependences of the transport coefficients, and the radial flow of the expanding fireball. Our model calculations show good descriptions of the D meson suppression and elliptic flow observed at the LHC and RHIC. The prediction for the Pb-Pb collisions at √ sNN=5.02 TeV is provided.

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

We construct a theoretical framework to describe the evolution of heavy flavors produced in relativistic heavy-ion collisions. The in-medium energy loss of heavy quarks is described using our modified Langevin equation that incorporates both quasi-elastic scatterings and the medium-induced gluon radiation. The space-time profiles of the fireball is described by a (2+1)-dimensional hydrodynamics simulation. A hybrid model of fragmentation and coalescence is utilized for heavy quark hadronization, after which the produced heavy mesons together with the soft hadrons produced from the bulk QGP are fed into the hadron cascade UrQMD model to simulate the subsequent hadronic interactions. We find that the medium-induced gluon radiation contributes significantly to heavy quark energy loss at high pT; heavy-light quark coalescence enhances heavy meson production at intermediate pT; and scatterings inside the hadron gas further suppress the D meson RAA at large pT and enhance its v2. Our calculations provide good descriptions of heavy meson suppression and elliptic flow observed at both the LHC and RHIC.