s-Tetrazine is the essential candidate of many energetic compounds due to its high nitrogen content, enthalpy of formation and thermal stability. The present study explores the design of s-tetrazine derivatives in which different $-NO_2$, $-NH_2$ and $-N_3$ substituted azoles are attached to the tetrazine ring via C-N linkage. The density functional theory (DFT) is used to predict the geometries, heats of formation (HOFs) and other energetic properties. The predicted results show that azide group plays a very important role in increasing HOF values of the s-tetrazine derivatives. The densities for designed molecules were predicted by using the crystal packing calculations. The introduction of $-NO_2$ group improves the density as compared to $-N_3$, and $-NH_2$ groups and hence the detonation performance. Bond dissociation energy analysis and insensitivity correlations revealed that amino derivatives are better candidates considering insensitivity and stability.