We investigated the heat of formation (Delta(f)H) of polynitrocubanes using density functional theory B3LYP and HF methods with 6-31G*, 6-311 + G**, and cc-pVDZ basis sets. The results indicate that Delta(f)H firstly decreases (nitro number m = 0-2) and then increases (m = 4-8) with each additional nitro group being introduced to the cubane skeleton. Delta(f)H of octanitrocubane is predicted to be 808.08 kJ mol(-1) at the B3LYP/6-311 + G** level. The Gibbs free energy of formation (Delta(f)G) increases by about 40-60 kJ mol(-1) with each nitro group being added to the cubane when the substituent number is fewer than 4, then Delta(f)G increases by about 100-110 kJ mol(-1) with each additional group being attached to the Cubic skeleton. Both the detonation velocity and the pressure for polynitrocubanes increase as the number of substituents increases. Detonation velocity and pressure of octanitrocubane are substantially larger than the famous widely used explosive cyclotetramethylenetetranitramine (HMX).