Synchrotron-Dominated Cascades Explained

• Synchrotron-dominated cascades are electromagnetic regimes where secondary e± lose energy primarily via synchrotron radiation in strong magnetic fields, producing distinctive photon spectra.
• They occur in varied environments—from UHECR-seeded blazar halos and compact gamma-ray sources to magnetar magnetospheres—each showing unique spectral slopes and temporal behaviors.
• Numerical models hybridize Monte Carlo sampling with kinetic transport equations to simulate rapid pair injection and cooling, clarifying cascade closure and observable diagnostics.

Searching arXiv for relevant papers on synchrotron-dominated cascades and related cascade formalisms. Synchrotron-dominated cascades are electromagnetic cascade regimes in which the secondary $e^\pm$ produced by high-energy photons or hadrons lose energy primarily through synchrotron radiation rather than inverse Compton scattering. Across blazar environments, compact $\gamma$-ray–opaque sources, intense magnetic fields, and magnetar magnetospheres, this regime arises when magnetic energy density and kinematic conditions force rapid synchrotron cooling of cascade pairs, yielding distinctive photon spectra, angular signatures, temporal behavior, and, in some cases, strong polarization. In the literature summarized here, the term encompasses several closely related but physically distinct settings: ultra-high-energy cosmic-ray–induced cascades in magnetised large-scale structure around blazars (Oikonomou et al., 2014), time-dependent synchrotron-dominated electromagnetic cascades in relativistic jets (Asano et al., 2020), linear internal cascades in compact opaque sources (Fiorillo et al., 29 Aug 2025), one-photon pair cascades in strong magnetic fields (Wang et al., 2018), and pair-synchrotron cascades in ultra-strong magnetar fields (Harding et al., 15 Aug 2025).

1. Defining physical regime

The common structural element is a feedback loop between pair creation and synchrotron emission. A high-energy primary photon, proton, or electron injects secondary $e^\pm$; these pairs radiate synchrotron photons; and those photons may in turn generate additional pairs if the local opacity permits it. The cascade is therefore controlled by the competition among pair injection, radiative cooling, escape, and further conversion.

In the blazar-cosmic-web scenario, ultra-high-energy cosmic rays (UHECRs; primarily protons) accelerated in a blazar jet escape into a magnetised region of size $R\sim$ a few Mpc and field strength $B\sim10^{-7}$ G, where they interact with CMB and EBL photons through Bethe–Heitler pair production and photomeson production (Oikonomou et al., 2014). The immediately produced secondary electrons and positrons of energy $E\gg m_e c^2$ then lose energy almost exclusively by synchrotron radiation in the local field. This establishes a cascade whose observable output is a synchrotron pair halo or echo.

In compact-source internal-cascade theory, the defining assumption is a $\gamma$-ray–opaque region of radius $R$ containing a dense target photon field $n_t(\epsilon_t)$ and a magnetic field $B$ such that electron losses are dominated by synchrotron rather than inverse Compton scattering. The source is “linear,” meaning cascade particles interact only with the fixed background $\gamma$0, and one requires $\gamma$1 so that the cascade is fully developed (Fiorillo et al., 29 Aug 2025).

In intense magnetic fields, the regime changes from ordinary $\gamma$2-driven cascades to one-photon pair production plus quantum synchrotron emission. A high-energy photon of energy $\gamma$3 propagating across a magnetic field $\gamma$4 can be absorbed if

$\gamma$5

with $\gamma$6 G, producing an $\gamma$7 pair whose members radiate high-energy photons via synchrotron emission and continue the cascade (Wang et al., 2018).

In magnetar magnetospheres, the cascade is initiated by resonant inverse Compton scattering (RICS) of thermal X-rays by relativistic particles on closed magnetic loops. The resulting photons are attenuated by one-photon pair production and photon splitting, with the synchrotron emission of the secondary pairs dominating the emergent sub-MeV spectrum for most observer angles (Harding et al., 15 Aug 2025).

A general regime criterion appears explicitly in compact-source treatments: synchrotron dominates if

$\gamma$8

and, in jet applications, Klein–Nishina suppression of inverse Compton cooling can strengthen synchrotron dominance at high $\gamma$9 (Asano et al., 2020, Fiorillo et al., 29 Aug 2025).

2. Transport equations and cascade closure

The electron distribution in synchrotron-dominated cascades is usually described by kinetic transport equations. For UHECR-induced blazar cascades, the differential number density $e^\pm$0 obeys

$e^\pm$1

where $e^\pm$2 is the total energy-loss rate, $e^\pm$3 is the $e^\pm$4 injection rate from $e^\pm$5, and $e^\pm$6 is the escape timescale or diffusion time. In steady state and for slow escape,

$e^\pm$7

This form makes explicit that the steady lepton population is set by a balance between continuous synchrotron cooling and injection (Oikonomou et al., 2014).

The corresponding synchrotron loss rate is

$e^\pm$8

and the characteristic synchrotron photon energy is

$e^\pm$9

For first-generation $R\sim$0 from photomeson production with $R\sim$1 eV in $R\sim$2 G, one finds $R\sim$3 eV, i.e. in the TeV band (Oikonomou et al., 2014).

In time-dependent jet models, the evolution is written as a Fokker–Planck equation including stochastic acceleration: $R\sim$4 Here $R\sim$5 is the energy-diffusion coefficient, $R\sim$6 the primary injection, and $R\sim$7 the secondary-pair source term from photon-photon annihilation. For the stochastic-acceleration “hard-sphere” case,

$R\sim$8

so reacceleration applies equally to primary and secondary pairs (Asano et al., 2020).

In the linear compact-source theory, the steady-state electron and photon equations are given, in the $R\sim$9-function approximation, by

$B\sim10^{-7}$0

and

$B\sim10^{-7}$1

These equations codify, energy bin by energy bin, the balance between synchrotron cooling and pair injection on the lepton side, and between synchrotron photon production and pair-production losses on the photon side (Fiorillo et al., 29 Aug 2025).

This suggests that “cascade closure” can be formulated in a unified way across environments: one needs sufficiently rapid pair creation to replenish the radiating lepton pool and sufficiently rapid synchrotron cooling to transfer energy downward before escape or competing channels remove it.

3. Numerical methodologies and approximations

The numerical treatment depends strongly on the astrophysical setting. In UHECR-seeded blazar cascades, propagation is handled with CRPropa 2 in 3D mode. Interaction lengths for $B\sim10^{-7}$2 processes on the CMB, EBL, and radio backgrounds are sampled using SOPHIA for photo-pion production and the Kelner & Aharonian parametrisation for Bethe–Heitler pair production. Deflections in the structured magnetic field, with coherence $B\sim10^{-7}$3 Mpc, are tracked analytically. Every time a secondary photon or $B\sim10^{-7}$4 is produced, its position and momentum are recorded, and particles pointing within the observer cone ($B\sim10^{-7}$5) are collected. The collected leptons are then propagated by solving the kinetic equation on the fly, including continuous synchrotron losses, inverse Compton scattering when relevant, and further pair production above threshold; the code of Murase & Dermer implements the fully coupled cascade (Oikonomou et al., 2014).

In the 3C 279 flare model, the simulations are explicitly time-dependent and one-zone. The synchrotron-dominated electromagnetic cascade appears in a stochastic-acceleration scenario in which secondary pairs from $B\sim10^{-7}$6 annihilation are reaccelerated in the same emitting region. Klein–Nishina suppression is incorporated in the inverse Compton loss term through either $B\sim10^{-7}$7 or a suppression factor

$B\sim10^{-7}$8

so that

$B\sim10^{-7}$9

This treatment is essential because the synchrotron-dominated state is reached precisely where inverse Compton losses are suppressed at high Lorentz factor (Asano et al., 2020).

In magnetar calculations, Monte Carlo methods propagate “weighted” photons emitted from RICS spectral bins in the tangent direction to the local field. Along each photon path, opacities to one-photon pair production and photon splitting are accumulated, and a rejection-sampling scheme determines whether the photon splits, pair produces, or escapes. Secondary photons, whether split or synchrotron, are followed recursively until all weight escapes. For synchrotron emission from pairs created in excited Landau states, the code uses the classical form

$E\gg m_e c^2$0

in the “circular” frame when appropriate, but switches to full QED rates in supercritical fields $E\gg m_e c^2$1 (Harding et al., 15 Aug 2025).

In strong-field one-dimensional magnetic cascades, a dedicated Monte Carlo code is used to study the spectrum after the cascade reaches saturated status, when almost all the energy of the primary particles transfers to photons. The resulting spectral energy distribution is then fit analytically as a function only of the product $E\gg m_e c^2$2 (Wang et al., 2018).

A recurring methodological theme is the hybridization of event-based Monte Carlo interaction sampling with deterministic or semi-deterministic kinetic evolution of secondaries. This is especially pronounced in the blazar-halo framework, where proton interactions are stochastic while lepton cooling is treated through transport equations (Oikonomou et al., 2014).

4. Spectral structure and universal behaviors

Several distinct universal or quasi-universal spectral patterns appear in synchrotron-dominated cascades, but they are not identical across regimes.

For linear internal cascades in compact opaque sources, a key analytic result is the emergence of two photon power-law regimes when the injection is effectively monochromatic at $E\gg m_e c^2$3. At high photon energies $E\gg m_e c^2$4,

$E\gg m_e c^2$5

where $E\gg m_e c^2$6 and $E\gg m_e c^2$7. At lower photon energies $E\gg m_e c^2$8, in the cooling-only regime,

$E\gg m_e c^2$9

The approximate cascade spectrum can therefore be written as

$\gamma$0

The reported simulations confirm the analytic indices $\gamma$1 above $\gamma$2 and $\gamma$3 below it in AGN coronae, GRB internal shocks, some compact blazar models, and TDE environments (Fiorillo et al., 29 Aug 2025).

By contrast, saturated cascades in intense magnetic fields exhibit a different universal broken-power-law structure. Defining the normalized SED

$\gamma$4

the shape below the cutoff depends only on $\gamma$5 and is fit by

$\gamma$6

with

$\gamma$7

The break is at $\gamma$8. Below $\gamma$9, one has $R$0; between $R$1 and $R$2, $R$3; and above $R$4 the spectrum displays a very sharp cutoff, much steeper than a simple exponential (Wang et al., 2018).

In UHECR-driven structured-region cascades, the spectrum is described as a broad bump peaking at $R$5–10 TeV with the high-energy cutoff set by EBL absorption. Below the peak, the slope is governed by the $R$6 injection spectrum and synchrotron cooling, with $R$7 and $R$8–2 (Oikonomou et al., 2014). In the 3C 279 stochastic-cascade model, the resulting electron distribution is described as very hard and nearly flat, $R$9, feeding back into a flat $n_t(\epsilon_t)$0 continuum from X-rays to GeV $n_t(\epsilon_t)$1 rays, with a cooling break at $n_t(\epsilon_t)$2 MeV and a sharp cutoff above $n_t(\epsilon_t)$3 GeV from internal $n_t(\epsilon_t)$4 absorption and EBL attenuation (Asano et al., 2020).

These results caution against treating “synchrotron-dominated cascade spectrum” as a single universal template. The asymptotic form depends on whether the cascade is linear or nonlinear, whether it is driven by $n_t(\epsilon_t)$5 absorption or one-photon magnetic conversion, whether the synchrotron process is classical or QED, and whether the photon field is external, internal, or negligible.

5. Source classes and phenomenology

The best-developed astrophysical application in the provided material concerns extreme TeV blazars. For 1ES 0229+200, RGB J0710+591, and 1ES 1218+304, synchrotron emission of UHECR secondaries produced in blazars embedded in structured regions with magnetic field strengths of order $n_t(\epsilon_t)$6 G is presented as an alternative to both standard leptonic scenarios and inverse-Compton cascade channels (Oikonomou et al., 2014). The observed specific flux at Earth is written as

$n_t(\epsilon_t)$7

with $n_t(\epsilon_t)$8 and an obvious generalization to radially varying $n_t(\epsilon_t)$9. The total synchrotron luminosity in the halo is approximated as

$B$0

where $B$1 is the fraction of proton energy converted through $B$2 interactions within the few-Mpc region (Oikonomou et al., 2014).

The angular extent and time-delay phenomenology are specific and potentially diagnostic. The parent pairs are deflected over their cooling length by an angle

$B$3

so the halo is $B$4 for sources beyond $B$5 Mpc, implying a point-like image for current IACTs but a resolvable source for CTA at TeV energies. The delay of the synchrotron echo is

$B$6

which yields month-to-year echoes (Oikonomou et al., 2014).

In the 3C 279 June 2015 flare, the synchrotron-dominated electromagnetic cascade appears in a distinct jet context. The stochastic-acceleration solution uses $B$7, $B$8 cm, $B$9 G, $\gamma$00, primary monoenergetic injection at $\gamma$01, and one-sided jet luminosities $\gamma$02 erg s$\gamma$03, $\gamma$04 erg s$\gamma$05, and $\gamma$06 erg s$\gamma$07 (Asano et al., 2020). The observed maximum synchrotron photon energy follows from

$\gamma$08

and

$\gamma$09

with $\gamma$10 in Model C, implying $\gamma$11 GeV (Asano et al., 2020).

In magnetars, the relevant observables shift to hard X-rays and polarization. For a single loop with maximum radius $\gamma$12, surface field $\gamma$13–$\gamma$14, and $\gamma$15–$\gamma$16, the secondary synchrotron component dominates at 5–200 keV when RICS losses are included, with spectral index $\gamma$17 and cutoff at $\gamma$18–$\gamma$19 MeV, while split photons fill the MeV window up to $\gamma$20 MeV. The synchrotron photons are predicted to yield net polarization 40%–80% in $\gamma$21 mode below the SR cutoff, depending on viewing angle, with the possibility of a mode reversal when split photons dominate at higher energies (Harding et al., 15 Aug 2025).

6. Parameter dependence, robustness, and comparisons with inverse-Compton cascades

A central result of the blazar-halo scenario is that the generated synchrotron pair halo or echo flux at the peak energy is not sensitive to variations in the overall IGMF strength. The reason given is geometrical and environmental: filaments occupy only $\gamma$22 of void volume, so cascade photons produced far outside the structured region suffer large deflections in void IGMFs and are isotropised, suppressing their contribution to the beamed signal. The TeV synchrotron emission is therefore determined almost entirely by the few-Mpc, $\gamma$23 G region around the source and is insensitive to weaker fields in voids (Oikonomou et al., 2014).

The same work reports several robustness tests. Varying the filament magnetic field between $\gamma$24 and $\gamma$25 G shifts the synchrotron bump energy only within the observed TeV range. The absolute flux at peak depends mainly on the isotropic UHECR luminosity, $\gamma$26 erg s$\gamma$27, but not on the void field strength. Changing the UHECR injection index from $\gamma$28 to $\gamma$29 requires only modest $\gamma$30 factor-few adjustments in $\gamma$31 to fit the same data. Scans in maximum energy from $\gamma$32 to $\gamma$33 eV show that $\gamma$34 eV is required to reproduce the highest-energy tail above $\gamma$35 TeV after EBL attenuation, and the fit is only weakly sensitive to the choice among current EBL models (Oikonomou et al., 2014).

In the 3C 279 jet model, the defining comparison is instead with prompt-injection or reconnection-based synchrotron scenarios. A too strong magnetic field yields a too bright synchrotron X-ray flux due to secondary electron–positron pairs, whereas the stochastic-cascade scenario operates in a very low-magnetization regime, with $\gamma$36 at $\gamma$37, and uses KN suppression of IC to force synchrotron-dominated cooling of the highest-energy particles (Asano et al., 2020). This indicates that synchrotron-dominated does not necessarily imply magnetically dominated in the dynamical sense; it may instead arise from radiative suppression of the IC channel.

In generalized linear-cascade theory, universality is contingent rather than automatic. Under sufficiently high $\gamma$38 and $\gamma$39, a synchrotron cascade with $\gamma$40 can emerge in the optical–X-ray band in compact blazar environments. However, in many blazar models Bethe–Heitler pair injection over a broad range or partial $\gamma$41 attenuation spoils the universality (Fiorillo et al., 29 Aug 2025). A plausible implication is that observationally inferred departures from the idealized $\gamma$42 and $\gamma$43 indices need not falsify the synchrotron-cascade mechanism itself; they may instead indicate non-monochromatic injection, incomplete opacity, mixed cooling, or additional pair-production channels.

The comparison with inverse-Compton cascades is recurrent. In voids, inverse-Compton pair echoes are highly sensitive to the poorly known intergalactic magnetic field: for $\gamma$44 G delays are $\gamma$45 months, but for $\gamma$46 G delays become $\gamma$47 years and flux is suppressed. By contrast, the structured-region synchrotron signal is described as unavoidable and robust against void-field uncertainty (Oikonomou et al., 2014). In compact-source theory, inverse Compton losses must be subdominant throughout the cascade for the synchrotron-universal spectrum to appear (Fiorillo et al., 29 Aug 2025). In strong magnetic fields and magnetars, inverse Compton is not the competing closure process at all; instead, the competition is with magnetic pair creation and photon splitting (Wang et al., 2018, Harding et al., 15 Aug 2025).

7. Conceptual distinctions and common misconceptions

One common source of confusion is the use of the same label for physically different cascade architectures. In extragalactic blazar environments, the cascade may be seeded by UHECR interactions on the CMB and EBL and radiate in $\gamma$48 G structured regions, producing TeV synchrotron halos or echoes (Oikonomou et al., 2014). In compact jets, the cascade can instead be an internal $\gamma$49-driven process with reacceleration of secondary pairs in a one-zone emission region (Asano et al., 2020). In intense magnetic fields, the cascade is driven by one-photon pair production and quantum synchrotron emission rather than by $\gamma$50 absorption (Wang et al., 2018). In magnetars, photon splitting is a major participant in the cascade and shapes both spectrum and polarization (Harding et al., 15 Aug 2025). These are not merely different parameter choices of one theory.

A second misconception is that synchrotron-dominated cascades necessarily produce the same spectral slope. The material summarized here shows at least three distinct asymptotic families: the $\gamma$51 and $\gamma$52 branches of linear internal cascades (Fiorillo et al., 29 Aug 2025); the broken SED with indices $\gamma$53 and $\gamma$54 in the saturated strong-field QED cascade, expressed in $\gamma$55 rather than directly in $\gamma$56 (Wang et al., 2018); and the broad TeV bump or flat $\gamma$57 continuum characteristic of specific source applications (Oikonomou et al., 2014, Asano et al., 2020). The distinction reflects differences in kernel, opacity, geometry, and observable variable.

A third misconception is that synchrotron dominance requires extreme magnetic field strength in an absolute sense. The blazar-halo framework achieves TeV synchrotron emission at only $\gamma$58 G because the pairs are ultra-relativistic, with $\gamma$59 eV (Oikonomou et al., 2014). The 3C 279 stochastic-cascade model operates with $\gamma$60 G, but KN suppression renders inverse Compton inefficient at the relevant energies (Asano et al., 2020). By contrast, the strong-field and magnetar cases truly involve QED-scale fields where one-photon pair creation and photon splitting become central (Wang et al., 2018, Harding et al., 15 Aug 2025).

Finally, the observational appearance of synchrotron-dominated cascades is not uniformly extended or slowly varying. In structured extragalactic regions, the halo can be $\gamma$61 and point-like for current IACTs, with month-to-year echoes (Oikonomou et al., 2014). In the 3C 279 flare model, the same broad physical category accommodates minute-scale variability, with $\gamma$62 min (Asano et al., 2020). In magnetars, the distinguishing diagnostics are instead hard-X-ray tails and polarization fractions of 40%–80% in the keV band (Harding et al., 15 Aug 2025).

Taken together, these results establish synchrotron-dominated cascades as a broad class of radiation-mediated pair cascades whose detailed manifestation is controlled by magnetic geometry, opacity, injection channel, and radiative hierarchy. The unifying principle is rapid synchrotron cooling of secondary pairs; the diversity lies in how those pairs are generated, how the cascade closes, and which observables carry the imprint.