Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model interpretation”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Geothermal Exploration of Mount Baker Hot Springs Through Ground-Based Magnetic and Gravity Surveys

The Washington State Division of Geology and Earth Resources (WDGER), together with the USGS, Altarock Energy Inc., Colorado College, Temple University, Innovate Geothermal Ltd., University of British Columbia, and Western Washington University(WWU), conducted phase two of a geothermal play fairway exploration at three plays (Mount Baker, Mount St. Helens seismic zone, and Wind River valley) along the Washington Cascade Range. The Mount Baker (MB) play is an ideal location for further investigation as a potential geothermal resource due to the presence of thermal features, young volcanic centers, existing geothermal leases, accessibility, and its proximity to existing transmission lines at the Baker Lake Dam. The goal of phase II exploration at all three plays is to collect geologic and geophysical data to help constrain spatial changes in subsurface lithology and structures that may be associated with the geothermal systems. Furthermore, this paper discusses the geothermal potential of the southeast flank of Mount Baker, WA by interpreting newly collected ground-based gravity and magnetic data. The Mount Baker play was surveyed from July-September, 2016 by the USGS, WDGER, and WWU. Results of a detailed magnetic and gravity surveys of the region surrounding the Baker hot springs, the primary surface expression of the hydrothermal system, are presented here. Gravity data were collected at 495 stations in the study area amounting to an area of roughly 150 square km. Approximately 93 km of line magnetometer data were collected within the study area, over a dozen magnetic susceptibility measurements were taken at 50 locations and about 50 hand samples were collected for density measurements. Preliminary data and modeling suggest pronounced magnetic anomalies trending NE/SW, generally in line with previously mapped lineaments identified in LiDAR. We present a series of preliminary two-dimensional forward models of the gravity and magnetic data as well as a 3D inversion model of magnetic susceptibility and interpret these models in terms of the structural controls influencing the Baker hot springs and hydrothermal system.

Geophysics↗
DOE OSTI2018

Internal Variability and Disequilibrium Confound Estimates of Climate Sensitivity From Observations

An emerging literature suggests that estimates of equilibrium climate sensitivity (ECS) derived from recent observations and energy balance models are biased low because models project more positive climate feedback in the far future. Here we use simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) to show that across models, ECS inferred from the recent historical period (1979–2005) is indeed almost uniformly lower than that inferred from simulations subject to abrupt increases in CO2 radiative forcing. However, ECS inferred from simulations in which sea surface temperatures are prescribed according to observations is lower still. ECS inferred from simulations with prescribed sea surface temperatures is strongly linked to changes to tropical marine low clouds. However, feedbacks from these clouds are a weak constraint on long-term model ECS. One interpretation is that observations of recent climate changes constitute a poor direct proxy for long-term sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Search for π 0 decays to invisible particles

The NA62 experiment at the CERN SPS reports a study of a sample of 4 × 10 9 tagged π0 mesons from K + → π + π 0 (γ), searching for the decay of the π 0 to invisible particles. No signal is observed in excess of the expected background fluctuations. An upper limit of 4.4 × 10 –9 is set on the branching ratio at 90% confidence level, improving on previous results by a factor of 60. This result can also be interpreted as a model- independent upper limit on the branching ratio for the decay K + → π + X, where X is a particle escaping detection with mass in the range 0.110–0.155 GeV/c 2 and rest lifetime greater than 100 ps. Model-dependent upper limits are obtained assuming X to be an axion-like particle with dominant fermion couplings or a dark scalar mixing with the Standard Model Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Combined search for electroweak production of charginos and neutralinos in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A statistical combination of several searches for the electroweak production of charginos and neutralinos is presented. All searches use proton-proton collision data at $ \sqrt{s}=13 $ TeV, recorded with the CMS detector at the LHC in 2016 and corresponding to an integrated luminosity of 35.9 fb$^{−1}$. In addition to the combination of previous searches, a targeted analysis requiring three or more charged leptons (electrons or muons) is presented, focusing on the challenging scenario in which the difference in mass between the two least massive neutralinos is approximately equal to the mass of the Z boson. The results are interpreted in simplified models of chargino-neutralino or neutralino pair production. For chargino-neutralino production, in the case when the lightest neutralino is massless, the combination yields an observed (expected) limit at the 95% confidence level on the chargino mass of up to 650 (570) GeV, improving upon the individual analysis limits by up to 40 GeV. If the mass difference between the two least massive neutralinos is approximately equal to the mass of the Z boson in the chargino-neutralino model, the targeted search requiring three or more leptons obtains observed and expected exclusion limits of around 225 GeV on the second neutralino mass and 125 GeV on the lightest neutralino mass, improving the observed limit by about 60 GeV in both masses compared to the previous CMS result. In the neutralino pair production model, the combined observed (expected) exclusion limit on the neutralino mass extends up to 650–750 (550–750) GeV, depending on the branching fraction assumed. This extends the observed exclusion achieved in the individual analyses by up to 200 GeV. The combined result additionally excludes some intermediate gaps in the mass coverage of the individual analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Search for a new resonance decaying into two spin-0 bosons in a final state with two photons and two bottom quarks in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for a new boson X is presented using CERN LHC proton-proton collision data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV in 2016–2018, and corresponding to an integrated luminosity of 138 fb -1 . The resonance X decays into either a pair of Higgs bosons HH of mass 125 GeV or an H and a new spin-0 boson Y. One H subsequently decays to a pair of photons, and the second H or Y, to a pair of bottom quarks. The explored mass ranges of X are 260–1000 GeV and 300–1000 GeV, for decays to HH and to HY, respectively, with the Y mass range being 90–800 GeV. For a spin-0 X hypothesis, the 95% confidence level upper limit on the product of its production cross section and decay branching fraction is observed to be within 0.90–0.04 fb, depending on the masses of X and Y. The largest deviation from the background-only hypothesis with a local (global) significance of 3.8 (below 2.8) standard deviations is observed for X and Y masses of 650 and 90 GeV, respectively. The limits are interpreted using several models of new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Novel imaging techniques of rectal cancer: what do radiomics and radiogenomics have to offer? A literature review

Introduction: As computational capabilities have advanced, radiologists and their collaborators have looked for novel ways to analyze diagnostic images. This has resulted in the development of radiomics and radiogenomics as new fields in medical imaging. Radiomics and radiogenomics may change the practice of medicine, particularly for patients with colorectal cancer. Radiomics corresponds to the extraction and analysis of numerous quantitative imaging features from conventional imaging modalities in correlation with several endpoints, including the prediction of pathology, genomics, therapeutic response, and clinical outcome. In radiogenomics, qualitative and/or quantitative imaging features are extracted and correlated with genetic profiles of the imaged tissue. Thus far, several studies have evaluated the use of radiomics and radiogenomics in patients with colorectal cancer; however, there are challenges to be overcome before its routine implementation including challenges related to sample size, model design and interpretability, and the lack of robust multicenter validation set. Material and methods: In this article, we will review the concepts of radiomics and radiogenomics and their potential applications in rectal cancer. Conclusion: Radiologists should be aware of the basic concepts, benefits, pitfalls, and limitations of new radiomic and radiogenomics techniques to achieve a balanced interpretation of the results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Atmospheric humidity deficits tell us how soil moisture deficits down-regulate ecosystem evaporation

Knowing how actual evaporation is down-regulated from potential evaporation during periods with soil moisture deficits is one of the greatest challenges towards computing evaporation everywhere on a regular basis. We propose the hypothesis that vegetated landscapes transmit information on soil moisture deficits through its down-regulation of evaporation, which in turn affects the humidification and growth of the planetary boundary layer. To test this hypothesis, here, we examined how the evaporative fraction, defined as the ratio between actual and potential evaporation, corresponds with the soil moisture stress index, defined as the relative humidity with a power law dependence on vapor pressure deficit times a parameter, β. We tested the parameterized soil moisture stress index with direct eddy covariance measurements of actual evaporation and computations of potential evaporation based on meteorological conditions, averaged on monthly time steps. The analysis was conducted on a dataset obtained from 144 FLUXNET sites. These sites spanned much of the world's climates and biomes and contained over 5900 months of observations. Observations of the evaporative fraction fit the model of the soil moisture stress index best for semi-arid and arid ecosystems, and least for the humid and wet ecosystems. Consequently, a significant relationship between measurements of the evaporative fraction and soil moisture stress index held for about one-half of the population of sites. Under this condition, the median value of the parameter, β, was 1.41. To investigate the mechanism of this empirical soil moisture stress index, we diagnosed it with a coupled evaporation-planetary boundary layer model. The soil moisture stress index is strongly related to surface resistance, as defined by inverting Penman-Monteith equation. Consequently, this index provides an independent estimate of surface resistance based on easy to measure mean monthly weather conditions like relative humidity and temperature. Thereby, this soil moisture stress index has potential to be applied to weather, climate and biogeochemical models and the interpretation of satellite derived evaporation products, like the one provided by the ECOSTRESS mission.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Establishing nationwide power system vulnerability index across US counties using interpretable machine learning

Power outages have become increasingly frequent, intense, and prolonged in the US due to climate change, aging electrical grids, and rising energy demand. However, largely due to the absence of granular spatiotemporal outage data, we lack data-driven evidence and analytics-based metrics to quantify power system vulnerability. This limitation has hindered the ability to effectively evaluate and address vulnerability to power outages in US communities. Here, in this work, we collected ∼179 million power outage records at 15-min intervals across 3022 US contiguous counties (96.15 % of the area) from 2014 to 2023. We developed a power system vulnerability assessment framework based on three dimensions (intensity, frequency, and duration) and applied interpretable machine learning models (XGBoost and SHAP) to compute Power System Vulnerability Index (PSVI) at the county level. Our analysis reveals a consistent increase in power system vulnerability across the US counties over the past decade. We identified 318 counties across 45 states as hotspots for high power system vulnerability, particularly in the West Coast (California and Washington), the East Coast (Florida and the Northeast area), the Great Lakes megalopolis (Chicago-Detroit metropolitan areas), and the Gulf of Mexico (Texas). Our heterogeneity analysis indicates that urban counties and those located along regional transmission boundaries tend to exhibit significantly higher vulnerability. Our results highlight the significance of the proposed PSVI for evaluating the vulnerability of communities to power outages. The findings underscore the widespread and pervasive impact of power outages across the country and offer crucial insights to support infrastructure operators, policymakers, and emergency managers in formulating policies and programs aimed at enhancing the resilience of the US power infrastructure.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Structural characterization of high-protein system through ultra-small and small-angle X-ray scattering

High-protein systems exhibit a hierarchical structure consisting of interconnected multiscale length assemblies. Recent investigations are directed toward identifying these structural units under various processing and environmental conditions to establish structure–function correlations. Ultra-small and small-angle X-ray scattering (USAXS/SAXS) has become a crucial tool for characterizing the structures of proteins and their clusters/aggregates, ranging from nanometers to micrometers, with minimal disruption to their original state. Here, this review first describes the facilities and principles of X-ray scattering, followed by discussions on the analysis of scattering data, including the interpretation of fitting models. It then delves into the main applications of USAXS/SAXS in plant and dairy protein-rich systems. Future strategies to enhance the utilization of scattering techniques for elucidating the structure of high-protein systems are also included.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

A data-driven CO 2 leakage detection using seismic data and spatial–temporal densely connected convolutional neural networks

In carbon capture and sequestration (also known as carbon capture and storage, or CCS), developing effective monitoring methods is needed to detect and respond to CO 2 leakage. CO 2 leakage detection methods rely on geophysical observations and monitoring sensor network. However, traditional methods usually require the development of site-specific physical models and expert interpretation, and the effectiveness of these methods can be limited to different application locations, operational scenarios, and conditions. Here, we developed a novel data-driven leakage detection method based on densely connected convolutional neural networks. Our method is an end-to-end detection approach, that differs from conventional leakage monitoring methods by directly learning a mapping relationship between seismic data and the CO 2 leakage mass. To account for the spatial and temporal characteristics of seismic data, our novel networks architecture combines 1D and 2D convolutional neural networks. To overcome the computational expense of solving optimization problems, we apply a densely-connecting strategy in our network architecture that reduces the number of network parameters. Based on the features generated by our convolutional neural networks, we further incorporate a long short-term memory network to utilize time-sequential information, which further improves the detection accuracy. Lastly, we employ our detection method to synthetic seismic datasets generated based on flow simulations of a hypothetical CO 2 storage scenario with injection into a partially compartmentalized sandstone storage reservoir. To evaluate method performance, we conducted multiple experiments including a random leakage test, a sequential test, and a robustness test. Numerical results show that our CO 2 leakage detection method successfully detects the leakage and accurately predicts the leakage mass, suggesting that it has the potential for application in monitoring of real CO 2 storage sites.

58 GEOSCIENCES↗
DOE OSTI2022

Flow in additively manufactured super-rough channels

Metal additive manufacturing has enabled geometrically complex internal cooling channels for turbine and heat exchanger applications, but the process gives rise to large-scale roughness whose size is comparable to the channel height (which is 500 $\mathrm {\mu }$ m). These super-rough channels pose previously unseen challenges for experimental measurements, data interpretation and roughness modelling. First, it is not clear if measurements at a particular streamwise and spanwise location still provide accurate representation of the mean (time- and plane-averaged) flow. Second, we do not know if the logarithmic layer survives. Third, it is unknown how well previously developed rough-wall models work for these large-scale roughnesses. To answer the above practical questions, we conduct direct numerical simulations of flow in additively manufactured super-rough channels. Three rough surfaces are considered, all of which are obtained from computed tomography scans of additively manufactured surfaces. The roughness’ trough to peak sizes are 0.1 $h$ , 0.3 $h$ and 0.8 $h$ , respectively, where $h$ is the intended half-channel height. Each rough surface is placed opposite a smooth wall and the other two rough surfaces, leading to six rough-wall channel configurations. Two Reynolds numbers are considered, namely $Re_\tau =180$ and $Re_\tau =395$ . We show first that measurements at one streamwise and spanwise location are insufficient due to strong mean flow inhomogeneity across the entire channel, second that the logarithmic law of the wall survives despite the mean flow inhomogeneity and third that the established roughness sheltering model remains accurate.

Altland, Samuel↗
DOE OSTI2020

Understanding the Interplay between the Nonvalence and Valence States of the Uracil Anion upon Monohydration

In this work we present an ab initio investigation into the effect of monohydration on the interaction of uracil with low energy electrons. Electron attachment and photodetachment experimental studies have previously shown dramatic changes in uracil upon solvation with even a single water molecule, due to an inversion of the character of the ground state of the anion. Here we explore the interplay between the nonvalence and valence states of the uracil anion, as a function of geometry and site of solvation. Our model provides unambiguous interpretation of previous photoelectron studies, reproducing the binding energies and photoelectron images for bare uracil and a single isomer of the U•(H 2 O) 1 cluster. Furthermore, the results of this study provide insight into how electrons may attach to hydrated nucleobases. These results lay the foundations for further investigations into the effect of microhydration on the electronic structure and electron capture dynamics of nucleobases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Effects of Fulvic Acid on Uranium(VI) Sorption Kinetics

Many studies by Dr. James A. Davis and others have demonstrated the influence of natural organic matter (NOM) on metal sorption and transport behavior in subsurface environments. Overall, NOM effects on metal sorption have been ascribed to: (1) different sorption affinities of metal-ligand solution complexes compared to ‘free’ metals in solution, and (2) NOM sorption onto minerals leading to competitive sorption behavior and/or changes in mineral surface characteristics. All of these processes have the potential to decrease metal sorption onto mineral surfaces, and alter the pH-dependence of metal sorption characteristics. Maybe even more importantly, a change in metal sorption mechanisms in the presence of NOM may lead to differences in metal surface speciation and affect the potential for later metal desorption and remobilization. Furthermore, it is well-known that overall reaction rates change if the series of elemental reactions comprising a reaction pathway has been altered. In the presence of NOM, metal-ligand dissociation reactions occurring prior to the sorption of ‘free’ metal ions may slow down overall metal sorption kinetics. In addition, a competition between metals and organic ligands binding to the same reactive surface sites may potentially decrease metal sorption rates based on mass action principles. Last, the sorption kinetics of metal-ligand solution complexes may be controlled by the rates of NOM surface reactions, which could be either faster or slower than ‘free’ metal sorption depending on ligand sorption mechanisms. In this presentation, we will use a uranium(VI)-fulvic acid-silica sand system to demonstrate how kinetic sorption experiments – combined with modeling strategies – can provide a powerful tool to study changes in metal sorption mechanisms as a function of organic ligand concentrations. Based on our kinetic sorption data and modeling-based data interpretation, U(VI) sorption onto silica sand can be either faster or slower in the presence of fulvic acid compared to a ligand-free system, depending on the relative concentrations of metals, organic ligands and mineral surface sites. Supporting simulations of U(VI) solution speciation and mineral surface characterization data further allowed us to conclude that a shift occurred from competitive sorption behavior to strong influences of uranium(VI)-fulvic acid complexes with increasing ligand concentrations in this system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Response of the Atmospheric Boundary Layer to Submesoscale Sea Surface Temperature Fronts

Submesoscale sea surface temperature fronts are ubiquitous throughout much of the global ocean; however, the response of the marine atmospheric boundary layer (MABL) to the ocean submesoscale is not well understood. Here large-eddy simulation is used to explore the time-dependent response of the MABL to idealized submesoscale sea surface temperature fronts, with an emphasis on how the dynamics of the MABL determine the strength and position of gradients in wind speed and air temperature. Results suggest that horizontal mixing only becomes important in response to frontogenesis by horizontally convergent ageostrophic flows, contrary to the common assumption that the MABL response will be strongly dependent on horizontal turbulent mixing. The fronts that develop in the MABL are also associated with large vertical relative vorticity, suggesting the possibility that submesoscale fronts may induce inertial instability in the MABL. Finally, these results provide guidance for high-resolution ocean and atmosphere modeling and for interpreting observations.

58 GEOSCIENCES↗
DOE OSTI2018

The Importance of Capturing Topographic Features for Modeling Groundwater Flow and Transport in Mountainous Watersheds

Abstract Physics‐based distributed hydrological models that include groundwater are widely used to understand and predict physical and biogeochemical processes within watersheds. Typically, due to computational limitations, watershed modelers minimize the number of elements used in domain discretization, smoothing or even ignoring critical topographic features. We use an idealized model to investigate the implications of mesh refinement along streams and ridges for modeling three‐dimensional groundwater flow and transport in mountainous watersheds. For varying degrees of topographic complexity level (TCL), which increases with the level of mesh refinement, and geological heterogeneity, we estimate and compare steady state baseflow discharge, mean age, and concentration of subsurface weathering products. Results show that ignoring lower‐order streams or ridges diminishes flow through local flow paths and biases higher the contribution of intermediate and regional flow paths, and biases baseflow older. The magnitude of the bias increases for systems where permeability rapidly decreases with depth and is dominated by shallow flow paths. Based on a simple geochemical model, the concentration of weathering products is less sensitive to the TCL, partially due to the thermodynamic constraints on chemical reactions. Our idealized model also reproduces the observed emergent scaling relationship between the groundwater contribution to streamflow and drainage area, and finds that this scaling relationship is not sensitive to mesh TCL. The bias effects have important implications for the use of hydrological models in the interpretation of environmental tracer data and the prediction of biogeochemical evolution of stream water in mountainous watersheds.

Wang, Chao↗
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
DOE OSTI2024

In-beam spectroscopy reveals competing nuclear shapes in the rare isotope 62 Cr

In recent decades, rareisotope facilities have enabled the study of short-lived, neutron-rich nuclei. Their measured properties indicate that shell structure changes in the regime of unbalanced neutron-to-proton ratios compared to that of stable nuclei. In the so-called islands of inversion in the nuclear chart – around the neutron-rich nuclei 32 Mg, 42 Si and 64 Cr, for example – the textbook shell-model picture predicts the presence of spherical shapes due to the respective magic neutron numbers of 20, 28, and 40 of these nuclei. However, nuclei in these regions turn out to be deformed in their ground states. Another hallmark of these islands is shape coexistence, where a nucleus assumes different shapes with excitation energy. Furthermore, we present evidence for this phenomenon from the observation of an excited 0 + state in 62 Cr, two neutrons away from the heart of the island of inversion around neutron number N = 40. We use large-scale shell-model calculations to interpret the results and report extrapolations for the doubly-magic nucleus 60 Ca.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Examining the Chemical and Structural Properties that Influence the Sensitivity of Energetic Nitrate Esters

The sensitivity of explosives is controlled by factors that span from intrinsic chemical reactivity and chemical intramolecular effects to mesoscale structure and defects, and has been a topic of extensive study for over 50 years. Due to these complex competing chemical and physical elements, a unifying relationship between molecular framework, crystal structure, and sensitivity has yet to be developed. In order to move towards this goal, ideally experimental studies should be performed on systems with small, systematic structural modifications, with modeling utilized to interpret experimental results. Pentaerythritol tetranitrate (PETN) is a common nitrate ester explosive that has been widely studied due to its use in military and commercial explosives. We have synthesized PETN derivatives with modified sensitivity characteristics by substituting the CCH 2 ONO 2 moiety with other substituents, including CH, CNH 2 , CNH3X, CCH 3 , and PO. We relate the handling sensitivity properties of each PETN derivative to its structural properties, and discuss the potential roles of thermodynamic properties such as heat capacity and heat of formation, thermal stability, crystal structure, compressibility, and inter- and intramolecular hydrogen bonding on impact sensitivity. Reactive molecular dynamics (MD) simulations of the C/H/N/O-based PETN-derivatives have been performed under cook-off conditions that mimic those accessed in impact tests. As a result, these simulations infer how changes in chemistry affect the subsequent decomposition pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗