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

Molecular Tailoring Approach for the Direct Estimation of Individual Noncovalent Interaction Energies in Molecular Systems

The noncovalent interactions (NCIs) are omnipresent in chemistry, physics, and biology. The study of such interactions offers insights into various physicochemical phenomena. Some indirect approaches proposed in the literature for exploring the NCIs are briefly reviewed in Section 1 of this Perspective. These include: (i) Shift in the stretching frequency of an X-Y bond involved in X-Y···Z interaction. (ii) Topological analysis of molecular electron density. (iii) Empirical equations derived employing experimental and theoretical quantities. However, a direct method for estimating individual intramolecular/intermolecular interaction energies has been conspicuous by its absence from the literature. We have developed a molecular tailoring approach (MTA)-based method enabling a direct and reliable estimation of the energy of intra- as well as intermolecular interactions. This method offers a direct and reliable estimation of these interactions, in particular of the hydrogen bonds (HB) in molecules/weakly bound clusters along with the respective cooperativity contribution. In Section 2, the basis of our method is discussed, along with some illustrative examples. The application of this method to a variety of molecules and clusters, with a special emphasis on estimating the HB energy along with the energy of other NCIs is presented in Section 3. Section 4 discusses some computational strategies for applying our method to large molecular clusters. The last Section provides a summary and a discussion on future developments.

Analysis of bonding motifs in unusual molecules I: planar hexacoordinated carbon atoms

The bonding structures of CO3Li3+ and CS3Li3+ are studied by means of oriented quasi-atomic orbitals (QUAOs) to assess the possibility of these molecules being planar hexacoordinated carbon (phC) systems. CH3Li and CO32- are employed as reference molecules. It is found that the introduction of Li+ ions into the molecular environment of carbonate has a greater effect on the orbital structure of the O atoms than it does on the C atom. Partial charges computed from QUAO populations imply repulsion between the positively charged C and Li atoms in CO3Li3+. Upon the transition from CO3Li3+ to CS3Li3+, the analysis reveals that the substitution of O atoms by S atoms inverts the polarity of the carbon-chalcogen σ bond. This is linked to the difference in s- and p-fractions of the QUAOs of C and S, as element electronegativities do not explain the observed polarity of the CSσ bond. Partial charges indicate that the larger electron population on the C atom in CS3Li3+ makes C-Li attraction possible. Upon comparison with the C-Li bond in methyllithium, it is found that the C-Li covalent interactions in CO3Li3+ and CS3Li3+ have about 14% and 6% of the strength of the C-Li covalent interaction in CH3Li, respectively. Consequently, it is concluded that only CS3Li3+ may be considered to be a phC system.

Quasi-atomic orbital analysis of halogen bonding interactions

A quasi-atomic orbital analysis of the halogen bonded NH3⋯XF complexes (X = F, Cl, Br, and I) is performed to gain insight into the electronic properties associated with these σ-hole interactions. It is shown that significant sharing of electrons between the nitrogen lone pair of the ammonia molecule and the XF molecule occurs, resulting in a weakening of the X-F bond. In addition, the N-X bond shows increasing covalent character as the size of the halogen atom X increases. While the Mulliken outer complex NH3⋯XF appears to be overall the main species, the strength of the covalent interaction of the N-X bond becomes increasingly similar to that of the N-X bond in the [NH3X]+ cation as the size of X increases.

Intermolecular interactions in clusters of ethylammonium nitrate and 1-amino-1,2,3-triazole

The intermolecular interaction energies, including hydrogen bonds (H-bonds), of clusters of the ionic liquid ethylammonium nitrate (EAN) and 1-amino-1,2,3-triazole (1-AT) based deep eutectic propellants (DeEP) are examined. 1-AT is introduced as a neutral hydrogen bond donor (HBD) to EAN in order to form a eutectic mixture. The effective fragment potential (EFP) is used to examine the bonding interactions in the DeEP clusters. The resolution of the Identity (RI) approximated second order Møller-Plesset perturbation theory (RI-MP2) and coupled cluster theory (RI-CCSD(T)) are used to validate the EFP results. The EFP method predicts that there are significant polarization and charge transfer effects in the EAN:1-AT complexes, along with Coulombic, dispersion and exchange repulsion interactions. The EFP interaction energies are in good agreement with the RI-MP2 and RI-CCSD(T) results. The quasi-atomic orbital (QUAO) bonding and kinetic bond order (KBO) analyses are additionally used to develop a conceptual and semi-quantitative understanding of the H-bonding interactions as a function of the size of the system. The QUAO and KBO analyses suggest that the H-bonds in the examined clusters follow the characteristic hydrogen bonding three-center four electron interactions. The strongest H-bonding interactions between the (EAN)1:(1-AT)n and (EAN)2:(1-AT)n (n = 1-5) complexes are observed internally within EAN; that is, between the ethylammonium cation [EA]+ and the nitrate anion ([NO3]-). The weakest H-bonding interactions occur between [NO3]- and 1-AT. Consequently, the average strengths of the H-bonds within a given (EAN)x:(1-AT)n complex decrease as more 1-AT molecules are introduced into the EAN monomer and EAN dimer. The QUAO bonding analysis suggests that 1-AT in (EAN)x:(1-AT)n can act as both a HBD and a hydrogen bond acceptor simultaneously. It is observed that two 1-AT molecules can form H-bonds to each other. Although the KBOs that correspond to H-bonding interactions in [EA]+:1-AT, [NO3]-:1-AT and between two 1-AT molecules are weaker than the H-bonds in EAN, those weak H-bond networks with 1-AT could be important to form a stable DeEP.

Analysis of the bonding in tetrahedrane and phosphorus-substituted tetrahedranes

The bonding structures of tetrahedrane, phosphatetrahedrane, diphosphatetrahedrane and triphosphatetrahedrane are studied by employing an intrinsic quasi-atomic orbital analysis. Ethane, cyclopropane and tetrahedral P4 are employed as reference systems. The orbital analysis is paired with the computation of strain energies via isodesmic reactions. The results show that the increase in geometric strain upon transition from ethane to cyclopropane to tetrahedrane weakens the CC bonds, despite leading to shorter C-C interatomic distances. With the increase in strain, the orbitals centered on C and involved in the bonding of the cage structure are observed to have elevated p-character, and the orbital structure of C deviates from sp3 hybridization. The systematic substitution of CH groups by P atoms in the cage structure of tetrahedrane leads to stronger CC bonds, larger angles in the cage structures of the resulting phosphatetrahedranes, lower p-character in the orbitals involved in the bonding of the cages, and lower strain energies. It is found that P is more amenable to strained molecular arrangements than is C, and that the propensity of a given atom to hybridize s and p functions, or the lack thereof, has implications in the stability of molecules with strained geometries. The combination of the calculations presented here with the existing literature provides insight into the apparent propensity of tetrahedrane and P4 to 'break' their tetrahedral structures. Trends in the bonding interactions, such as bond strengths, s- and p-orbital characters and charge transfer are identified and related to the strain energy observed in each of the analyzed systems.

Carbon-[n]Triangulenes and Sila-[n]Triangulenes: Which Are Planar?

Using our recently suggested concept of a quasi-molecule ("tile") and, in the case of the planarity here at stake, its generalization to larger than tetratomics, we explain why carbon [n]triangulenes tend to be planar, while hybrids, where just a few or even all a- or b-type carbon atoms are silicon-substituted (sila-[n]triangulenes), tend to be planar/nonplanar when compared with the unsubstituted carbon-[n]triangulenes. Because other spin states of the parent carbon- and sila-[n]triangulenes tend to correlate with the same tiles, it is conjectured that no structural changes are expected to depend on their spin state. Other polycyclic and sila-compounds are also discussed.

Energy natural orbital characterization of nonadiabatic electron wavepackets in the densely quasi-degenerate electronic state manifold

Dynamics and energetic structure of largely fluctuating nonadiabatic electron wavepackets are studied in terms of Energy Natural Orbitals (ENOs) [K. Takatsuka and Y. Arasaki, J. Chem. Phys. 154, 094103 (2021)]. Such huge fluctuating states are sampled from the highly excited states of clusters of 12 boron atoms (B₁₂), which have densely quasi-degenerate electronic excited-state manifold, each adiabatic state of which gets promptly mixed with other states through the frequent and enduring nonadiabatic interactions within the manifold. Yet, the wavepacket states are expected to be of very long lifetimes. This excited-state electronic wavepacket dynamics is extremely interesting but very hard to analyze since they are usually represented in large time-dependent configuration interaction wavefunctions and/or in some other complicated forms. We have found that ENO gives an invariant energy orbital picture to characterize not only the static highly correlated electronic wavefunctions but also those time-dependent electronic wavefunctions. Hence, we first demonstrate how the ENO representation works for some general cases, choosing proton transfer in water dimer and electron-deficient multicenter chemical bonding in diborane in the ground state. We then penetrate with ENO deep into the analysis of the essential nature of nonadiabatic electron wavepacket dynamics in the excited states and show the mechanism of the coexistence of huge electronic fluctuation and rather strong chemical bonds under very random electron flows within the molecule. To quantify the intra-molecular energy flow associated with the huge electronic-state fluctuation, we define and numerically demonstrate what we call the electronic energy flux.

Quantum Chemical Modeling of Propellant Degradation

An ab initio quantum chemical approach for the modeling of propellant degradation is presented. Using state-of-the-art bonding analysis techniques and composite methods, a series of potential degradation reactions are devised for a sample hydroxyl-terminated-polybutadiene (HTPB) type solid fuel. By applying thermochemical procedures and isodesmic reactions, accurate thermochemical quantities are obtained using a modified G3 composite method based on the resolution of the identity. The calculated heats of formation for the different structures produced presents an ∼2 kcal/mol average error when compared against experimental values.

N-Electron Valence Perturbation Theory with Reference Wave Functions from Quantum Computing: Application to the Relative Stability of Hydroxide Anion and Hydroxyl Radical

Quantum simulations of the hydroxide anion and hydroxyl radical are reported, employing variational quantum algorithms for near-term quantum devices. The energy of each species is calculated along the dissociation curve, to obtain information about the stability of the molecular species being investigated. It is shown that simulations restricted to valence spaces incorrectly predict the hydroxyl radical to be more stable than the hydroxide anion. Inclusion of dynamical electron correlation from nonvalence orbitals is demonstrated, through the integration of the variational quantum eigensolver and quantum subspace expansion methods in the workflow of N-electron valence perturbation theory, and shown to correctly predict the hydroxide anion to be more stable than the hydroxyl radical, provided that basis sets with diffuse orbitals are also employed. Finally, we calculate the electron affinity of the hydroxyl radical using an aug-cc-pVQZ basis on IBM's quantum devices.

Maximally valent orbitals in systems with non-ideal bond-angles: atomic Wannier orbitals guided by the Mayer bond order

In pursuit of a directed minimal set of basis for systems with non-ideal bond angles, in this work we find the exact orientation of the major overlapping orbitals along the nearest neighbouring coordination segments in a given system such that they maximally represent the covalent interactions throughout the system. We compute Mayer's bond order, akin to Wiberg's bond index, on the basis of atomic Wannier orbitals with customizable non-degenerate hybridization leading to variable orientations, constructed from first principles, in a representative variety of molecules and layered systems. We put such orbitals in perspective with unbiased maximally localized descriptions of bonding and non-bonding orbitals, and energetics to tunneling of electrons through them between nearest neighbours, to describe the different physical aspects of covalent interactions, which are not necessarily represented using a single unique set of atomic or bonding orbitals.

Resolving the π-assisted U-N σf-bond formation using quantum information theory

We model the potential energy profiles of the UO₂ (NCO)Cl₂^- → NUOCl₂^- + CO₂ reaction pathway [Y. Gong, V. Vallet, M. del Carmen Michelini, D. Rios and J. K. Gibson, J. Phys. Chem. A, 2014, 118, 325-330] using different pair coupled-cluster doubles (pCCD) methods. Specifically, we focus on pCCD and pCCD-tailored coupled cluster models in predicting relative energies for the various intermediates and transition states along the reaction coordinate. Furthermore, we augment our study on energetics with an orbital-pair correlation analysis of the complete reaction pathway that features two distinct paths. Our analysis of orbital correlations sheds new light on the formation and breaking of respective bonds between the uranium, oxygen, and nitrogen atoms along the reaction coordinates where the "yl" bond is broken and a nitrido compound formed. Specifically, the strengthening of the U-N σ_f-bond is assisted by a π-type interaction that is delocalized over the C-N-U backbone of the UO₂ (NCO)Cl₂^- complex.

Nature of chemical bond and potential barrier in an invariant energy-orbital picture

Physical nature of the chemical bond and potential barrier is studied in terms of energy natural orbitals (ENOs), which are extracted from highly correlated electronic wavefunctions. ENO provides an objective one-electron picture about energy distribution in a molecule, just as the natural orbitals (NOs) represent one electron view about electronic charge distribution. ENO is invariant in the same sense as NO is, that is, ENOs converge to the exact ones as a series of approximate wavefunctions approach the exact one, no matter how the methods of approximation are adopted. Energy distribution analysis based on ENO can give novel insights about the nature of chemical bonding and formation of potential barriers, besides information based on the charge distribution alone. With ENOs extracted from full configuration interaction wavefunctions in a finite yet large enough basis set, we analyze the nature of chemical bonding of three low-lying electronic states of a hydrogen molecule, all being in different classes of the so-called covalent bond. The mechanism of energy lowering in bond formation, which gives a binding energy, is important, yet not the only concern for this small molecule. Another key notion in chemical bonding is whether a potential basin is well generated stiff enough to support a vibrational state(s) on it. Based on the virial theorem in the adiabatic approximation and in terms of the energy and force analyses with ENOs, we study the roles of the electronic kinetic energy and its nuclear derivative(s) on how they determine the curvature (or the force constant) of the potential basins. The same idea is applied to the potential barrier and, thereby, the transition states. The rate constant within the transition-state theory is formally shown to be described in terms of the electronic kinetic energy and the nuclear derivatives only. Thus, the chemical bonding and rate process are interconnected behind the scenes. Besides this aspect, we pay attention to (1) the effects of electron correlation that manifests itself not only in the orbital energy but also in the population of ENOs and (2) the role of nonadiabaticity (diabatic state mixing), resulting in double basins and a barrier on a single potential curve in bond formation. These factors differentiate a covalent bond into subclasses.

Orbital contraction and covalent bonding

According to Ruedenberg's classic treatise on the theory of chemical bonding [K. Ruedenberg, Rev. Mod. Phys. 34, 326-376 (1962)], orbital contraction is an integral consequence of covalent bonding. While the concept is clear, its quantification by quantum chemical calculations is not straightforward, except for the simplest of molecules, such as H₂ ⁺ and H₂. This paper proposes a new, yet simple, approach to the problem, utilizing the modified atomic orbital (MAO) method of Ehrhardt and Ahlrichs [Theor. Chim. Acta 68, 231 (1985)]. Through the use of MAOs, which are an atom-centered minimal basis formed from the molecular and atomic density operators, the wave functions of the species of interest are re-expanded, allowing the computation of the kinetic energy (and any other expectation value) of free and bonded fragments. Thus, it is possible to quantify the intra- and interfragment changes in kinetic energy, i.e., the effects of contraction. Computations are reported for a number of diatomic molecules H₂, Li₂, B₂, C₂, N₂, O₂, F₂, CO, P₂, and Cl₂ and the polyatomics CH₃-CH₃, CH₃-SiH₃, CH₃-OH, and C₂H₅-C₂H₅ (where the single bonds between the heavy atoms are studied) as well as dimers of He, Ne, Ar, and the archetypal ionic molecule NaCl. In all cases, it is found that the formation of a covalent bond is accompanied by an increase in the intra-fragment kinetic energy, an indication of orbital contraction and/or deformation.

Anatomy of Classical Boron-Boron Bonding: Overlap and sp Dissonance

Classical bonding is predominantly understood using the insipid spⁿ hybridization for σ-bonds as well as π bonds and their delocalized variants. Because hybridization ignores intricate differences in the energy and size of valence atomic orbitals, its naïve application to classically bonded boron atoms leads to numerous surprises in bond strengths, frontier MOs/bands, and even geometry. Here we show that the sp dissonance caused by size mismatch between the valence s and p orbitals of boron plays a crucial role in its bonding, subtly distinct from that of carbon and silicon. Unlike the heavier p block elements, boron prefers to actively engage its compact 2s orbitals in bonding. This leads to the overreach of p-p σ-type overlap that reduces its magnitude in the entire B─B bonding range. Consequently, the π-type overlap remains substantial, making its electronic structure visibly distinct in saturated and unsaturated regimes. The deltahedral frameworks offer a compromise by breaking this symmetry-enforced dichotomy of classical σ- and π-type bonding and following alternate electron counts that suit the electron deficiency of the boron. The pathological anatomy of classical B─B σ-bonding also explains the origins of puzzling metallic character and disorder in their classical boride networks even with ideal electron count, unlike deltahedral borides. The implications of sp dissonance are illustrated in classical boron networks of various hybridizations, explaining the unusual preference for unique sp³ lattice with strained four-membered rings in CrB₄, origins of observed σ holes in MgB₂ that lead to its superconducting nature, and the absence of Peierls distortion in LiB.

Metavalent bonding in chalcogenides: DFT-chemical pressure approach

Understanding the chemical bond nature has attracted considerable attention as it is crucial to analyze and comprehend the different physical and chemical properties of materials. This work is considered a complementary part of our previous work in studying the nature of different types of bonding interactions in a wide variety of molecules and materials using the DFT Chemical Pressure (CP) approach. Recently, a new type of chemical bond, the metavalent bond (MVB), has been defined. We show how the CP formalism can be used to analyze and study the establishment of MVB in two chalcogenides, GeSe and PbSe, in a similar fashion as the electron localization function (ELF) profiles. This is accomplished by analyzing the CP maps of these two chalcogenides at different pressures (up to 40 GPa for GeSe and 10 GPa for PbSe). The CP maps show distinctive features related to the MVB, providing insights into the existence of such chemical interaction in the crystal structure of the two compounds. Similar to ELF profiles, CP maps can visualize and track the strength of the MVB in GeSe and PbSe under pressure.

From six to eight Π-electron bare rings of group-XIV elements and beyond: can planarity be deciphered from the "quasi-molecules" they embed?

Ab initio molecular orbital theory is used to study the structures of six and eight π-electron bare rings of group-XIV elements, and even larger [n]annulenes up to C₁₈H₁₈, including some of their mono-, di-, tri-, and tetra-anions. While some of the above rings are planar, others are nonplanar. A much spotlighted case is cyclo-octatetraene (C₈H₈), which is predicted to be nonplanar together with its heavier group-XIV analogues Si₈H₈ and Ge₈H₈, with the solely planar members of its family having the stoichiometric formulas C₄Si₄H₈ and C₄Ge₄H₈. A similar situation arises with the six π-electron bare rings, where benzene and substituted ones up to C₃Si₃H₆ or so are planar, while others are not. However, the explanations encountered in the literature find support in ab initio calculations for such species, often rationalized from distinct calculated features. Using second-order Møller-Plesset perturbation theory and, when affordable (particularly tetratomics, which may allow even higher levels), the coupled-cluster method including single, double, and perturbative triple excitations, a common rationale is suggested based on a novel concept of quasi-molecules or the (3+4)-atom partition scheme. Any criticism of tautology is therefore avoided. The same analysis has also been successfully applied to even larger [n]annulenes, to their mixed family members involving silicon and germanium atoms, and to the C₁₈ carbon ring. Furthermore, it has been extended to annulene anions to check the criteria of the popular Hückel rule for planarity and aromaticity. Exploratory work on cycloarenes is also reported. Besides a partial study of the involved potential energy surfaces, equilibrium geometries and harmonic vibrational frequencies have been calculated anew, for both the parent and the actual prototypes of the quasi-molecules.

Robust Pipek-Mezey Orbital Localization in Periodic Solids

We describe a robust method for determining Pipek-Mezey (PM) Wannier functions (WF), recently introduced by Jónsson et al. (J. Chem. Theor. Chem. 2017, 13, 460), which provide some formal advantages over the more common Boys (also known as maximally-localized) Wannier functions. The Broyden-Fletcher-Goldfarb-Shanno-based PMWF solver is demonstrated to yield dramatically faster convergence compared to the alternatives (steepest ascent and conjugate gradient) in a variety of one-, two-, and three-dimensional solids (including some with vanishing gaps) and can be used to obtain Wannier functions robustly in supercells with thousands of atoms. Evaluation of the PM functional and its gradient in periodic linear combination of atomic orbital representation used a particularly simple definition of atomic charges obtained by Moore-Penrose pseudoinverse projection onto the minimal atomic orbital basis. An automated "canonicalize phase then randomize" method for generating the initial guess for WFs contributes significantly to the robustness of the solver.

Bonding analysis of water clusters using quasi-atomic orbitals

The quasi-atomic orbital (QUAO) bonding analysis introduced by Ruedenberg and co-workers is used to develop an understanding of the hydrogen bonds in small water clusters, from the dimer through the hexamer (bag, boat, book, cyclic, prism and cage conformers). Using kinetic bond orders as a metric, it is demonstrated that as the number of waters in simple cyclic clusters increases, the hydrogen bonds strengthen, from the dimer through the cyclic hexamer. However, for the more complex hexamer isomers, the strength of the hydrogen bonds varies, depending on whether the cluster contains double acceptors and/or double donors. The QUAO analysis also reveals the three-center bonding nature of hydrogen bonds in water clusters.

Hybrid Atomic Orbital Basis from First Principles: Bottom-Up Mapping of Self-Energy Correction to Large Covalent Systems

Construction of hybrid atomic orbitals is proposed as the approximate common eigenstates of finite first moment matrices. Their hybridization and orientation can be a priori tuned as per their anticipated neighborhood. Their Wannier function counterparts constructed from the Kohn-Sham (KS) single particle states constitute an orthonormal multiorbital tight binding (TB) basis resembling hybrid atomic orbitals locked to their immediate atomic neighborhood, while spanning the subspace of KS states. The proposed basis thus renders predominantly single TB parameters from first principles for each nearest neighbor bond involving no more than two orbitals irrespective of their orientation and also facilitates an easy route for the transfer of such TB parameters across isostructural systems exclusively through mapping of neighborhoods and projection of orbital charge centers. With hybridized 2s, 2p and 3s, 3p valence electrons, the spatial extent of the self-energy correction (SEC) to TB parameters in the proposed basis is found to be localized mostly within the third nearest neighborhood, thus allowing effective transfer of self-energy-corrected TB parameters from smaller reference systems to much larger target systems, with nominal additional computational cost beyond that required for explicit computation of SEC in the reference systems. The proposed approach promises inexpensive estimation of the quasi-particle structures of large covalent systems with workable accuracy.

The Valence-Bond (VB) Model and Its Intimate Relationship to the Symmetric or Permutation Group

VB and molecular orbital (MO) models are normally distinguished by the fact the first looks at molecules as a collection of atoms held together by chemical bonds while the latter adopts the view that each molecule should be regarded as an independent entity built up of electrons and nuclei and characterized by its molecular structure. Nevertheless, there is a much more fundamental difference between these two models which is only revealed when the symmetries of the many-electron Hamiltonian are fully taken into account: while the VB and MO wave functions exhibit the point-group symmetry, whenever present in the many-electron Hamiltonian, only VB wave functions exhibit the permutation symmetry, which is always present in the many-electron Hamiltonian. Practically all the conflicts among the practitioners of the two models can be traced down to the lack of permutation symmetry in the MO wave functions. Moreover, when examined from the permutation group perspective, it becomes clear that the concepts introduced by Pauling to deal with molecules can be equally applied to the study of the atomic structure. In other words, as strange as it may sound, VB can be extended to the study of atoms and, therefore, is a much more general model than MO.

Multiple Bonding in Rhodium Monoboride. Quasi-atomic Analyses of the Ground and Low-Lying Excited States

The bonding structures of the ground state and the lowest five excited states of rhodium monoboride are identified by determining the quasi-atomic orbitals in full valence space MCSCF wave functions and the interactions between these orbitals. A quadruple bond, namely two π-bonds and two σ-bonds, is identified and characterized for the X¹Σ⁺ ground state, in agreement with a previous report (Cheung J. Phys. Chem. Lett. 2020, 11, 659-663). However, in all excited states, the bonding is predicted to be weaker because, in these states, one of the σ-bonding interactions has a small magnitude. In the a³Δ and A¹Δ states, the bond order is between a triple and quadruple bond. In the b³Σ⁺ state, the Rh-B linkage is a triple bond. In the c³Π and B¹Π states, the atoms are linked by a double bond due to an additional weakening of the two π-bonds. The decreases in the predicted bond strengths are reflected in the decreases of the predicted binding energies and in the increases of the predicted bond lengths from the X¹Σ⁺ ground state to the c³Π and the B¹Π excited states. Notably, the 5pσ orbital of rhodium, which is vacant in the ground state of the atom, plays a significant role in the molecule.

Substituent Effects on the Quantum Interference of Two-Center One-Electron Bonds: [B2X6]- (X = H, F, Cl, CN, OH, CH3, and OCH3)

The interference energy analysis (IEA) provided by the generalized product function energy partitioning (GPF-EP) method was applied to investigate the influence of the neighboring atoms on the nature of the two-center one-electron (2c1e) bonds in the anion dimers of BX₃ species (X = H, F, Cl, CN, OH, CH₃, and OCH₃). The species were studied at the GVB-PP(6/12).SC(1,2)/6-31**G++ level of calculation. The IEA has revealed that there is a balance between two main factors determining the chemical stability of the species. Quantum interference acts as the sole stabilizing effect in the formation of the chemical bonds, particularly as the result of the drop in kinetic energy, and the electronegativity of the substituent has a direct influence on the magnitude of this effect. The quasi-classical energy is responsible for the destabilizing factors, mainly the group bulkiness, and the "electron-withdrawing" effect in the case of the cyano group.

Generalization of Intrinsic Orbitals to Kramers-Paired Quaternion Spinors, Molecular Fragments, and Valence Virtual Spinors

Localization of molecular orbitals finds its importance in the representation of chemical bonding (and antibonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed on this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic acid, and ferrocene molecules, and then demonstrate the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium, and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.

Stability and Dissociation of Ethylenedione (OCCO)

This work examines the electronic structure and apparent instability of ethylenedione (OCCO), including an analysis of the singlet and triplet potential energy surfaces along the bending vibrations. While the singlet state is inherently unstable due to the Renner-Teller effect, theory predicts the triplet state to have a stable minimum on the potential energy surface. The stability of the triplet state is examined in detail, taking into account spin-orbit interactions. Using multireference quantum chemical methods, the lifetime of the triplet state is estimated to be in the picosecond range, significantly lower than previously computed. A quasi-atomic molecular orbital (QUAO) analysis is also used to elucidate the nature of bonding along the potential energy surface in both the singlet and triplet states. These calculations confirm the transient nature of the OCCO molecule, although they do not fully explain the lack of experimental detection via spectroscopy, which is known have the capability to probe even shorter lifetimes.

Why is Si2H2 Not Linear? An Intrinsic Quasi-Atomic Bonding Analysis

The molecular energy of Si₂H₂ geometric structures increases in the order dibridged < trans-bent < linear, in contrast to acetylene, C₂H₂, for which the linear structure is the global minimum. In this study, the intra-atomic (antibonding) and bonding contributions to the total molecular energy of these valence isoelectronic molecules are computed by expressing the density matrices of the full valence space multiconfiguration self-consistent field wave function in terms of quasi-atomic orbitals. The analysis shows that the intra-atomic contributions to the molecular energy become less favorable in the order dibridged → trans-bent → linear for both C₂H₂ and Si₂H₂. By contrast, the inter-atomic bonding contributions become energetically more favorable in that order for both C₂H₂ and Si₂H₂. The two systems differ as follows. For Si₂H₂, the antibonding intra-atomic energy changes that occur when the dibridged molecule reconstructs into the trans-bent and linear structures prevail over the interatomic interactions that induce bond formation. In contrast, for C₂H₂, the interatomic interactions that create bonds prevail over the intra-atomic energy changes that occur when the dibridged molecule reconstructs into the trans-bent and linear structures. The intra-atomic energy changes that occur in these systems are related to the hybridization of the heavy atoms in an analogous manner to the hybridization of C in CH₄ from (2s)²(2p)² to sp³ hybrid orbitals.

Recent developments in the general atomic and molecular electronic structure system

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.

On the connection between probability density analysis, QTAIM, and VB theory

Classification of bonds is essential for understanding and predicting the reactivity of chemical compounds. This classification mainly manifests in the bond order and the contribution of different Lewis resonance structures. Here, we outline a first principles approach to obtain these orders and contributions for arbitrary wave functions in a manner that is both, related to the quantum theory of atoms in molecules and consistent with valence bond theory insight: the Lewis structures arise naturally as attractors of the all-electron probability density |Ψ|2. Doing so, we introduce a valence bond weight definition that does not collapse in the basis set limit.

Quasi-Atomic Bond Analyses in the Sixth Period: II. Bond Analyses of Cerium Oxides

The role of the 4f orbitals in bonding is examined for the molecules cerium monoxide and cerium dioxide that have cerium formally in the +2 and +4 oxidation states, respectively. It is shown that the 4f orbitals are used primarily for polarization of the 5d orbitals when cerium is in the lower oxidation state, while the 4f orbitals play a significant role in chemical bonding via 5d/4f hybridization when cerium is in the +4 oxidation state.

Scalable Electron Correlation Methods. 6. Local Spin-Restricted Open-Shell Second-Order Møller-Plesset Perturbation Theory Using Pair Natural Orbitals: PNO-RMP2

We present a (near) linear scaling implementation of high-spin open-shell Møller-Plesset perturbation theory using pair natural orbitals (PNO-RMP2). The theory is based on a new variant of open-shell MP2 which is fully spin-adapted and uses a single set of spin-free amplitudes, as in closed-shell MP2. This method, denoted SROMP2, is invariant to unitary orbital transformations within the closed, open, and virtual orbital subspaces. Accordingly, only a single set of PNOs per spatial orbital pair is needed, and the efficiency is similar to closed-shell calculations. The PNOs are obtained using a semicanonical approximation with large domains of projected atomic orbitals (PAOs). Linear scaling is achieved provided that the open-shell orbitals are local, and distant pairs are treated by multipole approximations. The method is efficiently parallelized. The convergence of ionization and reaction energies as a function of the PAO and PNO domain sizes is demonstrated and found to be very similar as for closed-shell calculations. The suitability of the PNOs for explicitly correlated PNO-RCCSD-F12 calculations is also tested. So far, this method is only simulated using a conventional program with appropriate projections to the PAO and PNO subspaces. It is demonstrated for radical stabilization energies as well as ionization potentials that the errors caused by the local domain approximations with our default thresholds are negligible.

A Quasi-Atomic Analysis of Three-Center Two-Electron Zr-H-Si Interactions

A comprehensive analysis of the bonding structure of the disilyl zirconocene amide cation {Cp₂Zr[N(SiHMe₂)₂]}⁺ is conducted by application of an intrinsic orbital localization method that yields quasi-atomic orbitals (QUAOs). An emphasis is placed on describing a previously characterized three-center two-electron interaction between zirconium, hydrogen, and silicon that presents structural and spectroscopic features similar to that of agostic bonds. Expressions of the first-order density matrix in terms of the QUAOs yields bond orders (BOs), kinetic bond orders (KBOs), and the extent of transfer of charge that are useful to determine the electronic nature of the Zr-H-Si bond. The interactions between the QUAOs demonstrate the importance of vicinal interactions in the stabilization of the molecule. In addition, the evolution of the QUAOs during reactions with Lewis bases reveals the role of the Zr-H-Si interaction in facilitating the reaction.

Proton Transfer in 1,2,4-Triazolium Dinitramide: Effect of Aqueous Microsolvation

The gas phase proton transfer process in 1,2,4-triazolium dinitramide (TD) was studied using second-order perturbation theory to determine how the presence of one and two water molecules modifies the potential energy surface that connects the ion pair to the neutral pair. The presence of one water molecule can introduce small proton transfer energy barriers that separate the ion pair from the lower-energy neutral pair. These energy barriers are easily surmounted. Reaction paths were determined for single proton transfers and double proton transfers via one water molecule. In the presence of two water molecules, the global minimum is an ion pair, as are most of the lower-energy local minima. Energy barriers for single, double, and triple proton transfers were also found for TD in the presence of two water molecules. One TD ion pair structure with two water molecules has no corresponding neutral pair energy minimum. A quasi-atomic orbital analysis is used to understand the nature of the bonding in the various species studied in this work.

Chemical Pressure Maps of Molecules and Materials: Merging the Visual and Physical in Bonding Analysis

The characterization of bonding interactions in molecules and materials is one of the major applications of quantum mechanical calculations. Numerous schemes have been devised to identify and visualize chemical bonds, including the electron localization function, quantum theory of atoms in molecules, and natural bond orbital analysis, whereas the energetics of bond formation are generally analyzed in qualitative terms through various forms of energy partitioning schemes. In this Article, we illustrate how the chemical pressure (CP) approach recently developed for analyzing atomic size effects in solid state compounds provides a basis for merging these two approaches, in which bonds are revealed through the forces of attraction and repulsion acting between the atoms. Using a series of model systems that include simple molecules (H₂, CO₂, and S₈), extended structures (graphene and diamond), and systems exhibiting intermolecular interactions (ice and graphite), as well as simple representatives of metallic and ionic bonding (Na and NaH, respectively), we show how CP maps can differentiate a range of bonding phenomena. The approach also allows for the partitioning of the potential and kinetic contributions to the interatomic interactions, yielding schemes that capture the physical model for the chemical bond offered by Ruedenberg and co-workers.

Covalent Bonding in the Hydrogen Molecule

This work addresses the continuing disagreement between two schools of thought concerning the mechanism of covalent bonding. According to Hellmann, Ruedenberg, and Kutzelnigg, covalent bonding is a quantum mechanical phenomenon whereby lowering of the kinetic energy associated with electron sharing, i.e., delocalization, is the key stabilization mechanism. The opposing view of Slater, Feynman, and Bader has maintained that the source of stabilization is electrostatic potential energy lowering due to electron density redistribution to binding regions between nuclei. Following our study of H₂⁺ we present an analogous detailed study of H₂ where bonding involves an electron pair with repulsion and correlation playing a significant role in its properties. We use a range of different computational approaches to study and reveal the relevant contributions to bonding as seen in the electron density and corresponding kinetic and potential energy distributions. The energetics associated with the more complex electronic structure of H₂, when examined in detail, clearly agrees with the analysis of Ruedenberg; i.e., covalent bonding is due to a decrease in the interatomic kinetic energy resulting from electronic delocalization. Our results support the view that covalent bonding is a quantum dynamical phenomenon requiring a properly quantized kinetic energy to be used in its description.

Identification and Characterization of Molecular Bonding Structures by ab initio Quasi-Atomic Orbital Analyses

The quasi-atomic analysis of ab initio electronic wave functions in full valence spaces, which was developed in preceding papers, yields oriented quasi-atomic orbitals in terms of which the ab initio molecular wave function and energy can be expressed. These oriented quasi-atomic orbitals are the rigorous ab initio counterparts to the conceptual bond forming atomic hybrid orbitals of qualitative chemical reasoning. In the present work, the quasi-atomic orbitals are identified as bonding orbitals, lone pair orbitals, radical orbitals, vacant orbitals and orbitals with intermediate character. A program determines the bonding characteristics of all quasi-atomic orbitals in a molecule on the basis of their occupations, bond orders, kinetic bond orders, hybridizations and local symmetries. These data are collected in a record and provide the information for a comprehensive understanding of the synergism that generates the bonding structure that holds the molecule together. Applications to a series of molecules exhibit the complete bonding structures that are embedded in their ab initio wave functions. For the strong bonds in a molecule, the quasi-atomic orbitals provide quantitative ab initio amplifications of the Lewis dot symbols. Beyond characterizing strong bonds, the quasi-atomic analysis also yields an understanding of the weak interactions, such as vicinal, hyperconjugative and radical stabilizations, which can make substantial contributions to the molecular bonding structure.

Are One-Electron Bonds Any Different from Standard Two-Electron Covalent Bonds?

The nature of the chemical bond is perhaps the central subject in theoretical chemistry. Our understanding of the behavior of molecules developed amazingly in the last century, mostly with the rise of quantum mechanics (QM) and QM-based theories such as valence bond theory and molecular orbital theory. Such theories are very successful in describing molecular properties, but they are not able to explain the origin of the chemical bond. This problem was first analyzed in the 1960s by Ruedenberg, who showed that covalent bonds are the direct result of quantum interference between one-electron states. The generality of this result and its quantification were made possible through the recent development of the generalized product function energy partitioning (GPF-EP) method by our group, which allows the partitioning of the electronic density and energy into their interference and quasi-classical (noninterference) contributions. Furthermore, with GPF wave functions these effects can be analyzed separately for each bond of a molecule. This interference energy analysis has been applied to a large variety of molecules, including diatomics and polyatomics, molecules with single, double, and triple bonds, molecules with different degrees of polarity, linear or branched molecules, cyclic or acyclic molecules, conjugated molecules, and aromatics, in order to verify the role played by quantum interference. In all cases the conclusion is exactly the same: for each bond in each of the molecules considered, the main contribution to its stability comes from the interference term. Two-center one-electron (2c1e) bonds are the simplest kind of chemical bonds. Yet they are often viewed as odd or unconventional cases of bonding. Are they any different from conventional (2c2e) bonds? If so, what differences can we expect in the nature of (2c1e) bonds relative to electron-pair bonds? In this Account, we extend the GPF-EP method to describe bonds involving N electrons in M orbitals (N < M) and show its application to (2c1e) bonds. As examples we chose the molecules H₂⁺, H₃C·CH₃⁺, B₂H₄^-, [Cu·BH₃(PH₃)₃], and an alkali-metal cation dimer, and we evaluated the components of the electronic energy and density, which account for the formation of the bond, and compared the results with those for the respective analogous molecules exhibiting the "conventional" two-electron bond. In all cases, it was verified that interference is the dominant effect for the one-electron bonds. The GPF-EP results clearly indicate that molecules exhibiting (2c1e) bonds should not be considered as special systems, since one- and two-electron bonds result from quantum interference and therefore there is no conceptual difference between them. Moreover, these results show that quantum interference provides a way to unify the chemical bond concept.

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.