Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular explosives”

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 109 records · Page 6

Spatial and temporal control of lysis by the lambda holin

The infection cycle of phage λ terminates in lysis mediated by three types of lysis proteins, each disrupting a layer in the bacterial envelope: the S105 holin, the R endolysin, and the Rz/Rz1 spanin complex targeting the inner membrane, cell wall or peptidoglycan, and the outer membrane, respectively. Video microscopy has shown that in most infections, lysis occurs as a sudden, explosive event at a cell pole, such that the initial product is a less refractile ghost that retains rod-shaped morphology. Here, we investigate the molecular basis of polar lysis using time-lapse fluorescence microscopy. The results indicate that the holin determines the morphology of lysis by suddenly forming two-dimensional rafts at the poles about 100 s prior to lysis. Given the physiological and biochemical similarities between the lambda holin and other class I holins, dynamic redistribution and sudden concentration may be common features of holins, probably reflecting the fitness advantage of all-or-nothing lysis regulation

Microbiology↗

Atomic and molecular supernovae

Atomic and molecular physics of supernovae is discussed with an emphasis on the importance of detailed treatments of the critical atomic and molecular processes with the best available atomic and molecular data. The observations of molecules in SN 1987A are interpreted through a combination of spectral and chemical modelings, leading to strong constraints on the mixing and nucleosynthesis of the supernova. The non-equilibrium chemistry is used to argue that carbon dust can form in the oxygen-rich clumps where the efficient molecular cooling makes the nucleation of dust grains possible. For Type Ia supernovae, the analyses of their nebular spectra lead to strong constraints on the supernova explosion models.

Liu, Weihong↗

Electron–Ion Covariance Reveals Geometry-Resolved Inner-Shell Spectra in CS 2 Photodissociation

The chemical shifts of inner-shell atomic orbitals are highly sensitive to chemical bonding and molecular structure. In evolving systems, however, the spectra of distinct molecular species and geometries overlap, obscuring the underlying chemical dynamics. Here we demonstrate the use of electron–ion covariance analysis to combine the structural sensitivity of Coulomb explosion imaging with inner-shell spectroscopy, yielding geometry-specific spectra of transient and product species. We apply this approach to the excited state dynamics of CS 2 probed by ionization above the S 2p edge. Electron–ion covariance with time- and momentum-selected S + and S 2+ ions isolates distinct S 2p photoelectron spectra for ground-state CS 2 , bent photoexcited CS 2 , the CS photoproduct, and bare atomic sulfur ─ species whose spectra overlap strongly in the channel-averaged measurement. Clear chemical shifts are observed in the covariance photoelectron spectrum for each of these species, all of which are consistent with high-level calculations. Here, by extracting the atomic S contribution to the photoelectron spectrum in a finely time-resolved manner we can disentangle this contribution to the overall time-resolved photoelectron spectrum as the photodissociation proceeds. These results demonstrate the promise of electron–ion covariance as a general approach to geometry-resolved inner-shell spectroscopy, opening a route to tracking structural evolution through chemical shifts in complex photoexcited molecules.

Ionization↗

Data Science Enabled Enabled Discovery of Superconductors (Final Progress Report)

This Final Technical Report describes efforts by 4 PIs at the University of Florida (Peter Hirschfeld, Richard Hennig, Greg Stewart and James Hamlin), over the period September 2019-August 2023, to use data science and machine learning techniques to discover new conventional superconductors. The PIs constructed a discovery loop with two theorists and two experimentalists to: develop algorithms to machine learn descriptors correlating strongly with the critical temperature Tc (PI's Peter Hirschfeld, UF Physics and Richard Hennig, UF Materials Science and En), synthesize and measure properties of promising materials, and feed back the knowledge gained into the prediction algorithm. This work was motivated by the theoretical prediction and experimental discovery of high-pressure, high-pressure hydride superconductors, and to find ways to recreate the high critical temperatures in these systems at ambient pressure. Highlights from the grant include: 1) a new equation for Tc in terms of moments of the electron-phonon spectral function, improving on the so-called Allen-Dynes equation (1975); 2) study of the metastable A15 superconductor Nb3Si, formed under explosive compression at ~1000GPa to determine the kinetic barrier to the ground state structure; 3) the development of ultra-fast machine-learned atomic potentials for molecular dynamics, and 4) the discovery of superconductivity at 19K in WB2 arising from metastable defect structures in the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗

Sulfur and halogen chemistry of the stratosphere and of volcanic eruption plumes

The major eruption of Volcan de Fuego caused an enhancement of 1.6 x 10 to the 9 kg of SO4(-) in the stratosphere 6 months after the eruption. Measurable changes in nitric acid vapor and particulate chloride were not observed. Subsequent data suggests that HCl was not injected into the stratosphere. The first simultaneous measurement of halogen and sulfur content of gases and particles in explosive eruption plumes is reported. The average percents of elements in particles were sulfur, 2.5 + or - 2.1; chlorine, 18 + or - 12; and fluorine, 38 + or - 29. The average molecular ratio of HCl/SO2 and HCl/HF were 0.41 + or - 0.26 and 14 + or - 12 in the Guatemalan eruption plumes.

Lazrus, A. L.↗

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗

Comprehensive Data Analysis and Analytic Method Development for PBX9501 Material (FY2019 Annual Report of Aging and Lifetimes Program)

In the polymer bonded explosive (PBX) 9501, a binder consisting of 2.5 wt% nitroplasticizer (NP) and 2.5 wt% Estane R 5703 (Estane) is combined with 94.9 wt% of the high-explosive HMX and 0.1 wt% stabilizer. Because of its flexibility and tensile strength, this binder lowers the sensitivity and improves the manufacturability of PBX 9501. However, like many plasticizers comprised of low molecular weight components, NP has a tendency to diffuse out of the PBX 9501 matrix and can decompose at moderate temperatures into reactive byproducts, such as NO, NO 2 , H 2 O, and HNO x . Through oxidation and hydrolysis, these molecules can further degrade NP and Estane, ultimately degrading the properties of PBX 9501. While gaseous molecules can readily diffuse out of the PBX 9501 charges, the intermediates with lower volatility are more likely trapped in the condensed phase inside the charges and, in a closed system, will co-exist with Estane for an extended period of time creating an ongoing reactive environment. Understanding the rates of formation for these volatiles and intermediates during the NP degradation is therefore a critical prerequisite for understanding the long-term stability of polymeric binders used in munition systems. To work toward this ultimate goal, in the past year, we have conducted comprehensive studies on the physical properties of NP and finished three-year long aging experiment to understand the aging behavior of NP under thermal treatment. Important results are summarized in this annual report. Although NP is widely used in the DOE complex, their physical properties are rather scattered and inconsistent in the open literature. For example, there are at least two widely different values, 14.5°C and -15°C, for the melting point of NP. Although it is known that NP is a eutectic 50:50 mixture of BDNPA and BDNPF, their eutectic phase diagram is not well documented. Furthermore, the effect of temperature on the miscibility between water and eutectic BDNPA/F mixture is rarely reported. To fill these knowledge gaps, a large set of BDNPA/F mixtures with BDNPA concentration ranging from 0 to 100 wt% was analyzed using DSC techniques. In addition to determining the eutectic melt point as -25°C, a phase diagram of the BDNPA/F system was constructed from -30°C to 45°C. With this phase diagram, the phase transition temperatures and composition can be readily found for the BDNPA/F mixtures with various mass ratios.

36 MATERIALS SCIENCE↗

Optical Microsystems for Molecular Sensing (CRADA Final Report)

As part of the Cyclotron Road program, NuSight Photonics sought to implement and investigate an optical Raman spectroscopy instrument for material and molecular analysis. Raman spectroscopy is widely used for quick and non-invasive material measurement when accurate results are needed – e.g., for monitoring manufacturing process in pharmaceutical plants, or identifying illicit substances and explosives by Law Enforcement. This technique is also growing rapidly in popularity as it is applied successfully to a broad range of applications, from pathology and non-invasive early cancer detection, to food and drug safety. The biggest challenge for Raman spectroscopy to achieve its full potential is the hardware cost, ranging from $\$ $25k to $\$ $250k, with the lower cost solutions significantly limited in sensitivity—or detection limit. NuSight Photonics aimed to tackle this obstacle through a novel hardware approach: Swept-Source Raman Spectroscopy.

36 MATERIALS SCIENCE↗

Trace explosive residue detection of HMX and RDX in post-detonation dust from an open-air environment

Explosives are often used in industry, geology, mining, and other applications, but it is not always clear what remains after a detonation or the fate and transport of any residual material. The goal of this study was to determine to what extent intact molecules of high explosive (HE) compounds are detectable and quantifiable from post-detonation dust and particulates in a field experiment with varied topography. We focused on HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane), which is less studied in field detonation literature, as the primary explosive material and RDX (1,3,5-Trinitroperhydro-1,3,5-triazine) as the secondary material. The experiment was conducted at Site 300, Lawrence Livermore National Laboratory’s Experimental Test Site, in California, USA. Two 20.4 kg and one 40.8 kg above ground explosions (primarily comprised of LX-14, an HMX-based polymer-bonded high explosive) were detonated on an open-air firing area on separate days. The complex terrain of the firing area (e.g., buildings, berm, low-height obstacles) was advantageous to study HE deposition in relation to plume dynamics. Three types of samples were collected up to 100 m away from each shot: surface swipes of aluminum plates, surface swipes of fixed objects, and filters from air samples. We used atmospheric flow tube-mass spectrometry (AFT-MS) to quantify picogram levels of molecular residue of HE material in the post-detonation dust. An aliquot of sample extract in methanol (e.g., 1 µL of 0.5 mL) was placed onto a resistive material and then thermally desorbed into the AFT-MS. We successfully detected and quantified both HMX and RDX in many of the samples. Based on mass (pg) detected and solution dilution, we back-calculated the mass collected on the swipe or filter (ng per sample). The aerial distribution of molecular residue was consistent with the path of the plume, which was strongly determined by wind speed and direction at the time of each shot. The quantity of material detected appeared to correlate more with distance from the shot and the wind conditions than with shot size. This study demonstrates that the picogram detection levels of AFT-MS are well-suited for quantification of analytes (e.g., HMX and RDX) in environmental samples.

atmospheric flow tube-mass spectrometry (AFT-MS), ↗

Hydrogen Molecules in SN 1987A

The observations of CO and SiO in the infrared spectrum of SN 1987A clearly indicate that molecules can form in the debris of a supernova explosion. Since H2 is not easily observable we compute its abundance theoretically. For conditions typical of the inner (v less than 2500 km/s) envelope of SN 1987A, the fraction of H that is in molecular form rises to approx. 1% by t approx. 800 days. For t less than 500 days the formation is dominated by the gas-phase reactions H + H(+) yields H2(+) + hv; H2(+) + H yields H2 + H(+). Thereafter, the formation is dominated by the reactions H + e yields H(-) + hv; H(-) + H yields H2 + e. At early times the H(-) may absorb approx. 10%-30% of visible photons, contributing to the apparent paucity of H alpha emission. For t greater than 1000 days the abundance of H2 'freezes out' due to the slowing of all reactions. The opacity of the supernova envelope in the range 912 less than lambda less than or approx. equal to 1400 A (the upper limit depending on temperature) is dominated by resonance scattering in the Lyman and Werner bands of H2. The resulting fluorescence emission bands of H2 in the range 1150 less than lambda less than 1650 A may be observable in the UV spectra of supernovae at late times.

Culhane, Michael↗

Visualizing the strong field–induced molecular breakup of C 60 via x-ray diffraction

Laser-driven dynamics in polyatomic molecules poses a complex many-body problem. Understanding intense light-matter interaction is crucial for steering intramolecular quantum dynamical processes. Here, we record time-resolved x-ray diffraction images of C 60 molecules during and after their interaction with intense near-infrared fields, giving direct access to structural changes of the molecules and their fragmentation in real time. Tuning the intensity of the excitation pulses, we uncover a transition from a weak-field regime of excited but stable molecules to a high-field regime dominated by Coulomb explosion. In the transition region, the molecules expand by up to 50% of their initial size within just 140 fs, with major fragmentation only setting in afterward. This work demonstrates that x-ray diffractive imaging is capable of retrieving time-resolved structural information of large molecules reshaped by intense laser fields. Laser-driven fragmentation is a first step toward observing molecular processes modified by laser fields of increasing intensity.

Schnorr, Kirsten [Paul Scherrer Inst. (PSI), Villi↗

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗

Airborne spectrophotometry of SN 1987A from 1.7 to 12.6 microns - Time history of the dust continuum and line emission

Spectrophotometric observations of SN 1987A from the Kuiper Airborne Observatory are presented for five epochs at 60, 260, 415, 615, and 775 days after the explosion. The low-resolution (lambda/Delta lambda = 50-100) spectra of SN 1987A are combined with data from other wavelengths to model the continuum, subtract the continuum from the spectra to determine line strengths and reveal molecular bands, separate the atomic continuum radiation from the dust continuum, and derive constraints on the grain temperatures and optical depths. A scenario for the evolution of SN 1987A and that of the ejecta from which it arises is obtained on the basis of the analysis of the continuum emission.

Wooden, Diane H.↗

Regolith Volatile Characterization (RVC) in RESOLVE

Resource investigation in the lunar poles is of importance to the potential impact of in-situ resource utilization (ISRU). The RESOLVE project developed a payload to investigate the permanently shadowed areas of the lunar poles and demonstrate ISRU technology. As a part of the RESOLVE project, the regolith volatile characterization (RVC) subsystem was designed to examine the release of volatiles from sample cores. The test sample was heated in the reactor to release the volatiles where they were analyzed with gas chromatography. Subsequently, the volatile sample was introduced into the lunar water resource demonstration (LWRD) subsystem where the released hydrogen and water were selectively captured. The objective of the Regolith Volatile Characterization (RVC) subsystem was to heat the crushed core sample and determine the desorption of volatile species of interest. The RVC subsystem encompasses the reactor and the system for volatile analysis. The system was designed to analyze H2, He, CO, CO2, N2, 02, CH4, H2S and H2O. The GC chosen for this work is a Siemens MicroSAM process GC with 3 columns and 8 TCD detectors. Neon was chosen as the carrier gas to enhance the analysis of hydrogen and helium.The limit of detection for the gases is approx.1000ppm for H2, CO. CO2 , N2, O2 and H2 S. The limit of detection for CH4 is approx.4000ppm and the water limit of detection is -10000 ppm with a sample analysis time of 2-3 minutes. These values (with the exception of water and H2S) were determined by dilution of a six gas mixture from Scott Gas (5% CO2, CO, O2, N2, 4% CH4 and H2) using mass flow controllers (MFC5). Water was calibrated at low levels using an in house relative humidity (RH) generator. H 2S and high concentrations of H2 were calibrated by diluting a pure stream of gas with MFCs. Higher concentrations of N2 and 02 were calibrated using Air again diluting with MFCs. There were three modification goals for the GC in EBU2 that would allow this process GC to be used in the field demo for RESOLVE. The first modification was to decrease the weight associated with the GC, this included eliminating the explosion proof case (Figure 1) and replacing it with a lightweight case as well as using an on board COPV tank for the neon carrier gas. The next goal was to add a second oven for the molecular sieve column to allow for dual temperature control during GC operation; the separation of hydrogen and helium is optimum at lower temperatures while the water analysis required higher temperatures creating a competing design requirement. The second oven also allows a lower limit of detection for water quantification and avoids the possibility of water condensing in the GC which could ruin the column characteristics. The final goal was to modify the column arrangement to optimize the system for our specific application. Figure 2 shows the internal details of the module optimized optimized for our field application. The modifications and performance of the gas analysis system will be discussed in detail.

Captain, Janine↗