6.18.1 Introduction

The methods described below are related to valence bond theory and are handled by the GVBMAN module. The following models are available:

Molecules where electron correlation is strong are characterized by small energy gaps between the nominally occupied orbitals (that would comprise the Hartree-Fock wave function, for example) and nominally empty orbitals. Examples include so-called diradicaloid molecules, ⁶¹⁸ Jung Y., Head-Gordon M.
ChemPhysChem
(2003), 4, pp. 522. Link or molecules with partly broken chemical bonds (as in some transition-state structures). Because the energy gap is small, electron configurations other than the reference determinant contribute to the molecular wave function with considerable amplitude, and omitting them leads to a significant error.

Including all possible configurations however, is a vast overkill. It is common to restrict the configurations that one generates to be constructed not from all molecular orbitals, but just from orbitals that are either “core” or “active”. In this section, we consider just one type of active space, which is composed of two orbitals to represent each electron pair: one nominally occupied (bonding or lone pair in character) and the other nominally empty, or correlating (it is typically anti-bonding in character). This is usually called the perfect pairing active space, and it clearly is well-suited to represent the bonding/anti-bonding correlations that are associated with bond-breaking.

The quantum chemistry within this (or any other) active space is given by a Complete Active Space SCF (CASSCF) calculation, whose exponential cost growth with molecule size makes it prohibitive for systems with more than about 14 active orbitals. One well-defined coupled cluster (CC) approximation based on CASSCF is to include only double substitutions in the valence space whose orbitals are then optimized. In the framework of conventional CC theory, this defines the valence optimized doubles (VOD) model, ⁶⁸⁶ Krylov A. I. et al.
J. Chem. Phys.
(1998), 109, pp. 10669. Link which scales as $𝒪(N^6)$ (see Section 6.13.2). This is still too expensive to be readily applied to large molecules.

The methods described in this section bridge the gap between sophisticated but expensive coupled cluster methods and inexpensive methods such as DFT, HF and MP2 theory that may be (and indeed often are) inadequate for describing molecules that exhibit strong electron correlations such as diradicals. The coupled cluster perfect pairing (PP), ²⁷⁰ Cullen J.
Chem. Phys.
(1996), 202, pp. 217. Link ^{, 103} Beran G. J. O. et al.
J. Phys. Chem. A
(2005), 109, pp. 9183. Link imperfect pairing ¹²⁹³ Van Voorhis T., Head-Gordon M.
Chem. Phys. Lett.
(2000), 317, pp. 575. Link (IP) and restricted coupled cluster ¹²⁹⁶ Van Voorhis T., Head-Gordon M.
J. Chem. Phys.
(2001), 115, pp. 7814. Link (RCC) models are local approximations to VOD that include only a linear and quadratic number of double substitution amplitudes respectively. They are close in spirit to generalized valence bond (GVB)-type wave functions, ⁴²⁹ Goddard III W. A., Harding L. B.
Annu. Rev. Phys. Chem.
(1978), 29, pp. 363. Link because in fact they are all coupled cluster models for GVB that share the same perfect pairing active space. The most powerful method in the family, the Coupled Cluster Valence Bond (CCVB) method, ¹¹⁸⁶ Small D. W., Head-Gordon M.
J. Chem. Phys.
(2009), 130, pp. 084103. Link ^{, 1187} Small D. W., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2011), 13, pp. 19285. Link ^{, 1188} Small D. W., Head-Gordon M.
J. Chem. Phys.
(2012), 137, pp. 114103. Link is a valence bond approach that goes well beyond the power of GVB-PP and related methods, as discussed below in Section 6.18.3.