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 127 records · Page 7

The behavior of antioxidant irganox 1010 during the thermal degradation of a plastic bonded explosive

The effect of water concentration on the aging behavior of blend components in plastic bonded explosive (PBX) 9501 is investigated when samples were aged up to 24 months under various conditions. Additionally, the blend components studied here are: poly(urethane ester) (Estane®5703) (Estane), nitroplasticizer (NP), and antioxidant Irganox 1010 (Irg1010). The experimental results reveal that NP is prone to thermally degrading and producing H 2 O, NO x , and HNO x species, which are the predominant species to consume Irg1010 during PBX 9501 aging under inert environment. As Irg1010 is completely consumed, Estane degrades through oxidation and NP addition, in addition to well anticipated hydrolysis. The competition among hydrolysis, oxidation, and NP addition results in non-monotonical changes in the molecular weight of Estane over the aging process.

36 MATERIALS SCIENCE↗

A Hotspot’s Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material’s microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. Furthermore, the origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

36 MATERIALS SCIENCE↗

S&TR September 2025: Computing Grand Challenge Turns 20

Livermore’s Computing Grand Challenge Program enters its 20th year with more unclassified high-performance computing (HPC) power than ever before. This unique, peer-reviewed competition awards HPC allocations on top supercomputers to multidisciplinary teams with high-impact projects. The Grand Challenge encourages researchers to innovate, pushes scientific discovery to new heights, improves the Laboratory’s HPC capabilities, and extends HPC accessibility to collaborators. Awardees must adapt to successive generations of HPC hardware and learn to run simulations at scale. The feature article spotlights three Grand Challenge teams whose research broke new ground in key scientific pursuits—the essence of dark matter, explosion-generated seismic waves, and protein interactions linked to cancer—while underscoring the importance of academic partnerships and considering the program’s future.

07 ISOTOPE AND RADIATION SOURCES↗

High pressure suppression of plasticity due to an overabundance of shear embryo formation

Abstract High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is known to occur in a variety of these materials, the governing dynamics of the mechanisms are not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands, “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high-strain rate application materials such as those for high explosives and hypersonic systems.

36 MATERIALS SCIENCE↗

Nitrogen‐Nitrogen Bond Breaking in Irradiation Products of Hexanitrohexaazaisowurtzitane (CL‐20)

While the primary result of an interaction of ionizing radiation with an organic material is the ejection of electrons from molecular orbitals producing cations, the free electrons can further interact, yielding excited states; and, potentially, anions. In this work, we computed reaction barriers and energies for N−N bond breaking in the CL-20 cation, anion, and excited state to investigate if CL-20 degradation is accelerated after radiation exposure. We focused on N−N bond dissociation because it is the rate-determining step in the thermal degradation of CL-20. We found that N−N cleavage rapidly takes place in the CL-20 cation and anion with greatly reduced reaction energies as compared to CL-20. We also outlined a potential path for photodissociation of the N−N bond. While CL-20 is kinetically stable, initial degradation occurs readily in its irradiation products. In conclusion, this is of importance for the performance of the explosive after high-dose exposure and influences aging if CL-20 is irradiated at a low dose over extended periods of time.

CL-20 anion↗

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dependence of the Elastic Stiffness Tensors of PETN, α‐RDX, γ‐RDX, ϵ‐RDX, ϵ‐CL‐20, DAAF, FOX‐7, and β‐HMX on Hydrostatic Compression

Abstract The dependence of the components of the elastic stiffness tensors (or elastic constants) of the organic explosives PETN, RDX, CL‐20, DAAF, FOX‐7, and HMX on hydrostatic pressure up to 10 GPa have been computed using dispersion‐corrected density functional theory. We report the evolution of lattice parameters and the non‐zero stiffnesses for the tetragonal, orthorhombic, and monoclinic crystal symmetries. Linear and quadratic dependencies of the components of the elastic stiffness tensors on volumetric compression and hydrostatic pressure are tabulated for use in single crystal plasticity models.

36 MATERIALS SCIENCE↗

The complex polyploid genome architecture of sugarcane

Sugarcane, the world’s most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype (‘monoploid’) representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗

BioTransformer 3.0 – A Web Server for Accurately Predicting Metabolic Transformation Products

BioTransformer 3.0 is a freely available web server that supports accurate, rapid and comprehensive in silico metabolism prediction. It combines machine learning approaches with a rule-based system to predict small-molecule metabolism in human tissues, the human gut as well as the external environment (soil and water microbiota). Simply stated, BioTransformer takes a molecular structure as input (SMILES or SDF) and outputs an interactively viewable/sortable table of the predicted metabolites or transformation products (SMILES, PNG images) along with the enzymes that are predicted to be responsible for those reactions and richly annotated downloadable files (CSV and JSON). The entire process typically takes a few seconds. Previous versions of BioTransformer focused exclusively on predicting the metabolism of xenobiotics (such as plant natural products, drugs, cosmetics and other synthetic compounds) using a limited number of pre-defined steps and somewhat limited rule-based methods. BioTransformer 3.0, uses much more sophisticated methods and incorporates new databases, new constraints and new prediction modules to not only more accurately predict the metabolic transformation products of exogenous xenobiotics but also the transformation products of endogenous metabolites, such as amino acids, peptides, carbohydrates, organic acids, and lipids. BioTransformer 3.0 can also support customized sequential combinations of these transformations along with multiple iterations to simulate multi-step human and/or environmental biotransformation events. Performance tests indicate that BioTransformer 3.0 is 40-50% more accurate, much less prone to combinatorial “explosions” and far more comprehensive in terms of metabolite coverage/capabilities than previous versions of BioTransformer.

59 BASIC BIOLOGICAL SCIENCES↗

X-ray multiphoton-induced Coulomb explosion images complex single molecules

Following structural dynamics in real time is a fundamental goal towards a better understanding of chemical reactions. Recording snapshots of individual molecules with ultrashort exposure times is a key ingredient towards this goal, as atoms move on femtosecond (10 -15 s) timescales. For condensed-phase samples, ultrafast, atomically resolved structure determination has been demonstrated using X-ray and electron diffraction. Pioneering experiments have also started addressing gaseous samples. However, they face the problem of low target densities, low scattering cross sections and random spatial orientation of the molecules. Therefore, obtaining images of entire, isolated molecules capturing all constituents, including hydrogen atoms, remains challenging. Here we demonstrate that intense femtosecond pulses from an X-ray free-electron laser trigger rapid and complete Coulomb explosions of 2-iodopyridine and 2-iodopyrazine molecules. We obtain intriguingly clear momentum images depicting ten or eleven atoms, including all the hydrogens, and thus overcome a so-far impregnable barrier for complete Coulomb explosion imaging—its limitation on molecules consisting of three to five atoms. In combination with state-of-the-art multi-coincidence techniques and elaborate theoretical modelling, this allows tracing ultrafast hydrogen emission and obtaining information on the result of intramolecular electron rearrangement. Our work represents an important step towards imaging femtosecond chemistry via Coulomb explosion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Strong-field ionization of water: Nuclear dynamics revealed by varying the pulse duration

Polyatomic molecules in strong laser fields can undergo substantial nuclear motion within tens of femtoseconds. Ion imaging methods based on dissociation or Coulomb explosion therefore have difficulty faithfully recording the geometry dependence of the field ionization that initiates the dissociation process. Here we compare the strong-field double ionization and subsequent dissociation of water (both H 2 O and D 2 O) in 10-fs and 40-fs 800-nm laser pulses. We find that 10-fs pulses turn off before substantial internuclear motion occurs, whereas rapid internuclear motion can take place during the double ionization process for 40-fs pulses. The short-pulse measurements are consistent with a simple tunnel ionization picture whose predictions help interpret the motion observed in the long-pulse measurements.

74 ATOMIC AND MOLECULAR PHYSICS↗

The promising role of proteomes and metabolomes in defining the single-cell landscapes of plants

The plant community has a strong track-record of RNA sequencing technology deployment, which combined with the recent advent of spatial platforms (e.g., 10x genomics), has resulted in an explosion of outstanding single cell and nuclei datasets that can be put in an in situ context within tissues (e.g., a cell atlas)1. In the genomics era, application of proteomics technologies in the plant sciences has always trailed behind that of RNA sequencing technologies, largely due to accessibility, ease-of-use and access to expertise along with depth of analysis benefits. On the other hand, the use of early analytical tools for characterizing small molecules (metabolites) from plant systems predates nucleic acid sequencing and proteomics analysis2, as the search for plant-based natural products has played a significant role in improving human health throughout history. However, the employment of proteomics and metabolomics assays for characterizing plant cell processes now remains significantly behind transcriptional approaches, even though both provide a direct functional readout of cell states and phenotypes.

Anderton, Christopher R. [BATTELLE (PACIFIC NW LAB↗

The mechanical properties of Kel-F 800 (FK-800) as a function of crystallinity

Kel-F 800 is a copolymer of chlorotrifluoroethylene PTFE (75 wt. %) and vinylidene fluoride PVDF (25 wt. %). It has previously been used as a PBX binder for insensitive explosives such as PBX 9502 and LX-17. 3M started production of Kel-F 800 in 1957 and small-scale batches continued to be made until 2002 when production ceased due to environmental concerns regarding one of the emulsifiers used during production. Around 2000 the Kel-F 800 name was changed to FK-800 to avoid trademark concerns because rights to produce another polymer with a similar tradename (Kel-F 81) had been sold to another manufacturer. The Kel designation came from the original manufacturer of PCTFE (Kel-F 81), the Kellog company. To avoid confusion this document will only refer to Kel-F 800. In 2006, production of small-scale batches of Kel-F 800 was started again by 3M in response to customer enquiries. This new material, the first blended batch is referred to as LOT 1, was produced with a different emulsifier than used previously. Because Kel-F 800 is made in a small batch reactor, considerable variation in crystallinity can be expected from lot to lot and year to year. In many ways, this is not significant since the material is dissolved in a solvent (often MEK, ethylmethyl ketone or ethyl acetate) for PBX production purposes. This destroys the as received crystallinity and the resulting crystallinity in the processed material is a function of polymer molecular weight and thermal history. Producing large billets of Kel-F 800 from solvent extraction is not practical and so a compression molding technique has been used above the melting temperature. This method also removes residual crystallinity from the supplied granules. The molecular weight of a polymer can be estimated by several techniques, the most common being gel permittivity chromatography (GPC), size exclusion chromatography (SEC) and shear rheometry measurements of polymer/solvent solutions. Changes in molecular weight will affect the crystallization rate and the maximum crystallinity reached for a specific thermal history. Both references agree that the new LOT 1 material molecular weight falls within the deviation found from averaging previous historical lots of Kel-F 800.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fate of multiparticle resonances: From Q-balls to 3 He droplets

We study a system of N nonrelativistic particles which form a near-threshold resonance. Assuming no subset of these particles can form a bound state, the resonance can only decay through an “explosion” into N particles. We find that the decay width of the resonance scales as E Δ –5/2 in the limit when the energy E of the resonance goes to zero, where Δ is the ground-state energy of a system of N particles in a spherical harmonic trap with unit frequency. Here, the formula remains valid when some pairs of final particles have zero-energy s-wave resonance, but the Efimov effect is not present. In the limit of large N, we show that the final particles follow a Maxwell-Boltzmann distribution if they are bosons and a semicirclelike law if they are fermions. We expect our general result to be applicable to various systems that exist in nature. In particular, we argue that metastable 3 He droplets exist with the lifetime varying over many orders of magnitude ranging from a fraction of a nanosecond to values greatly exceeding the age of the Universe.

74 ATOMIC AND MOLECULAR PHYSICS↗

Preshot Report for Carbonara: Dynamics of a Boron Carbide Sandwich

The main purpose of the Carbonara experiment is to provide much-needed data to inform physics modeling of boron carbide (B 4 C) subjected to a sweeping high explosive (HE) detonation. The proposed experiment will include three separate shots, each employing a rectangular slab of HE (4” × 3” × 1”) sandwiched between two wedges of B 4 C, see Fig.1. In each case the top B 4 C wedge tapers from a height of 2.5 cm to 0.5 cm, while the bottom wedge tapers from 1.5 cm to 0.5 cm. The wedges are 4” in the direction of the proton beam, and 3” in the direction transverse to the beam. Three shots are currently planned, each one using a different HE: PBX 9501, PBX 9502, and PBX 9701. The HE will be initiated using a line wave generator (LWG) to preserve planar symmetry, in conjunction with a Comp. B booster to smooth out the profile for the initial condition on the main charge and to aid in initiation of the PBX 9502 shot. In each case the booster spans the entire length and height of the main HE slab and extends 0.5” in toward the LWG.

36 MATERIALS SCIENCE↗

Multiparticle cumulant mapping for Coulomb explosion imaging: Calculations and algorithm

We present a versatile cumulant mapping algorithm for analyzing correlated particle emission, offering insights into complex electronic and nuclear dynamics. Recently, we have demonstrated the use of cumulant mapping to extract information-rich correlations between the momenta of multiple fragments produced in Coulomb explosion imaging experiments [C. Cheng et al., Phys. Rev. Lett. 130, 093001 (2023)]. We define cumulant mapping in terms of histograms, enabling fast computation of linear (additive) observables. However, applying the same algorithm to nonlinear (nonadditive) observables poses challenges, as the computation time of conventional estimators scales nonlinearly with data size. To overcome this, we develop estimators and an accompanying algorithm to enable computationally efficient estimation of the cumulant of interest. Comparisons of computation times and signal-to-noise ratios reveal the superior performance of our approach. This method is demonstrated on the (D+, D+, C+, O+) dissociation channel of CD 2 ⁢O 4+ produced in a strong-field ionization experiment. Additionally, Poisson statistics are used to simulate the two methods and provide insights into the efficiency of our algorithm. The proposed methodology unlocks efficient computation of cumulant mapping for a broader range of complex systems and observables, such as the laser pulse dependence of ionization dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Initial decomposition mechanisms of 2,4,6-triamino-1,3,5-trinitrobenzene (TATB) and their kinetic isotope effect

2,4,6-triamino-1,3,5-trinitrobenzene (TATB) is an insensitive High Explosive (HE) that is widely studied to better understand the physical properties of safety and sensitivity of HE. A dominant initial decomposition mechanism of TATB is believed to be a dehydration reaction that forms mono- and di-furazans, although other mechanisms have been reported. In this work, seven initial decomposition mechanisms were modeled with ab initio simulations to calculate its free energy barriers, decomposition rates, and kinetic isotope effects. The energy barrier for mono-benzofurazan mechanisms was found to be high, >61 kcal/mol in the gas phase; however, the reaction energy can decrease significantly in a disordered condensed state. The predicted kinetic isotope effect ratio of the furazan mechanism was found to be k h /k d ≈ 1.41 at 60K, in agreement with the experiment. The NO 2 scission mechanism was found to be an entropy-driven mechanism because the free energy barrier decreased significantly with temperature, making it the most energetically favorable mechanism at high temperatures in the gas phase. So the results provide a better understanding of the atomistic decomposition mechanisms of TATB and may be useful for improving models of safety and sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Constraining neutron-star matter with microscopic and macroscopic collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not probed only in astrophysical observations, but also in terrestrial heavy-ion collision experiments. Here we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars and from heavy-ion collisions of gold nuclei at relativistic energies with microscopic nuclear theory calculations to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent observations by the Neutron Star Interior Composition Explorer mission. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗