Type Ia Supernova Remnants Research Articles

Chandrasekhar and Sub-Chandrasekhar Models for the X-Ray Emission of Type Ia Supernova Remnants. I. Bulk Properties

Abstract Type Ia supernovae originate from the explosion of carbon–oxygen white dwarfs in binary systems, but the exact nature of their progenitors remains elusive. The bulk properties of Type Ia supernova remnants, such as the radius and the centroid energy of the Fe Kα blend in the X-ray spectrum, are determined by the properties of the supernova ejecta and the ambient medium. We model the interaction between Chandrasekhar and sub-Chandrasekhar models for Type Ia supernova ejecta and a range of uniform ambient medium densities in one dimension up to an age of 5000 years. We generate synthetic X-ray spectra from these supernova remnant models and compare their bulk properties at different expansion ages with X-ray observations from Chandra and Suzaku. We find that our models can successfully reproduce the bulk properties of most observed remnants, suggesting that Type Ia SN progenitors do not modify their surroundings significantly on scales of a few pc, although more detailed models are required to establish quantitative limits on the density of any such surrounding circumstellar material. Ambient medium density and expansion age are the main contributors to the diversity of the bulk properties in our models. Chandrasekhar and sub-Chandrasekhar progenitors make similar predictions for the bulk remnant properties, but detailed fits to X-ray spectra have the power to discriminate explosion energetics and progenitor scenarios.

Constraints on Cosmic-ray Acceleration Efficiency in Balmer Shocks of Two Young Type Ia Supernova Remnants in the Large Magellanic Cloud

Abstract We present results from an optical study of two young Balmer-dominated remnants of SNIa in the Large Magellanic Cloud, 0509−67.5 and 0519−69.0, in an attempt to search for signatures of efficient cosmic-ray (CR) acceleration. We combine proper motion measurements with corresponding optical spectroscopic measurements of the Hα line at multiple rim positions from VLT/FORS2 and SALT/RSS and compare our results to published Balmer shock models. Analyses of the optical spectra result in broad Hα widths in the range of 1800–4000 km s−1 for 12 separate Balmer-dominated filaments that show no evidence for forbidden line emission; the corresponding shock speeds from proper motion measurements from HST span a range of 1700–8500 km s−1. Our measured values of shock speeds and broad Hα widths in 0509−67.5 and 0519−69.0 are fit well with a Balmer shock model that does not include effects of efficient CR acceleration. We determine an upper limit of 7%/χ (95% confidence) on the CR acceleration efficiency for our ensemble of data points, where χ is the ionization fraction of the pre-shock gas. The upper limits on the individual remnants are 6%/χ (0509−67.5) and 11%/χ (0519−69.0). These upper limits are below the integrated CR acceleration efficiency in the Tycho supernova remnant, where the shocks predominantly show little Hα emission, indicating that Balmer-dominated shocks are not efficient CR accelerators.

Asymmetric Type-Ia supernova origin of W49B as revealed from spatially resolved X-ray spectroscopic study

The origin of the asymmetric supernova remnant (SNR) W49B has been a matter of debate: is it produced by a rare jet-driven core-collapse (CC) supernova, or by a normal supernova that is strongly shaped by its dense environment? Aiming to uncover the explosion mechanism and origin of the asymmetric, centrally filled X-ray morphology of W49B, we have performed spatially resolved X-ray spectroscopy and a search for potential point sources. We report new candidate point sources inside W49B. The Chandra X-ray spectra from W49B are well-characterized by two-temperature gas components (~0.27 keV + 0.6–2.2 keV). The hot component gas shows a large temperature gradient from the northeast to the southwest and is over-ionized in most regions with recombination timescales of 1–10 × 1011 cm−3 s. The Fe element shows strong lateral distribution in the SNR east, while the distribution of Si, S, Ar, Ca is relatively smooth and nearly axially symmetric. Asymmetric Type-Ia explosion of a Chandrasekhar-mass white dwarf (WD) well-explains the abundance ratios and metal distribution of W49B, whereas a jet-driven explosion and normal CC models fail to describe the abundance ratios and large masses of iron-group elements. A model based on a multi-spot ignition of the WD can explain the observed high MMn∕MCr value (0.8–2.2). The bar-like morphology is mainly due to a density enhancement in the center, given the good spatial correlation between gas density and X-ray brightness. The recombination ages and the Sedov age consistently suggest a revised SNR age of 5–6 kyr. This study suggests that despite the presence of candidate point sources projected within the boundary of this SNR, W49B is likely a Type-Ia SNR, which suggests that Type-Ia supernovae can also result in mixed-morphology SNRs.

A search for a surviving companion in SN 1006

Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surveys, however, have only looked for stars as faint as $\approx 0.1 L_\odot$ and thus would have missed a surviving white dwarf companion. In this work, we present very deep DECAM imaging $(u, g, r, z)$ of the Type Ia supernova remnant SN 1006 specifically to search for a potential surviving white dwarf companion. We find no object within the inner third of the SN 1006 remnant that is consistent with a relatively young cooling white dwarf. We find that if there is a companion white dwarf, it must have formed long ago and cooled undisturbed for $> 10^8$ yr to be consistent with the redder objects in our sample. We conclude that our findings are consistent with the complete destruction of the secondary (such as in a merger) or an unexpectedly cool and thus very dim surviving companion white dwarf.

Radio Evolution of Supernova Remnants Including Nonlinear Particle Acceleration: Insights from Hydrodynamic Simulations

Abstract We present a model for the radio evolution of supernova remnants (SNRs) obtained by using three-dimensional hydrodynamic simulations coupled with nonlinear kinetic theory of cosmic-ray (CR) acceleration in SNRs. We model the radio evolution of SNRs on a global level by performing simulations for a wide range of the relevant physical parameters, such as the ambient density, supernova (SN) explosion energy, acceleration efficiency, and magnetic field amplification (MFA) efficiency. We attribute the observed spread of radio surface brightnesses for corresponding SNR diameters to the spread of these parameters. In addition to our simulations of Type Ia SNRs, we also considered SNR radio evolution in denser, nonuniform circumstellar environments modified by the progenitor star wind. These simulations start with the mass of the ejecta substantially higher than in the case of a Type Ia SN and presumably lower shock speed. The magnetic field is understandably seen as very important for the radio evolution of SNRs. In terms of MFA, we include both resonant and nonresonant modes in our large-scale simulations by implementing models obtained from first-principles, particle-in-cell simulations and nonlinear magnetohydrodynamical simulations. We test the quality and reliability of our models on a sample consisting of Galactic and extragalactic SNRs. Our simulations give Σ − D slopes between −4 and −6 for the full Sedov regime. Recent empirical slopes obtained for the Galactic samples are around −5, while those for the extragalactic samples are around −4.

The surviving companions in type Ia supernova remnants

The single-degenerate (SD) model is one of the most popular progenitor models of type Ia supernovae (SNe Ia), in which the companion star can survive after an SN Ia explosion and show peculiar properties. Therefore, searching for the surviving companion in type Ia supernova remnants (SNRs) is a potential method to verify the SD model. In the SN 1604 remnant (Kepler’s SNR), although Chandra X-ray observation suggests that the progenitor is most likely a WD+AGB system, a the surviving companion has not been found. One possible reason is rapid rotation of the white dwarf (WD), causing explosion of the WD to be delayed for a spin-down timescale, and then the companion evolved into a WD before the supernova explosion, so the companion is too dim to be detected. We aim to verify this possible explanation by carrying out binary evolution calculations. In this paper, we use Eggleton’s stellar evolution code to calculate the evolution of binaries consisting of a WD+red giant (RG). We assume that the rapidly rotating WD can continuously increase its mass when its mass exceeds the Chandrasekhar mass limit () until the mass-transfer rate decreases to be lower than a critical value. Eventually, we obtain the final masses of a WD in the range 1.378 M⊙ to 2.707 M⊙. We also show that if the spin-down time is less than 106 yr, the companion star will be very bright and easily observed; but if the spin-down time is as long as ∼ 107 yr, the luminosities of the surviving companion would be lower than the detection limit. Our simulation provides guidance in hunting for the surviving companion stars in SNRs, and the fact that no surviving companion has been found in Kepler’s SNR may not be definite evidence disfavoring the SD origin of Kepler’s SN.

Observational Evidence for High Neutronization in Supernova Remnants: Implications for Type Ia Supernova Progenitors

Abstract The physical process whereby a carbon–oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using Ca and S lines in the X-ray spectra of Type Ia supernova remnants (SNRs). Applying this method to Suzaku data of Tycho, Kepler, 3C 397, and G337.2−0.7 in the Milky Way, and N103B in the Large Magellanic Cloud, we find that the neutronization of the ejecta in N103B is comparable to that of Tycho and Kepler, which suggests that progenitor metallicity is not the only source of neutronization in SNe Ia. We then use a grid of SN Ia explosion models to infer the metallicities of the stellar progenitors of our SNRs. The implied metallicities of 3C 397, G337.2−0.7, and N103B are major outliers compared to the local stellar metallicity distribution functions, indicating that progenitor metallicity can be ruled out as the origin of neutronization for these SNRs. Although the relationship between ejecta neutronization and equivalent progenitor metallicity is subject to uncertainties stemming from the 12C + 16O reaction rate, which affects the Ca/S mass ratio, our main results are not sensitive to these details.

Physical Structures of the Type Ia Supernova Remnant N103B

N103B is a Type Ia supernova remnant (SNR) projected in the outskirt of the superbubble around the rich cluster NGC 1850 in the Large Magellanic Cloud (LMC). We have obtained H$\alpha$ and continuum images of N103B with the $\textit{Hubble Space Telescope}$ ($\textit{HST}$) and high-dispersion spectra with 4m and 1.5m telescopes at Cerro Tololo Inter-American Observatory. The $\textit{HST}$ H$\alpha$ image exhibits a complex system of nebular knots inside an incomplete filamentary elliptical shell that opens to the east where X-ray and radio emission extends further out. Electron densities of the nebular knots, determined from the [S II] doublet, reach 5300 cm$^{-3}$, indicating an origin of circumstellar medium, rather than interstellar medium. The high-dispersion spectra reveal three kinematic components in N103B: (1) a narrow component with [N II]6583/H$\alpha$ $\sim$ 0.14 from the ionized interstellar gas associated with the superbubble of NGC 1850 in the background, (2) a broader H$\alpha$ component with no [N II] counterpart from the SNR's collisionless shocks into a mostly neutral ambient medium, and (3) a broad component, $\Delta V$ $\sim$ 500 km s$^{-1}$, in both H$\alpha$ and [N II] lines from shocked material in the nebular knots. The Balmer-dominated filaments can be fitted by an ellipse, and we adopt its center as the site of SN explosion. We find that the star closest to this explosion center has colors and luminosity consistent with a 1 $M_\odot$ surviving subgiant companion as modelled by Podsiadlowski. Follow-up spectroscopic observations are needed to confirm this star as the SN's surviving companion.

Gamma-ray Observations of Tycho’s Supernova Remnant with VERITAS and Fermi

Abstract High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho’s SNR is a particularly good target because it is a young, type Ia SNR that has been well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho’s SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho’s SNR with 147 hr of VERITAS and 84 months of Fermi-LAT observations, which represent about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is 2.92 ± 0.42stat ± 0.20sys. It is also softer than the spectral index in the GeV energy range, 2.14 ± 0.09stat ± 0.02sys, measured in this study using Fermi-LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi-LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously.

Searching for Young Type Ia Supernova Remnants in M83

AbstractYoung Type Ia supernova remnants (SNRs) are characterized by Balmer-dominated optical spectra, well-defined shell morphologies, > 1036 ergs s−1 X-ray luminosities, and a lack of massive stars and dense interstellar gas in their vicinity. Applying these characteristics and using archival deep HST and Chandra observations of M83, we search for young Type Ia SNRs in this spiral galaxy. This is a very difficult task!

THE ORIGIN OF THE IRON-RICH KNOT IN TYCHO’S SUPERNOVA REMNANT

ABSTRACT X-ray observations of supernova remnants (SNRs) allow us to investigate the chemical inhomogeneity of ejecta, offering unique insight into the nucleosynthesis in supernova explosions. Here we present detailed imaging and spectroscopic studies of the “Fe knot” located along the eastern rim of the Type Ia SNR Tycho (SN 1572) using Suzaku and Chandra long-exposure data. Surprisingly, the Suzaku spectrum of this knot shows no emission from Cr, Mn, or Ni, which is unusual for the Fe-rich regions in this SNR. Within the framework of the canonical delayed-detonation models for SN Ia, the observed mass ratios , , and (at 90% confidence) can only be achieved for a peak temperature of (5.3–5.7) K and a neutron excess of . These constraints rule out the deep, dense core of a Chandrasekhar-mass white dwarf as the origin of the Fe knot and favor either incomplete Si burning or an α-rich freeze-out regime, probably close to the boundary. An explosive He burning regime is a possible alternative, although this hypothesis is in conflict with the main properties of this SNR.

EXPANDING MOLECULAR BUBBLE SURROUNDING TYCHO’S SUPERNOVA REMNANT (SN 1572) OBSERVED WITH THE IRAM 30 m TELESCOPE: EVIDENCE FOR A SINGLE-DEGENERATE PROGENITOR

ABSTRACT Whether the progenitors of SNe Ia are single-degenerate or double-degenerate white dwarf (WD) systems is a highly debated topic. To address the origin of Tycho’s Type Ia supernova remnant (SNR), SN 1572, we have carried out a 12CO J = 2–1 mapping and a 3 mm line survey toward the remnant using the IRAM 30 m telescope. We show that Tycho is surrounded by a clumpy molecular bubble at a local standard of rest velocity of , which expands at a speed of and has a mass of (at the distance of 2.5 kpc). Enhanced 12CO J = 2–1 line emission relative to 12CO J = 1–0 emission and possible line broadenings (in velocity range −64 to −60 km s−1) are found at the northeastern boundary of the SNR, where the shell is deformed and decelerated. These features, combined with the morphological correspondence between the expanding molecular bubble and Tycho, suggest that the SNR is associated with the bubble at the velocity range −66 to −57 km s−1. The most plausible origin for the expanding bubble is the fast outflow (with velocity of hundreds km s−1) driven from the vicinity of a WD as it accreted matter from a nondegenerate companion star. The SNR has been expanding in the low-density wind-blown bubble, and the shock wave has just reached the molecular cavity wall. This is the first unambiguous detection of an expanding bubble driven by the progenitor of a Type Ia SNR, which constitutes evidence for a single-degenerate progenitor for this SN Ia.

On the brightness of surviving companions in type Ia supernova remnants

Abstract The progenitor systems for type Ia supernovae are still controversial. One of the methods used to test the proposed scenario for the progenitor systems is to identify companions that are supposed to survive according to the so-called single degenerate scenario. These companions might be affected by supernova ejecta. We present several numerical simulations of surviving red giant companions whose envelopes were stripped and heated. We find that red giants with less-massive helium cores (≲0.30 M⊙) can be so faint after the supernovae that we cannot detect them. In addition, we apply the results to the case of SNR 0509−67.5, and put constraints on the helium core mass, envelope stripping, and energy injection under the single degenerate scenario for type Ia supernovae.

Suzaku study on the ejecta of the supernova remnant G272.2−3.2

Abstract We report reanalyses of the Suzaku observations of the Galactic supernova remnant (SNR) G272.2−3.2, for which previous studies were limited below 3 keV. With careful data reduction and background subtraction, we discover the K-shell lines of Ar, Ca, and Fe above 3 keV. The X-ray spectrum of G272.2−3.2 consists of two components, a low-temperature collisional ionization equilibrium (CIE) plasma (kTe ∼ 0.2 keV) and a high-temperature non-equilibrium ionization (NEI) plasma (kTe = 0.6–3 keV). The CIE plasma has solar abundances over the entire area, hence it originates from the interstellar medium. On the other hand, the abundances of the NEI plasma increase toward the inner region, suggesting ejecta origin. The line center energy of the Fe K-shell emission (∼6.4 keV) suggests that the ejecta are recently heated by reverse shock, a common feature in Type Ia SNRs.

Type Ia supernova remnants: shaping by iron bullets

Using 2D numerical hydrodynamical simulations of Type Ia supernova remnants (SNR Ia) we show that iron clumps few times denser than the rest of the SN ejecta might form protrusions in an otherwise spherical SNR. Such protrusions exist in some SNR Ia, e.g. SNR 1885 and Tycho. Iron clumps are expected to form in the deflagration to detonation explosion model. In SNR Ia where there are two opposite protrusions, termed ‘ears’, such as Kepler's SNR and SNR G1.9+0.3, our scenario implies that the dense clumps, or iron bullets, were formed along an axis. Such a preferred axis can result from a rotating white dwarf progenitor. If our claim holds, this offers an important clue to the SN Ia explosion scenario.

X-RAY SPECTROSCOPY OF POTENTIAL SMALL MAGELLANIC CLOUD TYPE Ia SUPERNOVA REMNANTS AND THEIR ENVIRONMENTS

We examine three supernova remnants in the SMC, IKT 5 (supernova remnant (SNR) 0047-73.5), IKT 25 (SNR 0104-72.3), and DEM S 128 (SNR 0103-72.4), which are designated as Type Ia in the literature due to their spectra and morphology. This is troublesome because of their asymmetry, a trait not usually associated with young Type Ia remnants. We present Chandra X-ray Observatory data on these three remnants and perform a maximum likelihood analysis on their spectra. We find that the X-ray emission is dominated by interactions with the interstellar medium. In spite of this, we find a significant Fe overabundance in all three remnants. Through examination of radio, optical, and infrared data, we conclude that these three remnants are likely not Type Ia SNRs. We detect potential point sources that may be members of the progenitor systems of both DEM S 128 and IKT 5, which could suggest these could be Fe-rich core-collapse remnants.

A CHANDRASEKHAR MASS PROGENITOR FOR THE TYPE Ia SUPERNOVA REMNANT 3C 397 FROM THE ENHANCED ABUNDANCES OF NICKEL AND MANGANESE

Despite decades of intense efforts, many fundamental aspects of Type Ia supernova (SNe Ia) remain elusive. One of the major open questions is whether the mass of the exploding white dwarf (WD) is close to the Chandrasekhar limit. Here we report the detection of strong K-shell emission from stable Fe-peak elements in the Suzaku X-ray spectrum of the Type Ia supernova remnant (SNR) 3C 397. The high Ni/Fe and Mn/Fe mass ratios (0.11-0.24 and 0.018-0.033, respectively) in the hot plasma component that dominates the K-shell emission lines indicate a degree of neutronization in the SN ejecta which can only be achieved by electron captures in the dense cores of exploding WDs with a near-Chandrasekhar mass. This suggests a single-degenerate origin for 3C 397, since Chandrasekhar mass progenitors are expected naturally if the WD accretes mass slowly from a companion. Together with other results supporting the double-degenerate scenario, our work adds to the mounting evidence that both progenitor channels make a significant contribution to the SN Ia rate in star-forming galaxies.

ASYMMETRY IN THE OBSERVED METAL-RICH EJECTA OF THE GALACTIC TYPE IA SUPERNOVA REMNANT G299.2–2.9

We have performed a deep Chandra observation of Galactic Type Ia supernova remnant G299.2-2.9. Here we report the initial results from our imaging and spectral analysis. The observed abundance ratios of the central ejecta are in good agreement with those predicted by delayed-detonation Type Ia supernovae models. We reveal inhomogeneous spatial and spectral structures of metal-rich ejecta in G299.2-2.9. The Fe/Si abundance ratio in the northern part of the central ejecta is higher than that in the southern part. An elongation of ejecta material extends out to the western outermost boundary of the remnant. In this western elongation, both the Si and Fe are enriched with a similar abundance ratio to that in the southern part of the central nebula. These structured distributions of metal-rich ejecta material suggest that this Type Ia supernova might have undergone a significantly asymmetric explosion and/or has been expanding into a structured medium.

SEARCH FOR SURVIVING COMPANIONS IN TYPE Ia SUPERNOVA REMNANTS

The nature of the progenitor systems of type~Ia supernovae is still unclear. One way to distinguish between the single-degenerate scenario and double-degenerate scenario for their progenitors is to search for the surviving companions. Using a technique that couples the results from multi-dimensional hydrodynamics simulations with calculations of the structure and evolution of main-sequence- and helium-rich surviving companions, the color and magnitude of main-sequence- and helium-rich surviving companions are predicted as functions of time. The surviving companion candidates in Galactic type~Ia supernova remnants and nearby extragalactic type~Ia supernova remnants are discussed. We find that the maximum detectable distance of main-sequence surviving companions (helium-rich surviving companions) is $0.6-4$~Mpc ($0.4-16$~Mpc), if the apparent magnitude limit is 27 in the absence of extinction, suggesting that the Large and Small Magellanic Clouds and the Andromeda Galaxy are excellent environments in which to search for surviving companions. However, only five Ia~SNRs have been searched for surviving companions, showing little support for the standard channels in the singe-degenerate scenario. To better understand the progenitors of type Ia supernovae, we encourage the search for surviving companions in other nearby type Ia supernova remnants.

REVERSE AND FORWARD SHOCK X-RAY EMISSION IN AN EVOLUTIONARY MODEL OF SUPERNOVA REMNANTS UNDERGOING EFFICIENT DIFFUSIVE SHOCK ACCELERATION

We present new models for the forward and reverse shock thermal X-ray emission from core-collapse and Type Ia supernova remnants (SNRs) which include the efficient production of cosmic rays via non-linear diffusive shock acceleration (DSA). Our CR-hydro-NEI code takes into account non-equilibrium ionization (NEI), hydrodynamic effects of efficient CR production on the SNR evolution, and collisional temperature equilibration among heavy ions and electrons in both the shocked supernova (SN) ejecta and the shocked circumstellar material. While X-ray emission is emphasized here, our code self-consistently determines both thermal and non-thermal broadband emission from radio to TeV energies. We include Doppler broadening of the spectral lines by thermal motions of the ions and by the remnant expansion. We study, in general terms, the roles which the ambient environment, progenitor models, temperature equilibration, and processes related to DSA have on the thermal and non-thermal spectra. The study of X-ray line emission from young SNRs is a powerful tool for determining specific SN elemental contributions, and for providing critical information that helps to understand the type and energetics of the explosion, the composition of the ambient medium in which the SN exploded, and the ionization and dynamics of the hot plasma in the shocked SN ejecta and interstellar medium. With the approaching launch of the next-generation X-ray satellite Astro-H, observations of spectral lines with unprecedented high resolution will become a reality. Our self-consistent calculations of the X-ray spectra from various progenitors will help interpret future observations of SNRs.