Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Mathematical Physics

Structure of Fermion Density Matrices. II. Antisymmetrized Geminal Powers

Alan Coleman

Published Sep 1, 1965• 269 citations
Publisher

Abstract

Functions which can be expressed as an antisymmetrized power of a single two-particle function (a geminal) occur in the BCS ansatz. They constitute a comparatively tractable generalization of a single Slater determinant. The second-order reduced density matrix and bounds on its eigenvalues are obtained for the most general such antisymmetrized geminal power (AGP). These eigenvalues are physically of great importance, being analogs, in an analysis of the wavefunction into pair states, of the one-particle occupation numbers of the independent-particle model. A necessary and sufficient condition is given that an alleged first-order reduced density matrix be derivable from an N-particle AGP function. It is shown that only for exceptional functions (called extreme) is the generating geminal of the AGP function an eigenfunction of the 2-matrix. The extreme 2-matrix is diagonalized explicitly and provides a vivid example of the fact that one-particle occupation numbers alone are unable to convey certain information of decisive significance about the wavefunction. The extreme 2-matrix satisfies the assumption basic to current theories of superconductivity.

DOI
10.1063/1.1704794

Related papers

A solvable model with an extreme AGP ground state: Relationships among fermion pairs, pairons, and natural spin geminals

Everett G. Larson

A model many-fermion Hamiltonian is presented for which the ground state is asymptotically an Antisymmetrized Geminal Powers (AGP) wave function with largest possible greatest eigenvalue for its two-particle reduced density matrix. Closed analytical expressions and plane-wave expansions are presented for the generating geminal of the AGP ground state and for its one-particle reduced density matrix. The natural orbitals for this generating geminal are plane waves. The generating geminal shows intensely local character in its intracule and corresponds to the formation of a quasi-boson from two fermions. One may appropriately modify this generating geminal to introduce zero occupation numbers of its one-particle reduced density matrix and to make all the nonzero occupation numbers of its one-particle reduced density matrix equal, thus making this geminal a generator of an extreme AGP wave function, with an extreme large eigenvalue for its two-particle reduced density matrix. Closed analytical expressions are also given for this modified geminal and for its one-particle reduced density matrix. The modified generating geminal develops anomolously long “tails” in the intracule function. Nevertheless, both generating geminals introduced here reduce to a renormalized Dirac δ-function of the intracule coordinate, in the appropriate limit. In this limit, both of these generating geminals produce an absolute extreme AGP wave function with the absolute extreme value, N/2, for the Löwdin-normalized occupation number of the generating geminal, and vanishing occupation numbers for all other natural geminals, in the two-particle reduced density matrix of this wavefunction. A demonstration is given showing the formation of quasi-bosons in this ground state and showing their relationship to the generating geminal and to the other natural geminals of the two-particle reduced density matrix of the extreme AGP wave function. The similarities and differences of the features of this model and the accepted models of the superconducting ground state of electrons in metals, and the superfluid ground state of liquid He4 are mentioned.

International Journal of Quantum Chemistry•Mar 10, 1986•1 citation

The AGP model for fermion systems

Alan Coleman

A mathematically rigorous treatment of the antisymmetrized geminal power (AGP) model is given. The model is advocated as an extremely flexible tool for studying a wide variety of fermion systems, characterized by a single Slater determinant, on the one hand, or Yang's most highly correlated wave function, on the other. When the first-order reduced density matrix of the actual system has eigenvalues which are evenly degenerate, the model can treat one-particle properties exactly. A formula is obtained for the precise energy of an AGP state in the Born-Oppenheimer approximation. An estimate is given of the error which occurs when the second-order density matrix of an AGP state is approximated by the widely used expression involving the so-called anomalous Green's function. © 1997 John Wiley & Sons, Inc.

International Journal of Quantum Chemistry•Jan 1, 1997•42 citations

Pseudo-BCS wave function from density matrix decomposition: Application in auxiliary-field quantum Monte Carlo

Zhi-Yu Xiao, Hao Shi, Shiwei Zhang

We present a method to construct pseudo-BCS wave functions from the one-body density matrix.The resulting many-body wave function, which can be produced for any fermion systems, including those with purely repulsive interactions, has the form of a number-projected BCS form, or antisymmetrized geminal power (AGP).Such wave functions provide a better ansatz for correlated fermion systems than a single Slater determinant, and often better than a linear combination of Slater determinants (for example from a truncated active space calculation).We describe a procedure to build such a wave function conveniently from a given reduced density matrix of the system, rather than from a mean-field solution (which gives a Slater determinant for repulsive interactions).The pseudo-BCS wave function thus obtained reproduces the density matrix or minimizes the difference between the input and resulting density matrices.One application of the pseudo-BCS wave function is in auxiliary-field quantum Monte Carlo (AFQMC) calculations as the trial wave function to control the sign/phase problem.AFQMC is often among the most accurate general methods for correlated fermion systems.We show that the pseudo-BCS form further reduces the constraint bias and leads to improved accuracy compared to the usual Slater determinant trial wave functions, using the two-dimensional Hubbard model as an example.Furthermore, the pseudo-BCS trial wave function allows a new systematically improvable self-consistent approach, with pseudo-BCS trial wave function iteratively generated by AFQMC via the one-body density matrix.

arXiv:2010.06055•Physical Review Research•Jan 21, 2021•9 citations

Excitation operators associated with antisymmetrized geminal-power states

Brian Weiner, Osvaldo Goscinski

Various expressions for antisymmetrized geminal-power (AGP) states are reviewed. With the use of a simple relationship for creation and annihilation operators of BCS (Bardeen, Cooper, and Schrieffer)-type states, such operators are characterized for AGP and GAGP (generalized AGP) schemes. The effect of degeneracy in the first-order reduced density matrix on these operators is displayed.

Physical Review A•Jan 1, 1983•29 citations

Cumulant expansion of the reduced density matrices

Werner Kutzelnigg, Debashis Mukherjee

k -particle cumulants λk (for 2⩽k⩽n) corresponding to the k-particle reduced density matrices γk for an n-fermion system are defined via a generating function. The two-particle cumulant λ2 describes two-particle correlations (excluding exchange), λ3 genuine three-particle correlations etc. The properties of these cumulants are analyzed. Conditions for vanishing of certain λk are formulated. Necessary and sufficient for λ2=0 is the well-known idempotency condition γ2=γ for γ≡γ1. For λ3=0 to hold, a general necessary condition is Tr{2γ3−3γ2+γ}=0, for three special forms of the wave function (arbitrary two-electron state, antisymmetrized product of strongly orthogonal geminals on antisymmetrized geminal power wave function of extreme type) 2γ3−3γ2+γ=0 turns out to be necessary and sufficient. For a multiconfiguration self-consistent field wave function the only nonvanishing matrix elements of the cumulants are those where all labels refer to active (partially occupied) spin orbitals. Spin-free cumulants Λk corresponding to the spin-free reduced density matrices Γk are also defined and analyzed. The main interest in the density cumulants is in connection with the recently formulated normal ordering and the corresponding Wick theorem for arbitrary reference functions, but they are also useful for an analysis of electron correlation.

The Journal of Chemical Physics•Feb 8, 1999•291 citations

Density matrix embedding in an antisymmetrized geminal power bath

Takashi Tsuchimochi, Matthew Welborn, Troy Van Voorhis

Density matrix embedding theory (DMET) has emerged as a powerful tool for performing wave function-in-wave function embedding for strongly correlated systems. In traditional DMET, an accurate calculation is performed on a small impurity embedded in a mean field bath. Here, we extend the original DMET equations to account for correlation in the bath via an antisymmetrized geminal power (AGP) wave function. The resulting formalism has a number of advantages. First, it allows one to properly treat the weak correlation limit of independent pairs, which DMET is unable to do with a mean-field bath. Second, it associates a size extensive correlation energy with a given density matrix (for the models tested), which AGP by itself is incapable of providing. Third, it provides a reasonable description of charge redistribution in strongly correlated but non-periodic systems. Thus, AGP-DMET appears to be a good starting point for describing electron correlation in molecules, which are aperiodic and possess both strong and weak electron correlation.

The Journal of Chemical Physics•Jul 13, 2015•56 citations
Functions which can be expressed as an antisymmetrized power of a single two-particle function (a geminal) occur in the BCS ansatz. They constitute a comparatively tractable generalization of a single Slater determinant. The second-order reduced density matrix and bounds on its eigenvalues are obtained for the most general such antisymmetrized geminal power (AGP). These eigenvalues are physically of great importance, being analogs, in an analysis of the wavefunction into pair states, of the one-particle occupation numbers of the independent-particle model. A necessary and sufficient condition is given that an alleged first-order reduced density matrix be derivable from an N-particle AGP function. It is shown that only for exceptional functions (called extreme) is the generating geminal of the AGP function an eigenfunction of the 2-matrix. The extreme 2-matrix is diagonalized explicitly and provides a vivid example of the fact that one-particle occupation numbers alone are unable to convey certain information of decisive significance about the wavefunction. The extreme 2-matrix satisfies the assumption basic to current theories of superconductivity.
A model many-fermion Hamiltonian is presented for which the ground state is asymptotically an Antisymmetrized Geminal Powers (AGP) wave function with largest possible greatest eigenvalue for its two-particle reduced density matrix. Closed analytical expressions and plane-wave expansions are presented for the generating geminal of the AGP ground state and for its one-particle reduced density matrix. The natural orbitals for this generating geminal are plane waves. The generating geminal shows intensely local character in its intracule and corresponds to the formation of a quasi-boson from two fermions. One may appropriately modify this generating geminal to introduce zero occupation numbers of its one-particle reduced density matrix and to make all the nonzero occupation numbers of its one-particle reduced density matrix equal, thus making this geminal a generator of an extreme AGP wave function, with an extreme large eigenvalue for its two-particle reduced density matrix. Closed analytical expressions are also given for this modified geminal and for its one-particle reduced density matrix. The modified generating geminal develops anomolously long “tails” in the intracule function. Nevertheless, both generating geminals introduced here reduce to a renormalized Dirac δ-function of the intracule coordinate, in the appropriate limit. In this limit, both of these generating geminals produce an absolute extreme AGP wave function with the absolute extreme value, N/2, for the Löwdin-normalized occupation number of the generating geminal, and vanishing occupation numbers for all other natural geminals, in the two-particle reduced density matrix of this wavefunction. A demonstration is given showing the formation of quasi-bosons in this ground state and showing their relationship to the generating geminal and to the other natural geminals of the two-particle reduced density matrix of the extreme AGP wave function. The similarities and differences of the features of this model and the accepted models of the superconducting ground state of electrons in metals, and the superfluid ground state of liquid He4 are mentioned.
A mathematically rigorous treatment of the antisymmetrized geminal power (AGP) model is given. The model is advocated as an extremely flexible tool for studying a wide variety of fermion systems, characterized by a single Slater determinant, on the one hand, or Yang's most highly correlated wave function, on the other. When the first-order reduced density matrix of the actual system has eigenvalues which are evenly degenerate, the model can treat one-particle properties exactly. A formula is obtained for the precise energy of an AGP state in the Born-Oppenheimer approximation. An estimate is given of the error which occurs when the second-order density matrix of an AGP state is approximated by the widely used expression involving the so-called anomalous Green's function. © 1997 John Wiley & Sons, Inc.
We present a method to construct pseudo-BCS wave functions from the one-body density matrix.The resulting many-body wave function, which can be produced for any fermion systems, including those with purely repulsive interactions, has the form of a number-projected BCS form, or antisymmetrized geminal power (AGP).Such wave functions provide a better ansatz for correlated fermion systems than a single Slater determinant, and often better than a linear combination of Slater determinants (for example from a truncated active space calculation).We describe a procedure to build such a wave function conveniently from a given reduced density matrix of the system, rather than from a mean-field solution (which gives a Slater determinant for repulsive interactions).The pseudo-BCS wave function thus obtained reproduces the density matrix or minimizes the difference between the input and resulting density matrices.One application of the pseudo-BCS wave function is in auxiliary-field quantum Monte Carlo (AFQMC) calculations as the trial wave function to control the sign/phase problem.AFQMC is often among the most accurate general methods for correlated fermion systems.We show that the pseudo-BCS form further reduces the constraint bias and leads to improved accuracy compared to the usual Slater determinant trial wave functions, using the two-dimensional Hubbard model as an example.Furthermore, the pseudo-BCS trial wave function allows a new systematically improvable self-consistent approach, with pseudo-BCS trial wave function iteratively generated by AFQMC via the one-body density matrix.
k -particle cumulants λk (for 2⩽k⩽n) corresponding to the k-particle reduced density matrices γk for an n-fermion system are defined via a generating function. The two-particle cumulant λ2 describes two-particle correlations (excluding exchange), λ3 genuine three-particle correlations etc. The properties of these cumulants are analyzed. Conditions for vanishing of certain λk are formulated. Necessary and sufficient for λ2=0 is the well-known idempotency condition γ2=γ for γ≡γ1. For λ3=0 to hold, a general necessary condition is Tr{2γ3−3γ2+γ}=0, for three special forms of the wave function (arbitrary two-electron state, antisymmetrized product of strongly orthogonal geminals on antisymmetrized geminal power wave function of extreme type) 2γ3−3γ2+γ=0 turns out to be necessary and sufficient. For a multiconfiguration self-consistent field wave function the only nonvanishing matrix elements of the cumulants are those where all labels refer to active (partially occupied) spin orbitals. Spin-free cumulants Λk corresponding to the spin-free reduced density matrices Γk are also defined and analyzed. The main interest in the density cumulants is in connection with the recently formulated normal ordering and the corresponding Wick theorem for arbitrary reference functions, but they are also useful for an analysis of electron correlation.
Density matrix embedding theory (DMET) has emerged as a powerful tool for performing wave function-in-wave function embedding for strongly correlated systems. In traditional DMET, an accurate calculation is performed on a small impurity embedded in a mean field bath. Here, we extend the original DMET equations to account for correlation in the bath via an antisymmetrized geminal power (AGP) wave function. The resulting formalism has a number of advantages. First, it allows one to properly treat the weak correlation limit of independent pairs, which DMET is unable to do with a mean-field bath. Second, it associates a size extensive correlation energy with a given density matrix (for the models tested), which AGP by itself is incapable of providing. Third, it provides a reasonable description of charge redistribution in strongly correlated but non-periodic systems. Thus, AGP-DMET appears to be a good starting point for describing electron correlation in molecules, which are aperiodic and possess both strong and weak electron correlation.