---
title: Proton Synchrotron Emission Physics
url: https://www.emergentmind.com/topics/proton-synchrotron-emission
type: topic
---

# Proton Synchrotron Emission Physics

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 [2304.13893][1912.02185].

## 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}\) [2304.13893].

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 \((m_p/m_e)^3\) suppression in the cooling efficiency [2304.13893].

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
\[
t^{\rm obs}_{\rm cool,S,p} = t^{\rm obs}_{\rm cool,S,e}\left(\frac{m_p}{m_e}\right)^{5/2},
\]
which is the central reason proton synchrotron has been proposed for incomplete-cooling prompt spectra [1912.02185]. In TXS 0506+056, the jet-frame proton energy required to produce synchrotron photons at observed critical frequency \(\nu_c\) is written as
\[
E'_p \simeq 4.38\times 10^{19}\ \text{eV}\, \left(\nu_{c,25} B_3\right)^{1/2}\frac{(1+z)}{\delta},
\]
and the observer-frame cooling time as
\[
t^{\rm obs}_{\rm syn,p}\simeq 7.52\ B_{2.4}^{-2}\ \delta_{16}^{-1}\ E_{p,19}^{\prime-1}\ \text{days},
\]
showing directly how long-lived VHE proton synchrotron activity can emerge when the magnetic field is of order a few gauss [2211.02493].

Analytical treatments usually supplement these relations with a Hillas-type confinement condition and a transparency condition against internal \(\gamma\gamma\) absorption. In one-zone blazar applications, the allowed parameter space in the \((B,\delta_{\rm D})\) plane is therefore carved out simultaneously by the Eddington constraint, the acceleration requirement, and the condition \(\tau_{\gamma\gamma}<1\) [2304.13893].

## 2. Energetics in blazars

The most systematic energetic critique of proton synchrotron models in blazars is the study of steady \(\gamma\)-ray emission in a sample of 145 sources. In that framework the absolute power of a two-sided jet is written as
\[
P_j = 2\pi r^{\prime 2} \beta \Gamma^2 c \sum_{i=B,e,p} \left(u_i^\prime + p_i^\prime\right) + P_j^{r} + P_j^{c},
\]
and the quantity \(P_{j,\min}\) 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 \(P_{j,\min}\) is already a conservative lower limit [2003.10460].

For the 145-source sample, the observables were Monte Carlo sampled with \(10^4\) realizations per source, using uncertainties of \(0.5\) dex in luminosities, \(0.3\) dex in peak frequencies, and a variability-timescale distribution with
\[
\mu = 10^5\,{\rm s}, \qquad \sigma = 3\times10^4\,{\rm s}.
\]
The central result is the so-called energy crisis: for most sources,
\[
P_{j,\min}\sim 10^2 L_d,
\]
and the same lower bound is typically about two orders of magnitude above \(L_{\rm Edd}\) and the optimistic Blandford–Znajek estimate \(P_{\rm BZ}\). The derived magnetic fields imply either local amplification by a factor of about \(30\) or a \(\gamma\)-ray production site at sub-parsec scales; the former corresponds to a median amplification factor \(f_{\rm amp}\approx 27\), while the latter yields a median \(z_{\rm em}\sim 0.03\) pc and places the emission well inside the BLR. The predicted neutrino emission peaks at \(\sim 0.1-10\) EeV, with typical muon-neutrino peak fluxes \(\sim 10^{-4}\) of the peak \(\gamma\)-ray flux [2003.10460].

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 \(R\sim10^{16}\)–\(10^{17}\,\mathrm{cm}\), \(B\lesssim 10\) G, and \(\delta_{\rm D}\lesssim 30\). For humps peaking in the \(0.1\)–\(10\) GeV range, however, the model typically requires super-Eddington jet power and, if \(\delta_{\rm D}\lesssim 30\), magnetic fields that can exceed \(10^3\) G. For humps peaking in the \(10\)–\(100\) keV band, the outcome depends strongly on luminosity: Mrk 421 admits an allowed region, whereas 3C 279 does not [2304.13893].

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 [2003.10460][2304.13893].

## 3. Gamma-ray bursts and other transients

Prompt-GRB applications were motivated by spectra that exhibit
\[
F_\nu \propto \nu^{1/3}
\]
below a break and
\[
F_\nu \propto \nu^{-1/2}
\]
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 \(\sim 100\) keV with a cooling time of order \(1\) s, the proton-synchrotron interpretation was argued to work with \(B'\sim 10^6\) G and \(R\sim 10^{13}\) cm, whereas the electron-synchrotron alternative would require \(B'\sim 1\) G and \(R\sim 10^{16}\) cm, in tension with prompt variability [1912.02185].

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
\[
B \sim (1-20)\times10^6~{\rm G},
\]
and for \(\Gamma=300\) the median jet Poynting luminosity is
\[
L_{B,j}\approx 5\times10^{54}\ {\rm erg\,s^{-1}},
\]
with strong low-energy spectral contamination from secondary pairs unless \(\Gamma \gtrsim 10^3\). On that basis the model is described as strongly disfavoured for the low-energy spectral breaks of prompt GRBs [2102.02501]. A related study of Bethe–Heitler pair production finds two regimes: at high \(\Gamma\), 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 \(\Gamma\), small radius, and high luminosity, Bethe–Heitler cooling drives the spectrum toward a single fast-cooling power law \(F_\nu\propto \nu^{-1/2}\) across the entire GBM/Swift band, which is incompatible with observations [2112.07231].

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 \(0.2\)–\(1\) TeV MAGIC emission reproduces the data with isotropic explosion energy \(\sim 10^{54.5}\) erg, ambient density \(\sim 10\)–\(100~{\rm cm^{-3}}\), a few-percent accelerated fractions, and protons reaching a few \(10^{20}\) eV [2210.02363]. For GRB 221009A, reverse-shock proton synchrotron has been used to explain \(\gtrsim\)TeV photons, including the possibility of \(\sim 18\) TeV photons with reasonable EBL models, and a structured-jet analysis associates the rise of the reverse-shock proton-synchrotron component with the \(\sim 500\)–\(800\) s hardening interval and the arrival of the \(\sim 13\) TeV photon [2211.05754][2311.13671].

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 \(B\simeq 2.4\) G, \(\delta=16\), \(R'=1.23\times10^{16}\) cm, and a proton luminosity \(L'_p \simeq 10^{47}\) erg s\(^{-1}\), with proton synchrotron cooling times in the \(\sim 11\) to \(44\) day interval matching the observed \(\sim 45\)-day VHE episode [2211.02493].

## 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 \(R=1\) kpc, \(B=30\) mG, \(\Gamma_j=3\), \(\delta_D=1\), and requires a total luminosity of about \(10^{46}\,\mathrm{erg\,s^{-1}}\). In the escape-dominated case, the parameters become \(R=0.8\) kpc and \(B=10\) mG, with powers \(L_{j,B}=5.9\times10^{44}\,\mathrm{erg\,s^{-1}}\) and \(L_{j,p}=1.0\times10^{44}\,\mathrm{erg\,s^{-1}}\), making the escape-dominated solution the more favorable one [1406.5978].

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 \(10^{20}\)–\(10^{21}\,\mathrm{eV}\), but the luminosity budgets are comparatively modest: for the combined knots of PKS 0637-752 the total power is about \(6\times 10^{43}\,\mathrm{erg\,s^{-1}}\), or about \(0.6\%\) of Eddington, and for 3C 273 the required power is around \(5\times 10^{43}\,\mathrm{erg\,s^{-1}}\), or about \(0.5\%\) of Eddington [1511.00258]. AP Librae represents the opposite extreme. There the extended-jet proton-synchrotron interpretation of the VHE emission requires \(B\sim 1\) mG, \(E_{p,\max}\approx 3.98\times 10^{21}\) eV, and
\[
P_{\rm jet}\sim 4.86\times 10^{48}\ \text{erg s}^{-1},
\]
which is more than \(100\) times the Eddington luminosity of AP Librae, leading that scenario to be described as unlikely [1706.04895].

Compact radio lobes introduce another characteristic proton-synchrotron signature. In parsec-scale mini lobes with \(R=2\) pc and \(B=0.1\) G, the proton synchrotron characteristic photon energy is
\[
\nu_{p,\rm syn}\approx 0.7 \left(\frac{E_p}{10^{18}\,\mathrm{eV}}\right)^2 \left(\frac{B}{0.1\,\mathrm{G}}\right) \mathrm{MeV},
\]
placing the direct proton-synchrotron bump in the sub-MeV band. Its visibility depends strongly on the primary electron synchrotron luminosity: for \(L_e<10^{43}\,\mathrm{erg\,s^{-1}}\) the sub-MeV bump appears clearly, whereas for \(L_e\sim10^{45}\,\mathrm{erg\,s^{-1}}\) leptonic emission hides the hadronic signature. In intermediate cases, synchrotron from secondary \(e^\pm\) produced in the photopion cascade can emerge in the GeV–TeV range, yielding a double-bump hadronic spectrum [1006.3384].

## 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,
\[
\left[ \gamma_\mu (i\partial^\mu - eA^\mu) - m_N - \frac{e\kappa_p}{2m_N}\sigma_{\mu\nu} (\partial^\mu A^\nu - \partial^\nu A^\mu) \right]\psi(x)=0,
\]
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 \(10^3\) [1503.05635].

A later relativistic quantum treatment extends this framework to photon, pion, and \(\rho\)-meson production and derives a scaling rule for transitions between very large Landau levels. There the proton curvature parameter is
\[
\chi_p=\frac{eBE_i}{M_p^3},
\]
the transition probabilities can be extrapolated to \(n_i\gtrsim 10^{15}\), recoil is included explicitly, and the total widths scale linearly with \(B\),
\[
\Gamma_A \propto B.
\]
At fixed \(\chi_p\), the luminosities become universal functions of \(\chi_p\), and the photon width is roughly two orders of magnitude below the mesonic widths [2509.24366].

Not all proton-controlled synchrotron emission is direct proton radiation. In hadronic \(\pi^0\)-decay scenarios, secondary electrons and positrons generated in \(pp\) interactions can radiate synchrotron X-rays. For HESS J1641-463, new NuSTAR data with \(82\) ks exposure and archival Chandra data yield upper limits of \(\sim 6\times 10^{-13}\) erg cm\(^{-2}\) s\(^{-1}\) in \(2\)–\(10\) keV and \(\sim 3\times 10^{-13}\) erg cm\(^{-2}\) s\(^{-1}\) in \(10\)–\(20\) keV, but those limits are not yet deep enough to constrain the primary proton spectrum [2404.11012].

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 \(e^\pm\). In that framework standard equipartition assumptions underestimate the true emission-weighted magnetic field by factors of \(2\)–\(3\), and spectral evolution is crucial near galactic centers, where neglecting it can overpredict the emission by factors of \(\sim 6\) or even \(\sim 10\)–\(50\) in the most extreme central regions [2309.04526].

## 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
\[
\varepsilon_{\rm syn}^{\pm}\approx 24~{\rm MeV},
\]
proton synchrotron peaks at
\[
\varepsilon_{\rm syn}^{i}\approx 43~{\rm GeV},
\]
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 \(5\%\) to \(20\%\) of the total dissipation power [2507.14002].

The same trend appears in fast surrogate modeling. A CNN-based hadronic blazar framework trained on \(7\times 10^6\) SOPRANO spectra covers both proton-synchrotron and hybrid lepto-hadronic regimes, evolving the system to equilibrium at \(t=4t_{\rm dyn}\) and reproducing both electromagnetic and neutrino emission. Applied to TXS 0506+059, a Gaussian neutrino-flux fit favors a hybrid solution with \(B\simeq 0.03\) G, \(R\sim 5.8\times 10^{17}\) cm, \(\delta\sim 15.5\), and \(L_p\sim 1.5\times 10^{51}\,\mathrm{erg\,s^{-1}}\), whereas a Poisson event-count fit favors a compact proton-synchrotron solution with \(B\sim 5.3\times10^2\) G, \(R\sim 3.4\times10^{14}\) cm, \(\delta\sim 23.2\), and \(\gamma_{p,\max}\sim 10^8\). For PKS 0735+178, the posterior remains bimodal, with a proton-synchrotron-compatible best fit near \(B\sim 7.5\times10^2\) G, \(R\sim 3.6\times10^{14}\) cm, and \(\delta\sim 26\) [2506.23885].

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 \(10^6\)–\(10^7\) G fields and extreme Poynting fluxes; and extended-jet or VHE interpretations can become implausible when the required proton energies approach or exceed \(10^{21}\) eV [2003.10460][2102.02501][1706.04895]. At the same time, specific source classes remain viable: high-peaking blazar humps above \(\sim 10\) 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 [2304.13893][2211.02493][2210.02363][2507.14002].

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.

Source: https://www.emergentmind.com/topics/proton-synchrotron-emission