Pavlyukh Y. and Berakdar J.

We present a procedure for the construction of accurate Kohn-Sham potentials of quasispherical molecules starting from the first-principles valence densities. The method is demonstrated for the case of icosahedral C₂₀²⁺ and C₆₀ molecules. Provided the density is N representable the Hohenberg-Kohn theorem guarantees the uniqueness of the obtained potentials. The potential is iteratively built following the suggestion of R. van Leeuwen and E. J. Baerends [Phys. Rev. A 49, 2421 (1994)]. The high symmetry of the molecules allows a parametrization of the angular dependence of the densities and the potentials using a small number of symmetry-adapted spherical harmonics. The radial behavior of these quantities is represented on a grid and the density is reconstructed from the approximate potential by numerically solving the coupled-channel Kohn-Sham equations. Subsequently, the potential is updated and the procedure is continued until convergence is achieved.

PDF (794kB)

(For personal use only. Please send an email to our secretary if you need the paper.)