Signatures of a Maxwellian Component in Shock-Accelerated Electrons in GRBs

Recent particle-in-cell simulations suggest that a large fraction of the energy dissipated in a relativistic shock is deposited into a Maxwellian distribution of electrons that is connected to the high-energy power-law tail. Here, we explore the observational implications of such a mixed thermal-nonthermal particle distribution for the afterglow and prompt emission of gamma-ray bursts. When the Maxwellian component dominates the energy budget, the afterglow lightcurves show a very steep decline phase followed by a more shallow decay when the characteristic synchrotron frequency crosses the observed band. The steep decay appears in the X-rays at ~100 sec after the burst and is accompanied by a characteristic hard-soft-hard spectral evolution that has been observed in a large number of early afterglows. If internal shocks produce a similar mixed electron distribution, a bump is expected at the synchrotron peak of the nu*f_nu spectrum.

💡 Research Summary

The paper investigates the observational consequences of a mixed thermal‑non‑thermal electron distribution that emerges from recent particle‑in‑cell (PIC) simulations of relativistic shocks. These simulations indicate that, in a weakly magnetized electron‑ion shock, only a few percent of electrons are accelerated into a high‑energy power‑law tail while the bulk of the dissipated energy is deposited into a relativistic Maxwellian (thermal) component whose temperature is of order the bulk Lorentz factor Γ of the shocked plasma.

The authors model the downstream electron spectrum as a continuous function that joins a relativistic Maxwellian, (N_{\rm th}(\gamma)\propto \gamma^{2}\exp(-\gamma/\Theta)), to a power‑law, (N_{\rm pl}(\gamma)\propto \gamma^{-p}), at a transition Lorentz factor (\gamma_{\rm nth}). Two key parameters control the shape: (i) (\delta), the fraction of the total electron energy residing in the non‑thermal tail, and (ii) (\epsilon_{e}), the fraction of the shock‑dissipated energy transferred to electrons. For a given bulk Lorentz factor (\Gamma), the temperature (\Theta) and transition Lorentz factor (\gamma_{\rm nth}) are obtained by solving the energy‑budget equations. When (\delta=1) the distribution reduces to the classic pure power‑law; when (\delta\ll1) the Maxwellian dominates and (\Theta\simeq \epsilon_{e}\Gamma m_{p}/(3m_{e})).

Using standard synchrotron formulas (magnetic energy fraction (\epsilon_{B}=10^{-2})), the authors compute the emitted spectra for both slow‑cooling ((\nu_{c}>\nu_{\rm ch})) and fast‑cooling ((\nu_{c}<\nu_{\rm ch})) regimes, varying (\delta) over several orders of magnitude. The main spectral signature of a dominant Maxwellian is a pronounced break at the characteristic synchrotron frequency (\nu_{\rm ch}). Below (\nu_{\rm ch}) the emission is set by the thermal electrons and is relatively flat; just above (\nu_{\rm ch}) the flux drops sharply because the Maxwellian contributes little, and only at higher frequencies does the power‑law tail re‑emerge, producing a hardening of the spectrum. This “sharp decline‑then‑hardening” is not a simple power‑law break and is a direct imprint of the mixed distribution.

To translate these spectral features into light curves, the authors adopt the standard blast‑wave dynamics for a spherical outflow of isotropic equivalent energy (E=10^{53}) erg expanding into either a constant‑density interstellar medium ((n_{\rm ext}=1) cm(^{-3})) or a wind‑like medium ((n\propto R^{-2})). The bulk Lorentz factor evolves as (\Gamma\propto R^{-3/2}) (ISM) or (\Gamma\propto R^{-1/2}) (wind). With fiducial parameters (\epsilon_{e}=0.3), (\epsilon_{B}=0.01), and electron index (p=2.5), they calculate multi‑band light curves (optical, X‑ray at 1 keV, and radio).

For (\delta\gtrsim0.3) the light curves show the usual single break when (\nu_{\rm ch}) passes through the observing band. However, for (\delta\lesssim0.1) a distinct two‑stage behavior appears: as (\nu_{\rm ch}) sweeps across the X‑ray band (typically a few × 10² s after trigger), the flux undergoes a very steep decline (temporal index (\alpha\sim3–5)), accompanied by a characteristic “hard‑soft‑hard” spectral evolution (photon index softening then re‑hardening). After the steep phase, the non‑thermal tail begins to dominate the observed band, leading to a shallower decay. This pattern matches the early steep X‑ray decays observed by Swift/XRT in many bursts, which often display spectral evolution inconsistent with the high‑latitude emission model.

The authors also discuss implications for the prompt emission. If internal shocks generate a similar mixed electron distribution, the synchrotron (\nu f_{\nu}) spectrum would exhibit a bump or “bump‑like” excess at the peak frequency, deviating from the smooth Band‑function shape. Detection of such a feature would indicate efficient electron heating (large (\Theta)) and relatively weak non‑thermal acceleration.

In summary, the paper provides a quantitative framework linking the microphysics of relativistic shock acceleration (parameterized by (\delta) and (\epsilon_{e})) to observable signatures in GRB afterglows and prompt spectra. It shows that a dominant Maxwellian component naturally explains the early steep X‑ray decay and its associated spectral evolution, while also predicting a possible synchrotron bump in the prompt phase. This work extends the standard afterglow model by incorporating realistic electron heating found in modern PIC simulations, offering a more physically motivated interpretation of several puzzling GRB observations.