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,
(2003), 4, pp. 522. 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,
J. Chem. Phys.
(1998), 109, pp. 10669. which scales as (see Section 6.12.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),
(1996), 202, pp. 217. , 87 J. Phys. Chem. A
(2005), 109, pp. 9183. imperfect pairing 1114 Chem. Phys. Lett.
(2000), 317, pp. 575. (IP) and restricted coupled cluster 1117 J. Chem. Phys.
(2001), 115, pp. 7814. (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, 368 Annu. Rev. Phys. Chem.
(1978), 29, pp. 363. 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, 1015 J. Chem. Phys.
(2009), 130, pp. 084103. , 1016 Phys. Chem. Chem. Phys.
(2011), 13, pp. 19285. , 1017 J. Chem. Phys.
(2012), 137, pp. 114103. is a valence bond approach that goes well beyond the power of GVB-PP and related methods, as discussed below in Sec. 6.18.2.