Source author record

Stephen P. Reynolds

Stephen P. Reynolds appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Evidence for a Dense, Inhomogeneous Circumstellar Medium in the Type Ia SNR 0519-69.0

We perform an expansion study of the Balmer dominated outer shock of the SNR 0519$-$69.0 in the LMC by using a combination of new HST WFC3 imagery obtained in 2020 and archival ACS images from 2010 and 2011. Thanks to the very long time baseline, our proper motion measurements are of unprecedented accuracy. We find a wide range of shock velocities, with the fastest shocks averaging 5280 km/s and the slowest grouping of shocks averaging just 1670 km/s. We compare the H_alpha images from HST with X-ray images from Chandra and mid-IR images from Spitzer, finding a clear anti-correlation between the brightness of the remnant in a particular location and the velocity of the blast wave at that location, supporting the idea that the bright knots of X-ray and IR emission result from an interaction with a dense inhomogeneous circumstellar medium. We find no evidence for X-ray emission, thermal or nonthermal, associated with the fastest shocks, as expected if the fastest velocities are the result of the blast wave encountering the lower density ambient medium of the LMC. We derive an age of the remnant of 670 +/- 70 yr, consistent with results derived from previous investigations.

preprint2022arXiv

The Eel Pulsar Wind Nebula: a PeVatron-Candidate Origin for HAWC J1826-128 and HESS J1826-130

HAWC J1826-128 is one of the brightest Galactic TeV gamma-ray sources detected by the High Altitude Water Cherenkov (HAWC) Observatory, with photon energies extending up to nearly $\sim$100 TeV. This HAWC source spatially coincides with the H.E.S.S. TeV source HESS J1826-130 and the "Eel" pulsar wind nebula (PWN), which is associated with the GeV pulsar PSR J1826-1256. In the X-ray band, Chandra and XMM-Newton revealed that the Eel PWN is composed of both a compact nebula ($\sim$15") and diffuse X-ray emission ($\sim$6'$\times$2') extending away from the pulsar. Our NuSTAR observation detected hard X-ray emission from the compact PWN up to $\sim$20 keV and evidence of the synchrotron burn-off effect. In addition to the spatial coincidence between HESS J1826-130 and the diffuse X-ray PWN, our multi-wavelength spectral energy distribution (SED) analysis using X-ray and gamma-ray data establishes a leptonic origin of the TeV emission associated with the Eel PWN. Furthermore, our evolutionary PWN SED model suggests (1) a low PWN B-field of $\sim$1 $μ$G, (2) a significantly younger pulsar age ($t \sim5.7$ kyr) than the characteristic age ($τ= 14.4$ kyr) and (3) a maximum electron energy of $E_{max} = 2$ PeV. The low B-field as well as the putative supersonic motion of the pulsar may account for the asymmetric morphology of the diffuse X-ray emission. Our results suggest that the Eel PWN may be a leptonic PeVatron particle accelerator powered by the $\sim$6-kyr-old pulsar PSR J1826-1256 with a spin-down power of $3.6 \times 10^{36}$ erg s$^{-1}$.

preprint2021arXiv

Type Ia Supernova Models: Asymmetric Remnants and Supernova Remnant G1.9+0.3

The youngest Galactic supernova remnant G1.9+0.3, probably the result of a Type Ia supernova, shows surprising anomalies in the distribution of its ejecta in space and velocity. In particular, high-velocity shocked iron is seen in several locations far from the remnant center, in some cases beyond prominent silicon and sulfur emission. These asymmetries strongly suggest a highly asymmetric explosion. We present high-resolution hydrodynamic simulations in two and three dimensions of the evolution from ages of 100 seconds to hundreds of years of two asymmetric Type Ia models, expanding into a uniform medium. At the age of G1.9+0.3 (about 100 years), our 2D model shows almost no iron shocked to become visible in X-rays. Only in a much higher-density environment could significant iron be shocked, at which time the model's expansion speed is completely inconsistent with the observations of G1.9+0.3. Our 3D model, evolving the most asymmetric of a suite of Type Ia SN models from Seitenzahl et al.~(2013), shows some features resembling G1.9+0.3. We characterize its evolution with images of composition in three classes: C and O, intermediate-mass elements (IMEs), and iron-group elements (IGEs). From ages of 13 to 1800 years, we follow the evolution of the highly asymmetric initial remnant as the explosion asymmetries decrease in relative strength to be replaced by asymmetries due to evolutionary hydrodynamic instabilities. At an age of about 100 years, our 3D model has comparable shocked masses of C+O, IMEs, and IGEs, with about 0.03 $M_\odot$ each. Evolutionary changes appear to be rapid enough that continued monitoring with the Chandra X-ray Observatory may show significant variations.

preprint2020arXiv

Fast Blast Wave and Ejecta in the Young Core-Collapse Supernova Remnant MSH 15-52/RCW 89

One of the youngest known remnants of a core-collapse supernova (SN) in our Galaxy is G320.4$-$1.2/MSH 15-52 containing an energetic pulsar with a very short (1700 yr) spindown age and likely produced by a stripped-envelope SN Ibc. Bright X-ray and radio emission north of the pulsar overlaps with an H$α$ nebula RCW 89. The bright X-rays there have a highly unusual and quite puzzling morphology, consisting of both very compact thermally emitting knots and much more diffuse emission of nonthermal origin. We report new X-ray observations of RCW 89 in 2017 and 2018 with Chandra that allowed us to measure the motions of many knots and filaments on decade-long time baselines. We identify a fast blast wave with a velocity of $(4000 \pm 500)d_{5.2}$ km/s ($d_{5.2}$ is the distance in units of 5.2 kpc) with a purely nonthermal spectrum, and without any radio counterpart. Many compact X-ray emission knots are moving vary fast, with velocities as high as 5000 km/s, predominantly radially away from the pulsar. Their spectra show that they are Ne- and Mg-rich heavy-element SN ejecta. They have been significantly decelerated upon their recent impact with the dense ambient medium north of the pulsar. We see fast evolution in brightness and morphology of knots in just a few years. Ejecta knots in RCW 89 resemble those seen in Cas A at optical wavelengths in terms of their initial velocities and densities. They might have the same origin, still not understood but presumably related to stripped-envelope SN explosions themselves.

preprint2019arXiv

On the Expansion, Age, and Origin of the Puzzling Shell/Pulsar Wind Nebula G310.6-1.6

We present a 142-ks Chandra observation of the enigmatic combination supernova remnant G310.6-1.6 consisting of a bright pulsar-wind nebula driven by an energetic pulsar, surrounded by a highly circular, very faint shell with a featureless, probably synchrotron, spectrum. Comparison with an observation 6 years earlier shows no measurable expansion of the shell, though some features in the pulsar-wind nebula have moved. We find an expansion age of at least 2500 yr, implying a current shock velocity less than about 1000 km/s. We place severe upper limits on thermal emission from the shell; if the shell locates the blast wave, a Sedov interpretation would require the remnant to be very young, about 1000 yr, and to have resulted from a dramatically sub-energetic supernova, ejecting << 0.02 M_sun with energy E < 3 x 10^47 erg. Even a merger-induced collapse of a white dwarf to a neutron star, with a low-energy explosion, is unlikely to produce such an event. Other explanations seem equally unlikely.

preprint2016arXiv

A Compact X-ray Source in the Radio Pulsar-Wind Nebula G141.2+5.0

We report the results of a 50 ks Chandra observation of the recently discovered radio object G141.2+5.0, presumed to be a pulsar-wind nebula. We find a moderately bright unresolved X-ray source which we designate CXOU J033712.8 615302 coincident with the central peak radio emission. An absorbed power-law fit to the 241 counts describes the data well, with absorbing column $N_H = 6.7 (4.0, 9.7) \times 10^{21}$ cm$^{-2}$ and photon index $Γ= 1.8 (1.4, 2.2)$. For a distance of 4 kpc, the unabsorbed luminosity between 0.5 and 8 keV is $ 1.7^{+0.4}_{-0.3} \times 10^{32}$ erg s$^{-1}$ (90\% confidence intervals). Both $L_X$ and $Γ$ are quite typical of pulsars in PWNe. No extended emission is seen; we estimate a conservative $3 σ$ upper limit to the surface brightness of any X-ray PWN near the point source to be $3 \times 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ arcsec$^{-2}$ between 0.5 and 8 keV, assuming the same spectrum as the point source; for a nebula of diameter $13"$, the flux limit is 6\% of the flux of the point source. The steep radio spectrum of the PWN ($α\sim -0.7$), if continued to the X-ray without a break, predicts $L_X\ \rm{(nebula)} \sim 1 \times 10^{33}$ erg s$^{-1}$, so additional spectral steepening between radio and X-rays is required, as is true of all known PWNe. The high Galactic latitude gives a $z$-distance of 350 pc above the Galactic plane, quite unusual for a Population I object.

preprint2016arXiv

An X-ray and Radio Study of the Varying Expansion Velocities in Tycho's Supernova Remnant

We present newly obtained X-ray and radio observations of Tycho's supernova remnant using {\it Chandra} and the Karl G. Jansky Very Large Array in 2015 and 2013/14, respectively. When combined with earlier epoch observations by these instruments, we now have time baselines for expansion measurements of the remnant of 12-15 year in the X-rays and 30 year in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. We find, consistent with earlier measurements, a clear gradient in the expansion velocity of the remnant, despite its round shape. The proper motions on the western and southwestern sides of the remnant are about a factor of two higher than those in the east and northeast. We showed in an earlier work that this is related to an offset of the explosion site from the geometric center of the remnant due to a density gradient in the ISM, and using our refined measurements reported here, we find that this offset is $\sim 23"$ towards the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

preprint2016arXiv

G11.2-0.3: The Young Remnant of a Stripped-Envelope Supernova

We present results of a 400-ks Chandra observation of the young shell supernova remnant (SNR) G11.2-0.3, containing a pulsar and pulsar-wind nebula (PWN). We measure a mean expansion rate for the shell since 2000 of 0.0277+/-0.0018% per yr, implying an age between 1400 and 2400 yr, and making G11.2-0.3 one of the youngest core-collapse SNRs in the Galaxy. However, we find very high absorption ($A_V \sim 16^m \pm 2^m$), confirming near-IR determinations and ruling out a claimed association with the possible historical SN of 386 CE. The PWN shows strong jets and a faint torus within a larger, more diffuse region of radio emission and nonthermal X-rays. Central soft thermal X-ray emission is anticorrelated with the PWN; that, and more detailed morphological evidence, indicates that the reverse shock has already reheated all ejecta and compressed the PWN. The pulsar characteristic energy-loss timescale is well in excess of the remnant age, and we suggest that the bright jets have been produced since the recompression. The relatively pronounced shell and diffuse hard X-ray emission in the interior, enhanced at the inner edge of the shell, indicate that the immediate circumstellar medium into which G11.2-0.3 is expanding was quite anisotropic. We propose a possible origin for G11.2-0.3 in a stripped-envelope progenitor that had lost almost all its envelope mass, in an anisotropic wind or due to binary interaction, leaving a compact core whose fast winds swept previously lost mass into a dense irregular shell, and which exploded as a Type cIIb or Ibc supernova.

preprint2016arXiv

The distribution of radioactive $^{44}$Ti in Cassiopeia A

The distribution of elements produced in the inner-most layers of a supernova explosion is a key diagnostic for studying the collapse of massive stars. Here we present the results of a 2.4 Ms \textit{NuSTAR} observing campaign aimed at studying the supernova remnant Cassiopeia A (Cas A). We perform spatially-resolved spectroscopic analyses of the $^{44}$Ti ejecta which we use to determine the Doppler shift and thus the three-dimensional (3D) velocities of the $^{44}$Ti ejecta. We find an initial $^{44}$Ti mass of 1.54 $\pm$ 0.21 $\times 10^{-4}$ M$_{\odot}$ which has a present day average momentum direction of 340$^{\circ}$ $\pm$ 15$^{\circ}$ projected on to the plane of the sky (measured clockwise from Celestial North) and tilted by 58$^{\circ}$ $\pm$ 20$^{\circ}$ into the plane of the sky away from the observer, roughly opposite to the inferred direction of motion of the central compact object. We find some $^{44}$Ti ejecta that are clearly interior to the reverse shock and some that are clearly exterior to the reverse shock. Where we observe $^{44}$Ti ejecta exterior to the reverse shock we also see shock-heated iron; however, there are regions where we see iron but do not observe $^{44}$Ti. This suggests that the local conditions of the supernova shock during explosive nucleosynthesis varied enough to suppress the production of $^{44}$Ti in some regions by at least a factor of two, even in regions that are assumed to be the result of processes like $α$-rich freezeout that should produce both iron and titanium.

preprint2015arXiv

A Spatially Resolved Study of the Synchrotron Emission and Titanium in Tycho's Supernova Remnant with NuSTAR

We report results from deep observations (~750 ks) of Tycho's supernova remnant (SNR) with NuSTAR. Using these data, we produce narrow-band images over several energy bands to identify the regions producing the hardest X-rays and to search for radioactive decay line emission from 44Ti. We find that the hardest (>10 keV) X-rays are concentrated in the southwest of Tycho, where recent Chandra observations have revealed high emissivity "stripes" associated with particles accelerated to the knee of the cosmic-ray spectrum. We do not find evidence of 44Ti, and we set limits on its presence and distribution within the SNR. These limits correspond to a upper-limit 44Ti mass of M44 < 2.4x10^-4 M_sun for a distance of 2.3 kpc. We perform spatially resolved spectroscopic analysis of sixty-six regions across Tycho. We map the best-fit rolloff frequency of the hard X-ray spectra, and we compare these results to measurements of the shock expansion and ambient density. We find that the highest energy electrons are accelerated at the lowest densities and in the fastest shocks, with a steep dependence of the roll-off frequency with shock velocity. Such a dependence is predicted by models where the maximum energy of accelerated electrons is limited by the age of the SNR rather than by synchrotron losses, but this scenario requires far lower magnetic field strengths than those derived from observations in Tycho. One way to reconcile these discrepant findings is through shock obliquity effects, and future observational work is necessary to explore the role of obliquity in the particle acceleration process.

preprint2015arXiv

Energy Dependence of Synchrotron X-Ray Rims in Tycho's Supernova Remnant

Several young supernova remnants exhibit thin X-ray bright rims of synchrotron radiation at their forward shocks. Thin rims require strong magnetic field amplification beyond simple shock compression if rim widths are only limited by electron energy losses. But, magnetic field damping behind the shock could produce similarly thin rims with less extreme field amplification. Variation of rim width with energy may thus discriminate between competing influences on rim widths. We measured rim widths around Tycho's supernova remnant in 5 energy bands using an archival 750 ks Chandra observation. Rims narrow with increasing energy and are well described by either loss-limited or damped scenarios, so X-ray rim width-energy dependence does not uniquely specify a model. But, radio counterparts to thin rims are not loss-limited and better reflect magnetic field structure. Joint radio and X-ray modeling favors magnetic damping in Tycho's SNR with damping lengths ~1--5% of remnant radius and magnetic field strengths ~50--400 $μ$G assuming Bohm diffusion. X-ray rim widths are ~1% of remnant radius, somewhat smaller than inferred damping lengths. Electron energy losses are important in all models of X-ray rims, suggesting that the distinction between loss-limited and damped models is blurred in soft X-rays. All loss-limited and damping models require magnetic fields $\gtrsim$ 20 $μ$G, affirming the necessity of magnetic field amplification beyond simple compression.

preprint2015arXiv

Locating the most energetic electrons in Cassiopeia A

We present deep ($>$2.4 Ms) observations of the Cassiopeia A supernova remnant with {\it NuSTAR}, which operates in the 3--79 keV bandpass and is the first instrument capable of spatially resolving the remnant above 15 keV. We find that the emission is not entirely dominated by the forward shock nor by a smooth "bright ring" at the reverse shock. Instead we find that the $>$15 keV emission is dominated by knots near the center of the remnant and dimmer filaments near the remnant's outer rim. These regions are fit with unbroken power-laws in the 15--50 keV bandpass, though the central knots have a steeper ($Γ\sim -3.35$) spectrum than the outer filaments ($Γ\sim -3.06$). We argue this difference implies that the central knots are located in the 3-D interior of the remnant rather than at the outer rim of the remnant and seen in the center due to projection effects. The morphology of $>$15 keV emission does not follow that of the radio emission nor that of the low energy ($<$12 keV) X-rays, leaving the origin of the $>$15 keV emission as an open mystery. Even at the forward shock front we find less steepening of the spectrum than expected from an exponentially cut off electron distribution with a single cutoff energy. Finally, we find that the GeV emission is not associated with the bright features in the {\it NuSTAR} band while the TeV emission may be, suggesting that both hadronic and leptonic emission mechanisms may be at work.

preprint2015arXiv

Second Epoch Hubble Space Telescope Observations of Kepler's Supernova Remnant: The Proper Motions of Balmer Filaments

We report on the proper motions of Balmer-dominated filaments in Kepler's supernova remnant using high resolution images obtained with the Hubble Space Telescope at two epochs separated by about 10 years. We use the improved proper motion measurements and revised values of shock velocities to derive a distance to Kepler of 5.1 [+0.8, -0.7] kpc. The main shock around the northern rim of the remnant has a typical speed of 1690 km/s and is encountering material with densities of about 8 cm^-3. We find evidence for the variation of shock properties over small spatial scales, including differences in the driving pressures as the shock wraps around a curved cloud surface. We find that the Balmer filaments ahead of the ejecta knot on the northwest boundary of the remnant are becoming fainter and more diffuse. We also find that the Balmer filaments associated with circumstellar material in the interior regions of the remnant are due to shocks with significantly lower velocities and that the brightness variations among these filaments trace the density distribution of the material, which may have a disk-like geometry.

preprint2014arXiv

High-Energy X-ray Imaging of the Pulsar Wind Nebula MSH~15-52: Constraints on Particle Acceleration and Transport

We present the first images of the pulsar wind nebula (PWN) MSH 15-52 in the hard X-ray band (>8 keV), as measured with the Nuclear Spectroscopic Telescope Array (NuSTAR). Overall, the morphology of the PWN as measured by NuSTAR in the 3-7 keV band is similar to that seen in Chandra high-resolution imaging. However, the spatial extent decreases with energy, which we attribute to synchrotron energy losses as the particles move away from the shock. The hard-band maps show a relative deficit of counts in the northern region towards the RCW 89 thermal remnant, with significant asymmetry. We find that the integrated PWN spectra measured with NuSTAR and Chandra suggest that there is a spectral break at 6 keV which may be explained by a break in the synchrotron-emitting electron distribution at ~200 TeV and/or imperfect cross calibration. We also measure spatially resolved spectra, showing that the spectrum of the PWN softens away from the central pulsar B1509-58, and that there exists a roughly sinusoidal variation of spectral hardness in the azimuthal direction. We discuss the results using particle flow models. We find non-monotonic structure in the variation with distance of spectral hardness within 50" of the pulsar moving in the jet direction, which may imply particle and magnetic-field compression by magnetic hoop stress as previously suggested for this source. We also present 2-D maps of spectral parameters and find an interesting shell-like structure in the NH map. We discuss possible origins of the shell-like structure and their implications.

preprint2014arXiv

Magnetic-Field Amplification in the Thin X-ray Rims of SN1006

Several young supernova remnants (SNRs), including SN1006, emit synchrotron X-rays in narrow filaments, hereafter thin rims, along their periphery. The widths of these rims imply 50 to 100 $μ$G fields in the region immediately behind the shock, far larger than expected for the interstellar medium compressed by unmodified shocks, assuming electron radiative losses limit rim widths. However, magnetic-field damping could also produce thin rims. Here we review the literature on rim width calculations, summarizing the case for magnetic-field amplification. We extend these calculations to include an arbitrary power-law dependence of the diffusion coefficient on energy, $D \propto E^μ$. Loss-limited rim widths should shrink with increasing photon energy, while magnetic-damping models predict widths almost independent of photon energy. We use these results to analyze Chandra observations of SN 1006, in particular the southwest limb. We parameterize the full widths at half maximum (FWHM) in terms of energy as FWHM $\propto E^{m_E}_γ$. Filament widths in SN1006 decrease with energy; $m_E \sim -0.3$ to $-0.8$, implying magnetic field amplification by factors of 10 to 50, above the factor of 4 expected in strong unmodified shocks. For SN 1006, the rapid shrinkage rules out magnetic damping models. It also favors short mean free paths (small diffusion coefficients) and strong dependence of $D$ on energy ($μ\ge 1$).

preprint2014arXiv

NuSTAR study of Hard X-Ray Morphology and Spectroscopy of PWN G21.5-0.9

We present NuSTAR high energy X-ray observations of the pulsar wind nebula (PWN)/supernova remnant G21.5-0.9. We detect integrated emission from the nebula up to ~40 keV, and resolve individual spatial features over a broad X-ray band for the first time. The morphology seen by NuSTAR agrees well with that seen by XMM-Newton and Chandra below 10 keV. At high energies NuSTAR clearly detects non-thermal emission up to ~20 keV that extends along the eastern and northern rim of the supernova shell. The broadband images clearly demonstrate that X-ray emission from the North Spur and Eastern Limb results predominantly from non-thermal processes. We detect a break in the spatially integrated X-ray spectrum at ~9 keV that cannot be reproduced by current SED models, implying either a more complex electron injection spectrum or an additional process such as diffusion compared to what has been considered in previous work. We use spatially resolved maps to derive an energy-dependent cooling length scale, $L(E) \propto E^{m}$ with $m = -0.21 \pm 0.01$. We find this to be inconsistent with the model for the morphological evolution with energy described by Kennel & Coroniti (1984). This value, along with the observed steepening in power-law index between radio and X-ray, can be quantitatively explained as an energy-loss spectral break in the simple scaling model of Reynolds (2009), assuming particle advection dominates over diffusion. This interpretation requires a substantial departure from spherical magnetohydrodynamic (MHD), magnetic-flux-conserving outflow, most plausibly in the form of turbulent magnetic-field amplification.

preprint2014arXiv

Spitzer Observations of the Type Ia Supernova Remnant N103B: Kepler's Older Cousin?

We report results from Spitzer observations of SNR 0509-68.7, also known as N103B, a young Type Ia supernova remnant in the Large Magellanic Cloud that shows interaction with a dense medium in its western hemisphere. Our images show that N103B has strong IR emission from warm dust in the post-shock environment. The post-shock gas density we derive, 45 cm$^{-3}$, is much higher than in other Type Ia remnants in the LMC, though a lack of spatial resolution may bias measurements towards regions of higher than average density. This density is similar to that in Kepler's SNR, a Type Ia interacting with a circumstellar medium. Optical images show H$α$ emission along the entire periphery of the western portion of the shock, with [O III] and [S II] lines emitted from a few dense clumps of material where the shock has become radiative. The dust is silicate in nature, though standard silicate dust models fail to reproduce the "18 $μ$m" silicate feature that peaks instead at 17.3 $μ$m. We propose that the dense material is circumstellar material lost from the progenitor system, as with Kepler. If the CSM interpretation is correct, this remnant would become the second member, along with Kepler, of a class of Type Ia remnants characterized by interaction with a dense CSM hundreds of years post-explosion. A lack of N enhancement eliminates symbiotic AGB progenitors. The white dwarf companion must have been relatively unevolved at the time of the explosion.

preprint2014arXiv

The Hard X-Ray View of the Young Supernova Remnant G1.9+0.3

NuSTAR observed G1.9+0.3, the youngest known supernova remnant in the Milky Way, for 350 ks and detected emission up to $\sim$30 keV. The remnant's X-ray morphology does not change significantly across the energy range from 3 to 20 keV. A combined fit between NuSTAR and CHANDRA shows that the spectrum steepens with energy. The spectral shape can be well fitted with synchrotron emission from a power-law electron energy distribution with an exponential cutoff with no additional features. It can also be described by a purely phenomenological model such as a broken power-law or a power-law with an exponential cutoff, though these descriptions lack physical motivation. Using a fixed radio flux at 1 GHz of 1.17 Jy for the synchrotron model, we get a column density of N$_{\rm H}$ = $(7.23\pm0.07) \times 10^{22}$ cm$^{-2}$, a spectral index of $α=0.633\pm0.003$, and a roll-off frequency of $ν_{\rm rolloff}=(3.07\pm0.18) \times 10^{17}$ Hz. This can be explained by particle acceleration, to a maximum energy set by the finite remnant age, in a magnetic field of about 10 $μ$G, for which our roll-off implies a maximum energy of about 100 TeV for both electrons and ions. Much higher magnetic-field strengths would produce an electron spectrum that was cut off by radiative losses, giving a much higher roll-off frequency that is independent of magnetic-field strength. In this case, ions could be accelerated to much higher energies. A search for $^{44}$Ti emission in the 67.9 keV line results in an upper limit of $1.5 \times 10^{-5}$ $\,\mathrm{ph}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$ assuming a line width of 4.0 keV (1 sigma).

preprint2014arXiv

X-ray Spectral Components in the Afterglow of GRB 130925A

We have identified spectral features in the late-time X-ray afterglow of the unusually long, slow-decaying GRB 130925A using NuSTAR, Swift-XRT, and Chandra. A spectral component in addition to an absorbed power-law is required at $>4σ$ significance, and its spectral shape varies between two observation epochs at $2\times10^5$ and $10^6$ seconds after the burst. Several models can fit this additional component, each with very different physical implications. A broad, resolved Gaussian absorption feature of several keV width improves the fit, but it is poorly constrained in the second epoch. An additive black body or second power-law component provide better fits. Both are challenging to interpret: the blackbody radius is near the scale of a compact remnant ($10^8$ cm), while the second powerlaw component requires an unobserved high-energy cutoff in order to be consistent with the non-detection by Fermi-LAT.

preprint2013arXiv

Azimuthal Density Variations Around the Rim of Tycho's Supernova Remnant

{\it Spitzer} images of Tycho's supernova remnant in the mid-infrared reveal limb-brightened emission from the entire periphery of the shell and faint filamentary structures in the interior. As with other young remnants, this emission is produced by dust grains, warmed to $\sim 100$ K in the post-shock environment by collisions with energetic electrons and ions. The ratio of the 70 to 24 $μ$m fluxes is a diagnostic of the dust temperature, which in turn is a sensitive function of the plasma density. We find significant variations in the 70/24 flux ratio around the periphery of Tycho's forward shock, implying order-of-magnitude variations in density. While some of these are likely localized interactions with dense clumps of the interstellar medium, we find an overall gradient in the ambient density surrounding Tycho, with densities 3-10 times higher in the NE than in the SW. This large density gradient is qualitatively consistent with the variations in the proper motion of the shock observed in radio and X-ray studies. Overall, the mean ISM density around Tycho is quite low ($\sim 0.1-0.2$ cm$^{-3}$), consistent with the lack of thermal X-ray emission observed at the forward shock. We perform two-dimensional hydrodynamic simulations of a Type Ia SN expanding into a density gradient in the ISM, and find that the overall round shape of the remnant is still easily acheivable, even for explosions into significant gradients. However, this leads to an offset of the center of the explosion from the geometric center of the remnant of up to 20%, although lower values of 10% are preferred. The best match with hydrodynamical simulations is achieved if Tycho is located at a large (3-4 kpc) distance in a medium with a mean preshock density of $\sim 0.2$ cm$^{-3}$. Such preshock densities are obtained for highly ($\ga 50$%) porous ISM grains.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2012arXiv

Dust in a Type Ia Supernova Progenitor: Spitzer Spectroscopy of Kepler's Supernova Remnant

Characterization of the relatively poorly-understood progenitor systems of Type Ia supernovae is of great importance in astrophysics, particularly given the important cosmological role that these supernovae play. Kepler's Supernova Remnant, the result of a Type Ia supernova, shows evidence for an interaction with a dense circumstellar medium (CSM), suggesting a single-degenerate progenitor system. We present 7.5-38 $μ$m infrared (IR) spectra of the remnant, obtained with the {\it Spitzer Space Telescope}, dominated by emission from warm dust. Broad spectral features at 10 and 18 $μ$m, consistent with various silicate particles, are seen throughout. These silicates were likely formed in the stellar outflow from the progenitor system during the AGB stage of evolution, and imply an oxygen-rich chemistry. In addition to silicate dust, a second component, possibly carbonaceous dust, is necessary to account for the short-wavelength IRS and IRAC data. This could imply a mixed chemistry in the atmosphere of the progenitor system. However, non-spherical metallic iron inclusions within silicate grains provide an alternative solution. Models of collisionally-heated dust emission from fast shocks ($>$ 1000 km s$^{-1}$) propagating into the CSM can reproduce the majority of the emission associated with non-radiative filaments, where dust temperatures are $\sim 80-100$ K, but fail to account for the highest temperatures detected, in excess of 150 K. We find that slower shocks (a few hundred km s$^{-1}$) into moderate density material ($n_{0} \sim 50-250$ cm$^{-3}$) are the only viable source of heating for this hottest dust. We confirm the finding of an overall density gradient, with densities in the north being an order of magnitude greater than those in the south.

preprint2012arXiv

The First Reported Infrared Emission from the SN 1006 Remnant

We report results of infrared imaging and spectroscopic observations of the SN 1006 remnant, carried out with the Spitzer Space Telescope. The 24 micron image from MIPS clearly shows faint filamentary emission along the northwest rim of the remnant shell, nearly coincident with the Balmer filaments that delineate the present position of the expanding shock. The 24 micron emission traces the Balmer filaments almost perfectly, but lies a few arcsec within, indicating an origin in interstellar dust heated by the shock. Subsequent decline in the IR behind the shock is presumably due largely to grain destruction through sputtering. The emission drops far more rapidly than current models predict, however, even for a higher proportion of small grains than would be found closer to the Galactic plane. The rapid drop may result in part from a grain density that has always been lower -- a relic effect from an earlier epoch when the shock was encountering a lower density -- but higher grain destruction rates still seem to be required. Spectra from three positions along the NW filament from the IRS instrument all show only a featureless continuum, consistent with thermal emission from warm dust. The dust-to-gas mass ratio in the pre-shock interstellar medium is lower than that expected for the Galactic ISM -- as has also been observed in the analysis of IR emission from other SNRs but whose cause remains unclear. As with other SN Ia remnants, SN 1006 shows no evidence for dust grain formation in the supernova ejecta.

preprint2012arXiv

X-ray Emission from Strongly Asymmetric Circumstellar Material in the Remnant of Kepler's Supernova

Kepler's supernova remnant resulted from a thermonuclear explosion, but is interacting with circumstellar material (CSM) lost from the progenitor system. We describe a statistical technique for isolating X-ray emission due to CSM from that due to shocked ejecta. Shocked CSM coincides well in position with 24 $μ$m emission seen by {\sl Spitzer}. We find most CSM to be distributed along the bright north rim, but substantial concentrations are also found projected against the center of the remnant, roughly along a diameter with position angle $\sim 100^\circ$. We interpret this as evidence for a disk distribution of CSM before the SN, with the line of sight to the observer roughly in the disk plane. We present 2-D hydrodynamic simulations of this scenario, in qualitative agreement with the observed CSM morphology. Our observations require Kepler to have originated in a close binary system with an AGB star companion.

preprint2012arXiv

X-ray Proper Motions and Shock Speeds along the Northwest Rim of SN 1006

We report the results of an X-ray proper motion measurement for the NW rim of SN1006, carried out by comparing Chandra observations from 2001 and 2012. The NW limb has predominantly thermal X-ray emission, and it is the only location in SN1006 with significant optical emission: a thin, Balmer-dominated filament. For most of the NW rim, the proper motion is about 0.30 arcsec/yr, essentially the same as has been measured from the H-alpha filament. Isolated regions of the NW limb are dominated by nonthermal emission, and here the proper motion is much higher, 0.49 arcsec/yr, close to the value measured in X-rays along the much brighter NE limb, where the X-rays are overwhelmingly nonthermal. At the 2.2 kpc distance to SN1006, the proper motions imply shock velocities of about 3000 km/s and 5000 km/s in the thermal and nonthermal regions, respectively. A lower velocity behind the H-alpha filament is consistent with the picture that SN1006 is encountering denser gas in the NW, as is also suggested by its overall morphology. In the thermally-dominated portion of the X-ray shell, we also see an offset in the radial profiles at different energies; the 0.5-0.6 keV peak dominated by O VII is closer to the shock front than that of the 0.8-3 keV emission--due to the longer times for heavier elements to reach ionization states where they produce strong X-ray emission.

preprint2011arXiv

Dusty Blastwaves of Two Young LMC Supernova Remnants: Constraints on Postshock Compression

We present results from mid-IR spectroscopic observations of two young supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) done with the {\it Spitzer Space Telescope}. We imaged SNRs B0509-67.5 and B0519-69.0 with {\it Spitzer} in 2005, and follow-up spectroscopy presented here confirms the presence of warm, shock heated dust, with no lines present in the spectrum. We use model fits to {\it Spitzer} IRS data to estimate the density of the postshock gas. Both remnants show asymmetries in the infrared images, and we interpret bright spots as places where the forward shock is running into material that is several times denser than elsewhere. The densities we infer for these objects depend on the grain omposition assumed, and we explore the effects of differing grain porosity on the model fits. We also analyze archival {\it XMM-Newton} RGS spectroscopic data, where both SNRs show strong lines of both Fe and Si, coming from ejecta, as well as strong O lines, which may come from ejecta or shocked ambient medium. We use model fits to IRS spectra to predict X-ray O line strengths for various grain models and values of the shock compression ratio. For 0509-67.5, we find that compact (solid) grain models require nearly all O lines in X-ray spectra to originate in reverse-shocked ejecta. Porous dust grains would lower the strength of ejecta lines relative to those arising in the shocked ambient medium. In 0519-69.0, we find significant evidence for a higher than standard compression ratio of 12, implying efficient cosmic-ray acceleration by the blast wave. A compact grain model is favored over porous grain models. We find that the dust-to-gas mass ratio of the ambient medium is significantly lower than what is expected in the ISM.

preprint2011arXiv

Pulsar Wind Nebulae with Thick Toroidal Structure

We investigate a class of pulsar wind nebulae that show synchrotron emission from a thick toroidal structure. The best studied such object is the small radio and X-ray nebula around the Vela pulsar, which can be interpreted as the result of interaction of a mildly supersonic inward flow with the recent pulsar wind. Such a flow near the center of a supernova remnant can be produced in a transient phase when the reverse shock reaches the center of the remnant. Other nebulae with a thick toroidal structure are G106.6+2.9 and G76.9+1.0. Their structure contrasts with young pulsar nebulae like the Crab Nebula and 3C 38, which show a more chaotic, filamentary structure in the synchrotron emission. In both situations, a torus-jet structure is present where the pulsar wind passes through a termination shock, indicating the flow is initially toroidal. We suggest that the difference is due to the Rayleigh-Taylor instability that operates when the outer boundary of the nebula is accelerating into freely expanding supernova ejecta. The instability gives rise to mixing in the Crab and related objects, but is not present in the nebulae with thick toroidal regions.

preprint2011arXiv

RCW 86: A Type Ia Supernova in a Wind-Blown Bubble

We report results from a multi-wavelength analysis of the Galactic SNR RCW 86, the proposed remnant of the supernova of 185 A.D. We report new infrared observations from {\it Spitzer} and {\it WISE}, where the entire shell is detected at 24 and 22 $μ$m. We fit the infrared flux ratios with models of collisionally heated ambient dust, finding post-shock gas densities in the non-radiative shocks of 2.4 and 2.0 cm$^{-3}$ in the SW and NW portions of the remnant, respectively. The Balmer-dominated shocks around the periphery of the shell, large amount of iron in the X-ray emitting ejecta, and lack of a compact remnant support a Type Ia origin for this remnant. From hydrodynamic simulations, the observed characteristics of RCW 86 are successfully reproduced by an off-center explosion in a low-density cavity carved by the progenitor system. This would make RCW 86 the first known case of a Type Ia supernova in a wind-blown bubble. The fast shocks ($> 3000$ km s$^{-1}$) observed in the NE are propagating in the low-density bubble, where the shock is just beginning to encounter the shell, while the slower shocks elsewhere have already encountered the bubble wall. The diffuse nature of the synchrotron emission in the SW and NW is due to electrons that were accelerated early in the lifetime of the remnant, when the shock was still in the bubble. Electrons in a bubble could produce gamma-rays by inverse-Compton scattering. The wind-blown bubble scenario requires a single-degenerate progenitor, which should leave behind a companion star.

preprint2010arXiv

Radioactive Scandium in the Youngest Galactic Supernova Remnant G1.9+0.3

We report the discovery of thermal X-ray emission from the youngest Galactic supernova remnant G1.9+0.3, from a 237-ks Chandra observation. We detect strong K-shell lines of Si, S, Ar, Ca, and Fe. In addition, we detect a 4.1 keV line with 99.971% confidence which we attribute to 44Sc, produced by electron capture from 44Ti. Combining the data with our earlier Chandra observation allows us to detect the line in two regions independently. For a remnant age of 100 yr, our measured total line strength indicates synthesis of $(1 - 7) \times 10^{-5}$ solar masses of 44Ti, in the range predicted for both Type Ia and core-collapse supernovae, but somewhat smaller than the $2 \times 10^{-4}$ solar masses reported for Cas A. The line spectrum indicates supersolar abundances. The Fe emission has a width of about 28,000 km/s, consistent with an age of about 100 yr and with the inferred mean shock velocity of 14,000 km/s deduced assuming a distance of 8.5 kpc. Most thermal emission comes from regions of lower X-ray but higher radio surface brightness. Deeper observations should allow more detailed spatial mapping of scandium, with significant implications for models of nucleosynthesis in Type Ia supernovae.

preprint2010arXiv

Steady X-Ray Synchrotron Emission in the Northeastern Limb of SN 1006

We investigate time variations and detailed spatial structures of X-ray synchrotron emission in the northeastern limb of SN 1006, using two Chandra observations taken in 2000 and 2008. We extract spectra from a number of small (about 10") regions. After taking account of proper motion and isolating the synchrotron from the thermal emission, we study time variations in the synchrotron emission in the small regions. We find that there are no regions showing strong flux variations. Our analysis shows an apparent flux decline in the overall synchrotron flux of about 4% at high energies, but we suspect that this is mostly a calibration effect, and that flux is actually constant to about 1%. This is much less than the variation found in other remnants where it was used to infer magnetic-field strengths up to 1 mG. We attribute the lack of variability to the smoothness of the synchrotron morphology, in contrast to the small-scale knots found to be variable in other remnants. The smoothness is to be expected for a Type Ia remnant encountering uniform material. Finally we find a spatial correlation between the flux and the cut-off frequency in synchrotron emission. The simplest interpretation is that the cut-off frequency depends on the magnetic-field strength. This would require that the maximum energy of accelerated electrons is not limited by synchrotron losses, but by some other effect. Alternatively, the rate of particle injection and acceleration may vary due to some effect not yet accounted for, such as a dependence on shock obliquity.

preprint2009arXiv

Deep Chandra Observations of the Crab-like Pulsar Wind Nebula G54.1+0.3 and Spitzer Spectroscopy of the Associated Infrared Shell

G54.1+0.3 is a young pulsar wind nebula (PWN), closely resembling the Crab, for which no thermal shell emission has been detected in X-rays. Recent Spitzer observations revealed an infrared (IR) shell containing a dozen point sources arranged in a ring-like structure, previously proposed to be young stellar objects. An extended knot of emission located in the NW part of the shell appears to be aligned with the pulsar's X-ray jet, suggesting a possible interaction with the shell material. Surprisingly, the IRS spectrum of the knot resembles the spectrum of freshly formed dust in Cas A, and is dominated by an unidentified dust emission feature at 21 microns. The spectra of the shell also contain various emission lines and show that some are significantly broadened, suggesting that they originate in rapidly expanding supernova (SN) ejecta. We present the first evidence that the PWN is driving shocks into expanding SN ejecta and we propose an alternative explanation for the origin of the IR emission in which the shell is composed entirely of SN ejecta. In this scenario, the freshly formed SN dust is being heated by early-type stars belonging to a cluster in which the SN exploded. Simple dust models show that this interpretation can give rise to the observed shell emission and the IR point sources.

preprint2007arXiv

Spitzer Space Telescope Observations of Kepler's Supernova Remnant: A Detailed Look at the Circumstellar Dust Component

We present 3.6 - 160 micron infrared images of Kepler's supernova remnant (SN1604) obtained with the IRAC and MIPS instruments on the Spitzer Space Telescope. We also present MIPS SED low resolution spectra in the 55 - 95 micron region. The observed emission in the MIPS 24 micron band shows the entire shell. Emission in the MIPS 70 micron and IRAC 8 micron bands is seen only from the brightest regions of 24 micron emission, which also correspond to the regions seen in optical Halpha images. Shorter wavelength IRAC images are increasingly dominated by stars, although faint filaments are discernible. The SED spectrum of shows a faint continuum dropping off to longer wavelengths and confirms that strong line emission does not dominate the mid-IR spectral region. The emission we see is due primarily to warm dust emission from dust heated by the primary blast wave; no excess infrared emission is observed in regions where supernova ejecta are seen in X-rays. We use models of the dust to interpret the observed 70/24 micron ratio and constrain the allowed range of temperatures and densities. We estimate the current mass of dust in the warm dust component to be 5.4 x 10^{-4} solar masses, and infer an original mass of about 3 x 10^{-3} solar masses before grain sputtering. The MIPS 160 micron band shows no emission belonging to the remnant. We place a conservative but temperature dependent upper limit on any cold dust component roughly a factor of 10 below the cold dust mass inferred from SCUBA observations. Finally, we comment on issues relevant to the possible precursor star and the supernova type.