Particle acceleration signatures in the time-dependent one-zone synchrotron self-Compton model of blazar flares

The study of multiwavelength flux and spectral variations during rapid flares from blazars provides strong constraints on the physical parameters of the compact emission regions responsible for these still poorly understood events. Although a full description of the continuous and transient emission from blazars seems to require more sophisticated scenarios, particle acceleration and loss mechanisms can be approximately described within the simple leptonic one-zone framework, enabling a systematic study of their impact on the observable properties of multiwavelength flare light curves. Our goal is to identify characteristic signatures in these light curve profiles that permit one to discriminate between the main physical processes situated inside the relativistic jet and commonly invoked to explain blazar flares. The present study exclusively focuses on modeling rapid flares from BL Lac type objects, which can be described within the synchrotron self-Compton (SSC) emission scenario. Combinations of several commonly employed mechanisms to describe the gain and loss of energetic particles in onezone models during flaring events are studied in a systematic way: particle injection; diffusive shock and stochastic acceleration and reacceleration; particle escape; adiabatic losses; radiative losses through synchrotron and inverse-Compton radiation. The current study is limited to the case of “hard-sphere” scattering. A large variety of light curve shapes arises from the different scenarios under study. Characteristic signatures, in particular energy-dependent time delays and differences in the shapes of the rising part of the flare, should allow the distinction to be made between different injection and acceleration scenarios, given the availability of sufficiently high-quality multiwavelength data sets. This is illustrated with a simplified application to a flare event from the blazar Mrk 421.

💡 Research Summary

This paper presents a systematic investigation of rapid flares in BL Lac type blazars using a time‑dependent, one‑zone synchrotron self‑Compton (SSC) framework. The authors employ the EMBLEM code to evolve the electron distribution Nₑ(γ,t) inside a homogeneous spherical blob, solving a Fokker‑Planck equation that includes synchrotron and inverse‑Compton cooling, adiabatic expansion, particle escape, and two forms of acceleration: diffusive shock (Fermi I) and stochastic (Fermi II). All processes are treated under the “hard‑sphere” scattering assumption, i.e., energy‑independent diffusion and escape coefficients, which simplifies the parameter space while retaining the essential physics.

The study explores a suite of scenarios combining different particle injection and (re)acceleration mechanisms: (i) a simple power‑law injection, (ii) shock‑driven injection with a time‑dependent maximum Lorentz factor, (iii) stochastic acceleration with a constant acceleration timescale, and (iv) re‑acceleration of an already injected power‑law population by either Fermi I or Fermi II. Each scenario is further examined with or without adiabatic expansion and with varying escape timescales, yielding more than ten distinct model configurations.

For each configuration the authors compute the evolving spectral energy distribution (SED) and multi‑wavelength light curves (LCs) in radio, optical, X‑ray, and γ‑ray bands. They focus on four observable diagnostics: (a) the asymmetry between rise and decay times (τ_rise vs τ_decay), (b) energy‑dependent time delays between the peaks in different bands (Δt(E)), (c) the ratio of peak fluxes across bands, and (d) the presence of a plateau or extended tail in the high‑energy decay. The results reveal clear, model‑specific signatures. Pure injection flares produce a fast low‑energy rise followed by a delayed high‑energy peak, with roughly symmetric decay. Shock (Fermi I) acceleration leads to high‑energy peaks that precede the low‑energy ones (negative Δt), while stochastic (Fermi II) acceleration yields simultaneous rises across the spectrum and a slower, more gradual decay. Re‑acceleration hardens the electron spectrum, shifting both synchrotron and SSC peaks to higher energies and often generating a pronounced high‑energy plateau when adiabatic expansion is included. Short escape times produce steep, rapid decays; long escape times generate lingering emission.

The authors condense these findings into a “signature matrix” that maps observable LC features to underlying physical processes. To demonstrate practical applicability, they apply the matrix to a well‑studied 2013 flare of Mrk 421. By fitting the observed multi‑band light curves, they find that a model combining stochastic (Fermi II) acceleration with moderate adiabatic expansion best reproduces the measured ~0.5 h inter‑band delay and the asymmetric rise/decay profiles.

In conclusion, even within the simplicity of a one‑zone SSC model, distinct acceleration and loss mechanisms imprint recognizable patterns on flare light curves. The paper argues that forthcoming high‑cadence, broadband monitoring campaigns (e.g., CTA, LSST, IXPE) will be able to exploit these signatures to discriminate between competing flare scenarios, thereby constraining jet microphysics such as magnetic field geometry, shock properties, and turbulence spectra. The work thus provides both a theoretical toolkit and a roadmap for interpreting rapid blazar variability in the era of next‑generation observatories.