Anisotropic hydrogen bond structures and orientation dependence of shock sensitivity in crystalline 1,3,5-tri-amino-2,4,6-tri-nitrobenzene (TATB)

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.

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

Engineering the Microstructure and Morphology of Explosive Films via Control of Interfacial Energy

Physical vapor deposition of organic explosives enables growth of polycrystalline films with a unique microstructure and morphology compared to the bulk material. This study demonstrates the ability to control crystal orientation and porosity in pentaerythritol tetranitrate films by varying the interfacial energy between the substrate and the vapor-deposited explosive. Variation in density, porosity, surface roughness, and optical properties is achieved in the explosive film, with significant implications for initiation sensitivity and detonation performance of the explosive material. Various surface science techniques, including angle-resolved X-ray photoelectron spectroscopy and multiliquid contact angle analysis, are utilized to characterize interfacial characteristics between the substrate and explosive film. Optical microscopy and scanning electron microscopy of pentaerythritol tetranitrate surfaces and fracture cross sections illustrate the difference in morphology evolution and the microstructure achieved through surface energy modification. X-ray diffraction studies with the Tilt-A-Whirl three-dimensional pole figure rendering and texture analysis software suite reveal that high surface energy substrates result in a preferred (110) out-of-plane orientation of pentaerythritol tetranitrate crystallites and denser films. Low surface energy substrates create more randomly textured pentaerythritol tetranitrate and lead to nanoscale porosity and lower density films. This work furthers the scientific basis for interfacial engineering of polycrystalline organic explosive films through control of surface energy, enabling future study of dynamic and reactive detonative phenomena at the microscale. Results of this study also have potential applications to active pharmaceutical ingredients, stimuli-responsive polymer films, organic thin film transistors, and other areas.

Anisotropic response of the co-crystal of CL-20/TNT under shock loading: molecular dynamics simulation

Anisotropic response of the co-crystal of CL-20/TNT under shock loading is studied by analyzing the changes of thermodynamic parameters, product evolution and cluster evolution.