Non-empirical pseudopotentials in the HF-LCAO approach to crystalline solids

Abstract

The core pseudopotentials of Durand and Barthelat have been introduced in the periodic Hartree-Fock program CRYSTAL. This new option allows the treatment of crystals involving heavy atoms and therefore extends the applicability of the method. A comparison between all-electron and pseudopotential results is presented. Binding energies, optimized cell parameters, bulk moduli, band structures and electron density curves calculated with the two methods are presented for covalent (diamond and silicon), metallic (lithium) and ionic (MgO) systems. As expected the pseudopotential results for covalent compounds are in good agreement with those obtained from all-electron calculations; for MgO larger discrepancies are observed, possibly due to core polarization effects.