The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures.

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.

Atoms and bonds in molecules as synergisms of interactions between electrons and nuclei.

The molecular electronic wave functions are expressed in quasi-atomic form. The corresponding global energy expression exhibits the unified resolution of the various modes of chemical bonding in terms of physical interactions. The bonding patterns in several molecules is elucidated.

Atoms and interatomic bonding synergism inherent in molecular electronic wave functions.

The chemical model of matter consists of atoms held together by bonds. The success of this model implies that the physical interactions of the electrons and nuclei in molecules combine into compound interactions that create the bonding. In the quantum mechanical description, the modified atoms in molecules and the bonding synergism are contained in the molecular electronic wave function. So far, only part of this information has been recovered from the wave function. Notably, the atoms have remained unidentified in the wave function. One reason is that conventional energy decomposition analyses formulate separate model wave functions, independent of the actual wave function, to represent "prepared atoms" and preconceived interactions and, then, intuitively catenate the parts. In the present work, the embedded modified atoms and the inherent physical synergisms between them are determined by a unified derivation entirely from the actual molecular valence space wave function. By means of a series of intrinsic orbital and configurational transformations of the wave function, the energy of formation of a molecule is additively resolved in terms of intra-atomic energy changes, interference energies, and quasi-classical, non-classical, and charge-transfer Coulombic interactions. The analysis furnishes an algorithm for the quantitative resolution of the energy of formation, which enables analyses elucidating reaction energies.

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 X1Σ+ 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 a3Δ and A1Δ states, the bond order is between a triple and quadruple bond. In the b3Σ+ state, the Rh-B linkage is a triple bond. In the c3Π and B1Π 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 X1Σ+ ground state to the c3Π and the B1Π 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.

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

The molecular energy of Si2H2 geometric structures increases in the order dibridged < trans-bent < linear, in contrast to acetylene, C2H2, 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 C2H2 and Si2H2. By contrast, the inter-atomic bonding contributions become energetically more favorable in that order for both C2H2 and Si2H2. The two systems differ as follows. For Si2H2, 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 C2H2, 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 CH4 from (2s)2(2p)2 to sp3 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.

Quasi-Atomic Bond Analyses in the Sixth Period: I. Relativistic Accurate Atomic Minimal Basis Sets for the Elements Cesium to Radon.

Full-valence relativistic accurate atomic minimal basis set (AAMBS) orbitals are developed for the sixth-row elements from cesium to radon, including the lanthanides. Saturated primitive atomic basis sets are developed and subsequently used to form the AAMBS orbitals. By virtue of the use of a saturated basis, properties computed based on the AAMBS orbitals are basis set independent. In molecules, the AAMBS orbitals can be used to construct valence virtual orbitals (VVOs) that provide chemically meaningful ab initio lowest unoccupied molecular orbitals (LUMOs) with basis set independent orbital energies. The optimized occupied molecular orbitals complemented with the VVOs form a set of full-valence molecular orbitals. They can be transformed into a set of oriented quasi-atomic orbitals (QUAOs) that provide information on intramolecular bonding via an intrinsic density analysis. In the present work, the development of the AAMBS for the sixth row is presented.

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.

The Virial Theorem and Covalent Bonding.

A long-held view of the origin of covalent binding is based on the notion that electrostatic forces determine the stability of a system of charged particles and that, therefore, potential energy changes drive the stabilization of molecules. A key argument advanced for this conjecture is the rigorous validity of the virial theorem. Rigorous in-depth analyses have however shown that the energy lowering of covalent bonding is due to the wave mechanical drive of electrons to lower their kinetic energy through expansion. Since the virial theorem applies only to systems with Coulombic interaction potentials, its relevance as a foundation of the electrostatic view is tested here by calculations on analogues of the molecules H2+ and H2, where all 1/ r interaction potentials are replaced by Gaussian-type potentials that yield one-electron "atoms" with realistic stability ranges. The virial theorem does not hold in these systems, but covalent bonds are found to form nonetheless, and the wave mechanical bonding analysis yields analogous results as in the case of the Coulombic potentials. Notably, the key driving feature is again the electron delocalization that lowers the interatomic kinetic energy component. A detailed discussion of the role of the virial theorem in the context of covalent binding is given.

Quasi-Atomic Bonding Analysis of Xe-Containing Compounds.

The origin of bonding in the rare-gas-containing molecules HXeCCH, HXeCCXeH, and HXeOXeH is explored using a quasi-atomic orbital (QUAO) analysis. The QUAOs provide qualitative and quantitative data about bonding through transformations of the density matrix. Bond orders, kinetic bond orders, and the extent of transfer of charge are analyzed. Localized molecular orbitals formed from the QUAOs provide additional insights about the relative polarity of bonds formed by xenon. The analysis suggests significant covalent bonding for Xe-Y (Y = C, O) as well as Xe-H, with both bonds using the same pσ-type orbital on Xe. These covalent interactions are established by substantial charge shifts from Xe to Y as well as to H. Accordingly, a covalent three-center four-electron bond links the atoms H-Xe-Y. On the basis of the analysis, electrostatic interactions do not play a significant role in the Xe-Y or Xe-H bonding.

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.

Dispersion Interactions in QM/EFP.

The dispersion energy term between quantum-mechanical (QM) and classical (represented by effective fragment potentials, EFP) subsystems is developed and implemented. A new formulation is based on long-range perturbation theory and uses dynamic polarizability tensors of the effective fragments and electric field integrals and orbital energies of the quantum-mechanical subsystem. No parametrization is involved. The accuracy of the QM-EFP dispersion energy is tested on a number of model systems; the average mean unsigned error is 0.8 kcal/mol or 13% with respect to the symmetry adapted perturbation theory on the S22 data set of noncovalent interactions. The computational cost of the dispersion energy computation is low compared to the self-consistent field calculation of the QM subsystem. The dispersion energy is sensitive to the level of theory employed for the QM part and to the electrostatic interactions in the system. The latter means that the dispersion interactions in the QM/EFP method are not purely two-body but have more complex many-body behavior.

Relativistic ab Initio Accurate Atomic Minimal Basis Sets: Quantitative LUMOs and Oriented Quasi-Atomic Orbitals for the Elements Li-Xe.

Valence virtual orbitals (VVOs) are a quantitative and basis set independent method for extracting chemically meaningful lowest unoccupied molecular orbitals (LUMOs). The VVOs are formed based on a singular value decomposition (SVD) with respect to precomputed and internally stored ab initio accurate atomic minimal basis sets (AAMBS) for the atoms. The occupied molecular orbitals and VVOs together form a minimal basis set that can be transformed into orthogonal oriented quasi-atomic orbitals (OQUAOs) that provide a quantitative description of the bonding in a molecular environment. In the present work, relativistic AAMBS are developed that span the full valence orbital space. The impact of using full valence AAMBS for the formation of the VVOs and OQUAOs and the resulting bonding analysis is demonstrated with applications to the cuprous chloride, scandium monofluoride, and nickel silicide diatomic molecules.

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 C2 molecule.

Correlation Energy Extrapolation by Many-Body Expansion.

Accounting for electron correlation is required for high accuracy calculations of molecular energies. The full configuration interaction (CI) approach can fully capture the electron correlation within a given basis, but it does so at a computational expense that is impractical for all but the smallest chemical systems. In this work, a new methodology is presented to approximate configuration interaction calculations at a reduced computational expense and memory requirement, namely, the correlation energy extrapolation by many-body expansion (CEEMBE). This method combines a MBE approximation of the CI energy with an extrapolated correction obtained from CI calculations using subsets of the virtual orbitals. The extrapolation approach is inspired by, and analogous to, the method of correlation energy extrapolation by intrinsic scaling. Benchmark calculations of the new method are performed on diatomic fluorine and ozone. The method consistently achieves agreement with CI calculations to within a few mhartree and often achieves agreement to within ∼1 millihartree or less, while requiring significantly less computational resources.

The transition from the open minimum to the ring minimum on the ground state and on the lowest excited state of like symmetry in ozone: A configuration interaction study.

The metastable ring structure of the ozone 1(1)A1 ground state, which theoretical calculations have shown to exist, has so far eluded experimental detection. An accurate prediction for the energy difference between this isomer and the lower open structure is therefore of interest, as is a prediction for the isomerization barrier between them, which results from interactions between the lowest two (1)A1 states. In the present work, valence correlated energies of the 1(1)A1 state and the 2(1)A1 state were calculated at the 1(1)A1 open minimum, the 1(1)A1 ring minimum, the transition state between these two minima, the minimum of the 2(1)A1 state, and the conical intersection between the two states. The geometries were determined at the full-valence multi-configuration self-consistent-field level. Configuration interaction (CI) expansions up to quadruple excitations were calculated with triple-zeta atomic basis sets. The CI expansions based on eight different reference configuration spaces were explored. To obtain some of the quadruple excitation energies, the method of Correlation Energy Extrapolation by Intrinsic Scaling was generalized to the simultaneous extrapolation for two states. This extrapolation method was shown to be very accurate. On the other hand, none of the CI expansions were found to have converged to millihartree (mh) accuracy at the quadruple excitation level. The data suggest that convergence to mh accuracy is probably attained at the sextuple excitation level. On the 1(1)A1 state, the present calculations yield the estimates of (ring minimum-open minimum) ∼45-50 mh and (transition state-open minimum) ∼85-90 mh. For the (2(1)A1-(1)A1) excitation energy, the estimate of ∼130-170 mh is found at the open minimum and 270-310 mh at the ring minimum. At the transition state, the difference (2(1)A1-(1)A1) is found to be between 1 and 10 mh. The geometry of the transition state on the 1(1)A1 surface and that of the minimum on the 2(1)A1 surface nearly coincide. More accurate predictions of the energy differences also require CI expansions to at least sextuple excitations with respect to the valence space. For every wave function considered, the omission of the correlations of the 2s oxygen orbitals, which is a widely used approximation, was found to cause errors of about ±10 mh with respect to the energy differences.

A Comprehensive Analysis in Terms of Molecule-Intrinsic, Quasi-Atomic Orbitals. II. Strongly Correlated MCSCF Wave Functions.

A methodology is developed for the quantitative identification of the quasi-atomic orbitals that are embedded in a strongly correlated molecular wave function. The wave function is presumed to be generated from configurations in an internal orbital space whose dimension is equal to (or slightly larger) than that of the molecular minimal basis set. The quasi-atomic orbitals are found to have large overlaps with corresponding orbitals on the free atoms. They separate into bonding and nonbonding orbitals. From the bonding quasi-atomic orbitals, localized bonding and antibonding molecular orbitals are formed. The resolution of molecular density matrices in terms of these orbitals furnishes a basis for analyzing the interatomic bonding patterns in molecules and the changes in these bonding patterns along reaction paths. A new bond strength measure, the kinetic bond order, is introduced.

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.

Seniority number description of potential energy surfaces: Symmetric dissociation of water, N2, C2, and Be2.

The present study further explores the concept of the seniority number (Ω) by examining different configuration interaction (CI) truncation strategies in generating compact wave functions in a systematic way. While the role of Ω in addressing static (strong) correlation problem has been addressed in numerous previous studies, the usefulness of seniority number in describing weak (dynamic) correlation has not been investigated in a systematic way. Thus, the overall objective in the present work is to investigate the role of Ω in addressing also dynamic electron correlation in addition to the static correlation. Two systematic CI truncation strategies are compared beyond minimal basis sets and full valence active spaces. One approach is based on the seniority number (defined as the total number of singly occupied orbitals in a determinant) and another is based on an excitation-level limitation. In addition, molecular orbitals are energy-optimized using multiconfigurational-self-consistent-field procedure for all these wave functions. The test cases include the symmetric dissociation of water (6-31G), N2 (6-31G), C2 (6-31G), and Be2 (cc-pVTZ). We find that the potential energy profile for H2O dissociation can be reasonably well described using only the Ω = 0 sector of the CI wave function. For the Be2 case, we show that the full CI potential energy curve (cc-pVTZ) is almost exactly reproduced using either Ω-based (including configurations having up to Ω = 2 in the virtual-orbital-space) or excitation-based (up to single-plus-double-substitutions) selection methods, both out of a full-valence-reference function. Finally, in dissociation cases of N2 and C2, we shall also consider novel hybrid wave functions obtained by a union of a set of CI configurations representing the full valence space and a set of CI configurations where seniority-number restriction is imposed for a complete set (full-valence-space and virtual) of correlated molecular orbitals, simultaneously. We discuss the usefulness of the seniority number concept in addressing both static and dynamic electron correlation problems along dissociation paths.

A Comprehensive Analysis in Terms of Molecule-Intrinsic Quasi-Atomic Orbitals. IV. Bond Breaking and Bond Forming along the Dissociative Reaction Path of Dioxetane.

The quantitative analysis of molecular density matrices in terms of oriented quasi-atomic orbitals (QUAOs) is shown to yield detailed conceptual insight into the dissociation of dioxetane on the basis of ab initio wave functions. The QUAOs persist and can be followed throughout the reaction path. The kinetic bond orders and the orbital populations of the QUAOs quantitatively reveal the changes of the bonding interactions along the reaction path. At the transition state the OO bond is broken, and the molecule becomes a biradical. After the transition state the reaction path bifurcates. The minimum energy path gently descends from the transition state via a valley-ridge inflection point to a second saddle point, from which two new minimum energy paths lead to two equivalent formaldehyde dimers. The CC bond breaks, and the π-bonds of the formaldehyde fragments form in close vicinity of the second saddle point. The changes of the interactions in this region are elucidated by the analysis of the rearrangements of the QUAOs.