X-ray Fourier lag-frequency spectra modulated by stochastic turbulent acceleration in the jets of high-frequency-peaked BL Lac

Technical Summary: X-ray Fourier Lag-Frequency Spectra Modulated by Stochastic Turbulent Acceleration in HBL Jets

Problem Statement
X-ray interband time lags serve as critical diagnostics for jet physics in High-frequency-peaked BL Lac (HBL) objects. Observations frequently reveal diverse lag signatures, including hard (positive) lags where high-energy variations lag low-energy ones, and soft (negative) lags where the reverse occurs. While previous theoretical models based on Fourier-transformed electron transport equations could reproduce soft lags or required generalized transport equations to explain hard lags, they often neglected nonlinear Synchrotron Self-Compton (SSC) cooling. Furthermore, the role of Stochastic Turbulent Acceleration (STA) in shaping these lag-frequency spectra, particularly in conjunction with nonlinear cooling effects, remained to be fully explored. This work addresses the need for a unifying quantitative framework that accounts for STA, radiative cooling, particle escape, and nonlinear SSC cooling to interpret the complex time-lag signatures observed in HBL flares.

Methodology
The authors employ a generic one-zone leptonic model assuming a spatially homogeneous, optically thin spherical plasma blob moving with a relativistic Lorentz factor. The intrinsic temporal evolution of the electron population is governed by a time-dependent Fokker-Planck (FP) equation. This equation incorporates:

1. Injection: A power-law distribution of electrons injected impulsively.
2. Acceleration: Stochastic turbulent acceleration (STA) modeled via a diffusion coefficient in energy space ($D_p \propto \gamma'^q$), with "hard-sphere" turbulence ($q=2$) as the primary case.
3. Cooling: Radiative losses from synchrotron emission and SSC scattering, the latter including nonlinear effects and Klein-Nishina corrections.
4. Escape: An energy-independent escape timescale.

The FP equation is solved numerically using an implicit Crank-Nicholson scheme. The resulting time-dependent electron distributions are used to calculate synchrotron and SSC emission spectra. To account for causality, light-crossing effects are incorporated by integrating flux contributions from slices of the emitting region.

Time lags between soft (0.5–2 keV) and hard (2–10 keV) X-ray light curves (LCs) are analyzed using three methods:

• Discrete Fourier Transform (DFT): To derive Fourier lag-frequency spectra.
• Normalized Flux Matching Delay (NFMD): To track instantaneous time delays.
• Cross-Correlation Function (CCF): To measure centroid lags over the full LC.

The study systematically varies parameters such as the maximum injected Lorentz factor ($\gamma'_{i,max}$), the spectral index ($p$), the diffusion coefficient ($D_0$), and the injection rate to simulate three distinct regimes: hard lags, soft lags, and a transitional state.

Key Contributions and Results

1. Unification of Lag Regimes: The model successfully reproduces three distinct time-lag regimes:

   • Hard/Positive Lags: Occur when the maximum injected energy is significantly below the equilibrium energy ($\gamma'_{i,max} \ll \gamma'_{eq}$), indicating acceleration-dominated dynamics.
   • Soft/Negative Lags: Occur when $\gamma'_{i,max} \gg \gamma'_{eq}$, where cooling dominates.
   • Transitional Regime: A smooth transition from positive to negative lags across frequencies is observed when $\gamma'_{i,max} \approx \gamma'_{eq}$. This transition reflects the coexistence of acceleration-dominated low-energy electrons and cooling-dominated high-energy electrons.

2. Role of STA: The competition between STA, radiative cooling, and escape processes shapes the lag-frequency spectra. Specifically, STA suppresses the radiative cooling of high-energy electrons by continuously re-accelerating lower-energy electrons. This suppression amplifies time lags in the transitional and soft-lag regimes compared to pure cooling models.

3. Impact of Nonlinear SSC Cooling:

   • SSC cooling further amplifies time lags, particularly in the transitional and soft-lag regimes.
   • The amplification arises because STA sustains soft X-ray emission (via re-accelerated lower-energy electrons) while SSC cooling accelerates the decay of hard X-ray emission.
   • The model predicts a positive correlation between flare duration and lag amplitude: longer flares allow for greater accumulation of electrons, strengthening SSC cooling and resulting in larger lags. This explains observations of large lags in TeV-bright flares.

4. Methodological Comparison: The study highlights that CCF lags, which average over the full light curve, may obscure transitional behaviors and frequency-dependent variations that are clearly resolved by Fourier lag-frequency spectra. In asymmetric light curves (soft-lag cases), CCF lags can be significantly smaller than low-frequency Fourier lags.

5. Robustness: The qualitative features of the lag-frequency spectra remain consistent when varying the turbulence index from "hard-sphere" ($q=2$) to Kolmogorov ($q=5/3$), though quantitative differences in lag amplitude and characteristic frequency arise. The inclusion of Fermi-I (systematic) acceleration was found to have negligible effects in the hard-lag case but can modify spectra in other regimes; however, the main conclusions regarding STA and SSC remain robust.

Significance
The paper claims to provide a unifying quantitative framework for interpreting the diverse time-lag signatures observed in HBL X-ray flares. By integrating STA and nonlinear SSC cooling, the model explains:

• The existence of both hard and soft lags and the transition between them.
• The amplification of lags in TeV-bright flares.
• The observed trend where larger flare durations correspond to larger time lags.

The authors suggest that the transitional regime, characterized by a smooth shift from positive to negative lags, may serve as evidence for the re-acceleration of relativistic electrons. They conclude that Fourier lag-frequency spectra, particularly when combined with Spectral Energy Distribution (SED) modeling, offer a powerful diagnostic tool for constraining the extreme physical conditions (acceleration efficiency, magnetic field strength, and cooling rates) in the compact emission regions of relativistic jets. The work anticipates that future missions like Athena, with superior time and energy resolution, will enable precise measurements of these spectra to further test these theoretical predictions.