Abstract:Computationally inexpensive particle-based coarse-grained (CG) models are essential for use in molecular dynamics (MD) simulations of mesoscopically slow cooperative phenomena, such as plastic deformations in solids. Molecular crystals possessing complex symmetry present enormous practical challenges for particle-based coarse-graining at molecularly resolved scales, when each molecule is in a single-site representation, and beyond. Presently, there is no published pairwise non-bonded single-site CG potential that is able to predict the space group and structure of a molecular crystal. In this paper, we present a successful coarse-graining at a molecular level from first principles of an energetic crystal, hexahydro-1,3,5-trinitro-s-triazine (RDX) in the alpha phase, using the force-matching-based multiscale coarse-graining (MSCG/FM) approach. The new MSCG/FM model, which implements an optimal pair decomposition of the crystal Helmholtz free energy potential in molecular center-of-mass coordinates, was obtained by force-matching atomistic MD simulations of liquid, amorphous, and crystalline states and in a wide range of pressures (up to 20 GPa). The MSCG/FM potentials for different pressures underwent top-down optimization to fine-tune the mechanical and thermodynamic properties, followed by consolidation into a transferable density-dependent model referred to as RDX-TC-DD (RDX True-Crystal Density-Dependent). The RDX-TC-DD model predicts accurately the crystal structure of α-RDX at room conditions and reproduces the atomistic reference system under isothermal (300 K) hydrostatic compression up to 20 GPa, in particular, the Pbca symmetry of α-RDX in the elastic regime. The RDX-TC-DD model was then used to simulate the plastic response of uniaxially ([100]) compressed α-RDX resulting in nanoscale shear banding, a key mechanism for plastic deformation and defect-free detonation initiation proposed for many molecular crystalline explosives. Additionally, a comparative analysis of the effect of core-softening of the RDX-TC-DD potential and the degree of molecular rigidity in the all-atom treatment suggests a stress-induced short-range softening of the effective intermolecular interaction as a fundamental cause of plastic instability in α-RDX. The reported RDX-TC-DD model and overall workflow to develop it open up possibilities to perform high quality simulation studies of molecular energetic materials under thermal and mechanical stimuli, including extreme conditions.

High-Fidelity Hugoniots of Α Phase RDX Solid from High-Quality Force Field with Thermal, Zero-Point Vibration, and Anharmonic Effects

Properties and Phase Transitions of the Solid Β-HMX: Different Force Fields

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

The first principle studies of the structural and vibrational properties of solid β-HMX under compression

The Pressure-Induced Phase Transition of the Solid -Hmx

Ab Initioand Molecular Dynamics Studies of Solid Β-Hmx: Effects of Hydrostatic Pressure and High Temperature

First-Principles Study of High-Pressure Behavior of Solid β-HMX

Accurately Predicting the Density and Hydrostatic Compression of Hexahydro-1,3,5-Trinitro-1,3,5-Triazine from First Principles

A unified model of anisotropy, thermoelasticity, inelasticity, phase transition and reaction for high-pressure ramp-loaded RDX single crystal

All-atom, Non-Empirical, and Tailor-Made Force Field for Α-Rdx from First Principles

Isentropic Compression Loading Locus of Α-Rdx Based on Its Self-consistent Force Field

Phase Transition, Thermodynamics Properties and IR Spectrum of Α- and Γ-Rdx: First Principles and MD Studies

Development and Validation of an All-Atom Force Field for the Energetic Materials TATB, RDX and HMX

First Principles Predict the Structural Properties, Electronic Properties and Phase Diagrams of RDX under High Temperature and Pressure

The Pressure Effects And Vibrational Properties Of Energetic Material: Hexahydro-1,3,5-Trinitro-1,3,5-Triazine (Alpha-Rdx)

Exploring the High-Pressure Polymorphs of Flexible Molecules by Crystal Structure Prediction: A Case of 1,3,5-Trinitro-1,3,5-triazinane

Many-Body Expansion-Based Quantum Mechanical Force Field for Cyclotrimethylene Trinitramine under High Pressure

Energetic Potential of Hexogen Constructed by Machine Learning

First-principles Study of Structural, Hydrogen Bonds and Mechanical Properties of Α-Rdx under Hydrostatic Compression

Md Simulation Study on the Mechanical Properties of Hmx Crystals and Hmx/F-2311 Pbxs

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