Search NASA⌕ Search

Engineering topics

Williams, Brian J.

Publications and source records attributed to Williams, Brian J..

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multi-phenomenology Yield Characterization

This report serves as the first delivery of a four-year applied science effort to transform and advance the error bounds for the yield estimate of an explosion. Each year’s delivery will be in this form, culminating in the submission of this work for peer review to a scientific journal. Importantly, the yearly progress reports can then also be viewed as expanding drafts working towards a formal journal article submission. For the first tranche of funding, we collaborated with Air Force Technical Applications Center (AFTAC) scientists to identify unclassified real-world data that demonstrate and validate our advanced error propagation methods. Collaboration includes visits to AFTAC and telecons. For this development, we illustrate the fusion of seismic, acoustic, optical, and surface effect signatures from an explosion. The mathematics and code being adapted to this specific application (Williams et al., 2021) involves physics models of multiple sensor signatures. We have also identified related physics models and have integrated them into code. Current methods of underground explosion yield estimation for the Threshold Test Ban Treaty (TTBT) have served the US treaty monitoring mission well for decades. A research objective of the Defense Nuclear Nonproliferation Research and Development (DNN R&D) office of the National Nuclear Security Administration (NNSA) has always been to provide new technical capabilities for monitoring lower thresholds. The general model and error propagation code to be developed in this project is based on significant advances in error modeling and propagation needed to analyze data at lower detection thresholds. The second tranche of funding for this project began on May 1, 2022, and planned work for the second tranche includes: i) completing the integration of physical model code into the general error model framework; this code accommodates a wide range of linear/nonlinear source models, fixed/ random effects, and frequentist/Bayesian analyses (the purpose of which is not to dictate to users how to analyze data, but instead to allow users the maximum flexibility in their work); ii) illustrative application of code to identified data, and; iii) initial planning with AFTAC researchers on delivery of code to the Common Development Environment at AFTAC, and continued writing of the planned final journal article submission (year two of this progress report), with particular emphasis on descriptions of data identified for this effort.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Advanced Error Modeling for Prompt Yield Estimation [Slides]

Multi-Phenomenology Explosion Monitoring (MultiPEM) Toolbox provides a capability for post-detonation characterization focusing on yield estimation; Synthesize information from multiple sensor types; Leverage historical databases to inform statistical models; Allows for rapid or comprehensive analysis; Implemented in R code.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2018

Expansion and Age of the X-Ray Synchrotron-Dominated Supernova Remnant G330.2+1.0

We report new Chandra observations of one of the few Galactic supernova remnants whose X-ray spectrum is dominated by nonthermal synchrotron radiation, G330.2+1.0. We find that between 2006 and 2017, some parts of the shell have expanded by about 1%, giving a free-expansion (undecelerated) age of about 1000 yr, and implying shock velocities there of 9000 km/s for a distance of 5 kpc. Somewhat slower expansion is seen elsewhere around the remnant periphery, in particular in compact knots. Because some deceleration must have taken place, we infer that G330.2+1.0 is less than about 1000 yr old. Thus, G330.2+1.0 is one of only four Galactic core-collapse remnants of the last millennium. The large size, low brightness, and young age require a very low ambient density, suggesting expansion in a stellar-wind bubble. We suggest that in the east, where some thermal emission is seen and expansion velocities are much slower, the shock has reached the edge of the cavity. The high shock velocities can easily accelerate relativistic electrons to X-ray-emitting energies. A few small regions show highly significant brightness changes by 10%–20%, both brightening and fading, a phenomenon previously observed in only two supernova remnants, indicating strong and/or turbulent magnetic fields.

Borkowski, Kazimierz J.↗
NASA NTRS2018

The Expansion of the Young Supernova Remnant 0509-68.7 (N103B)

We present a second epoch of Chandra observations of the Type Ia Large Magellanic Cloud supernova remnant(SNR) 0509-68.7 (N103B) obtained in 2017. When combined with the earlier observations from 1999, we have a17.4 year baseline with which we can search for evidence of the remnant's expansion. Although the lack of strongpoint source detections makes absolute image alignment at the necessary accuracy impossible, we can measure thechange in the diameter and the area of the remnant, and find that it has expanded by an average velocity of 4170(2860, 5450) km s1. This supports the picture of this being a young remnant; this expansion velocity correspondsto an undecelerated age of 850 years, making the real age somewhat younger, consistent with results from lightecho studies. Previous infrared observations have revealed high densities in the western half of the remnant, likelyfrom circumstellar material, so it is probable that the real expansion velocity is lower on that side of the remnantand higher on the eastern side. A similar scenario is seen in Kepler's SNR. N103B joins the rare class ofMagellanic Cloud SNRs with measured proper motions.

Williams, Brian J.↗
NASA NTRS2018

A Deep X-Ray View of the Synchrotron-Dominated Supernova Remnant G330.2+1.0

We present moderately deep (125 ks) XMM-Newton observations of supernova remnant G330.2+1.0. This remnant is one of only a few known that fall into the "synchrotron-dominated" category, with the emission almost entirely dominated by a nonthermal continuum. Previous X-ray observations could only characterize the spectra of a few regions. Here, we examine the spectra from 14 regions surrounding the entire rim, finding that the spectral properties of the nonthermal emission do not vary significantly in any systematic way from one part of the forward shock to another, unlike several other remnants of this class. We confirm earlier findings that the power-law index, Γ, ranges from about 2.1-2.5, while the absorbing column density is generally between (2.0-2.6) × 10(exp 22) cm(exp -2). Fits with the srcut model find values of the roll-off frequency in the range of 10(exp 17.1)-10(exp 17.5) Hz, implying energies of accelerated electrons of approximately 100 TeV. These values imply a high shock velocity of approximately 4600 km s(exp -1), favoring a young age of the remnant. Diffuse emission from the interior is nonthermal in origin as well, and fits to these regions yield similar values to those along the rim, also implying a young age. Thermal emission is present in the east, and the spectrum is consistent with approximately 650 km s(exp -1) shock wave encountering interstellar or circumstellar material with a density of approximately 1 cm(exp -3).

Williams, Brian J.↗
NASA NTRS2017

The Three-Dimensional Expansion of the Ejecta from Tycho's Supernova Remnant

We present the first 3D measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12 yr baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 tufts of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line-of-sight velocity, we use two different methods: a nonequilibrium ionization model fit to the strong Si and S lines in the 1.22.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the redshift or blueshift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km/s, with a mean of 4430 km/s. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km/s. Some SN Ia explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and we discuss our findings in light of various explosion models, favoring those delayed-detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

ISM: individual objects (Tycho’s SNR) – ISM: s↗
NASA NTRS2017

The Refined Shock Velocity of the X-Ray Filaments in the RCW 86 Northeast Rim

A precise measurement of shock velocities is crucial for constraining the mechanism and efficiency of cosmic-ray (CR) acceleration at supernova remnant (SNR) shock fronts. The northeastern rim of the SNR RCW 86 is thought to be a particularly efficient CR acceleration site, owing to the recent result in which an extremely high shock velocity of 6000 km s1 was claimed. Here, we revisit the same SNR rim with the Chandra X-ray Observatory, 11 years after the first observation. This longer baseline than previously available allows us to determine a more accurate proper motion of the nonthermal X-ray filament, revealing a much lower velocity of 3000 340 km s1 (and even slower at a brighter region). Although the value has dropped to one-half of that from the previous X-ray measurement, it is still higher than the mean velocity of the H filaments in this region (1200 km s1). This discrepancy implies that the filaments bright in nonthermal X-rays and H emission trace different velocity components, and thus a CR pressure constrained by combining the X-ray kinematics and the H spectroscopy can easily be overestimated. We also measure the proper motion of the thermal X-ray filament immediately to the south of the nonthermal one. The inferred velocity (720 360 km s1) is significantly lower than that of the nonthermal filament, suggesting the presence of denser ambient material, possibly a wall formed by a wind from the progenitor, which has drastically slowed down the shock.

acceleration of particles - ISM: individuals objec↗
NASA NTRS2016

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 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 years in the X-rays and 30 years in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. Consistent with earlier measurements, we find 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 approximately 23″ toward the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

dust↗
NASA NTRS2015

SOFIA Observations of SN 2010jl: Another Non-Detection of the 9.7 Micrometer Silicate Dust Feature

We present photometric observations from the Stratospheric Observatory for Infrared Astronomy (SOFIA) at 11.1 micrometers of the Type IIn supernova (SN IIn) 2010jl. The SN is undetected by SOFIA, but the upper limits obtained, combined with new and archival detections from Spitzer at 3.6 and 4.5 micrometers, allow us to characterize the composition of the dust present. Dust in other SN IIn has been shown in previous works to reside in a circumstellar shell of material ejected by the progenitor system in the few millennia prior to explosion. Our model fits show that the dust in the system shows no evidence for the strong, ubiquitous 9.7 micrometer feature from silicate dust, suggesting the presence of carbonaceous grains. The observations are best fit with 0.01−0.05 solar mass of carbonaceous dust radiating at a temperature of approximately 550−620 degrees Kelvin. The dust composition may reveal clues concerning the nature of the progenitor system, which remains ambiguous for this subclass. Most of the single star progenitor systems proposed for SNe IIn, such as luminous blue variables, red supergiants, yellow hypergiants, and B[e] stars, all clearly show silicate dust in their pre-SN outflows. However, this post-SN result is consistent with the small sample of SNe IIn with mid-infrared observations, none of which show signs of emission from silicate dust in their infrared spectra.

SOFIA Observations of SN 2010jl↗
NASA NTRS2015

A Chandrasekhar Mass Progenitor for the Type Ia Supernova Remnant 3C 397 from the Enhanced Abundances of Nickel and Manganese

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

G41.1-0.3); supernova remnants; nuclear reactions;↗
NASA NTRS2014

Magnetic Field Amplification in the Thin X-Ray Rims of SN 1006

Several young supernova remnants, including SN 1006, emit synchrotron X-rays in narrow filaments, hereafter thin rims, along their periphery. The widths of these rims imply 50–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 ∝ 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 FWHM in terms of energy as FWHM ∝ EmE γ. Filament widths in SN 1006 decrease with energy; mE ∼ −0.3 to −0.8, implying magnetic field amplification by factors of 10–50, above the factor of four 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 (μ ≥ 1).

MAGNETIC FIELD AMPLIFICATION IN THE THIN X-RAY RIM↗
NASA NTRS2014

A High-Resolution X-Ray and Optical Study of SN1006: Asymmetric Expansion and Small-Scale Structure in a Type Ia Supernova Remnant

We introduce a deep (670 ks) X-ray survey of the entire SN 1006 remnant from the Chandra X-Ray Observatory, together with a deep Hα image of SN 1006 from the 4 m Blanco telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of the X-ray proper motions around the entire periphery, carried out over a 9 yr baseline. We find that the expansion velocity varies significantly with azimuth. The highest velocity of approx.7400 km/s (almost 2.5 times that in the northwest (NW)) is found along the southeast (SE) periphery, where both the kinematics and the spectra indicate that most of the X-ray emission stems from ejecta that have been decelerated little, if at all. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Neon emission is strongest in a sharp filament just behind the primary shock along the NWrim, where the pre-shock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in Hα, features that we interpret as ejecta knots that have reached relatively dense regions of the surrounding interstellar medium, but that appear in the interior in projection. Recent three-dimensional hydrodynamic models for Type Ia supernovae display small-scale features that strongly resemble the ones seen in X-rays in SN 1006; an origin in the explosion itself or from subsequent hydrodynamic instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front, as expected from diffusive shock acceleration theory, to numerous regions along both the northeast and southwest rims of the shell. Our data require that a precursor be thinner than about 3, and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of seven or more in a narrow precursor region, promoting diffusive particle acceleration.

galaxies: active↗
NASA NTRS2014

Grain Destruction in a Supernova Remnant Shock Wave

Dust grains are sputtered away in the hot gas behind shock fronts in supernova remnants, gradually enriching the gas phase with refractory elements. We have measured emission in C IV (lambda)1550 from C atoms sputtered from dust in the gas behind a non-radiative shock wave in the northern Cygnus Loop. Overall, the intensity observed behind the shock agrees approximately with predictions from model calculations that match the Spitzer 24 micron and the X-ray intensity profiles. Thus these observations confirm the overall picture of dust destruction in SNR shocks and the sputtering rates used in models. However, there is a discrepancy in that the CIV intensity 10'' behind the shock is too high compared to the intensities at the shock and 25'' behind it. Variations in the density, hydrogen neutral fraction and the dust properties over parsec scales in the pre- shock medium limit our ability to test dust destruction models in detail.

dust↗
NASA NTRS2012

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

We report the results of an X-ray proper motion measurement for the NW rim of SN 1006, 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 SN 1006 with signi cant optical emission: a thin, Balmer-dominated lament. For most of the NW rim, the proper motion is approximately equal to 0.30"yr(exp -1), essentially the same as has been measured from the H alpha lament. Isolated regions of the NW limb are dominated by nonthermal emission, and here the proper motion is much higher, 0:49"yr(exp -1), 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 SN 1006, the proper motions imply shock velocities of approximately 3000 kms(exp -1) and approximately 5000 kms(exp -1) in the thermal and nonthermal regions, respectively. A lower velocity behind the H alpha filament is consistent with the picture that SN 1006 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 o set 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.

Katsuda, Satoru↗
NASA NTRS2012

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 micron IR spectra of the remnant, obtained with the Spitzer Space Telescope, dominated by emission from warm dust. Broad spectral features at 10 and 18 micron, 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) propagating into the CSM can reproduce the majority of the emission associated with non-radiative filaments, where dust temperatures are approx 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) into moderate density material (n(sub o) approx 50-100 / cubic cm) 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.

Williams, Brian J.↗
NASA NTRS2012

The First Reported Infrared Emission from the SN1006 Remnant

We report results of infrared imaging and spectroscopic observations of the SN 1006 remnant, carried out with the Spitzer Space Telescope. The 24 m 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 m 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 lowera relic effect from an earlier epoch when the shock was encountering a lower densitybut 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 SNIa remnants, SN1006 shows no evidence for dust grain formation in the supernova ejecta.

Winkler, P. Frank↗