Finite temperature string by K-means clustering sampling with order parameters as collective variables for molecular crystals: application to polymorphic transformation between β-CL-20 and ε-CL-20

Polymorphic transformation of molecular crystals is a fundamental phase transition process, and it is important practically in the chemical, material, biopharmaceutical, and energy storage industries. However, understanding of the transformation mechanism at the molecular level is poor due to the extreme simulating challenges in enhanced sampling and formulating order parameters (OPs) as the collective variables that can distinguish polymorphs with quite similar and complicated structures so as to describe the reaction coordinate. In this work, two kinds of OPs for CL-20 were constructed by the bond distances, bond orientations and relative orientations. A K-means clustering algorithm based on the Euclidean distance and sample weight was used to smooth the initial finite temperature string (FTS), and the minimum free energy path connecting β-CL-20 and ε-CL-20 was sketched by the string method in collective variables, and the free energy profile along the path and the nucleation kinetics were obtained by Markovian milestoning with Voronoi tessellations. In comparison with the average-based sampling, the K-means clustering algorithm provided an improved convergence rate of FTS. The simulation of transformation was independent of OP types but was affected greatly by finite-size effects. A surface-mediated local nucleation mechanism was confirmed and the configuration located at the shoulder of potential of mean force, rather than overall maximum, was confirmed to be the critical nucleus formed by the cooperative effect of the intermolecular interactions. This work provides an effective way to explore the polymorphic transformation of caged molecular crystals at the molecular level.

High-pressure structural studies and pressure-induced sensitisation of 3,4,5-trinitro-1H-pyrazole

Herein we report the first high-pressure study of the energetic material 3,4,5-trinitro-1H-pyrazole (3,4,5-TNP) using neutron powder diffraction and single-crystal X-ray diffraction. A new high-pressure phase, termed Form II, was first identified through a substantial change in the neutron powder diffraction patterns recorded over the range 4.6-5.3 GPa, and was characterised further by compression of a single crystal to 5.3 GPa in a diamond-anvil cell using X-ray diffraction. 3,4,5-TNP was found to be sensitive to initiation under pressure, as demonstrated by its unexpected and violent decomposition at elevated pressures in successive powder diffraction experiments. Initiation coincided with the sluggish phase transition from Form I to Form II. Using a vibrational up-pumping model, its increased sensitivity under pressure can be explained by pressure-induced mode hardening. These findings have potential implications for the safe handling of 3,4,5-TNP, on the basis that shock- or pressure-loading may lead to significantly increased sensitivity to initiation.

Ab Initio Crystal Structure Prediction of the Energetic Materials LLM-105, RDX, and HMX

Crystal structure prediction (CSP) is performed for the energetic materials (EMs) LLM-105 and α-RDX, as well as the α and β conformational polymorphs of 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), using the genetic algorithm (GA) code, GAtor, and its associated random structure generator, Genarris. Genarris and GAtor successfully generate the experimental structures of all targets. GAtor's symmetric crossover scheme, where the space group symmetries of parent structures are treated as genes inherited by offspring, is found to be particularly effective. However, conducting several GA runs with different settings is still important for achieving diverse samplings of the potential energy surface. For LLM-105 and α-RDX, the experimental structure is ranked as the most stable, with all of the dispersion-inclusive density functional theory (DFT) methods used here. For HMX, the α form was persistently ranked as more stable than the β form, in contrast to experimental observations, even when correcting for vibrational contributions and thermal expansion. This may be attributed to insufficient accuracy of dispersion-inclusive DFT methods or to kinetic effects not considered here. In general, the ranking of some putative structures is found to be sensitive to the choice of the DFT functional and the dispersion method. For LLM-105, GAtor generates a putative structure with a layered packing motif, which is desirable thanks to its correlation with low sensitivity. Our results demonstrate that CSP is a useful tool for studying the ubiquitous polymorphism of EMs and shows promise of becoming an integral part of the EM development pipeline.

Study on the Cocrystallization Mechanism of CL-20/HMX in a Propellant Aging Process through Theoretical Calculations and Experiments

Energetic materials undergo physical and chemical aging due to environmental effects, resulting in the degradation of safety and detonation performances. Therefore, studying the aging performance of energetic materials is of great importance for the efficient application of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20)-based solid propellants. In this paper, XRD and FTIR of the CL-20-based propellant and CL-20/1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX)-based propellant samples showed CL-20/HMX cocrystal formation according to appearance of new peaks. SEM and EDS analyses showed that pores and dehumidification in the propellant occurred with the cocrystallization of CL-20 and HMX during the aging process. Furthermore, molecular dynamics simulation was used to predict the crystal transformation of the CL-20- and HMX-based propellant under a long-term storage process. The stability of ε-CL-20 was obtained by analyzing the crystal transformation rate. The binding energy, radial distribution function between CL-20 and HMX, as well as mechanical properties of the CL-20/HMX cocrystal and the mixture were calculated to reveal the stronger binding between CL-20 and HMX in the cocrystal. Meanwhile, the inducer effect of a nitrate ester during the cocrystallization process was analyzed. The theoretical calculation shows that during aging, ε-CL-20 tends to exist stably, while CL-20/HMX tends to form cocrystals because of the strong bond. The present work on the transformation and cocrystallization of CL-20 and HMX during long-term storage is beneficial for understanding the degradation mechanism of the propellant performances, facilitating safe storage and life evaluation of propellants.

Strategies for Achieving Balance between Detonation Performance and Crystal Stability of High-Energy-Density Materials

Performance-stability contradiction of high-energy-density materials (HEDMs) is a long-standing puzzle in the field of chemistry and material science. Bridging the gap that exists between detonation performance of new HEDMs and their stability remains a formidable challenge. Achieving optimal balance between the two contradictory factors is of a significant demand for deep-well oil and gas drilling, space exploration, and other civil and defense applications. Herein, supercomputers and latest quantitative computational strategies were employed and high-throughput quantum calculations were conducted for 67 reported HEDMs. Based on statistical analysis of large amounts of physico-chemical data, in-crystal interspecies interactions were identified to be the one that provokes the performance-stability contradiction of HEDMs. To design new HEDMs with both good detonation performance and high stability, the proposed systematic and comprehensive strategies must be satisfied, which could promote the development of crystal engineering of HEDMs to an era of theory-guided rational design of materials.

Micronization of CL-20 using supercritical and liquefied gases

Liquid and supercritical CO₂ and TFE have been systematically studied as media for CL-20 recrystallization to estimate their potential for the selective preparation of ultrafine particles.

Modification of crystalline energetic salts through polymorphic transition: enhanced crystal density and energy performance

We presented a detailed investigation of polymorphic transition of energetic salts and explored a new path for modifying crystalline energetic salts.

A Study of the Shock Sensitivity of Energetic Single Crystals by Large-Scale Ab Initio Molecular Dynamics Simulations

Understanding the reaction initiation of energetic single crystals under external stimuli is a long-term challenge in the field of high energy density materials. Herewith, we developed an ab initio molecular dynamics method based on the multiscale shock technique (MSST) and reported the reaction initiation mechanism by performing large-scale simulations for the sensitive explosive benzotrifuroxan (BTF), insensitive explosive triaminotrinitrobenzene (TATB), four polymorphs of hexanitrohexaazaisowurtzitane (CL-20) pristine crystals and five novel CL-20 cocrystals. A theoretical indicator, t_initiation, the delay of decomposition reaction under shock, was proposed to characterize the shock sensitivity of energetic single crystal, which was proved to be reliable and satisfactorily consistent with experiments. We found that it was the coupling of heat and pressure that drove the shock reaction, wherein the vibrational spectra, the specific heat capacity, as well as the strength of the trigger bonds were the determinants of the shock sensitivity. The intermolecular hydrogen bonds were found to effectively buffer the system from heating, thereby delaying the decomposition reaction and reducing the shock sensitivity of the energetic single crystal. Theoretical rules for synthesizing novel energetic materials with low shock sensitivity were given. Our work is expected to provide a useful reference for the understanding, certifying and adjusting of the shock sensitivity of novel energetic materials.

Design, preparation, characterization and formation mechanism of a novel kinetic CL-20-based cocrystal

2,4,6,8,10,12-Hexanitrohexaazaisowurtzitane (CL-20)-based cocrystals have gained increasing attention as a means of obtaining insensitive high explosives. However, the design of ideal candidates for these cocrystals remains difficult. This work compares the crystal energies of the CL-20-dinitrobenzene (DNB) and CL-20-2,4,6-trinitrotoluene (TNT) cocrystals with those of the respective pure coformers. The results indicate that the cocrystal formation is driven by the differences in the energies of the cocrystals and the coformers. Furthermore, analysis via Hirshfeld surfaces and two-dimensional fingerprint plots confirms that the O...O, O...H, O...N and C...O interactions were the main force for stabilizing the CL-20-based cocrystal structure. Based on these findings, a novel energetic-energetic cocrystal of CL-20-2,4,6-trinitrophenol (TNP) was designed and prepared by means of a rapid method for solvent removal. The crystal structure was investigated via powder X-ray diffraction methods, solid-state nuclear magnetic resonance spectroscopy and Fourier transform infrared spectroscopy. The results revealed that the O-H...O hydrogen bonding interaction between the phenolic hydroxyl group of TNP and nitro groups of CL-20, as well as nitro...π, nitro...nitro and O_NO2...π(N)_NO2 interactions, based on the benzene ring and nitro groups, are the main interactions occurring in the cocrystal.

Host–guest energetic materials constructed by incorporating oxidizing gas molecules into an organic lattice cavity toward achieving highly-energetic and low-sensitivity performance

Powerful oxidizer N₂O was incorporated into an organic lattice cavity through aeration crystallization, and smart host–guest energetic materials with highly-energetic and low-sensitivity performance were obtained.

Boosting the performance of energetic materials through thermally-induced conformational transition

We presented an effective strategy to improve the performance of energetic materials through thermally-induced conformational transition.

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.

Discovery and recovery of delta p-aminobenzoic acid

A new high-pressure recoverable form has been observed in the model system, p-aminobenzoic acid.

High-pressure crystallography of periodic and aperiodic crystals

More than five decades have passed since the first single-crystal X-ray diffraction experiments at high pressure were performed. These studies were applied historically to geochemical processes occurring in the Earth and other planets, but high-pressure crystallography has spread across different fields of science including chemistry, physics, biology, materials science and pharmacy. With each passing year, high-pressure studies have become more precise and comprehensive because of the development of instrumentation and software, and the systems investigated have also become more complicated. Starting with crystals of simple minerals and inorganic compounds, the interests of researchers have shifted to complicated metal-organic frameworks, aperiodic crystals and quasicrystals, molecular crystals, and even proteins and viruses. Inspired by contributions to the microsymposium 'High-Pressure Crystallography of Periodic and Aperiodic Crystals' presented at the 23rd IUCr Congress and General Assembly, the authors have tried to summarize certain recent results of single-crystal studies of molecular and aperiodic structures under high pressure. While the selected contributions do not cover the whole spectrum of high-pressure research, they demonstrate the broad diversity of novel and fascinating results and may awaken the reader's interest in this topic.

The crystal structure and morphology of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) p-xylene solvate: a joint experimental and simulation study

The crystal structure of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaiso-wurtzitane (CL-20) p-xylene solvate, and the solvent effects on the crystal faces of CL-20 were studied through a combined experimental and theoretical method. The properties were analyzed by thermogravimetry-differential scanning calorimetry (TG-DSC), Fourier transform infrared spectroscopy (FTIR) and X-ray diffraction (XRD).The growth morphology of CL-20p-xylene solvate crystal was predicted with a modified attachment energy model. The crystal structure of CL-20p-xylene solvate belonged to the Pbca space group with the unit cell parameters, a = 8.0704(12) Å, b=13.4095(20) Å, c = 33.0817(49) Å, and Z = 4, which indicated that the p-xylene solvent molecules could enter the crystal lattice of CL-20 and thus the CL-20 p-xylene solvate is formed. According to the solvent-effected attachment energy calculations, (002) and (11−1) faces should not be visible at all, while the percentage area of the (011) face could be increased from 7.81% in vacuum to 12.51% in p-xylene solution. The predicted results from the modified attachment energy model agreed very well with the observed morphology of crystals grown from p-xylene solution.