The nitro group was introduced into the nitrogen heterocycle of 1H-benzotriazole to design a total of 31 derivatives. To estimate the thermal stability of these derivatives, the heat of formation (HOF) is calculated based on the isodesmic reaction. The bond dissociation energy (BDE) was also predicted based on the homolytic reaction to further evaluate the dynamic stability. To evaluate the possibility of utilizing as high energy density compounds (HEDCs), the detonation parameters including the detonation pressure (P), detonation velocity (D), and explosive heat (Q) are predicted by taking advantage of the Kamlet-Jacobs empirical equation. To measure the sensitivity to impact, both the characteristic height (H50) and free space in crystal (∆V) are considered in this paper. Based on our calculations, D-series and E are found to be the candidates for HEDCs. The Gaussian 09 software package was used in this paper. The B3PW91 hybrid function with the 6-311 + G(d,p) basis set was chosen to perform the structural optimization, frequency analysis, heat of formation, and bond dissociation energy. The detonation parameters were calculated following the Kamlet-Jacobs equation.