Energy Decomposition of the Intermolecular Interactions in Energetic Co-crystals
Energetic co-crystals (ECCs) are now thriving and becoming alternatives to energetic materials. Thereby, it is important to understand the intermolecular interactions to obtain the knowledge for ECC engineering. However, the physical sources of the interactions remain unclear, even though the interactions have already been understood as the three conventional basic interaction kinds, or the three main traditional engineering motifs of organic crystals, including hydrogen bonding, π-stacking and halogen bonding. Twelve typical molecular pairs extracted from five observed EECs covering all the three interaction kinds are selected to partition intermolecular interaction energy and discuss the physical sources, by density functional theory calculations and the block located wavefunction energy decomposition analyses (BLW-EDA). We find that, after carefully examining all observed ECCs, each conventional interaction motif in energetic-energetic molecular pairs is always weak, and dominated by a frozen effect, i.e., van de Waals and electrostatic interactions. The rather strong hydrogen bonding exists in the molecular pairs with one non-energetic molecule at least, and is predominated by the polarization and charge transfer effects. Meanwhile, we find small bond order variations caused by the crystal packing of energetic-energetic co-crystals (EECCs), thereby showing small molecular stability variation. It suggests that it is difficult to increase the molecular stability of energetics by cocrystallization to improve safety; while the safety of EECCs will benefit from the enhanced intermolecular interactions and the improved crystal packing mode favoring ready shear slip.