---
title: Proton Synchrotron & Hadronic Models
url: https://www.emergentmind.com/topics/proton-synchrotron-and-hadronic-models
type: topic
---

# Proton Synchrotron & Hadronic Models

A proton synchrotron is a theoretical and computational framework within high-energy astrophysics and astroparticle physics in which the high-energy photons in nonthermal astrophysical sources are produced by ultra-relativistic protons radiating synchrotron emission in strong magnetic fields. Hadronic models constitute a broad class of physical scenarios where, in addition to or instead of electrons, relativistic ions (typically protons) serve as the main energy carriers and radiators in compact objects such as blazar jets and gamma-ray bursts (GRBs). These models contrast with leptonic frameworks, where the emission is primarily due to electrons and positrons upscattering ambient photons. Integrating proton synchrotron processes with the full suite of hadronic interactions—such as photopion production, Bethe–Heitler pair creation, and the ensuing electromagnetic cascades—enables the prediction of both electromagnetic and high-energy neutrino emission, thereby forging a unified multimessenger paradigm for compact jet sources [2411.14218, 1207.5227].

## 1. Physical and Mathematical Foundation of Proton Synchrotron

In a strong magnetic field $B$, a relativistic proton ($\gamma_p \gg1$) emits synchrotron radiation with a characteristic power per particle:
\[
P_{\rm syn}(\gamma_p) = \frac{4}{3}\,\sigma_T\,c\,U_B\,\left( \frac{m_e}{m_p} \right)^2\,\gamma_p^2\,,
\]
where $U_B=B^2/(8\pi)$ is the comoving magnetic energy density, and $\sigma_T$ is the Thomson cross-section. The typical photon energy radiated is
\[
\varepsilon_{\rm syn}(\gamma_p) = \frac{3\,e\,B}{4\pi\,m_p\,c}\,\gamma_p^2\,.
\]
The synchrotron cooling time for these protons is
\[
\tau_{p,\rm syn} = \frac{6\pi\,m_p^3\,c}{\sigma_T\,m_e^2\,B^2\,\gamma_p}\,.
\]
For the observed high-energy $\gamma$-ray emission in jets, this necessitates either very large $B$ or extremely high $\gamma_p$, often approaching or exceeding $\gamma_{p,\max}\sim 10^9$ [1411.5968, 2411.14218].

Modeling these processes, particularly in the context of blazar spectral energy distributions (SEDs), further enforces constraints on the allowed parameter space due to requirements such as confinement (Hillas criterion) and avoidance of excessive internal photon-photon ($\gamma\gamma$) opacity at high $B$ [1207.5227].

## 2. Structure and Regimes of Hadronic Models

Hadronic models for high-energy astrophysical sources generally involve:

- **Pure proton-synchrotron regime:** Protons are the dominant ultrarelativistic species, and their synchrotron emission produces the observed high-energy hump. This typically requires $B$ in the range $1$–$100$ G, $\delta$ (Doppler factor) in the range $15$–$50$, and maximum proton energies $E_{p,\max}\sim 10^{17}$–$10^{19}$ eV, with often super-Eddington total energy requirements [1207.5227, 1411.5968, 2003.10460].
- **Cascade–plus–SSC regime:** At lower $B$ ($0.1$–$1$ G), proton $p\gamma$ and Bethe–Heitler interactions dominate proton energy losses, resulting in secondary $e^\pm$ and photon cascades, accompanied by electron SSC emission, collectively producing the GeV–TeV emission. In mixed (hybrid) lepto-hadronic models, both electron and proton populations contribute significantly to the high-energy component [1411.5968, 1912.07448].

Transitions between particle-dominated and magnetic-dominated regimes are set by the ratio $U_p/U_B$, which, along the minimum Doppler-factor curve $\delta_{\min}(B)$, can exceed $10^3$ (particle-dominated) or be $\ll1$ (magnetically dominated) [1207.5227].

## 3. Computational Formulation and Code Systematics

Contemporary implementations solve coupled, stationary or time-dependent kinetic equations for the distributions of protons, electrons, photons, and all secondary particles, incorporating injection, radiative losses, escape, and all relevant hadronic interactions. For instance, 
\[
\frac{\partial N_p(\gamma_p)}{\partial t} = -\frac{\partial}{\partial \gamma_p}\left[ \dot\gamma_{\rm tot} N_p(\gamma_p) \right] - \frac{N_p(\gamma_p)}{t_{\rm esc}} + Q_p(\gamma_p)
\]
with 
\[
\dot\gamma_{\rm tot}= \dot\gamma_{\rm syn} + \dot\gamma_{p\gamma} + \dot\gamma_{\rm BH}\,,
\]
and analogous terms for electrons and photons [1411.5968, 1111.0557, 2411.14218].

Comprehensive code comparisons [2411.14218] between numerical schemes (AM³, ATHENA, B13, LeHa-Paris, LeHaMoC) reveal excellent agreement in spectral shapes, with a normalization uncertainty envelope of $\pm40\%$, which should be systematically included in statistical-hadronic SED modeling.

Surrogate modeling via convolutional neural networks (CNNs), trained on extensive parameter grids from kinetic codes such as SOPRANO, is now used for rapid Bayesian sampling and observational fitting, enabling multidimensional parameter exploration on timescales orders of magnitude faster than traditional solvers [2506.23885].

## 4. Energetics, Constraints, and Observational Implications

The necessary conditions for matching observed high-energy emissions with proton-synchrotron models are highly restrictive. For blazars with observed $\gtrsim 100$ MeV–TeV $\gamma$-ray luminosities, minimum jet powers $P_{j, \min}$ calculated for $\sim 150$ sources are typically $\gtrsim 10^2$ times larger than the Eddington luminosity $L_{\rm Edd}$, accretion-disk luminosity $L_{\rm disk}$, or spin-extraction (Blandford–Znajek) power $P_{\rm BZ}$ [2003.10460]. Achieving this requires either extreme amplification ($f_{\rm amp}\sim 30$) of the jet magnetic field or a production site inside the broad-line region, often incompatible with VLBI or variability constraints.

A plausible implication is that, where relativistic hadrons are present, they can only contribute a radiatively subdominant high-energy (steady-state) component rather than accounting for the dominant $\gamma$-ray luminosity of powerful blazars [2003.10460].

The situation is somewhat more favorable for:
- Ultra-high-frequency-peaked BL Lacs (UHBLs), where sub-Eddington jet powers are achievable for $B\sim 1$–$100$ G, $\delta\simeq 30$, and $E_{p,\max} \lesssim 10^{19}$ eV [1411.5968].
- Gamma-ray burst afterglows, where proton-synchrotron can explain GeV–TeV afterglow emission if the allowed parameter space (large initial Lorentz factor, $B\sim 10^3$–$10^6$ G, high $\epsilon_B$) is realized [1004.3330].

## 5. Polarization as a Discriminant and Multimessenger Diagnostics

Polarization predictions provide an incisive test of proton-synchrotron models. Proton-synchrotron emission in a perfectly ordered field produces a frequency-independent maximum linear polarization,
\[
\Pi_{\rm pl} = \frac{p+1}{p+7/3}\,,
\]
reaching $70\%$–$75\%$ for typical proton indices [1307.4187]. Hadronic models predict substantially higher maximal X-ray and $\gamma$-ray polarization than leptonic SSC or external Compton scenarios, especially in low- and intermediate-synchrotron-peaked blazars. High values of observed high-energy polarization ($\gtrsim f\cdot 50\%$ with $f$ the field-order correction) strongly favor hadronic emission.

Multimessenger diagnostics further leverage the correlated prediction of high-energy neutrinos. Pure proton-synchrotron blazar fits predict neutrino SED peaks at $0.1$–$10$ EeV, with peak all-flavor fluxes $\sim 10^{-4} \,$ of the $\gamma$-ray flux. In the case of TXS 0506+056, only hybrid (SSC + hadronic cascade) scenarios reproduce the observed $E_\nu \sim 0.3$–$1$ PeV neutrino coincident with the $\gamma$-ray flare [1912.07448, 1807.04335, 2506.23885].

## 6. Generalizations and Applications across Source Classes

The proton-synchrotron and broader hadronic framework is applied across multiple source classes:
- **Blazars:** SED modeling for FSRQs, BL Lacs, and UHBLs, with code implementations (SOPRANO, AM³, etc.) allowing parameter studies over $B$, $\delta$, $R_{\rm blob}$, and particle injection rates [2411.14218, 2506.23885].
- **Gamma-Ray Bursts:** Hadronic models for the prompt and afterglow emission of GRBs invoke proton-synchrotron and photohadronic cascades to account for high-energy ($>$100 MeV) LAT data, with strong constraints from energy budget and neutrino upper limits [1210.7802, 2102.02501, 1004.3330].
- **Jet Simulation and Microphysics:** Integration of hadronic energy losses and proton-synchrotron cooling into particle-in-cell and general-relativistic MHD schemes is achieved using relativistic Boris-pusher and guiding-center algorithms, incorporating continuous and stochastic hadronic collision terms, validated against analytic benchmarks [2410.22781].

These methodologies enable not only electromagnetic SED fitting, but increasingly detailed predictions for neutrino event rates and polarimetric signatures, with systematic uncertainties in SED normalization ($\pm40\%$) now quantified for robust inference [2411.14218, 2506.23885].

## 7. Limitations, Systematics, and Outlook

The proton-synchrotron model, while conceptually attractive for linking $\gamma$-ray and neutrino astrophysics, faces severe constraints:
- Energetic requirements for steady blazar emission are typically orders of magnitude above the Eddington or disk luminosity [2003.10460].
- Viable parameter regimes require high $B$, large $\delta$, and ultrarelativistic proton populations that challenge standard acceleration scenarios and jet energetics.
- Secondary $e^\pm$ production via $p\gamma$ and $\gamma\gamma$ interactions is generally subdominant in the TeV band for high $B$, but can reshape the spectrum below the peak unless $\gamma$-ray emission zones are highly compact and magnetized [2102.02501].

These constraints suggest a radiatively subdominant role for protons in most steady blazar jets, though episodic or highly-magnetized environments may admit plausible hadronic dominance, particularly in orphan flares or narrow SED intervals [1207.5227, 1411.5968, 2411.14218]. The integration of deep learning surrogates and code cross-comparison establishes a reproducible, systematics-aware foundation for future multimessenger studies [2506.23885].

Source: https://www.emergentmind.com/topics/proton-synchrotron-and-hadronic-models