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⁻¹, while for some novel high-nitrogen molecules they reach even hundreds of kJ mol⁻¹. 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⁻¹, 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⁻¹. Therefore, we performed a reparametrization of the latter equation with simple chemical descriptors that reduces the mean error down to 30 kJ mol⁻¹. 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⁻¹.

In the present work, we studied in detail the thermochemistry, thermal stability, mechanical sensitivity, and detonation performance for 20 nitro-, cyano-, and methyl derivatives of 1,2,5-oxadiazole-2-oxide (furoxan), along with their bis-derivatives. For all species studied, we also determined the reliable values of the gas-phase formation enthalpies using highly accurate multilevel procedures W2-F12 and/or W1-F12 in conjunction with the atomization energy approach and isodesmic reactions with the domain-based local pair natural orbital (DLPNO) modifications of the coupled-cluster techniques. Apart from this, we proposed reliable benchmark values of the formation enthalpies of furoxan and a number of its (azo)bis-derivatives. Additionally, we reported the previously unknown crystal structure of 3-cyano-4-nitrofuroxan. Among the monocyclic compounds, 3-nitro-4-cyclopropyl and dicyano derivatives of furoxan outperformed trinitrotoluene, a benchmark melt-cast explosive, exhibited decent thermal stability (decomposition temperature >200 °C) and insensitivity to mechanical stimuli while having notable volatility and low melting points. In turn, 4,4′-azobis-dicarbamoyl furoxan is proposed as a substitute of pentaerythritol tetranitrate, a benchmark brisant high explosive. Finally, the application prospects of 3,3′-azobis-dinitro furoxan, one of the most powerful energetic materials synthesized up to date, are limited due to the tremendously high mechanical sensitivity of this compound. Overall, the investigated derivatives of furoxan comprise multipurpose green energetic materials, including primary, secondary, melt-cast, low-sensitive explosives, and an energetic liquid.

A reliable kinetic description of the thermal stability of energetic materials (EM) is very important for safety and storage-related problems. Among other pertinent issues, autocatalysis very often complicates the decomposition kinetics of EM. In the present study, the kinetics and decomposition mechanism of a promising energetic compound, 5-amino-3,4-dinitro-1H-pyrazole (5-ADP) were studied using a set of complementary experimental (e.g., differential scanning calorimetry in the solid state, melt, and solution along with advanced thermokinetic models, accelerating rate calorimetry, and evolved gas analysis) and theoretical techniques (CCSD(T)-F12 and DLPNO-CCSD(T) predictive quantum chemical calculations). The experimental study revealed that the strong acceleration of the decomposition rate of 5-ADP is caused by two factors: the progressive liquefaction of the sample directly observed using in situ optical microscopy, and the autocatalysis by reaction products. For the first time, the processing of the non-isothermal data was performed with a formal Manelis–Dubovitsky kinetic model that accounts for both factors. With the aid of quantum chemical calculations, we have rationalized the autocatalysis present in the formal kinetic models at the molecular level. Theory revealed an unusual primary decomposition channel of 5-ADP, viz., the two subsequent sigmatropic H-shifts in the pyrazole ring followed by the C–NO₂ bond scission yielding a pyrazolyl and nitrogen dioxide radicals as simple primary products. Moreover, we found the secondary reactions of the latter radical with the 5-ADP to be kinetically unimportant. On the contrary, the substituted pyrazolyl radical turned out to undergo a facile addition to 5-ADP, followed by a fast exothermic elimination of another ˙NO₂ species. We believe the latter process to contribute remarkably to the observed autocatalytic behavior of 5-ADP. Most importantly, the calculations provide detailed mechanistic evidence complementing the thermoanalytical experiment and formal kinetic models.

Novel energetic materials (EM) often combine two intrinsically counter trends, viz., a high energy density and mediocre safety parameters, like thermal stability and sensitivity toward mechanical stimuli. A rational design of promising EMs requires a proper understanding of their thermal stability at both macroscopic and molecular levels. In the present contribution, we studied in detail the thermal stability of 4,4′-dinitro-3,3′-diazenofuroxan (DDF), an ultrahigh-performance energetic material with a reliable experimental detonation velocity being very close to 10 km s⁻¹. To this end, we employed a set of complementary thermoanalytical (DSC and TGA in the solid state along with advanced thermokinetic models, optical microscopy, and gas products detection) and theoretical techniques (DLPNO-CCSD(T) quantum chemical calculations). According to the DSC measurements, the solid-state thermolysis of DDF turned out to be a complex three-step process. The decomposition commences at ∼85 °C and the most intense heat release occurs at ∼130 °C depending on the heating rate. In order to properly describe the kinetics of DDF thermolysis beyond the simple Kissinger and Friedman methods, we applied a “top-down” kinetic approach resulting in the formal model comprised of three independent stages. A flexible Kolmogorov–Johnson–Mehl–Avrami–Erofeev equation was applied for the first decomposition stage along with the extended Prout–Tompkins equation for the second and third processes, respectively. The formal exponent in the former equation turned out to be close to a second order, thus suggesting a two-dimensional nuclei-growth model for the first stage. We rationalized this fact with the aid of optical microscopy experiments tracking the changes in the morphology of a solid DDF sample. Then, we complemented the formal macroscopic kinetics with some mechanistic patterns of the primary decomposition channels from quantum chemical calculations. The three reactions involving all important moieties of the DDF molecule turned out to compete very closely: viz., the nitro-nitrite isomerization, radical C(heterocycle)–N(bridge) bond scission and molecular decomposition comprised of the consequent N–O and C–C bond scissions in a furoxane ring. The DLPNO-CCSD(T) activation barriers of all these reactions were close to ∼230 kJ mol⁻¹. Most importantly, the calculations provide some mechanistic details missing in thermoanalytical experiment and formal kinetic models. Apart from this, we also determined a mutually consistent set of thermochemical and phase change data for DDF.