The abundant dataset of detonation parameters for energetic materials is reported, empirical and thermodynamic calculation methods are benchmarked, and implications for design of novel promising compounds provided.

Thermal decomposition of a novel promising high-performance explosive dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate (TKX-50) was studied using a number of thermal analysis techniques (thermogravimetry, differential scanning calorimetry, and accelerating rate calorimetry, ARC). To obtain more comprehensive insight into the kinetics and mechanism of TKX-50 decomposition, a variety of complementary thermoanalytical experiments were performed under various conditions. Non-isothermal and isothermal kinetics were obtained at both atmospheric and low (up to 0.3 Torr) pressures. The gas products of thermolysis were detected in situ using IR spectroscopy, and the structure of solid-state decomposition products was determined by X-ray diffraction and scanning electron microscopy. Diammonium 5,5′-bistetrazole-1,1′-diolate (ABTOX) was directly identified to be the most important intermediate of the decomposition process. The important role of bistetrazole diol (BTO) in the mechanism of TKX-50 decomposition was also rationalized by thermolysis experiments with mixtures of TKX-50 and BTO. Several widely used thermoanalytical data processing techniques (Kissinger, isoconversional, formal kinetic approaches, etc.) were independently benchmarked against the ARC data, which are more germane to the real storage and application conditions of energetic materials. Our study revealed that none of the Arrhenius parameters reported before can properly describe the complex two-stage decomposition process of TKX-50. In contrast, we showed the superior performance of the isoconversional methods combined with isothermal measurements, which yielded the most reliable kinetic parameters of TKX-50 thermolysis. In contrast with the existing reports, the thermal stability of TKX-50 was determined in the ARC experiments to be lower than that of hexogen, but close to that of hexanitrohexaazaisowurtzitane (CL-20).

Thermal decomposition of solids often includes simultaneous occurrence of the overlapping processes with unequal contributions in a specific physical quantity variation during the overall reaction (e.g., the opposite effects of decomposition and evaporation on the caloric signal). Kinetic analysis for such reactions is not a straightforward, while the applicability of common kinetic calculation methods to the particular complex processes has to be justified. This study focused on the critical analysis of the available kinetic calculation methods applied to the mathematically simulated thermogravimetry (TG) and differential scanning calorimetry (DSC) data. Comparing the calculated kinetic parameters with true kinetic parameters (used to simulate the thermoanalytical curves), some caveats in the application of the Kissinger, isoconversional Friedman, Vyazovkin and Flynn–Wall–Ozawa methods, mathematical and kinetic deconvolution approaches and formal kinetic description were highlighted. The model-fitting approach using simultaneously TG and DSC data was found to be the most useful for the complex processes assumed in the study.