Binding energies of first row diatomics in the light of the interacting quantum atoms approach

Binding energies of first row diatomics are revisited within the interacting quantum atoms (IQA) approach. This is a formalism in chemical bonding theory based upon the quantum theory of atoms in molecules. It is characterized by the preservation of the energetic identity of atoms within molecules. Quantum mechanically computed binding energies are recovered in IQA as a sum of small atomic deformation energies and large pairwise interaction terms. We show how this partition responds faithfully to chemical intuition, and how the different evolution of deformations and interactions accounts in a unified manner for the subtle variations of the binding energy of these molecules.     

A Comparison of the Interacting Quantum Atoms (IQA) Analysis of the Two-Particle Density-Matrices of MP4SDQ and CCSD

Within the quantum topological energy partitioning method called Interacting Quantum Atoms (IQA) we transition from Møller-Plesset (MP4SDQ) to CCSD in calculating intra- and interatomic electron correlation energies for a set of hydrides, diatomics, a few simple molecules and non-covalently bonded complexes, using the uncontracted basis set 6-31++G(2d,2p). CCSD-IQA allows a more rigorous analysis of atomic electron correlation than that offered by Møller-Plesset, which returns IQA contributions that are identical to Hartree–Fock counterparts except for two-electron terms. The CCSD-IQA analysis returns bond and other interatomic correlation energies that are typically much larger in magnitude than the MP4SDQ values. Crisp patterns of energy transferability are detected in water clusters, both for intra-atomic and interatomic correlation energies. CCSD determines that the intra-atomic correlation energy of an oxygen drops by 15 kJ · mol^–1 for donating a hydrogen and by 25 kJ · mol^–1 for accepting a hydrogen.     

Hydrogen‐Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach

The cooperative effects of hydrogen bonding in small water clusters (H₂O)_n (n=3–6) have been studied by using the partition of the electronic energy in accordance with the interacting quantum atoms (IQA) approach. The IQA energy splitting is complemented by a topological analysis of the electron density (ρ( r )) compliant with the quantum theory of atoms‐in‐molecules (QTAIM) and the calculation of electrostatic interactions by using one‐ and two‐electron integrals, thereby avoiding convergence issues inherent to a multipolar expansion. The results show that the cooperative effects of hydrogen bonding in small water clusters arise from a compromise between: 1) the deformation energy (i.e., the energy necessary to modify the electron density and the configuration of the nuclei of the isolated water molecules to those within the water clusters), and 2) the interaction energy (E_int) of these contorted molecules in (H₂O)_n. Whereas the magnitude of both deformation and interaction energies is enhanced as water molecules are added to the system, the augmentation of the latter becomes dominant when the size of the cluster is increased. In addition, the electrostatic, classic, and exchange components of E_int for a pair of water molecules in the cluster (H₂O)_n?1 become more attractive when a new H₂O unit is incorporated to generate the system (H₂O)_n with the last‐mentioned contribution being consistently the most important part of E_int throughout the hydrogen bonds under consideration. This is opposed to the traditional view, which regards hydrogen bonding in water as an electrostatically driven interaction. Overall, the trends of the delocalization indices, δ(Ω,Ω′), the QTAIM atomic charges, the topology of ρ( r ), and the IQA results altogether show how polarization, charge transfer, electrostatics, and covalency contribute to the cooperative effects of hydrogen bonding in small water clusters. It is our hope that the analysis presented in this paper could offer insight into the different intra‐ and intermolecular interactions present in hydrogen‐bonded systems.     

Steric repulsions, rotation barriers, and stereoelectronic effects: a real space perspective

Widely used chemical concepts like Pauli repulsion or hyperconjugation, and their role in determining rotation barriers or stereoelectronic effects, are analyzed from the real space perspective of the interacting quantum atoms approach (IQA). IQA emerges from the quantum theory of atoms in molecules (QTAIM), but is free from the equilibrium geometry constraint of the former. A framework with both electronically unrelaxed and relaxed wavefunctions is presented that leads to an approximate correspondence between the IQA concepts and those used in the EDA (energy decomposition analysis) or NBO (natural bond orbital) procedures. We show that no net force acts upon the electrons in an electronically relaxed system, so that any reasonable definition of Pauli repulsion must involve unrelaxed state functions. Using antisymmetrized fragments clarifies that Pauli repulsions are energetically connected to the IQA deformation energies, leaving footprints in the finally relaxed states. Similarly, EDA or NBO hyperconjugative stabilizations are found to be naturally related to the IQA electron delocalization patterns. Applications to the rotation barrier of ethane and other simple systems are presented, and the very often forgotten role of electrostatic contributions in determining preferred conformations is highlighted.     

Efficient implementation of the interacting quantum atoms energy partition of the second-order Møller–Plesset energy

We describe an efficient implementation of the partition of the second-order Møller–Plesset (MP2) correlation energy within the interacting quantum atoms (IQA) energy decomposition. We simplify the IQA integration bottleneck by considering only the occupied to virtual elements of the second order reduced density matrix, a procedure that reduces substantially the size of the two-electron matrix, which has to be addressed. The algorithmic improvements described herein allow to perform the decomposition of the MP2 correlation energy for medium size molecular systems using moderate computational resources. We expect that the methods developed in this investigation will prove useful to understand electron correlation effects through a real space perspective.     

Reinvestigation of intramolecular hydrogen bond in malonaldehyde derivatives: An ab initio,AIM and NBO study

The RAHB systems in malonaldehyde and its derivatives at MP2/ 6‐311++G(d,p) level of theory were studied and their intramolecular hydrogen bond energies by using the related rotamers method was obtained. The topological properties of electron density distribution in O? H···O intramolecular hydrogen bond have been analyzed in term of quantum theory of atoms in molecules (QTAIM). Correlations between the H‐bond strength and topological parameters are probed. The results of QTAIM clearly showed that the linear correlation between the electron density distribution at HB critical point and RAHB ring critical point with the corresponding hydrogen bond energies was obtained. Moreover, it was found a linear correlation between the electronic potential energy density, V(r_cp), and hydrogen bond energy which can be used as a simple equation for evaluation of HB energy in complex RAHB systems. Finally, the similar linear treatment between the geometrical parameters, such as O···O or O? H distance, and Lp(O)→σ*_OH charge transfer energy with the intramolecular hydrogen bond energy is observed. © 2010 Wiley Periodicals, Inc., Int J Quantum Chem, 2011     

Hydrogen–hydrogen interaction in planar biphenyl: A theoretical study based on the interacting quantum atoms and Hirshfeld atomic energy partitioning methods

The nature of H‐H interaction between ortho‐hydrogen atoms in planar biphenyl is investigated by two different atomic energy partitioning methods, namely fractional occupation iterative Hirshfeld (FOHI) and interacting quantum atoms (IQA), and compared with the traditional virial‐based approach of quantum theory of atoms in molecules (QTAIM). In agreement with Bader's hypothesis of H? H bonding, partitioning the atomic energy into intra‐atomic and interatomic terms reveals that there is a net attractive interaction between the ortho‐hydrogens in the planar biphenyl. This falsifies the classical view of steric repulsion between the hydrogens. In addition, in contrast to the traditional QTAIM energy analysis, both FOHI and IQA show that the total atomic energy of the ortho‐hydrogens remains almost constant when they participate in the H‐H interaction. Although, the interatomic part of atomic energy of the hydrogens plays a stabilizing role during the formation of the H? H bond, it is almost compensated by the destabilizing effects of the intra‐atomic parts and consequently, the total energy of the hydrogens remains constant. The trends in the changes of intra‐atomic and interatomic energy terms of ortho‐hydrogens during H? H bond formation are very similar to those observed for the H₂ molecule. © 2014 Wiley Periodicals, Inc.     

Do one-step mechanisms always involve simultaneous evolution of electron density? QTAIM/IQA analysis of the Curtius rearrangement

The Curtius rearrangement reaction is studied by using quantum theory of atoms in molecules (QTAIM) analysis of the electron density and the interacting quantum atoms (IQA) formalism. Although the rearrangements take place in one stage, two phases are distinguished when the rearranged atom is H: the first one corresponds to the separation of N₂, and the second one to the N-H/C-H bond rearrangement. The transition state (TS) for the reaction does not represent an intermediate between reagent and product for the migration but for the isolation of the N₂ molecule. When the migration is undergone by a fluorine atom, no electronic phases can be distinguished and the process is really concerted. As the migration happens closer to the TS, the TS is more similar to the product. The IQA analysis reveals different electron density evolutions for H and F migrations, and the scarce relevance (in terms of energy) of the point where BCPs appear or disappear.     

The effects of higher orders of perturbation theory on the correlation energy of atoms and bonds in molecules

We examine, for the first time, the effects of higher orders of Møller–Plesset perturbation theory on the individual atoms within a molecule and the bonds between them, via the topological energy partitioning method of interacting quantum atoms. In real terms (i.e., not by absolute value) MP3 decreases the correlation energy of a bond, and MP4SDQ also decreases the energy of the atoms at either end of the bond. In addition, we investigated long‐range through‐space dispersive effects on a H₂ oligomer. Overall, MP3 is the largest correction to the correlation energy, and most of that energy is allocated to chemical bonds, reducing their values in actual terms. The MP4SDQ bond correlation correction, despite being relatively small, tends to have two effects: (i) for small or negative correlation energies MP4SDQ tends to decrease the bond correlation values even more, and (ii) for large (positive) bond correlation energies MP4SDQ tends to restore the bond correlation energies from the MP3 back toward the MP2 values. Furthermore, each individual part of a molecule or complex (atom or bond) has a specific convergence pattern for the MPn series: through‐space interactions converge at MP2 but bonds converge at MP3 level. The atomic correlation energy appears to head toward convergence at the MP4 level.     

Forced Bonding and QTAIM Deficiencies: A Case Study of the Nature of Interactions in He@Adamantane and the Origin of the High Metastability

Calculations within the framework of the interacting quantum atoms (IQA) approach have shown that the interactions of the helium atom with both tertiary, tC, and secondary, sC, carbon atoms in the metastable He@adamantane (He@adam) endohedral complex are bonding in nature, whereas the earlier study performed within the framework of Bader’s quantum theory of atoms in molecules (QTAIM) revealed that only He???tC interactions are bonding. The He???tC and He???sC bonding interactions are shown to be forced by the high pressure that the helium and carbon atoms exert upon each other in He@adam. The occurrence of a bonding interaction between the helium and sC atoms, which are not linked by a bond path, clearly shows that the lack of a bond path between two atoms does not necessarily indicate the lack of a bonding interaction, as is asserted by QTAIM. IQA calculations showed that not only the destabilization of the adamantane cage, but also a huge internal destabilization of the helium atom, contribute to the metastability of He@adam, these contributions being roughly equal. This result disproves previous opinions based on QTAIM analysis that only the destabilization of the adamantane cage accounts for the endothermicity of He@adam. Also, it was found that there is no homeomorphism of the ρ( r ) and ‐v( r ) fields of He@adam. Comparison of the IQA and QTAIM results on the interactions in He@adam exposes other deficiencies of the QTAIM approach. The reasons for the deficiencies in the QTAIM approach are analyzed.     

The Activation Strain Model in the Light of Real Space Energy Partitions

The distortion/interaction or activation strain model (ASM) of chemical reactivity is examined in real space through the interacting quantum atoms (IQA) approach. Attention is paid to the role that the geometrically constrained ASM structures of the fragments play in the chemical interpretation of the driving forces that lead to a given reaction channel. These fictitious intermediate states are necessary in the ASM, but IQA may or may not use them at will. Similarities and differences are highlighted by studying the endo/exo preference rules of simple [4+2] Diels–Alder cycloadditions. Although overall the agreement is reasonable, we warn about a blind use of the plain ASM if no further energy decomposition analyses of its interaction energy are done.     

Actinium coordination chemistry: A density functional theory study with monodentate and bidentate ligands

In the current study, the coordination chemistry of nine-coordinate Ac(III) complexes with 35 monodentate and bidentate ligands was investigated using density functional theory (DFT) in terms of their geometries, charges, reaction energies, and bonding interactions. The energy decomposition analysis with naturals orbitals for chemical valence (EDA-NOCV) and the quantum theory of atoms in molecules (QTAIM) were employed as analysis methods. Trivalent Ac exhibits the highest affinities toward hard acids (such as charged oxophilic donors, fluoride), so its classification as a hard acid is justified. Natural population analysis quantified the involvement of 5f orbitals on Ac to be about 30% of total valence electron natural configuration indicating that Ac is a member of the actinide series. Pearson correlation coefficients were used to study the pairwise correlations among the bond lengths, ΔG reaction energies, charges on Ac and donor atoms, and data from EDA-NOCV and QTAIM. Strong correlations and anticorrelations were found between Voronoi charges on donor atoms with ΔG, EDA-NOCV interaction energies and QTAIM bond critical point densities.     

Congested molecules. Where is the steric repulsion? An analysis of the electron density by the method of interacting quantum atoms

The computed electron density of several congested saturated hydrocarbons and halogenated derivatives has been analyzed by the method of interacting quantum atoms (IQA). For all the molecules studied, the calculations show the existence of a bond path between the congested atoms and which, according to the Quantum Theory of Atoms in Molecules, indicates that there is a stabilizing interaction between these atoms. The bond path is found to exist up to interatomic distances well‐beyond the sum of the van der Waals radii. The IQA results indicate that steric hindrance is not a repulsive force between the congested atoms but that is the result of an increase in the intra‐atomic or self‐energy of the congested atoms. This increase in self‐energy is caused by the deformation of the atomic basin of the congested atoms. © 2013 Wiley Periodicals, Inc.     

Performance of the density matrix functional theory in the quantum theory of atoms in molecules

The generalization to arbitrary molecular geometries of the energetic partitioning provided by the atomic virial theorem of the quantum theory of atoms in molecules (QTAIM) leads to an exact and chemically intuitive energy partitioning scheme, the interacting quantum atoms (IQA) approach, that depends on the availability of second-order reduced density matrices (2-RDMs). This work explores the performance of this approach in particular and of the QTAIM in general with approximate 2-RDMs obtained from the density matrix functional theory (DMFT), which rests on the natural expansion (natural orbitals and their corresponding occupation numbers) of the first-order reduced density matrix (1-RDM). A number of these functionals have been implemented in the promolden code and used to perform QTAIM and IQA analyses on several representative molecules and model chemical reactions. Total energies, covalent intra- and interbasin exchange-correlation interactions, as well as localization and delocalization indices have been determined with these functionals from 1-RDMs obtained at different levels of theory. Results are compared to the values computed from the exact 2-RDMs, whenever possible.     

On the origin of the relative stability of ZnIINTA and ZnIINTPA metal complexes. An insight from the IQA,IQF, and π‐FARMS methods

Relative stability of Zn^II complexes with nitrilotriacetic acid (ZnNTA) and nitrilotri‐3‐propionic acid (ZnNTPA) was investigated. Classical analysis of individual interactions using local indices failed to explain the preferential formation of ZnNTA. This work shows that the preferential formation of ZnNTA is not due to the size of coordination five‐membered rings or the absence of the steric CH‐‐HC contacts, as commonly considered. By combining interacting quantum atoms/fragments, IQA/IQF‐defined properties implemented in the π‐FARMS (preorganized‐interacting fragment attributed relative molecular stability) method, (i) several measures of Zn^II “affinity” to NTPA were shown to be consistently greater than to NTA and (ii) larger stability of ZnNTA was attributed to coordinated water molecules. Being smaller, NTA occupies less space around the metal center. This results in less destabilised Zn? OH₂ coordination bonds and preorganization energy of H₂O fragments being smaller in ZnNTA. Only by summing preorganization energies (of ligand and water fragments) and binding energy between fragments (using π‐FARMS method) we recovered the experimental trend. Importantly, the fundamental origin of all major energy components controlling relative stability of metal complexes was pin‐pointed using the π‐FARMS method.