The ease with which an energetic material can be initiated by mechanical impact is a critical parameter directing material safety and application. While impact sensitivity metrics are traditionally derived experimentally, recent developments have highlighted that the phenomenon is amenable to first principles simulation. In this chapter, we will outline a fully ab initio approach to predict the relative impact sensitivities of energetic materials based on the mechanochemical principles that link the impact event to vibrational energy transfer.
This mechanism is key to rationalizing how a mechanical impact—which deposits energy into the low-frequency lattice vibrations—results in a molecular response. By simulating the vibrational energy levels (the so-called phonon density of states, PDOS) using first-principles computational methods (typically dispersion-corrected plane-wave density functional theory, PW-DFT) we can calculate the relative rate of energy propagation from the delocalized low-energy lattice vibrations through to the localized molecular modes. The latter traps the energy, which eventually results in bond rupture through heightened vibrational excitation.
This method, based on vibrational up-pumping, offers a route toward predicting the impact sensitivities of a broad range of energetic materials, provided the crystal structure of the compound (or salt or co-crystal) is known. While it does not offer insight into the sensitizing roles undoubtedly played by crystal defects or grain boundaries, it does provide a level of understanding at the molecular and crystal packing levels. Correspondingly, this approach offers a feedback mechanism to chemists and materials scientists to guide the design of new materials with desired impact sensitivity behavior.
The reliable determination of gas-phase and solid-state heats of formation are important considerations in energetic materials research. Herein, the ability of PM7 to calculate the gas-phase heats of formation for CNHO-only and inorganic compounds has been critically evaluated, and for the former, comparisons drawn with isodesmic equations and Atom equivalence methods. Routes to obtain solid-state heats of formation for a range of singlecomponent molecular solids, salts, and co-crystals were also evaluated. Finally, local vibrational mode analysis has been used to calculate bond length/force constant curves for seven different chemical bonds occurring in CHNO-containing molecules, which allow for rapid identification of the weakest bond, opening up great potential to rationalise decomposition pathways. Both metrics are important tools in rationalising the design of new energetic materials through computational screening processes.
The impact sensitivity (IS) of FOX-7 polymorphs is predicted by phonon up-pumping to decrease as layers of FOX-7 molecules flatten. Experimental validation proved anomalous owing to a phase transition during testing, raising questions regarding Impact sensitivity measurement and highlighting the need for models to predict IS of polymorphic energetic materials.
Predicting the Impact Sensitivities of Energetic Materials through Zone-Center Phonon Up-Pumping
(2021)
The development of new energetic materials (EMs) is accompanied by significant hazards, prompting interest in their computational design. Before reliable in silico design strategies can be realized, however, approaches to understand and predict EM response to mechanical impact must be developed. We present here a fully ab initio model based in phonon up-pumping which successfully ranks the relative impact sensitivity of a series of organic EMs. The methodology depends only on the crystallographic unit cell and Brillouin zone center vibrational frequencies. We therefore expect this approach to become an integral tool in the large-scale screening of potential EMs.
The ease with which an energetic material (explosives, propellants, and pyrotechnics) can be initiated is a critical parameter to assess their safety and application. Impact sensitivity parameters are traditionally derived experimentally, at great cost and risk to safety. In this work we explore a fully ab initio Approach based on concepts of vibrational energy transfer to predict impact sensitivities for a series of chemically, structurally and energetically diverse molecular materials. The quality of DFT calculations is assessed for a subset of the materials by comparison with experimental inelastic neutron scattering spectra (INS). A variety of models are considered, including both qualitative and quantitative analysis of the vibrational spectra. Excellent agreement against experimental impact sensitivity is achieved by consideration of a multi-phonon ladder-type up-pumping mechanism that includes both overtone and combination pathways, and is improved further by the added consideration of temperature. This fully ab initio approach not only permits ranking of energetic materials in terms of their impact sensitivity but also provides a tool to guide the targeted design of advanced energetic compounds with tailored properties.
A high-pressure neutron diffraction study was conducted on polycrystalline samples of the two known polymorphs of 2,4,6-trinitrotoluene [monoclinic (m) and orthorhombic (o) TNT] under hydrostatic conditions. Isothermal equations of state were obtained for both polymorphic forms. Neither polymorph was observed to undergo a phase transition in the pressure region 0–5 GPa, with both polymorphs displaying smooth compression behavior across the pressure range. This differs somewhat from previous X-ray diffraction and Raman spectroscopy investigations in which discontinuities were observed in the P–V curves and spectral changes were reported at ∼2 GPa. The high-pressure response of these materials is supported by dispersion-corrected density functional theory calculations which, while overestimating the experimental bulk moduli, give excellent agreement with the observed smooth compression response of both phases.