Molecular conformational stability in cyclotrimethylene trinitramine crystals

Properties of Erythritol Tetranitrate from Molecular Dynamics Simulation

The nonpolarizable force field for alkyl nitrates developed by Borodin et al. [J. Phys. Chem. B, 2008, 112, 734-742] has been employed to calculate selected properties of crystalline and liquid erythritol tetranitrate (ETN). The set of partial charges proposed by Borodin for pentaerythritol tetranitrate (PETN) was used except for a small correction to the H atom charges to ensure charge neutrality owing to the absence of the neopentyl carbon in ETN. The force field was used to compute the isothermal compression curve, lattice parameters, heat capacity, thermal expansivity, single crystal elastic constants, and Gruneisen parameters of crystalline ETN. The density- and temperature-dependent viscosities of liquid ETN are also reported. We anticipate that these data will be of some utility to the development of equations of state and thermomechanical models for ETN.

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

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 theoretical study of the relaxation of a phenyl group chemisorbed to an RDX freestanding thin film

Energy relaxation from an excited phenyl group chemisorbed to the surface of a crystalline thin film of α-1,3,5-trinitro-1,3,5-triazacyclohexane (α-RDX) at 298 K and 1 atm is simulated using molecular dynamics. Two schemes are used to excite the phenyl group. In the first scheme, the excitation energy is added instantaneously as kinetic energy by rescaling momenta of the 11 atoms in the phenyl group. In the second scheme, the phenyl group is equilibrated at a higher temperature in the presence of static RDX geometries representative of the 298 K thin film. An analytical model based on ballistic phonon transport that requires only the harmonic part of the total Hamiltonian and includes no adjustable parameters is shown to predict, essentially quantitatively, the short-time dynamics of the kinetic energy relaxation (∼200 fs). The dynamics of the phenyl group for times longer than about 6 ps follows exponential decay and agrees qualitatively with the dynamics described by a master equation. Long-time heat propagation within the bulk of the crystal film is consistent with the heat equation.

Rotational defects in cyclotrimethylene trinitramine (RDX) crystals

Cyclotrimethylene trinitramine (RDX) crystalizes in the orthorhombic α-phase at the ambient pressure and temperature. In principle, the point defects commonly found in monatomic crystals, such as vacancies and interstitials, may exist in RDX as well. However, in molecular crystals one encounters additional point defects associated with the distortion of the molecules. A set of rotational defects are described in this article. These are molecules which are located in the proper positions in the crystal but are rotated relative to the molecules in the perfect crystal, and their ring is slightly puckered. The energetic barriers for defect formation and for their annealing back to the perfect crystal configuration are computed using an atomistic model. It is shown that the formation energy of rotational defects is smaller than the vacancy formation energy. Such defects are identified in the cores of dislocations in RDX and hence their concentration in the crystal is expected to increase during plastic deformation. The importance of such point defects is related to their role in phonon scattering and in dislocation-mediated plastic deformation.