Prediction of enthalpies of sublimation of high-nitrogen energetic compounds: Modified Politzer model

Enthalpy of formation of 6-phenyl-1,5-diazabicyclo[3.1.0]hexane by combustion calorimetry and theoretical approach for efficient prediction of thermochemistry of diaziridines

The combustion energy and standard molar enthalpy of formation of crystalline 6-phenyl-1,5-diazabicyclo[3.1.0]hexane (PDABH) were determined using an isoperibolic calorimeter with a static bomb. PDABH is the first diaziridine for which the experimental value of the enthalpy of formation was obtained. This value was validated by the theoretical values of gas phase enthalpy of formation and enthalpy of sublimation. The gas phase enthalpy of formation was calculated using the DLPNO-CCSD(T1)/CBS method in conjunction with isodesmic-type reactions. This method was chosen in comparison to another high quality evaluative method (G4), which has been shown to provide unreliable results for cyclic nitrogen containing compounds. The descriptors of the molecular electrostatic potential (MEP) were used to estimate the enthalpy of sublimation of PDABH. The proposed MEP model is based on experimental enthalpies of sublimation for 75 compounds structurally similar to PDABH. The high-level ab initio calculations of gas phase enthalpies of formation combined with enthalpies of sublimations estimated using descriptors of MEP allow predicting the enthalpies of formation of diaziridines in the solid phase.

Performance of Dispersion-Inclusive Density Functional Theory Methods for Energetic Materials

Molecular crystals of energetic materials (EMs) are denser than typical molecular crystals and are characterized by distinct intermolecular interactions between nitrogen-containing moieties. To assess the performance of dispersion-inclusive density functional theory (DFT) methods, we have compiled a data set of experimental sublimation enthalpies of 31 energetic materials. We evaluate the performance of three methods: the semilocal Perdew-Burke-Ernzerhof (PBE) functional coupled with the pairwise Tkatchenko-Scheffler (TS) dispersion correction, PBE with the many-body dispersion (MBD) method, and the PBE-based hybrid functional (PBE0) with MBD. Zero-point energy contributions and thermal effects are described using the quasi-harmonic approximation (QHA), including explicit treatment of thermal expansion, which we find to be non-negligible for EMs. The lattice energies obtained with PBE0+MBD are the closest to experimental sublimation enthalpies with a mean absolute error of 9.89 kJ/mol. However, the state-of-the-art treatment of vibrational and thermal contributions makes the agreement with experiment worse. Pressure-volume curves are also examined for six representative materials. For pressure-volume curves, all three methods provide reasonable agreement with experimental data with mean absolute relative errors of 3% or less. Most of the intermolecular interactions typical of EMs, namely nitro-amine, nitro-nitro, and nitro-hydrogen interactions, are more sensitive to the choice of the dispersion method than to the choice of the exchange-correlation functional. The exception is π-π stacking interactions, which are also very sensitive to the choice of the functional. Overall, we find that PBE+TS, PBE+MBD, and PBE0+MBD do not perform as well for energetic materials as previously reported for other classes of molecular crystals. This highlights the importance of testing dispersion-inclusive DFT methods for diverse classes of materials and the need for further method development.

Synthetic Strategies Toward Nitrogen-Rich Energetic Compounds Via the Reaction Characteristics of Cyanofurazan/Furoxan

The structural units of amino-/cyano-substituted furazans and furoxans played significant roles in the synthesis of nitrogen-rich energetic compounds. This account focused on the synthetic strategies toward nitrogen-rich energetic compounds through the transformations based on cyanofurazan/furoxan structures, including 3-amino-4-cyanofurazan, 4-amino-3-cyano furoxan, 3,4-dicyanofurazan, and 3,4-dicyanofuroxan. The synthetic strategies toward seven kinds of nitrogen-rich energetic compounds, such as azo (azoxy)-bridged, ether-bridged, methylene-bridged, hybrid furazan/furoxan-tetrazole–based, tandem furoxan–based, hybrid furazan-isofurazan–based, hybrid furoxan-isoxazole–based and fused framework–based energetic compounds were fully reviewed, with the corresponding reaction mechanisms toward the nitrogen-rich aromatic frameworks and examples of using the frameworks to create high energetic substances highlighted and discussed. The energetic properties of typical nitrogen-rich energetic compounds had also been compared and summarized.

Learning to fly: thermochemistry of energetic materials by modified thermogravimetric analysis and highly accurate quantum chemical calculations

The standard state enthalpy of formation and the enthalpy of sublimation are essential thermochemical parameters determining the performance and application prospects of energetic materials (EM). Direct experimental measurements of these properties are complicated by low volatility and high heat release in bomb calorimetry experiments. As a result, the uncertainties in the reported enthalpies of formation for a number of even well-known CHNO-containing compounds might amount up to tens kJ mol-1, while for some novel high-nitrogen molecules they reach even hundreds of kJ mol-1. The present study reports a facile approach to determining the solid-state formation enthalpies comprised of complementary high-level quantum chemical calculations of the gas-phase thermochemistry and advanced thermal analysis techniques yielding sublimation enthalpies. The thermogravimetric procedure for the measurement of sublimation enthalpy was modified by using low external pressures (down to 0.2 Pa). This allows for observing sublimation/vaporization instead of thermal decomposition of the compounds studied. Extensive benchmarking on nonenergetic and energetic compounds reveals the average and maximal absolute errors of the sublimation enthalpies of 3.3 and 11.0 kJ mol-1, respectively. The comparison of the results with those obtained from the widely used Trouton-Williams empirical equation shows that the latter underestimates the sublimation enthalpy up to 140 kJ mol-1. Therefore, we performed a reparametrization of the latter equation with simple chemical descriptors that reduces the mean error down to 30 kJ mol-1. Highly accurate multi-level procedures W2-F12 and/or W1-F12 in conjunction with the atomization energy approach were used to calculate theoretically the gas-phase formation enthalpies. In several cases, the DLPNO-CCSD(T) enthalpies of isodesmic reactions were also employed to obtain the gas-phase thermochemistry for medium-sized important EMs. Combining the obtained thermochemical properties, we determined the solid-state enthalpies of formation for nearly 60 species containing various important explosophoric groups, from common nitroaromatics, nitroethers, and nitramines to novel nitrogen-rich heterocyclic species (e.g., the derivatives of pyrazole, tetrazole, furoxan, etc.). The large-scale benchmarking against the available experimental solid-state enthalpies of formation yielded the maximal inaccuracy of the proposed method of 25 kJ mol-1.

Thermal and Sensitiveness Determination of Cubanes: Towards Cubane-Based Fuels for Infrared Countermeasures

As infrared seeking technology evolves, threats are better able to distinguish defensive infrared (IR) flares from true targets. Spectrally matched flares, which generally employ carbon-based fuels, are better able to decoy some advanced missiles by more closely mimicking the IR emission of the target. Cubane is a high-energy carbon-based scaffold which may be suitable for use as a fuel in spectrally matched flares. The enthalpy of formation and strain energy of a series of cubanes was predicted in silico, and their thermal and impact stability examined. All were found to undergo highly exothermic decomposition in sealed cell differential scanning calorimetry, and two cubanes subsequently underwent quantitative sensitiveness testing. Despite their F of I values being in the secondary explosive range, cubane-1,4-dicarboxylic acid (F of I=70) and 4-carbamoylcubane-1-carboxylic acid (F of I=90) were identified as potentially useful fuels for pyrotechnic infrared countermeasure flare formulations.

High pressure behavior of crystal [2,2'-bi(1,3,4-oxadiazole)]-5,5'-dinitramide: A DFT investigation

Density functional theory (DFT) computation was carried out to investigate the crystal, molecular and electronic structures of high energy crystal [2,2'-bi(1,3,4-oxadiazole)]-5,5'-dinitramide (BODN) with the pressure 0-120 GPa. The relaxed crystal structure by the GGA/PBE-TS functional matches well with the experimental data at ambient pressure condition. With the intensifying of pressure, the lattice parameters, volumes, bond lengths, H-bond energies, atomic charges, bond populations, band gaps and density of states of crystal BODN change gently. Under the pressure of 48, 104, and 107 GPa, three pressure-induced transformations occurred. The intramolecular six membered rings pose strong affect in stabilizing systems in the pressure range 0-120 GPa. Between O1 and H2 atoms, the H-bond interaction transforms into covalent interaction under the circumstance of 48 GPa. At 104 GPa, structural transformation occurs with the distortion of the intramolecular six membered ring. In addition, O1⋅⋅⋅H2 and O2⋅⋅⋅H1 have the largest H-bond energies in comparison with the others. When the pressure reaches 107 GPa, the H-bond O1⋅⋅⋅H2 is formed again with the deformation and non-coplanarity of two oxadiazoles in crystal BODN. The electrons can be moved easily based on the density of states and energy bands under high pressure. Helpful information will be conveyed by this work in the field of further analysis connected the pressure effect on molecular transformations.