Theoretical study on intra-molecule interactions in TKX-50

Abstract

To fully and deeply understand the weak interactions in the gaseous structure of the TKX-50 molecule, two conformations I and II of the TKX-50 molecule confirmed in a crystal cell were optimized at the B3LYP/6-311g(d,p) level in the gas state, and the single point energy of the optimized structure was calculated at the M06-2X/ma-TZVPP level. Analyzing methods for weak interactions such as the interaction region indicator (IRI), topological basin analysis, and the extended transition state-natural orbitals for chemical valence (ETS-NOCV) theory with the help of Multiwfn code were employed to reveal the corresponding intramolecular weak interactions. The results showed that there were 5 kinds of intramolecular weak interaction in both conformations. They are two types of H bond, two types of intra-ring weak interaction, and one type of O–N bond within the molecular fragment containing the bis-tetrazole ring. The combined effect of all these weak interactions holds the bis-tetrazole ring of TKX-50 retaining an almost coplanar configuration. Meanwhile, the strength of these weak interactions is significantly different in conformation I and conformation II. The most obvious difference is that conformation II has a significant H transfer between intramolecular fragments due to the mirror rotation of almost 180° of cations (NH₃OH)⁺ perpendicular to the N–O bond axis thereof as compared to the reference conformation I. This conformational difference not only makes the weak interaction between the two conformations very different but also forms a quasi-covalent bond in conformation II with much larger bonding energy than other H bonds, thus resulting in conformation II having lower electron energy and more stable geometry. In addition, the order of breaking various H bonds in the combustion decomposition process of TKX-50 is deduced by comparing various H bonds.