Quantum chemical calculations are used to compute the heats of formation (HOFs) for 24 nitroester (NE) energetic compounds in which only 5 nitroester energetic compounds have the available experimental heats of formation. The heats of formation of the five compounds are calculated from isodesmic reactions by employing the hybrid density functional theory (DFT) (B3LYP, B3PW91, and B3P86) methods with 6-31G** and 6-311G** basis sets. It is demonstrated that B3PW91/6-31G** method can yield reliable HOFs, which has the mean absolute deviation of 1.1 kcal/mol. The HOF's of other 19 nitroester energetic compounds are calculated by using B3PW91/6-31G** method. Through the analysis of the Calculated result, it is found that the computed heat of formation decreases when the number of methylene (CH2) group increases for normal chain nitroester compounds. The further study shows that Our results about gas-phase heats of formation of nitroester compounds are better than the results of Muthurajan et al. [H. Muthurajan, R. Sivabalan, M.B. Talawar, M. Anniyappan, S. Venugopalan, Prediction of heat of formation and related parameters of high energy materials, J. Hazard. Mater. A133 (2006) 30-45], and Byrd Edward and Rice Betsy [F.C. Byrd Edward, M. Rice Betsy, Improved prediction of heats of formation of energetic materials using quantum mechanical calculations. J. Phys. Chem. A 110 (2006) 1005-1013]. In addition, the condensed phase heat of formation of the nitroester compounds are computed through the same method of Byrd Edward and Rice Betsy. (C) 2008 Elsevier B.V. All rights reserved.