Can a shock-induced phonon up-pumping model relate to impact sensitivity of molecular crystals, polymorphs and cocrystals?

Impact sensitivity engineering of high-energy molecular crystals requires accurate predictive models. For this purpose, the promising multi-phonon based approach is selected, assessing a bit more its strengths and weaknesses. Presently used with high-quality phonon calculations of 22 molecular crystals, using a physics-based criterion to determine the phonon bath extent, the resulting intrinsic shock sensitivity index (SSI) is compared to the most common marker of impact sensitivity, h 50, as determined from drop-weight impact tests. Selecting a data subset from experiments performed under very similar conditions (2.5 kg hammer with grit and 30-40 mg samples), the model can predict h 50 values for mono-molecular crystals with very good accuracy, including the ability to discriminate the polymorphs of HMX and CL20. This very good agreement validates an initial indirect up-pumping mechanism occurring under these conditions, where the doorway modes also interact with the phonon bath. However, the phonon bath criterion for mono-molecular crystals does not transfer well to cocrystals. Owing to the vibrational coupling of the co-molecules, it seems a broader phonon bath should be considered. Additionally recalling experimental uncertainty and various experimental factors affecting h 50 values for a given compounds, we recommend that the density of the sample, granularity and morphology be systematically considered and reported along with measurements, which will in turn allow for more systematic data and predictive capabilities for sensitivity models.

Improved predictions of thermomechanical properties of molecular crystals from energy and dispersion corrected DFT

Thermal stability and pressure-dependent changes are key to molecular crystals and their properties. The determination of their thermal properties from ab initio methods is, however, a challenging task. While the low-frequency phonon spectrum related to intermolecular vibrations remains difficult to describe, the Quasi-Harmonic Approximation (QHA) also induces for molecular crystals a significant volume deviation, which makes their thermal behavior ill-determined. To overcome these difficulties, we consider a pragmatic energy correction (EC) that has long been used for atomic crystals, and we presently report the first ever use for molecular crystals. Applying the QHA in dispersion-corrected density functional theory (DFT-D) calculations with an ab initio parameterized EC, the resulting model can simultaneously and accurately derive thermal and mechanical properties of high-explosive molecular crystals. When compared to experiments, the mean absolute percent error of previous DFT-based thermomechanical models is 12% for mechanical and 31% for thermal properties. Our model performs significantly better and reduces these uncertainties to 4.1% and 9.8%, respectively. In particular, the agreement between our model and experiments for the thermal properties is three times better. This significant improvement greatly benefits the determination of thermomechanical properties such as the Grüneisen parameter and the shock properties. The method has been successfully applied to molecular crystals showing a large diversity of weak intermolecular interactions (β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), α-1,1-diamino-2,2-dinitroethylene (FOX-7), Triaminotrinitrobenzene (TATB), ε-Hexanitrohexaazaisowurtzitane (CL20), and Pentaerythritol tetranitrate (PETN)-I). Due to its accuracy and transferability, our model is expected to work for a large class of computationally designed molecular crystals and co-crystals, providing a basis for a predictive framework.

A flexible-molecule force field to model and study hexanitrohexaazaisowurtzitane (CL-20) – polymorphism under extreme conditions

The quantum-chemistry based force field (FF) developed for HMX by Smith and Bharadwaj (SB) [G. D. Smith and R. K. Bharadwaj, J. Phys. Chem. B, 1999, 103(18), 3570–3575] is transferred to another nitramine of different stoichiometry: hexanitrohexaazaisowurtzitane (CL-20 or HNIW). The modification of a single parameter alongside a very small number of add-ons related to carbon–carbon bonds, angles and dihedrals lead to two SB FF variants denoted SB-CL20 and SB-CL20 + CCNN. These flexible-molecule FFs should inherit the predictive capabilities of SB FF. For this purpose, we perform Molecular Dynamics simulations at ambient temperature and selected pressures. The modeled structures of the various CL-20 polymorphs are consistent with experimental data. Focusing on the ε-polymorph, we determine an equation of state which consolidates the general trend underpinned by most published results, and we confirm the increasing stiffness of the crystal under pressures up to 90 GPa. Moreover, we link some subtle pressure-induced changes of the elastic and structural properties to the flexibility and mobility of well-identified nitro groups. Finally, the simulations of the γ ↔ ζ phase transition suggest different multiple-step direct and reverse thermodynamic paths.

The quantum-chemistry based force field developed by Smith and Bharadwaj is transferred to hexanitrohexaazaisowurtzitane (CL20), revealing pressure-induced alterations of ε-CL20.

Granularity impact on hotspot formation and local chemistry in shocked nanostructured RDX

The nanostructuration of energetic materials results in interesting properties. In particular, the detonation of carbon-rich explosives leads to the formation of nanodiamonds, the size of which is linked to the initial size of the explosive nanograins. This correlation could come from the role played by the granularity in the shock properties, especially the local temperature, which could be enhanced when the shock front crosses the various interfaces and nanoporosities of the material. More generally, the granularity-dependent reactivity also concerns some aspects of the sensitivity of energetic materials and subsequent inhibition or safety issues. In order to investigate this hypothesis, we perform classical and reactive molecular dynamics shock-simulations on cyclotrimethylene trinitramine (RDX), a common high-explosive. We design various nanogranular structures using the body-centered-cubic stacking of nanoparticles shaped as Kelvin cells (truncated octahedra) with slightly-bumped faces. The nanograin size and the bump radius allow to independently control the porosity value and its nanostructure. We show that the global shock properties, especially the temperature, are sensitive to the porosity value but not to the nanograin size. The porosity nanostructure has a local impact, enhancing the temperature heterogeneities between the inter- and intra-grain regions by a few hundred degrees and slowing down the thermal homogenization. For a given porosity, the larger the nanograins, the larger and the hotter the hotspots. In such hotspots, the local chemistry is significantly modified, resulting in a larger reactivity with a quicker formation of some final products. We suggest that the quicker consumption of heteroatoms (namely, H, O, and N) along with higher local temperatures is likely to impact the formation process of solid carbonaceous phases.

Dynamic formation of nanodiamond precursors from the decomposition of carbon suboxide (C3O2) under extreme conditions-A ReaxFF study

We use molecular dynamics simulations with the ReaxFF-lg potential to model the high pressure pyrolysis of carbon suboxide (C₃O₂) in mixture with argon as a pressure bath. We show that the reactive simulations catch the experimental behavior of the low-pressure detonation of C₃O₂ (around 10 bars in shock tube experiments) and allow extrapolations to the high-pressure range of solid-state explosive detonation (up to 60 GPa). While at low pressure carbonaceous nanostructures are formed through the aggregation of species such as carbon dimers C₂, it appears that the high pressure deeply modifies the process, with the aggregation of growing C_xO_y heterostructures, in which the oxygen amount is driven by the pressure and the temperature. Pressures in the order of 60 GPa lead to high oxygen ratios, which prevent carbon atoms to get four carbon neighbors (the first condition to get a diamond structure). But a pressure lowering leads to a substantial carbon enrichment through CO₂/CO release and facilitates the formation of pure sp³-carbon phases where diamond precursors can form. These results give new insights on the conditions leading to nanodiamonds during the detonation of carbon-rich high explosives.