⚛️ phenomenology

Thermal photon production rate from QGP at finite density: Analytic cutoff independence

This paper analytically demonstrates for the first time that the thermal photon production rate from a hot and dense quark-gluon plasma is independent of the intermediate momentum cutoff even in the presence of a nonzero quark chemical potential, thereby resolving a long-standing issue and providing a generalized expression suitable for relativistic hydrodynamic simulations.

Original authors: Sourav Duari, Nilanjan Chaudhuri, Sourav Sarkar, Pradip Roy

Published 2026-10-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sourav Duari, Nilanjan Chaudhuri, Sourav Sarkar, Pradip Roy

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

In the extreme heat of a particle collision, matter behaves in ways that defy our everyday experience. When heavy atomic nuclei smash together at nearly the speed of light, they create a fleeting, super-hot soup known as a quark-gluon plasma. In this state, the protons and neutrons that usually make up atomic nuclei melt apart, freeing their internal constituents: quarks and gluons. For a brief moment, this plasma exists at temperatures trillions of degrees hotter than the center of the sun. Scientists study this primordial state of matter to understand how the universe looked just microseconds after the Big Bang. One of the most valuable tools for studying this invisible fireball is light itself. Because photons, the particles of light, do not interact strongly with the plasma, they escape the collision zone immediately, carrying a direct snapshot of the conditions inside. By measuring the energy and number of these escaping photons, researchers can reconstruct the temperature and density of the plasma. However, a complete picture requires understanding not just the heat, but also the density of matter within the plasma. In many modern experiments, the collision creates a region where there are more matter particles than antimatter particles, a condition described by a non-zero chemical potential. Understanding how this imbalance affects the light emitted is crucial for interpreting data from current and future experiments designed to map the properties of this exotic matter.

For years, physicists have been able to calculate how many photons are produced by this hot plasma, but a significant gap remained in their understanding when the plasma is dense. Theoretical calculations for photon production are naturally divided into two parts: a "hard" component involving high-energy collisions between particles, and a "soft" component involving lower-energy interactions. To make these calculations work, scientists must introduce an arbitrary dividing line, or cutoff, to separate the hard interactions from the soft ones. While it was known that the final answer should not depend on where this line is drawn, proving this mathematically had only been possible when the plasma was perfectly balanced between matter and antimatter. When a chemical imbalance was introduced, previous attempts to prove this independence analytically failed, leaving researchers to rely on numerical approximations that were not fully rigorous. This uncertainty meant that the theoretical tools used to interpret experimental data were not as robust as they needed to be, particularly for the dense environments created in recent beam-energy scan programs.

A team of researchers has now closed this gap by demonstrating, for the first time, that the mathematical independence of this dividing line holds true even when the plasma is dense and imbalanced. They achieved this by carefully re-evaluating the hard component of the photon production rate. In their approach, they treated the particles involved in the high-energy collisions using a simplified distribution that is valid when the photon energy is much larger than the temperature of the plasma. This allowed them to perform the complex integrations required to find the photon rate without getting stuck in mathematical dead ends that had plagued previous attempts. By combining this new, precise calculation for the hard part with the already known formula for the soft part, they showed that the arbitrary cutoff parameter cancels out perfectly. The result is a single, clean mathematical expression that describes the total rate of photon production without relying on any arbitrary choices or numerical shortcuts.

The researchers found that their new formula reveals a clear physical effect: the presence of a chemical imbalance, or a surplus of matter over antimatter, increases the number of photons produced. This happens because the imbalance changes the availability of quarks and antiquarks, effectively increasing the number of particles ready to collide and emit light. When they compared their new analytical result with older, approximate formulas that had been widely used to interpret data, they found that the older methods deviated significantly, especially at lower photon energies and higher densities. Their new formula provides a more accurate description of the physics, particularly in the regime where the chemical potential is comparable to the temperature. This level of precision is vital as experiments move toward exploring the densest regions of the phase diagram.

Because the final result is a closed-form analytical expression, it can be directly plugged into the computer simulations that model the entire evolution of a heavy-ion collision. This removes the need for computationally expensive numerical calculations or the systematic errors introduced by older approximations. The ability to use this precise formula in real-time simulations will allow physicists to extract more accurate information about the temperature and density of the quark-gluon plasma from experimental data. This advancement is particularly relevant for upcoming experiments at facilities like FAIR and NICA, which are specifically designed to probe matter at high densities. By providing a reliable theoretical foundation, this work ensures that the light emitted from these collisions can be read correctly, helping scientists to map the true nature of matter under the most extreme conditions in the universe.

Technical Summary: Thermal Photon Production Rate from QGP at Finite Density

Problem Statement
The production of thermal photons from a Quark-Gluon Plasma (QGP) serves as a critical probe for relativistic heavy-ion collisions, particularly in experiments where the net baryon density is non-vanishing (e.g., RHIC beam-energy scan, FAIR, NICA). While the thermal photon production rate has been well-established at zero chemical potential ( $\mu=0$ ), extending this framework to finite quark chemical potential ( $\mu \neq 0$ ) presents a significant theoretical challenge.

The calculation of the photon rate involves separating the process into "hard" (large momentum exchange) and "soft" (small momentum exchange) regimes, separated by an arbitrary intermediate momentum cutoff $k_c$ . While the mutual cancellation of the $k_c$ dependence between these two regimes is required for a physical result, previous studies at finite density [1, 10] could only demonstrate this cancellation numerically. A long-standing issue, highlighted in Ref. [18], was that attempts to derive a closed-form analytic expression for the hard contribution at finite $\mu$ resulted in terms containing branch points (specifically involving exponential integral functions $Ei$ ), which prevented a fully analytic demonstration of cutoff independence.

Methodology
The authors compute the thermal photon emission rate $E dR/d^3p$ using the imaginary part of the retarded photon self-energy. The calculation follows the standard weak-coupling framework:

Soft Contribution: Calculated using the Hard Thermal Loop (HTL) approximation for the soft quark propagator. The result depends on the thermal quark mass $M$ , which is derived as a function of temperature $T$ and chemical potential $\mu$ using polylogarithm functions ( $Li_2$ ).
Hard Contribution: Calculated from leading-order perturbative QCD processes: quark-antiquark annihilation ( $q\bar{q} \to g\gamma$ $q \overset{q}{ˉ} \to g γ$ ) and QCD Compton scattering ( $qg \to q\gamma$ $q g \to q γ$ and $\bar{q}g \to \bar{q}\gamma$ $\overset{q}{ˉ} g \to \overset{q}{ˉ} γ$ ).
- The authors employ the Maxwell-Boltzmann (MB) approximation for the distribution functions of incoming particles, valid in the kinematic regime where the photon energy $E$ is much larger than the temperature ( $E \gg T$ ).
- Crucially, the authors re-evaluate the integration over the Mandelstam variable $s$ for the Compton scattering process. They demonstrate that the derivation in Ref. [18] was inconsistent. By expanding the logarithmic terms in the integrand and utilizing series expansions of the polylogarithm functions, they derive a closed-form analytic expression for the hard contribution that remains well-behaved as $\mu \to 0$ .

Key Contributions and Results

Analytic Cutoff Independence: The primary contribution of this work is the first analytic demonstration that the dependence on the arbitrary cutoff $k_c$ cancels exactly between the soft and hard contributions, even in the presence of a non-zero quark chemical potential $\mu$ . This resolves the issue previously addressed only numerically.
Generalized Analytic Formula: The authors provide a closed-form expression for the total photon production rate (Eq. 24) that incorporates finite $\mu$ . This expression is a generalization of the results from Refs. [2, 3].
Small- $\mu$ Expansion: A compact expansion for small $\mu$ (Eq. 25) is derived, which isolates a new $O(\mu^2)$ term. This term represents the specific enhancement due to finite baryon density.
Validation: The analytic hard contribution (Eq. 23) was compared against a semi-analytical approach (analytic integration of $E_2, V, t$ but numerical integration of $s$ ) and fully numerical results. The analytic expression shows excellent agreement for physically relevant values where $k_c \ll T$ .
Comparison with Parametrizations: The new analytic result (Eq. 24) and its small- $\mu$ approximation (Eq. 25) were compared to the phenomenological parametrization by Traxler et al. [1]. The authors find that their analytic result agrees more systematically with the exact numerical rate, particularly at lower photon energies and moderate $\mu/T$ ratios, where the Traxler parametrization shows deviations up to 7%.

Significance and Claims
The paper claims to resolve a long-standing theoretical issue regarding the analytic treatment of thermal photon rates in dense QGP. By providing a closed-form, cutoff-independent analytic expression, the authors enable the direct implementation of these rates into relativistic hydrodynamic simulations of heavy-ion collisions without the computational overhead of full numerical integration or the systematic biases of older parametrizations.

The results indicate that a finite quark chemical potential enhances the thermal photon yield by modifying the phase-space density of quarks and antiquarks. This enhancement is monotonic with respect to $\mu$ and is particularly significant at lower photon energies. The authors state that their framework is immediately applicable to the extraction of direct-photon spectra and flow observables in baryon-rich matter, relevant for current and future experiments at RHIC, FAIR, and NICA. The work also sets the stage for future extensions, such as incorporating chirally imbalanced media and next-to-leading-order QCD corrections.

Technical Summary: Thermal Photon Production Rate from QGP at Finite Density

More like this