ABSTRACT:Haberlen and Rosch HR demonstrated Chem. Phys. Lett. 199, 491Ž.x Ž . 1992 the feasibility of performing scalar-relativistic, density functional theory DFT Ž . linear combination of Gaussian-type orbitals᎐fitting function LCGTO-FF calculations on clusters of atoms using an ''incomplete'' Douglas᎐Kroll transformation. Some of the approximations used in their multiatom calculations have not yet been fully explored for isolated atoms, especially the neglect of matrix elements involving vector products of the momentum and the use of fitting functions. In this investigation, scalar-relativistic Ž . LCGTO calculations without fitting functions were carried out at four levels of approximation, including that employed by HR, for four atoms: Ce, Au, Pb, and Pu. The results of these calculations suggest that the HR approximation should be applicable to atoms through the light actinides, so long as spin᎐orbit effects are unimportant. A comparison of very large basis-set results with results obtained using ordinary-size basis Ž . sets indicates that the valence bonding states from scalar-relativistic LCGTO-FF calculations should be no more sensitive to the orbital basis-set size than are their nonrelativistic counterparts.