Nonequilibrium Phenomenology of Identified Particle Spectra in Heavy-Ion Collisions at LHC Energies

Technical Summary: Nonequilibrium Phenomenology of Identified Particle Spectra in Heavy-Ion Collisions at LHC Energies

Problem Statement
The paper addresses a persistent discrepancy in the theoretical description of identified particle spectra in ultrarelativistic heavy-ion collisions (HIC) at CERN LHC energies. While statistical models successfully describe particle yield ratios via a chemical freeze-out scenario, they fail to accurately reproduce the transverse-momentum ($p_T$) spectra of pions and kaons at low momenta. Specifically, standard thermal distribution functions and conventional hydrodynamic codes systematically underestimate the pion yield in the low-$p_T$ region. Previous explanations attributed this enhancement to resonance decays and radial flow, but recent high-precision data from the ALICE Collaboration indicate that these mechanisms are insufficient to describe the observed overpopulation of low-energy states.

Methodology
The authors employ the Zubarev method of the Nonequilibrium Statistical Operator (NSO) to develop a systematic approach for describing non-equilibrium pion production. The core of the methodology involves:

1. Theoretical Framework: Utilizing the NSO formalism to derive distribution functions where long-living fluctuations (specifically pion and kaon numbers, $N_\pi$ and $N_K$) are treated as relevant observables. This leads to the introduction of effective, non-equilibrium chemical potentials ($\mu_\pi, \mu_K$) in the distribution function, distinct from equilibrium parameters.
2. Modeling Approaches: The study tests four distinct models (A–D) to describe particle freeze-out:
   • Model A: A standard blast-wave model without resonance decays or final-state interactions.
   • Model B: A blast-wave model incorporating non-equilibrium effective chemical potentials for pions and kaons within the Bose distribution.
   • Models C & D: A hybrid approach coupling a blast-wave-based thermal particle generator to the SMASH hadronic transport model. Model C assumes full chemical equilibrium, while Model D includes an effective pion chemical potential to account for overpopulated states.
3. Data Analysis: The models are fitted to high-precision transverse-momentum spectra of identified particles ($\pi^\pm, K^\pm, p/\bar{p}$) measured by ALICE in central Pb+Pb collisions at $\sqrt{s_{NN}} = 2.76$ TeV.
4. Statistical Inference: Bayesian inference methods, utilizing Gaussian process emulators and Markov chain Monte Carlo (MCMC) sampling, are applied to extract the most probable thermodynamic parameters (temperature, flow velocity, chemical potentials) and evaluate model probabilities via Bayes factors.

Key Results

• Model A (Standard): Replicates previous findings, yielding a freeze-out temperature of $\sim 97.7$ MeV and surface velocity $\sim 0.88c$. It significantly deviates from ALICE data at low $p_T$.
• Model B (Non-equilibrium Chemical Potential): The introduction of effective chemical potentials ($\mu_\pi \approx 126$ MeV, $\mu_K \approx 403$ MeV) allows the model to describe the low-$p_T$ enhancement without resonance decays. This model yields a higher temperature ($\sim 118$ MeV) and provides a significantly better fit ($\chi^2$ improvement by a factor of 6.7) over a broad $p_T$ range, including regions outside the fitted intervals.
• Models C & D (Transport Coupled): When resonance decays and final-state interactions are explicitly simulated via SMASH, the resulting temperature is $\sim 156$ MeV, consistent with Lattice QCD and statistical hadronization models.
  • Model C (equilibrium) describes the data well.
  • Model D (non-equilibrium $\mu_\pi$) yields a $\mu_\pi$ compatible with zero and offers no significant improvement in fit quality over Model C. The Bayes factor ($B^D_C \approx 0.43$) indicates that the data does not strongly favor the non-equilibrium extension when explicit resonance decays are included.

Key Contributions

1. Systematic NSO Application: The paper demonstrates the application of the Zubarev NSO method to HIC phenomenology, formally justifying the introduction of effective chemical potentials as Lagrange multipliers for conserved particle numbers in a non-equilibrium setting.
2. Alternative Explanation for Low-$p_T$ Enhancement: The study establishes that a non-equilibrium pion chemical potential (Model B) offers an efficient alternative to the conventional explanation of low-$p_T$ enhancement via resonance decays. It successfully describes the overpopulation of low-momentum states using a thermal framework extended by non-equilibrium parameters.
3. Validation of Resonance Decays: By coupling the blast-wave model to the SMASH transport code, the authors confirm that explicit treatment of resonance decays and thermalization (Model C) naturally accounts for the low-$p_T$ enhancement, rendering the explicit non-equilibrium chemical potential (Model D) redundant in that specific context.

Significance and Claims
The authors claim that their work provides a more fundamental approach to the problem of low-momentum pion enhancement, moving beyond phenomenological fixes. They demonstrate that the "overpopulation" of low-$p_T$ states can be described either by:

• Introducing a non-equilibrium chemical potential in a simplified thermal model (Model B), which serves as an efficient alternative to tracking microscopic thermalization processes.
• Or, by explicitly simulating resonance decays and final-state interactions (Model C).

The paper concludes that while Model B is a successful phenomenological tool for describing kinetic freeze-out spectra, it is an approximation. The authors note that all models presented (particularly C and D) remain approximations because they do not yet fully incorporate in-medium effects (such as the Mott effect or continuum correlations) which are essential at high densities. The ultimate goal stated is to extend the NSO approach to include flow observables and resonance formation dynamics to provide a consistent scenario that justifies the success and limitations of thermal statistical models. The work suggests that fugacity parameters used in thermal models may be rigorously derived as Lagrange multipliers within the NSO formalism.