Reactive Molecular Dynamics Simulations of the Thermal Decomposition Mechanism of 1,3,3-Trinitroazetidine

Investigation of dislocation and twinning behavior in HMX under high-velocity impact employing molecular dynamics simulations

CONTEXT

Under the ReaxFF/lg force field, the multiscale shock technique (MSST) was employed to investigate the decomposition behavior of perfect, dislocated, and twinned octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) at a velocity of 11 km/s. This study aimed to analyze the changes in system temperature, bond formation and breaking, variations in the number of small molecules, and the number of clusters. The results indicated that the sensitivity of dislocated HMX was the lowest, while the sensitivity of twinned and perfect HMX was comparable. Comparing the formation and breaking of bonds in HMX during the shock process, it was found that the change in the number of bonds in dislocated HMX was similar to that in perfect HMX, whereas twinning accelerated the breaking of bonds. By analyzing the changes in small molecular fragments (CO, CO2, H, H2, H2O, N2, N2H, NH2, NO, NO2, and O) during the shock process of HMX, it was found that dislocation had a relatively minor effect on the small molecular fragments, while twinning promoted the generation of CO, H, NO, and O and accelerated the decomposition of NO. A comparison of the number, weight, and atomic ratio of clusters under perfect, dislocated, and twinned conditions revealed that under the influence of shock, the number of clusters initially increased sharply and then decreased slowly. Meanwhile, compared to the perfect and dislocated explosives, the number of clusters under the twinned structure was significantly fewer, indicating that the twinned structure could reduce cluster formation. The proportion of oxygen to carbon in the twinned HMX was lower than that in the perfect and dislocated explosives, possibly due to the higher content of small molecular fragment O in twinned HMX.

METHODS

Different structures of HMX crystals were constructed, including twinned defect structure (with a supercell containing 6458 atoms), dislocation defect structure (with a supercell containing 2352 atoms), and perfect structure (also with a supercell containing 2352 atoms). The modeling of defect crystal structures was carried out using the Atomsk software. For the twinned defect structure, we first constructed a mirror symmetric structure of the original configuration and then merged these two structures together. For the dislocation defect structure, we shifted a segment of the originally ordered perfect crystal structure by a certain distance using Atomsk. Before conducting the simulations, we performed geometric optimization of the models using the conjugate gradient (CG) algorithm and carried out 10 ps of NVT and NPT simulations to equilibrate the energy, temperature, and other parameters within the system. Finally, a 50-ps MSST impact simulation was performed using Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS).

Investigation of the decomposition mechanism of MTNP melt-cast explosive at different temperatures and pressures through ReaxFF/lg molecular dynamics simulations

CONTEXT

Thermal decomposition of 1-methyl-3,4,5-trinitropyrazole (MTNP), a melt-cast explosive, was investigated at different temperatures (2500, 2750, 3000, 3250, and 3500 K) and pressures (3000 K/0.5 GPa, 3000 K/1 GPa) using the ReaxFF/lg force field. The study aimed to analyze the changes in reactant quantities, initial reaction pathways, and final product yields. The results demonstrated that complete decomposition of MTNP molecules occurred within a timeframe of 200 ps, with shorter decomposition times observed as the temperature increased. The high-temperature thermal decomposition of MTNP was found to follow two primary reaction pathways. Reaction 1 involved denitration, while reaction 2 proceeded with nitro group isomerization. DFT calculations indicated that nitro group isomerization was the most favorable reaction. During the initial stages, higher quantities of NO2, NO, and N2 were observed compared to other species. This can be attributed to the relatively higher nitrogen and oxygen content in the MTNP structure. Among the five reaction temperatures, it was observed that the quantities of small molecules followed the order of NO2 > NO > N2 > CO. Moreover, with increasing temperature, the quantities of all four small molecules increased, indicating that higher temperatures promoted the progression of the reactions. However, as the pressure increased, there was a trend of initially increasing and then decreasing to zero for the quantities of NO2 and NO. This suggests that high temperature accelerated the high-temperature thermal decomposition of NO2 and NO, leading to a significant increase in the content of N2.

METHODS

A 3 × 5 × 5 supercell model of MTNP was constructed in Materials Studio, consisting of 75 unit cells and 300 MTNP molecules. The model was then subjected to a 20 ps geometric optimization using the Polak-Ribiere version of the conjugate gradient (CG) algorithm in the large-scale atomic/molecular massively parallel simulator (LAMMPS) under the isothermal-isobaric (NPT) ensemble at 1 atm pressure and 300 K temperature. Following the optimization, molecular dynamics simulations were performed on the model at five temperatures (2500, 2750, 3000, 3200, and 3500 K) under 1 atm using the NPT ensemble for a total duration of 1 ns. During the simulations, atomic trajectories, as well as information on atomic and molecular species, were output every 500 steps. Subsequently, a custom script was utilized to analyze the thermal decomposition pathways and products. A time step of 0.1 fs was employed for the calculations, and periodic boundary conditions were applied to eliminate boundary effects.

Coupling Effect of Non-Ignition Impact and Heat on the Decay of FOX-7

Non-ignition impact and heat stimuli are the most common external stimuli loaded on energetic materials. Nevertheless, there is thereby an urgent need, but it is still a significant challenge to comprehend their coupling effects on the decay and safety mechanisms of energetic materials. Then, reactive molecular dynamics simulation was employed to mimic practical situations and reveal the impact heat coupling effect on the decay mechanism of FOX-7. The temperature and the degree of compression of the crystal caused by the impact are considered variables in the simulation. Both increasing the degree of compression and elevating the temperature promotes the decay of FOX-7. However, their underlying response mechanism is not the same. The acceleration of decomposition is due to the elevated potential energy of the FOX-7 molecules because of elevating the temperature. In addition to the elevated potential energy of the molecule, the main contribution to the decomposition from the compression is to change the decomposition path. The results of the analysis show that compression reduces the stability of the C=C bond, so that chemical reactions related to the double bond occur. In addition, interestingly, the compression along the c direction has an almost equal effect on the final product as the compression along the b direction. Finally, the decay reaction networks are proposed to provide insights into the decomposition mechanism on atomic level. All these findings are expected to pave a way to understand the underlying response mechanism for the FOX-7 against external stimuli.

The Anisotropic Chemical Reaction Mechanism of 1,3,3-trinitroazetidine (TNAZ) under Different Shock Wave Directions by ReaxFF Reactive Molecular Dynamics Simulations

1,3,3-Trinitroazetidine (TNAZ) has good thermal stability and low shock sensitivity, among other properties, and it has broad prospects in insensitive ammunition applications. In this study, a molecular dynamics calculation based on the ReaxFF-lg force field and multiscale shock technique (MSST) was used to simulate the shock-induced chemical reaction of TNAZ with different shock wave directions. The results showed that the shock sensitivity of TNAZ was in the order of [100] > [010] > [001]. There were significant differences in molecular arrangements in different shock directions, which affected the reaction rate and reaction path in different directions. The molecular arrangement in the [010] and [001] directions formed a “buffer” effect. The formation and cleavage of bonds, formation of small molecules and growth of clusters were analyzed to show the effect of the “buffer”. The polymerization reactions in the [010] and [001] directions appeared later than that in the [100] direction, and the cluster growth in the [010] and [001] directions was slower than that in the [100] direction. In different shock loading directions, the formation and cleavage mechanisms of the N-O bonds of the TNAZ molecules were different, which resulted in differences in the initial reaction path and reaction rate in the three directions