Analytical Techniques in Molecular Simulation and Its Application in Energetic Materials

Efficient design on the molecular and crystal levels is an urgent need to accelerate the development of energetic materials (EMs), and the performance analysis of microstructures is the most important thing in the research and design of EMs. Although molecular simulation methods are widely used in various research fields, there are few comprehensive reviews on analytical techniques. It is urgent to understand the basic principles of various analytical methods in the research of EMs. In this article, the characterization/analysis methods in quantum mechanics and molecular mechanics simulation technology are summarized, and their applications in the field of EMs are listed. At the molecular level, energy, composition, geometric structure, and electronic structure are all related to macroscopic properties, and most of them have been widely used as variables in numerical models to predict and compare the properties of EMs. In addition, this paper emphasizes that the correlation between theoretical calculation and confirmatory experiment needs to be further verified because the current experimental characterization can also be dealt with by molecular simulation, which is helpful to the popularization and application of theoretical research.

Theoretical mechanistic insights on the thermal and acid-catalyzed rearrangements of N-methyl-N-nitroanilines

The thermal and acid-catalyzed rearrangement mechanisms of N-methyl-N-nitroanilines were theoretically investigated via density functional theory (DFT) calculations for all possible proposed mechanisms. The results indicate that the thermal rearrangement of N-methyl-N-nitroanilines undergoes a radical pair complex mechanism through the homolysis of their N-N bond to generate a radical pair complex and the recombination of the radical pairs followed by aromatization. For the acid-catalyzed rearrangements, N-methyl-N-nitroanilines are first protonated on the nitrogen atom of their aniline moiety and then generate protonated N-methyl-O-nitroso-N-phenylhydroxylamines through a three-membered spirocyclic oxadiaziridine transition state. The N-protonated N-methyl-O-nitroso-N-phenylhydroxylamines favor homolytic dissociation to generate N-methylaniline cationic radical and nitrogen dioxide complexes, which further combine together and aromatize to afford protonated N-methyl-o-nitroanilines and N-methyl-p-nitroanilines, respectively. The radical pair complexes are more stable than the corresponding solvent-caged radical pairs. The thermal rearrangements require higher activation energy than the corresponding acid-catalyzed rearrangements.

Revealing the Ultrafast Energy Transfer Pathways in Energetic Materials: Time-Dependent and Quantum State-Resolved

Intramolecular vibrational energy transfer is gaining tremendous attention as a regulator of condensed-phase behavior and reactions. In polyatomic molecules, this transfer is an ultrafast process involving multiple modes with numerous quantum states. The inherent complexity and rapid evolution of these processes pose significant challenges to experimental observation, and the high computational costs make full quantum mechanical calculations impractical with current technology. In the intramolecular energy transfer process, whether the doorway modes are intermediaries for transferring energy from lattice phonons to high-frequency intramolecular vibrational modes has been a controversial issue. However, the broad range of doorway modes complicates the experimental identification of a specific doorway in the transfer process corresponding to a specific end point. Here, for the first time, we utilize a combination of vibrational projection, statistical analysis, and the local quantum vibrational embedding (LQVE) method to elucidate the ultrafast energy transfer pathways that upconvert energy from lattice phonons to intramolecular modes in the typical energetic material β-HMX. This approach enables us to resolve the coupled vibrational mode groups, identify the most probable energy transfer pathways corresponding to the different final modes, and clearly confirm that the doorway region is a mandatory pathway for energy transfer. The LQVE method's time-dependent and quantum state-resolved advantages are leveraged to reveal the microscopic mechanism of the energy transfer process. The time scale of these processes is determined at about 1 ps, and the first theoretical two-dimensional infrared spectroscopy evidence is provided, which is confirmed by the experimental results. These findings deliver important insights into the fundamental mechanisms of ultrafast energy transfer in energetic materials, providing theoretical support for controlling explosive behavior and designing new explosives. The methodologies developed in this work can be extended to other condensed phase materials and used to evaluate the coupling between multiple vibrational modes.

Theoretical advances in understanding and enhancing the thermostability of energetic materials

The quest for thermally stable energetic materials is pivotal in advancing the safety of applications ranging from munitions to aerospace. This perspective delves into the role of theoretical methodologies in interpreting and advancing the thermal stability of energetic materials. Quantum chemical calculations offer an in-depth understanding of the molecular and electronic structure properties of energetic compounds related to thermal stability. It is also essential to incorporate the surrounding interactions and their impact on molecular stability. Ab initio molecular dynamics (AIMD) simulations provide detailed theoretical insights into the reaction pathways and the key intermediates during thermal decomposition in the condensed phase. Analyzing the kinetic barrier of rate-determining steps under various temperature and pressure conditions allows for a comprehensive assessment of thermal stability. Recent advances in machine learning have demonstrated their utility in constructing potential energy surfaces and predicting thermal stability for newly designed energetic materials. The machine learning-assisted high-throughput virtual screening (HTVS) methodology can accelerate the discovery of novel energetic materials with improved properties. As a result, the newly identified and synthesized energetic molecule ICM-104 revealed excellence in performance and thermostability. Theoretical approaches are pivotal in elucidating the mechanisms underlying thermal stability, enabling the prediction and design of enhanced thermal stability for emerging EMs. These insights are instrumental in accelerating the development of novel energetic materials that optimally balance performance and thermal stability.

Interpretable and Physicochemical-Intuitive Deep Learning Approach for the Design of Thermal Resistance of Energetic Compounds

Thermal resistance of energetic materials is critical due to its impact on safety and sustainability. However, developing predictive models remains challenging because of data scarcity and limited insights into quantitative structure-property relationships. In this work, a deep learning framework, named EM-thermo, was proposed to address these challenges. A data set comprising 5029 CHNO compounds, including 976 energetic compounds, was constructed to facilitate this study. EM-thermo employs molecular graphs and direct message-passing neural networks to capture structural features and predict thermal resistance. Using transfer learning, the model achieves an accuracy of approximately 97% for predicting the thermal-resistance property (decomposition temperatures above 573.15 K) in energetic compounds. The involvement of molecular descriptors improved model prediction. These findings suggest that EM-thermo is effective for correlating thermal resistance from the atom and covalent bond level, offering a promising tool for advancing molecular design and discovery in the field of energetic compounds.

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.

Predicting the thermal decomposition temperature of energetic materials from a simple model

CONTEXT

The key factor in designing heat-resistant energetic materials is their thermal sensitivity. Further research and prediction of thermal sensitivity remains a great challenge for us. This study is based on first-principles calculations and establishes a theoretical model, which comprehensively considers band gap, density of states, and Young's modulus to obtain a empirical parameter Ψ. A quantitative relationship was established between the new parameter and the thermal decomposition temperature. The value of Ψ is calculated for 10 energetic materials and is found to have a strong correlation with the experimental thermal decomposition temperature. This further proves the reliability of our model. Specifically, the larger the value of Ψ, the higher the thermal decomposition temperature, and the more stable the energetic material will be. Therefore, to some extent, we can use the new parameter Ψ calculated by the model to predict thermal sensitivity.

METHODS

Based on first-principles, this paper used the Cambridge Serial Total Energy Package (CASTEP) module of Materials Studio (MS) for calculations. The Perdew-Burke-Ernzerhof (PBE) functionals in Generalized Gradient Approximation (GGA) method as well as the Grimme dispersion correction was used in this paper.

High energy barrier hydroxyl radical dissociation mechanism of a low shock sensitivity dihydroxylammonium 5,5'-bistetrazole-1,1'-diolate (TKX-50) explosive

As a representative of the new generation of high-energy explosives, TKX-50 has attracted widespread attention due to its remarkably low sensitivity toward shock. However, the reported decomposition barriers of TKX-50 (∼37 kcal mol-1) are comparable to those of commonly used explosives. The mechanism of its low shock sensitivity remains unclear. In this study, using an ab initio molecular dynamics method combined with a multiscale shock simulation technique and transition state calculations (at the B2PLYP-D3/Def2TZVP level), we discovered an unconventional reaction pathway of TKX-50 under shock, and its rate-controlling step is the dissociation of the hydroxyl radical (OH) from the anion ring after proton transfer, followed by ring rupture and the production of H2O and N2. The barrier for this OH dissociation reaction is as high as 51.9 kcal mol-1. In contrast, under thermal stimuli, TKX-50 prefers to open rings directly after proton transfer without losing the OH. The corresponding barrier is 35.4 kcal mol-1, which is in good agreement with previous studies. The reason for the unconventional reaction pathway of TKX-50 under shock may be the suppression of anion ring opening in thermal decomposition by steric hindrance upon shock compression. In addition, the dominant N2 generation pathway under shock releases less energy than pyrolysis which further explains the low shock sensitivity of TKX-50. This study comprehensively elucidates the different reaction mechanisms of TKX-50 under thermal and shock conditions and proposes a crucial reaction pathway leading to its low shock sensitivity. These findings will contribute to the understanding and application of tetrazole anionic energetic salts.

Molecular dynamics simulation of sensitivity of HMX, FOX-7, and TATB crystals

METHODS

This study used molecular dynamics (MD) to simulate three materials (HMX, FOX-7, and TATB) under the NVT ensemble. Six temperatures (100 K, 200 K, 300 K, 400 K, 500 K, and 600 K) were simulated. In addition, the trigger bond lengths, energy bands, and density of states of three materials were obtained at different temperatures and compared with the calculated results at 0 K.

CONTEXT

The results indicate that the trigger bond lengths of the three materials are very close to the experimental values. Overall, the maximum and average bond lengths of the trigger bonds increase with increasing temperature. The band gap value decreases with increasing temperature. The changes in trigger bond length and band gap value are consistent with the experimental fact that sensitivity increases with increasing temperature. And Eg > 1 eV is consistently found within the temperature range of 0-600 K, indicating that all three materials are non-metallic compounds.

Exploring the thermal decomposition and detonation mechanisms of 2,4-dinitroanisole by TG-FTIR-MS and molecular simulations

2,4-dinitroanisole (DNAN), an insensitive explosive, has replaced trinitrotoluene (TNT) in many melt-cast explosives to improve the safety of ammunition and becomes a promising material to desensitize novel explosives of high sensitivity. Here, we combine thermogravimetric-Fourier transform infrared spectrometry-Mass spectrometry (TG-FTIR-MS), density functional theory (DFT), and ReaxFF molecular dynamics (MD) to investigate its thermal decomposition and detonation mechanisms. As revealed by TG-FTIR-MS, the thermal decomposition of DNAN starts at ca. 453 K when highly active NO2 is produced and quickly converted to NO resulting in the formation of a large amount of Ph(OH)(OH2)OCH3+. DFT calculations show that the activation energy of DNAN is higher than that of TNT due to the lack of α-H. Further steps in both thermal decomposition and detonation reactions of the DNAN are dominated by bimolecular O-transfers. ReaxFF MD indicates that DNAN has a lower heat of explosion than TNT, in accordance with the observation that the activation energies of polynitroaromatic explosives are inversely proportional to their heat of explosion. The inactive -OCH3 group and less nitro groups also render DNAN higher thermal stability than TNT.

Bistriazolotriazole-tetramine: commendable energetic moiety and cation

CONTEXT

Various 7H,7'H-[6,6'-bi[1,2,4]triazolo[4,3-b][1,2,4]triazole]-3,3',7,7'-tetramine (A) based nitrogen-rich energetic salts were designed and their properties explored. All energetic salts possess relatively high nitrogen content (> 48%), positive heats of formation (> 429 kJ/mol) and stability owing to a significant contribution from fused backbone. The cationic component shows a very high heat of formation (2516 kJ/mol); therefore, it is highly suitable for enthalpy enhancement in new energetic salts. The cation was paired with the energetic anions nitrate (NO3-), perchlorate (ClO4-), dinitromethanide (CH(NO2)2-), trinitromethanide (C(NO2)3-), nitroamide (NHNO2-), and dinitroamide (N(NO2)2-) to improve oxygen balance and detonation performance. Designed salts show moderate detonation velocities (7.9-8.7 km/s) and pressures (23.8 - 33.1 GPa). The distribution of frontier molecular orbitals, molecular electrostatic surface potentials, QTAIM topological properties, and noncovalent interactions of designed salts were simulated to understand the electronic structures, charge distribution in molecules, hydrogen bonding, and other nonbond interactions. The predicted safety factor (SF) and impact sensitivity (H50) of designed salts suggest their insensitivity to mechanical stimuli. This work explored the 7H,7'H-[6,6'-bi[1,2,4]triazolo[4,3-b][1,2,4]triazole]-3,3',7,7'-tetramine as a suitable cationic component which could be promising and serve exemplarily in energetic materials.

METHODS

The optimization and energy calculations of all the designed compounds were carried out at the B3LYP/6-311 +  + G(d,p) and M06-2X/def2-TZVPP levels, utilizing the Gaussian software package. The molecular surface electrostatic potential, quantum theory of atoms in molecules (QTAIM), reduced density gradient (RDG), and noncovalent interaction (NCI) analysis were performed by employing Multiwfn software. The EXPLO5 (v 7.01) thermochemical code and PILEM web application were used to predict the detonation properties.

Probing into the theory of impact sensitivity: propelling the understanding of phonon-vibron coupling coefficients

The determination of impact sensitivity of energetic materials traditionally relies on expensive and safety-challenged experimental means. This has instigated a shift towards scientific computations to gain insights into and predict the impact response of energetic materials. In this study, we refine the phonon-vibron coupling coefficients ζ in energetic materials subjected to impact loading, building upon the foundation of the phonon up-pumping model. Considering the full range of interactions between high-order phonon overtones and molecular vibrational frequencies, this is a pivotal element for accurately determining phonon-vibron coupling coefficients ζ. This new coupling coefficient ζ relies exclusively on phonon and molecular vibrational frequencies within the range of 0-700 cm-1. Following a regression analysis involving ζ and impact sensitivity (H50) of 45 molecular nitroexplosives, we reassessed the numerical values of damping factors, establishing a = 2.5 and b = 35. This coefficient is found to be a secondary factor in determining sensitivity, secondary to the rate of decomposition propagation and thermodynamic factor (heat of explosion). Furthermore, the relationship between phonon-vibron coupling coefficients ζ and impact sensitivity was studied in 16 energetic crystalline materials and eight nitrogen-rich energetic salts. It was observed that as the phonon-vibron coupling coefficient increases, the tendency for reduced impact sensitivity H50 still exists.

A physical organic strategy to predict and interpret stabilities of chemical bonds in energetic compounds for the discovery of thermal-resistant properties

CONTEXT

The in-depth understanding about the stability of chemical bonds in energetic compounds plays a central role for molecular design and safety-related evaluations. Most energetic compounds contain nitro as explosophores, and nitro cleavage is fundamental for thermal and mechanical stability. However, the quantum chemistry approach to accurately predict energy and temperature properties related to bond stability is challenging, due to the tradeoff between computational costs and deviations. Herein, the bond orders are proposed as accurate and computational-cost efficient descriptors for predicting the chemical bond stability and thermal-resistant properties. The intrinsic bond strength index (IBSI) demonstrates the best prediction for experimental homolytic bond dissociation energies (R2 > 0.996), which is on par with the results from high-precision quantum chemistry methods. The effects from bond connectivity and steric hindrance hierarchy were analyzed to reveal underlying mechanisms. Additionally, the IBSI descriptors are successfully applied to predict the thermal decomposition temperatures of 24 heat-resistant energetic compounds (R2 = 0.995), thus validating the effectiveness for the prediction and interpretation of chemical bond stability in energetic compounds via a physical organic approach.

METHODS

All DFT calculations were performed with Gaussian 09 software. To investigate the dependence of the method on functionals and basis sets, 9 DFT methods were considered (B3LYP/6-31G(d,p), B3LYP/6-311G(d,p), B3LYP/def2-TZVP, M062X/6-31G(d,p), M062X/6-311G(d,p), M062X/def2-TZVP, ωB97XD/6-31G(d,p), ωB97XD/6-311G(d,p), and ωB97XD/def2-TZVP). The bond order descriptors LBO and IBSI are obtained through the bond order analysis module in the Multiwfn software.

Azo-Bridged Triazole Macrocycles: Computational Design, Energy Content, Performance, and Stability Assessment

Oxadiazole and triazole are extensively investigated heterocyclic scaffolds in the development of energetic materials. New energetic molecules were designed by replacing 1,2,5-oxadiazole with 2H-1,2,3-triazole in the reported conjugated macrocyclic systems to assess the influence on the energetic properties and stability. In addition, nitro groups were introduced in triazole units (N-functionalization) to improve the energetic performance. Energetic properties, including heat of formation, oxygen balance, density, detonation pressure and velocity, and impact sensitivity, were estimated for these triazole-based macrocycles. The replacement of 1,2,5-oxadiazole with 2H-1,2,3-triazole and 2-nitro-1,2,3-triazole significantly enhances the energy content, detonation performance, and noncovalent interactions. The theoretically computed energetic properties of triazole-based macrocycles reveal high positive heats of formation (1507-2761 kJ/mol), oxygen balance (-88.8 to -22.8%), high densities (1.87-1.90 g/cm3), superior detonation velocities (8.41-9.52 km/s), pressures (26.64-40.55 GPa), acceptable impact sensitivity (27-40 cm), and safety factor (51-290). The overall energetic assessment highlights triazole-based macrocycles as a potential framework that will be useful for developing advanced energetic materials.

Theoretical study on the correlation between the structure, excess energy, surface energy, electronic structure, nitro charge, and friction sensitivity of N,N'-dinitroethylenediamine (EDNA)

The sensitivity of energetic materials along different crystal directions is not the same and is anisotropic. In order to explore the difference in friction sensitivity of different surfaces, we calculated the structure, excess energy, surface energy, electronic structure, and the nitro group along (1 1 1), (1 1 0), (1 0 1), (0 1 1), (0 0 1), (0 1 0), and (1 0 0) surfaces of EDNA based on density functional theory. The analysis results showed that relative to other surfaces, the (0 0 1) surface has the shortest N-N average bond length, largest N-N average bond population, smallest excess energy and surface energy, widest band gap, and the largest nitro group charge value, which indicates that the (0 0 1) surface has the lowest friction sensitivity compared to other surfaces. Furthermore, the conclusions obtained by analyzing the excess energy are consistent with the results of the N-N bond length and bond population, band gap, and nitro charge. Therefore, we conclude that the friction sensitivity of different surfaces of EDNA can be evaluated using excess energy.

Optimization of performance and sensitivity: preparation of two Ag(I)-based ECPs by using isomeric ligands

For energetic compounds, their structure determines their performance, and even minor variations in their structure can have a significant impact on their performance. The application scenarios for energetic materials are diverse, and their performance requirements vary as well. To investigate the influence of different substituent positions on the performance of primary explosives, we prepared two Ag(I)-based complexes, [Ag(2-IZCA)ClO4]n (ECPs-1) and [Ag(4-IZCA)ClO4]n (ECPs-2), using structurally isomeric ligands, 1H-imidazole-2-carbohydrazide (2-IZCA) and 1H-imidazole-4-carbohydrazide (4-IZCA). The structures were confirmed using infrared, elemental analysis, and single-crystal X-ray diffraction. Experimental results demonstrate that both ECPs exhibit good thermal stability. However, compared to ECPs-1, ECPs-2 exhibits a lower thermal initial decomposition temperature (Td = 210 °C), lower mechanical sensitivity (IS = 27 J, FS = 84 N), and more concentrated energy output. Although theoretical predictions suggest similar detonation velocities and pressures for both compounds, actual detonation performance tests indicate that ECPs-2 has stronger explosive power and initiating capability, with potential for use as a laser initiator (E = 126 mJ). The simple preparation method and inexpensive starting materials enrich the research on primary explosives.

Fluorine-Containing Functional Group-Based Energetic Materials

Molecules featuring fluorine-containing functional groups exhibit outstanding properties with high density, low sensitivity, excellent thermal stability, and good energetic performance due to the strong electron-withdrawing ability and high density of fluorine. Hence, they play a pivotal role in the field of energetic materials. In light of current theoretical and experimental reports, this review systematically focuses on three types of energetic materials possessing fluorine-containing functional groups F- and NF2 - substituted trinitromethyl groups (C(NO2 )2 F, C(NO2 )2 NF2 ), trifluoromethyl group (CF3 ), and difluoroamino and pentafluorosulfone groups (NF2 , SF5 ) and investigates the synthetic methods, physicochemical parameters, and energetic properties of each. The incorporation of fluorine-containing functional moieties is critical for the development of novel high energy density materials, and is rapidly being adopted in the design of energetic materials.

Halogenated PETN derivatives: interplay between physical and chemical factors in explosive sensitivity

Determining the factors that influence and can help predict energetic material sensitivity has long been a challenge in the explosives community. Decades of literature reports identify a multitude of factors both chemical and physical that influence explosive sensitivity; however no unifying theory has been observed. Recent work by our team has demonstrated that the kinetics of "trigger linkages" (i.e., the weakest bonds in the energetic material) showed strong correlations with experimental drop hammer impact sensitivity. These correlations suggest that the simple kinetics of the first bonds to break are good indicators for the reactivity observed in simple handling sensitivity tests. Herein we report the synthesis of derivatives of the explosive pentaerythritol tetranitrate (PETN) in which one, two or three of the nitrate ester functional groups are substituted with an inert group. Experimental and computational studies show that explosive sensitivity correlates well with Q (heat of explosion), due to the change in the number of trigger linkages removed from the starting material. In addition, this correlation appears more significant than other observed chemical or physical effects imparted on the material by different inert functional groups, such as heat of formation, heat of explosion, heat capacity, oxygen balance, and the crystal structure of the material.

Synthesis, Characterization, and Properties of High-Energy Fillers Derived from Nitroisobutylglycerol

Herein we report a comprehensive laboratory synthesis of a series of energetic azidonitrate derivatives (ANDP, SMX, AMDNNM, NIBTN, NPN, 2-nitro-1,3-dinitro-oxypropane) starting from the readily available nitroisobutylglycerol. This simple protocol allows obtaining the high-energy additives from the available precursor in yields higher than those reported using safe and simple operations not presented in previous works. A detailed characterization of the physical, chemical, and energetic properties including impact sensitivity and thermal behavior of these species was performed for the systematic evaluation and comparison of the corresponding class of energetic compounds.

A Physical Organic Approach towards Statistical Modeling of Tetrazole and Azide Decomposition

Nitrogen atom-rich heterocycles and organic azides have found extensive use in many sectors of modern chemistry from drug discovery to energetic materials. The prediction and understanding of their energetic properties are thus key to the safe and effective application of these compounds. In this work, we disclose the use of multivariate linear regression modeling for the prediction of the decomposition temperature and impact sensitivity of structurally diverse tetrazoles and organic azides. We report a data-driven approach for property prediction featuring a collection of quantum mechanical parameters and computational workflows. The statistical models reported herein carry predictive accuracy as well as chemical interpretability. Model validation was successfully accomplished via tetrazole test sets with parameters generated exclusively in silico. Mechanistic analysis of the statistical models indicated distinct divergent pathways of thermal and impact-initiated decomposition.

Which molecular properties determine the impact sensitivity of an explosive? A machine learning quantitative investigation of nitroaromatic explosives

We decomposed density functional theory charge densities of 53 nitroaromatic molecules into atom-centered electric multipoles using the distributed multipole analysis that provides a detailed picture of the molecular electronic structure. Three electric multipoles, (the charge of the nitro groups), (the total dipole, i.e., polarization, of the nitro groups), (the total electron delocalization of the C ring atoms), and the number of explosophore groups (#NO₂) were selected as features for a comprehensive machine learning (ML) investigation. The target property was the impact sensitivity h₅₀ (cm) values quantified by drop-weight measurements, with a large h₅₀ (e.g., 150 cm) indicating that an explosive is insensitive and vice versa. After a preliminary screening of 42 ML algorithms, four were selected based on the lowest root mean square errors: Extra Trees, Random Forests, Gradient Boosting, and AdaBoost. Compared to experimental data, the predicted h₅₀ values of molecules having very different sensitivities for the four algorithms have differences in the range 19-28%. The most important properties for predicting h₅₀ are the electron delocalization in the ring atoms and the polarization of the nitro groups with averaged weights of 39% and 35%, followed by the charge (16%) and number (10%) of nitro groups. A significant result is how the contribution of these properties to h₅₀ depends on their actual sensitivities: for the most sensitive explosives (h₅₀ up to ∼50 cm), the four properties contribute to reducing h₅₀, and for intermediate ones (∼50 cm ≲ h₅₀ ≲ 100 cm) #NO₂ and contribute to increasing it and the other two properties to reducing it. For highly insensitive explosives (h₅₀ ≳ 200 cm), all four properties essentially contribute to increasing it. These results furnish a consistent molecular basis of the sensitivities of known explosives that also can be used for developing safer new ones.

Heuristic-Based Alkaline Hydrolysis Mechanism of Nitrate Ester (Nitrocellulose Monomer) and Nitroamine (Hexogen) Compounds: Electrostatic Attraction Effect of the Nitro Group

The alkaline hydrolysis reaction of energetic materials is important and complex. With improved performance, AMK_Mountain was used to systematically study the alkaline hydrolysis of the nitrocellulose monomer and hexogen. The reaction pathways showed that the nitrocellulose monomer produces the nitrate anion and nitrite anion differently, while hexogen only produces the nitrite anion. Electronic structure results at the M06-2X/6-311G(d,p)/PCM(Pauling) level showed that the nitrocellulose monomer and hexogen have a similar pathway in their main energy-releasing process (nitrite anion production): with electrostatic attraction effects after proton transfer, the nitrite anion dissociates from the original structure with a low barrier. Moreover, during the alkaline hydrolysis of the nitrocellulose monomer, the metastable intermediates after proton transfer may be directly generated following transition states that, structurally, tend to produce nitrite anions "proximal" to the proton transfer site and produce nitrate anions "distal" to the proton transfer site. Electronic structure analysis showed that representative metastable intermediates revealed that the charge transfer caused by electrostatic attraction may be the direct cause of these reactions.

Triggering the mechanism of the initial reaction of energetic materials under pressure based on Raman intensity analysis

The Raman intensity and other stoichiometric calculations of nitromethane (NM) and 2-nitrimino-5-nitro-hexahydro-1,3,5-triazine (NNHT) have been made by using first-principles density functional theory. We propose a method to judge the initial reaction mechanism of NM and NNHT under pressure based on the Raman intensity. Both the resulting NM and NNHT undergo hydrogen transfer and conventional trigger bond cleavage. And the results obtained from the Raman peak intensities infer a reaction path that is not inferior to the traditional C-NO₂ and N-NO₂ bond cleavage, thus verifying our results.

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.

Effects of pressure on structural, electronic, optical, and mechanical properties of nitrogen-rich energetic material: 6-azido-8-nitrotetrazolo[1,5-b]pyridazine-7-amine (3at)

CONTEXT AND RESULTS

6-Azido-8-nitrotetrazolo[1,5-b]pyridazine-7-amine (3at) is a promising green energetic material, which meets the development requirements of environment-friendly explosives. By discussing the relationship between lattice parameters and pressure, it is found that the compression ratio indicates anisotropy of compressibility. And bond lengths get shorter under pressure, resulting in stronger intermolecular bonds. The N₃ group rotates under pressure. And then, the optical properties basically change regularly with the change of pressure. As the pressure increases, the absorption range widens. In the low energy interval, it shows transparency, and then with the increase of energy and pressure, it shows better optical activity. With the increase of pressure and energy, the absorption coefficient increases, representing that the optical activity becomes high. Finally, according to the analysis of mechanical properties, 3at exhibited brittle behavior at 0 GPa and 100 GPa, while at 10 to 90 GPa, the values of ν and B/G are malleable.

COMPUTATIONAL AND THEORETICAL TECHNIQUES

Based on density functional theory, the crystal parameters, electronic properties, optical properties, and elastic and mechanical properties of 3at under different pressures were studied theoretically. The GGA-PW91+OBS method was used to calculate the physical parameters under pressure, such as lattice parameters, energy band structures, dielectric function, refractive index, absorption coefficient, and elastic constants. Physical properties under (3at) pressure are predicted.

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.

Spun of improvised cis-1,3,4,6-tetranitrooctahydroimidazo-[4,5-D]-Imidazole (BCHMX) in polystyrene nanofibrous membrane by electrospinning techniques

Development of ultra-fine fiber technology and nano-sized materials are widely taking place to enhance the characteristic of different materials. In our study, a newly developed technique was used to produce improvised nano energetic fibers with the exploitation of cis-1,3,4,6-Tetranitrooctahydroimidazo-[4,5-d] imidazole (BCHMX) to spin in a polystyrene nanofiber membrane. Scanning electron microscopy (SEM) showed the synthesized nanofibrous polystyrene (PS)/BCHMX sheets with clear and continual fiber were imaged with scanning electron microscopy (SEM). Characterization of the produced nanofiber was examined by Fourier Transform Infrared (FTIR), and X-ray diffractometer (XRD). Explosive sensitivity was also evaluated by both BAM impact and friction apparatus. Thermal behavior for the synthesized PS/BCHMX fiber and the pure materials were also investigated by thermal gravimetric analysis (TGA). The results show enhancement in the fabrication of nano energetic fibers with a size of 200-460 nm. The TG confirms the high weight percentage of BCHMX which reaches 60% of the total mass. PS/BCHMX fiber was confirmed with the XRD, FTIR spectrum. Interestingly, XRD sharp peaks showed the conversion of amorphous PS via electrospinning into crystalline shape regarding the applied high voltage. The synthesized PS/BCHMX nanofiber was considered insensitive to the mechanical external stimuli; more than 100 J impact energy and  > 360 N initiation force as friction stimuli. PS/BCHMX is considering a candidate tool to deal with highly sensitive explosives safely and securely for explosives detection training purposes.

Unimolecular Decomposition Reactions of Picric Acid and Its Methylated Derivatives─A DFT Study

To handle energetic materials safely, it is important to have knowledge about their sensitivity. Density functional theory (DFT) has proven a valuable tool in the study of energetic materials, and in the current work, DFT is employed to study the thermal unimolecular decomposition of 2,4,6-trinitrophenol (picric acid, PA), 3-methyl-2,4,6-trinitrophenol (methyl picric acid, mPA), and 3,5-dimethyl-2,4,6-trinitrophenol (dimethyl picric acid, dmPA). These compounds have similar molecular structures, but according to the literature, mPA is far less sensitive to impact than the other two compounds. Three pathways believed important for the initiation reactions are investigated at 0 and 298.15 K. We compare the computed energetics of the reaction pathways with the objective of rationalizing the unexpected sensitivity behavior. Our results reveal a few if any significant differences in the energetics of the three molecules, and thus do not reflect the sensitivity deviations observed in experiments. These findings point toward the potential importance of crystal structure, crystal morphology, bimolecular reactions, or combinations thereof on the impact sensitivity of nitroaromatics.

Co-agglomerated crystals of cyclic nitramines with sterically crowded molecules

Co-crystallization of cyclic nitramines with sterically crowded molecules, i.e., ε-2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (ε-CL-20) with cis-1,3,4,6-tetranitrooctahydroimidazo-[4,5-d]imidazole, is difficult, but co-agglomeration given good results.

Synthesis of Energetic 7-Nitro-3,5-dihydro-4H-pyrazolo[4,3-d][1,2,3]triazin-4-one Based on a Novel Hofmann-Type Rearrangement

Rearrangement reactions are efficient strategies in organic synthesis and contribute enormously to the development of energetic materials. Here, we report on the preparation of a fused energetic structure of 7-nitro-3,5-dihydro-4H-pyrazolo[4,3-d][1,2,3]triazin-4-one (NPTO) based on a novel Hofmann-type rearrangement. The 1,2,3-triazine unit was introduced into the fused bicyclic skeleton from a pyrazole unit for the first time. The new compound of NPTO was fully characterized using multinuclear NMR and IR spectroscopy, elemental analysis as well as X-ray diffraction studies. The thermal behaviors and detonation properties of NPTO were investigated through a differential scanning calorimetry (DSC-TG) approach and EXPLO5 program-based calculations, respectively. The calculation results showed similar detonation performances between NPTO and the energetic materials of DNPP and ANPP, indicating that NPTO has a good application perspective in insensitive explosives and propellants.

Nature of the Trigger Linkage in Explosive Materials Is a Charge-Shift Bond

Explosion begins by rupture of a specific bond, in the explosive, called a trigger linkage. We characterize this bond in nitro-containing explosives. Valence-bond (VB) investigations of X-NO₂ linkages in alkyl nitrates, nitramines, and nitro esters establish the existence of Pauli repulsion that destabilizes the covalent structure of these bonds. The trigger linkages are mainly stabilized by covalent-ionic resonance and are therefore charge-shift bonds (CSBs). The source of Pauli repulsion in nitro explosives is unique. It is traced to the hyperconjugative interaction from the oxygen lone pairs of NO₂ into the σ_(X-N)^* orbital, which thereby weakens the X-NO₂ bond, and depletes its electron density as X becomes more electronegative. Weaker trigger bonds have higher CSB characters. In turn, weaker bonds increase the sensitivity of the explosive to impacts/shocks which lead to detonation. Application of the analysis to realistic explosives supports the CSB character of their X-NO₂ bonds by independent criteria. Furthermore, other families of explosives also involve CSBs as trigger linkages (O-O, N-O, Cl-O, N-I, etc. bonds). In all of these, detonation is initiated selectively at the CSB of the molecule. A connection is made between the CSB bond-weakening and the surface-electrostatic potential diagnosis in the trigger bonds.

How aromatic system size affects the sensitivities of highly energetic molecules?

Positive values of electrostatic potentials above the central regions of the molecular surface are strongly related to the high sensitivities of highly energetic molecules. The influence of aromatic system size on the positive values of electrostatic potentials and bond dissociation energies of C-NO2 bonds was studied by Density Functional Theory (DFT) calculations on a series of polycyclic nitroaromatic molecules. Calculations performed at PBE/6-311G** level showed that with the increase of the aromatic system size, values of positive electrostatic potential above the central areas of selected energetic molecules decrease from 32.78 kcal mol-1 (1,2,4,5-tetranitrobenzene) to 15.28 kcal mol-1 (2,3,9,10-tetranitropentacene) leading to the decrease in the sensitivities of these molecules towards detonation. Results of the analysis of electrostatic potential maps were in agreement with the trends in bond dissociation energies calculated for C-NO2 bonds of studied nitroaromatic molecules. Bond dissociation energies values indicate that the C-NO2 bond in the molecule of 1,2,4,5-tetranitrobenzene (56.72 kcal mol-1) is weaker compared to the nitroaromatic molecules with the additional condensed aromatic rings and with a similar arrangement of -NO2 groups (59.75 kcal mol-1 in the case of 2,3,9,10-tetranitropentacene). The influence of the mutual arrangement of -NO2 groups on the sensitivity of nitroaromatic molecules was also analyzed. Results obtained within this study could be of great importance for the development of new classes of highly energetic molecules with lower sensitivity towards detonation.

Fabrication of Energetic Composites with 91% Solid Content by 3D Direct Writing

Direct writing is a rapidly developing manufacturing technology that is convenient, adaptable and automated. It has been used in energetic composites to manufacture complex structures, improve product safety, and reduce waste. This work is aimed at probing the formability properties and combustion performances of aluminum/ammonium perchlorate with a high solid content for direct writing fabrication. Four kinds of samples with different solid content were successfully printed by adjusting printing parameters and inks formulas with excellent rheological behavior and combustion properties. A high solid content of 91% was manufactured and facile processed into complex structures. Micromorphology, rheology, density, burning rate, heat of combustion and combustion performance were evaluated to characterized four kinds of samples. As the solid content increases, the density, burning rate and heat of combustion are greatly enhanced. Based on 3D direct writing technology, complex energetic 3D structures with 91% solid content are shaped easier and more flexibly than in traditional manufacturing process, which provides a novel way for the manufacture of complicated structures of energetic components.

Fault Detection and Diagnosis for Plasticizing Process of Single-Base Gun Propellant Using Mutual Information Weighted MPCA under Limited Batch Samples Modelling

Aiming at the lack of reliable gradual fault detection and abnormal condition alarm and evaluation ability in the plasticizing process of single-base gun propellant, a fault detection and diagnosis method based on normalized mutual information weighted multiway principal component analysis (NMI-WMPCA) under limited batch samples modelling was proposed. In this method, the differences of coupling correlation among multi-dimensional process variables and the coupling characteristics of linear and nonlinear relationships in the process are considered. NMI-WMPCA utilizes the generalization ability of a multi-model to establish an accurate fault detection model in limited batch samples, and adopts fault diagnosis methods based on a multi-model SPE statistic contribution plot to identify the fault source. The experimental results demonstrate that the proposed method is effective, which can realize the rapid detection and diagnosis of multiple faults in the plasticizing process.

Low-Power Laser Ignition of an Antenna-Type Secondary Energetic Copper Complex: Synthesis, Characterization, Evaluation, and Ignition Mechanism Studies

In recent years, development of new energetic compounds and formulations, suitable for ignition with relatively low-power lasers, is a highly active and competitive field of research. The main goal of these efforts is focused on achieving and providing much safer solutions for various detonator and initiator systems. In this work, we prepared, characterized, and studied thermal and ignition properties of a new laser-ignitable compound, based on the 5,6-bis(ethylnitroamino)-N'2,N'3-dihydroxypyrazine-2,3-bis(carboximidamide) (DS3) proligand. This new energetic proligand was prepared in three steps, starting with 5,6-bis(ethylamino)-pyrazine-2,3-dicarbonitrile. Crystallography studies of the DS3-derived Cu(II) complex (DS4) revealed a unique stacked antenna-type structure of the latter compound. DS4 has an exothermal temperature of 154.5 °C and was calculated to exhibit a velocity of detonation of 6.36 km·s^-1 and a detonation pressure of 15.21 GPa. DS4 showed properties of a secondary explosive, having sensitivity to impact, friction, and electrostatic discharge of 8 J, 360 N, and 12 mJ, respectively. In order to study the mechanism of ignition by a laser (using a diode laser, 915 nm), we conducted a set of experiments that enabled us to characterize a photothermal ignition mechanism. Furthermore, we found that a single pulse, with a time duration of 1 ms and with a total energy of 4.6 mJ, was sufficient for achieving a consistent and full ignition of DS4. Dual-pulse experiments, with variable time intervals between the laser pulses, showed that DS4 undergoes ignition via a photothermal mechanism. Finally, calculating the chemical mechanism of the formation of the complex DS4 and modeling its anhydrous and hydrated crystal structures (density functional theory calculations using Gaussian and HASEM software) allowed us to pinpoint a more precise location of water molecules in experimental crystallographic data. These results suggest that DS4 has potential for further development to a higher technology readiness level and for integration into small-size safe detonator systems as for many civil, aerospace, and defense applications.