Proton Synchrotron Emission Physics

• Proton synchrotron emission is a mechanism where ultra-relativistic protons emit radiation in strong magnetic fields, requiring extreme conditions like high Lorentz factors and intense fields.
• The process is invoked to explain high-energy components in blazars, GRBs, and AGN jets, with models constrained by cooling times, energy budgets, and compact emission regions.
• Recent research integrates proton synchrotron emission into multimessenger frameworks, evaluating its role in complex astrophysical environments and its compatibility with observational data.

Searching arXiv for recent and foundational papers on proton synchrotron emission relevant to the provided corpus. Searching arXiv for "proton synchrotron emission blazar GRB jet". Proton synchrotron emission is synchrotron radiation produced by ultra-relativistic protons in magnetic fields. It has been invoked as a direct radiative channel for the high-energy component of blazars, gamma-ray bursts, large-scale AGN jets, compact radio lobes, and magnetized reconnection layers, while in other contexts proton populations influence synchrotron emission indirectly through secondary leptons rather than through direct proton radiation. Because the proton mass strongly suppresses synchrotron losses relative to electrons, the mechanism generally requires very high proton Lorentz factors, strong magnetic fields, compact emitting regions, or some combination thereof, and much of the literature is organized around whether those requirements are physically acceptable in a given source class (Xue et al., 2023, Ghisellini et al., 2019).

1. Physical basis and characteristic scales

In the standard hadronic interpretation of a blazar or GRB spectrum, the low-energy hump is usually attributed to electron synchrotron radiation, whereas the high-energy hump is assigned to synchrotron radiation from ultra-relativistic protons in the comoving magnetic field. A common one-zone setup assumes a compact region of radius $R$, bulk Lorentz factor $\Gamma$, Doppler factor $\delta_{\rm D}\approx \Gamma$, and uniform magnetic field $B$, with proton synchrotron required to reproduce an observed peak energy $E_{\rm peak}^{\rm obs}$ and luminosity $(\nu L_\nu)^{\rm obs}_{\rm peak}$ (Xue et al., 2023).

A useful monochromatic approximation relates the proton Lorentz factor to the observed peak energy,

$$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$

while the proton synchrotron cooling efficiency is written as

$$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$

These expressions make explicit why proton synchrotron emission is difficult to realize: efficient radiation requires large $B$, large $\gamma_{\rm p,max}$, or both, because of the $\Gamma$0 suppression in the cooling efficiency (Xue et al., 2023).

The same scaling appears in prompt-GRB studies. For a fixed observed synchrotron frequency, the observed proton synchrotron cooling time is longer than the electron case by

$\Gamma$1

which is the central reason proton synchrotron has been proposed for incomplete-cooling prompt spectra (Ghisellini et al., 2019). In TXS 0506+056, the jet-frame proton energy required to produce synchrotron photons at observed critical frequency $\Gamma$2 is written as

$\Gamma$3

and the observer-frame cooling time as

$\Gamma$4

showing directly how long-lived VHE proton synchrotron activity can emerge when the magnetic field is of order a few gauss (Sunanda et al., 2022).

Analytical treatments usually supplement these relations with a Hillas-type confinement condition and a transparency condition against internal $\Gamma$5 absorption. In one-zone blazar applications, the allowed parameter space in the $\Gamma$6 plane is therefore carved out simultaneously by the Eddington constraint, the acceleration requirement, and the condition $\Gamma$7 (Xue et al., 2023).

2. Energetics in blazars

The most systematic energetic critique of proton synchrotron models in blazars is the study of steady $\Gamma$8-ray emission in a sample of 145 sources. In that framework the absolute power of a two-sided jet is written as

$\Gamma$9

and the quantity $\delta_{\rm D}\approx \Gamma$0 is defined as the minimum total jet power allowed by the proton synchrotron model after minimizing over the unknown magnetic field for fixed observables and jet kinematics. Because the calculation assumes monoenergetic particles, high proton radiative efficiency, and adopted Doppler factors, the resulting $\delta_{\rm D}\approx \Gamma$1 is already a conservative lower limit (Liodakis et al., 2020).

For the 145-source sample, the observables were Monte Carlo sampled with $\delta_{\rm D}\approx \Gamma$2 realizations per source, using uncertainties of $\delta_{\rm D}\approx \Gamma$3 dex in luminosities, $\delta_{\rm D}\approx \Gamma$4 dex in peak frequencies, and a variability-timescale distribution with

$\delta_{\rm D}\approx \Gamma$5

The central result is the so-called energy crisis: for most sources,

$\delta_{\rm D}\approx \Gamma$6

and the same lower bound is typically about two orders of magnitude above $\delta_{\rm D}\approx \Gamma$7 and the optimistic Blandford–Znajek estimate $\delta_{\rm D}\approx \Gamma$8. The derived magnetic fields imply either local amplification by a factor of about $\delta_{\rm D}\approx \Gamma$9 or a $B$0-ray production site at sub-parsec scales; the former corresponds to a median amplification factor $B$1, while the latter yields a median $B$2 pc and places the emission well inside the BLR. The predicted neutrino emission peaks at $B$3 EeV, with typical muon-neutrino peak fluxes $B$4 of the peak $B$5-ray flux (Liodakis et al., 2020).

A more selective view emerges when the proton synchrotron parameter space is mapped directly against peak energy and luminosity. In that analysis, proton synchrotron can fit the high-energy hump when it peaks beyond tens of GeV without violating basic observations and theories, especially for $B$6–$B$7, $B$8 G, and $B$9. For humps peaking in the $E_{\rm peak}^{\rm obs}$0–$E_{\rm peak}^{\rm obs}$1 GeV range, however, the model typically requires super-Eddington jet power and, if $E_{\rm peak}^{\rm obs}$2, magnetic fields that can exceed $E_{\rm peak}^{\rm obs}$3 G. For humps peaking in the $E_{\rm peak}^{\rm obs}$4–$E_{\rm peak}^{\rm obs}$5 keV band, the outcome depends strongly on luminosity: Mrk 421 admits an allowed region, whereas 3C 279 does not (Xue et al., 2023).

Taken together, these blazar results separate steady luminous GeV-peaked sources from high-peaking VHE sources. The former are described as energetically problematic, whereas the latter remain viable in larger emission regions and at moderate luminosity (Liodakis et al., 2020, Xue et al., 2023).

3. Gamma-ray bursts and other transients

Prompt-GRB applications were motivated by spectra that exhibit

$E_{\rm peak}^{\rm obs}$6

below a break and

$E_{\rm peak}^{\rm obs}$7

up to the peak, a shape interpreted as synchrotron emission from particles that are in fast cooling but do not cool completely during the dynamical time. Proton synchrotron was proposed because it can preserve compact emission regions and strong magnetic fields while lengthening the cooling time relative to the electron case. For prompt emission observed at $E_{\rm peak}^{\rm obs}$8 keV with a cooling time of order $E_{\rm peak}^{\rm obs}$9 s, the proton-synchrotron interpretation was argued to work with $(\nu L_\nu)^{\rm obs}_{\rm peak}$0 G and $(\nu L_\nu)^{\rm obs}_{\rm peak}$1 cm, whereas the electron-synchrotron alternative would require $(\nu L_\nu)^{\rm obs}_{\rm peak}$2 G and $(\nu L_\nu)^{\rm obs}_{\rm peak}$3 cm, in tension with prompt variability (Ghisellini et al., 2019).

Detailed semi-analytical and numerical calculations, however, make the prompt-GRB proton-synchrotron hypothesis highly constrained. In the marginally fast-cooling formulation, typical magnetic fields are

$(\nu L_\nu)^{\rm obs}_{\rm peak}$4

and for $(\nu L_\nu)^{\rm obs}_{\rm peak}$5 the median jet Poynting luminosity is

$(\nu L_\nu)^{\rm obs}_{\rm peak}$6

with strong low-energy spectral contamination from secondary pairs unless $(\nu L_\nu)^{\rm obs}_{\rm peak}$7. On that basis the model is described as strongly disfavoured for the low-energy spectral breaks of prompt GRBs (Florou et al., 2021). A related study of Bethe–Heitler pair production finds two regimes: at high $(\nu L_\nu)^{\rm obs}_{\rm peak}$8, large radius, and low luminosity, proton synchrotron can dominate and may leave a subdominant pair-synchrotron power law extending to tens or hundreds of MeV; at low $(\nu L_\nu)^{\rm obs}_{\rm peak}$9, small radius, and high luminosity, Bethe–Heitler cooling drives the spectrum toward a single fast-cooling power law $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$0 across the entire GBM/Swift band, which is incompatible with observations (Bégué et al., 2021).

By contrast, afterglow and reverse-shock applications are considerably more favorable. For GRB 190114C, a two-component synchrotron model with electron synchrotron for the X-rays and proton synchrotron for the $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$1–$$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$2 TeV MAGIC emission reproduces the data with isotropic explosion energy $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$3 erg, ambient density $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$4–$$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$5, a few-percent accelerated fractions, and protons reaching a few $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$6 eV (Isravel et al., 2022). For GRB 221009A, reverse-shock proton synchrotron has been used to explain $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$7TeV photons, including the possibility of $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$8 TeV photons with reasonable EBL models, and a structured-jet analysis associates the rise of the reverse-shock proton-synchrotron component with the $$\gamma_{\rm p,max} = \left( \frac{E_{\rm peak}^{\rm obs}(1+z)} {1.53\times10^3\,h\,B\,\delta_{\rm D}} \right)^{1/2},$$9–$$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$0 s hardening interval and the arrival of the $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$1 TeV photon (Zhang et al., 2022, Zhang et al., 2023).

Transient blazar activity provides an intermediate case. In TXS 0506+056, the 2017 sequence has been modeled with electron synchrotron plus SSC for the HE flare and proton synchrotron for the delayed VHE activity. In that fit the source parameters are $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$2 G, $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$3, $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$4 cm, and a proton luminosity $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$5 erg s$$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$6, with proton synchrotron cooling times in the $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$7 to $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$8 day interval matching the observed $$f_{\rm p,syn}=\min\left\{\frac{t_{\rm dyn}}{t_{\rm p,syn}},\,1\right\} =\min\left\{ \frac{\sigma_{\rm T}B^2R\gamma_{\rm p,max}} {6\pi m_e c^2 (m_p/m_e)^3}, \,1 \right\}.$$9-day VHE episode (Sunanda et al., 2022).

4. Extended jets, knots, and compact radio lobes

Large-scale AGN jets provide a distinct proton-synchrotron environment because the emission region is spatially extended and only mildly beamed. In the large-scale jet of 3C 273, the X-ray and GeV emission from knot A were modeled with a broken power-law proton distribution after the IC/CMB interpretation was judged inconsistent with Fermi-LAT constraints. Two regimes were considered. In the cooling-dominated case, the fit uses $B$0 kpc, $B$1 mG, $B$2, $B$3, and requires a total luminosity of about $B$4. In the escape-dominated case, the parameters become $B$5 kpc and $B$6 mG, with powers $B$7 and $B$8, making the escape-dominated solution the more favorable one (Kundu et al., 2014).

A broader application to PKS 0637-752 and 3C 273 treated proton synchrotron as an alternative to IC/CMB for extended quasar jets. The characteristic requirements are milligauss magnetic fields and proton energies of order $B$9–$\gamma_{\rm p,max}$0, but the luminosity budgets are comparatively modest: for the combined knots of PKS 0637-752 the total power is about $\gamma_{\rm p,max}$1, or about $\gamma_{\rm p,max}$2 of Eddington, and for 3C 273 the required power is around $\gamma_{\rm p,max}$3, or about $\gamma_{\rm p,max}$4 of Eddington (Bhattacharyya et al., 2015). AP Librae represents the opposite extreme. There the extended-jet proton-synchrotron interpretation of the VHE emission requires $\gamma_{\rm p,max}$5 mG, $\gamma_{\rm p,max}$6 eV, and

$\gamma_{\rm p,max}$7

which is more than $\gamma_{\rm p,max}$8 times the Eddington luminosity of AP Librae, leading that scenario to be described as unlikely (Basumallick et al., 2017).

Compact radio lobes introduce another characteristic proton-synchrotron signature. In parsec-scale mini lobes with $\gamma_{\rm p,max}$9 pc and $\Gamma$00 G, the proton synchrotron characteristic photon energy is

$\Gamma$01

placing the direct proton-synchrotron bump in the sub-MeV band. Its visibility depends strongly on the primary electron synchrotron luminosity: for $\Gamma$02 the sub-MeV bump appears clearly, whereas for $\Gamma$03 leptonic emission hides the hadronic signature. In intermediate cases, synchrotron from secondary $\Gamma$04 produced in the photopion cascade can emerge in the GeV–TeV range, yielding a double-bump hadronic spectrum (Kino et al., 2010).

5. Strong magnetic fields and indirect hadronic synchrotron signatures

In sufficiently strong magnetic fields, proton synchrotron emission becomes a quantum problem rather than a semiclassical one. The quantum-field-theoretic treatment of pion production in a uniform magnetic field starts from the Dirac equation with the proton anomalous magnetic moment,

$\Gamma$05

and computes the emission rate from the imaginary part of the proton self-energy. In this framework the fully quantum pion-synchrotron rate is much smaller than semiclassical estimates unless the anomalous magnetic moment is included; with the anomalous magnetic moment, the pion synchrotron decay width can be enhanced by about two orders of magnitude, and in favorable spin-flip channels the enhancement can approach $\Gamma$06 (Maruyama et al., 2015).

A later relativistic quantum treatment extends this framework to photon, pion, and $\Gamma$07-meson production and derives a scaling rule for transitions between very large Landau levels. There the proton curvature parameter is

$\Gamma$08

the transition probabilities can be extrapolated to $\Gamma$09, recoil is included explicitly, and the total widths scale linearly with $\Gamma$10,

$\Gamma$11

At fixed $\Gamma$12, the luminosities become universal functions of $\Gamma$13, and the photon width is roughly two orders of magnitude below the mesonic widths (Maruyama et al., 29 Sep 2025).

Not all proton-controlled synchrotron emission is direct proton radiation. In hadronic $\Gamma$14-decay scenarios, secondary electrons and positrons generated in $\Gamma$15 interactions can radiate synchrotron X-rays. For HESS J1641-463, new NuSTAR data with $\Gamma$16 ks exposure and archival Chandra data yield upper limits of $\Gamma$17 erg cm$\Gamma$18 s$\Gamma$19 in $\Gamma$20–$\Gamma$21 keV and $\Gamma$22 erg cm$\Gamma$23 s$\Gamma$24 in $\Gamma$25–$\Gamma$26 keV, but those limits are not yet deep enough to constrain the primary proton spectrum (Tsuji et al., 2024).

Galactic synchrotron studies make the same distinction explicitly. In FIRE simulations with spectrally resolved CR-MHD, synchrotron emission is produced by electrons and positrons rather than by protons directly, while protons dominate the CR energy budget, set the proton-to-electron ratio, and supply secondary $\Gamma$27. In that framework standard equipartition assumptions underestimate the true emission-weighted magnetic field by factors of $\Gamma$28–$\Gamma$29, and spectral evolution is crucial near galactic centers, where neglecting it can overpredict the emission by factors of $\Gamma$30 or even $\Gamma$31–$\Gamma$32 in the most extreme central regions (Ponnada et al., 2023).

6. Contemporary modeling and astrophysical interpretation

Recent work has increasingly embedded proton synchrotron emission in broader multimessenger and plasma-kinetic frameworks rather than treating it as an isolated spectral component. In reconnection-driven flare models for M87*, a mixed pair–proton plasma in a MAD current sheet yields a three-component radiative partition: pair synchrotron peaks at

$\Gamma$33

proton synchrotron peaks at

$\Gamma$34

and pair inverse Compton produces the TeV component. In that scenario protons are subdominant in number but can dominate the energy budget, and proton synchrotron accounts for approximately $\Gamma$35 to $\Gamma$36 of the total dissipation power (Hakobyan et al., 18 Jul 2025).

The same trend appears in fast surrogate modeling. A CNN-based hadronic blazar framework trained on $\Gamma$37 SOPRANO spectra covers both proton-synchrotron and hybrid lepto-hadronic regimes, evolving the system to equilibrium at $\Gamma$38 and reproducing both electromagnetic and neutrino emission. Applied to TXS 0506+059, a Gaussian neutrino-flux fit favors a hybrid solution with $\Gamma$39 G, $\Gamma$40 cm, $\Gamma$41, and $\Gamma$42, whereas a Poisson event-count fit favors a compact proton-synchrotron solution with $\Gamma$43 G, $\Gamma$44 cm, $\Gamma$45, and $\Gamma$46. For PKS 0735+178, the posterior remains bimodal, with a proton-synchrotron-compatible best fit near $\Gamma$47 G, $\Gamma$48 cm, and $\Gamma$49 (Sahakyan et al., 30 Jun 2025).

The literature represented here repeatedly reaches a restricted rather than universal verdict. Steady proton-synchrotron models for luminous GeV blazars are commonly found to require super-Eddington power or magnetic-field configurations that are difficult to reconcile with other constraints; prompt-GRB implementations often demand $\Gamma$50–$\Gamma$51 G fields and extreme Poynting fluxes; and extended-jet or VHE interpretations can become implausible when the required proton energies approach or exceed $\Gamma$52 eV (Liodakis et al., 2020, Florou et al., 2021, Basumallick et al., 2017). At the same time, specific source classes remain viable: high-peaking blazar humps above $\Gamma$53 GeV, delayed VHE activity in TXS 0506+056, reverse shocks in GRB afterglows, some extended quasar knots, sub-MeV bumps in faint mini lobes, and proton-synchrotron GeV emission from reconnection layers near black holes (Xue et al., 2023, Sunanda et al., 2022, Isravel et al., 2022, Hakobyan et al., 18 Jul 2025).

Within that landscape, proton synchrotron emission is best understood not as a generic hadronic default but as a tightly constrained mechanism whose relevance depends on the detailed balance among acceleration, synchrotron cooling, escape, opacity, and global energetics.