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.

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

Characterization of Nano-Scale Parallel Lamellar Defects in RDX and HMX Single Crystals by Two-Dimension Small Angle X-ray Scattering

Electron Microscopy and Diffraction of Crystalline Dendrimers and Macrocycles

Mechanical and thermophysical properties of energetic crystals: evaluation methods and recent achievements

Exposed Facets and Surface Properties of Highly Explosive RDX Studied by the Crystal-Face-Indexing Method and Theoretical Calculations

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

Generation and Evolution Mechanism of Cracks in HMX Single Crystals Around Phase Transition Temperature

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

Experimental Study of the Crystal Habit of High Explosive Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) in Acetone and Dimethyl Sulfoxide

Research on the Mechanism of Threefold Discrepancy in Thermal Conductivity Between Room Temperature and Phase Transition Temperature for 1,3,5,7-Tetranitro-1,3,5,7-tetrazocine (HMX) Crystals.

Influence of the shape of RDX grains on the creation of hot spots by mesoscale modeling

Effects of Grain Refinement and Thermal Aging on Atomic Scale Local Structures of Ultra-Fine Explosives by X-ray Total Scattering

Researches and Progresses of Novel Energetic Materials

Molecular Dynamics Study of RDX/AMMO Propellant

Crystal Quality and Properties of Spherical HMX

Evaluation of Crystal Properties and Initiation Characteristics of Decreased Sensitivity RDX

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

Molecular dynamics simulations of void defects in the energetic material HMX

Insights into the Mechanical-Structural-Energetic Responses of RDX Nanoparticles from Molecular Dynamics Simulation

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

The effect of crystal anisotropy and plastic response on the dynamic fracture of energetic materials