We quantitatively estimate properties of the quark-gluon plasma created in ultra-relativistic heavy-ion collisions utilizing Bayesian statistics and a multi-parameter model-to-data comparison. The study is performed using a recently developed parametric initial condition model, T R ENTo, which interpolates among a general class of particle production schemes, and a modern hybrid model which couples viscous hydrodynamics to a hadronic cascade. We calibrate the model to multiplicity, transverse momentum, and flow data and report constraints on the parametrized initial conditions and the temperature-dependent transport coefficients of the quark-gluon plasma. We show that initial entropy deposition is consistent with a saturation-based picture, extract a relation between the minimum value and slope of the temperature-dependent specific shear viscosity, and find a clear signal for a nonzero bulk viscosity.

The iEBE-VISHNU code package performs event-by-event simulations for relativistic heavy-ion collisions using a hybrid approach based on (2+1)-dimensional viscous hydrodynamics coupled to a hadronic cascade model. We present the detailed model implementation, accompanied by some numerical code tests for the package. iEBE-VISHNU forms the core of a general theoretical framework for model-data comparisons through large scale Monte-Carlo simulations. A numerical interface between the hydrodynamically evolving medium and thermal photon radiation is also discussed. This interface is more generally designed for calculations of all kinds of rare probes that are coupled to the temperature and flow velocity evolution of the bulk medium, such as jet energy loss and heavy quark diffusion. It is impossible to use external probes to study the properties of the quark-gluon plasma (QGP), a novel state of matter created during the collisions. Experiments can only measure the momentum information of stable hadrons, who are the remnants of the collisions. In order to extract the thermal and transport properties of the QGP, one needs to rely on Monte-Carlo event-by-event model simulations, which reverse-engineer the experimental measurements to the early time dynamics of the relativistic heavy-ion collisions. Solution method: Relativistic heavy-ion collisions contain multiple stages of evolution. The physics that governs each stage is implemented into individual code component. A general driver script glues all the modular packages as a whole to perform large-scale Monte-Carlo simulations. The final results are stored into SQLite database, which supports standard querying for massive data analysis. By tuning transport coefficients of the QGP as free parameters, e.g. the specific shear viscosity η/s, we can constrain various transport properties of the QGP through model-data comparisons. Keywords Running time:The following running time is tested on a laptop computer with a 2.4 GHz Intel Core i5 CPU, 4GB memory. All the C++ and Fortran codes are compiled with the GNU Compiler Collection (GCC) 4.9.2 and -O3 optimization.

No abstract

We introduce T R ENTo, a new parametric initial condition model for high-energy nuclear collisions based on eikonal entropy deposition via a "reduced thickness" function. The model simultaneously describes experimental proton-proton, proton-nucleus, and nucleus-nucleus multiplicity distributions, and generates nucleus-nucleus eccentricity harmonics consistent with experimental flow constraints. In addition, the model is compatible with ultra-central uranium-uranium data unlike existing models that include binary collision terms.

I develop and apply a Bayesian method for quantitatively estimating properties of the quark-gluon plasma (QGP), an extremely hot and dense state of fluid-like matter created in relativistic heavy-ion collisions.The QGP cannot be directly observed-it is extraordinarily tiny and ephemeral, about 10 −14 meters in size and living 10 −23 seconds before freezing into discrete particles-but it can be indirectly characterized by matching the output of a computational collision model to experimental observations. The model, which takes the QGP properties of interest as input parameters, is calibrated to fit the experimental data, thereby extracting a posterior probability distribution for the parameters.In this dissertation, I construct a specific computational model of heavyion collisions and formulate the Bayesian parameter estimation method, which is based on general statistical techniques. I then apply these tools to estimate fundamental QGP properties, including its key transport coefficients and characteristics of the initial state of heavy-ion collisions.Perhaps most notably, I report the most precise estimate to date of the temperature-dependent specific shear viscosity η/s, the measurement of which is a primary goal of heavy-ion physics. The estimated minimum value is η/s = 0.085 +0.026 −0.025 (posterior median and 90% uncertainty), remarkably close to the conjectured lower bound of 1/4π 0.08. The analysis also shows that η/s likely increases slowly as a function of temperature.Other estimated quantities include the temperature-dependent bulk viscosity ζ/s, the scaling of initial state entropy deposition, and the duration of the pre-equilibrium stage that precedes QGP formation.

We posit a unified hydrodynamic and microscopic description of the quark-gluon plasma (QGP) produced in ultrarelativistic p-Pb and Pb-Pb collisions at √ sNN = 5.02 TeV and evaluate our assertion using Bayesian inference. Specifically, we model the dynamics of both collision systems using initial conditions with parametric nucleon substructure, a pre-equilibrium free streaming stage, event-by-event viscous hydrodynamics, and a microscopic hadronic afterburner. Free parameters of the model which describe the initial state and QGP medium are then simultaneously calibrated to fit charged-particle yields, mean pT , and flow cumulants. We argue that the global agreement of the calibrated model with the experimental data strongly supports the existence of hydrodynamic flow in small collision systems at ultrarelativistic energies, and that the flow produced develops at length scales smaller than a single proton. Posterior estimates for the model's input parameters are obtained, and new insights into the temperature dependence of the QGP transport coefficients and event-by-event structure of the proton are discussed.

International audienceBy applying a Bayesian model-to-data analysis, we estimate the temperature and momentum dependence of the heavy quark diffusion coefficient in an improved Langevin framework. The posterior range of the diffusion coefficient is obtained by performing a Markov chain Monte Carlo random walk and calibrating on the experimental data of D-meson RAA and v2 in three different collision systems at the Relativistic Heavy-Ion Collidaer (RHIC) and the Large Hadron Collider (LHC): Au-Au collisions at 200 GeV and Pb-Pb collisions at 2.76 and 5.02 TeV. The spatial diffusion coefficient is found to be consistent with lattice QCD calculations and comparable with other models' estimation. We demonstrate the capability of our improved Langevin model to simultaneously describe the RAA and v2 at both RHIC and the LHC energies, as well as the higher order flow coefficient such as D meson v3. We show that by applying a Bayesian analysis, we are able to quantitatively and systematically study the heavy flavor dynamics in heavy-ion collisions

The iEBE-VISHNU code package performs event-by-event simulations for relativistic heavy-ion collisions using a hybrid approach based on (2+1)-dimensional viscous hydrodynamics coupled to a hadronic cascade model. We present the detailed model implementation, accompanied by some numerical code tests for the package. iEBE-VISHNU forms the core of a general theoretical framework for model-data comparisons through large scale Monte-Carlo simulations. A numerical interface between the hydrodynamically evolving medium and thermal photon radiation is also discussed. This interface is more generally designed for calculations of all kinds of rare probes that are coupled to the temperature and flow velocity evolution of the bulk medium, such as jet energy loss and heavy quark diffusion.