The frequently-used Zk and Yk functions in the framework of numerical Hartree–Fock method are formulated for atomic systems embedded in weakly-coupled plasmas modeled by Debye–Hückel potential. The Gegenbauer expansion is used to decompose the screened interelectronic Coulomb potential. Two numerical methods, the integration and differentiation ones, are developed to accurately and efficiently calculate the Zk and Yk integrals with each one showing special advantages in the implementation of numerical Hartree–Fock programs. Comparing to the analytical results specialized for Slater-type radial functions, it is generally concluded that the integration method possesses more accuracy and robustness than the differentiation one. Calculations on the two-electron C4+ ion in three different screening models are performed to demonstrate the applicability and stability of present methods in a wider range of screening parameters than existing Hartree–Fock calculations. With benchmark Hartree–Fock energies obtained in this work and accurate non-relativistic energies calculated based on Hylleraas-CI basis functions, it is found that the correlation energies are smallest in the fully screening model.
Read full abstract