A New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Matthew Kroonblawd,Paul Lafourcade,Laurence Fried,Jean-Bernard Maillet,Tommy Sewell

DOI: https://doi.org/10.26434/chemrxiv-2023-rql38

2023-12-13

Abstract:Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics (MD) modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of the single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature-pressure-dependent equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the existing nonreactive FFs are more accurate in describing TATB's physical properties than the ReaxFF-based variants. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. This work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

Chemistry

What problem does this paper attempt to address?

The paper aims to address the following issues: 1. **Development of a new non-reactive force field**: For molecular dynamics simulations of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under extreme conditions, existing force field models are not accurate enough in describing properties under high-pressure conditions. Therefore, researchers have developed a new non-reactive force field (CEA-LLNL-Missouri FF or CLM FF) to improve the accurate modeling of TATB material properties. 2. **Validation of existing force field models**: A series of experimental data and electronic structure calculation results were collected as a benchmark dataset, and these data were used to evaluate four existing TATB force field models as well as the newly developed CLM FF model. The study found that existing non-reactive force fields are more accurate in describing the physical properties of TATB compared to variants based on ReaxFF. 3. **Improvement of TATB's elasticity and equation of state**: The new CLM FF shows significant improvements in describing the elastic properties and equation of state of TATB, which helps to enhance the accuracy of complex dynamic response based on force field models, thereby better understanding and predicting the safety and performance characteristics of this material.

A New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

A New Method for Calculation of Elastic Properties of Anisotropic Material by Constant Pressure Molecular Dynamics

Molecular dynamics informed calibration of crystal plasticity critical shear stresses for the mesoscopic mechanical modeling of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) single crystal

Towards the Low-Sensitive and High-Energetic Co-Crystal Explosive CL-20/TNT: from Intermolecular Interactions to Structures and Properties

Ab Initio Neural Network MD Simulation of Thermal Decomposition of a High Energy Material CL-20/TNT.

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

Explosive chemistry: Simulating the chemistry of energetic materials at extreme conditions

Bottom-up coarse-grain modeling of nanoscale shear bands in shocked α-RDX

TATB thermal decomposition: An improved kinetic model for explosive safety analysis

The Reactivity of Energetic Materials Under High Pressure and Temperature

Bottom-up coarse-grain modeling of plasticity and nanoscale shear bands in α-RDX

Evaluating shock sensitivity and decomposition of energetic materials by ReaxFF molecular dynamics

First principles reactive simulation for equation of state prediction

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

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

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

Layers can be deceiving: A hopping model for small molecule diffusion in TATB crystal

Non-adiabatic Dynamics Modeling Framework for Materials in Extreme Conditions

A Study of the Shock Sensitivity of Energetic Single Crystals by Large-Scale Ab Initio Molecular Dynamics Simulations