Force-Free Pulsar Magnetospheres

A force-free pulsar magnetosphere is a plasma-filled electromagnetic solution in which charges and currents reshape the rotating dipole while their inertia is neglected. Closed field lines corotate inside a separatrix, open flux carries energy and angular momentum away, and an equatorial current sheet completes the circuit. The light cylinder is a regularity surface, not an outer boundary. Computing this structure requires explicit prescriptions for the stellar conductor, current sheet, outer boundary, and constraint control; this page gives that rerun contract without claiming an execution.

Required background. Frames, Coframes, and Anholonomy supplies locally measured electric and magnetic components and drift frames.

Helpful background. Flat-Background Hilbert Stress-Energy from Metric Variation supplies electromagnetic energy and angular-momentum flux; The Einstein Equation supplies relativistic corrections near the star; Relativistic Matter and Compact Objects supplies the conducting stellar boundary; Computational Relativity supplies hyperbolic evolution and constraint checks.

Force-free closure and its domain

In Gaussian electromagnetic units with $c=1$, Maxwell evolution obeys

$$\partial_t\boldsymbol B=-\boldsymbol\nabla\times\boldsymbol E, \qquad \partial_t\boldsymbol E=\boldsymbol\nabla\times\boldsymbol B-4\pi\boldsymbol J,$$

with $\boldsymbol\nabla\mathbin{\cdot}\boldsymbol B=0$ and $\boldsymbol\nabla\mathbin{\cdot}\boldsymbol E=4\pi\rho_e$. Force-free closure sets the Lorentz-force density to zero,

$$\rho_e\boldsymbol E+\boldsymbol J\times\boldsymbol B=0, \qquad \boldsymbol E\mathbin{\cdot}\boldsymbol B=0, \qquad B^2-E^2\gt0.$$

The last inequality guarantees a local frame in which the electric field vanishes. A formal solution that develops $E^2\gt B^2$ is outside ideal force-free validity; clipping $E$ can stabilize a code but acts as an unresolved dissipation model.

The corotation field is $\boldsymbol E=-(\boldsymbol\Omega\times \boldsymbol r)\times\boldsymbol B$. Applying Gauss’s law and $\boldsymbol\nabla\times(\boldsymbol\Omega\times\boldsymbol r)=2\boldsymbol\Omega$ gives

$$4\pi\rho_e=-2\boldsymbol\Omega\mathbin{\cdot}\boldsymbol B +(\boldsymbol\Omega\times\boldsymbol r)\mathbin{\cdot} (\boldsymbol\nabla\times\boldsymbol B).$$

This derivation shows the approximation in the familiar charge density: neglect the current-dependent second term well inside the light cylinder. Corotation fixes a charge density, while force-free closure also fixes a current; meeting one requirement does not automatically meet the other.

Near a corotating aligned region, the charge-density scale reduces to the Goldreich–Julian result $\rho_{\rm GJ}\simeq-\boldsymbol\Omega\mathbin{\cdot}\boldsymbol B/(2\pi)$ in Gaussian units. Its derivation and charge-separated interpretation appear in Goldreich and Julian 1969, §III, Eqs. (5)–(9). Force-free theory assumes a carrier supply dense enough to provide the required charge and current, including both signs wherever a one-sign flow cannot meet that demand. Charge-separated force-free configurations are possible in restricted current regimes. The closure cannot decide where pairs are created or what happens in a charge-starved acceleration gap.

Open flux and light-cylinder regularity

In cylindrical coordinates $(\varpi,\phi,z)$, define the flux function and current function by

$$\boldsymbol B_p=\frac{\boldsymbol\nabla\Psi\times\hat{\boldsymbol\phi}}{\varpi}, \qquad B_\phi=\frac{A(\Psi)}{\varpi}, \qquad\boldsymbol E=-\Omega\boldsymbol\nabla\Psi.$$

Thus $2\pi\Psi$ is poloidal magnetic flux relative to the axis and the enclosed current is $I=A/2$ in the stated Gaussian $c=1$ convention. For constant field-line angular velocity $\Omega$, the stationary aligned force balance reduces to

$$\bigl(1-\Omega^2\varpi^2\bigr) \left(\Psi_{,\varpi\varpi}-\frac{\Psi_{,\varpi}}{\varpi} +\Psi_{,zz}\right) -2\Omega^2\varpi\Psi_{,\varpi}=-A\frac{dA}{d\Psi}.$$

It is elliptic on either side of the light cylinder and degenerate at $\varpi=R_{\rm LC}=\Omega^{-1}$. A differentiable crossing therefore requires

$$A\frac{dA}{d\Psi}=\frac{2}{R_{\rm LC}}\Psi_{,\varpi} \qquad(\varpi=R_{\rm LC}).$$

The current function and flux geometry must be solved together: an arbitrary prescribed current generally fails this condition. These definitions and the pulsar equation are those of Contopoulos, Kazanas, and Fendt 1999, §2, Eqs. (2)–(7), written here with dimensional cylindrical coordinates rather than coordinates scaled by $R_{\rm LC}$. The solution contains a closed zone, a separatrix terminating near a Y point, open polar flux, and an equatorial current sheet. Contopoulos, Kazanas, and Fendt obtained the canonical self-consistent axisymmetric solution by iterating the light-cylinder regularity condition in Contopoulos, Kazanas, and Fendt 1999, §§2–4, Eqs. (1)–(11).

Time-dependent three-dimensional force-free evolution gives, for a dipole moment convention matched to the vacuum page,

$$L_{\rm FFE}\simeq\mu^2\Omega^4(1+\sin^2\alpha).$$

Spitkovsky reports the oblique solutions, current sheet, and luminosity fit in Spitkovsky 2006, pp. L51–L54, Eq. (3), and Fig. 3. The nonzero aligned luminosity is the essential contrast with a point vacuum dipole. The fitted coefficient is a benchmark for the stated force-free setup, not a kinetic pair-plasma theorem.

Inputs and environment boundary

Continuum inputs are the stellar radius $R_*/R_{\rm LC}$, rotation rate, dipole moment and inclination, initial magnetic field, and flat or specified relativistic background. At the conducting stellar surface, impose the normal magnetic field and tangential corotation electric field consistently. Record how surface charge and outgoing current are supplied. The outer domain must extend several light-cylinder radii and use an outgoing boundary whose reflection level is measured.

The numerical record includes coordinates, dimensionality, grid or mesh refinement, polar treatment, reconstruction and Riemann solver, time integrator, divergence control, force-free current implementation, floors or limiters, the rule enforcing $E\mathbin{\cdot}B=0$ and $E^2\lt B^2$, current-sheet dissipation, and stopping time in rotations. It also includes software commit, compiler, precision, libraries, input files, and deterministic or stochastic choices. Outputs should include field snapshots, flux function, charge and current, open magnetic flux, Y-point location, and energy and angular-momentum fluxes on several radii.

No force-free executable or frozen input set is registered on this page, and no simulation has been rerun. A clean reproduction would provide a versioned environment and one command that regenerates benchmark time series and machine-readable diagnostics rather than only a final image.

Benchmarks, convergence, and uncertainty

Begin with a stationary monopole or split-monopole whose flux and luminosity are known, then recover the axisymmetric dipole solution and its light-cylinder smoothness. In three dimensions, compare aligned and oblique luminosities with Spitkovsky’s fit. As a separate vacuum benchmark, disable the complete force-free current closure, the $\boldsymbol E\mathbin{\cdot}\boldsymbol B$ projection, and the magnetic-dominance limiter before recovering the rotating vacuum dipole; setting only the plasma current to zero is not vacuum Maxwell evolution if either projection remains active. A single match in total luminosity does not validate current closure.

Use at least three resolutions and repeat with larger outer domains. Track $\boldsymbol\nabla\mathbin{\cdot}\boldsymbol B$, the volume and current-sheet distribution of $\boldsymbol E\mathbin{\cdot}\boldsymbol B$, the measure of $E^2\gt B^2$, energy balance, angular-momentum balance, open flux, Y-point position, and luminosity versus extraction radius. Current-sheet thickness often follows the grid; reconnection rate and dissipated fraction should not be called converged unless they approach a stable limit under both grid and scheme changes.

Separate truncation and outer-boundary uncertainty from closure uncertainty. Ideal force-free electrodynamics has no physical resistivity or pair injection scale, so different limiters can agree on large-scale luminosity while disagreeing on dissipation and high-energy emission.

Failure modes and interpretation

A numerical failure mode can be triggered by $E^2\gt B^2$ in the current sheet, inconsistent conductor boundary data, reflected waves, loss of magnetic divergence control, polar-axis artifacts, or excessive numerical resistivity. Another is false stationarity: the global luminosity can settle while the Y point or current sheet continues to drift. Run duration must be compared with the crossing time of the full domain.

Force-free validity breaks down in acceleration gaps, reconnection layers, and charge-starved regions where particle inertia and radiation matter. Such regions require kinetic or resistive models. A force-free current sheet can locate where dissipation is plausible, but it cannot predict a photon spectrum without additional particle physics.

Check your understanding

Using the force-free luminosity fit, what is the ratio of the orthogonal rotator luminosity at $\alpha=\pi/2$ to the aligned luminosity at $\alpha=0$? How does the aligned result differ from the vacuum point-dipole model?

Solution

The force-free fit gives

$$\frac{L_{\perp}}{L_{\parallel}} =\frac{1+\sin^2(\pi/2)}{1+\sin^2(0)}=2.$$

The aligned force-free luminosity is $L_{\parallel}=\mu^2\Omega^4$, whereas the vacuum point-dipole formula is proportional to $\sin^2\alpha$ and vanishes at alignment. Open flux and plasma current supply the aligned force-free torque.

References

• Contopoulos, Ioannis, Demosthenes Kazanas, and Christian Fendt. “The Axisymmetric Pulsar Magnetosphere.” The Astrophysical Journal 511 (1999): 351–358. DOI.
• Goldreich, Peter, and William H. Julian. “Pulsar Electrodynamics.” The Astrophysical Journal 157 (1969): 869–880. DOI.
• Spitkovsky, Anatoly. “Time-Dependent Force-Free Pulsar Magnetospheres: Axisymmetric and Oblique Rotators.” The Astrophysical Journal Letters 648 (2006): L51–L54. DOI.