On the basis-set dependence of local and integrated electron density properties: Application of a new computer program for quantum-chemical density analysis

A new computer program for post-processing analysis of quantum-chemical electron densities is described. The code can work with Slater- and Gaussian-type basis functions of arbitrary angular momentum. It has been applied to explore the basis-set dependence of the electron density and its Laplacian in terms of local and integrated topological properties. Our analysis, including Gaussian/Slater basis sets up to sextuple/quadruple-zeta order, shows that these properties considerably depend on the choice of type and number of primitives utilized in the wavefunction expansion. Basis sets with high angular momentum (l = 5 or l = 6) are necessary to achieve convergence for local properties of the density and the Laplacian. In agreement with previous studies, atomic charges defined within Bader's Quantum Theory of Atoms in Molecules appear to be much more basis-set dependent than the Hirshfeld's stockholder charges. The former ones converge only at the quadruple-zeta/higher level with Gaussian/Slater functions.