Pair Cascades and Acceleration Gaps

A pulsar needs plasma to carry the current demanded by its global magnetosphere, yet the same plasma screens the electric field that accelerates it. An acceleration gap is the self-consistent failure of complete screening: an electric field parallel to the magnetic field accelerates a primary particle, its photons convert into electron–positron pairs, and the new charges modify the field. Where that cycle closes—polar cap, slot gap, outer gap, separatrix, or current sheet—depends on charge supply, current demand, spacetime, magnetic geometry, and radiation physics. No single “death line” is independent of those choices.

Required background. Frames, Coframes, and Anholonomy supplies local electric and magnetic components, particle energies, and curvature scales.

Helpful background. Flat-Background Hilbert Stress-Energy from Metric Variation supplies electromagnetic energy flow; The Einstein Equation supplies redshift and frame dragging; Relativistic Matter and Compact Objects supplies the conducting surface and magnetic geometry; Computational Relativity supplies kinetic resolution and convergence checks.

From charge deficit to an accelerator

In Gaussian electromagnetic units with $c=\hbar=1$, the total electric field obeys $\boldsymbol\nabla\mathbin{\cdot}\boldsymbol E=4\pi\rho$. A reduced flux-tube gap model subtracts the corotation or global background field. If $\boldsymbol\nabla\mathbin{\cdot}\boldsymbol E_{\rm cor}=4\pi\rho_{\rm req}$, the residual accelerator field $\boldsymbol E_{\rm gap}=\boldsymbol E-\boldsymbol E_{\rm cor}$ obeys

$$\boldsymbol\nabla\mathbin{\cdot}\boldsymbol E_{\rm gap} =4\pi(\rho-\rho_{\rm req}).$$

$\rho_{\rm req}$ is the charge density needed by the corotating or global current solution; in the simplest flat aligned limit it becomes the Goldreich–Julian scale. If the surface and incoming plasma cannot supply the required charge and current, $E_\parallel$ develops. A primary electron or positron reaches Lorentz factor $\gamma$ until radiation reaction, pair screening, or the finite voltage ends acceleration.

For a field-line curvature radius $\rho_c$, curvature photons have a characteristic energy

$$\varepsilon_c\simeq\frac{3\gamma^3}{2\rho_c},$$

with $\hbar=c=1$. They may convert through one-photon magnetic pair creation, whose rate depends steeply on energy and the field perpendicular to the photon, or through photon–photon collisions. For two photons of energies $\varepsilon_1$ and $\varepsilon_2$ meeting at angle $\vartheta$, the vacuum threshold is

$$\varepsilon_1\varepsilon_2(1-\cos\vartheta)\geq2m_e^2.$$

The daughters radiate and may make further pairs. Pair multiplicity, spectrum, and screening length are outputs of the cascade, not free synonyms for plasma density.

Current demand changes the polar-cap result

A gap cannot be solved from the local magnetic field alone because the global magnetosphere sets the current. In one-dimensional space-charge-limited polar flow, Timokhin and Arons found mild acceleration without pairs for some sub-Goldreich–Julian currents, but intermittent discharges for super- Goldreich–Julian and return currents. Their current-controlled formulation, Monte Carlo photons, and boundary assumptions are documented in Timokhin and Arons 2013, §§2–6, Figs. 2–23, and Appendix A. The result demonstrates time dependence and current sensitivity; its one-dimensional flux tube cannot determine global current closure.

Global particle-in-cell calculations connect local supply to open flux and a current sheet, but their physical scale separations and pair prescriptions are necessarily compressed. Three-dimensional oblique simulations with threshold pair injection found force-free-like structure and inclination-dependent pair regions in Philippov, Spitkovsky, and Cerutti 2015, pp. L19-1–L19-5 and Figs. 1–4. Their photons and secondary injection implement a specified effective cascade; they do not resolve every surface, QED, and radio-emission scale of a real pulsar.

Competing gap locations and general relativity

Evidence cutoff: 2026-08-18. Polar-cap models emphasize strong near-surface fields and magnetic conversion. Slot-gap models extend an unscreened region along the last open lines. Outer-gap models use photon–photon creation near the null-charge surface, while global kinetic calculations often place strong acceleration and radiation in the return-current layer and equatorial current sheet. These are competing formulations of where charge starvation persists; more than one region may operate in one magnetosphere.

The comparison below fixes the evidence date while keeping unlike models from being mistaken for measurements of one universal gap. Each row asks a different closure question and compresses different parts of the cascade.

Model	Field geometry and assumptions	Seed photons and radiative channels	Pair multiplicity and global-current boundary	Method and evidence type	Result	Limitation	Exact source locator	Evidence date
Polar-cap kinetic discharge	Local/rescaled one-dimensional tube with constant $\rho_{\rm GJ}$ in the main pair runs, prescribed curvature radius $\rho_c=10^6\,\mathrm{cm}$, and $B=B_0$ below $x_B$ but $B=0$ above it	Curvature photons undergo one-photon magnetic conversion; Compton seeds are excluded from the cited runs	Multiplicity is time dependent, while the current ratio is imposed by an external global solution	Kinetic discharge calculation with Monte Carlo photon conversion	Sub-Goldreich–Julian flow can remain quiet; super-Goldreich–Julian and return currents are generally intermittent, with a narrow fine-tuned quasi-stationary exception	The one-dimensional geometry over-coheres the discharge; transverse structure and continued photon transport are essential for a physical pair yield	Timokhin and Arons 2013, § 4, § 6 including Eq. (15) and Figs. 11–23, and § 7	2026-08-18
Slot gap	General-relativistic space-charge-limited flow in a prescribed static dipole, including the lapse and frame-dragging contribution, with narrow equipotential slot boundaries	Curvature photons set the screening pair-formation front; the cascade code also tracks inverse-Compton scattering of thermal X-rays, but those pairs alone do not close the slot gap	The pair-formation front and gap width are imported from earlier cascade calculations; this paper does not evolve multiplicity or the global return current	Analytic GR electrodynamics plus a particle-and-photon cascade calculation	A small residual parallel field maintains high-altitude acceleration and produces caustic curvature emission	The imported gap width, static dipole, and local boundary conditions are not a global kinetic closure	Muslimov and Harding 2004, § 1, § 2 Eqs. (1)–(4), § 3 Eqs. (40)–(59), and § 3.2.2	2026-08-18
Outer gap	Charge starvation near a null-charge surface in a prescribed outer-magnetosphere geometry	Curvature or inverse-Compton gamma rays interact with an assumed soft-photon field through photon–photon conversion	Pair closure bounds the gap, but multiplicity and boundary current are closure assumptions rather than a global output	Semianalytic gap electrodynamics and radiative closure	An accelerator can persist far from the polar cap	Target photons, geometry, and current boundaries are imposed; the model predates current-consistent global kinetics	Cheng, Ho, and Ruderman 1986, part I, pp. 500–521	2026-08-18
Return-current layer and current sheet	Three-dimensional oblique magnetosphere with a compressed kinetic scale hierarchy	Effective threshold injection stands in for the photon cascade; separate curvature and Compton photon transport is not resolved	Pair supply follows the injection prescription, while the global current distribution emerges from the simulation	Global particle-in-cell calculation	Pair formation and acceleration occur in the return-current layer and current sheet; quantified current-sheet dissipation within $2R_{\rm LC}$ decreases from about 20% for the aligned rotator to 3% for the orthogonal rotator	The imposed injection threshold does not predict a physical multiplicity, surface cascade, or coherent radio emission	Philippov, Spitkovsky, and Cerutti 2015, §§ 4.1–4.3, Figs. 2–5, and § 5	2026-08-18

General relativity can change the competition. Frame dragging alters the charge density required for corotation and can turn a flat-spacetime weak polar cap into a spacelike-current accelerator. Axisymmetric aligned GR-PIC simulations of compact millisecond-star models find this activation in Torres et al. 2024, §§3–5, Figs. 3–6, and Appendix A. The credible contrary evidence is not that flat calculations were numerically wrong: they address a different spacetime and current ratio. The GR result is itself limited to aligned axisymmetry, chosen compactness and dipole geometry, finite particle number, and compressed kinetic scales.

An observed radio cutoff or gamma-ray light curve does not uniquely choose a gap. Beaming, viewing angle, magnetic multipoles, surface temperature, and selection against faint or intermittent pulsars can move an empirical boundary. Coherent radio emission remains an unresolved step beyond producing pairs; pair supply is necessary in many models but is not a derivation of a radio spectrum.

Comparing results and possible new tests

A comparison must record spacetime and compactness; dipole or multipolar field; obliquity; surface work function and charge extraction; imposed global current; curvature, synchrotron, inverse-Compton, magnetic-conversion, and two-photon processes; dimensionality; boundary conditions; grid and particles per cell; stellar-to-light-cylinder ratio; skin-depth and gyroradius resolution; pair injection rule; duration; and diagnostics for voltage, multiplicity, current, and energy conservation.

Evidence for a self-sustained cascade requires that pairs generated from the modeled photons maintain the demanded current after arbitrary seed plasma is removed. Evidence for screening requires convergence of the residual $E_\parallel$, not only a force-free-looking field image. A death boundary requires a selection-aware population comparison rather than the absence of one simulated discharge.

An independent budget check begins with the electromagnetic work $\dot E_{\rm acc}=\int_{\rm gap}\boldsymbol J\mathbin{\cdot}\boldsymbol E\,dV$. For a quasi-steady one-dimensional channel in which a primary current $I_{\rm prim}=e\dot N_{\rm prim}$ crosses one potential drop $\Delta\Phi$, this reduces to $\dot E_{\rm acc}=\lvert I_{\rm prim}\Delta\Phi\rvert$. If $\kappa_{\rm pair}$ denotes created pairs per injected primary and a fraction $\eta_\gamma$ of the work reaches pair-producing photons, then

$$\kappa_{\rm pair} \leq\frac{\eta_\gamma\dot E_{\rm acc}} {2m_e\dot N_{\rm prim}} =\frac{\eta_\gamma e\lvert\Delta\Phi\rvert}{2m_e}.$$

This bound applies only to that restricted channel. Photon escape, daughter kinetic energy, and radiation make it stricter. Charge continuity supplies a separate test: with $I_{\rm in}$ and $I_{\rm out}$ defined as positive entering and leaving current magnitudes, $\dot Q_{\rm gap}=I_{\rm in}-I_{\rm out}$. A calculation that reports a large multiplicity while violating either budget has not demonstrated pair closure.

The main unresolved issues are coherent radio emission, realistic multipolar surface fields, ion extraction, QED rates in compressed simulations, coupling among polar and outer gaps, and how kinetic dissipation converges in the current sheet. The interpretation could change with a global three-dimensional GR-PIC suite with physical radiation and converged scale ratios, an independently reproduced pair-closure boundary, or an official population data release with a forward selection model that changes the empirical death region.

Check your understanding

A $1\,\mathrm{GeV}$ photon collides head-on with a soft photon. What minimum soft-photon energy reaches the two-photon pair threshold?

Solution

For a head-on collision, $1-\cos\vartheta=2$, so the threshold reduces to $\varepsilon_1\varepsilon_2=m_e^2$. Hence

$$\varepsilon_2 =\frac{(0.511\,\mathrm{MeV})^2}{1000\,\mathrm{MeV}} =2.61\times10^{-4}\,\mathrm{MeV} =0.261\,\mathrm{keV}.$$

At a smaller collision angle the required energy is larger. The calculation is purely kinematic; a substantial optical depth also needs sufficient target photon density and the angle-dependent cross section.

References

• Cheng, K. S., C. Ho, and M. Ruderman. “Energetic Radiation from Rapidly Spinning Pulsars. I. Outer Magnetosphere Gaps.” The Astrophysical Journal 300 (1986): 500–521. DOI.
• Muslimov, Alex G., and Alice K. Harding. “High-Altitude Particle Acceleration and Radiation in Pulsar Slot Gaps.” The Astrophysical Journal 606 (2004): 1143–1153. DOI.
• Philippov, Alexander A., Anatoly Spitkovsky, and Benoît Cerutti. “Ab Initio Pulsar Magnetosphere: Three-Dimensional Particle-in-Cell Simulations of Oblique Pulsars.” The Astrophysical Journal Letters 801 (2015): L19. DOI.
• Timokhin, Andrey N., and Jonathan Arons. “Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars’ Force-Free Magnetospheres: Space Charge Limited Flow.” Monthly Notices of the Royal Astronomical Society 429 (2013): 20–54. DOI.
• Torres, Rui, Thomas Grismayer, Fábio Cruz, Ricardo A. Fonseca, and Luís O. Silva. “OSIRIS-GR: General Relativistic Activation of the Polar Cap of a Compact Neutron Star.” New Astronomy 112 (2024): 102261. DOI.