Cage strain in nitro-substituted 1,3,5,7-tetraazacubanes

Homodesmotic reactions were designed for the computation of strain energies (SE) for four nitro-substituted 1,3,5,7-tetraazacubane derivatives. Total energies of the optimized geometric structures at the DFT-B3LYP/6-31G* and DFT-B3LYP/6-311G** levels were used to derive the SE. The variation of SE with respect to the number of substituents is similar with both basis sets. The SE value is 237.32 kcal/mol at the B3LYP/6-311G** level for 2,4,6,8-tetranitro-1,3,5,7-tetraazacubane, which is unexpectedly much larger than that of its cubane analogue. The SE increases remarkably with more nitro groups being attached to the cage skeleton of tetraazacubane. The ‘bending’ of the bonds within the cubic skeleton attributes to the increase of strains as the attached number of nitro groups increases.