Hadron Resonance Gas Model and Extensions

• Hadron Resonance Gas model is a statistical framework that describes the thermodynamics of strongly interacting hadronic matter using a partition sum over known resonances.
• It incorporates quantum statistical corrections, two-body correlations via the Beth–Uhlenbeck approach, and repulsive effects like excluded volume and van der Waals interactions to refine thermodynamic predictions.
• Precise lattice QCD fits reveal both the model's success in capturing equilibrium properties and its limitations at higher temperatures and baryochemical potentials, spurring further research on multi-body interactions.

The hadron resonance gas (HRG) model provides a statistical framework for describing the thermodynamics of strongly interacting matter in the hadronic phase of quantum chromodynamics (QCD). At temperatures and baryochemical potentials below the deconfinement transition, and up to the vicinity of the chiral crossover, the equation of state and related observables are accurately described by a partition sum over all known hadrons and resonances, treated as non-interacting or weakly interacting species. Modern extensions of the HRG include quantum statistical corrections, explicit two-body correlations, van der Waals–type repulsive and attractive forces, surface tension effects, and empirical constraints from lattice QCD. These developments result in a highly flexible and physically grounded model, relevant for interpreting equilibrium properties and fluctuations in heavy-ion collisions, as well as for benchmarking QCD thermodynamics against lattice results.

1. Theoretical Foundations and Partition Function

The foundational construct is the grand-canonical partition function, which for the ideal HRG reads

$$\ln Z(T,V,\{\mu_i\}) = \sum_i \pm \frac{Vg_i}{2\pi^2} \int_0^\infty p^2 dp \ln\bigl[ 1 \pm \lambda_i e^{-E_i(p)/T} \bigr],$$

where $g_i$ is the degeneracy, $E_i(p) = \sqrt{p^2 + m_i^2}$, $$\lambda_i = \exp\bigl[(B_i\mu_B + S_i\mu_S + Q_i\mu_Q)/T\bigr]$$, with signs for Bose/Fermi statistics. The summation extends over all relevant hadron and resonance species. For thermodynamic quantities, pressure is given by $P = (T/V)\ln Z$, with energy and number densities obtained via temperature and chemical potential derivatives. The inclusion of quantum statistics becomes essential below $\sim 150$ MeV but can be neglected at higher temperatures.

2. Extensions: Interactions Beyond the Ideal Gas

2.1 Quantum Statistical Correlations: Beth–Uhlenbeck Approach

Quantum interactions enter via the Beth–Uhlenbeck (BU) formalism, which expresses the second virial coefficient as an integral over two-body scattering phase shifts,

$$b_2(T) \propto \sum_\ell(2\ell+1)\int dE\, e^{-E/T} \frac{d\delta_\ell(E)}{dE}.$$

In the Boltzmann limit, pair correlations introduce a multiplicative modification to the occupation probability, parameterized by a correlation length $r$, such that a species $i$ receives a factor $1\pm e^{-4\pi^2 m_i T r^2}$ (plus for bosons, minus for fermions), and all thermodynamic functions become $r$-dependent (Hanafy et al., 2020). This construction encodes both attractive and repulsive quantum mechanical effects in a minimal form.

2.2 Excluded Volume and van der Waals Interactions

Repulsion at short distances is traditionally modeled by an excluded volume (EV) prescription, shifting the chemical potentials,

$\mu_i^* = \mu_i - v_i P,$

with eigenvolume $v_i$ for each species, leading to pressure and densities subject to transcendental equations. The van der Waals (vdW) approach incorporates both repulsion and mean-field attraction (for baryons),

$$p_B = \sum_{i\in B}p_i^{\rm id}(T,\mu_i^*) - a n_B^2, \quad \mu_i^* = \mu_i - b p_B - a b n_B^2 + 2a n_B,$$

with empirical parameters $a$ and $b$ fixed to nuclear matter properties (Vovchenko et al., 2017, Vovchenko, 2020).

2.3 Surface Tension and Multicomponent Generalizations

The induced surface tension (IST) model further refines the treatment of repulsive forces by introducing a surface tension coefficient $\Sigma$, yielding coupled equations for pressure and $\Sigma$, and capturing higher-order virial contributions accurately. This approach produces a causal EoS up to higher packing fractions compared to vdW/EVM schemes and manages arbitrarily many hadron radii with only two equations (Sagun et al., 2017, Bugaev et al., 2016).

3. Correlations, Fluctuations, and Constraints from Lattice QCD

Generalized susceptibilities (cumulants) of conserved charges ($\chi_n^B,\,\chi_n^Q,\,\chi_n^S$) are accessible in the HRG via analytic derivatives of the pressure with respect to chemical potentials, and can be compared directly to lattice QCD and experiment. Extensions to HRG, including BU correlations, EV/vdW corrections, and additional "missing" resonances, are empirically constrained by systematic comparison with lattice data at both $\mu_B=0$ and finite baryochemical potential.

• The inclusion of a finite correlation length $r$ (BU-type) improves agreement with lattice thermodynamics at moderate $\mu_B$, but ideal HRG ($r=0$) remains optimal at $\mu_B=0$ (Hanafy et al., 2020).
• vdW interactions systematically improve susceptibilities such as $\chi_2^B$, $\chi_{11}^{BQ}$ (Vovchenko et al., 2017).
• Adding unconfirmed (QM/PDG+) baryon resonances further adjusts fluctuations and helps achieve better correspondence for $\chi_{11}^{BQ}/\chi_2^B$, but some tension persists, especially in ratios where vdW corrections cancel analytically (Vovchenko et al., 2017, Karthein et al., 2021).
• At high T or $\mu_B$, all classical and quantum extended HRG variants underpredict the rapid rise of pressure and energy above $T_c$, indicating the breakdown of the hadronic description and the onset of deconfinement.

4. Numerical Fits and Parameter Dependence

The efficacy of each HRG extension is quantified by fits ($\chi^2/$dof) to lattice QCD observables across a grid of $T$ and $\mu_B$:

μ_B (MeV)	Optimal r (fm)	χ²/dof (P/T⁴)	Notes
0	0	~0.02	Ideal HRG
170	0–0.05	0.002–0.01	BU improves fit
340	0.1–0.15	~0.005	Larger r optimal
425	—	—	No convincing fit

The best-fit correlation length increases with baryon density, suggesting increased quantum correlations in baryon-rich environments. However, at high temperature above $T \approx 170$ MeV, all BU/EV/vdW/IST extensions fail to match the sharp rise of lattice observables, consistent with the emergence of partonic degrees of freedom (Hanafy et al., 2020).

5. Physical Implications for Strongly Interacting Matter

The extended HRG paradigm robustly supports the resonance dominance of equilibrium thermodynamics up to the chiral crossover. The success of adding two-body quantum correlations, surface tension corrections, and flavor-dependent repulsion underlines the importance of both repulsive and attractive forces in the hot hadronic phase. Quantum and mean-field repulsions become increasingly relevant at finite baryon densities, while the necessity of an ever-richer hadron spectrum hints at incipient deconfinement. Systematic studies of fluctuations, higher-order susceptibilities, and fits to precise lattice data continue to test the limits and physical content of the HRG framework.

6. Limitations and Outlook

Contemporary HRG models employing BU corrections or van der Waals/IST features are limited by several factors:

• Only two-body virial corrections are included; three-body (or higher) cluster effects are neglected.
• No explicit use of S-matrix phase shifts in BU models, precluding detailed separation of dynamical resonance dynamics, channel coupling, or a systematic treatment of repulsive/attractive forces via empirical scattering data.
• Validity is typically limited to $T \lesssim 200$ MeV, $\mu_B \lesssim 340$ MeV; deviation from lattice data occurs beyond these domains as partonic degrees of freedom dominate.
• Full treatment of higher-order susceptibilities and explicit inclusion of empirical scattering phases is an active direction for future refinement (Hanafy et al., 2020).

Potential extensions include complete S-matrix implementations, density-dependent interaction parameters, systematic inclusion of three- and four-body interactions, and mapping to heavy-ion observables at lower beam energies and higher densities.

---

In summary, the hadron resonance gas model, and its quantum-statistical and repulsive-interaction extensions, form the foundation for understanding the equilibrium thermodynamics of QCD matter in heavy-ion phenomenology and for benchmarking the hadronic phase against first-principles lattice QCD (Hanafy et al., 2020, Vovchenko et al., 2017, Sagun et al., 2017, Karthein et al., 2021).