Search NASA⌕ Search

SEARCH · Search NASA

Results for “NUCLEAR EXPLOSION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Influence of background sources and topographic resolution in the Weather Research and Forecasting Model on xenon plume characteristics at monitoring stations

For many atmospheric monitoring applications, networks of measurement sites—such as the radionuclide stations of the International Monitoring System—can be sparse. With measurement locations potentially hundreds to thousands of kilometers from a release it is important to quantify the effects of physical processes on transport and dispersion of plumes between source and measurement locations. This study addresses the effects of background sources and topography resolution near the release location of radionuclides. We use the Weather Research and Forecasting (WRF) model with inline chemistry to investigate (1) how an additional, time-varying source of 133 Xe, such as an operational medical isotope production facility, contributes to activity concentration measurements at monitoring sites, and (2) how complex topography influences on atmospheric conditions near emission sources impact plume concentrations at varying distances from the source. Two 133 Xe emission sources, including (1) a high flux rate of short duration representing an explosive event, and (2) a variable and continuous background source, are simulated. The continuous background source contributes significantly to total 133 Xe concentrations at several monitoring stations. Further, a WRF simulation at 9 km horizontal resolution is compared with a nested grid simulation, where the innermost domain has a resolution of 1 km. Increased topographic resolution leads to an improved representation of plume responses to local winds, with topographic influences greatest at locations closest to the sources. Differences between the two domain resolutions decrease at greater distances from the sources, as plumes have time to spread and mix and are influenced by synoptic scale circulation patterns that are represented similarly in both simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Analysis of measurements from an array of radioxenon samplers near to Hartlepool Nuclear Power Station

As part of a scientific research and development project, the radionuclide fingerprint of an operating advanced gas-cooled reactor 25 (AGR) has been studied across several facets (Goodwin et al., 2024). One part of this project was to deploy an array of 26 radioxenon samplers to the region for a period of around 1 year, to measure any radioxenon emissions from the reactors of the 27 Hartlepool nuclear power station at a range of tens of kilometres away. The array of 3 sensors was operational for around 12 28 months from March 2022 and detected many occurrences of isotopes of radioxenon. Here we provide a detailed analysis and 29 interpretation of the data and where possible, attribution of detections to a source or region. Whilst a key part of this work is to 30 measure any emissions of radioxenon from Hartlepool, this work presents one of the most comprehensive efforts to determine the 31 source of a great number of (mostly) 133Xe detections. A combination of different types of atmospheric dispersion modelling 32 techniques, including the use of stack monitoring data from nearby civil radioxenon-emitting nuclear facilities has enabled the 33 majority of detections to be attributed to one or more possible sources. Whilst emissions from Hartlepool have been detected on 34 the systems, the majority of detections are associated with a medical isotope production facility in Fleurus, Belgium (IRE).

Advanced gas-cooled reactor↗
DOE OSTI2025

Development and validation of a software for simulating γ-γ coincidence emission and detection probabilities

Gamma-gamma coincidence spectrometers have the potential to significantly enhance detection sensitivity for ultra-trace radionuclide measurements. The implementation of these spectrometers, however, is limited by the complexity of acquisition hardware, data processing and quantification. This work reports development of a novel radionuclide quantification software for γ-γ coincidence measurements. For any radionuclide, the software parses the Evaluated Nuclear Structure Data File (ENSDF) database, recursively simulating all possible γ-γ coincidence signatures and their respective emission and detection probabilities. Implemented using Python programming language, the software employs several strategies to boost overall computational performance. Since coincidence-based spectrometers are of notable interest in monitoring compliance for the Comprehensive Nuclear-Test-Ban Treaty (CTBT), the software’s execution was tested for 84 CTBT-relevant radionuclides. To date, the software has been experimentally validated for 15 radionuclides using the Advanced Radionuclide Gamma spectrOmeter (ARGO) at Pacific Northwest National Laboratory, USA (PNNL). Notably, the software can be operated in convergence mode, whereby coincidence detection efficiency’s convergence behavior can help avoid unreliable radionuclide activity estimates. With growing number of coincidence spectrometers worldwide, this paper aims to assist the radiation metrology community in developing similar software for their system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Novel application of neutrinos to evaluate U.S. nuclear weapons performance

There is a growing realization that neutrinos can be used as a diagnostic tool to better understand the inner workings of a nuclear weapon. Robust estimates demonstrate that an Inverse Beta Decay (IBD) neutrino scintillation detector built at the Nevada Test Site with a 1000-ton active target mass at a standoff distance of 500 m would detect thousands of antineutrino events per nuclear test. This would provide less than 4% statistical error on the measured antineutrino rate and 5% error on antineutrino energy. Extrapolating this to an error on the test device explosive yield requires knowledge from evaluated nuclear databases, non-equilibrium fission rates, and assumptions on internal neutron fluxes. Initial calculations demonstrate that the total number of neutrinos emitted per fission in the first 10 3 s after a short pulse of 239 Pu fission is about a factor of two less than that from Pu fissioning under steady state conditions. Furthermore, there are significant energy spectral differences as a function of time after the pulse that must be considered. These and other model dependencies will be discussed in the paper. In the absence of nuclear weapons testing, many of the technical and theoretical challenges of a full nuclear test could be mitigated with a low cost smaller scale 20 ton fiducial mass IBD demonstration detector placed near a pulsed reactor. Potential reactors include the Texas A&M University TRIGA 1 GW–10 ms pulsed facility or the Sandia Annular Core Research Reactor. The short duty cycle and repeatability of pulses would provide critical real environment testing and measurements, which would be valuable for planning a possible real test shot in the future. Furthermore, the antineutrino rate as a function of time data would provide unique constraints on fission databases and model assumptions. Finally, there are impactful science drivers such as sensitive searches for ∼1 eV 2 sterile neutrinos and ∼MeV scale axions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗
DOE OSTI2026

An Exploration of the Equation of State Dependence of Core-collapse Supernova Explosion Outcomes and Signatures

We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.

Rusakov, Aleksandr [Princeton University, NJ (Unit↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2026

Single-Nucleon Transfer on Unstable 59 Cu Probes the NiCu Cycle in Astrophysical X-Ray Bursts

Recent models of the rapid proton (𝑟⁢𝑝) capture process indicate that a competition between the 59 Cu⁢(𝑝,𝛾) ⁢60 Zn and 59 Cu⁢(𝑝,𝛼) ⁢56 Ni reactions may result in the formation of a nickel-copper (NiCu) cycle that traps the flux of material between 56 Ni and 60 Zn . Here, in this work, we report the identification of 15 proton-unbound levels in 60 Zn , populated via 59 Cu⁢(𝑑,𝑛) transfer, which govern the rate of the 59 Cu⁢(𝑝,𝛾)⁢ 60 Zn reaction in XRBs. Precise excitation energies for levels in 60 Zn were obtained from observed 𝛾 decays, and spectroscopic factors were determined from angle-integrated cross sections. Incorporating these results into stellar-model calculations, we find that with experimentally constrained uncertainties a NiCu cycle in XRBs is indeed possible, though we limit its branching strength to less than 38%. While modest, such a branching has significant impact on the light curve, motivating further studies of the relevant rates. Our calculations also indicate that a significant NiCu cycle leads to an increase in the amount of odd-𝐴 nuclei in the burst ashes, which may affect Urca cooling processes in neutron star crusts.

direct reactions↗
DOE OSTI2025

Measurements of Gamow-Teller transitions from 59 Co via the 59 Co ⁢(𝑡, 3 He +𝛾) charge-exchange reaction and its application to the stellar electron-capture rates

Electron-capture reactions on iron-group nuclei play a crucial role in the late stages of massive star evolution. Since stellar evolution simulations depend on accurate electron-capture rates—which are highly sensitive to the detailed Gamow-Teller (GT) strength distributions—reliable theoretical models are essential. However, experimental data on GT strength distributions are scarce. High-resolution measurements are therefore vital for benchmarking and improving these theoretical calculations. To provide high-resolution data on Gamow-Teller strength distributions of iron-group nuclei and to compare these results with theoretical calculations within this mass region. Differential cross sections for the 59 Co ⁢(𝑡, 3 He)⁢ 59 Fe charge-exchange reaction at 115 MeV/u were measured using the S800 spectrometer. Furthermore, to resolve individual levels that are not distinguishable in the S800 particle singles data, coincident 𝛾 rays from the 59 Fe residual nucleus were detected by using the Gamma-Ray Energy Tracking In-beam Nuclear Array 𝛾-ray tracking array. Here, the Gamow-Teller transition strength distribution from the ground state of 59 Co to 59 Fe was extracted up to an excitation energy of 10 MeV. Additionally, transition strengths for several low-lying states were determined from coincident 𝛾-ray measurements. Electron-capture rates calculated using the present data indicate that these low-lying states contribute significantly to the overall rates in relevant stellar environments. The experimental results show reasonable agreement with theoretical predictions based on both shell-model and projected shell-model calculations. High-resolution data on Gamow-Teller strength distributions—particularly for individual low-lying states—are essential for accurately determining electron-capture rates in iron-group nuclei. Coincident 𝛾-ray measurements provide a powerful tool for obtaining such detailed information. While the present work demonstrates that shell-model calculations successfully reproduce the experimental results, such comparisons are scarce and more experimental data are desirable.

59 ≤ A ≤ 89↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2025

PySolate : A Python‐Based Thresholding Tool to Denoise or Designal Seismic Waveforms Based on the Continuous Wavelet Transform

PySolate is a Python‐based toolset that implements the continuous wavelet transform and nonlinear thresholding operations to denoise or designal seismic data, following Langston and Mousavi (2019). This filtering approach can remove microseismic noise to isolate intermediate‐period seismic signals that are key to enabling full‐waveform modeling and analysis of smaller‐magnitude regional events. This approach is best for the application to signals with frequency or time separation of signal and noise, in contrast to Fourier analysis, which is effective when signal and noise are separated in frequency. We demonstrate the Python toolset using the six announced Democratic People’s Republic of Korea declared nuclear tests, showing the effectiveness of isolating the seismic signal compared to standard bandpass filtering. In conclusion, we also demonstrate the ease of using the toolset with any Python processing tools.

Asia↗
DOE OSTI2025

Shear-Wave Radiation Patterns from Explosive and Earthquake Sources in Scattering, Heterogeneous Media

Distinguishing whether a seismic event is an earthquake or an explosion is a core problem in explosion monitoring. For simplistic models, earthquakes produce a predictable S-wavefield based on their radiation pattern, while an explosion produces no S-wavefield. However, observations from nuclear tests show that explosions can produce significant S-wave energy on both horizontal components. We perform numerical experiments using SPECFEM to constrain when an S-wavefield generated from an explosion might differ from one produced by an earthquake. We generate the S-wavefield for our explosions by placing the source location within a small region consisting of large velocity heterogeneities. For the 2D cases investigated here, a crucial condition is that the size of the region of heterogeneities must be much smaller (approximately 1/12th) than the minimum wavelength of the wavefield for our simulation setups. We investigate both isotropic heterogeneities and anisotropic heterogeneities. Initial results demonstrate that the heterogeneous region lowers the P/S amplitude ratios of an explosion. The presence of strong anisotropic heterogeneities near the source can also produce an S-wavefield with similar amplitudes to the P-wavefield. Lastly, we show that S-wave polarization angles vary as a function of azimuth and distance for explosive sources in the presence of small-scale heterogeneities.

58 GEOSCIENCES↗
DOE OSTI2025

A case study in contrastive learning information combination: Application to technical forensics of additive manufacturing filament source identification

Combination of information from disparate data sources into a single decision is a core challenge in many fields, including the field of technical forensics. Technical forensics (TF) utilizes technical characterization of questioned samples to determine properties of that sample; these properties are then used to infer information of forensic interest, such as provenance, age, or attribution. TF is utilized in traditional forensic applications, such as the attribution of material fragments from an explosive, and in nuclear forensic applications, such as the attribution of actinides which have been interdicted out of regulatory control. The challenge of combining information from disparate sources, described alternately by many terms including “Data Fusion” and “Data Integration”, is exacerbated in the technical forensics domain due to at least two factors: the challenge of interpreting each information source singularly, and the relatively small data set sizes available. Extensive literature exists attempting to combine technical forensics information sources, both in manual and automated processes. These attempts are often bespoke to the specific information sources (such as the bi-, tri-, or quad-isotope chart (Moody, Grant, and Hutcheon 2005)), with some emerging examples of simple early- and late- fusion (, respectively). Simultaneous to the information combination efforts described in the previous paragraph, the field of natural language processing attempted (and largely succeeded) in combining information from multiple non-technical information sources. The ecosystem of “multi-modal” language models, which can take text and images as input, and generate text and images as output, became large and diverse by 2025 (Khan et al. 2025). In a generalized sense, many of these methods are trained by learning neural networks which can convert raw text or images into a vector of numbers describing the text or image, hereafter called “embeddings” and the neural networks performing the conversion are called “embedders”. By using a separate embedder for text and images, finding coincident text and images (such as images with their captions), and optimizing the parameters of the embedders such that the embeddings for the text and the image are similar, the field has found a bridge between text and images (Girdhar et al. 2023). It is the contention of the authors of this report that this insight is not limited to text and images but instead can be extended to any modality which can be found coincidently. The subject of the rest of this report is the application of this method to example multi-modal technical forensic data. Some details about the data used in this report are not appropriate for this report, and are included in a companion report (PNNL-38669).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

First Proton-Induced Cross Sections on a Stored Rare Ion Beam: Measurement of Te 118 ( p , γ ) for Explosive Nucleosynthesis

We present the first nuclear cross-section measurements of (p,$\gamma$) and (p,n) reactions on 118 Te at energies relevant for the $\gamma$-process nucleosynthesis. Absolute cross-section values for center-of-mass energies of 6, 7 and 10 MeV are provided, together with a theoretical extrapolation to the Gamow window. This experiment marks the first time that direct proton-induced reactions have been measured on a radioactive ion beam at the Experimental Storage Ring (ESR) at GSI, Darmstadt. This paves the way for a large variety of measurements, delivering new constraints for explosive nucleosynthesis and for physics beyond nuclear stability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A flow-through gamma detection system for sub-surface radioxenon tracer measurements in LYNM PE1

A continuous flow-through, real-time radioxenon sensor array has been designed to monitor either 127 Xe or 133 Xe tracer gas as part of a series of underground chemical explosions in the Low-Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1). The sensors are pairs of shielded thallium-doped sodium iodide gamma-ray detectors that surround a gas measurement chamber. Each measurement chamber is connected to a sampling location in the geology surrounding the chemical explosion or in the access tunnel. The first experiment in the series, PE1 A, occurred in 2023 and the system quantified the 127 Xe tracer gas in the gas sampling lines for four weeks following the experiment execution. The array generated a set of radioxenon concentration time series with a 15-minute resolution that will be used to refine subsurface gas transport models.

Gamma-ray coincidence detection↗
DOE OSTI2025

Quantum physics of stars

Stars are slowly developing objects; the lifetimes of the different burning phases are determined by the strength of nuclear reactions, which in turn are defined by the quantum structure of the associated nuclei at the threshold and the respective reaction mechanisms. Stars, from the nuclear physics perspective, are cold environments where only a few of the key nuclear reactions have been measured at the actual stellar plasma temperatures. This is also the case for more dynamic astrophysical phenomena from the big bang to stellar explosions. Most of the nuclear reaction rates are therefore based on theoretical extrapolations. A number of discrepancies between these predictions and the associated stellar signatures have been observed, and many may be due to low-energy or near-threshold quantum effects. These effects need to be understood in order to reliably model nuclear reaction processes, not only for stars but also for low-temperature plasma environments such as controlled magnetic or inertial confinement fusion systems, which operate in similar temperature regimes. This review summarizes the various theoretical techniques presently used for deriving reaction rates and discusses possible quantum effects that may impact the reaction cross section near the reaction threshold. These resemble enhanced single-particle and cluster structures near threshold and associated interference effects. New experimental techniques such as deep-underground accelerators or the study of transfer reactions to mimic the quantum-mechanical transition strength, the so-called Trojan horse method, provide ways to directly or indirectly probe the reaction features that determine the reaction rates at stellar energies. Furthermore, this is demonstrated on a number of key nuclear reactions for different nucleosynthesis environments. Finally, current inconsistencies between experimental predictions and observations are discussed.

Models & methods for nuclear reactions↗