A comprehensive analysis of molecule-intrinsic quasi-atomic, bonding, and correlating orbitals. I. Hartree-Fock wave functions

Accurately extracting the signature of intermolecular interactions present in the NCI plot of the reduced density gradient versus electron density

An electron density (ED)-based methodology is developed for the automatic identification of intermolecular interactions using pro-molecular density. The expression of the ED gradient in terms of atomic components furnishes the basis for the Independent Gradient Model (IGM). This model leads to a density reference for non interacting atoms/fragments where the atomic densities are added whilst their interaction turns off. Founded on this ED reference function that features an exponential decay also in interference regions, IGM model provides a way to identify and quantify the net ED gradient attenuation due to interactions. Using an intra/inter uncoupling scheme, a descriptor (δg^inter) is then derived that uniquely defines intermolecular interaction regions. An attractive feature of the IGM methodology is to provide a workflow that automatically generates data composed solely of intermolecular interactions for drawing the corresponding 3D isosurface representations.

Intrinsic Resolution of Molecular Electronic Wave Functions and Energies in Terms of Quasi-atoms and Their Interactions

A general intrinsic energy resolution has been formulated for strongly correlated wave functions in the full molecular valence space and its subspaces. The information regarding the quasi-atomic organization of the molecular electronic structure is extracted from the molecular wave function without introducing any additional postulated model state wave functions. To this end, the molecular wave function is expressed in terms of quasi-atomic molecular orbitals, which maximize the overlap between subspaces of the molecular orbital space and the free-atom orbital spaces. As a result, the molecular wave function becomes the superposition of a wave function representing the juxtaposed nonbonded quasi-atoms and a wave function describing the interatomic electron migrations that create bonds through electron sharing. The juxtaposed nonbonded quasi-atoms are shown to consist of entangled quasi-atomic states from different atoms. The binding energy is resolved as a sum of contributions that are due to quasi-atom formation, quasiclassical electrostatic interactions, and interatomic interferences caused by electron sharing. The contributions are further resolved according to orbital interactions. The various transformations that generate the analysis are determined by criteria that are independent of the working orbital basis used for calculating the molecular wave function. The theoretical formulation of the resolution is quantitatively validated by an application to the C₂ molecule.

A Comprehensive Analysis in Terms of Molecule-Intrinsic, Quasi-Atomic Orbitals. III. The Covalent Bonding Structure of Urea

The analysis of molecular electron density matrices in terms of quasi-atomic orbitals, which was developed in previous investigations, is quantitatively exemplified by a detailed application to the urea molecule. The analysis is found to identify strong and weak covalent bonding interactions as well as intramolecular charge transfers. It yields a qualitative as well as quantitative ab initio description of the bonding structure of this molecule, which raises questions regarding some traditional rationalizations.

Covalent bonds are created by the drive of electron waves to lower their kinetic energy through expansion

An analysis based on the variation principle shows that in the molecules H2 (+), H2, B2, C2, N2, O2, F2, covalent bonding is driven by the attenuation of the kinetic energy that results from the delocalization of the electronic wave function. For molecular geometries around the equilibrium distance, two features of the wave function contribute to this delocalization: (i) Superposition of atomic orbitals extends the electronic wave function from one atom to two or more atoms; (ii) intra-atomic contraction of the atomic orbitals further increases the inter-atomic delocalization. The inter-atomic kinetic energy lowering that (perhaps counter-intuitively) is a consequence of the intra-atomic contractions drives these contractions (which per se would increase the energy). Since the contractions necessarily encompass both, the intra-atomic kinetic and potential energy changes (which add to a positive total), the fact that the intra-atomic potential energy change renders the total potential binding energy negative does not alter the fact that it is the kinetic delocalization energy that drives the bond formation.