XMM-Newton Large Program on SN1006 - I: Methods and Initial Results of Spatially-Resolved Spectroscopy

TL;DR

This paper introduces advanced spatially-resolved X-ray spectroscopy methods applied to SN1006, mapping multiple physical parameters across the remnant and revealing complex ejecta and ISM structures.

Contribution

It develops new analysis techniques for detailed spectral mapping of supernova remnants, enabling comprehensive characterization of thermal, ionization, and abundance variations.

Findings

Mapped temperature, density, ionization, and metal abundances across SN1006.

Identified a low-abundance 'O Hole' with reverse shocked ejecta features.

Revealed asymmetric metal distributions suggesting explosion or ISM asymmetry.

Abstract

Based on our newly developed methods and the XMM-Newton large program of SN1006, we extract and analyze the spectra from 3596 tessellated regions of this SNR each with 0.3-8 keV counts $>10^4$. For the first time, we map out multiple physical parameters, such as the temperature ($kT$), electron density ($n_e$), ionization parameter ($n_et$), ionization age ($t_{ion}$), metal abundances, as well as the radio-to-X-ray slope ($\alpha$) and cutoff frequency ($\nu_{cutoff}$) of the synchrotron emission. We construct probability distribution functions of $kT$ and $n_et$, and model them with several Gaussians, in order to characterize the average thermal and ionization states of such an extended source. We construct equivalent width (EW) maps based on continuum interpolation with the spectral model of each regions. We then compare the EW maps of OVII, OVIII, OVII K$\delta-\zeta$, Ne, Mg, SiXIII, SiXIV, and S lines constructed with this method to those constructed with linear interpolation. We further extract spectra from larger regions to confirm the features revealed by parameter and EW maps, which are often not directly detectable on X-ray intensity images. For example, O abundance is consistent with solar across the SNR, except for a low-abundance hole in the center. This "O Hole" has enhanced OVII K$\delta-\zeta$ and Fe emissions, indicating recently reverse shocked ejecta, but also has the highest $n_et$, indicating forward shocked ISM. Therefore, a multi-temperature model is needed to decompose these components. The asymmetric metal distributions suggest there is either an asymmetric explosion of the SN or an asymmetric distribution of the ISM.