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.

A study of the solvent effect on the morphology of RDX crystal by molecular modeling method

Molecular dynamics investigation of the effect of solvent adsorption on crystal habits of hexogen

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

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

Effect of the Casting Solvent on the Crystal Characteristics and Pervaporative Separation Performances of P(Vdf-Co-Hfp) Membranes

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

Study on the Effect of Solvent on Cocrystallization of CL-20 and HMX Through Theoretical Calculations and Experiments.

Exploring the mechanisms of benzanilide crystal growth and morphology: Crystal surface structure, solvent diffusion and solvent-interface interactions

Theoretical study on the influence of different solvents on the growth morphology of 3,4-bis(3-nitrofurazan-4-yl)furoxan (DNTF)

Morphology prediction of 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) crystal in dimethyl sulfoxide (DMSO) solvent with different models using molecular dynamics simulation

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

Theoretical Study on the Effect of Solvent Behavior on Ammonium Dinitramide (ADN)/1,4,7,10,13,16‐Hexaoxacyclooctadecane (18‐Crown‐6) Cocrystal Growth Morphology at Different Temperatures

The effect of solvents on morphology of sulfamerazine crystals and its molecular mechanism

Experimental and simulation study on hydrogen-bond-induced crystallization of spherical ammonium dinitramide

The Effect of Solvents on the Crystal Morphology of Pyriproxyfen

Theoretical study on initial decomposition paths of energetic materials featuring RDX ring

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

Theoretical calculation into the effect of molar ratio on the structures, stability, mechanical properties and detonation performance of 1,3,5,7-tetranitro-1,3,5,7-tetrazocane/ 1,3,5-trinitro-1,3,5-triazacyco-hexane cocrystal

Theoretical simulation study on crystal property and hygroscopicity of ADN doping with nitramine explosives (RDX, HMX, and CL-20)

Solubility determination and comparison of β-HMX and RDX in two binary mixed solvents (acetonitrile + water, nitric acid + water)

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