How does low temperature coupled with different pressures affect initiation mechanisms and subsequent decompositions in nitramine explosive HMX?

Exploring the Interfacial Reaction of Nano Al/CuO Energetic Films through Thermal Analysis and Ab Initio Molecular Dynamics Simulation

The effect of the interface layer on energy release in nanoenergetic composite films is important and challenging for the utilization of energy. Nano Al/CuO composite films with different modulation periods were prepared by magnetron sputtering and tested by differential scanning calorimetry. With the increase in the modulation period of the nano Al/CuO energetic composite films, the interface layer contained in the energetic composite film decreased meaningfully, increasing the total heat release meaningfully. Ab initio molecular dynamics (AIMD) simulation were carried out to study the preparation process changes and related properties of the nano Al/CuO energetic composite films under different configurations at 400 K. The results showed that the diffusion of oxygen atoms first occurred at the upper and lower interfaces of CuO and Al, forming AlO_x and Cu_xAl_yO_z. The two-modulation-period structure changed more obviously than the one-modulation-period structure, and the reaction was faster. The propagation rate and reaction duration of the front end of the diffusion reaction fronts at the upper and lower interfaces were different. The Helmholtz free energy loss of the nano Al/CuO composite films with a two-modulation-period configuration was large, and the number of interfacial layers had a great influence on the Helmholtz free energy, which was consistent with the results of the thermal analysis. Current molecular dynamics studies may provide new insights into the nature and characteristics of fast thermite reactions in atomic detail.

Pressure and Polymorph Dependent Thermal Decomposition Mechanism of Molecular Materials: A Case of 1,3,5,-Trinitro-1,3,5,-triazine

1,3,5,-Trinitro-1,3,5,-triazine (RDX) serves as an important energetic material and is widely used as various solid propellants and explosives. Understanding the thermal decomposition behaviors of various polymorphs of RDX at high pressure and high temperature is significantly important for safe storage and handling. The present work reveals the early thermal decay mechanisms of two polymorphs (α- and ε-forms) of RDX at high pressure and high temperature by ReaxFF reactive molecular dynamic simulations and climbing image nudged elastic band (CI-NEB) static calculations. It is found that the thermal decomposition rate has positive and negative effects on the pressure for α- and ε-RDX, respectively. This difference originates from the difference of pressure effect on the intermolecular H transfer of the two polymorphs, as we confirm that the bimolecular H transfer rather than the NO₂ partition initiates the decay with a significantly lower energy barrier therein. This finding may be beneficial to understand the pressure and polymorph dependent effect on the decay of RDX and to develop a kinetic model for the combustion of solid RDX.

Pressure effects on the thermal decomposition of the LLM-105 crystal

Thermal mechanical responses under high temperature and high pressure are basic information to understand the performance of energetic materials. In this work, the pressure effects on the thermal decay of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) are explored. Up to the initial pressure of 4.6 GPa, the pressure dependent decomposition boundary is built and no phase transition occurs until the decomposition of the LLM-105 crystal. The decomposition temperature is significantly lifted via a weak loading pressure. The experimental measurement confirms the decomposition products, including NO₂, CO₂ and NH₃, which are predicted by the density functional tight-binding molecular dynamics method. The calculation described the details of thermal decay in the initial stages under high pressure. The sudden drop in the shifts of the Raman modes associated with hydrogen bonds under high pressure indicates the strengthening of the intermolecular hydrogen bonds and the occurrence of intermolecular hydrogen transfer prior to crystal decomposition. The simulation supported the existence of intermolecular hydrogen transfer and provided the transfer path and decomposition mechanism. All of these jobs not only contribute significantly to the understanding of thermal decomposition of energetic materials as a function of pressure, but also contribute to the understanding of sensitivity mechanisms and safety issues.

Predicting the Initial Thermal Decomposition Path of Nitrobenzene Caused by Mode Vibration at Moderate-Low Temperatures: Temperature-Dependent Anti-Stokes Raman Spectra Experiments and First-Principals Calculations

The lack of understanding of the initial decomposition micromechanism of energetic materials subjected to external stimulation has hindered its safe storage, usage, and development. The initial thermal decomposition path of nitrobenzene triggered by molecular thermal motion is investigated using temperature-dependent anti-Stokes Raman spectra experiments and first-principles calculations to clarify the initial thermal decomposition micromechanism. The experiment shows that the symmetric nitro stretching, antisymmetric nitro stretching, and phenyl ring stretching vibration modes are active as increasing temperature below 500 K. The DFT method is used to examine the effects of the three mode vibrations on the initial decomposition of nitrobenzene by relaxed scan for each relevant change in bond lengths and bond angles to obtain the optimal reaction channel leading to initial thermal decomposition of nitrobenzene. The results demonstrate that the initial thermal decomposition is the isomerization of nitrobenzene to phenyl nitrite. The optimal reaction channel leading to the initial isomerization is the increase or decrease of angle O-N-C from the antisymmetric nitro stretching vibration, which causes the torsion of nitro group and the subsequent oxygen atom attacking carbon atom. The scanning energy barrier related to angle O-N-C is about 62.1 kcal/mol, which is very consistent with the calculated activation barrier of isomerization of nitrobenzene. This proves the reliability of our conclusions.

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.

Thermal decomposition of isolated and crystal 4,10-dinitro-2,6,8,12-tetraoxa-4,10-diazaisowurtzitane according to ab initio molecular dynamics simulations

We performed ab initio molecular dynamics simulations to study the initiation chemical reaction and subsequent decomposition mechanism of a 4,10-dinitro-2,6,8,12-tetraoxa-4,10-diazaisowutrzitane (TEX) crystal at 2160 K.

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.

Periodic density functional theory study of the high-pressure behavior of crystalline L-serine-L-ascorbic acid

A detailed study of the structural, electronic and absorption properties of crystalline L-serine-L-ascorbic acid (SAA) in the pressure range of 0-300 GPa was performed by density-functional theory (DFT) calculations in this work. Our results show that the compressible crystal of SAA is anisotropic. Furthermore, specific analysis of the variation tendencies of bond lengths and bond angles under different pressures show that the main structural transformations occur at pressures of 40, 50, 70, 100, 130 and 150 GPa, accompanied by repeated formations and disconnections of covalent bonds between O2(P1) and C2(P2) as well as C3(P1) and O1(P2), and a newly formed five-atom ring at 100 GPa. In addition, from 40 to 230 GPa, complex hydrogen bond transformations occur in SAA under compression, while from 240 to 300 GPa, the curve of lattice constants, bond lengths and bond angles of SAA barely changes, suggesting structural stability after 230 GPa. Then, by analyzing the band gap and density of states of SAA, it was found that the crystal undergoes a phase transformation from insulator to semiconductor at 150 GPa and it becomes more sensitive under compression. In addition, a relatively high optical activity with the pressure increases of SAA was seen from the absorption spectra, and two obvious changes of absorption coefficients were also observed at 50 GPa and 130 GPa, respectively, indicating that structural transformations occur here. Graphical abstract Structural formation and breaking of the five-atom ring O1(P2)-C2(P2)-O2(P1)-C2(P1)-C3(P1) with increasing pressure.

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.