Thermal Radiation from an Analytic Hydrodynamic Model with Hadronic and QGP Sources in Heavy-Ion Collisions

This paper presents a fully analytic hydrodynamic model incorporating a lattice-QCD-consistent equation of state to describe thermal photon production across the quark-hadron transition, demonstrating good agreement with PHENIX experimental data for Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV and enabling the investigation of initial temperature centrality dependence.

The Gist

Imagine the universe just a fraction of a second after the Big Bang. It wasn't filled with stars, planets, or even atoms. Instead, it was a super-hot, super-dense soup of the most fundamental building blocks of matter: quarks and gluons. Physicists call this state Quark-Gluon Plasma (QGP).

Today, scientists recreate this primordial soup in massive particle accelerators by smashing heavy gold atoms together at nearly the speed of light. The goal? To understand how the universe began and how matter behaves under extreme conditions.

This paper is about a new way to "listen" to that collision and figure out how hot it got. Here is the story, broken down into simple concepts.

1. The Problem: The "Flashlight" vs. The "Fog"

When these atoms collide, they create a tiny, expanding fireball.

• The Hadrons (The Fog): As the fireball cools down, the quarks and gluons freeze together to form normal particles (like protons and pions). These are like the fog that forms after a hot shower. They tell us what the fireball looked like at the very end of its life, but they are too slow to tell us what happened at the very beginning.
• The Photons (The Flashlight): Real light (photons) is different. Once created, they don't interact with the "fog." They zip straight out of the collision zone without getting stuck. They carry a perfect snapshot of the temperature and conditions at the exact moment they were born.

The challenge is that the collision produces light from two different "eras":

1. The Early Era (QGP): The super-hot, early stage where quarks are free.
2. The Late Era (Hadronic): The cooler, later stage where particles have already formed.

Scientists want to know: How hot was the fireball at the very start? To find out, they need to separate the "early light" from the "late light."

2. The Old Way vs. The New Way

Previously, scientists tried to model this using simple math, often assuming the fireball was just one uniform thing. It was like trying to describe a complex orchestra by only listening to the drums. It gave a rough idea, but the details were fuzzy, and the estimated starting temperatures were often wildly uncertain.

The New Model (The "Two-Stage" Recipe):
The author, Gábor László Kasza, built a new mathematical model that treats the fireball as having two distinct phases:

1. The QGP Phase: The hot, early soup.
2. The Hadronic Phase: The cooling, particle-forming soup.

He didn't just guess how they mix; he used advanced math (relativistic hydrodynamics) to describe exactly how the fireball expands and cools, switching from one phase to the other. Crucially, he made the math analytic, meaning he solved the equations with exact formulas rather than relying on slow, brute-force computer simulations.

3. The "1D" Shortcut

You might wonder: "Collisions happen in 3D space. Why does this model only look at 1 dimension (length)?"

The Analogy: Imagine you are trying to figure out how hot a loaf of bread is by looking at the crust. You don't need to know the exact shape of every crumb inside to know the oven temperature.

• The "heat" (temperature) of the fireball is the most important factor for the light it emits.
• While the sideways expansion (transverse flow) does affect the light, the temperature history is the dominant signal.
• By simplifying the model to just the "length" of the expansion (ignoring the sideways flow for a moment), the author could keep the math solvable and elegant. It's a "skeleton" model that captures the essential bones of the physics without getting bogged down in the muscle.

4. The Experiment: Checking the Recipe

The author took his new "Two-Stage" recipe and compared it to real data from the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). They looked at gold-gold collisions at 200 GeV.

The Results:

• The Fit: The model matched the experimental data very well. It successfully separated the "early light" (QGP) from the "late light" (hadrons).
• The Temperature: When the model included the "late light" (hadrons), the estimated starting temperature of the fireball was higher (around 430–490 MeV) and more stable across different collision sizes.
• The Surprise: If you ignore the late light (as older models sometimes did), the math suggests the starting temperature drops and changes wildly depending on how big the collision was. This proves that you cannot ignore the "late light" if you want to know the true starting temperature. The "late light" acts like a background noise that, if not accounted for, distorts the picture of the "early light."

5. Why This Matters

• A New Benchmark: This paper provides a clean, mathematical "yardstick" for future studies. Before, scientists had to rely on massive, complex computer simulations. Now, they have a precise formula to check those simulations against.
• Understanding the Transition: It helps us understand the exact moment matter switches from a "free quark soup" to "frozen particles."
• Simplicity with Power: It shows that you don't always need a supercomputer to get deep insights. Sometimes, a clever, simplified mathematical model can reveal the core truth of a complex physical phenomenon.

The Takeaway

Think of this paper as a detective story. The universe left behind a clue (the light from the collision). Previous detectives tried to read the clue with a blurry lens. This author built a new, sharper lens that separates the "early clues" from the "late noise." By doing so, he gave us a much clearer picture of how hot the universe was in its very first moments.

The model says: "To know how hot the fire was at the start, you have to account for the smoke that came out at the end."

Technical Summary

Here is a detailed technical summary of the paper "Thermal Radiation from an Analytic Hydrodynamic Model with Hadronic and QGP Sources in Heavy-Ion Collisions" by Gábor László Kasza.

1. Problem Statement

High-energy heavy-ion collisions create a strongly coupled quark-gluon plasma (QGP) that evolves into hadronic matter. While hadronic observables reflect conditions at kinetic freeze-out, direct photons serve as unique probes of the entire evolution history because they escape the medium without strong interaction.

• The Challenge: Experimental data (specifically from PHENIX) shows that the transverse momentum ($p_T$) spectrum of non-prompt direct photons cannot be adequately described by a simple Boltzmann distribution with a constant inverse slope. This implies that the spectrum is a superposition of sources with different temperatures and flow profiles (QGP vs. hadronic).
• The Gap: Previous analytic hydrodynamic models often assumed a single-component medium or a constant equation of state (EoS) parameter ($\kappa = \epsilon/p$). These models struggled to extract precise initial temperatures ($T_0$) and often overestimated them due to the lack of a realistic treatment of the quark-hadron transition and the temperature dependence of the EoS. Furthermore, numerical models that handle these complexities well often lack the transparency of analytic solutions.

2. Methodology

The author develops a completely analytic model based on relativistic perfect fluid hydrodynamics in 1+1 dimensions (longitudinal expansion only).

A. New Hydrodynamic Solutions

• Generalized Equation of State: Unlike standard analytic solutions that assume a constant $\kappa$, this model introduces a temperature-dependent $\kappa(T)$. The author derives a condition (Eq. 13) allowing for closed-form analytic solutions where $\kappa(T)$ varies.
• Two-Phase Structure: The model explicitly separates the evolution into two phases:
  1. QGP Phase ($T > T_{tr}$): Modeled with an EoS approaching the conformal limit ($\epsilon = 3p$) at high temperatures.
  2. Hadronic Phase ($T < T_{tr}$): Modeled with a distinct EoS behavior.
• Matching: The two phases are joined at a transition temperature $T_{tr}$ (fixed at 156 MeV, consistent with lattice QCD crossover). The velocity field is assumed to be locally accelerating, described by a fluid rapidity $\Omega(\eta_z)$.

B. Thermal Photon Emission Calculation

• Source Function: The model uses a collective source function integrating over the 4D space-time volume, rather than a freeze-out hypersurface (Cooper-Frye), acknowledging that photons escape continuously.
• Two-Component Spectrum: The total photon yield is calculated as the sum of contributions from the QGP phase and the hadronic phase.
• Approximations:
  • 1+1D Framework: Transverse flow is neglected. The model assumes transverse homogeneity, embedding the 1+1D solution into a 1+3D spacetime for the spectrum calculation.
  • Saddle-Point Approximation: Used to evaluate the integrals analytically, resulting in closed-form expressions for the $p_T$ spectrum involving incomplete Gamma functions.
  • Boltzmann Weight: A simplified Boltzmann distribution is used for the source, neglecting complex microscopic rate prefactors (like logarithmic terms) to maintain analytic tractability.

C. Data Comparison and Validation

• Datasets: The model is fitted to PHENIX data for non-prompt direct photons in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV across four centrality classes (0–20%, 20–40%, 40–60%, 60–93%).
• Hadronic Validation: The same parameter set is tested against PHOBOS pseudorapidity density ($dN/d\eta$) data to ensure consistency with hadronic observables.
• Equation of State Reconstruction: Using the fitted parameters, the author reconstructs the $\kappa(T)$ function and compares it to Lattice QCD simulations.

3. Key Contributions

1. Analytic Two-Component Model: The first fully analytic hydrodynamic model that explicitly decomposes thermal photon production into distinct QGP and hadronic components, accounting for the quark-hadron transition.
2. Temperature-Dependent EoS in Analytic Form: Derivation of exact hydrodynamic solutions where the ratio of energy density to pressure, $\kappa(T)$, is a function of temperature, bridging the gap between simple analytic models and realistic Lattice QCD EoS.
3. Resolution of Initial Temperature Ambiguity: Demonstrates that previous overestimations of $T_0$ in single-component models were partly due to neglecting the hadronic contribution, which dominates the low-$p_T$ region.
4. Consistency Check: Validates that a single parameter set can simultaneously describe thermal photon spectra and hadronic pseudorapidity distributions within the 1+1D framework.

4. Results

• Fit Quality: The two-component model provides a good description of the non-prompt photon spectra (Confidence Level > 0.1% for all centrality classes). The fit is generally better for central collisions (0–20%, 20–40%) when both components are included.
• Initial Temperature ($T_0$):
  • With Hadronic Component: The extracted $T_0$ values are higher (e.g., ~489 MeV for 0–20% centrality) and show no strong monotonic trend with centrality, remaining roughly constant above 20% centrality within uncertainties.
  • Without Hadronic Component: When the hadronic contribution is neglected, $T_0$ drops significantly (e.g., ~432 MeV) and exhibits a clearer decreasing trend with centrality, which contradicts some theoretical expectations of a constant initial temperature.
• Hadronic Channel: The model successfully reproduces PHOBOS $dN/d\eta$ data for central and mid-central collisions using parameters consistent with the photon fits.
• Equation of State: The reconstructed $\kappa(T)$ function matches Lattice QCD data well at low temperatures ($T \lesssim T_{tr}$). At high temperatures, the model reproduces the qualitative behavior and the conformal limit ($\kappa \to 3$), though quantitative agreement for $p/T^4$ degrades slightly above 200 MeV.
• Centrality Dependence: The study suggests that the initial temperature is likely independent of centrality, but this conclusion is only robust when the hadronic contribution to the photon spectrum is properly accounted for.

5. Significance and Limitations

Significance:

• Benchmarking: Provides a transparent, analytic benchmark for future numerical 3D hydrodynamic simulations.
• Physical Insight: Clarifies the role of the hadronic phase in the low-$p_T$ photon spectrum, resolving discrepancies in previous single-component analyses.
• Theoretical Bridge: Offers a tractable analytic framework to study the interplay between the QGP and hadronic phases without the computational cost of full numerical transport models.

Limitations:

• Dimensionality: The model is strictly 1+1 dimensional. It neglects radial flow and elliptic flow ($v_2$). Consequently, it cannot address the "direct photon puzzle" (the unexpectedly large $v_2$ of direct photons).
• Viscosity: The model assumes a perfect fluid. While the good fit suggests viscous effects may not drastically alter the shape of the spectrum, they are not explicitly calculated.
• Quantitative Precision: Due to simplifying assumptions (e.g., neglecting logarithmic prefactors, transverse homogeneity), the extracted parameters (like $T_0$) should be viewed as qualitatively indicative rather than quantitatively precise.

Future Outlook:
The author suggests extending the solution to 1+3 dimensions to include transverse dynamics and elliptic flow, and developing viscous generalizations of the analytic solutions to address the direct photon puzzle more rigorously.