The five-membered fused six-membered nitrogen heteroaromatic ring system is a crucial skeleton in the design and synthesis of energetic compounds. Based on this skeleton, many high-performance energetic compounds have been synthesized. However, to date, no one has conducted a systematic study on the characteristics of this skeleton itself. To assess how the number and position of nitrogen atoms affect the energy and stability of this type of skeleton, one to four nitrogen-substituted skeleton molecules were analyzed using Density Functional Theory (DFT) calculations. Natural population analysis (NPA), Laplacian bond order (LBO) analysis, aromaticity studies, and enthalpy of formation calculations were performed. Patterns observed in the computational results were summarized, and their potential correlations were analyzed. Based on these findings, design recommendations for derivatives of these skeletons in energetic compounds were proposed to serve as a reference for energetic material chemists.