🔭 astrophysics

Synchrotron Emission from Cooled Particle Distributions

This paper derives accurate and computationally efficient analytic fitting functions for synchrotron emission and absorption coefficients that account for electron cooling effects, extending previous local treatments to improve modeling of astrophysical sources like gamma-ray bursts and tidal disruption events.

Original authors: Ross Ferguson, Ben Margalit

Published 2026-07-16

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: Ross Ferguson, Ben Margalit

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 the universe as a giant, chaotic dance floor where particles are constantly being spun up to incredible speeds. When these particles, specifically electrons, zoom through invisible magnetic fields, they don't just move; they scream. This scream is a type of light called synchrotron radiation, which we can see as radio waves, X-rays, or visible light depending on how fast the electrons are going. This happens everywhere, from the violent collisions of stars to the explosive deaths of massive stars in events like gamma-ray bursts.

To understand what we see from these cosmic fireworks, scientists have to figure out exactly how the electrons are moving. Usually, they assume the electrons are like a crowd of runners: some are sprinting, some are jogging, and the number of runners at each speed follows a predictable pattern. But there's a catch: these electrons are tired. As they run, they lose energy. They radiate light (which takes energy away) and the space they are running through expands, stretching them out like taffy. This "cooling" changes the crowd's speed distribution. If you don't account for this tiredness, your prediction of the light we see will be wrong. For a long time, scientists had to use slow, heavy computer simulations to track every single tired electron, which made modeling these cosmic events very difficult and computationally expensive.

This paper by Ross Ferguson and Ben Margalit offers a clever shortcut. Instead of running a slow, heavy simulation for every single electron, the authors have derived a set of "cheat codes"—mathematical formulas that act like a fast, accurate map of how these cooling electrons behave. They took two common types of electron crowds (one that looks like a standard power-law distribution and another that looks like a hot, thermal soup) and figured out exactly how their energy changes as they cool down. They then created simple, easy-to-use fitting functions that describe the light these cooled crowds emit and absorb.

The main finding is that these new formulas are just as accurate as the slow, heavy computer simulations but are much faster to calculate. The authors tested these formulas by applying them to a model of a gamma-ray burst afterglow and found that their results matched the complex, established models almost perfectly, with errors usually less than 1%. They explicitly show that a simpler, older way of thinking about these electrons—where you just chop off the high-energy tail of the crowd abruptly—fails to capture the smooth, gradual way electrons actually cool down, especially when the cooling is strong. While the paper doesn't claim to have solved every mystery of the universe, it provides a highly reliable and efficient tool that allows astronomers to model the light from cosmic shocks with much greater speed and consistency, helping us better understand the violent events that light up our sky.

Technical Summary: Synchrotron Emission from Cooled Particle Distributions

Problem Statement
Synchrotron emission from relativistic electrons is a primary diagnostic tool for high-energy astrophysical sources, including gamma-ray bursts (GRBs), luminous fast blue optical transients (LFBOTs), and jetted tidal disruption events (TDEs). In these environments, electrons are accelerated at shock fronts and subsequently cool via radiative losses (synchrotron and inverse-Compton) and adiabatic expansion. This cooling shifts electrons in energy space, fundamentally altering the electron distribution function and, consequently, the synchrotron emission and absorption coefficients.

While "one-zone" models often assume uniform cooling and simplified distributions, "full-volume" models treat the three-dimensional structure of the shock and downstream fluid self-consistently. In full-volume models, the downstream electron distribution becomes highly spatially dependent due to varying cooling histories. Previous approaches to modeling these local effects have relied on numerical integration to compute the cooled distribution function and associated radiation coefficients at every point behind the shock. This reliance on numerical methods significantly increases the computational cost of simulating these systems, hindering the efficiency of full-volume modeling.

Methodology
The authors extend the "local" treatment of synchrotron emission by deriving analytic expressions for the electron distribution function and the resulting emission and absorption coefficients, accounting for both adiabatic and synchrotron cooling.

General Formalism: The paper begins by deriving the evolution of an arbitrary, impulsively injected electron distribution function, $(\partial n / \partial \gamma)_{\text{inj}}$ , as it is advected downstream. Assuming the total number of electrons in a fluid element is conserved, the authors relate the distribution at time $t$ to the injection distribution via the ratio of number densities and the mapping of Lorentz factors $\gamma_{\text{inj}} \to \gamma(t)$ . The cooling rate includes terms for synchrotron losses and adiabatic expansion (PdV work), governed by characteristic timescales $t_B$ and $t_{\text{ad}}$ .
Specific Cases: The authors specialize this general formalism to two injected distributions:
- A power-law distribution: $(\partial n / \partial \gamma)_{\text{inj, pl}} \propto \gamma^{-p}$ .
- A relativistic Maxwellian (thermal) distribution.
  Analytic expressions for the "cooled" distribution functions are derived for both cases. These expressions include cutoff factors proportional to $(1 - \gamma/\gamma_\infty)$ , where $\gamma_\infty$ represents the Lorentz factor an electron injected with infinite energy would possess at time $t$ .
Derivation of Fitting Functions: The authors calculate the pitch-angle-averaged emissivity ( $j_\nu$ $j_{ν}$ ) and absorption coefficient ( $\alpha_\nu$ $α_{ν}$ ) by integrating the cooled distributions against the synchrotron kernel. Because these integrals cannot be solved in closed form for all frequency regimes, the authors derive analytic fitting functions.
- Power-Law: The fitting functions for the power-law case are constructed by joining asymptotic limits: low frequencies ( $x_1 \ll 1$ ), intermediate frequencies ( $1 \le x_1 < \eta^2$ ), and high frequencies ( $x_\infty \gg 1$ ). The authors use the method of steepest descent for the high-frequency regime and hypergeometric functions for the low-frequency regime. The transition between these regimes is handled by sigmoid functions with coefficients fitted to numerical integrations for $2 < p \le 5$ .
- Thermal: Similarly, fitting functions for the thermal distribution are derived by interpolating between low-frequency limits (involving confluent hypergeometric functions) and high-frequency limits (using a steepest descent approximation for the saddle point of the integral).

Key Contributions and Results

Analytic Distribution Functions: The paper provides exact analytic forms for cooled power-law (Eq. 11) and thermal (Eq. 12) electron distributions. These forms explicitly account for the "pile-up" of electrons near the maximum Lorentz factor $\gamma_\infty$ due to cooling, a feature often neglected in sharp-cutoff models.
Analytic Fitting Functions: The primary contribution is the set of analytic fitting functions for the synchrotron emissivity and absorption coefficients for both power-law and thermal distributions.
- For power-law distributions, the fitting functions (Eqs. 33, 52) reproduce numerical integrations with typical relative errors at the percent level, with maximum errors generally below 15% (peaking around $\eta \approx 1.14$ ).
- For thermal distributions, the fitting functions (Eqs. 79, 80) achieve mean relative errors $\lesssim 10\%$ , with peak errors of 10–25% occurring near the turnover frequency where low- and high-frequency asymptotes meet.
Comparison to Sharp-Cutoff Models: The authors demonstrate that approximating the cooled distribution as a "chopped-off" power-law (sharp cutoff at $\gamma_\infty$ ) introduces significant errors, particularly at high frequencies and when the cooling range $\eta = \gamma_\infty/\gamma_1$ is small. The smooth cutoff derived from the cooling ODE is shown to be necessary for accurate local modeling.
Application to GRB Afterglows: The authors test the power-law fitting functions by inserting them into a full-volume numerical framework for GRB afterglows (Ferguson & Margalit 2026) and comparing the results to the analytic spectrum of Granot & Sari (2002). The comparison shows that the new fitting functions reproduce the standard "slow cooling" spectrum with relative errors consistently below 10% (typically $\sim 1\%$ ), validating their use as a substitute for computationally expensive numerical integrals.

Significance and Claims
The paper claims that these analytic fitting functions offer a computationally efficient and accurate alternative to numerical integration for calculating synchrotron radiation coefficients in full-volume models. By avoiding the need to solve the cooling integrals numerically at every spatial point, these functions allow for the consistent treatment of cooling in synchrotron spectra without the prohibitive computational cost of previous methods.

The authors emphasize that while one-zone models often approximate cooling by enforcing a steeper spectral index, full-volume models require the averaging of contributions from fluid elements with varying degrees of cooling. The derived fitting functions enable this averaging process, facilitating the calibration of global one-zone models against more precise local radiation properties. The authors note that their formalism is applicable to any impulsively injected electron population cooling over time, though they explicitly restrict their current analysis to regimes where electrons remain relativistic ( $\gamma \gg 1$ ) and the power-law index $p > 2$ . They acknowledge that effects such as inverse-Compton cooling can be trivially included, while other processes like Coulomb cooling or self-absorption heating would require further modification to the formalism.

More like this