Search

We argue that when one divides a molecular property into atom‐in‐a‐molecule contributions, one should perform the division based on the property density of the quantity being partitioned. This is opposition to the normal approach, where the electron density is given a privileged role in defining the properties of atoms‐in‐a‐molecule. Because partitioning each molecular property based on its own property density is inconvenient, we design a reference‐free approach that does not (directly) refer atomic property densities. Specifically, we propose a stockholder partitioning method based on relative influence of a molecule's atomic nuclei on the electrons at a given point in space. The resulting method does not depend on an “arbitrary” choice of reference atoms and it has some favorable properties, including the fact that all of the electron density at an atomic nucleus is assigned to that nucleus and the fact all the atoms in a molecule decay at a uniform asymptotic rate. Unfortunately, the resulting model is not easily applied to spatially degenerate ground states. Furthermore, the practical realizations of this strategy that we tried here gave disappointing numerical results. © 2017 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

Atoms in molecules methods that rely on reference promolecular densities typically require that one define, or otherwise determine, the densities of unbound atomic anions. Whereas the isolated atomic polyanions are always physically and computationally unbound, monoanions can be either physically bound but computationally unbound (like the oxygen anion at the Hartree-Fock level of theory), or physically unbound but computationally bound (like the nitrogen anion using many DFT methods with a basis set including diffuse functions). Depending on the level of theory and basis set used, the densities of negatively charged atomic ions can decay very slowly and even be nonmonotonically decreasing. These delocalized anionic densities induce ill-behaved atomic properties for compounds containing highly reduced atoms. To treat the problem of unphysical proatom densities in iterative Hirshfeld methods, we compute the smallest (typically fractional) nuclear charge to bind all electrons, called the effective nuclear charge $$ {Z}_{\mathrm{A}}^{\mathrm{eff}} $$ of an atom A. When $$ {Z}_{\mathrm{A}}^{\mathrm{eff}}>{Z}_{\mathrm{A}} $$ at a given level of theory, the scaled density corresponding to the effective nuclear charge is used as the negatively charged proatom density. This novel approach dramatically improves the computational robustness of the iterative Hirshfeld partitioning scheme. [ABSTRACT FROM AUTHOR]