Correlation‐consistent configuration interaction: Accurate bond dissociation energies from simple wave functions

A configuration interaction correction on top of pair coupled cluster doubles

Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit a posteriori corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the N₂ and F₂ dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between those of the linearized frozen pCCD and frozen pCCD variants.

Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry

In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.

Describing Strong Correlation with Block-Correlated Coupled Cluster Theory

A block-correlated coupled cluster (BCCC) method based on the generalized valence bond (GVB) wave function (GVB-BCCC in short) is proposed and implemented at the ab initio level, which represents an attractive multireference electronic structure method for strongly correlated systems. The GVB-BCCC method is demonstrated to provide satisfactory descriptions for multiple bond breaking in small molecules, although the GVB reference function is qualitatively wrong for the studied processes. For a challenging prototype of strongly correlated systems, tridecane with all 12 single C-C bonds at various distances, our calculations have shown that the GVB-BCCC2b method can provide highly comparable results as the density matrix renormalization group method for potential energy surfaces along simultaneous dissociation of all C-C bonds.

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.

Orbitals and the Interpretation of Photoelectron Spectroscopy and (e,2e) Ionization Experiments

Electron momentum spectroscopy, scanning tunneling microscopy, and photoelectron spectroscopy provide unique information about electronic structure, but their interpretation has been controversial. This essay discusses a framework for interpretation. Although this interpretation is not new, we believe it is important to present this framework in light of recent publications. The key point is that these experiments provide information about how the electron distribution changes upon ionization, not how electrons behave in the pre-ionized state. Therefore, these experiments do not lead to a "selection of the correct orbitals" in chemistry and do not overturn the well-known conclusion that both delocalized molecular orbitals and localized molecular orbitals are useful for interpreting chemical structure and dynamics. The two types of orbitals can produce identical total molecular electron densities and therefore molecular properties. Different types of orbitals are useful for different purposes.

Description of ground and excited electronic states by ensemble density functional method with extended active space

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units. We demonstrate that the new method correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena.

Using the GVB Ansatz to develop ensemble DFT method for describing multiple strongly correlated electron pairs

Ensemble density functional theory (DFT) furnishes a rigorous theoretical framework for describing the non-dynamic electron correlation arising from (near) degeneracy of several electronic configurations. Ensemble DFT naturally leads to fractional occupation numbers (FONs) for several Kohn-Sham (KS) orbitals, which thereby become variational parameters of the methodology. The currently available implementation of ensemble DFT in the form of the spin-restricted ensemble-referenced KS (REKS) method was originally designed for systems with only two fractionally occupied KS orbitals, which was sufficient to accurately describe dissociation of a single chemical bond or the singlet ground state of biradicaloid species. To extend applicability of the method to systems with several dissociating bonds or to polyradical species, more fractionally occupied orbitals must be included in the ensemble description. Here we investigate a possibility of developing the extended REKS methodology with the help of the generalized valence bond (GVB) wavefunction theory. The use of GVB enables one to derive a simple and physically transparent energy expression depending explicitly on the FONs of several KS orbitals. In this way, a version of the REKS method with four electrons in four fractionally occupied orbitals is derived and its accuracy in the calculation of various types of strongly correlated molecules is investigated. We propose a possible scheme to ameliorate the partial size-inconsistency that results from perfect spin-pairing. We conjecture that perfect pairing natural orbital (NO) functionals of reduced density matrix functional theory (RDMFT) should also display partial size-inconsistency.

Adaptive multiconfigurational wave functions

A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N2 and the potential energy curves for the first three singlet states of C2. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu2O2(2+) core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.

Formation of Carbon-Carbon Triply Bonded Molecules from Two Free Carbyne Radicals via a Conical Intersection

The recent proposal (Bogoslavsky, B.; Levy, O.; Kotlyar, A.; Salem, M.; Gelman, F.; Bino, A. Angew. Chem., Int. Ed.2012, 51, 90-94) that metallo-alkylidyne complexes decompose in aqueous solution and give rise to free carbynes, which couple to yield acetylenes, is examined here theoretically. On the basis of the known marker reactions of carbynes in the doublet and quartet state, it is concluded that most of the reactivity patterns observed in the Bino experiment arose from quartet carbynes. Indeed, theory shows that quartet carbynes can be funneled to acetylene via a conical intersection. Moreover, many of the minor products are also identified as markers of the quartet carbynes. Carbynes formation in their doublet state is a minor channel that branches from the conical intersection and leads to the formation of dienes and olefins in the Bino experiment. Thus, we show that conical intersections are important also in thermally initiated reactions. Coupled to the experimental approach, the study opens a window to studies of carbynes under mild conditions.

An approximate diatomics in molecules formulation of generalized valence bond theory

The slow computational speed of the generalized valence bond perfect pairing method (GVB-PP) has been an impediment to its routine use. We have addressed this problem by employing a diatomics in molecules Hamiltonian derived from a second quantization perturbation approach. This results in all three- and four-centered two-electron integrals being dropped from the traditional GVB-PP calculation. For moderate sized molecules, as for example C20 computed with a double zeta + polarization basis, there is on average a fifty-fold decrease in computational times. In this article, we present the theory behind our approach and analyze the accuracy and speed of this approximate GVB-PP method for several cases where density functional methods have produced ambivalent results.

Second-order correction to perfect pairing: An inexpensive electronic structure method for the treatment of strong electron-electron correlations

We have formulated a second-order perturbative correction for perfect-pairing wave functions [PP2] based on similarity-transformed perturbation techniques in coupled cluster theory. The perfect-pairing approximation is used to obtain a simple reference wave function which can qualitatively describe bond breaking, diradicals, and other highly correlated systems, and the perturbative correction accounts for the dynamical correlation. An efficient implementation of this correction using the resolution of the identity approximation enables PP2 to be computed at a cost only a few times larger than that of canonical MP2 for systems with hundreds of active electrons and tens of heavy atoms. PP2 significantly improves on MP2 predictions in various systems with a challenging electronic structure.

Bond electron pair: Its relevance and analysis from the quantum chemistry point of view

This paper first comments on the surprisingly poor status that Quantum Chemistry has offered to the fantastic intuition of Lewis concerning the distribution of the electrons in the molecule. Then, it advocates in favor of a hierarchical description of the molecular wave-function, distinguishing the physics taking place in the valence space (in the bond and between the bonds), and the dynamical correlation effects. It is argued that the clearest pictures of the valence electronic population combine two localized views, namely the bond (and lone pair) Molecular Orbitals and the Valence Bond decomposition of the wave-function, preferably in the orthogonal version directly accessible from the complete active space self consistent field method. Such a reading of the wave function enables one to understand the work of the nondynamical correlation as an enhancement of the weight of the low-energy VB components, i.e. as a better compromise between the electronic delocalization and the energetic preferences of the atoms. It is suggested that regarding the bond building, the leading dynamical correlation effect may be the dynamical polarization phenomenon. It is shown that most correlation effects do not destroy the bond electron pairs and remain compatible with Lewis' vision. A certain number of free epistemological considerations have been introduced in the development of the argument.

A survey of recent developments in ab initio valence bond theory

Starting from the 1980s and onwards, Valence Bond theory has been enjoying renaissance that is characterized by the development of a growing number of ab initio methods, and by many applications to chemical reactivity and to the central paradigms of chemistry. Owing the increase of computational power of modern computers and to significant advances in the methodology, valence bond theory begins to offer a sound and attractive alternative to Molecular Orbital theory. This review aims at summarizing the most important developments of ab initio valence bond methods during the last two or three decades, and is primarily devoted to a description of what the various methods can actually achieve within their specific scopes and limitations. Key available softwares are surveyed.

Size extensive modification of local multireference configuration interaction

We recently developed a reduced scaling multireference configuration interaction (MRCI) method based on local correlation in the internal (occupied) and external (virtual) orbital spaces. This technique can be used, e.g., to predict bond dissociation energies in large molecules with reasonable accuracy. However, the inherent lack of size extensivity of truncated CI is a disadvantage that in principle worsens as the system size grows. Here we implement an a priori size-extensive modification of local MRCI known as the averaged coupled pair functional (ACPF) method. We demonstrate that local MR-ACPF recovers more correlation energy than local MRCI, in keeping with trends observed previously for nonlocal ACPF. We test the size extensivity of local ACPF on noninteracting He atoms and a series of hydrocarbons. Basis set and core correlation effects are explored, as well as bond breaking in a variety of organic molecules. The local MR-ACPF method proves to be a useful tool for investigating large molecules and represents a further improvement in predictive accuracy over local MRCI.