⚛️ nuclear theory

An expanding spherical fireball model for light hadron production at RHIC ( $\sqrt{s_{\rm NN}}=7.7$ --$39$ GeV)

This paper presents an expanding spherical fireball model that successfully describes the transverse momentum spectra and rapidity distributions of light hadrons produced in Au+Au collisions at RHIC energies between 7.7 and 39 GeV across various centralities.

Original authors: Ashutosh Dwibedi, Anupam Panja, Sabyasachi Ghosh

Published 2026-07-07

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Ashutosh Dwibedi, Anupam Panja, Sabyasachi Ghosh

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a high-energy physics experiment not as a complex machine, but as a giant, high-speed collision of two heavy trains (gold nuclei) smashing into each other. When they hit, they create a tiny, super-hot, super-dense "fireball" of matter. This paper is about understanding how that fireball expands and cools down, eventually turning into a spray of ordinary particles like pions, kaons, and protons that we can detect.

Here is a simple breakdown of what the authors did and found, using everyday analogies:

1. The Model: A Balloon Instead of a Cylinder

Most scientists usually imagine this fireball expanding like a cylinder (like a soda can being crushed and stretching out). However, the authors of this paper decided to try a different shape: a perfect sphere, like an inflating balloon.

The Analogy: Think of the fireball as a balloon being blown up. As it grows, the surface moves outward. The authors tracked how fast the radius of this "balloon" grew over time.
The Flow: Inside the balloon, the air (or in this case, the hot matter) moves faster near the surface and slower near the center, similar to how water flows faster at the edges of a widening river. They used a mathematical recipe (called a "blast-wave" profile) to describe this movement.

2. The Experiment: The "Freeze-Out" Moment

The fireball is incredibly hot at first. As it expands, it cools down. Eventually, it gets cold enough that the particles stop interacting with each other and fly off in straight lines to be detected. This moment is called "freeze-out."

The Analogy: Imagine a crowd of people dancing wildly in a hot room. As the room cools, they stop dancing and start walking out the door in specific directions. The "freeze-out" is the exact moment they stop dancing and start walking.
The Goal: The authors wanted to figure out the temperature and speed of the "dancers" (particles) at the exact moment they stopped interacting.

3. The Method: Fitting the Puzzle Pieces

The researchers had real data from the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). They looked at how many particles came out at different speeds (momentum) and at different angles (rapidity).

The Strategy: They treated the pion (a type of light particle) data as the "master key." They adjusted their balloon model until it perfectly matched the pion data for different collision "centrality" (how head-on the crash was).
- Head-on collisions (Central): Like two trains hitting dead center. The fireball is huge and lasts longer.
- Glancing blows (Peripheral): Like the trains just clipping each other. The fireball is smaller and cools down faster.
The Rule: Once they found the perfect settings for the pions, they kept those settings fixed for the other particles (protons and kaons). The only thing they were allowed to change was a specific "chemical potential" (a number that adjusts for how many of each particle type are produced).

4. The Results: What They Found

The model worked surprisingly well with very few adjustable knobs (only five main parameters).

Temperature vs. Size: They found that in big, head-on collisions, the fireball stays hot for a longer time but cools down to a lower final temperature before freezing. In small, glancing collisions, it cools down faster but ends up slightly hotter.
Speed: The "balloon" expands much faster in the center of the collision because the pressure pushing it outward is stronger there.
The Shape of the Spray: When they predicted how the particles spread out forward and backward (rapidity distribution), the model predicted a smooth, bell-curve shape (like a Gaussian distribution). This matched the real data, suggesting that even though the collision is chaotic, the overall expansion is very orderly.

5. Why This Matters (According to the Paper)

The authors argue that you don't need a super-complex, 8-parameter computer simulation to understand these collisions. A simple, spherical "balloon" model with just a few physical rules can explain:

How fast the particles are moving.
How many of each type of particle are made.
How they are distributed in space.

It's like showing that you can predict the pattern of water splashing out of a bursting pipe just by knowing the pipe's size and water pressure, without needing to simulate every single water molecule.

In Summary:
The paper says, "We built a simple spherical balloon model to describe the explosion of matter in heavy-ion collisions. By tuning it to match the data for pions, we successfully predicted the behavior of protons and kaons across many different collision energies. The model shows that the fireball expands like a balloon, cools down systematically, and produces a smooth, predictable spray of particles."

Technical Summary: An Expanding Spherical Fireball Model for Light Hadron Production at RHIC

Problem Statement
Relativistic hydrodynamics has been highly successful in describing the space-time evolution of quark-gluon plasma (QGP) produced in heavy-ion collisions (HICs) at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). However, extracting specific freeze-out properties from full hydrodynamic simulations is computationally intensive and requires precise knowledge of initial conditions and transport coefficients. Consequently, "blast-wave" (BW) parameterizations are frequently employed to describe particle spectra and collective flow. While standard BW models often assume cylindrical geometries, non-central collisions at RHIC exhibit initial spatial anisotropy that evolves into momentum anisotropy (elliptic flow). The authors note that while azimuthally integrated transverse momentum ( $p_T$ ) spectra and rapidity distributions are largely insensitive to these azimuthal anisotropies, existing models that attempt to describe both spectra and anisotropic flow simultaneously often require a large number of free parameters (e.g., eight parameters in elliptic fire-cylinder models). The challenge addressed in this work is to determine whether a simpler, spherically symmetric fireball model, employing fewer free parameters, can adequately describe the $p_T$ spectra and rapidity distributions of identified light hadrons across a broad range of collision energies ( $\sqrt{s_{NN}} = 7.7–39$ GeV) and centralities.

Methodology
The authors propose an expanding spherical fireball model to describe the medium produced in Au+Au collisions. The methodology involves the following key components:

Fireball Geometry and Evolution: The medium is modeled as a spherically symmetric fireball with a time-dependent radius $r_B(t)$ . The radius evolution is parameterized as $r_B(t) = r_0 + v_\infty [t - (1 - e^{-At})/A]$ , where $r_0$ is the initial effective radius derived from the initial elliptic overlap volume of the colliding nuclei, $v_\infty$ is the asymptotic expansion velocity, and $A$ characterizes the rate of radial flow development.
Flow Profile: Unlike the standard cylindrical blast-wave, the fluid velocity is defined radially. The radial fluid rapidity $\eta_r$ is determined such that the velocity at the surface matches the expansion rate of the radius, while the interior follows a blast-wave-like profile ( $v_r = (r/r_B(t))\dot{r}_B(t)$ ).
Freeze-out Prescription: Particle spectra are calculated using the Cooper–Frye freeze-out prescription. The model assumes a local equilibrium distribution function (Bose-Einstein or Fermi-Dirac) integrated over a kinetic freeze-out hypersurface defined at a constant time $t_f$ in the lab frame.
Parameter Extraction Strategy:
- The model parameters—kinetic freeze-out temperature ( $T_{kin}$ ), asymptotic velocity ( $v_\infty$ ), flow development parameter ( $A$ ), and freeze-out time ( $t_f$ )—are first fitted to the mid-rapidity $p_T$ spectra of pions ( $\pi^\pm$ ) for various centralities.
- These collective flow parameters are then held fixed for other hadron species (protons, antiprotons, kaons).
- The only additional free parameter for other species is the particle-dependent kinetic freeze-out chemical potential ( $\mu_{kin}$ ), which is adjusted to reproduce experimental data.
- Two fitting schemes are tested: Set-1 uses the full available $p_T$ range (0.25–2.0 GeV), while Set-2 restricts the fit to intervals where resonance feed-down and hard processes are minimized (similar to standard STAR BW analyses).

Key Contributions and Results
The study confronts the model with STAR collaboration data for $\pi^\pm$ , $K^\pm$ , $p$ , and $\bar{p}$ across six centrality bins and five collision energies.

Spectral Description: The model provides a good description of the experimental $p_T$ $p_{T}$ spectra for all identified light hadrons. The fits yield systematic trends in the extracted parameters:
- Centrality Dependence: As collisions become more central, the freeze-out time ( $t_f$ ) increases, reflecting the longer lifetime of the larger fireball. Conversely, the kinetic freeze-out temperature ( $T_{kin}$ ) decreases, indicating more substantial cooling before decoupling.
- Flow Development: The radial flow development parameter ( $A$ ) is larger in central collisions, suggesting faster flow development due to stronger pressure gradients.
- Chemical Potentials: The extracted effective kinetic chemical potentials ( $\mu_{kin}$ ) vary by particle species and centrality, effectively absorbing contributions from resonance feed-down not explicitly included in the model.
Rapidity Distributions: Using the parameters extracted from the $p_T$ spectra, the model predicts rapidity distributions ( $dN/dy_p$ ) for all hadron species. The results consistently show Gaussian-like distributions across the considered energy range and centralities, aligning with observations from lower-energy experiments (CERN SPS, GSI) and extending up to $\sqrt{s_{NN}} = 39$ GeV.
Parameter Efficiency: The model achieves these results with only five free parameters (two for radius evolution, one for freeze-out time, and temperature/chemical potential), which is fewer than more elaborate elliptic blast-wave parametrizations that require eight parameters to describe spectra and anisotropic flow simultaneously.

Significance and Claims
The authors claim that their work demonstrates the viability of a simple, spherically symmetric dynamical fireball model in describing complex heavy-ion collision observables. Specifically, they assert that:

A model with a small number of physically motivated parameters can simultaneously account for transverse momentum spectra and predict corresponding rapidity distributions in a consistent framework.
The assumption of spherical symmetry is sufficient for describing azimuthally integrated spectra and rapidity distributions, despite the initial geometric anisotropy of non-central collisions.
The extracted freeze-out parameters exhibit physically intuitive systematic dependencies on collision centrality and energy, validating the model's ability to capture the bulk evolution and cooling of the fireball.
The approach offers a computationally efficient alternative to full hydrodynamic simulations for extracting freeze-out properties, particularly for the Beam Energy Scan energies where full hydrodynamic descriptions can be challenging.

The paper concludes that the expanding spherical fireball model successfully reproduces the STAR data for light hadrons, supporting the hypothesis that the bulk properties of the medium in this energy range can be effectively captured by a simplified, spherically symmetric expansion scenario.