Thermodynamically consistent treatment of repulsive corrections in HRG

This paper presents a thermodynamically consistent reformulation of the hadron resonance gas model using excluded-volume corrections and a liquid-drop inspired mass-radius relation, which successfully reproduces lattice QCD results for conserved-charge susceptibilities with only two adjustable parameters.

The Gist

The Big Picture: A Cosmic Dance Floor

Imagine the universe just after the Big Bang, or the very core of a neutron star. It's a place so hot and dense that it's like a massive, chaotic dance floor packed with billions of particles (hadrons) bumping into each other.

Physicists want to understand how this "soup" behaves. They use a model called the Hadron Resonance Gas (HRG). Think of this model as a way to predict how crowded the dance floor is and how much energy the dancers have.

However, there's a problem. In this model, the particles aren't just tiny, invisible points; they have size. They take up space. When you pack too many people into a room, they start pushing against each other. This is called repulsion.

The Problem: The "Pushy" Math Glitch

For a long time, scientists tried to account for this "pushing" by simply saying, "Okay, let's pretend the chemical potential (a fancy term for the 'desire' of particles to be in the system) changes based on how crowded it is."

But this approach had a glitch. It was like trying to calculate the total weight of a crowd by weighing everyone individually, but then accidentally adding the weight of the floor twice because the math got confused about who was pushing whom. This led to inconsistent results, especially when trying to predict complex behaviors like how the crowd fluctuates (susceptibilities). The math was "thermodynamically inconsistent"—it broke the fundamental rules of energy and heat.

The Solution: A New Way to Count

Somenath Pal proposes a clever fix. Instead of trying to calculate the unique "push" for every single type of particle (which is a nightmare because there are hundreds of different hadrons), he suggests a simplified, auxiliary view.

The Analogy: The "Universal Discount" vs. The "Crowded Room"

Imagine a concert where everyone wants to get in.

• The Old Way (Model II): You try to calculate exactly how much the person in the front row pushes the person behind them, and how the person in the back pushes the one in front. It's a chaotic, individual calculation for every single person. It gets messy and leads to errors.
• The New Way (Model I): Instead of tracking every shove, we imagine a "Universal Discount" (an energy shift) applied to everyone's ticket price equally.

How do we figure out what that discount should be? We use a simple rule: The total number of people in the room must stay the same.

1. We look at the real, messy quantum world where particles push each other.
2. We create a "fake" classical world where everyone is just a simple ball, but we adjust the "ticket price" (energy) for everyone by the same amount.
3. We tweak that adjustment until the number of people in the "fake" room matches the "real" room.

Once the numbers match, we use this simple "fake" world to do the heavy math. Because the math is simpler, it doesn't break the rules of thermodynamics. It's like using a smooth, flat map to navigate a mountainous terrain; the map isn't the terrain, but if you calibrate it right, it tells you exactly where you are without getting lost.

The "Liquid Drop" Trick: Guessing the Size of Particles

To make this work, the author needs to know how big the particles are. But nobody knows the exact size of every single hadron (like a proton or a neutron).

The Analogy: The Water Balloon

The author uses a "Liquid Drop" model. Imagine a water balloon.

• If you have a tiny drop of water, it's small.
• If you have a giant water balloon, it's big.
• The paper suggests that the size of a particle is related to its mass (how heavy it is) in a specific way, just like how the volume of a water balloon scales with the amount of water inside.

By assuming particles are like little liquid drops, the author creates a formula where the size of a heavy particle is just a scaled-up version of the size of a light particle (the pion). This reduces the problem to just two adjustable knobs:

1. The size of the smallest particle (the pion).
2. A scaling factor (how fast the size grows as mass increases).

The Results: A Perfect Match

The author tested this new method against Lattice QCD, which is like the "gold standard" supercomputer simulation of the strong force.

• The Test: They looked at how the "fluctuations" (wiggles and jiggles) of the crowd changed at different temperatures.
• The Outcome: The new method (Model I) matched the supercomputer data almost perfectly for almost every measurement. The old method (Model II) started to drift away and give wrong answers.

Why This Matters

This paper is like fixing the blueprint for a building. The old blueprint worked okay for the ground floor, but if you tried to build the 10th floor, the math would collapse.

By reformulating how we handle the "pushing" of particles and using a clever "universal shift" trick, this paper gives us a more reliable, consistent way to understand the hot, dense matter that existed at the birth of the universe and exists inside neutron stars today. It proves that sometimes, to understand the complex quantum world, you just need to find the right way to simplify the math without losing the truth.

Technical Summary

1. Problem Statement

The Hadron Resonance Gas (HRG) model is a standard effective approach for describing strongly interacting matter at high temperatures and low-to-moderate chemical potentials. While the model successfully incorporates attractive interactions by treating unstable resonances as stable particles, incorporating repulsive interactions remains challenging.

Current implementations of repulsion, primarily via the Excluded Volume (EV) method, introduce density-dependent shifts to the chemical potential ($\mu \to \mu - \nu$). The paper identifies two critical issues with traditional EV treatments:

1. Thermodynamic Inconsistency: The density-dependent chemical potential shifts are often phenomenological and not derived from the S-matrix framework. This leads to incorrect contributions to higher-order derivatives of the pressure (susceptibilities), making fluctuation observables highly sensitive to the specific implementation of repulsion.
2. Uncertainty in Hadron Sizes: The excluded volume depends on the radii of hadrons. While the pion radius is estimated around 0.2 fm, radii for other hadrons are largely unknown, hindering precise model calculations.

2. Methodology

The author proposes a reformulation of the EV treatment to ensure thermodynamic consistency and reduce model dependence.

A. Auxiliary Classical Representation

Instead of applying density-dependent shifts directly to the quantum statistical distribution functions (which complicates higher-order derivatives), the author constructs an auxiliary classical representation (Maxwell-Boltzmann statistics).

• Core Concept: A single, common energy shift $E$ is introduced for all hadron species.
• Constraint: This shift $E$ is determined by requiring the scalar number density in the auxiliary classical picture to equal the scalar number density in the actual quantum statistical picture.
• Mathematical Formulation:
  • The actual quantum number density ($n_{actual}$) is calculated using the standard Fermi-Dirac/Bose-Einstein integrals with individual species shifts $\nu_i$.
  • The auxiliary classical density ($n_{aux}$) is calculated using Maxwell-Boltzmann statistics with a common shift $E$:
    $$n_{aux} = \sum_i \frac{g_i}{2\pi^2} m_i^2 T K_2\left(\frac{m_i}{T}\right) e^{(\mu_i + E)/T}$$
  • $E$ is solved numerically by equating $n_{actual} = n_{aux}$.
• Thermodynamic Consistency: Once $E$ is determined, the pressure and all higher-order susceptibilities are calculated using total derivatives. This ensures that the implicit dependence of $E$ on the chemical potential $\mu$ is correctly accounted for, preserving thermodynamic consistency.

B. Hadron Radius Parametrization (Liquid-Drop Model)

To address the lack of known hadron radii, the paper introduces a phenomenological mass-radius relation inspired by the liquid-drop model:

• Assumption: Hadrons are treated as spherical liquid drops with constant inner pressure.
• Derivation: By integrating the pressure work required to expand a liquid drop from zero mass to a hadron of mass $m$, a power-law relation between radius ($R_h$) and mass ($m_h$) is derived:
  $$R_h = R_\pi \left( \frac{m_\pi}{m_h} \right)^A$$
  where $R_\pi$ is the pion radius and $A$ is a scaling exponent.
• Effective Radius: The mechanical radius used in the excluded volume calculation is defined as $r_i = \sqrt{3/5} R_h$.
• Parameters: The model relies on only two adjustable parameters:
  1. Pion radius ($R_\pi$).
  2. Scaling exponent ($A$).

3. Key Contributions

1. Reformulation of EV Corrections: The paper moves away from ad-hoc density-dependent chemical potential shifts. By mapping the problem to an auxiliary classical frame with a common energy shift, it resolves the thermodynamic inconsistencies associated with higher-order derivatives in standard EV models.
2. Minimal Parameterization: The framework successfully reproduces complex lattice QCD data using only two parameters ($R_\pi$ and $A$), avoiding the need for individual radii for every hadron species.
3. Unified Treatment: It provides a unified approach that handles both the statistical mechanics of repulsion and the geometric scaling of hadron sizes within a single consistent framework.

4. Results

The model (referred to as Model I) was tested against Lattice QCD (LQCD) data for conserved charge susceptibilities at zero chemical potential ($\mu_B = 0$).

• Second-Order Susceptibilities ($\chi_2$): The model reproduces LQCD data for baryon number ($\chi_B^2$), electric charge ($\chi_Q^2$), and strangeness ($\chi_S^2$) with high accuracy.
• Fourth-Order Susceptibilities ($\chi_4$):
  • Good agreement is observed for $\chi_B^4$ and $\chi_S^4$.
  • $\chi_Q^4$ shows some deviation, which the author attributes to larger lattice uncertainties in the LQCD data for this specific observable.
• Crossed Susceptibilities: Mixed susceptibilities (e.g., $\chi_{BQ}^{11}$, $\chi_{BS}^{11}$) also show good agreement with lattice data.
• Comparison with Model II: A comparison with a traditional treatment (Model II, where no equivalent effective particle picture is used) demonstrates that the proposed reformulation (Model I) yields significantly better agreement with lattice data, particularly for higher-order fluctuations.
• Parameters Used: The best fit was achieved with a pion radius $R_\pi = 0.2$ fm and a scaling exponent $A = 3$.

5. Significance and Limitations

• Significance: This work provides a robust, thermodynamically consistent method for including repulsive interactions in the HRG model. It bridges the gap between phenomenological excluded-volume models and rigorous statistical mechanics, offering a reliable tool for interpreting heavy-ion collision data and neutron star physics where fluctuation observables are crucial.
• Limitations:
  • The model assumes hadrons are spherical liquid drops, which is an oversimplification of their complex internal structure.
  • Discrepancies in higher-order strange sector susceptibilities suggest the potential existence of undiscovered high-mass strange hadrons not included in the current hadron list (QMHRG2020).
  • The model is currently optimized for zero chemical potential; its performance at high baryon densities (relevant for neutron stars) requires further investigation.

In conclusion, Somenath Pal presents a streamlined, physically motivated framework that resolves long-standing thermodynamic inconsistencies in the HRG model, enabling precise reproduction of Lattice QCD results with minimal phenomenological input.