Quantum van der Waals Hadron Resonance Gas

• Quantum van der Waals Hadron Resonance Gas is a model that extends the ideal HRG by incorporating species-dependent repulsive (excluded volume) and attractive (mean-field) interactions.
• The model self-consistently solves a system of quantum-statistical equations to determine effective chemical potentials and particle densities, capturing key nuclear phase transitions.
• It accurately reproduces nuclear liquid–gas transitions and lattice QCD susceptibilities, providing insights into heavy-ion collision phenomenology and the QCD phase diagram.

The Quantum van der Waals Hadron Resonance Gas (QvdW-HRG) model is a thermodynamically consistent extension of the ideal hadron resonance gas framework, incorporating both short-range repulsive (excluded-volume) and intermediate-range attractive (mean-field) interactions among hadrons, specifically baryons and antibaryons. It utilizes quantum statistical mechanics for all species, parametrizes the repulsive and attractive forces by species-dependent matrices, and self-consistently determines the equation of state (EoS) for a multi-component system relevant for nuclear and hadronic physics (Vovchenko et al., 2017).

1. Quantum-Statistical Foundation and Core Equations

The QvdW-HRG pressure in the grand-canonical ensemble at temperature $T$ and chemical potentials $\{\mu_i\}$ is given by

$$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$

where $P_i^{\mathrm{id}}(T, \mu_i^*)$ is the quantum ideal-gas pressure of species $i$ (Fermi-Dirac for baryons and anti-baryons, Bose-Einstein for mesons), evaluated at an effective chemical potential $\mu_i^*$. The densities $n_i$ are also determined self-consistently via

$$n_i = \frac{n_i^{\mathrm{id}}(T, \mu_i^*)}{1 + \sum_j \tilde{b}_{ji} n_j^{\mathrm{id}}(T, \mu_j^*)},$$

where $n_i^{\mathrm{id}}$ is the ideal quantum-gas density, with upper/lower signs for fermions and bosons, respectively,

$$n_i^{\mathrm{id}}(T, \mu) = \frac{g_i}{2\pi^2} \int_0^\infty k^2 dk \: [\exp((\sqrt{k^2 + m_i^2} - \mu)/T)\pm1]^{-1}.$$

The effective chemical potentials $\{\mu_i\}$0 solve the coupled, transcendental system

$\{\mu_i\}$1

with $\{\mu_i\}$2 the (possibly “mixed”) excluded-volume coefficient, often set to $\{\mu_i\}$3 (Vovchenko et al., 2017).

The grand potential is $\{\mu_i\}$4.

2. Parameterization of Repulsive and Attractive Forces

• Repulsive ('excluded volume') matrix: $\{\mu_i\}$5, derived from hard-sphere radii. In practice, one often uses a common value for all (anti)baryons, e.g., $\{\mu_i\}$6 for non-strange baryons, with smaller values for strange baryons (e.g., $\{\mu_i\}$7), motivated by nuclear matter fits and lattice constraints.
• Attractive matrix: $\{\mu_i\}$8, symmetric, typically constructed via $\{\mu_i\}$9, with typical values $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$0, $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$1.
• For mixtures, $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$2 may be replaced by $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$3 for simplicity; cross-terms are defined as geometric means for $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$4 and functions of component radii for $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$5.

These parameters can be refined for different baryon sectors and explicitly tuned using lattice QCD EoS benchmarks or nuclear matter properties (Vovchenko et al., 2017).

3. Treatment of Resonances and Mesons

In the QvdW-HRG, all hadrons and resonances up to a cutoff mass (typically 2–2.5 GeV) are included as distinct, quantum-statistical species. The van der Waals interactions are implemented among baryon–baryon and antibaryon–antibaryon pairs; mesons and meson–baryon (or meson–antibaryon) cross-terms are typically neglected ($$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$6 for such pairs).

Resonance widths can be incorporated by folding spectral mass distributions (e.g., relativistic Breit–Wigner), especially when comparing to susceptibility data or for detailed fluctuation studies (Vovchenko et al., 2017). Resonance decays are usually only included when computing thermal yields, not for bulk thermodynamic quantities (Vovchenko et al., 2017).

4. Thermodynamic Consistency and Approximations

The model employs only two-body (mean-field) van der Waals terms; higher virial corrections beyond excluded volume are neglected. Bose condensation for mesons is commonly ignored (valid for $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$7 MeV and sub-threshold chemical potentials).

Thermodynamic quantities such as entropy density and energy density are given by

$$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$8

with reduction factors $$P(T,\{\mu_i\}) = \sum_{i=1}^N P_i^{\mathrm{id}}(T, \mu_i^*) - \sum_{i,j=1}^N a_{ij}\,n_i\,n_j,$$9 denoting the effective packing fraction.

Maxwell constructions in $P_i^{\mathrm{id}}(T, \mu_i^*)$0 or $P_i^{\mathrm{id}}(T, \mu_i^*)$1 are used to identify first-order phase transitions and critical points, yielding mean-field van der Waals exponents and structures (Vovchenko et al., 2017).

5. Comparison with Lattice QCD and Experimental Phenomenology

• Reproduction of the nuclear liquid–gas transition: At $P_i^{\mathrm{id}}(T, \mu_i^*)$2 and the nuclear saturation point, the model yields a first-order transition ending at a critical point, e.g., $P_i^{\mathrm{id}}(T, \mu_i^*)$3 (Vovchenko et al., 2017).
• Crossover thermodynamics: QvdW-HRG accurately describes baryon number susceptibilities $P_i^{\mathrm{id}}(T, \mu_i^*)$4, with enhanced agreement above $P_i^{\mathrm{id}}(T, \mu_i^*)$5 compared to the ideal HRG, and characteristic suppression of higher-order cumulants $P_i^{\mathrm{id}}(T, \mu_i^*)$6 in the range $P_i^{\mathrm{id}}(T, \mu_i^*)$7–$P_i^{\mathrm{id}}(T, \mu_i^*)$8, aligning with lattice QCD (Vovchenko et al., 2017).
• Strangeness and off-diagonal correlators: The baryon–strangeness correlator $P_i^{\mathrm{id}}(T, \mu_i^*)$9 and net-strangeness susceptibility $i$0 are sensitive to the size of $i$1; taking these smaller than their non-strange counterparts aligns model predictions with lattice observations, though residual discrepancies for $i$2 indicate the likely importance of missing high-mass strange states.
• Energy and entropy densities: QvdW-HRG tracks lattice QCD up to $i$3 MeV.

6. Numerical Procedure and Model Variants

Numerical implementation requires iterative solution of the transcendental system for $i$4 and $i$5 at each $i$6. The model admits several extensions:

• Multi-component generalizations with full matrices $i$7 (Vovchenko et al., 2017).
• Inclusion of resonance widths and additional predicted (Hagedorn or quark model) states, which modify fit values for $i$8 and shift the phase diagram (Sarkar, 2023).
• Induced surface tension (IST) variants that introduce an extra pressure equation for surface contributions, yielding improved reproduction of quantum virial coefficients to higher order (Bugaev et al., 2017).

7. Principal Results and Physical Significance

The QvdW-HRG provides a minimal yet powerful extension to the ideal HRG that:

• Dynamically generates a first-order nuclear liquid–gas transition and associated critical exponents.
• Suppresses baryonic fluctuations, resulting in improved matching to lattice QCD susceptibilities (particularly for $i$9, $\mu_i^*$0).
• Reproduces nuclear matter saturation properties without additional tuning when using empirical $\mu_i^*$1, $\mu_i^*$2.
• Offers a framework to analyze the onset of non-ideal hadronic effects in heavy-ion collisions and to study the approach to the QCD critical (endpoint) region.
• Allows quantification of the impact of strangeness, extra resonances, and further refinements via a straightforward prescription.

The model's success lies in its ability to unify the nuclear matter equation of state with lattice QCD and experimental heavy-ion phenomenology under a single, quantum-statistically consistent formalism with a minimal set of effective interaction parameters (Vovchenko et al., 2017, Vovchenko et al., 2017, Sarkar, 2023).