In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.
The energy of a molecule is not the sum of the energies of its atomic components. Our aim is to compute an accurate approximation of this quantity based on an electron pair model, that is to say by using an antisymmetric product of two-electron wave functions, called “geminals”. In this model, the total wave function can also be described in a more practical way, by a set of matrices, one for each geminal. However, without further restrictions, such a model has a factorial computational complexity with the number of electrons and its applicability is therefore limited to small systems.We have introduced generalized orthogonality constraints between geminals to reduce the computational effort to a polynomial cost, without sacrificing the indistinguishability of the electrons. In other words, the electron pairs corresponding to the geminals are not fully distinguishable, so that our model remains complex, in the sense that geminal products have still to be antisymmetrized according to the Pauli principle to form bona fide electronic wave functions.Our geometrical constraints translate, in terms of geminal matrices, into simple equations involving the traces of products of these matrices. In the simplest non-trivial model, a set of solutions is given by block-diagonal matrices where each block is of size 2x2 and consists of a Pauli matrix or the Identity, multiplied by a complex parameter to be optimized. With this simplified ansatz for geminals, the number of terms in the calculation of the matrix elements of quantum observables is considerably reduced.
In previous research, we tested the wave-function format of a linear combination of several antisymmetrized geminal power states. A numerical problem in the geminal matrices was noted, which made the total energies of electronic systems with large numbers of electrons unstable. The underlying cause was found to be the large cancellation term in the geminal power series. We have obtained a format to resolve this problem for the case of total energies and partly for the first-order derivatives within the antisymmetrized geminal power states. By using this formalism, we have calculated the ground-state energies for several electronic systems, including the usage of a larger chemical basis set. The results are, in some cases, very close to the exact result, especially for one-dimensional Hubbard systems. Our result for a water molecule with the Dunning double zeta basis set is better than the configuration-interaction singles and doubles energy and approaches the coupled cluster singles and doubles energy.
The Lewis electron-pair picture pervades all of chemistry and a substantialdomain of materials science. “Geminals” are electron-pair quantum states whichare used to express approximate solutions of the Schrödinger equation of polyelec-tronic systems. In quantum chemistry, geminal-based methods are potentially moreeffective than the traditional ones based on the orbital picture: electrons occupyorbitals and experience only the average effect of each other. This is especially truefor strongly correlated-systems. However, the computational cost of using arbitrarygeminal wave functions is non-polynomial in the system size.It is therefore an important challenge to find additional constraints on geminalmodels which preserve the accuracy of the wavefunction description while reducingthe calculations to polynomial cost. The Ayers’ group at Mc Master has developedone such successful geminal model called AP1roG because the wavefunction is anantisymmetric product of geminals with one distinct reference orbital occupied ineach geminal [1]. This ansatz relaxes the so-called “strong-orthogonality” constraintbetween geminals. The latter makes wavefunction matrix element computationsvery easy but appears too drastic from the physical point of view, since the geminalelectron pairs are then distinguishable.We have followed a different path to relax the strong-orthogonality constraint. Wehave defined the geometrical concept of n-orthogonality [2] which can be viewed asa graded indistinguishability measure for electronic states: 1-orthogonality coincideswith strong orthogonality, the larger n the less distinguishable the n-orthogonal elec-tronic states will be. We have studied geminal models constrained by n-orthogonalityrelations for different values of n, imposed to different combinations of geminal prod-uct states [3]. However, in this work the scaling of matrix element computations wasstill unsatisfactory. We will present new ideas developed since then to remedy thisproblem.References[1] P. A Limacher, P. W. Ayers, P. A. Johnson, S. De Baerdemacker, D. Van Neck and P. Bultinck, J. of Chem. Theory and Computation 9, 1394-1401 (2013).[2] P. Cassam-Chenaı̈, Phys. Rev. A77, 032103 (2008).[3] P. Cassam-Chenaı̈, V. Rassolov, Chem. Phys. Lett. 487, 147-152 (2010).
We propose an approach to the electronic structure problem based on noninteracting electron pairs that has similar computational cost to conventional methods based on noninteracting electrons. In stark contrast to other approaches, the wave function is an antisymmetric product of nonorthogonal geminals, but the geminals are structured so the projected Schrödinger equation can be solved very efficiently. We focus on an approach where, in each geminal, only one of the orbitals in a reference Slater determinant is occupied. The resulting method gives good results for atoms and small molecules. It also performs well for a prototypical example of strongly correlated electronic systems, the hydrogen atom chain.