Time Evolution of the Reverse Shock in SN 1006

The Schweizer-Middleditch star, located behind the SN 1006 remnant and near its center in projection, provides the opportunity to study cold, expanding ejecta within the SN 1006 shell through UV absorption. Especially notable is an extremely sharp red edge to the Si II 1260 Angstrom feature, which stems from the fastest moving ejecta on the far side of the SN 1006 shell–material that is just encountering the reverse shock. Comparing HST far-UV spectra obtained with COS in 2010 and with STIS in 1999, we have measured the change in this feature over the intervening 10.5-year baseline. We find that the sharp red edge of the Si II feature has shifted blueward by 0.19 +/- 0.05 Angstroms, which means that the material hitting the reverse shock in 2010 was moving slower by 44 +/- 11 km/s than the material that was hitting it in 1999, a change corresponding to - 4.2 +/- 1.0 km/s/yr. This is the first observational confirmation of a long-predicted dynamic effect for a reverse shock: that the shock will work its way inward through expanding supernova ejecta and encounter ever slower material as it proceeds. We also find that the column density of shocked Si II (material that has passed through the reverse shock) has decreased by 7 +/- 2% over the ten-year period. The decrease could indicate that in this direction the reverse shock has been ploughing through a dense clump of Si,leading to pressure and density transients.

💡 Research Summary

The paper presents a direct observational test of a fundamental prediction of the reverse‑shock model in young supernova remnants: as the reverse shock propagates inward through freely expanding ejecta, it should encounter material with progressively lower expansion velocities. The authors exploit the Schweizer‑Middleditch (SM) star, an OB sub‑dwarf located just behind the SN 1006 shell, as a bright ultraviolet background source. By comparing far‑UV spectra of the SM star taken with HST/STIS in 1999 and with HST/COS in 2010, they track the evolution of the broad Si II λ1260 Å absorption feature that originates in the supernova ejecta.

The STIS observation (1150–1700 Å, resolution ≈ 4.6 km s⁻¹) and the COS observation (≈ 1170–1470 Å, resolution ≈ 18 km s⁻¹) were carefully wavelength‑registered. The authors used numerous narrow interstellar and stellar lines near the Si II edge to cross‑correlate the two data sets, correcting a systematic offset of ~40 mÅ (≈ 9.6 km s⁻¹) with a linear transformation and achieving a final registration accuracy of ±6 mÅ (±1.5 km s⁻¹).

The Si II 1260 Å line exhibits an extremely sharp red edge at +7026 km s⁻¹, interpreted as the velocity of the unshocked ejecta just before it encounters the reverse shock. Between 1999 and 2010 this edge shifted blueward by 0.19 ± 0.05 Å, corresponding to a decrease in the line‑of‑sight velocity of 44 ± 11 km s⁻¹, i.e. a deceleration of –4.2 ± 1.0 km s⁻¹ yr⁻¹. This is the first quantitative, spectroscopic confirmation that the reverse shock is indeed moving into slower‑moving ejecta.

In addition to the edge shift, the column density of shocked Si II (the component that has passed through the reverse shock) declined by 7 ± 2 % over the same interval. The authors interpret this as evidence that, along this particular line of sight, the reverse shock has been traversing a dense Si clump, producing transient pressure and density variations that affect the post‑shock Si II column.

To extract physical parameters, the authors model the Si II profile as a sum of an unshocked component (characterized by the free‑expansion velocity rₛ/t, where rₛ is the reverse‑shock radius and t the remnant age) and a shocked component (described by a mean velocity v and a one‑dimensional velocity dispersion σ). They enforce the shock‑jump condition Δv = √3 σ, which holds when virtually all shock energy goes into ion heating, with negligible electron heating or cosmic‑ray acceleration. The fitted parameters satisfy this condition, supporting the view that SN 1006’s reverse shock is ion‑dominated.

The measured shocked Si II column density is Nₛₕₖ ≈ 7.4 × 10¹⁴ cm⁻². Using the updated collisional ionization rate for Si II (σᵢ ≈ 4.4 × 10⁻⁸ cm³ s⁻¹ from Clark & Abdallah 2003), the authors compute the expected steady‑state column density and find good agreement with the observation. This resolves earlier discrepancies noted by Hamilton et al. (1997) that arose from older ionization rates.

Overall, the work provides the first spectroscopic, time‑domain verification of reverse‑shock dynamics in a Type Ia supernova remnant. It demonstrates that high‑resolution UV absorption spectroscopy of background sources can probe the interior kinematics of SNRs, complementing imaging studies such as those of Cas A. The observed decrease in shocked Si II column density also highlights that the ejecta are not homogeneous; small‑scale density or compositional clumps can modulate the shock’s progress.

Future work suggested by the authors includes expanding the sample of background UV sources (e.g., additional quasars behind SN 1006) to map the three‑dimensional structure of the reverse shock, and obtaining higher‑resolution, higher‑signal‑to‑noise UV spectra to track other ionic species (Fe II, Si IV, etc.). Such observations would enable a more complete picture of electron heating, ionization balance, and possible particle acceleration at the reverse shock, deepening our understanding of how young supernova remnants evolve and mix their nucleosynthetic products into the interstellar medium.