Orbital-free effective embedding potential at nuclear cusps

A strategy to construct approximants to the kinetic-energy-functional dependent component (v(t)[rho(A),rho(B)](r)) of the effective potential in one-electron equations for orbitals embedded in a frozen-density environment [Eqs. (20) and (21) in Wesolowski and Warshel, J. Phys. Chem. 97, (1993) 8050] is proposed. In order to improve the local behavior of the orbital-free effective embedding potential near nuclei in the environment, the exact behavior of v(t)[rho(A),rho(B)](r) at rho(A)-> 0 and integral rho(B)dr=2 is taken into account. As a result, the properties depending on the quality of this potential are invariably improved compared to the ones obtained using conventional approximants which violated the considered exact condition. The approximants obtained following the proposed strategy and especially the simplest one constructed in this work are nondecomposable, i.e., cannot be used to obtain the analytic expression for the functional of the total kinetic energy. (C) 2008 American Institute of Physics.