Betsy M. Rice

For many experimentally measured chemical properties that cannot be directly computed from first-principles, the existing physics-based models do not extrapolate well to out-of-sample molecules, and experimental datasets themselves are too small for traditional machine learning (ML) approaches. To overcome these limitations, we apply a transfer learning approach, whereby we simultaneously train a multi-target regression model on a small number of molecules with experimentally measured values and a large number of molecules with related computed properties. We demonstrate this methodology on predicting the experimentally measured impact sensitivity of energetic crystals, finding that both characteristics of the computed dataset and model architecture are important to prediction accuracy of the small experimental dataset. Our directed-message passing neural network (D-MPNN) ML model using transfer learning outperforms direct-ML and physics-based models on a diverse test set, and the new methods described here are widely applicable to modeling many other structure-property relationships.

Availability of heat capacity as a function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5‐trinitro‐1,3,5‐triazinane (RDX) and 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant‐volume and constant‐pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

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 new approach to identify chemical species from molecular dynamics (MD) simulations of reacting materials under extreme temperatures and pressures is presented. The approach is based on bond-distance and vibrational criteria, derived from the examination of atomic behavior during a density functional theory MD simulation of an overdriven shock of the explosive pentaerythritol tetranitrate. For comparison, the trajectory was analyzed using popular bonding criteria commonly used in analysis of reactive MD simulations, including distance, distance–time, and bond-order criteria. Cluster analyses using the new time-dependent bond definition approach presented here and a bond-order approach revealed that species and their corresponding lifetimes were strongly dependent on the chosen approach, indicating significant implications for the development of chemical mechanisms and chemical kinetics models using the results of reactive MD simulations.

The response of an energetic material to insult is perhaps one of the most difficult processes to model due to concurrent chemical and physical phenomena occurring over scales ranging from atomistic to continuum. Unraveling the interdependencies of these complex processes across the scales through modeling can only be done within a multiscale framework. In this paper, I will describe progress in the development of a predictive, experimentally validated multiscale reactive modeling capability for energetic materials at the Army Research Laboratory. I will also describe new challenges and research opportunities that have arisen in the course of our development which should be pursued in the future.