Ranking the Drop-Weight Impact Sensitivity of Common Explosives Using Arrhenius Chemical Rates Computed from Quantum Molecular Dynamics Simulations

Understanding Trigger Linkage Dynamics in Energetic Materials Using Mixed Picramide Nitrate Ester Explosives

The ability to predict the handling sensitivity of new organic energetic materials has been a longstanding goal. We report the synthesis and characterization of six new nitropicramide energetic materials with mixed functional groups that mimic known explosives such as nitroglycerin, erythritol tetranitrate (ETN), and pentaerythritol tetranitrate (PETN). The molecules have been studied theoretically using quantum molecular dynamics (QMD) simulations and density functional theory (DFT) calculations to identify the weakest bond in the reactants - the trigger-linkages - which control handling sensitivity, and to quantify their specific enthalpies of explosion. In good accord with the drop weight impact sensitivity data, our calculations predict that the sensitivities of the molecules are very similar owing to the small variations of the energy output and rates of trigger linkage rupture. In addition, both the QMD and DFT calculations point to the nitropicramide N-NO2 bonds as the trigger linkages rather than the more typical O-NO2 bonds. We propose that the switch of the trigger linkage from the nitrate esters to the nitramine groups arises from the strongly electron withdrawing character of the adjacent trinitrobenzene groups.

An Integrated Experimental and Modeling Approach for Assessing High-Temperature Decomposition Kinetics of Explosives

We present a new integrated experimental and modeling effort that assesses the intrinsic sensitivity of energetic materials based on their reaction rates. The High Explosive Initiation Time (HEIT) experiment has been developed to provide a rapid assessment of the high-temperature reaction kinetics for the chemical decomposition of explosive materials. This effort is supported theoretically by quantum molecular dynamics (QMD) simulations that depict how different explosives can have vastly different adiabatic induction times at the same temperature. In this work, the ranking of explosive initiation properties between the HEIT experiment and QMD simulations is identical for six different energetic materials, even though they contain a variety of functional groups. We have also determined that the Arrhenius kinetics obtained by QMD simulations for homogeneous explosions connect remarkably well with those obtained from much longer duration one-dimensional time-to-explosion (ODTX) measurements. Kinetic Monte Carlo simulations have been developed to model the coupled heat transport and chemistry of the HEIT experiment to explicitly connect the experimental results with the Arrhenius rates for homogeneous explosions. These results confirm that ignition in the HEIT experiment is heterogeneous, where reactions start at the needle wall and propagate inward at a rate controlled by the thermal diffusivity and energy release. Overall, this work provides the first cohesive experimental and first-principles modeling effort to assess reaction kinetics of explosive chemical decomposition in the subshock regime and will be useful in predictive models needed for safety assessments.

Pressure and temperature effects on the Raman spectra of LLM-105

Currently, only one crystal structure of LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide) (P21/n) has been discovered, and there are still debates on its phase transition point and phase diagram. Based on previous work, we performed crystal structure, Raman spectra, and vibrational properties calculations on LLM-105 crystal. Our results indicate that the crystal structure of LLM-105 remains stable until compressed to 49 GPa, beyond which it may undergo two phase transitions at pressure intervals of 49.0-49.1 GPa and 51.4-51.5 GPa, respectively. Analysis of Raman shift results suggests that these two phase transitions may be reversible, with an intermediate phase possibly acting as a transition phase. Additionally, based on the quasi-harmonic approximation, we fitted the experimental data of LLM-105 lattice expansion state, obtaining the volume at zero pressure and using it for Raman spectra calculations. The results demonstrated the accuracy of this quasi-harmonic approximation method in describing the redshift of Raman peaks during the heating process and the excitation ratio of Raman peaks in different wavenumber ranges.

Insensitive High-Energy Density Materials Based on Azazole-Rich Rings: 1,2,4-Triazole N-Oxide Derivatives Containing Isomerized Nitro and Amino Groups

It is an arduous and meaningful challenge to design and develop new energetic materials with lower sensitivity and higher energy. How to skillfully combine the characteristics of low sensitivity and high energy is the key problem in designing new insensitive high-energy materials. Taking a triazole ring as a framework, a strategy of N-oxide derivatives containing isomerized nitro and amino groups was proposed to answer this question. Based on this strategy, some 1,2,4-triazole N-oxide derivatives (NATNOs) were designed and explored. The electronic structure calculation showed that the stable existence of these triazole derivatives was due to the intramolecular hydrogen bond and other interactions. The impact sensitivity and the dissociation enthalpy of trigger bonds directly indicated that some compounds could exist stably. The crystal densities of all NATNOs were larger than 1.80 g/cm3, which met the requirement of high-energetic materials for crystal density. Some NATNOs (9748 m/s for NATNO, 9841 m/s for NATNO-1, 9818 m/s for NATNO-2, 9906 m/s for NATNO-3, and 9592 m/s for NATNO-4) were potential high detonation velocity energy materials. These study results not only indicate that the NATNOs have relatively stable properties and excellent detonation properties but also prove that the strategy of nitro amino position isomerization coupled with N-oxide is an effective means to develop new energetic materials.

Chemical Descriptors for a Large-Scale Study on Drop-Weight Impact Sensitivity of High Explosives

The drop-weight impact test is an experiment that has been used for nearly 80 years to evaluate handling sensitivity of high explosives. Although the results of this test are known to have large statistical uncertainties, it is one of the most common tests due to its accessibility and modest material requirements. In this paper, we compile a large data set of drop-weight impact sensitivity test results (mainly performed at Los Alamos National Laboratory), along with a compendium of molecular and chemical descriptors for the explosives under test. These data consist of over 500 unique explosives, over 1000 repeat tests, and over 100 descriptors, for a total of about 1500 observations. We use random forest methods to estimate a model of explosive handling sensitivity as a function of chemical and molecular properties of the explosives under test. Our model predicts well across a wide range of explosive types, spanning a broad range of explosive performance and sensitivity. We find that properties related to explosive performance, such as heat of explosion, oxygen balance, and functional group, are highly predictive of explosive handling sensitivity. Yet, models that omit many of these properties still perform well. Our results suggest that there is not one or even several factors that explain explosive handling sensitivity, but that there are many complex, interrelated effects at play.

Radiolytic degradation of dodecane substituted with common energetic functional groups†

Explosives exist in and are expected to withstand a variety of harsh environments up to and including ionizing radiation, though little is known about the chemical consequences of exposing explosives to an ionizing radiation field. This study focused on the radiation-induced chemical changes to a variety of common energetic functional groups by utilizing a consistent molecular backbone. Dodecane was substituted with azide, nitro, nitrate ester, and nitramine functional groups and γ-irradiated with ⁶⁰Co in order to study how the functional group degraded along with what the relative stability to ionizing radiation was. Chemical changes were assessed using a combination of analysis techniques including: nuclear magnetic resonance (NMR) spectroscopy, gas chromatography of both the condensed and gas phases, Raman spectroscopy, and Fourier transform infrared (FTIR) spectroscopy. Results revealed that much of the damage to the molecules was on the energetic functional group and often concentrated on the trigger linkage, also known as the weakest bond in the molecule. The general trend from most to least susceptible to radiolytic damage was found to be D–ONO₂ → D–N₃ → D–NHNO₂ → D–NO₂. These results also appear to be in line with the relative stability of these functional groups to things such as photolysis, thermolysis, and explosive insults.

The relative radiolytic stability of dodecane functionalized with common energetic functional groups was explored with gamma irradiation and probed by various analytical techniques.

Identifying the Molecular Properties that Drive Explosive Sensitivity in a Series of Nitrate Esters

Energetic materials undergo hundreds of chemical reactions during exothermic runaway, generally beginning with the breaking of the weakest chemical bond, the "trigger linkage." Herein we report the syntheses of a series of pentaerythritol tetranitrate (PETN) derivatives in which the energetic nitrate ester groups are systematically substituted by hydroxyl groups. Because all the PETN derivatives have the same nitrate ester-based trigger linkages, quantum molecular dynamics (QMD) simulations show very similar Arrhenius kinetics for the first reactions. However, handling sensitivity testing conducted using drop weight impact indicates that sensitivity decreases precipitously as nitrate esters are replaced by hydroxyl groups. These experimental results are supported by QMD simulations that show systematic decreases in the final temperatures of the products and the energy release as the nitrate ester functional groups are removed. To better interpret these results, we derive a simple model based only on the specific enthalpy of explosion and the kinetics of trigger linkage rupture that accounts qualitatively for the decrease in sensitivity as nitrate ester groups are removed.

Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes

The photolytic stability of explosives and energetic functional groups is of importance for those who regularly handle or are exposed to explosives in typical environmental conditions. This study models the photolytic degradation of dodecane substituted with various energetic functional groups: azide, nitro, nitrate ester, and nitramine. For the studied molecules, it was found that excitons localize on the energetic functional group, no matter where they were initially formed, and thus, the predominant degradation pathway involves the degradation of the energetic functional group. The relative trends for both 4 and 8 eV excitation energies followed with what is expected from the relative stability of the energetic functional groups to thermal and sub-shock degradation. The one notable exception was the azide functional group; more work should be done to further understand the photolytic effects on the azide functional group.

Understanding Explosive Sensitivity with Effective Trigger Linkage Kinetics

We present a simple linear model for ranking the drop weight impact sensitivity of organic explosives that is based explicitly on chemical kinetics. The model is parameterized to specific heats of explosion, Q, and Arrhenius kinetics for the onset of chemical reactions that are obtained from gas-phase Born-Oppenheimer molecular dynamics simulations for a chemically diverse set of 24 molecules. Reactive molecular dynamics simulations sample all possible decomposition pathways of the molecules with the appropriate probabilities to provide an effective reaction barrier. In addition, the calculations of effective trigger linkage kinetics can be accomplished without prior physical intuition of the most likely decomposition pathways. We found that the specific heat of explosion tends to reduce the effective barrier for decomposition in accordance with the Bell-Evans-Polanyi principle, which accounts naturally for the well-known correlations between explosive performance and sensitivity. Our model indicates that sensitive explosives derive their properties from a combination of weak trigger linkages that react at relatively low temperatures and large specific heats of explosion that further reduce the effective activation energy.

Design and Computational Insight on Two Novel CL-20 Analogues, BNMTNIW and BNIMTNIW: High Performance Energetic Materials

In this study, a theoretical Insight into two newly designed novel CL-20-based high performance energetic compounds, namely bis(nitromethyl)-tetranitrohexaaza-isowurtzitane (BNMTNIW) and bis(nitratomethyl)-tetranitrohexaaza-isowurtzitane BNIMTNIW), is reported. The title compounds are expected to...

Analysis of Ignition Sites for the Explosives 3,3'-Diamino-4,4'-azoxyfurazan (DAAF) and 1,3,5,7-Tetranitro-1,3,5,7-tetrazoctane (HMX) Using Crush Gun Impact Testing

The handling safety characteristics of energetic materials must be measured in order to ensure the safe transport and use of explosives. Drop-weight impact sensitivity measurements are one of the first standardized tests performed for energetics. They utilize a small amount of the explosive sample and a standard weight, which is dropped on the material from various heights to determine its sensitivity. While multiple laboratories have used the impact sensitivity test as an initial screening tool for explosive sensitivity for the past 60 years, variability exists due to the use of different instruments, different methods to determine the initiation, and the scatter commonly associated with less-sensitive explosives. For example, standard explosives such as 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX) initiate reliably and consistently on the drop-weight impact test, whereas insensitive explosives such as 3,3'-diamino-4,4'-azoxyfurazan (DAAF) exhibit variability in sound levels and the expended material. Herein we investigate the impact sensitivity of DAAF and HMX along with a more detailed investigation of ignition sites using a novel "crush gun" apparatus: a pneumatically powered drop-weight tower with advanced diagnostics, including high-speed visual and infrared cameras. Using this crush gun assembly, the ignition sites in HMX and DAAF were analyzed with respect to the effects of particle size and the presence of a source of grit. The formation of ignition sites was observed in both explosives; however, only HMX showed ignition sites that propagated to a deflagration at lower firing speeds. Finally, the presence of grit particles was shown to increase the occurrence of ignition sites in DAAF at lower firing speeds, though propagation to a full reaction was not observed on the time scale of the test. These results enable a better understanding of how ignition and propagation occurs during the impact testing of DAAF.

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005-1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2-3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. These backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si-O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.