Modeling the spectral evolution of PWNe inside SNRs

We present a new model for the spectral evolution of Pulsar Wind Nebulae inside Supernova Remnants. The model couples the long-term dynamics of these systems, as derived in the 1-D approximation, with a 1-zone description of the spectral evolution of the emitting plasma. Our goal is to provide a simplified theoretical description that can be used as a tool to put constraints on unknown properties of PWN-SNR systems: a piece of work that is preliminary to any more accurate and sophisticated modeling. In the present paper we apply the newly developed model to a few objects of different ages and luminosities. We find that an injection spectrum in the form of a broken-power law gives a satisfactory description of the emission for all the systems we consider. More surprisingly, we also find that the intrinsic spectral break turns out to be at a similar energy for all sources, in spite of the differences mentioned above. We discuss the implications of our findings on the workings of pulsar magnetospheres, pair multiplicity and on the particle acceleration mechanism(s) that might be at work at the pulsar wind termination shock.

💡 Research Summary

This paper presents a semi‑analytic framework for the long‑term evolution of pulsar wind nebulae (PWNe) embedded within supernova remnants (SNRs). The authors combine a one‑dimensional (1‑D) dynamical description of the PWN–SNR system, based on the Truelove & McKee (1999) model for SNR expansion, with a one‑zone (1‑zone) treatment of the particle distribution and radiative processes inside the nebula. The dynamical part assumes spherical symmetry, a stationary pulsar, and neglects radiative losses in the pressure balance, allowing analytic expressions for the nebular radius and internal pressure through the three evolutionary phases: free expansion, interaction with the SNR reverse shock, and subsequent reverberation/compression.

The spectral component treats the injection of electron‑positron pairs as a broken power‑law, (Q(E)\propto E^{-\alpha_1}) below a break energy (E_b) and (Q(E)\propto E^{-\alpha_2}) above it. The authors adopt typical indices (\alpha_1\approx1.5) and (\alpha_2\approx2.5) and determine (E_b) by fitting the full multi‑wavelength spectral energy distribution (SED) of each nebula. The particle evolution equation includes adiabatic expansion, synchrotron cooling, and inverse‑Compton scattering on the cosmic microwave background and local infrared/optical photon fields. The resulting particle spectrum is then used to compute synchrotron and IC emission, producing model SEDs that can be directly compared with observations.

Applying the model to four well‑studied PWNe—Crab, 3C 58, B1509‑58, and Kes 75—the authors find that a single broken‑power‑law injection reproduces the radio‑to‑γ‑ray data for all sources. Remarkably, the break energy (E_b) is found to lie in a narrow range (∼10⁴–10⁵ MeV) despite the wide spread in age, spin‑down power, and ambient density. This suggests a universal acceleration mechanism operating at the pulsar wind termination shock, largely independent of the pulsar’s global properties.

From the fitted normalisation the authors infer pair multiplicities (\kappa = \dot N / \dot N_{\rm GJ}) of order (10^4)–(10^5), far exceeding the Goldreich‑Julian minimum. Such high values challenge existing magnetospheric pair‑creation models (polar‑cap, slot‑gap, outer‑gap), which typically predict lower rates for the observed spin parameters. The paper discusses possible resolutions, including additional low‑energy particle sources (thermal pick‑up, early‑time non‑relativistic acceleration) and the role of a heavy‑ion component that could drive cyclotron‑resonant acceleration, but argues that the observed smooth spectral continuity favours a single, pulsar‑originated injection spectrum.

The dynamical fit also yields wind Lorentz factors (\gamma_w\sim10^6) and magnetisation parameters (\sigma_w\ll1) (∼10⁻³–10⁻²), consistent with earlier 1‑D MHD models (Kennel & Coroniti 1984) and recent multidimensional simulations that find the nebular flow to be particle‑dominated downstream of the termination shock.

Overall, the study demonstrates that a coupled 1‑D dynamical + 1‑zone spectral model can capture the essential physics of PWN evolution while remaining computationally inexpensive. This makes it a valuable tool for exploring the large parameter space of older, fainter PWNe where full 3‑D magnetohydrodynamic simulations are impractical. The authors conclude that the apparent universality of the broken‑power‑law injection and the high pair multiplicities provide strong constraints on pulsar magnetosphere physics and on the microphysics of particle acceleration at relativistic termination shocks.