Substituent effects on the thermal stability and initial decomposition of methylene bis(dinitropyrazole) derivatives

Elucidating the role of functional-group substitution in the thermal stability of energetic materials is of considerable importance for the rational development of advanced energetic compounds. In this study, density functional theory (DFT) calculations, together with electronic-structure, weak-interaction, crystal-packing, and transition-state analyses, were performed to systematically explore the influence of amino (A), hydroxyl (B), and azido (C) substituents on methylene-bridged bis(3,5-dinitropyrazole) (MBDNP) derivatives. The results reveal that functional-group substitution induces distinct electronic redistribution, leading to significant differences in key bond stability, intra- and intermolecular interactions, and crystal packing. Among the three derivatives, A exhibits the strongest electron delocalization, the most favorable interaction network, and the strongest intermolecular interactions, whereas C shows the weakest molecular and crystal stability. Transition-state calculations further demonstrate that different substituents alter the preferred initial decomposition pathways, and the calculated activation free energies follow the order A > B > C, consistent with the experimentally observed decomposition temperatures. These findings demonstrate that the thermal stability of MBDNP derivatives is governed by the synergistic effects of electronic structure, molecular interactions, crystal packing, and initial decomposition kinetics, providing molecular-level insights into the substituent-dependent thermal decomposition mechanism of energetic materials.

Density functional theory (DFT) calculations at the M06-2X/def2-TZVPP level were performed to investigate the electronic structures and initial decomposition behavior of compounds A-C. Natural bond orbital (NBO), electrostatic potential (ESP), localized bond order (LBO), interaction region indicator (IRI), and quantum theory of atoms in molecules (QTAIM) analyses, together with Hirshfeld surface, energy framework, and transition-state analyses, were further performed to clarify how different substituents influence the electronic characteristics, intermolecular interactions, crystal-packing behavior, and thermal stability of the studied compounds.