P.W. Ayers

We propose a framework to construct the ground-state energy and density matrix of an N-electron system by solving a self-consistent set of single-particle equations. The method can be viewed as a nontrivial extension of the Kohn-Sham scheme (which is embedded as a special case). It is based on separating the Green’s function into a quasiparticle part and a background part, and expressing only the background part as a functional of the density matrix. The calculated single-particle energies and wave functions have a clear physical interpretation as quasiparticle energies and orbitals.

A previous study of diatomic molecules revealed that variational second-order density matrix theory has serious problems in the dissociation limit when the N-representability is imposed at the level of the usual two-index (P,Q,G) or even three-index (T1,T2) conditions [ H. Van Aggelen et al., Phys. Chem. Chem. Phys. 11, 5558 (2009) ]. Heteronuclear molecules tend to dissociate into fractionally charged atoms. In this paper we introduce a general class of N-representability conditions, called subsystem constraints, and show that they cure the dissociation problem at little additional computational cost. As a numerical example the singlet potential energy surface of Be B+ is studied. The extension to polyatomic molecules, where more subsystem choices can be identified, is also discussed.

A double-atom partitioning of the molecular one-electron density matrix is used to describe atoms and bonds. All calculations are performed in Hilbert space. The concept of atomic weight functions (familiar from Hirshfeld analysis of the electron density) is extended to atomic weight matrices. These are constructed to be orthogonal projection operators on atomic subspaces, which has significant advantages in the interpretation of the bond contributions. In close analogy to the iterative Hirshfeld procedure, self-consistency is built in at the level of atomic charges and occupancies. The method is applied to a test set of about 67 molecules, representing various types of chemical binding. A close correlation is observed between the atomic charges and the Hirshfeld-I atomic charges.

A previously introduced partitioning of the molecular one-electron density matrix over atoms and bonds [D. Vanfleteren et al., J. Chem. Phys. 133, 231103 (2010)] is investigated in detail. Orthogonal projection operators are used to define atomic subspaces, as in Natural Population Analysis. The orthogonal projection operators are constructed with a recursive scheme. These operators are chemically relevant and obey a stockholder principle, familiar from the Hirshfeld-I partitioning of the electron density. The stockholder principle is extended to density matrices, where the orthogonal projectors are considered to be atomic fractions of the summed contributions. All calculations are performed as matrix manipulations in one-electron Hilbert space. Mathematical proofs and numerical evidence concerning this recursive scheme are provided in the present paper. The advantages associated with the use of these stockholder projection operators are examined with respect to covalent bond orders, bond polarization, and transferability.

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.