Cosmic Ray Electron Evolution in the Supernova Remnant RX J1713.7-3946

A simple formalism to describe nonthermal electron acceleration, evolution, and radiation in supernova remnants (SNRs) is presented. The electron continuity equation is analytically solved assuming that the nonthermal electron injection power is proportional to the rate at which the kinetic energy of matter swept up in an adiabatically expanding SNR shell. We apply this model to \fermi\ and HESS data from the SNR \rxj, and find that a one-zone leptonic model with Compton-scattered cosmic microwave background (CMB) and interstellar infrared photons has difficulty providing a good fit to its spectral energy distribution, provided the source is at a distance $\sim 1\ \kpc$ from the Earth. However, the inclusion of multiple zones, as hinted at by recent {\em Chandra} observations, does provide a good fit, but requires a second zone of compact knots with magnetic fields $B\sim 16\ \mu$G, comparable to shock-compressed fields found in the bulk of the remnant.

💡 Research Summary

The paper presents a concise yet comprehensive framework for describing non‑thermal electron acceleration, evolution, and radiation in supernova remnants (SNRs), with a specific application to RX J1713.7‑3946. Starting from the electron continuity equation, the authors assume that the power injected into relativistic electrons is a fixed fraction of the rate at which kinetic energy is transferred to the expanding shell as it sweeps up ambient material. This proportionality yields a time‑dependent injection term that can be integrated analytically together with loss processes—synchrotron cooling, inverse‑Compton scattering, and particle escape—producing a closed‑form solution for the electron energy distribution throughout both the free‑expansion and Sedov‑Taylor phases.

Using this electron population, the authors compute synchrotron emission (radio to X‑ray) and inverse‑Compton emission from scattering of the cosmic microwave background (CMB) and interstellar infrared (IR) photon fields. When the model is constrained to a single homogeneous zone (single magnetic field, uniform volume), it fails to simultaneously fit the GeV data from Fermi‑LAT and the TeV data from HESS. The discrepancy arises because the required electron spectrum that reproduces the TeV flux either overproduces the GeV flux or demands an unphysically low magnetic field, indicating that a simple one‑zone leptonic scenario is insufficient for a source at the commonly adopted distance of ~1 kpc.

Motivated by recent Chandra observations that reveal compact, bright knots within the remnant, the authors introduce a two‑zone model. The first, extended zone represents the bulk of the shell with a magnetic field of roughly 10 µG and a relatively modest electron injection rate. The second, compact zone corresponds to the knots, characterized by a higher magnetic field (~16 µG) and a significantly enhanced injection rate, reflecting shock‑compressed conditions. The knot zone dominates the inverse‑Compton output at TeV energies, while the extended zone supplies the bulk of the synchrotron X‑ray emission and the lower‑energy gamma‑rays. By adjusting the relative volumes and injection efficiencies, the combined spectrum reproduces the observed SED across the entire GeV–TeV range without invoking extreme parameters.

The study concludes that spatially inhomogeneous magnetic fields and electron acceleration efficiencies are essential to explain the high‑energy phenomenology of RX J1713.7‑3946. It also demonstrates that the analytic formalism, while simple, can be readily extended to multi‑zone configurations and applied to other young SNRs. The authors advocate for higher‑resolution X‑ray and gamma‑ray observations, as well as detailed three‑dimensional simulations, to further test the multi‑zone hypothesis and to refine our understanding of particle acceleration in shock‑dominated astrophysical environments.