A Systematic Search for Spectral Hardening in Blazar Flares with the Fermi-Large Area Telescope

TL;DR

This work conducts the first population-level search for spectral hardening in γ-ray blazar flares using Fermi-LAT data, focusing on HSPs and EHSPs to test whether the GeV spectrum hardens during flares. By identifying flares with Bayesian Block Analysis and comparing broken and single power-law fits via a likelihood-ratio TS, the study finds only three marginally significant events across 4063 flares, implying a frequency below $0.1\%$. The rarity of such events supports a picture in which blazar emission is dominated by smoothly varying power-law spectra in the LAT band, while remaining open to rare, exceptional physical conditions such as multi-zone emission, magnetic reconnection, or hadronic contributions. The results establish population-level constraints and motivate targeted, multi-wavelength follow-ups, including VHE observations, to clarify the physical origins of spectral hardening in blazars.

Abstract

Blazars are a subclass of active galactic nuclei (AGN) that emit non-thermal radiation through relativistic jets, characterized by rapid flux and polarization variability. High synchrotron-peaked blazars (HSPs) and extreme high synchrotron-peaked blazars (EHSPs), with synchrotron peaks exceeding $10^{15}$ Hz and $10^{17}$ Hz, respectively, are crucial for understanding the full range of blazar phenomena and testing models of jet physics. Yet, their understanding remains challenging. This work aims to systematically identify and characterize the most extreme $γ$-ray blazars using data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. The focus is on spectral hardening, where the $γ$-ray spectrum becomes harder at higher energies, particularly during flaring episodes. This represents the first dedicated analysis of spectral hardening across a population of EHSPs, as previous studies explored it only in individual sources. We analyze 138 blazars selected from the 4FGL-DR2 catalog with high synchrotron peak frequencies and well-sampled light curves. Flaring periods are automatically identified, and each flare is analyzed, with the significance of spectral hardening assessed through a test statistic based on the likelihood ratio of two spectral models. We identify two flaring episodes with indications of spectral hardening, in 4FGL J0238.4$-$3116 and PKS 2155$-$304, the latter detected independently by both methods but referring to the same period. These events are consistent with expectations from statistical fluctuations, suggesting that spectral hardening is a rare occurrence (< 0.1 %). These results constrain its frequency and support a smoothly varying power-law blazar emission model, motivating future multi-wavelength studies to clarify whether these rare flares reflect distinct physical processes within blazar jets.

Figures (3)

• Figure 1: Light curve of the blazar 4FGL J0021.9$-$5140, as an example, from Fermi-LAT data (black flux measurements with error bars) with flare detection methods using the Baseline (top panel) and Flip (bottom panel) methods. We show the Bayesian blocks (blue line) and flares (shaded boxes, with different colors only for visualization convenience).
• Figure 2: Light curves of the sources where a spectral hardening feature is detected with TS$_{\mathrm{hardening}} \geq 12$. Insets show the Bayesian Blocks segmentation (blue line), with identified flares shaded alternately in purple and green for clarity. Flares showing spectral hardening are highlighted in red, with the 4FGL source name and corresponding time range labeled above each inset.
• Figure 3: Gamma-ray spectra for the flaring states with TS$_{\mathrm{hardening}} \geq 12$. The Fermi-LAT data (black dots), with the best-fit BPL model and its uncertainties (red). The PKS 2155$-$304 flare identified by the Flip method (top right) partially overlaps with that selected by the Baseline method (lower left) but spans twice the duration, resulting in a more significant SED for the Flip detection. Low-energy data are more uncertain, though we minimized systematics by using PSF event types and by identifying spectral hardening through likelihood ratios rather than individual flux bins.