, which controls how closely the truncated expansion approximates the full AGP-CI state. Benchmark applications to the Hubbard model and to the small molecules H

_2

O and N

_2

demonstrate that the proposed wavefunctions achieve consistently high accuracy and outperform the LC-AGP, particularly for systems with more electrons and in strongly correlated regimes.

DOI: 10.48550/arxiv.2604.14115

Related papers

Configuration interaction extension of AGP for incorporating inter-geminal correlations

Airi Kawasaki, Fei Gao, Gustavo E. Scuseria

In this paper, we develop a class of antisymmetrized geminal power configuration interaction (AGP-CI) wave functions that extend the AGP framework by incorporating inter-geminal correlations through a CI expansion. To make these wavefunctions computationally tractable, we evaluate them by rewriting the AGP-CI ansatz as a linear combination of AGPs (LC-AGP), for which overlaps and Hamiltonian matrix elements can be computed with standard AGP machinery. Motivated by border-rank decompositions, we further reorganize this ansatz into a compact linear combination of AGPs depending on a small deformation parameter τ, which controls how closely the truncated expansion approximates the full AGP-CI state. Benchmark applications to the Hubbard model and to the small molecules H2O and N2 demonstrate that the proposed wavefunctions achieve consistently high accuracy and outperform the LC-AGP, particularly for systems with more electrons and in strongly correlated regimes.In this paper, we develop a class of antisymmetrized geminal power configuration interaction (AGP-CI) wave functions that extend the AGP framework by incorporating inter-geminal correlations through a CI expansion. To make these wavefunctions computationally tractable, we evaluate them by rewriting the AGP-CI ansatz as a linear combination of AGPs (LC-AGP), for which overlaps and Hamiltonian matrix elements can be computed with standard AGP machinery. Motivated by border-rank decompositions, we further reorganize this ansatz into a compact linear combination of AGPs depending on a small deformation parameter τ, which controls how closely the truncated expansion approximates the full AGP-CI state. Benchmark applications to the Hubbard model and to the small molecules H2O and N2 demonstrate that the proposed wavefunctions achieve consistently high accuracy and outperform the LC-AGP, particularly for systems with more electrons and in strongly correlated regimes.

The Journal of Chemical PhysicsJun 17, 20260 citations

Configuration interaction with antisymmetrized geminal powers

Wataru Uemura, Shusuke Kasamatsu, Osamu Sugino

To avoid the combinatorial computational cost of configuration interaction (CI), we previously introduced the symmetric tensor decomposition CI (STD-CI) method, which takes advantage of the antisymmetric nature of the electronic wave function and expresses the CI coefficients compactly as a series of Kronecker product states (STD series) [W. Uemura and O. Sugino, Phys. Rev. Lett. 109, 253001 (2012)]. Here we extend the variational degrees of freedom by using different molecular orbitals for different terms in the STD series. This scheme is equivalent to the linear combination of the Hartree-Fock-Bogoliubov state or the antisymmetrized geminal powers (AGPs). The total energy converges very rapidly within 0.72 $\ensuremath{\mu}\mathrm{hartree}$ taking only 10 terms for the water molecule, and the convergence is likewise fast for Hubbard tetramers. The computational cost scales as the fifth power of the number of electrons and the square of the number of terms in the STD series, indicating the promise of this AGP-based scheme for highly accurate and efficient computation of quantum systems.To avoid the combinatorial computational cost of configuration interaction (CI), we previously introduced the symmetric tensor decomposition CI (STD-CI) method, which takes advantage of the antisymmetric nature of the electronic wave function and expresses the CI coefficients compactly as a series of Kronecker product states (STD series) [W. Uemura and O. Sugino, Phys. Rev. Lett. 109, 253001 (2012)]. Here we extend the variational degrees of freedom by using different molecular orbitals for different terms in the STD series. This scheme is equivalent to the linear combination of the Hartree-Fock-Bogoliubov state or the antisymmetrized geminal powers (AGPs). The total energy converges very rapidly within 0.72 $\ensuremath{\mu}\mathrm{hartree}$ taking only 10 terms for the water molecule, and the convergence is likewise fast for Hubbard tetramers. The computational cost scales as the fifth power of the number of electrons and the square of the number of terms in the STD series, indicating the promise of this AGP-based scheme for highly accurate and efficient computation of quantum systems.

arXiv:1504.06250Physical Review AJun 16, 201513 citations

Geminal-based configuration interaction

Thomas M. Henderson, Gustavo E. Scuseria

The antisymmetrized geminal power (AGP) wave function has a long history and considerable conceptual appeal, but in many situations its accuracy is wanting. Here, we consider a form of configuration interaction (CI) based upon the AGP wave function and take advantage of its killing operators to construct an excitation manifold. Our geminal CI reduces to standard single-determinant-based CI in the limit in which AGP reduces to a single determinant. It substantially improves upon AGP in the reduced BCS Hamiltonian, which serves as a prototype for the kinds of strong pairing correlations relevant in Bardeen-Cooper-Schrieffer-style superconductivity. Moreover, our geminal CI naturally generalizes to add correlation to more general geminal-based wave functions than AGP.The antisymmetrized geminal power (AGP) wave function has a long history and considerable conceptual appeal, but in many situations its accuracy is wanting. Here, we consider a form of configuration interaction (CI) based upon the AGP wave function and take advantage of its killing operators to construct an excitation manifold. Our geminal CI reduces to standard single-determinant-based CI in the limit in which AGP reduces to a single determinant. It substantially improves upon AGP in the reduced BCS Hamiltonian, which serves as a prototype for the kinds of strong pairing correlations relevant in Bardeen-Cooper-Schrieffer-style superconductivity. Moreover, our geminal CI naturally generalizes to add correlation to more general geminal-based wave functions than AGP.

arXiv:1906.11361The Journal of Chemical PhysicsAug 2, 201947 citations

Geminal Replacement Models Based on AGP

Rishab Dutta, Thomas M. Henderson, Gustavo E. Scuseria

The antisymmetrized geminal power (AGP) wave function has a long history and is known by different names in various chemical and physical problems. There has been recent interest in using AGP as a starting point for strongly correlated electrons. Here, we show that in a seniority-conserving regime, different AGP-based correlator representations based on generators of the algebra, killing operators, and geminal replacement operators are all equivalent. We implement one representation that uses number operators as correlators and has linearly independent curvilinear metrics to distinguish the regions of Hilbert space. This correlation method called J -CI provides excellent accuracy in energies when applied to the pairing Hamiltonian.The antisymmetrized geminal power (AGP) wave function has a long history and is known by different names in various chemical and physical problems. There has been recent interest in using AGP as a starting point for strongly correlated electrons. Here, we show that in a seniority-conserving regime, different AGP-based correlator representations based on generators of the algebra, killing operators, and geminal replacement operators are all equivalent. We implement one representation that uses number operators as correlators and has linearly independent curvilinear metrics to distinguish the regions of Hilbert space. This correlation method called J -CI provides excellent accuracy in energies when applied to the pairing Hamiltonian.

arXiv:2008.00552Journal of Chemical Theory and ComputationSep 1, 202035 citations

Low-rank decomposition on the antisymmetric product of geminals for strongly correlated electrons

Airi Kawasaki, Naoki Nakatani

We investigated some variational methods to compute a wavefunction based on antisymmetric product of geminals (APG). The Waring decomposition on the APG wavefunction leads a finite sum of antisymmetrized geminal power (AGP) wavefunctions, each for which the variational principle can be applied. We call this as AGP-CI method which provides a variational solution of the APG wavefunction efficiently. However, number of AGP wavefunctions in the exact AGP-CI formalism become exponentially large in case of many-electron systems. Therefore, we also investigate the low-rank APG wavefunction, in which the coefficient matrices of geminals are factorized by the Schur decomposition. Interestingly, only a few non-zero eigenvalues (up to half number of electrons) were found from the Schur decomposition on the APG wavefunction. We developed some methods to approximate the APG wavefunction by lowering the ranks of coefficient matrices of geminals, and demonstrate their performance. Our new geminal method based on the low-rank decomposition can drastically reduce the number of variational parameters, although no efficient formalism has been elaborated so far, due to the difficulty of optimization caused by the inability to compute analytical derivatives.We investigated some variational methods to compute a wavefunction based on antisymmetric product of geminals (APG). The Waring decomposition on the APG wavefunction leads a finite sum of antisymmetrized geminal power (AGP) wavefunctions, each for which the variational principle can be applied. We call this as AGP-CI method which provides a variational solution of the APG wavefunction efficiently. However, number of AGP wavefunctions in the exact AGP-CI formalism become exponentially large in case of many-electron systems. Therefore, we also investigate the low-rank APG wavefunction, in which the coefficient matrices of geminals are factorized by the Schur decomposition. Interestingly, only a few non-zero eigenvalues (up to half number of electrons) were found from the Schur decomposition on the APG wavefunction. We developed some methods to approximate the APG wavefunction by lowering the ranks of coefficient matrices of geminals, and demonstrate their performance. Our new geminal method based on the low-rank decomposition can drastically reduce the number of variational parameters, although no efficient formalism has been elaborated so far, due to the difficulty of optimization caused by the inability to compute analytical derivatives.

arXiv:2410.06666Molecular PhysicsJan 21, 20255 citations

Four-body correlation embedded in antisymmetrized geminal power wave function

Airi Kawasaki, Osamu Sugino

We extend the Coleman's antisymmetrized geminal power (AGP) to develop a wave function theory that can incorporate up to four-body correlation in a region of strong correlation. To facilitate the variational determination of the wave function, the total energy is rewritten in terms of the traces of geminals. This novel trace formula is applied to a simple model system consisting of one dimensional Hubbard ring with a site of strong correlation. Our scheme significantly improves the result obtained by the AGP-configuration interaction scheme of Uemura et al. and also achieves more efficient compression of the degrees of freedom of the wave function. We regard the result as a step toward a first-principles wave function theory for a strongly correlated point defect or adsorbate embedded in an AGP-based mean-field medium.We extend the Coleman's antisymmetrized geminal power (AGP) to develop a wave function theory that can incorporate up to four-body correlation in a region of strong correlation. To facilitate the variational determination of the wave function, the total energy is rewritten in terms of the traces of geminals. This novel trace formula is applied to a simple model system consisting of one dimensional Hubbard ring with a site of strong correlation. Our scheme significantly improves the result obtained by the AGP-configuration interaction scheme of Uemura et al. and also achieves more efficient compression of the degrees of freedom of the wave function. We regard the result as a step toward a first-principles wave function theory for a strongly correlated point defect or adsorbate embedded in an AGP-based mean-field medium.

arXiv:1609.01438The Journal of Chemical PhysicsDec 28, 20169 citations

τ

, which controls how closely the truncated expansion approximates the full AGP-CI state. Benchmark applications to the Hubbard model and to the small molecules H

_2

O and N

_2

demonstrate that the proposed wavefunctions achieve consistently high accuracy and outperform the LC-AGP, particularly for systems with more electrons and in strongly correlated regimes.