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.

Impact Sensitivity in Respect of the Crystal Lattice Free Volume and the Characteristics of Plasticity of Some Nitramine Explosives

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

Initial Decomposition of the Co-crystal of CL-20/TNT: Sensitivity Decrease under Shock Loading

ReaxFF Molecular Dynamics Simulation of Hydrostatic and Uniaxial Compression of Nitrate Energetic Materials

Synthesis, Performance, and Thermal Behavior of a Novel Insensitive EDNA/DAT Co‐crystal

Molecular simulation of the influence of interface faceting on the shock sensitivity of a model plastic bonded explosive

Sensitivity of Energetic Materials: Theoretical Relationships to Detonation Performance and Molecular Structure

Investigation of the microstructure of energetic crystals by means of X-ray powder diffraction

Mechanism of the Impact-Sensitivity Reduction of Energetic CL-20/TNT Cocrystals: A Nonequilibrium Molecular Dynamics Study

Johnson–Cook yield functions for cyclotetramethylene-tetranitramine (HMX) and cyclotrimethylene-trinitramine (RDX) derived from single crystal plasticity models

Evaluation of the compatibility of energetic nitrocellulose/plasticizer blends through molecular dynamics simulation

Polymorphism in a Nonsensitive-High-Energy Material: Discovery of a New Polymorph and Crystal Structure of 4,4′,5,5′-Tetranitro-1 H ,1′ H -[2,2′-biimidazole]-1,1′-diamine

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

Computational analysis of mesoscale thermomechanical ignition behavior of impacted LLM-105 based explosives

Investigation of Ethylenedinitramine as a Versatile Building Block in Energetic Salts, Cocrystals, and Coordination Compounds

Theoretical investigations on a novel CL-20/LLM-105 cocrystal explosive by molecular dynamics method

The influence of temperature and component proportion on stability, sensitivity, and mechanical properties of LLM-105/HMX co-crystals via molecular dynamics simulation

Co-agglomerated crystals of cyclic nitramines with the nitrogen rich 3,6-bis(1H-1,2,3,4-tetrazol-5-ylamino)-1,2,4,5-tetrazine (BTATz)

Study on Mechanical Properties and Penetration Stability of DNAN Based Explosives Under Additive Effects

Rheological Behavior of DNP/HMX Melt-Cast Explosives with Bimodal and Trimodal Particle-Size Distributions

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