Atomic charges,dipole moments,and Fukui functions using the Hirshfeld partitioning of the electron density

In the Hirshfeld partitioning of the electron density, the molecular electron density is decomposed in atomic contributions, proportional to the weight of the isolated atom density in the promolecule density, constructed by superimposing the isolated atom electron densities placed on the positions the atoms have in the molecule. A maximal conservation of the information of the isolated atoms in the atoms-in-molecules is thereby secured. Atomic charges, atomic dipole moments, and Fukui functions resulting from the Hirshfeld partitioning of the electron density are computed for a large series of molecules. In a representative set of organic and hypervalent molecules, they are compared with other commonly used population analysis methods. The expected bond polarities are recovered, but the charges are much smaller compared to other methods. Condensed Fukui functions for a large number of molecules, undergoing an electrophilic or a nucleophilic attack, are computed and compared with the HOMO and LUMO densities, integrated over the Hirshfeld atoms in molecules.     

Atoms-in-molecules from the Stockholder Partition of the Molecular Two-electron Distribution

The effective one-electron distributions of bonded atoms obtained from the “stockholder” partition of the molecular two-electron density are reported. These two-electron stockholder (S) atoms are compared with their one-electron analogs represented by the corresponding Hirshfeld (H) one-electron stockholder pieces of the molecular electron density. The influence of the exchange (Fermi) and Coulomb correlation between electrons on the resultant shapes of bonded atoms is investigated The vertical (for the fixed molecular electron density) and horizontal (involving the electron density displacement) correlation influences on the two-electron stockholder atoms are examined. The two sets of bonded stockholder atoms in the near-dissociation bond-elongated diatomics are compared for different approximations of the electron correlation effects. The cluster components in atomic resolution of the S-partitioning scheme are investigated for illustrative homonuclear and heteronuclear diatomics: H₂, LiH, HF, LiF, and N₂. This framework facilitates an understanding of the origins of the observed differences between the S and H variants of Atoms-in-Molecules. With the exception of hydrogen atoms, especially in light molecules, the two sets of bonded atoms were found to be practically identical. For H₂ and LiH the S atoms were shown to exhibit a distinctly higher degree of the bonding character, compared to their H analogs. The main electron correlation effects have been found to be well represented already at the exchange-only level, e.g., in the unrestricted Hartree–Fock (UHF) theory. An inclusion of the extra vertical Coulomb correlation exerts a marginal moderating influence on the ionic/covalent composition of the chemical bond already predicted by the UHF approximation, in the direction of a slightly more covalent (less ionic) bond character. The horizontal shifts of the molecular density due to Coulomb correlation, relative to the UHF reference, often act in the opposite direction.     

Kinetic energy decomposition scheme based on information theory

We proposed a novel kinetic energy decomposition analysis based on information theory. Since the Hirshfeld partitioning for electron densities can be formulated in terms of Kullback–Leibler information deficiency in information theory, a similar partitioning for kinetic energy densities was newly proposed. The numerical assessments confirm that the current kinetic energy decomposition scheme provides reasonable chemical pictures for ionic and covalent molecules, and can also estimate atomic energies using a correction with viral ratios.     

Tracing the Fingerprint of Chemical Bonds within the Electron Densities of Hydrocarbons: A Comparative Analysis of the Optimized and the Promolecule Densities

The equivalence of the molecular graphs emerging from the comparative analysis of the optimized and the promolecule electron densities in two hundred and twenty five unsubstituted hydrocarbons was recently demonstrated [Keyvani et al. Chem. Eur. J. 2016 , 22, 5003]. Thus, the molecular graph of an optimized molecular electron density is not shaped by the formation of the C?H and C?C bonds. In the present study, to trace the fingerprint of the C?H and C?C bonds in the electron densities of the same set of hydrocarbons, the amount of electron density and its Laplacian at the (3, ?1) critical points associated with these bonds are derived from both optimized and promolecule densities, and compared in a newly proposed comparative analysis. The analysis not only conforms to the qualitative picture of the electron density build up between two atoms upon formation of a bond in between, but also quantifies the resulting accumulation of the electron density at the (3, ?1) critical points. The comparative analysis also reveals a unified mode of density accumulation in the case of 2318 studied C?H bonds, but various modes of density accumulation are observed in the case of 1509 studied C?C bonds and they are classified into four groups. The four emerging groups do not always conform to the traditional classification based on the bond orders. Furthermore, four C?C bonds described as exotic bonds in previous studies, for example the inverted C?C bond in 1,1,1‐propellane, are naturally distinguished from the analysis.     

A self-consistent Hirshfeld method for the atom in the molecule based on minimization of information loss

Based on the so-called Hirshfeld atom in the molecule scheme, a new AIM method is presented. The method is similar to the Hirshfeld-I scheme, with the AIM weight function being constructed by minimizing the information loss upon formation of the molecule, but now requiring explicitly that the promolecular densities integrate to the same number of electrons as the AIM densities. This new weight function leads to a new iterative AIM scheme, and the resulting operative scheme is examined and discussed. The final results indicate that the newly proposed method does not perform as well as the Hirshfeld-I method.     

Comment on “Extending Hirshfeld‐I to bulk and periodic materials”

In recent years, several methods have been developed that partition the electron density among atoms using spherically symmetric atomic weights. D. E. P. Vanpoucke, P. Bultinck, and I. Van Driessche (J. Comput. Chem. 2012, doi: 10.1002/jcc.23088) recently reported a periodic implementation of the Hirshfeld‐I method that uses a combination of Becke‐style and uniform integration grids and modified atomic reference densities to compute net atomic charges in periodic materials. Herein, this method is discussed in the context of earlier periodic implementations of the Hirshfeld‐I method, the Iterated Stockholder Atoms method, and the density derived electrostatic and chemical method.     

To What Extent are “Atoms in Molecules” Structures of Hydrocarbons Reproducible from the Promolecule Electron Densities?

The “atoms in molecules” structures of 225 unsubstituted hydrocarbons are derived from both the optimized and the promolecule electron densities. A comparative analysis demonstrates that the molecular graphs derived from these two types of electron densities at the same geometry are equivalent for almost 90 % of the hydrocarbons containing the same number and types of critical points. For the remaining 10 % of molecules, it is demonstrated that by inducing small perturbations, through the variation of the used basis set or slight changes in the used geometry, the emerging molecular graphs from both densities are also equivalent. Interestingly, the (3, ?1) critical point between two “non‐bonded” hydrogen atoms, which triggered “H?H bonding” controversy is also observed in the promolecule densities of certain hydrocarbons. Evidently, the topology of the electron density is not dictated by chemical bonds or strong interactions and deformations induced by the interactions of atoms in molecules have a quite marginal role, virtually null, in shaping the general traits of the topology of molecular electron densities of the studied hydrocarbons, whereas the key factor is the underlying atomic densities.     

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. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Critical analysis and extension of the Hirshfeld atoms in molecules

The computational approach to the Hirshfeld [Theor. Chim. Acta 44, 129 (1977)] atom in a molecule is critically investigated, and several difficulties are highlighted. It is shown that these difficulties are mitigated by an alternative, iterative version, of the Hirshfeld partitioning procedure. The iterative scheme ensures that the Hirshfeld definition represents a mathematically proper information entropy, allows the Hirshfeld approach to be used for charged molecules, eliminates arbitrariness in the choice of the promolecule, and increases the magnitudes of the charges. The resulting "Hirshfeld-I charges" correlate well with electrostatic potential derived atomic charges.     

Performance of 3D‐space‐based atoms‐in‐molecules methods for electronic delocalization aromaticity indices

Several definitions of an atom in a molecule (AIM) in three‐dimensional (3D) space, including both fuzzy and disjoint domains, are used to calculate electron sharing indices (ESI) and related electronic aromaticity measures, namely, I_ring and multicenter indices (MCI), for a wide set of cyclic planar aromatic and nonaromatic molecules of different ring size. The results obtained using the recent iterative Hirshfeld scheme are compared with those derived from the classical Hirshfeld method and from Bader's quantum theory of atoms in molecules. For bonded atoms, all methods yield ESI values in very good agreement, especially for C–C interactions. In the case of nonbonded interactions, there are relevant deviations, particularly between fuzzy and QTAIM schemes. These discrepancies directly translate into significant differences in the values and the trends of the aromaticity indices. In particular, the chemically expected trends are more consistently found when using disjoint domains. Careful examination of the underlying effects reveals the different reasons why the aromaticity indices investigated give the expected results for binary divisions of 3D space. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011.