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

Transport Properties of Liquid Pentaerythritol Tetranitrate (PETN) from Molecular Dynamics Simulations

We have used molecular dynamics simulations to determine the transport properties of liquid pentaerythritol tetranitrate (PETN), an important energetic material. The density, ρ, self-diffusion coefficient, D, thermal conductivity, κ, and shear viscosity, μ, have been computed over pressures and temperatures relevant to the subshock regime (up to 1000 K and a few GPa), where PETN is known to melt prior to initiation. We find that the thermal conductivity κ(P, T) can be represented by a simple analytical function that fits the data points with very good accuracy, even beyond the subshock regime, up to 2000 K and 20 GPa. The self-diffusion coefficient, D, exhibits nonmonotonic behavior, with notably the temperature-independent prefactor decreasing by several orders of magnitude between 0 and 2 GPa before remaining nearly constant after, and the activation energy varying little in the subshock regime before increasing linearly beyond. Lastly, the viscosity, μ, is well described by Nahme's law, which is fitted to the MD results and allows us to predict μ(P, T) for temperatures and pressures corresponding to the subshock regime. These results can be used to model the response of PETN to low-velocity impacts, where the material melts prior to the first reactions, and thermal conduction and viscosity play a crucial role.

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.

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.

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.

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.

Effects of Low-Level Gamma Radiation on Common Nitroaromatic, Nitramine, and Nitrate Ester Explosives

The aging of high explosives in an ionizing radiation field is not well understood, and little work has been done in the low dose and low dose rate regime. In this study, four explosives were exposed to low-level gamma irradiation from a 137Cs source: PETN, PATO, and PBX 9501 both with and without the Irganox 1010 stabilizer. Post-irradiation analysis included GC-MS of the headspace gas, SEM of the pellets and powder, NMR spectroscopy, DSC analysis, impact sensitivity tests, and ESD sensitivity tests. Overall, no significant change to the materials was seen for the dose and dose rate explored in this study. A small change in the 1H NMR spectrum of PETN was observed and SEM and ESD results suggest a surface energy change in PATO, but these differences are minor and do not appear to have a substantial impact on the handling safety.