Two-Zone Leptonic Model: Acceleration & Cooling

• The Two-Zone Leptonic Model is a framework that divides the emission region into a compact acceleration zone and a larger cooling zone to overcome the limits of one-zone models.
• It decouples acceleration from cooling by assigning different magnetic fields, particle distributions, and photon environments to each zone, thereby fitting hard X-ray and γ-ray spectra.
• The model explains energy-dependent lags, variability patterns, and distinct spectral features in blazars and supernova remnants by leveraging differing local conditions.

A two-zone leptonic model is a radiative-emission framework in which the observed broadband spectral energy distribution is produced not in a single homogeneous region, but by the superposition of radiation from two physically distinct zones with different magnetic fields, particle distributions, transport histories, or external photon environments. In the cited applications, the two zones are typically organized as a compact acceleration or high-field region and a larger cooling or downstream region, or as a quasi-stationary background component plus a compact variable component. This class of models has been developed most explicitly for blazars and supernova remnants, where it is used to relax the tensions that arise in one-zone interpretations of hard $\gamma$-ray spectra, thin X-ray rims, or energy-dependent timing behavior (Boula et al., 2021, Hu et al., 2024, Atoyan et al., 2011).

1. Definition and physical motivation

The basic motivation for a two-zone leptonic description is the repeated failure of a single homogeneous emitting region to account for all observed bands with the same particle population and the same local conditions. In extreme BL Lacs, the hardness of the $\gamma$-ray spectrum, with photon index $<2$, “demands extreme parameters” in standard one-zone SSC, including unusual Doppler factor, magnetic field, electron energies, and source compactness. In Tycho, a standard single-zone reading of the thin hard X-ray rim leads to $B \gtrsim 400\,\mu{\rm G}$ everywhere, which strongly suppresses inverse Compton emission and makes the observed GeV–TeV flux difficult to explain leptonically. In Mrk 421, a one-zone interpretation has difficulty accounting simultaneously for quasi-steady emission, rapid X-ray flares, and the systematic increase of lag magnitude with energy separation (Aguilar-Ruiz et al., 2021, Atoyan et al., 2011, Hu et al., 2024).

In response, the two-zone approach distributes the radiative burden across two regions with different efficiencies. One zone can retain the high magnetic field or rapid acceleration required by X-rays, while the other can provide a larger particle reservoir, weaker synchrotron losses, stronger bremsstrahlung or inverse Compton output, or a distinct variability timescale. This suggests that the central role of the model is not merely geometric subdivision, but the decoupling of acceleration, cooling, and target-photon conditions.

2. Canonical architectures and transport structure

In blazar implementations aimed at the blazar sequence, the architecture is an acceleration zone plus a downstream cooling zone. Zone I is a compact acceleration region near the central engine; Zone II is a larger region farther downstream, with

$R_{II} \simeq 10\,R_I.$

Electrons accelerated in Zone I escape into Zone II, where they continue radiating. The magnetic field decreases with distance as

$B \propto \frac{1}{z},$

and injection into the second zone is supplied by escape from the first,

$Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$

The external radiation field is attributed to disk photons scattered by an accretion-disk MHD wind, with

$$U_{\rm sc} \simeq \frac{L_{\rm disk}\,\tau_T}{4\pi R_2^2 c},\qquad \tau_T = n_0 \sigma_T R_s \ln\!\left(\frac{R_2}{R_1}\right),\qquad U_{\rm ext} = \Gamma^2 U_{\rm sc}.$$

In this formulation, the relative strengths of synchrotron, SSC, and external Compton are tied to the accretion environment rather than assigned independently (Boula et al., 2021, Boula et al., 16 Dec 2025).

A different but still explicitly leptonic architecture appears in time-dependent HBL modeling. There the two zones are a quasi-stationary spherical blob and a highly variable compact blob, both moving relativistically along the jet. The compact zone carries the time dependence, while the stationary zone reproduces the persistent baseline X-ray/TeV SED. The electron evolution includes stochastic acceleration, first-order Fermi acceleration, radiative cooling, particle escape, and time-dependent injection. For the stochastic term, the adopted hard-sphere turbulence prescription is

$$D_p(\gamma') = D_0 \gamma'^2,\qquad \dot{\gamma}'_{\rm sto} = 4D_0\gamma'^2,$$

while the shock term is written as

$$\dot{\gamma}'_{\rm sh}=A_0\gamma',\qquad t'_{\rm esc} = \frac{t'_{\rm acc}}{\eta_{\rm esc}}.$$

This formulation is explicitly designed to connect spectral curvature and interband lags (Hu et al., 2024).

A supernova-remnant realization uses a thin acceleration region immediately behind the forward shock and a larger downstream shell interior. In Tycho, Zone 1 is the thin rim with enhanced magnetic turbulence; Zone 2 is the broader, lower-$\gamma$0 shell interior into which electrons from Zone 1 advect and diffuse. This separation is essential because synchrotron emissivity scales strongly with $\gamma$1, whereas the soft photon field is nearly the same in both zones, allowing Zone 2 to remain X-ray dim but $\gamma$2-ray efficient (Atoyan et al., 2011).

3. Blazar applications: sequence models and time-dependent HBLs

A self-consistent two-zone leptonic model for blazars links the zone structure to the blazar sequence by varying essentially one physical control parameter, the mass accretion rate $\gamma$3. As $\gamma$4 increases, the disk luminosity rises, the wind density rises, the scattered external photon energy density $\gamma$5 rises strongly, and the phenomenological acceleration timescale increases. The resulting trend is stronger Compton dominance, softer $\gamma$6-ray spectra, lower synchrotron peak frequencies at high $\gamma$7, and high-energy output dominated by external Compton. In this framework, Flat Spectrum Radio Quasars have strong external Compton emission from the extended zone, whereas BL Lac objects are dominated by synchrotron and synchrotron self-Compton emission from the compact acceleration region (Boula et al., 16 Dec 2025, Boula et al., 2021).

The illustrative fit to 3C273 makes the zoning explicit. The acceleration zone is placed at $\gamma$8 with $\gamma$9 and $<2$0, while the cooling zone has $<2$1 and $<2$2. The external field geometry is taken as $<2$3 and $<2$4, with $<2$5, $<2$6, and $<2$7. In the paper’s interpretation, SSC primarily explains the X-ray emission from the acceleration zone, while external Compton from the cooling zone dominates the $<2$8-ray output (Boula et al., 16 Dec 2025).

The time-dependent Mrk 421 model addresses a different phenomenology. The source exhibited various intensity states differing by close to an order of magnitude in flux, with the fractional variability amplitude increasing with energy through the X-ray band, and Bayesian power spectral density analysis gave colored-noise indices ranging from $<2$9 to $B \gtrsim 400\,\mu{\rm G}$0. The two-zone model reproduces both X-ray spectra and time lags by combining a background zone with a rapidly varying compact zone. For the 10 XMM-Newton observations, the compact-zone parameters are approximately $B \gtrsim 400\,\mu{\rm G}$1–53, average $B \gtrsim 400\,\mu{\rm G}$2, $B \gtrsim 400\,\mu{\rm G}$3–0.23 G, average $B \gtrsim 400\,\mu{\rm G}$4 G, and $B \gtrsim 400\,\mu{\rm G}$5 cm, average $B \gtrsim 400\,\mu{\rm G}$6 cm. The lag interpretation uses

$B \gtrsim 400\,\mu{\rm G}$7

for synchrotron-cooling soft lags and

$B \gtrsim 400\,\mu{\rm G}$8

for acceleration-driven hard lags. The fitted values imply that the active zone is much smaller than the background zone and that the jet is strongly particle dominated in the active zone (Hu et al., 2024).

4. Hard-spectrum BL Lacs and the boundary with hybrid models

The extreme-BL-Lac literature is central to the contemporary discussion of two-zone modeling, but it also marks an important conceptual boundary. In these works, the outer zone is explicitly leptonic, whereas the inner zone is hadronic or lepto-hadronic. The geometry is therefore genuinely two-zone, but not strictly a two-zone leptonic model in the narrow sense. The outer blob is placed outside the BLR but inside the dusty torus and produces synchrotron, SSC, and external inverse Compton emission; the inner blob, located very close to the SMBH, produces the hardest TeV component through photo-hadronic interactions with a pair-plasma annihilation line at $B \gtrsim 400\,\mu{\rm G}$9 (Aguilar-Ruiz et al., 2022, Aguilar-Ruiz et al., 2021).

The leptonic outer zone is nevertheless a key component. In 1ES 0229+200, the outer blob explains the X-ray hump and the sub-TeV $R_{II} \simeq 10\,R_I.$0-ray emission. The X-ray synchrotron peak is taken at about

$R_{II} \simeq 10\,R_I.$1

which constrains the electron break Lorentz factor to $R_{II} \simeq 10\,R_I.$2 for the quoted normalization. The same electrons produce SSC emission and external inverse Compton emission on BLR/DT photons, but the BLR contribution is negligible for Compton scattering and the DT dominates the EIC channel. The model estimates

$R_{II} \simeq 10\,R_I.$3

and for electrons with $R_{II} \simeq 10\,R_I.$4 the scattering remains in the Thomson regime. The corresponding outer-zone EIC peak is

$R_{II} \simeq 10\,R_I.$5

This places the leptonic EIC component just below the $R_{II} \simeq 10\,R_I.$6-ray peak from the inner blob, so the combined spectrum appears hard (Aguilar-Ruiz et al., 2021).

The 2022 extension to six extreme BL Lacs emphasizes the same division of labor. Typical outer-blob parameters are $R_{II} \simeq 10\,R_I.$7, $R_{II} \simeq 10\,R_I.$8, $R_{II} \simeq 10\,R_I.$9, and $B \propto \frac{1}{z},$0, with $B \propto \frac{1}{z},$1 in the outer blob. The paper’s point is that the outer blob can fit the data with much more typical conditions than a one-zone SSC model, while the inner blob supplies the hardest TeV component. Without the pair-plasma or inner blob, the source can still look EHSP-like, but the TBL/hard-TeV behavior is lost. A common misconception is therefore that every “two-zone leptonic” hard-TeV model is fully leptonic; in the extreme-BL-Lac papers most often cited, the decisive hard-TeV tail is not (Aguilar-Ruiz et al., 2022).

5. Supernova-remnant realization

The Tycho model is the clearest supernova-remnant example of a two-zone leptonic source. Zone 1 is the thin, high-$B \propto \frac{1}{z},$2 acceleration region at the forward shock, with rim width inferred from the X-ray filament width and volume filling factor

$B \propto \frac{1}{z},$3

Zone 2 is the much larger downstream shell interior, with

$B \propto \frac{1}{z},$4

The numerical example adopts

$B \propto \frac{1}{z},$5

with electron injection only in Zone 1,

$B \propto \frac{1}{z},$6

and no significant in-situ electron acceleration in Zone 2, $B \propto \frac{1}{z},$7 (Atoyan et al., 2011).

This geometry redistributes the emission by band. In the radio/sub-mm, the emission is dominated by Zone 2 because it occupies much more volume. In non-thermal X-rays, Zone 1 dominates because the higher magnetic field boosts synchrotron emissivity. The gamma-ray emission below $B \propto \frac{1}{z},$8 GeV is dominated by bremsstrahlung, mostly from Zone 2, with adopted density $B \propto \frac{1}{z},$9. At TeV energies the dominant leptonic process is inverse Compton scattering, again mainly from Zone 2, using the CMB and FIR photons from thermal dust in Tycho’s NE and NW shells; the FIR field can contribute up to about 30% of the total Compton flux at TeV energies. The reported total electron energies are

$Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$0

The model therefore shows how a purely leptonic two-zone interpretation can reproduce the full Fermi + VERITAS spectrum (Atoyan et al., 2011).

At the same time, the Tycho study is explicit that a hadronic model can also fit the radiation spectrum. The key discriminator is the low-energy $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$1-ray behavior below about $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$2 MeV, especially the $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$3-decay cutoff or turnover in the range $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$4. This is a recurring feature of two-zone modeling more generally: the two-zone structure may be strongly motivated, yet the exact leptonic-versus-hadronic content of one zone can remain degenerate.

6. Observational diagnostics, polarization, and limitations

A generalized multi-zone leptonic framework has sharpened the observational consequences of two-zone thinking by treating the source as $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$5 incoherent patches with uncorrelated magnetic fields, with

$Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$6

In that picture, the two-zone case is only a special case of a broader patchy-source model. For synchrotron emission, the standard random-walk suppression is the primary depolarization factor. For SSC, the paper identifies an additional secondary depolarization because each zone’s SSC output is itself a sum over seed photons arriving from many zones with different polarization orientations and internal light-travel delays. The resulting leptonic polarization degree is expected to be much lower than the one-zone prediction, with typical values around $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$7 for realistic patchiness. Detectability with IXPE within $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$8 ksec requires intrinsic optical synchrotron polarization $Q_{\rm inj} = \frac{n_{e,I}}{t_{\rm esc,I}}.$9 and a stable polarization angle during the exposure (Zhang et al., 2024).

This has two immediate implications. First, low high-energy polarization is a natural multi-zone SSC signature rather than a failure of the leptonic interpretation. Second, timing and polarization need not track each other one-to-one, because internal light-crossing delays and energy stratification decorrelate the seed synchrotron field from the observed SSC output. A plausible implication is that two-zone leptonic models are best tested through joint use of SED fitting, lag spectra, variability amplitude as a function of energy, and polarization limits, rather than through any single observable.

The limitations are equally explicit in the cited literature. The model is not unique: densities, filling factors, magnetic fields, and transport parameters are not uniquely constrained in the Tycho application; a more realistic remnant may contain additional zones; and in blazars, generalized multi-zone radiative transfer may be more appropriate than a literal two-zone toy model (Atoyan et al., 2011, Zhang et al., 2024). The term “two-zone leptonic model” therefore names a productive modeling strategy rather than a single universal construction. Its common content is the use of two radiatively distinct leptonic environments to separate acceleration from cooling, background from flare, or X-ray efficiency from $$U_{\rm sc} \simeq \frac{L_{\rm disk}\,\tau_T}{4\pi R_2^2 c},\qquad \tau_T = n_0 \sigma_T R_s \ln\!\left(\frac{R_2}{R_1}\right),\qquad U_{\rm ext} = \Gamma^2 U_{\rm sc}.$$0-ray efficiency, thereby avoiding the stronger assumptions of one-zone homogeneity.