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 Dynamics Simulations of the Collapse of a Cylindrical Pore in the Energetic Material α-RDX

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

Continuum and molecular dynamics simulations of pore collapse in shocked β-tetramethylene tetranitramine (β-HMX) single crystals

Pressure-Thresholded Response in Cylindrically Shocked Cyclotrimethylene Trinitramine (RDX)

Pore collapse in an energetic material from the micro-scale to the macro-scale

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

Influence of Polymer on Shock-Induced Pore Collapse: Hotspot Criticality through Reactive Molecular Dynamics

Reactive molecular dynamics simulations on the decomposition process of 1,3,5-trinitro-1,3,5-triazine crystal under high temperatures and pressure

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

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

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

Energy Localization Efficiency in TATB Pore Collapse Mechanisms

Atomistics-consistent continuum models and pore-collapse-generated hotspot temperatures in energetic crystals I: beta-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (beta-HMX)

Full-Field Quantitative Visualization of Shock-Driven Pore Collapse and Failure Modes in PMMA

Pressure-induced metallization of condensed-phase RDX: molecular dynamic simulations in conjunction with MSST method

High Pressure Suppression of Plasticity due to Over-Nucleation of Shear Strain

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

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

Molecular dynamics analysis of the transient temperature increase at void locations in shocked materials: RDX and Cu

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

The Potential Energy Hotspot: Effects from Impact Velocity, Defect Geometry, and Crystallographic Orientation