TeV pion bumps in the gamma-ray spectra of flaring blazars

TL;DR

This paper proposes a novel model explaining narrow TeV spectral features in flaring blazars as resulting from neutral pion decay, which could be detected with upcoming gamma-ray observatories like CTA.

Contribution

It introduces a new model for narrow VHE spectral features in blazars based on pion decay from proton-photon interactions during X-ray flares.

Findings

Narrow $$ TeV features can originate from $^0$ decay in proton-photon interactions.

The model fits the Mrk 501 flare spectrum well with a combined SSC and pion decay component.

Simulations suggest CTA can detect $^0$ bumps with 5$$ significance in 30 minutes.

Abstract

(Abridged) Very high-energy (VHE, $E>100$ GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions ($\pi^0$ bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. We explored the conditions needed for the emergence of narrow $\pi^0$ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies $\sim 100$ keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. We find that synchrotron-dominated flares with peak energies $\gtrsim100$ keV can be ideal periods for the search of $\pi^0$ bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct $\pi^0$ bump centered at 3 TeV. Spectral fitting of simulated SSC+$\pi^0$ spectra for the Cherenkov Telescope Array (CTA) show that a $\pi^0$ bump could be detected at a 5$\sigma$ significance level with a 30-min exposure.