Density matrices of seniority-zero geminal wavefunctions

Scalar products and density matrix elements of closed-shell pair geminal wavefunctions are evaluated directly in terms of the pair amplitudes, resulting in an analog of Wick’s theorem for fermions or bosons. This expression is, in general, intractable, but it is shown how it becomes feasible in three distinct ways for Richardson–Gaudin (RG) states, the antisymmetrized geminal power, and the antisymmetrized product of strongly orthogonal geminals. Dissociation curves for hydrogen chains are computed with off-shell RG states and the antisymmetrized product of interacting geminals. Both are near exact, suggesting that the incorrect results observed with ground state RG states (a local maximum rather than smooth dissociation) may be fixable using a different RG state.

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement─a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene → benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

Valence Bond Theory-Its Birth, Struggles with Molecular Orbital Theory, Its Present State and Future Prospects

This essay describes the successive births of valence bond (VB) theory during 1916–1931. The alternative molecular orbital (MO) theory was born in the late 1920s. The presence of two seemingly different descriptions of molecules by the two theories led to struggles between the main proponents, Linus Pauling and Robert Mulliken, and their supporters. Until the 1950s, VB theory was dominant, and then it was eclipsed by MO theory. The struggles will be discussed, as well as the new dawn of VB theory, and its future.

Spin-Coupled Generalized Valence Bond Theory: New Perspectves on the Electronic Structure of Molecules and Chemical Bonds

Spin-Coupled Generalized Valence Bond (SCGVB) theory provides the foundation for a comprehensive theory of the electronic structure of molecules. SCGVB theory offers a compelling orbital description of the electronic structure of molecules as well as an efficient and effective zero-order wave function for calculations striving for quantitative predictions of molecular structures, energetics, and other properties. The orbitals in the SCGVB wave function are usually semilocalized, and for most molecules, they can be interpreted using concepts familiar to all chemists (hybrid orbitals, localized bond pairs, lone pairs, etc.). SCGVB theory also provides new perspectives on the nature of the bonds in molecules such as C₂, Be₂ and SF₄/SF₆. SCGVB theory contributes unparalleled insights into the underlying cause of the first-row anomaly in inorganic chemistry as well as the electronic structure of organic molecules and the electronic mechanisms of organic reactions. The SCGVB wave function accounts for nondynamical correlation effects and, thus, corrects the most serious deficiency in molecular orbital (RHF) wave functions. Dynamical correlation effects, which are critical for quantitative predictions, can be taken into account using the SCGVB wave function as the zero-order wave function for multireference configuration interaction or coupled cluster calculations.

Ab initio valence bond theory: A brief history, recent developments, and near future

This Perspective presents a survey of several issues in ab initio valence bond (VB) theory with a primary focus on recent advances made by the Xiamen VB group, including a brief review of the earlier history of the ab initio VB methods, in-depth discussion of algorithms for nonorthogonal orbital optimization in the VB self-consistent field method and VB methods incorporating dynamic electron correlation, along with a concise overview of VB methods for complex systems and VB models for chemical bonding and reactivity, and an outlook of opportunities and challenges for the near future of the VB theory.

Evaluation of full valence correlation energies and gradients

Complete-active-space self-consistent field (CASSCF) wave functions are central to understanding strongly correlated molecules as they capture the entirety of electronic interactions within a subset of the orbital space. The most interesting case for CASSCF is the full valence limit, where all bonding and an equal number of virtual orbitals are included in the active space, and no approximation is made in selecting the important valence orbitals or electrons. While conventional algorithms require exponential computational time to evaluate full valence CASSCF, this article shows that the method of increments can do the same with polynomial effort, in a new method denoted iCASSCF. The method of increments can also provide density matrices and other necessary ingredients for the construction of the nuclear gradient. These goals are met through a many-body expansion that breaks the problem into smaller pieces that are subsequently reassembled to form close approximations of conventional CAS results. Practical demonstrations on a number of medium-sized molecules, with up to 116 valence electrons correlated in 116 orbitals, show the power of this methodology.

Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer

In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO–LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the ground-state energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.

Spin-Coupled Generalized Valence Bond Description of Group 14 Species: The Carbon, Silicon and Germanium Hydrides, XH n ( n = 1-4)

Although elements in the same group in the Periodic Table tend to behave similarly, the differences in the simplest Group 14 hydrides-XH _n (X = C, Si, Ge; n = 1-4)-are as pronounced as their similarities. Spin-coupled generalized valence bond (SCGVB) as well as coupled cluster [CCSD(T)] calculations are reported for all of the molecules in the CH _n/SiH _n/GeH _n series to gain insights into the factors underlying these differences. It is found that the relative weakness of the recoupled pair bonds of SiH and GeH gives rise to the observed differences in the ground state multiplicities, molecular structures, and bond energies of SiH _n and GeH _n. A number of factors that influence the strength of the recoupled pair bonds in CH, SiH, and GeH were examined. Two factors were identified as potential contributors to the decrease in the strengths of these bonds from CH to SiH and GeH: (i) a decrease in the overlap between the orbitals involved in the bond and (ii) an increase in Pauli repulsion between the electrons in the two lobe orbitals centered on the X atoms. Finally, an analysis of the hybridization of the SCGVB orbitals in XH₄ indicates that they are closer to sp hybrids than sp³ hybrids, which implies that the underlying cause of the tetrahedral structure of the XH₄ molecules is not a direct result of the hybridization of the X atom orbitals.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Extending density functional embedding theory for covalently bonded systems

Quantum embedding theory aims to provide an efficient solution to obtain accurate electronic energies for systems too large for full-scale, high-level quantum calculations. It adopts a hierarchical approach that divides the total system into a small embedded region and a larger environment, using different levels of theory to describe each part. Previously, we developed a density-based quantum embedding theory called density functional embedding theory (DFET), which achieved considerable success in metals and semiconductors. In this work, we extend DFET into a density-matrix-based nonlocal form, enabling DFET to study the stronger quantum couplings between covalently bonded subsystems. We name this theory density-matrix functional embedding theory (DMFET), and we demonstrate its performance in several test examples that resemble various real applications in both chemistry and biochemistry. DMFET gives excellent results in all cases tested thus far, including predicting isomerization energies, proton transfer energies, and highest occupied molecular orbital-lowest unoccupied molecular orbital gaps for local chromophores. Here, we show that DMFET systematically improves the quality of the results compared with the widely used state-of-the-art methods, such as the simple capped cluster model or the widely used ONIOM method.

The Quadruple Bonding in C2 Reproduces the Properties of the Molecule

Ever since Lewis depicted the triple bond for acetylene, triple bonding has been considered as the highest limit of multiple bonding for main elements. Here we show that C2 is bonded by a quadruple bond that can be distinctly characterized by valence-bond (VB) calculations. We demonstrate that the quadruply-bonded structure determines the key observables of the molecule, and accounts by itself for about 90% of the molecule's bond dissociation energy, and for its bond lengths and its force constant. The quadruply-bonded structure is made of two strong π bonds, one strong σ bond and a weaker fourth σ-type bond, the bond strength of which is estimated as 17-21 kcal mol(-1). Alternative VB structures with double bonds; either two π bonds or one π bond and one σ bond lie at 129.5 and 106.1 kcal mol(-1), respectively, above the quadruply-bonded structure, and they collapse to the latter structure given freedom to improve their double bonding by dative σ bonding. The usefulness of the quadruply-bonded model is underscored by "predicting" the properties of the (3)Σ+u state. C2's very high reactivity is rooted in its fourth weak bond. Thus, carbon and first-row main elements are open to quadruple bonding!

Description of Polar Chemical Bonds from the Quantum Mechanical Interference Perspective

The Generalized Product Function Energy Partitioning (GPF-EP) method has been applied to a set of molecules, AH (A = Li, Be, B, C, N, O, F), CO and LiF with quite different dipole moments, in order to investigate the role played by the quantum interference effect in the formation of polar chemical bonds. The calculations were carried out with GPF wave functions treating all the core electrons as a single Hartree-Fock group and the bonding electrons at the Generalized Valence Bond Perfect-Pairing (GVB-PP) level, with the cc-pVTZ basis set. The results of the energy partitioning into interference and quasi-classical contributions along the respective Potential Energy Surfaces (PES) show that the main contribution to the depth of the potential wells comes from the interference term, which is an indication that all the molecules mentioned above form typical covalent bonds. In all cases, the stabilization promoted by the interference term comes from the kinetic contribution, in agreement with previous results. The analysis of the effect of quantum interference on the electron density reveals that while polarization effects (quasi-classical) tend to displace electronic density from the most polarizable atom toward the less polarizable one, interference (quantum effects) counteracts by displacing electronic density to the bond region, giving rise to the right electronic density and dipole moment.

Active Thermochemical Tables: Sequential Bond Dissociation Enthalpies of Methane, Ethane, and Methanol and the Related Thermochemistry

Active Thermochemical Tables (ATcT) thermochemistry for the sequential bond dissociations of methane, ethane, and methanol systems were obtained by analyzing and solving a very large thermochemical network (TN). Values for all possible C-H, C-C, C-O, and O-H bond dissociation enthalpies at 298.15 K (BDE298) and bond dissociation energies at 0 K (D0) are presented. The corresponding ATcT standard gas-phase enthalpies of formation of the resulting CHn, n = 4-0 species (methane, methyl, methylene, methylidyne, and carbon atom), C2Hn, n = 6-0 species (ethane, ethyl, ethylene, ethylidene, vinyl, ethylidyne, acetylene, vinylidene, ethynyl, and ethynylene), and COHn, n = 4-0 species (methanol, hydroxymethyl, methoxy, formaldehyde, hydroxymethylene, formyl, isoformyl, and carbon monoxide) are also presented. The ATcT thermochemistry of carbon dioxide, water, hydroxyl, and carbon, oxygen, and hydrogen atoms is also included, together with the sequential BDEs of CO2 and H2O. The provenances of the ATcT enthalpies of formation, which are quite distributed and involve a large number of relevant determinations, are analyzed by variance decomposition and discussed in terms of principal contributions. The underlying reasons for periodic appearances of remarkably low and/or unusually high BDEs, alternating along the dissociation sequences, are analyzed and quantitatively rationalized. The present ATcT results are the most accurate thermochemical values currently available for these species.

Stereoelectronic switching in single-molecule junctions

A new intersection between reaction chemistry and electronic circuitry is emerging from the ultraminiaturization of electronic devices. Over decades chemists have developed a nuanced understanding of stereoelectronics to establish how the electronic properties of molecules relate to their conformation; the recent advent of single-molecule break-junction techniques provides the means to alter this conformation with a level of control previously unimagined. Here we unite these ideas by demonstrating the first single-molecule switch that operates through a stereoelectronic effect. We demonstrate this behaviour in permethyloligosilanes with methylthiomethyl electrode linkers. The strong σ conjugation in the oligosilane backbone couples the stereoelectronic properties of the sulfur-methylene σ bonds that terminate the molecule. Theoretical calculations support the existence of three distinct dihedral conformations that differ drastically in their electronic character. We can shift between these three species by simply lengthening or compressing the molecular junction, and, in doing so, we can switch conductance digitally between two states.

The first row anomaly and recoupled pair bonding in the halides of the late p-block elements

The dramatic differences between the properties of molecules formed from the late p-block elements of the first row of the periodic table (N–F) and those of the corresponding elements in subsequent rows is well recognized as the first row anomaly. Certain properties of the atoms, such as the relative energies and spatial extents of the ns and np orbitals, can explain some of these differences, but not others.

In this Account, we summarize the results of our recent computational studies of the halides of the late p-block elements. Our studies point to a single underlying cause for many of these differences: the ability of the late p-block elements in the second and subsequent rows of the periodic table to form recoupled pair bonds and recoupled pair bond dyads with very electronegative ligands. Recoupled pair bonds form when an electron in a singly occupied ligand orbital recouples the pair of electrons in a doubly occupied lone pair orbital on the central atom, leading to a central atom-ligand bond. Recoupled pair bond dyads occur when a second ligand forms a bond with the orbital left over from the initial recoupled pair bond.

Recoupled pair bonds and recoupled pair bond dyads enable the late p-block elements to form remarkably stable hypervalent compounds such as PF₅ and SF₆ and lead to unexpected excited states in smaller halides of the late p-block elements such as SF and SF₂. Recoupled pair bonding also causes the F_n–1X–F bond energies to oscillate dramatically once the normal valences of the central atoms have been satisfied. In addition, recoupled pair bonding provides a lower-energy pathway for inversion in heavily fluorinated compounds (PF₃ and PF₂H, but not PH₂F and PH₃) and leads to unusual intermediates and products in reactions involving halogens and late p-block element compounds, such as (CH₃)₂S + F₂.

Although this Account focuses on the halides of the second row, late p-block elements, recoupled pair bonds and recoupled pair bond dyads are important in the chemistry of p-block elements beyond the second row (As, Se, and Br) and for compounds of these elements with other very electronegative ligands, such as OH and O. Knowledge of recoupled pair bonding is thus critical to understanding the properties and reactivity of molecules containing the late p-block elements beyond the first row.

Quadruple bonding in C2 and analogous eight-valence electron species

Triple bonding is conventionally considered to be the limit for multiply bonded main group elements, despite higher metal-metal bond orders being frequently observed for transition metals and lanthanides/actinides. Here, using high-level theoretical methods, we show that C(2) and its isoelectronic molecules CN(+), BN and CB(-) (each having eight valence electrons) are bound by a quadruple bond. The bonding comprises not only one σ- and two π-bonds, but also one weak 'inverted' bond, which can be characterized by the interaction of electrons in two outwardly pointing sp hybrid orbitals. A simple way of assessing the energy of the fourth bond is proposed and is found to be ~12-17 kcal mol(-1) for the isoelectronic species studied, and thus stronger than a hydrogen bond. In contrast, the analogues of C(2) that contain higher-row elements, such as Si(2) and Ge(2), exhibit only double bonding.

Chemical bonding and the equilibrium composition of Grignard reagents in ethereal solutions

A thorough analysis of the electronic structure and thermodynamic aspects of Grignard reagents and its associated equilibrium composition in ethereal solutions is performed. Considering methylmagnesium halides containing fluorine, chlorine, and bromine, we studied the neutral, charged, and radical species associated with their chemical equilibrium in solution. The ethereal solvents considered, tetrahydrofuran (THF) and ethyl ether (Et(2)O), were modeled using the polarizable continuum model (PCM) and also by explicit coordination to the Mg atoms in a cluster. The chemical bonding of the species that constitute the Grignard reagent is analyzed in detail with generalized valence bond (GVB) wave functions. Equilibrium constants were calculated with the DFT/M06 functional and GVB wave functions, yielding similar results. According to our calculations and existing kinetic and electrochemical evidence, the species R(•), R(-), (•)MgX, and RMgX(2)(-) must be present in low concentration in the equilibrium. We conclude that depending on the halogen, a different route must be followed to produce the relevant equilibrium species in each case. Chloride and bromide must preferably follow a "radical-based" pathway, and fluoride must follow a "carbanionic-based" pathway. These different mechanisms are contrasted against the available experimental results and are proven to be consistent with the existing thermodynamic data on the Grignard reagent equilibria.