$S$-matrix calculation of $BQ$ correlation at finite baryon density

This paper calculates the baryon number–electric charge susceptibility at finite baryon density using an S-matrix formalism for pion-nucleon interactions within a hadron gas model, revealing a significant increase with chemical potential and evaluating its behavior along the chemical freeze-out line and in a cooling fireball under Partial Chemical Equilibrium.

The Gist

Imagine the universe, just a fraction of a second after the Big Bang, was a super-hot, super-dense soup of tiny particles. Physicists call this "quark-gluon plasma." As the universe cooled down, this soup froze into a "gas" of particles we know today: protons, neutrons, pions, and their many cousins.

This paper is about trying to understand the rules of that "freezing" process, specifically when the soup is heavy with matter (baryons) rather than just empty energy. The authors are trying to build a better map of how these particles interact, which helps us understand if there is a hidden "critical point" in the universe's phase diagram—a place where matter changes state in a dramatic, explosive way.

Here is the breakdown of their work using simple analogies:

1. The Problem: A Crowded Dance Floor

Think of the hot matter created in particle colliders (like the Large Hadron Collider) as a massive, crowded dance floor.

• The Old Model (HRG): Previously, scientists treated this dance floor like a room full of people just standing around or bumping into each other randomly. They counted the dancers but ignored how they actually held hands or danced together. This is called the "Hadron Resonance Gas" model.
• The New Model (S-Matrix): The authors say, "Wait, these particles aren't just standing there; they are interacting!" They use a mathematical tool called the S-matrix (Scattering Matrix). Think of this as a detailed choreography guide. It doesn't just count the dancers; it calculates exactly how they spin, bounce, and form temporary pairs (resonances) before breaking apart.

2. The Goal: Measuring the "Baryon-Charge" Connection

The authors are calculating something called $\chi_{BQ}$ (Baryon-Number–Electric-Charge Susceptibility).

• The Analogy: Imagine you are trying to predict how the crowd moves. If you know how many people are wearing red shirts (Baryon number) and how many are wearing blue hats (Electric charge), can you predict if they will clump together?
• The "Susceptibility": This number tells you how strongly the "red shirts" and "blue hats" are linked. If the number is high, it means the presence of a proton (red shirt) strongly influences the presence of an electric charge nearby.
• Why it matters: If this number spikes or changes suddenly as the temperature drops, it might be a sign that the matter is hitting a Critical Point—a special spot in the universe's history where physics behaves strangely.

3. The Twist: Adding "Strangeness" and "Three-Body" Chaos

The authors improved their calculation in two major ways:

• The "Strangeness" Rule: In the real universe, the total amount of "strangeness" (a property of certain heavy particles) must stay neutral, like a balanced scale. The authors made sure their dance floor simulation kept this balance, even when they added more "heavy" dancers (baryons).
• The "Three-Person" Dance: Most models only look at two particles bumping into each other. But sometimes, three particles interact at once (like a $\pi\pi N$ interaction).
  • The Metaphor: Imagine two people dancing, and a third person jumps in to change the rhythm. The authors had to invent three different ways to estimate how much this "three-person dance" changes the overall crowd behavior. They found that while it's hard to calculate exactly, it adds a significant "systematic uncertainty" (about 20%) to their results.

4. The Journey: From Chemical Freeze-Out to Kinetic Freeze-Out

The paper tracks what happens as the "fireball" (the hot soup) cools down.

• Chemical Freeze-Out: This is the moment the "recipe" is set. The number of protons, pions, and other particles stops changing. It's like the baker taking the cake out of the oven; the ingredients are fixed.
• Kinetic Freeze-Out: This is later, when the particles stop bumping into each other entirely and fly off to the detectors. The cake has cooled down completely.
• The PCE Model: The authors used a "Partial Chemical Equilibrium" model to simulate the time between these two moments.
  • The Metaphor: Imagine the cake is cooling. The ingredients (protons, pions) are no longer baking into new things, but they are still jostling around. The authors found that as the temperature drops from the "baking" point to the "cooling" point, the connection between baryons and charge ($\chi_{BQ}$) actually decreases.

5. The Big Discovery

When they crunched the numbers for a universe with a lot of matter (high baryon density, like in lower-energy collisions):

1. The Link Gets Stronger: As the density of matter increases, the connection between baryon number and electric charge gets much stronger.
2. The Cooling Effect: However, as the matter cools down after the chemical freeze-out, this connection weakens significantly (dropping to about 60% of its peak value).

Why Should You Care?

If you are an experimentalist looking for the "Critical Point" of the universe, you need a baseline. You need to know what the "normal" behavior looks like so you can spot the "weird" behavior.

This paper says: "Don't just look at the moment the particles stop reacting chemically. You have to account for the cooling process afterward, because the signal you are looking for might get diluted as the soup cools down."

By using a more realistic "choreography" (S-matrix) instead of a simple "standing crowd" (HRG), they provide a more accurate map. This helps physicists distinguish between normal physics and the exotic physics of the universe's critical point.

Technical Summary

1. Problem Statement

The paper addresses the calculation of the baryon number–electric charge susceptibility ($\chi_{BQ}$) in a hot, strongly interacting hadron gas at finite baryon density (non-vanishing baryo-chemical potential, $\mu_B$).

• Context: While the thermodynamics of hadron gas at $\mu_B = 0$ has been successfully described using the Hadron Resonance Gas (HRG) model improved by the S-matrix formalism (which accounts for interactions via phase shifts), previous calculations were limited to zero net-baryon density.
• Motivation: Fluctuations and correlations of conserved charges at finite $\mu_B$ are critical observables for locating the QCD critical point and understanding the phase diagram in relativistic heavy-ion collisions (e.g., at RHIC and FAIR).
• Gap: There was a need to extend the S-matrix formalism to finite $\mu_B$ to provide a reliable baseline for experimental data, distinguishing between standard hadronic interactions and potential critical point signals.

2. Methodology

The authors employ a statistical mechanics approach where the partition function ($Z$) is split into a free part ($Z_0$) and an interaction part ($Z_{int}$).

A. Partition Function Construction

The total partition function is defined as:
$$\ln Z = \ln Z_0 + \ln Z_{int}$$

1. Free Part ($\ln Z_0$): Includes all stable ground-state hadrons (meson octet, baryon octet, and decuplet) treated as zero-width particles.
2. Interaction Part ($\ln Z_{int}$): This is the core innovation, divided into three components:
   • Resonances ($\ln Z_R$): Standard zero-width resonances (analogous to HRG) for most channels.
   • Phase Shifts ($\ln Z_i$): Explicitly treats $\pi N$ (pion-nucleon) and $\eta N$ interactions using experimental phase shifts from the GWU/SAID partial wave analysis. This avoids the double-counting of resonances (like $N^*$ and $\Delta$) that are already included in the phase shifts.
   • Three-Body Interactions ($\ln Z_3$): Accounts for $\pi\pi N$ final states (decays of $N^*$ and $\Delta$) which are not fully captured by two-body phase shifts.

B. Treatment of Three-Body Interactions

Since a rigorous S-matrix treatment for three-body interactions is complex, the authors test three phenomenological models to estimate systematic uncertainty:

1. $\Delta(1232)$ Model: Modifies the mass of the $\Delta(1232)$ resonance to fit lattice QCD data at $T=155$ MeV, subtracting the original $\Delta$ to prevent double counting.
2. Shifted $N^*$ and $\Delta$ 1: Scales masses of $N^*$ and $\Delta$ resonances by factors fitted to lattice data at $T=135$ and $155$ MeV.
3. Shifted $N^*$ and $\Delta$ 2: Similar to model 1 but fitted over a wider temperature range ($135-160$ MeV).

C. Thermodynamic Quantities

The susceptibility is calculated as the second derivative of the partition function:
$$\chi_{BQ} = \frac{T}{V} \frac{\partial^2 \ln Z}{\partial \mu_B \partial \mu_Q}$$
The authors enforce strangeness neutrality ($n_S = 0$) to determine the strangeness chemical potential $\mu_S$ for a given $T$ and $\mu_B$.

D. Evolution Scenarios

The study evaluates $\chi_{BQ}$ under two distinct conditions:

1. Chemical Freeze-out: Using a parametrized freeze-out curve $T(\mu_B)$ derived from hadron yield fits.
2. Partial Chemical Equilibrium (PCE): Simulating the cooling of the fireball after chemical freeze-out. In this model, the effective numbers of stable hadrons (including decay products) are conserved as temperature drops, while resonances remain in equilibrium with their daughters.

3. Key Contributions

• Extension to Finite Density: Successfully extended the S-matrix formalism for $\chi_{BQ}$ from $\mu_B=0$ to the phenomenologically relevant region of finite baryon density.
• Refined Interaction Modeling: Explicitly separated the treatment of $\pi N$ interactions (via phase shifts) from other resonances, resolving inconsistencies in previous HRG models regarding resonance counting.
• Systematic Uncertainty Estimation: Provided a quantitative assessment of the uncertainty arising from the treatment of three-body ($\pi\pi N$) interactions by comparing three distinct models.
• PCE Evolution: Calculated the evolution of $\chi_{BQ}$ during the cooling phase of the fireball, moving beyond the static chemical freeze-out assumption.

4. Key Results

• Dependence on $\mu_B$:

  • At $\mu_B = 0$, the S-matrix results align well with recent Lattice QCD data.
  • As $\mu_B$ increases, $\chi_{BQ}$ grows significantly. The inclusion of interactions (specifically phase shifts) changes the temperature dependence from decreasing (in non-interacting models) to increasing.
  • The $I=3/2$ channel (dominated by $\Delta$ resonances) and three-body interactions contribute substantially to the enhancement of the susceptibility.

• Strangeness Neutrality:

  • Enforcing strangeness neutrality ($n_S=0$) enhances the $B-Q$ correlation compared to non-neutral matter. The dependence on temperature becomes smoother at high $\mu_B$.

• Comparison with HRG:

  • The S-matrix treatment yields slightly lower values of $\chi_{BQ}$ compared to the standard HRG model (which treats all resonances as stable particles without phase shifts), but the difference is not dramatic relative to the overall magnitude.

• PCE Cooling Effects:

  • As the fireball cools from chemical freeze-out ($T_{fo} \approx 150-160$ MeV) down to kinetic freeze-out ($T \approx 100$ MeV) under PCE conditions, the scaled susceptibility $\chi_{BQ}/T^2$ increases.
  • However, the absolute susceptibility $\chi_{BQ}$ decreases significantly (dropping to ~60% of its chemical freeze-out value at $T=100$ MeV for $\sqrt{s_{NN}} = 7.7$ GeV). This is a crucial finding for interpreting experimental data, as measurements occur at kinetic freeze-out.

• Systematic Uncertainty:

  • The variation between the three models for $\pi\pi N$ interactions results in a systematic uncertainty of approximately 20% in the final $\chi_{BQ}$ calculation.

5. Significance

• Baseline for Critical Point Search: The calculated $\chi_{BQ}$ serves as a robust "hadronic baseline." Any experimental deviation from these values in heavy-ion collisions could signal the presence of a critical point or a first-order phase transition.
• Interpretation of Experimental Data: The paper highlights that interpreting measured fluctuations requires accounting for the cooling evolution (PCE). The susceptibility at the time of particle emission (kinetic freeze-out) is significantly different from the value at chemical freeze-out.
• Model Validation: The work validates the S-matrix approach as a superior tool over the standard HRG model for describing interacting hadron matter, particularly in the context of finite baryon density where resonance structures play a dominant role.

In summary, this paper provides a comprehensive theoretical framework for understanding baryon-charge correlations in dense nuclear matter, bridging the gap between lattice QCD calculations at zero density and experimental heavy-ion collision data at finite density.