Picture change error correction in the radial distributions of canonical orbital densities and total electron density of radon atom: the effect of the size of nucleus and the basis set limit

The 2nd order Douglas-Kroll-Hess (DKH2) and the Infinite Order Two Component (IOTC) radial distributions of electron density of canonical Hartree-Fock (HF) orbitals of radon atom are presented. Furthermore, the total electron density is revisited. The picture change error (PCE) correction is investigated by analytical means. The point charge model of nucleus and the Gaussian nucleus model are employed. The basis set is extrapolated by means of including tight s and also p Gaussians within the original triple zeta basis set. It is found that the DKH1 PCE corrected DKH2 total electron and s orbital contact densities are negative for the point charge model of nucleus if tight enough s Gaussians are included in the basis set. It is shown that this failure is caused due to the missing terms of the second order Douglas-Kroll transformation for the DKH2 electron density. PCE is found the most striking in the DKH2/IOTC electron density of s orbitals close to the nucleus. The radial distributions of the 2-component p _1/2 orbital densities are considerably affected by PCE at the nucleus as well. Furthermore, the PCE corrected DKH2/IOTC scalar p orbital densities have a non-zero value of electron density at nucleus and can be considered as an spin-orbit (SO) average of the p _1/2 and p _3/2 orbitals. The d and f orbitals are affected by PCE in the vicinity of the nucleus only little. The PCE corrected DKH2 and IOTC radial distributions of orbital densities are nodeless, which is completely in agreement with the radial distribution of the analytic or numeric DCH orbital densities.