---
title: Quantum van der Waals Hadron Resonance Gas
url: https://www.emergentmind.com/topics/quantum-van-der-waals-hadron-resonance-gas-qvdw-hrg
type: topic
---

# Quantum van der Waals Hadron Resonance Gas

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 [1707.09215].

## 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^*\}$ solve the coupled, transcendental system
$$
\mu_i^* = \mu_i - \sum_{j}\tilde{b}_{ij}\,P_j^{\mathrm{id}}(T, \mu_j^*) + \sum_{j}(a_{ij}+a_{ji}) n_j,
$$
with $\tilde{b}_{ij}$ the (possibly “mixed”) excluded-volume coefficient, often set to $b_{ij} = \frac{2\pi}{3}(r_i + r_j)^3$ [1707.09215].

The grand potential is $\Omega(T,\{\mu_i\}) = -V P(T,\{\mu_i\})$.

## 2. Parameterization of Repulsive and Attractive Forces

- **Repulsive ('excluded volume') matrix**: $b_{ij}$, derived from hard-sphere radii. In practice, one often uses a common value for all (anti)baryons, e.g., $b_{\rm NS} = 3.42\,\text{fm}^3$ for non-strange baryons, with smaller values for strange baryons (e.g., $b_{\rm S} = b_{\rm NS}/8$), motivated by nuclear matter fits and lattice constraints.
- **Attractive matrix**: $a_{ij}$, symmetric, typically constructed via $a_{ij} = \sqrt{a_{ii} a_{jj}}$, with typical values $a_{\rm NS} = 329\,\text{MeV}\,\text{fm}^3$, $a_{\rm S} = a_{\rm NS}/8$.
- For mixtures, $\tilde b_{ij}$ may be replaced by $b_{ij}$ for simplicity; cross-terms are defined as geometric means for $a_{ij}$ and functions of component radii for $b_{ij}$.

These parameters can be refined for different baryon sectors and explicitly tuned using lattice QCD EoS benchmarks or nuclear matter properties [1707.09215].

## 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 ($a_{ij}=b_{ij}=0$ 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 [1711.09863]. Resonance decays are usually only included when computing thermal yields, not for bulk thermodynamic quantities [1707.09215].

## 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 $T > 50$ MeV and sub-threshold chemical potentials).

Thermodynamic quantities such as entropy density and energy density are given by
$$
S/V = \sum_i f_i s_i^{\rm id}(T,\mu_i^*), \qquad \varepsilon = \sum_i f_i \varepsilon_i^{\rm id}(T,\mu_i^*) - \sum_{ij}a_{ij} n_i n_j,
$$
with reduction factors $f_i = 1 - \sum_j \tilde{b}_{ji} n_j$ denoting the effective packing fraction.

Maxwell constructions in $P(V)$ or $T(\mu)$ are used to identify first-order phase transitions and critical points, yielding mean-field van der Waals exponents and structures [1707.09215].

## 5. Comparison with Lattice QCD and Experimental Phenomenology

- **Reproduction of the nuclear liquid–gas transition**: At $T=0$ and the nuclear saturation point, the model yields a first-order transition ending at a critical point, e.g., $(T_c \approx 19.7\,\text{MeV}, n_c \approx 0.072\,\text{fm}^{-3})$ [1707.09215].
- **Crossover thermodynamics**: QvdW-HRG accurately describes baryon number susceptibilities $\chi_2^B(T)$, with enhanced agreement above $T \sim 130\,\text{MeV}$ compared to the ideal HRG, and characteristic suppression of higher-order cumulants $\chi_4^B/\chi_2^B$ in the range $T \sim 150$–$165\,\text{MeV}$, aligning with lattice QCD [1707.09215]. 
- **Strangeness and off-diagonal correlators**: The baryon–strangeness correlator $C_{BS}$ and net-strangeness susceptibility $\chi_2^S$ are sensitive to the size of $a_S, b_S$; taking these smaller than their non-strange counterparts aligns model predictions with lattice observations, though residual discrepancies for $C_{BS}$ indicate the likely importance of missing high-mass strange states.
- **Energy and entropy densities**: QvdW-HRG tracks lattice QCD up to $T \sim 160$ MeV.

## 6. Numerical Procedure and Model Variants

Numerical implementation requires iterative solution of the transcendental system for $\{\mu_i^*\}$ and $\{n_i\}$ at each $(T,\{\mu_i\})$. The model admits several extensions:
- Multi-component generalizations with full matrices $a_{ij}, b_{ij}$ [1707.09215].
- Inclusion of resonance widths and additional predicted (Hagedorn or quark model) states, which modify fit values for $a,b$ and shift the phase diagram [2304.11914].
- Induced surface tension (IST) variants that introduce an extra pressure equation for surface contributions, yielding improved reproduction of quantum virial coefficients to higher order [1704.06846].

## 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 $\chi_2^B$, $\chi_4^B/\chi_2^B$).
- Reproduces nuclear matter saturation properties without additional tuning when using empirical $a$, $b$.
- 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 [1707.09215, 1711.09863, 2304.11914].

Source: https://www.emergentmind.com/topics/quantum-van-der-waals-hadron-resonance-gas-qvdw-hrg