SN 1987A at High Resolution

Handed the baton from ROSAT, early observations of SN 1987A with the Chandra HETG and the XMM-Newton RGS showed broad lines with a FWHM of 10^4 km/s: the SN blast wave was continuing to shock the H II region around SN 1987A. Since then, its picturesque equatorial ring (ER) has been shocked, giving rise to a growing, dominant narrow-lined component. Even so, current HETG and RGS observations show that a broad component is still present and contributes 20% of the 0.5--2 keV flux. SN 1987As X-ray behavior can be modeled with a minimum of free parameters as the sum of two simple 1D hydrodynamic simulations: i) an on-going interaction with H II region material producing the broad emission lines and most of the 3--10 keV flux, and ii) an interaction with the dense, clumpy ER material that dominates the 0.5--2 keV flux. Toward the future, we predict a continued growth of the broad component but a drop in the 0.5--2 keV flux, once no new dense ER material is being shocked. When? Time, and new data, will tell.

Evolution and Hydrodynamics of the Very-Broad X-ray Line Emission in SN 1987A

Observations of SN 1987A by the Chandra High Energy Transmission Grating (HETG) in 1999 and the XMM-Newton Reflection Grating Spectrometer (RGS) in 2003 show very broad (v-b) lines with a full-width at half-maximum (FWHM) of order 10^4 kms; at these times the blast wave was primarily interacting with the HII region around the progenitor. Since then, the X-ray emission has been increasingly dominated by narrower components as the blast wave encounters dense equatorial ring (ER) material. Even so, continuing v-b emission is seen in the grating spectra suggesting that interaction with HII region material is on-going. Based on the deep HETG 2007 and 2011 data sets, and confirmed by RGS and other HETG observations, the v-b component has a width of 9300 +/-2000 kms FWHM and contributes of order 20% of the current 0.5--2 keV flux. Guided by this result, SN 1987As X-ray spectra are modeled as the weighted sum of the non-equilibrium-ionization (NEI) emission from two simple 1D hydrodynamic simulations, this 2x1D model reproduces the observed radii, light curves, and spectra with a minimum of free parameters. The interaction with the HII region (rho_init sim 130 amu/cc, +/- 15 degrees opening angle) produces the very-broad emission lines and most of the 3-10 keV flux. Our ER hydrodynamics, admittedly a crude approximation to the multi-D reality, gives ER densities of order 10^4 amu/cc, requires dense clumps (x5.5 density enhancement in sim 30% of the volume), and it predicts that the 0.5-2 keV flux will drop at a rate of sim 17% per year once no new dense ER material is being shocked.

Modeling SN 1996crs X-ray lines at high-resolution: Sleuthing the ejecta/CSM geometry

SN 1996cr, located in the Circinus Galaxy (3.7 Mpc, z ~ 0.001) was non-detected in X-rays at ~ 1000 days yet brightened to ~ 4 x 10^{39} erg/s (0.5-8 keV) after 10 years (Bauer et al. 2008). A 1-D hydrodynamic model of the ejecta-CSM interaction produces good agreement with the measured X-ray light curves and spectra at multiple epochs. We conclude that the progenitor of SN 1996cr could have been a massive star, M > 30 M_solar, which went from an RSG to a brief W-R phase before exploding within its ~ 0.04 pc wind-blown shell (Dwarkadas et al. 2010). Further analysis of the deep Chandra HETG observations allows line-shape fitting of a handful of bright Si and Fe lines in the spectrum. The line shapes are well fit by axisymmetric emission models with an axis orientation ~ 55 degrees to our line-of-sight. In the deep 2009 epoch the higher ionization Fe XXVI emission is constrained to high lattitudes: the Occam-est way to get the Fe H-like emission coming from high latitude/polar regions is to have more CSM at/around the poles than at mid and lower lattitudes, along with a symmetric ejecta explosion/distribution. Similar CSM/ejecta characterization may be possible for other SNe and, with higher-throughput X-ray observations, for gamma-ray burst remnants as well.

Kinematics of Supernova Remnants: Status of X-Ray Observations

A supernova (SN) explosion drives stellar debris into the circumstellar material (CSM) filling a region on a scale of parsecs with X-ray emitting plasma. The velocities involved in supernova remnants (SNRs), thousands of km/s, can be directly measured with medium and high-resolution X-ray spectrometers and add an important dimension to our understanding of the last stages of the progenitor, the explosion mechanism, and the physics of strong shocks. After touching on the ingredients of SNR kinematics, I present a summary of the still-growing measurement results from SNR X-ray observations. Given the advances in 2D/3D hydrodynamics, data analysis techniques, and especially X-ray instrumentation, it is clear that our view of SNRs will continue to deepen in the decades ahead.

Simulation and Fitting of Multi-Dimensional X-ray Data

Astronomical data generally consists of 2 or more high-resolution axes, e.g., X,Y position on the sky or wavelength and position-along-one-axis (long-slit spectrometer). Analyzing these multi-dimension observations requires combining 3D source models (including velocity effects), instrument models, and multi-dimensional data comparison and fitting. A prototype of such a Beyond-XSPEC (Noble & Nowak, 2008) system is presented here using Chandra imag- ing and dispersed HETG grating data. Techniques used include: Monte Carlo event generation, chi-squared comparison, conjugate gradient fitting adapted to the Monte Carlo characteristics, and informative visualizations at each step. These simple baby steps of progress only scratch the surface of the computational potential that is available these days for astronomical analysis.

Chandra HETG Spectra of SN 1987A at 20 years

We have undertaken deep, high-resolution observations of SN 1987A at ~20 years after its explosion with the Chandra HETG and LETG spectrometers. Here we present the HETG X-ray spectra of SN 1987A having unprecedented spectral resolution and signal-to-noise in the 6 A to 20 A bandpass, which includes the H-like and He-like lines of Si, Mg, Ne, as well as O VIII lines and bright Fe XVII lines. In joint analysis with LETG data, we find that there has been a significant decrease from 2004 to 2007 in the average temperature of the highest temperature component of the shocked-plasma emission. Model fitting of the profiles of individual HETG lines yields bulk kinematic velocities of the higher-Z ions, Mg and Si, that are significantly lower than those inferred from the LETG 2004 observations.