Applying Bayesian parameter estimation to relativistic heavy-ion collisions: simultaneous characterization of the initial state and quark-gluon plasma medium

Abstract: 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, TRENTO, 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 paper applies Bayesian estimation to accurately constrain both the initial state and medium properties of the quark-gluon plasma in heavy-ion collisions.
• It reveals that entropy deposition scales with the geometric mean of participant densities and determines a nucleon width of approximately 0.5 fm.
• The study employs a comprehensive calibration using Gaussian process emulation and MCMC sampling, achieving experimental replication within 10% accuracy.

An Exploration of Bayesian Parameter Estimation in Relativistic Heavy-Ion Collisions

The study at hand presents a rigorous investigation into the characterization of quark-gluon plasma (QGP) properties formed in relativistic heavy-ion collisions via Bayesian parameter estimation. Utilizing a model-to-data comparison paradigm, this paper addresses a multi-faceted modeling approach incorporating a parametric initial condition model alongside a coupled viscous hydrodynamics and hadronic cascade framework.

Core Methodology

The authors advance the field of relativistic heavy-ion collision modeling by implementing Bayesian statistical methods to estimate critical QGP properties. The Bayesian framework relies on the careful calibration of a comprehensive parametric model with experimental data from high-energy collisions observed at facilities like the RHIC and LHC. The sophisticated parametric model interpolates a broad class of particle production schemes, enhancing the ability to constrain initial conditions and temperature-dependent transport coefficients.

A key attribute of the study is the flexibility of the parametric initial condition model, known as the -0.5ex\ model, which interpolates various theoretical models characteristic of initial QGP conditions. Such versatility allows for the effective reproduction of different entropy deposition schemes by varying a single entropy deposition parameter. The integration with viscous hydrodynamics facilitates a nuanced description of QGP dynamics, transitioning seamlessly into microscopic hadronic transport models for the post-QGP evolution.

Major Findings and Implications

The Bayesian analysis within the paper reveals several pertinent outcomes:

1. Initial Condition Characterization: The entropy deposition is found to be proportional to the geometric mean of participant nuclear densities, corresponding to a -0.5ex\ parameter $p$ approximately equal to zero. This suggests compatibility with established models such as EKRT and IP-Glasma.
2. Transport Coefficients: A relation between the shear viscosity minimum value and temperature slope is discerned, with uncertainty minimized at intermediate temperatures (~200-225 MeV). There is substantial evidence for a non-zero bulk viscosity, indicating its necessity for modeling transverse momentum and flow dynamics accurately.
3. Nucleon Width and Fluctuations: The Gaussian nucleon width is robustly estimated to be around 0.5 fm, aligning well with prior experimental and theoretical studies. However, multiplicity fluctuations play a limited role in Pb+Pb collisions due to larger position fluctuations.
4. Particlization Temperature: Significant divergence in posterior distributions when using identified versus charged particle yields highlights the influence of hadron chemistry in hybrid model fits.

Model Calibration and Verification

Utilizing a Gaussian process emulator interpolated from extensive model evaluations, the paper performs a comprehensive parameter space exploration via Markov chain Monte Carlo sampling. The calibration is performed against identified and charged particle yields, discerning the model's proficiency in replicating experimental observables within a 10% accuracy for most metrics.

Implications for Future Research

This Bayesian approach provides a structured pathway for extracting QGP properties by integrating comprehensive initial condition models with hydrodynamic simulation frameworks. As computational resources become increasingly available, coupling data from varied collision systems and energies could further illuminate the temperature dependence of transport coefficients. Nonetheless, the study also underscores the need for enhanced particle production modeling to reconcile discrepancies in hadronic chemical outputs (e.g., pion/kaon ratios).

The methodologies and insights delineated in this paper pave the way for more precise theoretical models of QGP dynamics, setting a benchmark for subsequent studies aiming to unravel the complexities of the QCD matter state characterized post-collision. As such, the fusion of advanced statistical methods with high-fidelity QCD modeling stands as a promising frontier for theoretical and computational high-energy nuclear physics.