The minimum parameterization of the wave function is derived for the time-independent many-body problem. In this parameterization all electrons and electron pairs are indistinguishable which eliminates the exponential scaling with the number of particles plaguing all other correlation methods. An exponential scaling shown to stem from the expansion of the wave function in one-particle basis sets. It is demonstrated that using a geminal basis with a skew-symmetric orthogonal coefficient matrix, which fulfill a Lie algebra, the parameterization of the exact wave function becomes independent of the number of particles and only scale quadratic with the number of basis functions in the optimized basis and quartic in an arbitrary basis. The resulting antisymmetrized geminal power wave function is shown to fulfill the necessary and sufficient conditions for the exact wave function, treat all electrons and electron pairs as indistinguishable and be the most compact representation of the exact wave function possible. The number of parameters in the antisymmetrized geminal power wave function is exactly the same number of parameters as in the two-particle reduced density matrix thereby showing that these representations are equivalent just like the representations for a one-particle Hamiltonian. A detailed discussion of how the most compact wave function can be derived in general is also presented along with different geminal ansatze and sufficient closure relations. Due to the breaking of the scaling wall for the exact wave function it is expected that, if a general set of constraints can be found, even systems of biological relevance can be treated exactly in the future.
The minimum parameterization of the wave function is derived for the time-independent many-body problem. In this parameterization all electrons and electron pairs are indistinguishable which eliminates the exponential scaling with the number of particles plaguing all other correlation methods. An exponential scaling shown to stem from the expansion of the wave function in one-particle basis sets. It is demonstrated that using a geminal basis with a skew-symmetric orthogonal coefficient matrix, which fulfill a Lie algebra, the parameterization of the exact wave function becomes independent of the number of particles and only scale quadratic with the number of basis functions in the optimized basis and quartic in an arbitrary basis. The resulting antisymmetrized geminal power wave function is shown to fulfill the necessary and sufficient conditions for the exact wave function, treat all electrons and electron pairs as indistinguishable and be the most compact representation of the exact wave function possible. The number of parameters in the antisymmetrized geminal power wave function is exactly the same number of parameters as in the two-particle reduced density matrix thereby showing that these representations are equivalent just like the representations for a one-particle Hamiltonian. A detailed discussion of how the most compact wave function can be derived in general is also presented along with different geminal ansatze and sufficient closure relations. Due to the breaking of the scaling wall for the exact wave function it is expected that, if a general set of constraints can be found, even systems of biological relevance can be treated exactly in the future.
It is virtually impossible to directly solve the Schrödinger equation for a many-electron wave function due to the exponential growth in degrees of freedom with increasing particle number. The two-body reduced density matrix (2-RDM) formalism reduces this coordinate dependence to that of four particles irrespective of the wave function's dimensionality, providing a promising path to solve the many-body problem. Unfortunately, errors arise in this approach because the 2-RDM cannot practically be constrained to guarantee that it corresponds to a valid wave function. Here we approach this so-called N-representability problem by expanding the 2-RDM in a complete basis of two-electron wave functions and studying the matrix formed by the expansion coefficients. This quantity, which we call the geminal density matrix (GDM), is found to evolve in time by a unitary transformation that preserves N-representability. This evolution law enables us to calculate eigenstates of strongly correlated systems by a fictitious adiabatic evolution in which the electron-electron interaction is slowly switched on. We show how this technique is used to diagonalize atomic Hamiltonians, finding that the problem reduces to the solution of ∼N(N−1)/2 two-electron eigenstates of the Helium atom on a grid of electron-electron interaction scaling factors.
Recent studies have considered the possibility that the exact ground-state wavefunction from any Hamiltonian with two-particle interactions may be generated from a single finite two-body exponential transformation acting on an arbitrary Slater determinant [Piecuch et al., Phys. Rev. Lett. 90, 113001 (2003)]. Using the Campbell-Baker-Hausdorff relation, we show that it is difficult for the variational minimum of this trial wave function to satisfy the contracted Schr\"odinger equation which is a necessary and sufficient condition for the wave function to satisfy the Schr\"odinger equation. A counterexample is presented through the Lipkin quasispin model with 4--50 fermions. When the number of fermions exceeds four, the wave function from a finite two-body exponential transformation is shown to be inexact. If the trial wave function ansatz is extended to include products of finite two-body exponential transformations acting on an arbitrary Slater-determinant reference, then we show that the ansatz includes the exact ground-state wave function from any Hamiltonian with only two-particle interactions. Connections between the two-body exponential transformation of the wave function and recent research on two-body exponential similarity transformations of the Hamiltonian [S.R. White, J. Chem. Phys. 117, 7472 (2002)] are discussed.
The energy of any atom or molecule with an even number N of electrons is shown to be an exact functional of a single geminal where the functionals for both the kinetic energy and the external potential are explicitly known. We derive the foundations for geminal functional theory (GFT) through a generalized constrained search and the use of two theorems which demonstrate that all one-particle properties of atoms and molecules with even N may be parametrized by a single geminal [A. J. Coleman, Int. J. Quantum Chem. 63, 23 (1997); D. W. Smith, Phys. Rev. 147, 896 (1966)]. By generalizing constrained search to optimize the universal functionals with respect to the 2-RDM (two particle reduced density matrix) rather than the wave function, we closely connect the one-density, the 1-RDM (one-particle reduced density matrix), and the geminal functional theories with 2-RDM minimization of the energy. Constrained search with the 2-RDM emphasizes that all energy functional methods must implicitly account for the N-representability of the 2-RDM within their universal functionals. An approximate universal functional for GFT, equivalent to a variational ansatz using the antisymmetrized geminal power wave function, yields energies that are significantly better than those from Hartree–Fock and yet rigorously above the exact energy.
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.
Abstract We show the transformation from a one-particle basis to a geminal basis, transformations between different geminal bases demonstrate the Lie algebra of a geminal basis. From the basis transformations, we express both the wave function and Hamiltonian in the geminal basis. The necessary and sufficient conditions of the exact wave function expanded in a geminal basis are shown to be a Brillouin theorem of geminals. The variational optimization of the geminals in the antisymmetrized geminal power (AGP), antisymmetrized product of geminals (APG) and the full geminal product (FGP) wave function ansätze are discussed. We show that using a geminal replacement operator to describe geminal rotations introduce both primary and secondary rotations. The secondary rotations rotate two geminals in the reference at the same time due to the composite boson nature of geminals. Due to the completeness of the FGP, where all possible geminal combinations are present, the FGP is exact. The number of parameters in the FGP scale exponentially with the number of particles, like the full configuration interaction (FCI). Truncation in the FGP expansion can give compact representations of the wave function since the reference function in the FGP can be either the AGP or APG wave function.
Progress toward the solution of the strongly correlated electron problem has been stymied by the exponential complexity of the wave function. Previous work established an exact two-body exponential product expansion for the ground-state wave function. By developing a reduced density-matrix analog of Dalgarno-Lewis perturbation theory, we prove here that (i) the two-body exponential product expansion is rapidly and globally convergent with each operator representing an order of a renormalized perturbation theory, (ii) the energy of the expansion converges quadratically near the solution, and (iii) the expansion is exact for both ground and excited states. The two-body expansion offers a reduced parametrization of the many-particle wave function as well as the two-particle reduced density matrix with potential applications on both conventional and quantum computers for the study of strongly correlated quantum systems. We demonstrate the result with the exact solution of the contracted Schr\"odinger equation for the molecular chains H4 and H5.