▪ Abstract  Energetic materials are chemical compounds or mixtures that store significant quantities of energy. In this review, we explore recent approaches to property prediction and new material synthesis. We show how the successful design of new energetic materials with tailored properties is becoming a practical reality.

We report the first quantum-based multiscale simulations to study the reactivity of shocked perfect crystals of the insensitive energetic material triaminotrinitrobenzene (TATB). Tracking chemical transformations of TATB experiencing overdriven shock speeds of 9 km/s for up to 0.43 ns and 10 km/s for up to 0.2 ns reveal high concentrations of nitrogen-rich heterocyclic clusters. Further reactivity of TATB toward the final decomposition products of fluid N(2) and solid carbon is inhibited due to the formation of these heterocycles. Our results thus suggest a new mechanism for carbon-rich explosive materials that precedes the slow diffusion-limited process of forming the bulk solid from carbon clusters and provide fundamental insight at the atomistic level into the long reaction zone of shocked TATB.

Raman spectroscopy in a laser heated diamond anvil cell and first principles molecular dynamics simulations have been used to study water in the temperature range 300 to 1500 K and at pressures to 56 GPa. We find a substantial decrease in the intensity of the O-H stretch mode in the liquid phase with pressure, and a change in slope of the melting line at 47 GPa and 1000 K. Consistent with these observations, theoretical calculations show that water beyond 50 GPa is "dynamically ionized" in that it consists of very short-lived (<10 fs) H2O, H3O+, and OH- species, and also that the mobility of the oxygen ions decreases abruptly with pressure, while hydrogen ions remain very mobile. We suggest that this regime corresponds to a superionic state.

Delivery of prebiotic compounds to early Earth from an impacting comet is thought to be an unlikely mechanism for the origins of life because of unfavourable chemical conditions on the planet and the high heat from impact. In contrast, we find that impact-induced shock compression of cometary ices followed by expansion to ambient conditions can produce complexes that resemble the amino acid glycine. Our ab initio molecular dynamics simulations show that shock waves drive the synthesis of transient C-N bonded oligomers at extreme pressures and temperatures. On post impact quenching to lower pressures, the oligomers break apart to form a metastable glycine-containing complex. We show that impact from cometary ice could possibly yield amino acids by a synthetic route independent of the pre-existing atmospheric conditions and materials on the planet.

A new multiscale simulation method is formulated for the study of shocked materials. The method combines molecular dynamics and the Euler equations for compressible flow. Treatment of the difficult problem of the spontaneous formation of multiple shock waves due to material instabilities is enabled with this approach. The method allows the molecular dynamics simulation of the system under dynamical shock conditions for orders of magnitude longer time periods than is possible using the popular nonequilibrium molecular dynamics approach. An example calculation is given for a model potential for silicon in which a computational speedup of 10(5) is demonstrated. Results of these simulations are consistent with the recent experimental observation of an anomalously large elastic precursor on the nanosecond time scale.