Abstract:We report molecular dynamics (MD) simulations using the first-principles-based ReaxFF reactive force field to study the thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) at various densities and temperatures. TATB is known to produce a large amount (15-30%) of high-molecular-weight carbon clusters, whereas detonation of nitramines such as HMX and RDX (1,3,5-trinitroperhydro-1,3,5-triazine) generate predominantly low-molecular-weight products. In agreement with experimental observation, these simulations predict that TATB decomposition quickly (by 30 ps) initiates the formation of large carbonaceous clusters (more than 4000 amu, or approximately 15-30% of the total system mass), and HMX decomposition leads almost exclusively to small-molecule products. We find that HMX decomposes readily on this time scale at lower temperatures, for which the decomposition rate of TATB is about an order of magnitude slower. Analyzing the ReaxFF MD results leads to the detailed atomistic structure of this carbon-rich phase of TATB and allows characterization of the kinetics and chemistry related to this phase and their dependence on system density and temperature. The carbon-rich phase formed from TATB contains mainly polyaromatic rings with large oxygen content, leading to graphitic regions. We use these results to describe the initial reaction steps of thermal decomposition of HMX and TATB in terms of the rates for forming primary and secondary products, allowing comparison to experimentally derived models. These studies show that MD using the ReaxFF reactive force field provides detailed atomistic information that explains such macroscopic observations as the dramatic difference in carbon cluster formation between TATB and HMX. This shows that ReaxFF MD captures the fundamental differences in the mechanisms of such systems and illustrates how the ReaxFF may be applied to model complex chemical phenomena in energetic materials. The studies here illustrate this for modestly sized systems and modest periods; however, ReaxFF calculations of reactive processes have already been reported on systems with approximately 10(6) atoms. Thus, with suitable computational facilities, one can study the atomistic level chemical processes in complex systems under extreme conditions.

Initial Chemical Events in Shocked Octahydro-1,3,5,7-tetranitro-1,3,5,7- Tetrazocine: A New Initiation Decomposition Mechanism

Initial decomposition mechanism of solid-state HMX for steady shock waves

Initial Decomposition of the Condensed-Phase β-HMX under Shock Waves: Molecular Dynamics Simulations

Shock Response of 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX): the C-N Bond Scission Studied by Molecular Dynamics Simulations

Initial Decomposition of the Co-crystal of CL-20/TNT: Sensitivity Decrease under Shock Loading

Anisotropic Responses and Initial Decomposition of Condensed-Phase Β-Hmx under Shock Loadings Via Molecular Dynamics Simulations in Conjunction with Multiscale Shock Technique.

Phase Transitions and Mechanical Properties of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine in Different Crystal Phases by Molecular Dynamics Simulation

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

Molecular dynamics simulation of primary detonation process of TATB crystal under shock loading

Initial Decomposition Mechanism of 3-Nitro-1,2,4-triazol-5-one (NTO) under Shock Loading: ReaxFF Parameterization and Molecular Dynamic Study

Reaction mechanism and kinetic modeling of gas-phase thermal decomposition of prototype nitramine compound HMX

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

Pressure-induced metallization of condensed phase β-HMX under shock loadings via molecular dynamics simulations in conjunction with multi-scale shock technique

Dynamic Responses and Initial Decomposition under Shock Loading: A DFTB Calculation Combined with MSST Method for β-HMX with Molecular Vacancy

Thermal Decomposition Mechanism of 1,3,5,7-Tetranitro-1,3,5,7-tetrazocane Accelerated by Nano-Aluminum Hydride (AlH 3 ): ReaxFF-Lg Molecular Dynamics Simulation

ReaxFF molecular dynamics simulations of shock induced reaction initiation in TNT

Carbon cluster formation during thermal decomposition of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine and 1,3,5-triamino-2,4,6-trinitrobenzene high explosives from ReaxFF reactive molecular dynamics simulations

Size-Dependent Shock Response Mechanisms in Nanogranular RDX: A Reactive Molecular Dynamics Study

Early Decay Mechanism of Shocked ε-CL-20: A Molecular Dynamics Simulation Study

Molecular Dynamics Simulations of an Initial Chemical Reaction Mechanism of Shocked CL-20 Crystals Containing Nanovoids

High Energy Barrier Hydroxyl Radical Dissociation Mechanism of Low Shock Sensitive Dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate (TKX-50) Explosive