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.

Molecular dynamic simulations on TKX-50/RDX cocrystal.

Molecular dynamic simulations on TKX-50/HMX cocrystal

Molecular Dynamics Simulations on Dihydroxylammonium 5,5′-Bistetrazole-1,1′-diolate/hexanitrohexaazaisowurtzitane Cocrystal

Reactive molecular dynamics simulations on the thermal decompositions and oxidations of TKX-50 and twinned TKX-50

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

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

Molecular Simulation on Co-crystal Structure of HMX/TATB

Comparative studies on structure, sensitivity and mechanical properties of CL-20/DNDAP cocrystal and composite by molecular dynamics simulation

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

Anisotropic Impact Sensitivity and Shock Induced Plasticity of TKX-50 (dihydroxylammonium 5,5′-Bis(tetrazole)-1,1′-diolate) Single Crystals: from Large-Scale Molecular Dynamics Simulations

Crystal Morphology of Dihydroxylammonium 5,5′-Bistetrazole-1,1′-diolate (TKX-50) under Solvents System with Different Polarity Using Molecular Dynamics

Molecular dynamics simulations for 5,5′-bistetrazole-1,1′-diolate (TKX-50) and its PBXs

Spheroidal Crystal Morphology of RDX in Mixed Solvent Systems Predicted by Molecular Dynamics

Additives Effects on Crystal Morphology of Dihydroxylammonium 5,5ʹ-Bistetrazole-1,1ʹ-diolate by Molecular Dynamics Simulations

Theoretical study on intra-molecule interactions in TKX-50

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

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

Research on Structures, Mechanical Properties, and Mechanical Responses of TKX-50 and TKX-50 Based PBX with Molecular Dynamics

Microscale Investigations of Mechanical Responses of TKX-50 Based Polymer Bonded Explosives Using MD Simulations

Morphology prediction of dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate (TKX-50) crystal in different solvent systems using modified attachment energy model

Directly Insight Into The Inter- And Intramolecular Interactions Of Cl-20/Tnt Energetic Cocrystal Through The Theoretical Simulations Of Thz Spectroscopy