Efficient cosmic ray acceleration, hydrodynamics, and Self-consistent Thermal X-ray Emission applied to SNR RX J1713.7-3946

We model the broad-band emission from SNR RX J1713.7-3946 including, for the first time, a consistent calculation of thermal X-ray emission together with non-thermal emission in a nonlinear diffusive shock acceleration (DSA) model. Our model tracks the evolution of the SNR including the plasma ionization state between the forward shock and the contact discontinuity. We use a plasma emissivity code to predict the thermal X-ray emission spectrum assuming the initially cold electrons are heated either by Coulomb collisions with the shock heated protons (the slowest possible heating), or come into instant equilibration with the protons. For either electron heating model, electrons reach >10^7 K rapidly and the X-ray line emission near 1 keV is more than 10 times as luminous as the underlying thermal bremsstrahlung continuum. Since recent Suzaku observations show no detectable line emission, this places strong constraints on the unshocked ambient medium density and on the relativistic electron to proton ratio. For the uniform circumstellar medium (CSM) models we consider, the low densities and high relativistic electron to proton ratios required to match the Suzaku X-ray observations definitively rule out pion-decay as the emission process producing GeV-TeV photons. We show that leptonic models, where inverse-Compton scattering against the cosmic background radiation dominates the GeV-TeV emission, produce better fits to the broad-band thermal and non-thermal observations in a uniform CSM.

💡 Research Summary

This paper presents a comprehensive model of the broadband emission from the supernova remnant (SNR) RX J1713.7‑3946 that, for the first time, self‑consistently incorporates thermal X‑ray line emission together with the non‑thermal radiation produced by nonlinear diffusive shock acceleration (DSA). The authors follow the dynamical evolution of the remnant from the explosion to the present, explicitly tracking the ionization state of the plasma between the forward shock and the contact discontinuity. A state‑of‑the‑art plasma emissivity code is employed to compute the thermal X‑ray spectrum under two extreme assumptions about electron heating: (1) slow heating by Coulomb collisions with shock‑heated protons, and (2) instantaneous equilibration with protons. In both scenarios electrons are rapidly heated to temperatures exceeding 10⁷ K, leading to strong line emission around 1 keV that is more than an order of magnitude brighter than the underlying thermal bremsstrahlung continuum.

Suzaku observations, however, reveal no detectable line features in the X‑ray band. This discrepancy forces the model to adopt a very low ambient density (≲0.01 cm⁻³) in order to suppress line emission. At such low densities the ratio of relativistic electrons to protons (K_ep) must be unusually high (≈10⁻²) to reproduce the observed synchrotron X‑ray flux. The combination of low density and high K_ep has a decisive implication for the origin of the GeV–TeV γ‑rays. In a uniform circumstellar medium (CSM) the pion‑decay channel, which requires a substantial density of target nuclei for proton–proton collisions, becomes inefficient; the predicted γ‑ray flux from hadronic interactions falls far short of the measurements.

Conversely, a leptonic scenario in which inverse‑Compton (IC) scattering of cosmic microwave background photons by the same relativistic electrons dominates the γ‑ray production provides an excellent fit to the entire spectral energy distribution. The IC component naturally accounts for the GeV–TeV flux without invoking high ambient densities, and it is compatible with the required high K_ep and low magnetic field (∼10 μG) inferred from the synchrotron modeling. The authors systematically explore the parameter space, demonstrating that the best‑fit models have ambient densities of order 10⁻³–10⁻² cm⁻³, electron‑to‑proton ratios near 10⁻², and magnetic fields modest enough to keep synchrotron losses subdominant to IC losses at TeV energies.

The study therefore concludes that, for RX J1713.7‑3946 embedded in a uniform CSM, the γ‑ray emission is most plausibly leptonic in origin, with IC scattering dominating over pion‑decay. This result challenges earlier interpretations that favored a hadronic origin and underscores the importance of incorporating self‑consistent thermal emission calculations when constraining particle acceleration efficiencies in young SNRs. The authors suggest that future high‑resolution X‑ray spectroscopy and deeper γ‑ray observations will be crucial for testing the low‑density, high‑K_ep scenario and for refining our understanding of cosmic‑ray acceleration in supernova remnants.