[Paper Review] Magnetohydrodynamic Models for the Structure of Pulsar-Wind Nebulae

Stephen P. Reynolds|arXiv (Cornell University)|Aug 27, 2003

Pulsars and Gravitational Waves Research3 citations

TL;DR

This paper extends the steady-state magnetohydrodynamic (MHD) model of Kennel & Coroniti to explore its applicability to pulsar-wind nebulae (PWNe) beyond the Crab Nebula. Using spherical symmetry, a purely toroidal magnetic field, and steady-state dynamics, it models electron energy evolution and synchrotron emissivity across varying $σ$ (magnetic-to-particle energy flux ratio), finding that predicted X-ray photon index variations with radius and nebular size evolution conflict with observations, indicating the need for non-steady or non-spherical flows and electron transport beyond convection.

ABSTRACT

The Crab Nebula is well-described at optical wavelengths and above by a steady-state magnetohydrodynamic model due to Kennel and Coroniti (1984). Can this class of model describe other pulsar-wind nebulae? I exhibit simple generalizations of KC models,for various values of sigma, the ratio of magnetic to particle flux input at the wind shock. I calculate the evolution of the electron spectrum and synchrotron emissivity in the nebula assuming spherical symmetry, a steady state, and a purely toroidal magnetic field. Emission profiles depend on the initial magnetic field B_0, the electron spectral index s, and the angle of the toroidal axis with the line of sight, phi. I show integrated spectra and radial profiles for various cases, along with predicted variations of photon index Gamma with radius in X-rays. Most models predict much smaller sizes in X-rays, and curves of Gamma(r) which do not resemble observations. Further elaborations of the dynamics of PWNe seem necessary.

Motivation & Objective

• To assess whether the steady-state MHD model of Kennel & Coroniti (KC84) can describe other pulsar-wind nebulae beyond the Crab Nebula.
• To investigate how variations in $σ$, initial magnetic field $B_0$, electron spectral index $s$, and viewing angle $φ$ affect synchrotron emission profiles and spectra.
• To determine whether the model's predictions for nebular size and photon index evolution with radius match X-ray observations of non-Crab PWNe.
• To identify limitations of the simple MHD framework in explaining observed PWN morphologies and spectral behaviors, especially in fainter, larger nebulae like 3C 58.

Proposed method

• Adopts a steady-state, spherically symmetric MHD model with a purely toroidal magnetic field, generalizing the KC84 model for varying $σ$.
• Uses the velocity law $u(r) = \frac{1}{r^2}(1 - u_\infty) + u_\infty$, where $u_\infty = \sigma / (1 + \sigma)$, to describe postshock flow dynamics.
• Applies the electron energy evolution equation $E(t) = \frac{E_0 \alpha^{1/3}}{1 + E_0 (a B_0^2 r_0 / v_0) \int z^{-8/3} u^{-10/3} dz}$ to track radiative and adiabatic losses.
• Calculates the electron distribution $N(E) = N(E_0) \left( \frac{E_0^2}{E^2} \right) \left( \frac{\rho}{\rho_0} \right)^{4/3}$ to determine synchrotron emissivity at each radius.
• Numerically integrates emissivity over solid angle to generate brightness profiles for different aspect angles $\phi$, with $\phi = 0^\circ$ indicating the magnetic field plane-on to the line of sight.
• Computes integrated spectra and radial profiles for various $\sigma$, $s$, $B_0$, and $\phi$, and evaluates predicted photon index $\Gamma(r)$ in X-rays.

Experimental results

Research questions

• RQ1 Can the steady-state MHD model of Kennel & Coroniti describe the structure and emission of pulsar-wind nebulae other than the Crab Nebula?
• RQ2 How do variations in $σ$, electron spectral index $s$, and viewing angle $\phi$ affect the predicted X-ray brightness profiles and spectral indices of PWNe?
• RQ3 Do the model predictions for the radial dependence of the X-ray photon index $\Gamma(r)$ match observed trends in PWNe like G21.5-0.9?
• RQ4 Is the predicted decrease in nebular size from radio to X-ray frequencies consistent with observations of PWNe such as 3C 58, which show extended X-ray emission matching radio extents?
• RQ5 Do the model-predicted HWHM radii for X-ray emission fall within observed ranges, or are they too small?

Key findings

• All MHD models predict a significant decrease in nebular size between radio and X-ray frequencies, with size dropping by at least a factor of 2, contradicting observations of PWNe like 3C 58 where X-ray emission extends to radio edges.
• The models predict X-ray brightness profiles with HWHM radii only 3–6 times the injection radius, which may be too small compared to observed nebulae.
• The predicted radial dependence of the X-ray photon index $\Gamma(r)$ is not linear with radius, in contrast to observations of PWNe such as G21.5-0.9, which show a roughly linear $\Gamma(r)$ trend.
• For flat-injected spectra ($s \sim 1$), the models produce a distinct energy-loss 'bump' in the spectrum, causing concave-up curvature over a broad frequency range.
• The aspect angle $\phi$ affects flux normalization but has minimal impact on the shape of the integrated spectrum and only minor effects on profile shape at energies where radiative losses are significant.
• The models fail to reproduce the observed lack of strong size evolution and linear $\Gamma(r)$ trend in X-rays, indicating that non-steady or non-spherical flows and/or electron diffusion are likely necessary to explain the data.