MOLECULAR-BOND-ENERGY CALCULATIONS BASED ON THE HARRIS-FUNCTIONAL APPROXIMATION COUPLED WITH THE GENERALIZED-GRADIENT APPROXIMATION

We have developed a molecular-orbital (MO) calculational method, based on the Harris-functional approximation coupled with the generalized-gradient approximation (Harris-GGA), in order to get reasonable molecular-bond energies of large systems for which the Kohn-Sham self-consistent-field (SCF) calculation is impractical. It has been applied to some diatomic molecules and three types of cage-shaped carbon cluster. For the diatomic molecules Li2, Be2, C2, N2, F2, Cu2, and CO, the bond-energy differences between the Harris-GGA and experimental values are 30- 50 % less than those between the Harris-functional approximation with a local-density approximation (Harris-LDA) and experimental values. For the carbon clusters, C-8 and C24, the bond energies calculated with use of the Harris-GGA are in agreement with those calculated with use of the SCF-GGA within about 10%. For the C60 fullerene, the calculated bond lengths and highest-occupied-MO-lowest-unoccupied-MO gap energy are comparable to experimental values.