Gaussian Geminals in Coupled Cluster and Many-Body Perturbation Theories. In: Rychlewski J, editor. Explicitly Correlated Wave Functions in Chemistry and Physics. Dordrecht: Springer Netherlands; 2003. pp. 185-248. [DOI: 10.1007/978-94-017-0313-0_4]

Second order explicitly correlated R12 theory revisited: A second quantization framework for treatment of the operators’ partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.

Accurate Pair Interaction Energies for Helium from Supermolecular Gaussian Geminal Calculations

Nonrelativistic clamped-nuclei pair interaction energy for ground-state helium atoms has been computed for 12 interatomic separations ranging from 3.0 to 9.0 bohr. The calculations applied the supermolecular approach. The major part of the interaction energy was obtained using the Gaussian geminal implementation of the coupled-cluster theory with double excitations (CCD). Relatively small contributions from single, triple, and quadruple excitations were subsequently included employing the conventional orbital coupled-cluster method with single, double, and noniterative triple excitations [CCSD(T)] and the full configuration interaction (FCI) method. For three distances, the single-excitation contribution was taken from literature Gaussian-geminal calculations at the CCSD level. The orbital CCSD(T) and FCI calculations used very large basis sets, up to doubly augmented septuple- and sextuple-zeta size, respectively, and were followed by extrapolations to the complete basis set limits. The accuracy of the total interaction energies has been estimated to be about 3 mK or 0.03% at the minimum of the potential well. For the attractive part of the well, the relative errors remain consistently smaller than 0.03%. In the repulsive part, the accuracy is even better, except, of course, for the region where the potential goes through zero. For interatomic separations smaller than 4.0 bohr, the relative errors do not exceed 0.01%. Such uncertainties are significantly smaller than the expected values of the relativistic and diagonal Born-Oppenheimer contributions to the potential.

Application of Gaussian-type geminals in local second-order Møller-Plesset perturbation theory

In this work Gaussian-type Geminals (GTGs) are applied in local second-order Moller-Plesset perturbation theory to improve the basis set convergence. Our implementation is based on the weak orthogonality functional of Szalewicz et al., [Chem. Phys. Lett. 91, 169 (1982); J. Chem. Phys. 78, 1420 (1983)] and a newly developed program for calculating the necessary many-electron integrals. The local approximations together with GTGs in the treatment of the correlation energy are introduced and tested. First results for correlation energies of H(2)O, CH(4), CO, C(2)H(2), C(2)H(4), H(2)CO, and N(2)H(4) as well as some reaction and activation energies are presented. More than 97% of the valence-shell correlation energy is recovered using aug-cc-pVDZ basis sets and six GTGs per electron pair. The results are compared with conventional calculations using correlation-consistent basis sets as well as with MP2-R12 results.