Multi-zone Modeling of Blazar Jets: Constraints from GeV-Optical Correlation and Short-Timescale Variability

Abstract: We have developed a multi-zone model of blazar jet emission, in which the emission region contains many cells with individual magnetic fields and electron energy distributions. Nonthermal emission from radio to $γ$-rays is generated by electrons accelerated by shocks passing through the region via synchrotron and inverse-Compton (IC) processes. The optical and GeV variability at days-to-months time-scale simulated from our model are strongly correlated with no significant time lag, as observed in most blazars and indicated by the standard shock-in-jet model. However, the mechanism of the shorter time-scale variability has been less explored, although such fluctuations at X-rays, $γ$-rays and optical bands have been observed regularly in recent years. In our model, the hr time-scale variability of the synchrotron radiation is due to the spatial fluctuation of the magnetic field in the emission region. We found that to reproduce the short-timescale variability of the observed synchrotron emission in blazars, the required fluctuations of the magnetic field are in the range $1-2\%$ to $25-30\%$. Similar variability of the IC emission, which does not depend on the magnetic field, may be reproduced in our model by implementing equipartition of energy between the magnetic field and particles. We found that orphan flares in the optical or GeV band, or optical-GeV correlation with a significant time delay, as observed occasionally, may be reproduced in certain special conditions related to the orientation of the magnetic field in the cells.

• The paper develops a multi-zone leptonic shock-in-jet model showing that 1–30% spatial magnetic-field fluctuations can reproduce hour-scale synchrotron variability, while equipartition transfers these fluctuations into EC-dominated GeV emission.
• The model reproduces observed variability trends, including red-noise power spectra, strong zero-lag optical–GeV correlations, and flux–flux slopes near 2 for SSC-dominated BL Lacs and 1.5 for mixed SSC–EC FSRQs.
• The paper finds that rare optical or GeV orphan flares and optical–GeV lags of up to tens of days require specific magnetic-field geometries, while limited light-travel treatment and injection-only equipartition remain important uncertainties.

Bala, Mitra and Chatterjee present a multi-zone leptonic model of blazar jet emission designed to reproduce multi-wavelength variability from months down to hour timescales, and use it to constrain the physical origin of short-timescale fluctuations and of the occasionally anomalous optical–GeV correlation behavior (2608.18707). The central result is that hour-scale synchrotron variability requires spatial magnetic field fluctuations of roughly 1–2% to 25–30%, and that hour-scale variability in external-Compton (EC) dominated GeV light curves arises indirectly through equipartition between magnetic and particle energy densities.

Model architecture

The emission region is a cylinder sliced into cells along the jet axis, each with its own magnetic field magnitude and orientation and its own electron energy distribution evolving independently under radiative cooling. A mildly relativistic shock ($\beta_{\rm shock}=0.2$, $\Gamma_{\rm shock}\approx1.02$) traverses the region, injecting electrons with a power-law distribution $N(\gamma)=N_0\gamma^{-s}$ with $s=2.05$. The magnetic field in each cell is decomposed into a smoothly varying component (linear between input endpoint values) plus a fluctuating component drawn from uniform random distributions, with separate normalization constants for inter-zone and intra-zone fluctuations; the field direction varies cell-to-cell by a controlled angle increment $\delta\theta$.

The model adopts fixed bulk parameters: $\Gamma=15$, viewing angle $0.05$ rad ($\delta=20$), jet radius 0.001 pc, and equipartition ($U_B/U_e=1$) applied at the moment of shock injection into each cell. The electron density normalization $N_0$ is then set by $\Gamma_{\rm shock}\approx1.02$0, $\Gamma_{\rm shock}\approx1.02$1, and $\Gamma_{\rm shock}\approx1.02$2. Cooling includes synchrotron losses and IC losses against BLR and torus photon fields parameterized following Hayashida et al., with disk luminosities spanning $\Gamma_{\rm shock}\approx1.02$3–$\Gamma_{\rm shock}\approx1.02$4 erg s$\Gamma_{\rm shock}\approx1.02$5 for FSRQs. Both SSC (with light-travel-time effects included only along the jet axis) and EC are computed using the full Klein-Nishina cross-section. Synchrotron self-absorption is omitted, justified by the focus on optical and GeV bands.

Origin of short-timescale variability

The key physical claim is twofold. First, hour-scale fluctuations of the synchrotron emission are driven directly by spatial fluctuations of the magnetic field among cells; reproducing the observed range of sub-day synchrotron variability requires fluctuation amplitudes of $\Gamma_{\rm shock}\approx1.02$61–2% up to $\Gamma_{\rm shock}\approx1.02$725–30%. Second, since EC emission does not depend explicitly on $\Gamma_{\rm shock}\approx1.02$8, the comparable short-timescale fluctuations seen in GeV light curves of LSP/ISP sources such as PKS 1510-089, 3C 454.3, 3C 273, and S5 0716+714 must arise through a different channel. The authors show that imposing equipartition transfers magnetic-field fluctuations into electron energy-density fluctuations at injection, thereby producing EC variability; they verify that the same trend holds for non-equipartition partitions (1:2, 1:5, 5:1). They conclude that equipartition is a necessary ingredient of blazar jet energetics for reproducing small-amplitude, short-timescale GeV fluctuations in EC-dominated sources — a nontrivial constraint, since it elevates equipartition from an energetic convenience to a variability-motivated requirement.

Simulated fractional rms variability amplitudes are broadly consistent with Fermi-LAT observations: simulated LSP sources show $\Gamma_{\rm shock}\approx1.02$9 on 5–10 day timescales versus 0.18 on 50–100 days, while HSP simulations yield 0.16 and 0.17 respectively, bracketing the observed FSRQ values (0.11, 0.26) and BL Lac values (0.26, 0.36). The computed PSDs are red-noise over the simulated frequency range, consistent with observations and with other time-dependent leptonic simulations. Notably, increasing the magnetic fluctuation amplitude raises high-frequency PSD power by approximately one order of magnitude while leaving low-frequency power nearly unchanged, implying that the relative amplitude of short-timescale variability determines whether a broadband PSD shows a break — though the authors caution that their simulated light curves span too narrow a frequency range for a robust PSD fit.

Inter-band correlation, orphan flares, and lags

For all LSP, ISP, and HSP parameter sets, the simulated optical and GeV light curves are strongly correlated with zero lag, matching the dominant observational pattern from large-sample studies. On log-log optical-versus-GeV flux plots, the model naturally produces slopes near 2.0 for BL Lac-type objects (SSC-dominated, quadratic in electron density) and an intermediate slope of $N(\gamma)=N_0\gamma^{-s}$01.5 for FSRQs (mixed SSC+EC), consistent with the observed dichotomy of slopes near 1.0 and 2.0 reported by Liodakis et al.

Orphan flares emerge under specific magnetic geometry conditions. An optical orphan flare occurs when turbulent field fluctuations modulate $N(\gamma)=N_0\gamma^{-s}$1 without significantly changing $N(\gamma)=N_0\gamma^{-s}$2: if the mean field angle to the line of sight exceeds $N(\gamma)=N_0\gamma^{-s}$3, large angular fluctuations alter $N(\gamma)=N_0\gamma^{-s}$4 (and hence synchrotron flux) while leaving $N(\gamma)=N_0\gamma^{-s}$5 and the equipartition-driven EC output nearly unchanged. A GeV orphan flare arises when $N(\gamma)=N_0\gamma^{-s}$6 is tied to the parallel field component via $N(\gamma)=N_0\gamma^{-s}$7, motivated by shock acceleration efficiency depending on field-shock orientation; in ISP sources, where the synchrotron peak lies near $N(\gamma)=N_0\gamma^{-s}$8 Hz, resulting $N(\gamma)=N_0\gamma^{-s}$9 excursions produce UV/X-ray and GeV flares with little optical counterpart.

Non-zero optical–GeV lags of up to tens of days appear when the angular variation per cell is large ($s=2.05$0), whereas variations of $s=2.05$1 per cell yield zero lag. When $s=2.05$2 and $s=2.05$3 fluctuate comparably but out of phase (approximately sinusoidally with a phase offset), the synchrotron and EC light curves inherit a lag whose sign depends on which correlation peak dominates. The authors emphasize that reproducing orphan flares and significant lags requires fine-tuned conditions, consistent with their rarity in observations.

Limitations and open questions

Several caveats bear directly on the quantitative claims. Light-travel-time effects are treated only along the jet axis; transverse delays are neglected, which the authors argue is acceptable for highly collimated regions and for the optical-synchrotron versus EC-GeV comparison, but this assumption would need revisiting for polarization or more compact geometries. The equipartition condition is imposed instantaneously at injection rather than evolved self-consistently, so the claimed link between equipartition and EC variability rests on the injection prescription. The PSD analysis is limited to the narrow Fourier band spanned by the simulations, leaving open whether the observed scale-invariant variability from hours to years can be reproduced by extending the model across a wider dynamic range. Finally, no observational data were fitted directly; constraints come from comparison with published statistical properties, so a systematic Bayesian exploration of the parameter space against individual well-sampled multi-wavelength data sets remains untested.

Conclusion

This work provides a physically motivated account of hour-timescale blazar variability within a shock-in-jet framework: magnetic field fluctuations of 1–30% drive synchrotron fluctuations directly, while equipartition propagates those fluctuations into the EC-dominated GeV band. The same magnetic geometry controls the rarer anomalies — orphan flares and non-zero optical–GeV lags — explaining their infrequency as a consequence of required field orientations. The model's main open question is whether its fluctuation prescription, extended over a broader range of timescales and tested against individual sources, can account for the apparent scale invariance of blazar PSDs from hours to years.