Model Wavefunction Forms to Describe Strong Correlation in Quantum Chemistry

Abstract

Standard approaches to approximating solutions to the Schroedinger equation rely upon the assumption that there is a single dominant electronic configuration with a small number of low-lying excited states. As such, they are referred to as single reference approaches and are accurate only when the band gap is small in comparison with the inter-electronic repulsion of the valence electrons. Systems for which this is the case are weakly-correlated. If the inter-electronic repulsion larger than the band gap, then a single electronic reference is not an appropriate ansatz for the physical wavefunction. As the strength of the interaction is increased many Slater determinants become near-degenerate, each providing non-negligible contributions. Such systems are strongly-correlated. The goal of the thesis was to employ wavefunctions which are: non-trivial, tractable, size-consistent and extensive, and not given in an active space. We have outlined approaches based on wavefunction forms more suited to describe strong correlation. Built on the Lie algebra su(2), our first trial wavefunction was the Bethe ansatz solution to the reduced BCS Hamiltonian derived by Richardson. This wavefunction is expressed as a product of quasi-particles each parametrized in terms of quasi-momenta. Consistency among the quasi-momenta is ensured by a set of non-linear equations to be solved numerically. This approach is feasible in both variational and projected flavours. We next considered a wavefunction of the same form without requiring it to be an eigenvector of the reduced BCS Hamiltonian. A variational approach in this case is only feasible provided there are two equivalent descriptions in terms of particles and holes. In any case, a bivariational principle and a projected approach are both feasible. Also based on su(2), there is AP1roG (coupled cluster pair doubles). Open-shell singlets and triplets were constructed with the Lie algebras sp(N) and so(4N). The resulting formulae appear to be infeasible unless it is possible to efficiently compute generalizations of determinants to three-dimensional objects. Treating each spatial orbital individually, all four possible occupations close the Lie superalgebra gl(2|2).