Cargando…

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

We report a reactive molecular dynamic (ReaxFF-MD) study using the newly parameterized ReaxFF-lg reactive force field to explore the initial decomposition mechanism of 3-Nitro-1,2,4-triazol-5-one (NTO) under shock loading (shock velocity >6 km/s). The new ReaxFF-lg parameters were trained from ma...

Descripción completa

Detalles Bibliográficos
Autores principales: Du, Lixiaosong, Jin, Shaohua, Nie, Pengsong, She, Chongchong, Wang, Junfeng
Formato: Online Artículo Texto
Lenguaje:English
Publicado: MDPI 2021
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8400475/
https://www.ncbi.nlm.nih.gov/pubmed/34443396
http://dx.doi.org/10.3390/molecules26164808
Descripción
Sumario:We report a reactive molecular dynamic (ReaxFF-MD) study using the newly parameterized ReaxFF-lg reactive force field to explore the initial decomposition mechanism of 3-Nitro-1,2,4-triazol-5-one (NTO) under shock loading (shock velocity >6 km/s). The new ReaxFF-lg parameters were trained from massive quantum mechanics data and experimental values, especially including the bond dissociation curves, valence angle bending curves, dihedral angle torsion curves, and unimolecular decomposition paths of 3-Nitro-1,2,4-triazol-5-one (NTO), 1,3,5-Trinitro-1,3,5-triazine (RDX), and 1,1-Diamino-2,2-dinitroethylene (FOX-7). The simulation results were obtained by analyzing the ReaxFF dynamic trajectories, which predicted the most frequent chain reactions that occurred before NTO decomposition was the unimolecular NTO merged into clusters ((C(2)H(2)O(3)N(4))(n)). Then, the NTO dissociated from (C(2)H(2)O(3)N(4))(n) and started to decompose. In addition, the paths of NO(2) elimination and skeleton heterocycle cleavage were considered as the dominant initial decomposition mechanisms of NTO. A small amount of NTO dissociation was triggered by the intermolecular hydrogen transfer, instead of the intramolecular one. For α-NTO, the calculated equation of state was in excellent agreement with the experimental data. Moreover, the discontinuity slope of the shock-particle velocity equation was presented at a shock velocity of 4 km/s. However, the slope of the shock-particle velocity equation for β-NTO showed no discontinuity in the shock wave velocity range of 3–11 km/s. These studies showed that MD by using a suitable ReaxFF-lg parameter set, could provided detailed atomistic information to explain the shock-induced complex reaction mechanisms of energetic materials. With the ReaxFF-MD coupling MSST method and a cheap computational cost, one could also obtain the deformation behaviors and equation of states for energetic materials under conditions of extreme pressure.