Kinetic Freeze-Out Conditions and Net Baryon Density in Au+Au Collisions at $\sqrt{s_{NN}} = 7.7$--$39$ GeV within a Collective Flow Fireball Model

This study utilizes a covariant statistical fireball model to analyze Au+Au collisions at RHIC energies, revealing that longitudinal flow induces a kinematic degeneracy that artificially elevates kinetic freeze-out temperatures above the QCD crossover limit, thereby suggesting large longitudinal velocities are physically disfavored, while simultaneously identifying a net baryon density maximum below 11.5 GeV enhanced by dynamical flow.

The Gist

Imagine you are a detective trying to figure out what happened inside a massive, high-speed car crash. But instead of cars, these are two gold atoms smashing into each other at nearly the speed of light. When they hit, they create a tiny, super-hot, super-dense ball of energy called a "fireball."

This fireball is made of the fundamental building blocks of matter (quarks and gluons) that usually live inside protons and neutrons. As this fireball expands and cools down, it eventually "freezes." At this moment, called Kinetic Freeze-Out, the particles stop bumping into each other and fly off in straight lines to be caught by detectors.

The goal of this paper is to figure out exactly when and where this freezing happens, and what the conditions (temperature and density) were at that exact moment.

Here is the breakdown of their investigation using simple analogies:

1. The "Static" vs. "Moving" Fireball

In previous studies, scientists treated this fireball like a stationary hot potato. They assumed it just sat there cooling down. They looked at the particles flying out and guessed the temperature based on how fast they were moving.

However, the authors of this paper realized the fireball isn't just sitting still; it's exploding outward in all directions.

• The Analogy: Imagine a balloon being blown up. If you put a speck of dust on the balloon, the dust moves away from the center not just because it's vibrating (heat), but because the balloon itself is stretching (flow).
• The Problem: If you ignore the stretching (the flow), you might think the dust is vibrating faster than it actually is. You would guess the temperature is higher than it really is.

2. The Two Types of Flow

The team looked at two ways the fireball expands:

• Transverse Flow (Sideways): Like the balloon expanding outward in all directions.
• Longitudinal Flow (Forward/Backward): Like the balloon stretching long and thin along the direction the cars were driving.

The researchers asked: "What if the fireball is stretching forward really fast? Does that change our temperature reading?"

3. The "Speed Trap" (The Main Discovery)

They ran their math models with three different assumptions for how fast the fireball was stretching forward:

1. Not stretching at all ($v_z = 0$).
2. Stretching a little ($v_z = 0.2$).
3. Stretching a lot ($v_z = 0.4$).

The Surprise:
When they assumed the fireball was stretching forward fast ($v_z = 0.4$), their math forced the temperature to skyrocket to values that don't make sense physically (over 200 MeV).

• The Metaphor: It's like trying to measure the speed of a runner. If you forget that the runner is on a moving treadmill, you might think they are running at 20 mph when they are actually only running at 10 mph, and the treadmill is doing the rest.
• The Result: The math showed a "degeneracy." The model couldn't tell the difference between "hot particles" and "particles being pushed by a fast-moving fireball." To make the math fit the data, the model had to invent a super-high temperature to compensate for the fast forward motion.

The Conclusion: Since a temperature of 200+ MeV is too hot for the "Hadron Gas" (the soup of particles) to exist—it would have already turned into a different state of matter (Quark-Gluon Plasma)—the authors concluded that the fireball probably isn't stretching forward that fast. The "fast stretch" scenario is likely unphysical.

4. The "Baryon Density" Peak

The paper also looked at how "crowded" the fireball is with protons and neutrons (baryons).

• The Analogy: Think of a crowd of people. At low-energy crashes, the people are packed tight in the middle. At high-energy crashes, they zip past each other so fast that the middle is almost empty.
• The Finding: The researchers found that the "crowdedness" (net baryon density) hits a maximum at a specific energy level (around 11.5 GeV).
• Why it matters: This is the "sweet spot" where the matter is most compressed. If scientists want to find a "Critical Point" (a special phase transition in the universe's history), this is the best place to look.

5. Why This Matters

This paper is like refining a map for future explorers.

• For the "Thermometer": It tells us that if we ignore the forward motion of the fireball, we get the wrong temperature. We need to be careful not to confuse "flow" with "heat."
• For the "Critical Point": It confirms that the best place to look for the mysterious "Critical Point" of the universe is at intermediate energies, where the matter is most compressed.
• For Future Experiments: It gives better guidelines for experiments at facilities like the RHIC (Relativistic Heavy Ion Collider) and FAIR in Germany, helping them tune their machines to catch the most interesting physics.

Summary in One Sentence

The authors discovered that if you assume the exploding fireball is stretching forward too fast, your math tricks you into thinking it's impossibly hot; by correcting for this, they found the true temperature and confirmed that the "crowdedness" of the collision peaks at a specific energy, guiding scientists to the best place to hunt for the secrets of the early universe.

Technical Summary

Here is a detailed technical summary of the paper "Kinetic Freeze-Out Conditions and Net Baryon Density in Au+Au Collisions at $\sqrt{s_{NN}} = 7.7–39$ GeV within a Collective Flow Fireball Model."

1. Problem Statement

The primary objective of heavy-ion collision experiments, such as the RHIC Beam Energy Scan (BES), is to map the Quantum Chromodynamics (QCD) phase diagram and locate the critical point. A crucial step in this process is determining the kinetic freeze-out conditions—the state of the system when elastic interactions cease and particle momentum distributions are fixed.

Previous studies (e.g., Randrup and Cleymans) identified a maximum in net baryon density ($\rho_B$) at intermediate collision energies using a static fireball model. However, this approach neglected collective flow (both transverse and longitudinal), which is a dynamical feature of the expanding fireball. Neglecting flow biases the extraction of thermodynamic parameters (temperature $T$ and chemical potential $\mu_B$). This paper addresses the gap by investigating how the inclusion of longitudinal collective flow ($v_z$) modifies the kinetic freeze-out trajectory and the inferred net baryon density in the $(\rho_B, \epsilon)$ plane.

2. Methodology

The authors employ a covariant statistical fireball model to analyze experimental data from the STAR Collaboration (0–5% central Au+Au collisions at $\sqrt{s_{NN}} = 7.7–39$ GeV).

• Theoretical Framework:

  • The model generalizes the static thermal source by introducing a dynamical expansion with four-velocity $u^\mu$.
  • The invariant particle spectrum is derived from the Lorentz-invariant distribution function $f \propto \exp(-p_\mu u^\mu / T)$, where $p_\mu$ is the particle four-momentum.
  • The distribution function incorporates both transverse flow velocity ($v_T$) and longitudinal flow velocity ($v_z$).
  • The resulting differential spectrum (Eq. 17) involves modified Bessel functions ($I_0, I_1$) arising from the azimuthal integration.

• Fitting Procedure:

  • The model fits the transverse momentum ($p_T$) spectra of protons and positive pions measured by STAR.
  • Free Parameters: Temperature ($T$), baryon chemical potential ($\mu_B$), transverse flow velocity ($v_T$), and normalization constant ($\lambda$).
  • Fixed Parameter: Longitudinal flow velocity ($v_z$) is fixed at three distinct values: 0.0, 0.2, and 0.4 to isolate its systematic effects.
  • Constraints: Fits are performed in the mid-rapidity window $|y| < 0.1$.

• Thermodynamic Reconstruction:

  • Once $T$ and $\mu_B$ are extracted, the authors calculate the energy density ($\epsilon$) and net baryon density ($\rho_B$) using a Hadron Resonance Gas (HRG) model with full quantum statistics.
  • Parametrizations for $\mu_B(\sqrt{s_{NN}})$ and $T(\mu_B)$ are re-fitted to ensure thermodynamic consistency with the flow-modified data.

3. Key Contributions

• Systematic Inclusion of Longitudinal Flow: Unlike previous static models, this work explicitly treats longitudinal flow as a variable to quantify its impact on extracted freeze-out parameters.
• Identification of Kinematic Degeneracy: The paper demonstrates a fundamental kinematic degeneracy between longitudinal velocity ($v_z$) and temperature ($T$) in the Lorentz-invariant distribution function.
• Physical Consistency Check: The study uses the extracted temperatures to test the validity of the HRG model against the QCD crossover temperature ($T_c \approx 155–160$ MeV), effectively ruling out unphysical flow scenarios.
• Refined $\rho_B$ Mapping: It reconstructs the freeze-out trajectory in the $(\rho_B, \epsilon)$ plane, quantifying how dynamical flow enhances the inferred baryon compression.

4. Key Results

A. Temperature and Chemical Potential

• Longitudinal Flow Effect on $T$: There is a systematic upward shift in the extracted freeze-out temperature as $v_z$ increases.
  • $v_z = 0.0$: $T \in [143, 171]$ MeV.
  • $v_z = 0.2$: $T \in [150, 180]$ MeV.
  • $v_z = 0.4$: $T \in [175, 209]$ MeV.
  • Mechanism: This shift is not due to spectral hardening (since longitudinal motion does not transfer momentum to the transverse plane). Instead, it arises because a non-zero $v_z$ increases the effective energy in the fluid rest frame ($p_\mu u^\mu$). To reproduce the observed $p_T$ spectra, the fit compensates by increasing $T$.
• Chemical Potential ($\mu_B$): $\mu_B$ is largely insensitive to $v_z$. It decreases monotonically from $\sim 420$ MeV at 7.7 GeV to $\sim 200$ MeV at 39 GeV, governed by initial baryon stopping dynamics rather than collective expansion.

B. Physical Validity and the QCD Crossover

• The $v_z = 0.4$ Problem: At $v_z = 0.4$, the extracted temperatures ($175–209$MeV) significantly exceed the lattice QCD crossover temperature ($T_c \approx 155–160$ MeV).
• Implication: Since the model assumes a Hadron Resonance Gas (HRG), temperatures above $T_c$ imply the system is in a deconfined Quark-Gluon Plasma (QGP) phase, rendering the HRG model invalid. This suggests that large longitudinal velocities ($v_z \approx 0.4$) are physically disfavored at these collision energies.
• Reliable Regime: The results for $v_z = 0.0$ and $v_z = 0.2$ yield temperatures consistent with $T_c$, making them the most physically reliable extractions.

C. Net Baryon Density ($\rho_B$) and Energy Density ($\epsilon$)

• Maximum Density: A maximum in net baryon density is confirmed at intermediate energies ($\sqrt{s_{NN}} \lesssim 11.5$ GeV).
• Flow Enhancement: Dynamical flow enhances the inferred baryon compression. Compared to the static limit ($v_z=0$), the dynamical scenario (specifically $v_z=0.4$) increases the inferred $\rho_B$ by up to ~20%.
• Yield Independence: The integrated proton yields are found to be essentially independent of $v_z$. This is because the net-baryon content at mid-rapidity is established during the pre-equilibrium stopping phase, and the narrow rapidity window ($|y|<0.1$) is insensitive to moderate rigid rapidity boosts.

5. Significance and Outlook

• Benchmark for Hydrodynamics: The extracted freeze-out coordinates ($T, \mu_B, \rho_B, \epsilon$) provide refined benchmarks for future 3+1 dimensional viscous hydrodynamic simulations, particularly for experiments at FAIR and low-energy RHIC runs.
• Critical Point Search: By constraining the freeze-out trajectory relative to the QCD crossover and critical point, this work helps define the "safe" parameter space for searching for the critical point.
• Future Directions: The authors suggest that future work should:
  1. Include strange hadrons ($K, \Lambda, \Xi$) to self-consistently determine $\mu_S$ and $\mu_Q$.
  2. Impose $T \leq T_c$ as a hard constraint in fitting to eliminate unphysical solutions.
  3. Utilize wider rapidity windows to break the kinematic degeneracy between $v_z$ and $T$.
  4. Couple the statistical model with full hydrodynamic evolution to dynamically determine the freeze-out surface.

In conclusion, this study demonstrates that neglecting longitudinal flow or assuming unphysically high flow velocities leads to erroneous thermodynamic interpretations. It establishes that while flow significantly impacts the inferred temperature and density, the physical constraints of the QCD phase transition ($T \leq T_c$) serve as a critical filter for validating flow models in the BES energy range.