The SCF calculation produces a density, P, and Fock
matrix, F. These, together with the one-electron matrix, H,
allow the total electronic energy to be calculated via
A more useful quantity is the heat of formation of the compound from its elements in their standard state. This is obtained when the energy required to ionize the valence electrons of the atoms involved (calculated using semiempirical parameters), Eisol(A), and heat of formation for a gaseous atom from its standard state, Eatom(A), are added to the electronic plus nuclear energy. This yields:
ΔHf = (Eelect + Enuc) - ΣAEisol(A) + ΣAEatom(A),
or
ΔHf = Etot
- ΣAEisol(A)
+ ΣAEatom(A).
This is the quantity which MOPAC calls the "Heat of Formation". An alternative but equivalent definition of ΔHf , more suited for comparison with experimental ΔHf, is:
" ΔHf is the calculated gas-phase heat of formation at 298K of one mole of a compound from its elements in their standard state."
Things to note about this definition: unlike ab initio methods, which yield the energy at 0K, semiempirical methods give ΔHf at 298K. This follows from the way in which semiempirical methods are parameterized: the reference ΔHf are conventionally given at 298K. This means that semiempirical methods will reproduce ΔHf for 298K. Secondly, note that ΔHf are for gas-phase systems. To calculate ΔHf in the liquid or solid phases, additional terms are necessary.