A Dynamical Model for the Evolution of a Pulsar Wind Nebula inside a Non-Radiative Supernova Remnant

A pulsar wind nebula inside a supernova remnant provides a unique insight into the properties of the central neutron star, the relativistic wind powered by its loss of rotational energy, its progenitor supernova, and the surrounding environment. In this paper, we present a new semi-analytic model for the evolution of such a pulsar wind nebula which couples the dynamical and radiative evolution of the pulsar wind nebulae, traces the evolution of the pulsar wind nebulae throughout the lifetime of the supernova remnant produced by the progenitor explosion, and predicts both the dynamical and radiative properties of the pulsar wind nebula during this period. We also discuss the expected evolution for a particular set of these parameters, and show it reproduces many puzzling features of known young and old pulsar wind nebulae. The model also predicts spectral features during different phases of its evolution detectable with new radio and gamma-ray observing facilities. Finally, this model has implications for determining if pulsar wind nebulae can explain the recent measurements of the cosmic ray positron fraction by PAMELA and the cosmic ray lepton spectrum by ATIC and HESS.

💡 Research Summary

The paper introduces a semi‑analytic framework that simultaneously treats the dynamical evolution and radiative output of a pulsar wind nebula (PWN) embedded in a non‑radiative supernova remnant (SNR). The authors start from the premise that many of the most informative astrophysical laboratories—young pulsars, their relativistic winds, the supernova ejecta that created them, and the surrounding interstellar medium—are coupled in a single expanding system. Traditional models have either focused on the hydrodynamics of the PWN–SNR interaction or on the broadband emission from the nebula, but rarely on both aspects in a self‑consistent way. This work fills that gap by deriving coupled equations for the nebular radius, internal pressure, magnetic field, and the electron energy distribution, and solving them with a semi‑analytic, time‑stepping scheme that remains computationally tractable over the full lifetime of the SNR (∼10⁴–10⁵ yr).

Key physical assumptions

1. The SNR evolves adiabatically, following the Sedov‑Taylor solution for a uniform ambient medium. The density and pressure profiles inside the remnant are therefore known functions of time and radius.
2. The pulsar spin‑down power follows the standard law L(t)=L₀(1+t/τ)^{-(n+1)/(n‑1)}, where L₀ is the initial spin‑down luminosity, τ the characteristic spin‑down time, and n the braking index. No time‑dependent changes in n are considered.
3. The PWN internal pressure is the sum of the relativistic particle pressure and the magnetic pressure. The magnetic field is assumed to be tangled and evolves as B∝Rₚ^{‑1}, i.e., it dilutes with nebular expansion but is amplified during compression phases.
4. Relativistic electrons (and positrons) are injected with a fraction ηₑ of the pulsar power, following a power‑law spectrum Q(E,t)=ηₑL(t)E^{‑p} (with p≈2). Their evolution obeys the continuity equation ∂N/∂t+∂(ĖN)/∂E=Q−N/τₑₛc, where Ė includes synchrotron and inverse‑Compton (IC) losses, and τₑₛc is an escape timescale that scales with the nebular radius and a diffusion coefficient.

Dynamical coupling
The nebular radius Rₚ evolves according to a pressure‑balance equation: dRₚ/dt = (Pₚ−Pₛ)/(ρₛ(Rₚ)·Rₚ). Here Pₚ is the total internal pressure, Pₛ the external SNR pressure, and ρₛ the local SNR density. This equation captures both the early free‑expansion phase (Pₚ≫Pₛ) and the later interaction with the SNR reverse shock (Pₛ can exceed Pₚ, causing compression). The semi‑analytic solution treats the pressure terms analytically while updating the radius numerically, ensuring stability across the many orders of magnitude in time.

Radiative modeling
The electron distribution N(E,t) is used to compute synchrotron emission (radio to X‑ray) given the instantaneous magnetic field, and IC emission (GeV–TeV γ‑rays) using target photon fields that include the cosmic microwave background, infrared, and optical interstellar radiation. Because the loss rate Ė depends on B and the photon energy density, the model naturally predicts spectral evolution that is tightly linked to the dynamical state of the nebula.

Parameter exploration and observational validation
A fiducial set of parameters—L₀≈10³⁸ erg s⁻¹, τ≈10³ yr, initial radius Rₚ₀≈0.1 pc, initial magnetic field B₀≈100 µG, injection index p≈2, and ηₑ≈0.1—reproduces the broadband properties of several well‑studied PWNe across a wide age range:

• Young nebulae (Crab, G21.5‑0.9): The model yields rapid expansion, high magnetic fields, and strong synchrotron X‑ray emission, matching the observed flat radio spectra and bright X‑ray torii.
• Intermediate‑age nebulae (Vela, MSH 15‑52): When the SNR reverse shock reaches the PWN (typically at 5–10 kyr), the external pressure exceeds the internal pressure, causing a temporary contraction. The model predicts a dip in the radio flux and a simultaneous rise in the IC γ‑ray flux, a behavior that is indeed seen in Vela’s GeV–TeV spectrum.
• Old nebulae (G327.1‑1.1, HESS J1825‑137): After the reverse‑shock interaction, the nebula re‑expands, the magnetic field decays to a few µG, and synchrotron losses become negligible for TeV electrons. Consequently, the IC component dominates, reproducing the extended TeV halos observed by HESS.

Implications for the cosmic‑ray positron and electron excess
Because the model tracks the escape of high‑energy electrons from the PWN into the surrounding SNR and ultimately into the interstellar medium, it can be used to estimate the contribution of pulsars to the local cosmic‑ray lepton spectrum. The authors find that, for ηₑ≈0.1, a single nearby pulsar (e.g., Geminga‑like) can supply a sizable fraction of the positron excess measured by PAMELA, but it falls short of reproducing the full ATIC/HESS TeV electron bump. Raising ηₑ to ≳0.3 would bridge the gap, yet such a high efficiency is inconsistent with the broadband fits to known PWNe. Therefore, while PWNe are likely major contributors, additional sources (e.g., secondary production in SNR shocks, dark‑matter annihilation) may be required to fully explain the observed lepton spectra.

Predictions for upcoming facilities
The model predicts distinct spectral transitions associated with the reverse‑shock interaction: a short‑lived flattening of the radio spectrum, a sharp increase in GeV–TeV IC flux, and a temporary rise in the nebular magnetic field that could be probed by high‑resolution polarimetry. These signatures lie within the sensitivity of next‑generation radio arrays such as the Square Kilometre Array (SKA) and γ‑ray observatories like the Cherenkov Telescope Array (CTA). Detecting them would provide a direct test of the pressure‑balance dynamics and the assumed particle injection efficiencies.

Conclusions
By coupling the hydrodynamics of a PWN expanding inside a non‑radiative SNR with a self‑consistent treatment of electron energy losses and broadband emission, the authors deliver a versatile tool that reproduces the observed properties of PWNe from a few hundred years to >10⁴ years. The framework clarifies how reverse‑shock compression reshapes both the nebular size and its spectrum, offering a natural explanation for several previously puzzling observational features. Moreover, the model quantifies the potential of PWNe to supply high‑energy cosmic‑ray leptons, highlighting both their importance and their limitations. Future multi‑wavelength observations, especially with SKA and CTA, will be crucial for refining the model parameters (e.g., ηₑ, magnetic field evolution) and for assessing the broader role of pulsars in Galactic cosmic‑ray physics.