Rapid Materials Degradation Induced by Surfaces and Voids: Ab Initio Modeling of β-Octatetramethylene Tetranitramine

Initial Response of Pentaerythritol Tetranitrate (PETN) under the Coupling Effect of Preheating, Shock and Defect via the Molecular Dynamics Simulations with the Multiscale Shock Technique Method

The initial response of PETN under the coupling of preheating, impact and defects was simulated by Multiscale Shock Technique (MSST) method and molecular dynamics. The temperature change of PETN during impact compression can be divided into three stages: (1) the elastoplastic change of the system caused by initial compression; (2) part of PETN decomposes and releases energy to raise temperature; (3) a secondary chemical reaction occurs, resulting in rapid temperature rise. Under the given conditions, a higher initial preheating temperature will lead to faster decomposition of PETN; The existence of defects will accelerate the decomposition of PETN molecules; Coupling the highest preheating temperature with defects will lead to the fastest decomposition of PETN molecules, while in the defect-free PETN system with a preheating temperature of 300 K, the decomposition of PETN molecules is the slowest. For the case of U_s = 8 km·s^-1, the effect of defects on the initial PETN reaction is greater than the initial preheating temperature; When the impact velocity is greater than 9 km·s^-1, the impact velocity is an important factor affecting the decomposition of PETN molecules. For U_s = 10 km·s^-1, NO₂ is the main initial product in the defective PETN crystal, while in the perfect PETN crystal, it is the combination of NO₂ and HONO. The chemical reaction kinetics analysis shows that the preheating temperature and defects will accelerate the decomposition of PETN. The higher the preheating temperature, the faster the decomposition of PETN. For the case of U_s = 7 km·s^-1, 8 km·s^-1 and 9 km·s^-1, the existence of defects will increase the decomposition rate by more than 50% regardless of the initial preheating temperature. In the case of Us = 10 km·s^-1, the improvement of decomposition rate by defects is not as significant as the initial preheating temperature.

Molecular dynamic insight into octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) and the nano-HMX decomposition mechanism

The entire decomposition reaction process of a 30 Å HMX nanoparticle at 2000 K by ReaxFF molecular dynamics.

Review of the molecular and crystal correlations on sensitivities of energetic materials

Highly efficient design on the levels of molecule and crystal, as well as formulation, is highly desired for accelerating the development of energetic materials (EMs). Sensitivity is one of the most important characteristics of EMs and should be compulsorily considered in the design. However, owing to multiple factors responsible for the sensitivity, it usually undergoes a low predictability. Thus, it becomes urgent to clarify which factors govern the sensitivity and what is the importance of these factors. The present article focuses upon the progress of the molecular and crystal correlations on the sensitivity, and the molecule-based numerical models for sensitivity prediction in the past decades. On the molecular level, composition, geometric structure, electronic structure, energy and reactivity can be correlated with the sensitivity; while the sensitivity can be also related with molecular packing pattern, intermolecular interaction, crystal morphology, crystal size and distribution, crystal surface/interface and crystal defect on the crystal level. And most of these factors, in particle on the crystal level, have been employed as variables in numerical models for predicting sensitivity of categorized EMs. Besides, we stress that more attention should be paid to the sensitivity correlations on the inherent structures of EMs, molecule and crystal packing, because they can be readily dealt by molecular simulations nowadays, facilitating to reveal the physical nature of sensitivity.

A mechanism for two-step thermal decomposition of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a representative of the new generation of low-sensitivity energetic materials and has been applied extensively in formulations as an insensitive high-energetic ingredient. Although the initial thermal decomposition mechanism of LLM-105 has been studied based on quantum chemical calculations, the internal mechanism of the two-step thermal decomposition still lacks experimental research. Thus, this study involves a detailed experimental study to reveal the mechanism of the two-step thermal decomposition of LLM-105. The results showed that LLM-105 decay was a consecutive reaction. The first-step reaction dominated the early stage of the LLM-105 decomposition, and its products participated in the reaction of the second step. The cleavage of NO2 and NH2 groups of LLM-105 mainly occurred in the first step, while gaseous products NO and C2N2 were released during the second reaction step. The first-step reaction had a higher oxygen consumption rate and a lower carbon consumption rate, producing more heat due to more extensive oxidation of the carbon backbone. The difference in the oxidative ability and reaction rate between the two steps resulted in a two-step exothermic and mass loss behavior. This study provides further insights into the entire reaction process of LLM-105 and would be helpful for its better application and for the design of new explosives.

Decomposition mechanism scenarios of CL-20 co-crystals revealed by ReaxFF molecular dynamics: similarities and differences

Understanding the similarities and differences of decomposition mechanisms of CL-20 and its cocrystals is of great interest for practical applications of CL-20 cocrystals. The responses of CL-20 cocrystals to thermal stimulus were investigated using ReaxFF molecular dynamics simulations of two representative cocrystals, CL-20/HMX and CL-20/TNT, under adiabatic conditions and comparing to the baseline system of pure CL-20. The comprehensive chemical details were revealed with the aid of the unique code of VARxMD. The three CL-20-involved reactive systems all exhibit a distinct three-stage character during adiabatic decomposition when using the double peaks of the major intermediate NO₂ amount as the boundary. By taking advantage of the three-stage classification, a clear scenario for the similar stimulus-response of the CL-20 cocrystals can be elucidated for the dominant primary decomposition of CL-20 in stage I and the transition of favored chemical mechanisms from the generation of intermediates/radicals in stage II into their consumption to form stable products in stage III. The similar chemical behaviors are rooted in the dominance of CL-20 chemistry in the initial response of its cocrystals to thermal stimulus. The prolonged reaction zone uncovers the slowed decomposition kinetics of CL-20/HMX and CL-20/TNT, which is associated with the altered kinetics of CL-20 decomposition specifically by N-NO₂ bond scission and CL-20 skeleton decay. The retarded CL-20 decomposition in its cocrystals consequently results in more moderate self-heating and less violent exothermic reactions that agrees with the experimental observations of improved stability and damaged detonation performance of CL-20 cocrystals, particularly for CL-20/TNT. The results obtained in this work suggest that ReaxFF MD simulations can provide useful insight for the modulated chemical properties of varied CL-20 cocrystals.

A systematic study of the surface structures and energetics of CH3NO2 surfaces by first-principles calculations

Density functional theory (DFT) has been employed within the generalized gradient approximation and Perdew-Burke-Ernzerhof functional (GGA-PBE) to study the structural and electronic properties of nitromethane (NM) surface models. Different surfaces, including (100), (001), (101), (110), and (111), are considered in this work. The corresponding properties of bulk crystal for NM were also calculated to form a contrast to the slab models. Results with anisotropic characteristics of different surfaces have been observed in this study. There was an obviously great anisotropy in electronic parameters, especially the band gaps of different surfaces, indicating the anisotropic impact sensitivity along different directions of NM. The band gap value for (111) surface, 2.687 eV, was smaller than that of other surfaces, showing a higher impact sensitivity for NM. The estimated anisotropy has been revealed in surface energies for different surfaces. Graphical Abstract The valence band minimum (VBM) and conduction band maximum (CBM) of the nitromethane (100), (001), (101), (110) and (111) surface models.

The solid phase thermal decomposition and nanocrystal effect of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) via ReaxFF large-scale molecular dynamics simulation

The solid phase thermal decomposition and nanocrystal effect are extremely important to understand the ignition, combustion, reaction growth and buildup to detonation under shock wave action.

Adiabatic and constant volume decomposition process of condensed phase δ-1,3,5,7-tetranitro-1,3,5,7-tetrazocane at high temperatures: Quantum molecular dynamics simulations

We performed quantum molecular dynamics simulations to investigate the initiation chemistry of condensed phase δ-HMX at high temperatures by maintaining constant energy and volume to model adiabatic initiation process. The decomposition of HMX began by the C-N bond breaking in one molecule and by the C-H bond cleavage in other HMX molecule at 2400 K. At 2700 K, HMX is triggered by only one path that the C-N bond broke and the ring opened. At 3000 K, the decomposition of HMX is triggered by the C-H bond and N-O bond fission in the branch chains. There are seven decomposition channels observed during the whole decomposition stage. The N-O bond cleavage is a dominant reaction pathway. The boat configuration of the HMX molecule caused a new reaction channel to be happened by forming a new N-N bond. Another new reaction channel took place to form a new N-C bond due to intermolecular effects.

A computational study of ANTA and NTO derivatives

This work is a study of 5-amino-3-nitro-1,2,4-triazole (ANTA), 3-nitro-1,2,4-triazol-5-one (NTO), and nitrated derivatives of ANTA and NTO. RDX and TNT were studied for comparison. ANTA and NTO are low-sensitive high explosives with detonation properties comparable to 2,4,6-trinitrotoluene (TNT) and 1,3,5-trinitroperhydro-1,3,5-triazine (RDX). We showed previously that nitrated NTO and ANTA compounds, when used in a glycidyl azide polymer (GAP) matrix in rocket propellants, could give impulses above 2600 m/s and that the oxygen balance is positive. If used in aluminized explosives, the heat of detonation may be increased to a practical level significantly above RDX/aluminum compositions. Here, we use two different methods for sensitivity and two density functional theory functionals, B3LYP and M06-2X with the 6-31G(d) basis set, together with the complete basis set method CBS-4M. Calculations indicate that most of the nitrated derivatives have nearly equal sensitivity to RDX. Significantly different bond dissociation energies in the nitrimino functional group are predicted, although most models give much the same result.

Molecular Theory of Detonation Initiation: Insight from First Principles Modeling of the Decomposition Mechanisms of Organic Nitro Energetic Materials

This review presents a concept, which assumes that thermal decomposition processes play a major role in defining the sensitivity of organic energetic materials to detonation initiation. As a science and engineering community we are still far away from having a comprehensive molecular detonation initiation theory in a widely agreed upon form. However, recent advances in experimental and theoretical methods allow for a constructive and rigorous approach to design and test the theory or at least some of its fundamental building blocks. In this review, we analyzed a set of select experimental and theoretical articles, which were augmented by our own first principles modeling and simulations, to reveal new trends in energetic materials and to refine known existing correlations between their structures, properties, and functions. Our consideration is intentionally limited to the processes of thermally stimulated chemical reactions at the earliest stage of decomposition of molecules and materials containing defects.

Special catalytic effects of intermediate-water for rapid shock initiation of β-HMX

The intermediate-water efficiently promoted the decomposition of β-HMX, corresponding to a rapid shock initiation process.

Formation and growth mechanisms of natural metastable twin boundary in crystalline β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine: a computational study

The twin boundary (TB), a typical planar defect, occurs naturally in molecular explosives and manipulates their sensitivities to external stimuli.

Impact sensitivity and crystal lattice compressibility/free space

There is considerable evidence, which we discuss, indicating that compressibility and available free space in the crystal lattice are among the factors that govern the sensitivity of an explosive compound. Expanding and extending earlier work, we demonstrate, for 25 explosives, that there is an overall general tendency for greater impact sensitivity as the estimated free space per molecule increases. More specific relationships can be discerned by looking at subgroups of the compounds. The nitramine sensitivities, most of which are quite high, increase nearly linearly but only very gradually with free space. The nitroaromatics cover a wide range of sensitivities but all have an approximately similar intermediate level of free space. The remaining types of compounds show a reasonable sensitivity-free space relationship with one outlier: FOX-7 (1,1-diamino-2,2-dinitroethylene).

Decomposition mechanisms and kinetics of novel energetic molecules BNFF-1 and ANFF-1: quantum-chemical modeling

Decomposition mechanisms, activation barriers, Arrhenius parameters, and reaction kinetics of the novel explosive compounds, 3,4-bis(4-nitro-1,2,5-oxadiazol-3-yl)-1,2,5-oxadiazole (BNFF-1), and 3-(4-amino-1,2,5-oxadiazol-3-yl)-4-(4-nitro-1,2,5-oxadiazol-3-yl)-1,2,5-oxadiazole (ANFF-1) were explored by means of density functional theory with a range of functionals combined with variational transition state theory. BNFF-1 and ANFF-1 were recently suggested to be good candidates for insensitive high energy density materials. Our modeling reveals that the decomposition initiation in both BNFF-1 and ANFF-1 molecules is triggered by ring cleavage reactions while the further process is defined by a competition between two major pathways, the fast C-NO₂ homolysis and slow nitro-nitrite isomerization releasing NO. We discuss insights on design of new energetic materials with targeted properties gained from our modeling.

Molecular dynamics simulations of void defects in the energetic material HMX

A molecular dynamics (MD) simulation was carried out to characterize the dynamic evolution of void defects in crystalline octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine (HMX). Different models were constructed with the same concentration of vacancies (10 %) to discuss the size effects of void. Energetic ground state properties were determined by annealing simulations. The void formation energy per molecule removed was found to be 55-63 kcal/mol(-1), and the average binding energy per molecule was between 32 and 34 kcal/mol(-1) according to the change in void size. Voids with larger size had lower formation energy. Local binding energies for molecules directly on the void surface decreased greatly compared to those in defect-free lattice, and then gradually increased until the distance away from the void surface was around 10 Å. Analysis of 1 ns MD simulations revealed that the larger the void size, the easier is void collapse. Mean square displacements (MSDs) showed that HMX molecules that had collapsed into void present liquid structure characteristics. Four unique low-energy conformers were found for HMX molecules in void: two whose conformational geometries corresponded closely to those found in HMX polymorphs and two, additional, lower energy conformers that were not seen in the crystalline phases. The ratio of different conformers changed with the simulated temperature, in that the ratio of α conformer increased with the increase in temperature.

Surface-Accelerated Decomposition of δ-HMX

Despite extensive efforts to study the explosive decomposition of HMX, a cyclic nitramine widely used as a solid fuel, explosive, and propellant, an understanding of the physicochemical processes, governing the sensitivity of condensed HMX to detonation initiation is not yet achieved. Experimental and theoretical explorations of the initiation of chemistry are equally challenging because of many complex parallel processes, including the β-δ phase transition and the decomposition from both phases. Among four known polymorphs, HMX is produced in the most stable β-phase, which transforms into the most reactive δ-phase under heat or pressure. In this study, the homolytic NO2 loss and HONO elimination precursor reactions of the gas-phase, ideal crystal, and the (100) surface of δ-HMX are explored by first principles modeling. Our calculations revealed that the high sensitivity of δ-HMX is attributed to interactions of surfaces and molecular dipole moments. While both decomposition reactions coexist, the exothermic HONO-isomer formation catalyzes the N-NO2 homolysis, leading to fast violent explosions.

Vibrational spectra of an RDX film over an aluminum substrate from molecular dynamics simulations and density functional theory

We report calculated vibrational spectra in the range of 0-3,500 cm(-1) of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) molecules adsorbed on a model aluminum surface. A molecular film was modeled using two approaches: (1) density functional theory (DFT) was used to optimize a single RDX molecule interacting with its periodic images, and (2) a group of nine molecules extracted from the crystal structure was deposited on the surface and interacted with its periodic images via molecular dynamics (MD) simulations. In both cases, the molecule was initialized in the AAA conformer geometry having the three nitro groups in axial positions, and kept that conformation in the DFT examination, but some molecules were found to change to the AAE conformer (two nitro groups in axial and one in equatorial position) in the MD analysis. The vibrational spectra obtained from both methods are similar to each other, except in the regions where collective RDX intermolecular interactions (captured by MD simulations) are important, and compare fairly well with experimental findings.