Search

Updated imaging and photometric results from Chandra observations of SN 1987A, covering the last 16 years, are presented. We find that the 0.5-2 keV light curve has remained constant at ~8x10^-12 erg s^-1 cm^-2 since 9500 days, with the 3-8 keV light curve continuing to increase until at least 10000 days. The expansion rate of the ring is found to be energy dependent, such that after day 6000 the ring expands faster in the 2-10 keV band than it does at energies <2 keV. Images show a reversal of the east-west asymmetry between 7000 and 8000 days after the explosion. The latest images suggest the southeastern side of the equatorial ring is beginning to fade. Consistent with the latest optical and infrared results, our Chandra analysis indicates the blast wave is now leaving the dense equatorial ring, which marks the beginning of a major change in the evolutionary phase of the supernova remnant 1987A.

The Galactic blue supergiant SBW1 with its circumstellar ring nebula represents the best known analog of the progenitor of SN 1987A. High-resolution imaging has shown H-alpha and IR structures arising in an ionized flow that partly fills the ring's interior. To constrain the influence of the stellar wind on this structure, we obtained an ultraviolet (UV) spectrum of the central star of SBW1 with the HST Cosmic Origins Spectrograph (COS). The UV spectrum shows none of the typical wind signatures, indicating a very low mass-loss rate. Radiative transfer models suggest an extremely low rate below 10$^{-10}$ Msun/yr, although we find that cooling timescales probably become comparable to or longer than the flow time below 10$^{-8}$ Msun/yr. We therefore adopt this latter value as a conservative upper limit. For the central star, the model yields $T_{\rm eff}$=21,000$\pm$1000 K, $L\simeq$5$\times$10$^4$ $L_{\odot}$, and roughly Solar composition except for enhanced N abundance. SBW1's very low mass-loss rate may hinder the wind's ability to shape the surrounding nebula. The very low mass-loss rate also impairs the wind's ability to shed angular momentum; the spin-down timescale for magnetic breaking is more than 500 times longer than the age of the ring. This, combined with the star's slow rotation rate, constrain merger scenarios to form ring nebulae. The mass-loss rate is at least 10 times lower than expected from mass-loss recipes, without any account of clumping. The physical explanation for why SBW1's wind is so weak presents an interesting mystery., Comment: 12 pages, 6 figs. submitted to MNRAS. comments welome

Both CO and SiO have been observed at early and late phases in SN 1987A. H_2 was predicted to form at roughly the same time as these molecules, but was not detected at early epochs. Here we report the detection of NIR lines from H_2 at 2.12 mu and 2.40 mu in VLT/SINFONI spectra obtained between days 6489 and 10,120. The emission is concentrated to the core of the supernova in contrast to H-alpha and approximately coincides with the [Si I]/[Fe II] emission detected previously in the ejecta. Different excitation mechanisms and power sources of the emission are discussed. From the nearly constant H_2 luminosities we favour excitation resulting from the 44Ti decay., Comment: 5 pages 3 figures, accepted for ApJ Letters

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots up to day ~8,000 (~2009), both this and the unshocked emission are now fading. From the radial positions of the hotspots we see an acceleration of these up to 500-1000 km/s, consistent with the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations (2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from the shocks or by direct shock interaction. All of these observations indicate that the interaction with the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed decay, that the inner ring will be destroyed by ~2025., Comment: 5 pages, 4 figures, accepted for ApJ Letters

We present new {\it Hubble Space Telescope} images of high-velocity H-$\alpha$ and Lyman-$\alpha$ emission in the outer debris of SN~1987A. The H-$\alpha$ images are dominated by emission from hydrogen atoms crossing the reverse shock. For the first time we observe emission from the reverse shock surface well above and below the equatorial ring, suggesting a bipolar or conical structure perpendicular to the ring plane. Using the H$\alpha$ imaging, we measure the mass flux of hydrogen atoms crossing the reverse shock front, in the velocity intervals ($-$7,500~$<$~$V_{obs}$~$<$~$-$2,800 km s$^{-1}$) and (1,000~$<$~$V_{obs}$~$<$~7,500 km s$^{-1}$), $\dot{M_{H}}$ = 1.2~$\times$~10$^{-3}$ M$_{\odot}$ yr$^{-1}$. We also present the first Lyman-$\alpha$ imaging of the whole remnant and new $Chandra$ X-ray observations. Comparing the spatial distribution of the Lyman-$\alpha$ and X-ray emission, we observe that the majority of the high-velocity Lyman-$\alpha$ emission originates interior to the equatorial ring. The observed Lyman-$\alpha$/H-$\alpha$ photon ratio, $\langle$$R(L\alpha / H\alpha)$$\rangle$ $\approx$~17, is significantly higher than the theoretically predicted ratio of $\approx$ 5 for neutral atoms crossing the reverse shock front. We attribute this excess to Lyman-$\alpha$ emission produced by X-ray heating of the outer debris. The spatial orientation of the Lyman-$\alpha$ and X-ray emission suggests that X-ray heating of the outer debris is the dominant Lyman-$\alpha$ production mechanism in SN 1987A at this phase in its evolution., Comment: 6 pages, 5 figures. ApJL - accepted

We present observations with VLT and HST of the broad emission lines from the inner ejecta and reverse shock of SN 1987A from 1999 until 2012 (days 4381 -- 9100 after explosion). We detect broad lines from H-alpha, H-beta, Mg I], Na I, [O I], [Ca II] and a feature at 9220 A. We identify the latter line with Mg II 9218, 9244,most likely pumped by Ly-alpha fluorescence. H-alpha, and H-beta both have a centrally peaked component, extending to 4500 km/s and a very broad component extending to 11,000 km/s, while the other lines have only the central component. The low velocity component comes from unshocked ejecta, heated mainly by X-rays from the circumstellar ring collision, whereas the broad component comes from faster ejecta passing through the reverse shock. The reverse shock flux in H-alpha has increased by a factor of 4-6 from 2000 to 2007. After that there is a tendency of flattening of the light curve, similar to what may be seen in soft X-rays and in the optical lines from the shocked ring. The core component seen in H-alpha, [Ca II] and Mg II has experienced a similar increase, consistent with that found from HST photometry. The ring-like morphology of the ejecta is explained as a result of the X-ray illumination, depositing energy outside of the core of the ejecta. The energy deposition in the ejecta of the external X-rays illumination is calculated using explosion models for SN 1987A and we predict that the outer parts of the unshocked ejecta will continue to brighten because of this. We finally discuss evidence for dust in the ejecta from line asymmetries., Comment: Submitted to ApJ. Comments welcome

We present a study of the morphology of the ejecta in Supernova 1987A based on images and spectra from the HST as well as integral field spectroscopy from VLT/SINFONI. The HST observations were obtained between 1994 - 2011 and primarily probe the outer hydrogen-rich zones of the ejecta. The SINFONI observations were obtained in 2005 and 2011 and instead probe the [Si I]/[Fe II] emission from the inner regions. We find a strong temporal evolution of the morphology in the HST images, from a roughly elliptical shape before ~5,000 days, to a more irregular, edge-brightened morphology thereafter. This transition is a natural consequence of the change in the dominant energy source powering the ejecta, from radioactive decay before ~5,000 days to X-ray input from the circumstellar interaction thereafter. The [Si I]/[Fe II] images display a more uniform morphology, which may be due to a remaining significant contribution from radioactivity in the inner ejecta and the higher abundance of these elements in the core. Both the H-alpha and the [Si I]/[Fe II] line profiles show that the ejecta are distributed fairly close to the plane of the inner circumstellar ring, which is assumed to define the rotational axis of the progenitor. The H-alpha emission extends to higher velocities than [Si I]/[Fe II] as expected. There is no clear symmetry axis for all the emission and we are unable to model the ejecta distribution with a simple ellipsoid model with a uniform distribution of dust. Instead, we find that the emission is concentrated to clumps and that the emission is distributed somewhat closer to the ring in the north than in the south. This north-south asymmetry may be partially explained by dust absorption. We compare our results with explosion models and find some qualitative agreement, but note that the observations show a higher degree of large-scale asymmetry., Comment: Accepted for publication in ApJ

We present the most sensitive ultraviolet observations of Supernova 1987A to date. Imaging spectroscopy from the Hubble Space Telescope-Cosmic Origins Spectrograph shows many narrow (dv \sim 300 km/s) emission lines from the circumstellar ring, broad (dv \sim 10 -- 20 x 10^3 km/s) emission lines from the reverse shock, and ultraviolet continuum emission. The high signal-to-noise (> 40 per resolution element) broad LyA emission is excited by soft X-ray and EUV heating of mostly neutral gas in the circumstellar ring and outer supernova debris. The ultraviolet continuum at \lambda > 1350A can be explained by HI 2-photon emission from the same region. We confirm our earlier, tentative detection of NV \lambda 1240 emission from the reverse shock and we present the first detections of broad HeII \lambda1640, CIV \lambda1550, and NIV] \lambda1486 emission lines from the reverse shock. The helium abundance in the high-velocity material is He/H = 0.14 +/- 0.06. The NV/H-alpha line ratio requires partial ion-electron equilibration (T_{e}/T_{p} \approx 0.14 - 0.35). We find that the N/C abundance ratio in the gas crossing the reverse shock is significantly higher than that in the circumstellar ring, a result that may be attributed to chemical stratification in the outer envelope of the supernova progenitor. The N/C abundance ratio may have been stratified prior to the ring expulsion, or this result may indicate continued CNO processing in the progenitor subsequent to the expulsion of the circumstellar ring., Comment: 12 pages, 8 figures. ApJ - accepted

We have been monitoring the supernova remnant (SNR) 1987A with {\it Chandra} observations since 1999. Here we report on the latest change in the soft X-ray light curve of SNR 1987A. For the last $\sim$1.5 yr (since day $\sim$8000), the soft X-ray flux has significantly flattened, staying (within uncertainties) at $f_{\rm X}$ $\sim$ 5.7 $\times$ 10$^{-12}$ erg cm$^{-2}$ s$^{-1}$ (corresponding to $L_{\rm X}$ $\sim$ 3.6 $\times$ 10$^{36}$ erg s$^{-1}$) in the 0.5--2 keV band. This remarkable change in the recent soft X-ray light curve suggests that the forward shock is now interacting with a decreasing density structure, after interacting with an increasing density gradient over $\sim$10 yr prior to day $\sim$8000. Possibilities may include the case that the shock is now propagating beyond a density peak of the inner ring. We briefly discuss some possible implications on the nature of the progenitor and the future prospects of our {\it Chandra} monitoring observations., Comment: 11 pages, 2 tables, 1 color figure, Accepted by ApJ Letters

Observations with the Hubble Space Telescope (HST), conducted since 1990, now offer an unprecedented glimpse into fast astrophysical shocks in the young remnant of supernova 1987A. Comparing observations taken in 2010 using the refurbished instruments on HST with data taken in 2004, just before the Space Telescope Imaging Spectrograph failed, we find that the Ly-a and H-a lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We observe broad blueshifted Ly-a, which we attribute to resonant scattering of photons emitted from hotspots on the equatorial ring. We also detect NV~\lambda\lambda 1239,1243 A line emission, but only to the red of Ly-A. The profiles of the NV lines differ markedly from that of H-a, suggesting that the N^{4+} ions are scattered and accelerated by turbulent electromagnetic fields that isotropize the ions in the collisionless shock., Comment: Science, accepted. Science Express, 02 Sept 2010. 5 figures. Supporting online material can be found at http://www.sciencemag.org/cgi/content/full/sci;science.1192134/DC1

We continue to explore the validity of the reflected shock structure (RSS) picture in SNR 1987A that was proposed in our previous analyses of the X-ray emission from this object. We used an improved version of our RSS model in a global analysis of 14 CCD spectra from the monitoring program with Chandra. In the framework of the RSS picture, we are able to match both the expansion velocity curve deduced from the analysis of the X-ray images and light curve. Using a simplified analysis, we also show that the X-rays and the non-thermal radio emission may originate from the same shock structure (the blast wave). We believe that using the RSS model in the analysis of grating data from the Chandra monitoring program of SNR 1987A that cover a long enough time interval, will allow us to build a more realistic physical picture and model of SNR 1987A., Comment: 14 pages, 1 Table, 8 figures, accepted for publication in MNRAS

We present the evolution of the radial expansion of SNR 1987A as measured using Chandra X-ray observations taken over the last 10 years. To characterize the complex structure of the remnant and isolate the expansion measurement, we fit the images to several empirical models including: a simple circular torus, a torus with bilateral lobes, and a torus with four tangentially extended lobes. We discuss the results of this measure in the context of the overall evolution of the supernova remnant, for which we believe we have measured the end of the free expansion phase and its transition to the adiabatic phase (at least along the equatorial ring). The timing of this event is in agreement with early predictions of the remnant evolution., Comment: Accepted for publication in ApJ, 21 pages, 4 figures, 3 tables

We present an extended analysis of the deep Chandra LETG and HETG observations of the supernova remnant 1987A (SNR 1987A) carried out in 2007. The global fits to the grating spectra show that the temperature of the X-ray emitting plasma in the slower shocks in this system has remained stable for the last three years, while that in the faster shocks has decreased. This temperature evolution is confirmed by the first light curves of strong X-ray emission lines and their ratios. On the other hand, bulk gas velocities inferred from the X-ray line profiles are too low to account for the post-shock plasma temperatures inferred from spectral fits. This suggests that the X-ray emission comes from gas that has been shocked twice, first by the blast wave and again by shocks reflected from the inner ring of SNR 1987A. A new model that takes these considerations into account gives support to this physical picture., Comment: 36 pages, 10 figures, Accepted for publication in ApJ

Balmer-dominated shocks in supernova remnants (SNRs) produce strong hydrogen lines with a two-component profile composed of a narrow contribution from cold upstream hydrogen atoms, and a broad contribution from hydrogen atoms that have undergone charge transfer reactions with hot protons. Observations of emission lines from edge-wise shocks in SNRs can constrain the gas velocity and collisionless electron heating at the shock front. Downstream hydrogen atoms engage in charge transfer, excitation and ionization reactions, defining an interaction region called the shock transition zone. The properties of hot hydrogen atoms produced by charge transfers (called broad neutrals) are critical for accurately calculating the structure and radiation from the shock transition zone. This paper is the third in a series describing the kinetic, fluid and emission properties of Balmer-dominated shocks, and is the first to properly treat the effect of broad neutral kinetics on shock transition zone structure. We use our models to extract shock parameters from observations of Balmer-dominated SNRs. We find that inferred shock velocities and electron temperatures are lower than those of previous calculations by <10% for v_s<1500 km/s, and by 10-30% for v_s>1500 km/s. This effect is primarily due to the fact that excitation by proton collisions and charge transfer to excited levels favor the high speed part of the neutral hydrogen velocity distribution. Our results have a strong dependence on the ratio of electron to proton temperatures, \beta=T_e/T_p, which allows us to construct a relation \beta(v_s) between the temperature ratio and shock velocity. We compare our calculations to previous results by Ghavamian et al. (2007)., Comment: 41 pages, 15 figures, 2 tables. Improved comparison to previous results, added discussion, and incorporated referee's suggestions. Submitted to ApJ

We explore a simple model for the high luminosity of SN 2006gy involving photon diffusion of shock-deposited thermal energy. The distinguishing property of the model is that the large ``stellar'' radius of 160 AU required to prevent adiabatic losses is not the true stellar radius, but rather, the radius of an opaque, unbound circumstellar envelope, created when 10 Msun was ejected in the decade before the supernova in an eruption analogous to that of eta Carinae. The supernova light is produced primarily by diffusion of thermal energy following the passage of the blast wave through this shell. This model differs from traditional models of supernova debris interacting with external CSM in that here the shell is optically thick and the escape of radiation is delayed. We show that any model attempting to account for SN2006gy's huge luminosity with radiation emitted by ongoing CSM interaction fails for the following basic reason: the CSM density required to achieve the observed luminosity makes the same circumstellar envelope opaque, forcing a thermal diffusion solution. In our model, the weaker CSM interaction giving rise to SN2006gy's characteristic Type IIn spectrum and soft X-rays is not linked to the power source of the visual continuum; instead, it arises after the blast wave breaks free of the opaque shell into the surrounding wind. While a simple diffusion model can explain the gross properties of the early light curve of SN2006gy, it predicts that the light curve must plummet rapidly at late-times, unless an additional power source is present., Comment: ApJ Letters, accepted

We present the first high-dispersion spectroscopy of He I 10830 from the hotspots in the ring around SN1987A, obtained at Gemini South, spatially resolving the near and far sides of the ring. We compare these line profiles to similar echelle spectra of H$\alpha$ and [N II] 6583 obtained at the Magellan Observatory. We find that the He I profiles are much broader than H-alpha or [N II], but the He I profiles also have different shapes -- they have enhanced emission at high speeds, with extra blueshifted emission on the north side of the ring, and extra redshifted emission on the south side. To explain this, we invoke a simple geometric picture where the extra He I emission traces hotter gas from faster shocks that strike the apex of the hotspots directly, while the H-alpha preferentially traces cooler lower-ionization gas from slower transverse shocks that penetrate into the sides of the ring., Comment: 3 pages. To appear in proceedings: "Supernova 1987A: 20 Years After: Supernovae and Gamma-Ray Bursters" AIP, New York, eds. S. Immler, K.W. Weiler, and R. McCray

We examine the structure of the post-shock region in supernova remnants (SNRs). The ``shock transition zone'' is set up by charge transfer and ionization events between atoms and ions, and has a width $\sim 10^{15}$ cm$^{-2}$ $n^{-1}_0$, where $n_0$ is the total pre-shock density (including both atoms and ions). For Balmer-dominated SNRs with shock velocity $v_s \gtrsim 1000$ km s$^{-1}$, the Rankine-Hugoniot conditions for ion velocity and temperature are obeyed instantly, leaving the full width at half-maximum (FWHM) of the broad H$\alpha$ line versus $v_s$ relation intact. However, the spatial variation in the post-shock densities is relevant to the problem of Ly$\alpha$ resonant scattering in young, core-collapse SNRs. Both two- (pre-shock atoms and ions) and three-component (pre-shock atoms, broad neutrals and ions) models are considered. We compute the spatial emissivities of the broad ($\xi_b$) and narrow ($\xi_n$) H$\alpha$ lines; a calculation of these emissivities in SN 1006 is in general agreement with the computed ones of Raymond et al. (2007). The (dimensionless) spatial shift, $\Theta_{\rm{shift}}$, between the centroids of $\xi_b$ and $\xi_n$ is unique for a given shock velocity and $f_{\rm{ion}}$, the pre-shock ion fraction. Measurements of $\Theta_{\rm{shift}}$ can be used to constrain $n_0$., Comment: 25 pages, 8 figures. Accepted by Astrophysical Journal

We have been monitoring Supernova (SN) 1987A with {\it Chandra X-Ray Observatory} since 1999. We present a review of previous results from our {\it Chandra} observations, and some preliminary results from new {\it Chandra} data obtained in 2006 and 2007. High resolution imaging and spectroscopic studies of SN 1987A with {\it Chandra} reveal that X-ray emission of SN 1987A originates from the hot gas heated by interaction of the blast wave with the ring-like dense circumstellar medium (CSM) that was produced by the massive progenitor's equatorial stellar winds before the SN explosion. The blast wave is now sweeping through dense CSM all around the inner ring, and thus SN 1987A is rapidly brightening in soft X-rays. At the age of 20 yr (as of 2007 January), X-ray luminosity of SN 1987A is $L_{\rm X}$ $\sim$ 2.4 $\times$ 10$^{36}$ ergs s$^{-1}$ in the 0.5$-$10 keV band. X-ray emission is described by two-component plane shock model with electron temperatures of $kT$ $\sim$ 0.3 and 2 keV. As the shock front interacts with dense CSM all around the inner ring, the X-ray remnant is now expanding at a much slower rate of $v$ $\sim$ 1400 km s$^{-1}$ than it was until 2004 ($v$ $\sim$ 6000 km s$^{-1}$)., Comment: 8 pages, 4 figures (including 2 color figs), 3 tables. To appear in the proceedings of "Supernova 1987A: 20 Years after Supernovae and Gamma-Ray Bursters", held in Aspen CO. USA, Feb 19-23, 2007

We present a new formalism to describe the ratios and profiles of emission lines from hydrogen in Balmer-dominated shocks. We use this model to interpret the measured widths and ratios of broad and narrow H-alpha, H-beta and Ly-alpha emission lines in supernova remnants (SNRs). Our model results agree fairly well with those obtained previously by Chevalier, Kirshner & Raymond (1980) and are consistent with observations of several SNRs. The same model fails to account for the ratio of broad to narrow line emission from the reverse shock in SNR 1987A as observed by Heng et al. (2006). We suggest that this discrepancy between theory and observation results from a faulty assumption that Balmer-dominated shocks can be treated as sharp discontinuities. If the spatial structure of the shock transition zone is taken into account, the predicted ratios of broad to narrow line emission in most SNRs will change by modest factors, but the ratio in SNR 1987A will increase substantially. Significantly greater shock velocities will be required to account for the observed full widths at half-maximum of the broad emission lines in most SNRs., Comment: 34 pages, 12 figures. Accepted by ApJ

$\sim$18 yr after the supernova explosion, the blast wave of SNR 1987A is entering the main body of the equatorial circumstellar material, which is causing a dramatic brightening of the remnant. We recently reported the observational evidence for this event from our {\it Chandra} data (Park et al. 2005b; P05 hereafter). We present here the temporal evolution of the X-ray emitting shock parameters and the detailed description of the spectral and image analysis of SNR 1987A, on which P05 was based. While the remnant becomes brighter, the softening of the overall X-ray spectrum continues and is enhanced on around day $\sim$6200 (since the explosion). The two-component shock model indicates that the electron temperatures have been changing for the last $\sim$6 yr. The X-ray spectrum is now described by $kT$ $\sim$ 0.3 keV and 2.3 keV thermal plasmas which are believed to characteristically represent the shock-heated density gradient along the boundary between the H{\small II} region and the dense inner ring. As the blast wave sweeps through the inner circumstellar ring shining in X-rays, we expect that the shock parameters continue to change, revealing the density and abundance structure of the inner ring. Follow-up {\it Chandra} observations will thus uncover the past history of the progenitor's stellar evolution. The origin of the relatively faint hard X-ray emission ($E$ $>$ 3 keV) from SNR 1987A is still unclear (thermal vs. nonthermal). Considering the continuous brightening of the hard band intensity, as well as the soft band flux, follow-up monitoring observations will also be essential to reveal the origin of the hard X-ray emission of SNR 1987A., Comment: 21 pages including 9 figures, Accepted by ApJ

We discuss the results from deep Chandra LETG observations of the supernova remnant 1987A (SNR 1987A). We find that a distribution of shocks, spanning the same range of velocities (from 300 to 1700 km/s) as deduced in the first part of our analysis (Zhekov et al. 2005, ApJL, 628, L127), can account for the entire X-ray spectrum of this object. The post-shock temperature distribution is bimodal, peaking at kT 0.5 and 3 keV. Abundances inferred from the X-ray spectrum have values similar to those for the inner circumstellar ring, except that the abundances of nitrogen and oxygen are approximately a factor of two lower than those inferred from the optical/UV spectrum. The velocity of the X-ray emitting plasma has decreased since 1999, apparently because the blast wave has entered the main body of the inner circumstellar ring., Comment: Accepted for publication in the ApJ

We present new (2004 July) G750L and G140L Space Telescope Imaging Spectrograph (STIS) data of the H-alpha and Ly-alpha emission from supernova remnant (SNR) 1987A. With the aid of earlier data, from Oct 1997 to Oct 2002, we track the local evolution of Ly-alpha emission and both the local and global evolution of H-alpha emission. In addition to emission which we can clearly attribute to the surface of the reverse shock, we also measure comparable emission, in both H-alpha and Ly-alpha, which appears to emerge from supernova debris interior to the surface. New observations taken through slits positioned slightly eastward and westward of a central slit show a departure from cylindrical symmetry in the H-alpha surface emission. Using a combination of old and new observations, we construct a light curve of the total H-alpha flux, F, from the reverse shock, which has increased by a factor ~ 4 over about 8 years. However, due to large systematic uncertainties, we are unable to discern between the two limiting behaviours of the flux - F ~ t (self-similar expansion) and F ~ t^5 (halting of the reverse shock). Such a determination is relevant to the question of whether the reverse shock emission will vanish in less than about 7 years (Smith et al. 2005). Future deep, low- or moderate-resolution spectra are essential for accomplishing this task., Comment: 28 pages, 12 figures. Accepted by ApJ

We report on the morphological and spectral evolution of SNR 1987A from the monitoring observations with the Chandra/ACIS. As of 2005, the X-ray-bright lobes are continuously brightening and expanding all around the ring. The softening of the overall X-ray spectrum also continues. The X-ray lightcurve is particularly remarkable: i.e., the recent soft X-ray flux increase rate is significantly deviating from the model which successfully fits the earlier data, indicating even faster flux increase rate since early 2004 (day ~6200). We also report results from high resolution spectral analysis with deep Chandra/LETG observations. The high resolution X-ray line emission features unambiguously reveal that the X-ray emission of SNR 1987A is originating primarily from a "disk" along the inner ring rather than from a spherical shell. We present the ionization structures, elemental abundances, and the shock velocities of the X-ray emitting plasma., Comment: 6 pages (A4 size) including 8 figures & 1 table. To appear in the proceedings of "The X-Ray Universe 2005" conference, Sept 26-30, 2005, El Escorial, Madrid, Spain

We use low-dispersion spectra obtained at the Magellan Observatory to study the broad H-alpha emission from the reverse shock of the infant supernova remnant SNR1987A. These spectra demonstrate that the spatio-kinematic structure of the reverse shock can be distinguished from that of the circumstellar ring and hotspots, even at ground-based spatial resolution. We measure a total dereddened H-alpha flux of 1.99(pm0.22)e-13 ergs/s/cm2 at an epoch 18.00 years after outburst. At 50 kpc, the total reverse shock luminosity in H-alpha is roughly 15 Lsun, which implies a total flux of neutral hydrogen atoms across the reverse shock of 8.9e46/s, or roughly 2.3e-3 Msun/yr. This represents an increase by a factor ~4 since 1997. Lyman continuum radiation from gas shocked by the forward blast wave can ionize neutral hydrogen atoms in the supernova debris before they reach the reverse shock. If the inward flux of ionizing photons exceeds the flux of hydrogen atoms approaching the reverse shock, this pre-ionization will shut off the broad Ly-alpha and H-alpha emission. The observed X-ray emission of SNR1987A implies that the ratio of ionizing flux to hydrogen atom flux across the reverse shock is presently at least 0.04. The X-ray emission is increasing much faster than the flux of atoms, and if these trends continue, we estimate that the broad Ly-alpha and H-alpha emission will vanish within about 7 years., Comment: 5 pgs, Fig 1 in color, accepted by ApJ Letters

We report an up-turn in the soft X-ray light curve of supernova remnant (SNR) 1987A in late 2003 (~6200 days after the explosion), as observed with the Chandra X-ray Observatory. Since early 2004, the rapid increase of the 0.5-2 keV band X-ray light curve can no longer be described by the exponential density distribution model with which we successfully fitted the data between 1990 and 2003. Around day ~6200, we also find that the fractional contribution to the observed soft X-ray flux from the decelerated shock begins to exceed that of the fast shock and that the X-ray brightening becomes "global" rather than "spotty". We interpret these results as evidence that the blast wave has reached the main body of the dense circumstellar material all around the inner ring. This interpretation is supported by other recent observations, including a deceleration of the radial expansion of the X-ray remnant, a significant up-turn in the mid-IR intensities, and the prevalence of the optical hot spots around the entire inner ring, all of which occur at around day 6000. In contrast to the soft X-ray light curve, the hard band (3-10 keV) X-ray light curve increases at a much lower rate which is rather similar to the radio light curve. The hard X-ray emission may thus originate from the reverse shock where the radio emission is likely produced. Alternatively, the low increase rate of the hard X-rays may simply be a result of the continuous softening of the overall X-ray spectrum., Comment: AASTex preprint style 12 pages including 1 table and 4 figures, Accepted by ApJL

We present an update on the results of our monitoring observations of the X-ray remnant of supernova (SN) 1987A with the {\it Chandra X-Ray Observatory}. As of 2002 December, we have performed a total of seven observations of SN 1987A. The high angular resolution images from the latest data reveal developments of new X-ray bright spots in the northwestern and the southwestern portions of the remnant as well as changes on the eastern side. The latest 0.5-2 keV band flux ($f_X$ $\sim$ 6 $\times$ 10$^{-13}$ ergs cm$^{-2}$ s$^{-1}$) is four times brighter than three years earlier. The overall X-ray emission is primarily from the blast wave shock with $kT$ $\sim$ 2.4 keV. As the blast wave approaches the dense circumstellar material, the contribution from the decelerated slow shock ($kT$ $\sim$ 0.22 keV) to the observed X-ray emission is becoming significant. The increase of this slow shock contribution over the last two years is particularly noticeable in the western half of the remnant. These results indicate that the shock front is now reaching the main body of the inner circumstellar ring. Based on the best-fit two-shock spectral model, we derive approximate densities of the X-ray-emitting regions ($n_e$ $\sim$ 235 cm$^{-3}$ for the fast shock and $n_e$ $\sim$ 7500 cm$^{-3}$ for the slow shock). We obtain an upper limit on the observed X-ray luminosity of any embedded point source ($L_X$ $\le$ 1.5 $\times$ 10$^{34}$ ergs s$^{-1}$) in the 2$-$10 keV band. The X-ray remnant continues to expand linearly at a rate of 4167 km s$^{-1}$., Comment: 22 pages (ApJ preprint style), 7 Figures, Accepted by ApJ (scheduled on July 20, 2004), for high-quality Fig 1 and Fig 2, please contact park@astro.psu.edu

We report and interpret HST/STIS long-slit observations of the optical and ultraviolet (1150 - 10270 Angstrom) emission-line spectra of the rapidly brightening Spot 1 on the equatorial ring of SN 1987A between 1997 September and 1999 October (days 3869 -- 4606 after outburst). The emission is caused by radiative shocks created where the supernova blast wave strikes dense gas protruding inward from the equatorial ring. We measure and tabulate line identifications, fluxes and, in some cases, line widths and shifts. We compute flux correction factors to account for substantial interstellar line absorption of several emission lines. Nebular analysis shows that optical emission lines come from a region of cool (T_e ~ 10^4 K) and dense (n_e ~ 10^6 cm^-3) gas in the compressed photoionized layer behind the radiative shock. The observed line widths indicate that only shocks with shock velocities V_s < 250 km/s have become radiative, while line ratios indicate that much of the emission must have come from yet slower (V_s < 135 k/ms) shocks. We are able to fit the UV fluxes with an idealized radiative shock model consisting of two shocks (V_s = 135 and 250 km/s). The observed UV flux increase with time can be explained by the increase in shock surface areas as the blast wave overtakes more of the protrusion. The observed flux ratios of optical to highly-ionized UV lines are greater by a factor of ~ 2 -- 3 than predictions from the radiative shock models and we discuss the possible causes. We also present models for the observed H-alpha line widths and profiles, which suggests that a chaotic flow exists in the photoionized regions of these shocks. We discuss what can be learned with future observations of all the spots present on the equatorial ring., Comment: 67 pages, 21 figures. Accepted for publication in the Astrophysical Journal

We discuss the X-ray emission observed from Supernova Remnant 1987A with the Chandra X-ray Observatory. We analyze a high resolution spectrum obtained in 1999 October with the high energy transmission grating (HETG). From this spectrum we measure the strengths and an average profile of the observed X-ray lines. We also analyze a high signal-to-noise ratio CCD spectrum obtained in 2000 December. The good statistics (~ 9250 counts) of this spectrum and the high spatial resolution provided by the telescope allow us to perform spectroscopic analyses of different regions of the remnant. We discuss the relevant shock physics that can explain the observed X-ray emission. The X-ray spectra are well fit by plane parallel shock models with post-shock electron temperatures of ~ 2.6 keV and ionization ages of ~ 6 x 10^10 cm^3/s. The combined X-ray line profile has a FWHM of ~ 5000 km/s, indicating a blast wave speed of ~ 3500 km/s. At this speed, plasma with a mean post-shock temperature of ~ 17 keV is produced. This is direct evidence for incomplete electron-ion temperature equilibration behind the shock. Assuming this shock temperature, we constrain the amount of collisionless electron heating at the shock front at T_e0 / T_s = 0.11 (+0.02,-0.01). We find that the plasma has low metallicity (abundances are ~ 0.1 -- 0.4 solar) and is nitrogen enriched (N/O ~ 0.8 by number), similar to that found for the equatorial ring. Analysis of the spectra from different regions of the remnant reveals slight differences in the parameters of the emitting plasma. The plasma is cooler near the optical Spot 1 (at position angle ~ 30 degrees) and in the eastern half of the remnant, where the bright optical spots are found, than in the western half, consistent with the presence of slower (~ 500 km/s) shocks entering denser ring material., Comment: 36 pages, 9 figures. Submitted to ApJ on Dec. 11, 2001; accepted for publication on Mar. 8, 2002. New version has text adjusted to referee's Comments

We report on the results of our monitoring program of the remnant of SN 1987A with the Advanced CCD Imaging Spectrometer (ACIS) on board the {\it Chandra X-ray Observatory}. Two new observations have been performed in AO2, bringing the total to four monitoring observations over the past two years. Over this time period, new techniques for correction of ``Charge Transfer Inefficiency (CTI)'' and for use of charge spreading to provide angular resolution somewhat better than the pixel size of the CCD detector have become available at Penn State. We have processed all four observations using sub-pixel resolution to obtain the highest possible angular resolution, and using our CTI correction software to provide more reliable spectral analysis and flux estimations. The high angular resolution images indicate that the X-ray bright knots are convincingly correlated with the optical spots, primarily at $\la$1 keV, while higher energy photons are very well correlated with radio images. Our data also provide marginal evidence for radial expansion of the X-ray remnant at a rate of 5200 $\pm$ 2100 km s$^{-1}$. The X-ray flux appears to linearly increase by $\sim$60% over the 18 month period of these observations. The spectrum is dominated by broad complexes of atomic emission lines and can be fit with a simple model of a plane-parallel shock with electron temperatures of $kT$ $\sim$ 2 $-$ 4 keV and a postshock electron density of $n_{e}$ $\sim$ 210 $-$ 420 cm$^{-3}$. The implied 0.5 $-$ 10 keV band luminosity in 2001 April is $\sim$1.3 $\times$ 10$^{35}$ ergs s$^{-1}$; as of that date, we still observe no direct evidence for the central point source, with an upper limit on the {\it observed} luminosity of $L_{X}$ $\sim$ 5.5 $\times$ 10$^{33}$ ergs s$^{-1}$ in the 2 $-$ 10 keV band., Comment: 29 pages, 5 figures (2 color images), Accepted for ApJ

We present high resolution Chandra observations of the remnant of SN1987A in the Large Magellanic Cloud. The high angular resolution of the Chandra X-ray Observatory (CXO) permits us to resolve the X-ray remnant. We find that the remnant is shell-like in morphology, with X-ray peaks associated with some of the optical hot spots seen in HST images. The X-ray light curve has departed from the linear flux increase observed by ROSAT, with a 0.5-2.0 keV luminosity of 1.5 x 10^35 erg/s in January 2000. We set an upper limit of 2.3 x 10^34 ergs/s on the luminosity of any embedded central source (0.5 - 2 keV). We also present a high resolution spectrum, showing that the X-ray emission is thermal in origin and is dominated by highly ionized species of O, Ne, Mg, and Si., Comment: 16 pages, 3 figures, Accepted for publication in ApJ Letters

We present ultraviolet and optical spectra of the first bright spot (PA = 29 degrees) on Supernova 1987A's equatorial circumstellar ring taken with the Space Telescope Imaging Spectrograph. We interpret this spot as the emission produced by radiative shocks that occur where the supernova blast wave strikes an inward protrusion of the ring. The observed line widths and intensity ratios indicate the presence of radiative shocks with velocities ranging from 100 to 250 km s^-1 entering dense (> 10^4 cm^-3) gas. These observations, and future observations of the development of the spectra and line profiles, provide a unique opportunity to study the hydrodynamics of radiative shocks., Comment: Accepted by ApJL on Aug. 3, 2000

We describe and model high velocity (~15,000 km/s) Lyman alpha and Balmer alpha emission from supernova remnant 1987A seen in September and October 1997 with the Space Telescope Imaging Spectrograph. Part of this emission comes from a reverse shock located at ~75% of the radius of the inner boundary of the inner circumstellar ring and confined within 30 degrees of the equatorial plane. Departure from axisymmetry in the Lyman alpha and Balmer alpha emission correlates with that seen in nonthermal radio emission and reveals an asymmetry in the circumstellar gas distribution. We also see diffuse high velocity Lyman alpha emission from supernova debris inside the reverse shock that may be due to excitation by nonthermal particles accelerated by the shock., Comment: 10 pages, 3 postscript figures. Submitted to ApJL on Sept 8, 1998; accepted Oct. 20, 1998. Text adjusted to referee's comments

This paper reports on a workshop hosted by the University of Minnesota, March 23-26, 1997. It addressed fundamental dynamical issues associated with the evolution of shell supernova remnants and the relationships between supernova remnants and their environments. The workshop considered, in addition to classical shell SNRs, dynamical issues involving X-ray filled composite remnants and pulsar driven shells, such as that in the Crab Nebula. Approximately 75 participants with wide ranging interests attended the workshop. An even larger community helped through extensive on-line debates prior to the meeting. Each of the several sessions, organized mostly around chronological labels, also addressed some underlying, general physical themes: How are SNR dynamics and structures modified by the character of the CSM and the ISM and vice versa? How are magnetic fields generated in SNRs and how do magnetic fields influence SNRs? Where and how are cosmic-rays (electrons and ions) produced in SNRs and how does their presence influence or reveal SNR dynamics? How does SNR blast energy partition into various components over time and what controls conversion between components? In lieu of a proceedings volume, we present here a synopsis of the workshop in the form of brief summaries of the workshop sessions. The sharpest impressions from the workshop were the crucial and under-appreciated roles that environments have on SNR appearance and dynamics and the critical need for broad-based studies to understand these beautiful, but enigmatic objects. \\, Comment: 54 pages text, no figures, Latex (aasms4.sty). submitted to the PASP

With the High Resolution Spectrograph aboard the Hubble Space Telescope we obtained high signal-to-noise (S/N > 200 - 600 per 17 km/s resolution element) spectra of narrow absorption lines toward the Wolf-Rayet star HD 50896. The ring nebula S308 that surrounds this star is thought to be caused by a pressure-driven bubble bounded by circumstellar gas (most likely from a red supergiant or luminous blue variable progenitor) pushed aside by a strong stellar wind. Our observation has shown for the first time that blueshifted (approximately 70 km/s relative to the star) absorption components of C IV and N V arise in a conduction front between the hot interior of the bubble and the cold shell of swept-up material. These lines set limits on models of the conduction front. Nitrogen in the shell appears to be overabundant by a factor ~10. The P Cygni profiles of N V and C IV are variable, possibly due to a suspected binary companion to HD 50896., Comment: 32 pages, Latex, to appear in the Astrophysical Journal, April, 1997

Analysis of Hubble Space Telescope Band R band images from 1994 to 2009 show that the optical luminosity of SN 1987A has transitioned from being powered by radioactive decay of Ti-44 to energy deposited by X-rays produced as the ejecta interacts with the surrounding material (Larsson et al. 2011, Nature, 474, 484). The B and R band flux from the densest, central parts of the ejecta followed the expected exponential decline until 2001 (about day 5000) when the flux in these bands started increasing, more than doubling by the end of 2009. This increase is the result of heat deposited by X-rays from the shock interaction of the fast-moving outer ejecta with the inner circumstellar ring. In time, the X-rays will penetrate farther into the ejecta, enabling us to analyze the structure and chemistry of the vanished star.

Multiwavelength observations of supernova remnant (SNR) 1987A show that its morphology and luminosity are rapidly changing at X-ray, optical, infrared, and radio wavelengths as the blast wave from the explosion expands into the circumstellar equatorial ring, produced by mass loss from the progenitor star. The observed infrared (IR) radiation arises from the interaction of dust grains that formed in mass outflow with the soft X-ray emitting plasma component of the shocked gas. Spitzer IRS spectra at 5 - 30 microns taken on day 6190 since the explosion show that the emission arises from approx. 1.1 x 10(exp -6) solar mass of silicate grains radiating at a temperature of approx. 180+/-(15-20) K. Subsequent observations on day 7137 show that the IR flux had increased by a factor of 2 while maintaining an almost identical spectral shape. The observed IR-to-X-ray flux ratio (IRX) is consistent with that of a dusty plasma with standard LMC dust abundances. This flux ratio has decreased by a factor of approx. 2 between days 6190 and 7137, providing the first direct observation of the ongoing destruction of dust in an expanding SN blast wave on dynamic time scales. Detailed models consistent with the observed dust temperature, the ionization fluence of the soft X-ray emission component, and the evolution of IRX suggest that the radiating si1icate grains are immersed in a 3.5 x 10(exp 6) K plasma with a density of (0.3 - 1) x 10(exp 4)/cu cm, and have a size distribution that is confined to a narrow range of radii between 0.02 and 0.2 microns. Smaller grains may have been evaporated by the initial UV flash from the supernova.

This grant was intended to support the development of theoretical models needed to interpret and understand the observations by the Hubble Space Telescope and the Chandra X-ray telescope of the rapidly developing remnant of Supernova 1987A. In addition, we carried out a few investigations of related topics. The project was spectacularly successful. The models that we developed provide the definitive framework for predicting and interpreting this phenomenon. Following is a list of publications based on our work. Some of these papers include results of both theoretical modeling supported by this project and also analysis of data supported by the Space Telescope Science Institute and the Chandra X-ray Observatory. We first list papers published in refereed journals, then conference proceedings and book chapters, and also an educational web site.

Following is final report on the study "The Impact of SN1987A with its Circumstellar Ring", which is now complete. In 1994, it was predicted that the blast wave from SN1987A would strike the circumstellar ring in AD1999+/-3, and that the ring would brighten by several hundreds optically as the shock entered the ring. It was also predicted that the emission lines from the shocked ring would have linewidths of a few hundred km/s. Today, we see a "hot spot" on the ring that first appeared in Hubble Space Telescope (HST) images in 1995 and has doubled in brightness between August 1997 and February 1998. Moreover, spectra from STIS show that the emission lines from the hot spot have widths and blue shifts of order 200 km/s, just as we predicted in 1994. Our guess of 1999+/-3 for the impact time was lucky, because we assumed that the hot gas in the bubble between the supernova and the ring had a low density (approximately 10 cm(exp -3)). But it was pointed out that the ROSAT observations of soft X-rays from SN 1987A implied that the intervening gas had a higher density, (approximately 100 cm(exp -3)), which would delay the impact until ca. 2007. Others developed hydrodynamic models to fit the ROSAT X-ray emission spectrum and came to the same conclusion. But these models were oversimplified in that they assumed that the ring was perfectly round. Now we see clearly that the hot spot is a peninsula that protrudes inward from the ring -- the first spot on the ring to be struck by the blast wave. No doubt there are other protrusions on the ring, which we may expect to light up in the next few years until they finally merge to set the entire ring ablaze, probably within 5 - 10 years. It was also predicted that the X-ray emission must be accompanied by optical and ultraviolet emission from atoms in the supernova debris and in the circumstellar gas that cross the reverse shock and the blast wave, respectively, and that the Lyman-alpha and NV emission lines should be bright enough to see with the STIS.