Computational models that use single reference wave function describe molecules
in terms of independent electrons interacting via mean Coulomb and exchange
fields. It is natural to improve this description by using correlated electron
pairs, or geminals, as building blocks for molecular wave functions.
Requirements of computational efficiency and size consistency constrain
geminals to have ,
J. Chem. Phys.
(2002), 117, pp. 5978. with each geminal spanning its own subspace of molecular orbitals. 43 J. Chem. Phys.
(1960), 33, pp. 95. Geminal wave functions were introduced into computational chemistry by Hurley, Lennard-Jones, and Pople. 496 Proc. Roy. Soc. London A
(1953), 220, pp. 446. An excellent review of the history and properties of geminal wave functions is given by Surjan. 1101 Topics Curr. Chem.
(1999), 203, pp. 63.
We implemented a size consistent model chemistry based on Singlet type Strongly orthogonal Geminals (SSG). In SSG, the number of molecular orbitals in each singlet electron pair is an adjustable parameter chosen to minimize total energy. Open-shell orbitals remain uncorrelated. The SSG wave function is computed by setting SSG $rem variable to 1. Both spin-restricted (RSSG) and spin-unrestricted (USSG) versions are available, chosen by the UNRESTRICTED $rem variable.
The wave function has the form
with the coefficients , , and subspaces chosen to minimize the energy
evaluated with the exact Hamiltonian . A constraint for all MO coefficients yields a spin-restricted version of SSG.
SSG model can use any orbital-based initial guess. It is often advantageous to compute Hartree-Fock orbitals and then read them as initial guess for SSG. The program distinguishes Hartree-Fock and SSG initial guess wave functions, and in former case makes preliminary assignment of individual orbital pairs into geminals. The verification of orbital assignments is performed every ten wave function optimization steps, and the orbital pair is reassigned if total energy is lowered.
The convergence algorithm consists of combination of three types of
minimization steps. Direct minimization steps
J. Chem. Phys.
(1976), 65, pp. 265. seek a minimum along the gradient direction, rescaled by the quantity analogous to the orbital energy differences in SCF theory. 953 J. Chem. Phys.
(2002), 117, pp. 5978. If the orbitals are nearly degenerate or inverted, a perturbative re-optimization of single geminal is performed. Finally, new set of the coefficients and is formed from a linear combination of previous iterations, in a manner similar to DIIS algorithm. 938 Chem. Phys. Lett.
(1980), 73, pp. 393. , 939 J. Comput. Chem.
(1982), 3, pp. 556. The size of iterative subspace is controlled by the DIIS_SUBSPACE_SIZE keyword.
After convergence is achieved, SSG reorders geminals based on geminal energy. The energy, along with geminal expansion coefficients, is printed for each geminal. Presence of any but the leading coefficient with large absolute value (value of 0.1 is often used for the definition of “large”) indicates the importance of electron correlation in the system. The Mulliken population analysis is also performed for each geminal, which enables easy assignment of geminals into such chemical objects as core electron pairs, chemical bonds, and lone electron pairs.
As an example, consider the sample calculation of ScH molecule with 6-31G basis
set at the experimental bond distance of 1.776 Å. In its singlet ground
state the molecule has 11 geminals. Nine of them form core electrons on Sc.
Two remaining geminals are:
Geminal 10 E = -1.342609
0.99128 -0.12578 -0.03563 -0.01149 -0.01133 -0.00398
Geminal 11 E = -0.757086
0.96142 -0.17446 -0.16872 -0.12414 -0.03187 -0.01227 -0.01204 -0.00435 -0.00416 -0.00098
Mulliken population analysis shows that geminal 10 is delocalized between Sc and H, indicating a bond. It is moderately correlated, with second expansion coefficient of a magnitude 0.126. The geminal of highest energy is localized on Sc. It represents electrons and describes their excitation into orbitals. Presence of three large expansion coefficients show that this effect cannot be described within GVB framework.