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

arXiv:2307.13690

Single reference treatment of strongly correlated H$_4$ and H$_{10}$ isomers with Richardson-Gaudin states

Paul A. Johnson, A. Eugene DePrince

Published Jul 25, 2023• 0 citations
arXivHTMLPDF

Abstract

Richardson-Gaudin (RG) states are employed as a variational wavefunction ansatz for strongly correlated isomers of H$_4$ and H$_{10}$. In each case a single RG state describes the seniority-zero sector quite well. Simple natural orbital functionals offer a cheap and reasonable approximation of the outstanding weak correlation in the seniority-zero sector, while systematic improvement is achieved by performing a configuration interaction (CI) in terms of RG states. Other pair theories (e.g. generalized valence bond and pair-coupled-cluster doubles) can provide a good description of many of the geometries considered, but, at short distances, the wavefunctions for the 2D and 3D structures of H$_{10}$ take the form of an RG state that cannot be described by these other theories.

DOI
10.48550/arxiv.2307.13690

Related papers

Single Reference Treatment of Strongly Correlated H 4 and H 10 Isomers with Richardson–Gaudin States

Paul A. Johnson, A. Eugene DePrince

Richardson–Gaudin (RG) states are employed as a variational wave function ansatz for strongly correlated isomers of H 4 and H 10 . In each case, a single RG state describes the seniority-zero sector quite well. Simple natural orbital functionals offer a cheap and reasonable approximation of the outstanding weak correlation in the seniority-zero sector, while systematic improvement is achieved by performing a configuration interaction in terms of RG states.

Journal of Chemical Theory and Computation•Nov 13, 2023•18 citations

Richardson–Gaudin states of non-zero seniority. II. Single-reference treatment of strong correlation

Paul A. Johnson

Strongly correlated systems are well described as a configuration interaction of Slater determinants classified by their number of unpaired electrons. This treatment is, however, unfeasible. In this paper, it is demonstrated that single reference methods built from Richardson-Gaudin states yield comparable results at polynomial cost.

arXiv:2506.09379•The Journal of Chemical Physics•Aug 22, 2025•4 citations

Seniority eigenstate configuration interaction

Thomas M. Henderson, Guo P. Chen, Gustavo E. Scuseria

Zero-seniority methods have shown great promise for the description of strongly correlated electronic systems. Other seniority sectors have been much less explored, and in particular, the maximal seniority sector and zero seniority have the same underlying algebraic structure. We introduce a seniority eigenstate configuration interaction in which the wave function is constrained to have good fixed local seniority for each paired orbital, by which we mean we partition orbitals into a pairing set with seniority zero, and a spin set with seniority one. We show how to build the effective Hamiltonian for this ansatz, and demonstrate that high-seniority wave functions have unexpectedly excellent accuracy for strongly correlated fermionic systems, with accuracy competitive with or better than seniority zero for the Hubbard model and for the dissociation of the nitrogen molecule.

The Journal of Chemical Physics•Jul 7, 2026•0 citations

Seniority Eigenstate Configuration Interaction

Thomas M. Henderson, Guo P. Chen, Gustavo E. Scuseria

Zero-seniority methods have shown great promise for the description of strongly-correlated electronic systems. Other seniority sectors have been much less explored, and in particular the maximal seniority sector and zero seniority have the same underlying algebraic structure. We introduce a seniority eigenstate configuration interaction in which the wave function is constrained to have good fixed local seniority for each paired orbital, by which we mean we partition orbitals into a pairing set with seniority zero, and a spin set with seniority one. We show how to build the effective Hamiltonian for this ansatz, and demonstrate that high-seniority wave functions have unexpectedly excellent accuracy for strongly-correlated fermionic systems, with accuracy competitive with or better than seniority zero for the Hubbard model and for the dissociation of the nitrogen molecule.

arXiv:2604.19063•Apr 21, 2026•0 citations

Model Wavefunction Forms to Describe Strong Correlation in Quantum Chemistry

Paul A. Johnson

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).

•Jan 6, 2015•3 citations

Richardson-Gaudin Configuration-Interaction for nuclear pairing correlations

Stijn De Baerdemacker, Pieter W. Claeys, Jean-Sébastien Caux +2 more

Background: The nuclear many-body system is a strongly correlated quantum system, posing serious challenges for perturbative approaches starting from uncorrelated reference states. The last decade has witnessed considerable progress in the accurate treatment of pairing correlations, one of the major components in medium-sized nuclei, reaching accuracies below the 1% level of the correlation energy. Purpose: Development of a quantum many-body method for pairing correlations that is (a) competitive in the 1% error range, and (b) can be systematically improved with a fast (exponential) convergence rate. Method: The present paper capitalizes upon ideas from Richardson-Gaudin integrability. The proposed method is a two-step approach. The first step consists of the optimization of a Richardson-Gaudin ground state as variational trial state. At the second step, the complete set of excited states on top of this Richardson-Gaudin ground state is used as an optimal basis for a Configuration Interaction method in an increasingly large effective Hilbert space. Results: The performance of the variational Richardson-Gaudin (varRG) and Richardson-Gaudin Configuration Interaction (RGCI) method is benchmarked against exact results using an effective $G$-matrix interaction for the Sn region. The varRG already reaches accuracies around the 1% level of the correlation energies, and the RGCI step sees an additional improvement scaling exponentially with the size of the effective Hilbert space. Conclusions: The Richardson-Gaudin models of integrability provide an optimized complete basis set for pairing correlations.

arXiv:1712.01673•Dec 5, 2017•2 citations
Single reference treatment of strongly correlated H4_44​ and H10_{10}10​ isomers with Richardson-Gaudin states
Richardson-Gaudin (RG) states are employed as a variational wavefunction ansatz for strongly correlated isomers of H4_44​ and H10_{10}10​. In each case a single RG state describes the seniority-zero sector quite well. Simple natural orbital functionals offer a cheap and reasonable approximation of the outstanding weak correlation in the seniority-zero sector, while systematic improvement is achieved by performing a configuration interaction (CI) in terms of RG states. Other pair theories (e.g. generalized valence bond and pair-coupled-cluster doubles) can provide a good description of many of the geometries considered, but, at short distances, the wavefunctions for the 2D and 3D structures of H10_{10}10​ take the form of an RG state that cannot be described by these other theories.
Zero-seniority methods have shown great promise for the description of strongly correlated electronic systems. Other seniority sectors have been much less explored, and in particular, the maximal seniority sector and zero seniority have the same underlying algebraic structure. We introduce a seniority eigenstate configuration interaction in which the wave function is constrained to have good fixed local seniority for each paired orbital, by which we mean we partition orbitals into a pairing set with seniority zero, and a spin set with seniority one. We show how to build the effective Hamiltonian for this ansatz, and demonstrate that high-seniority wave functions have unexpectedly excellent accuracy for strongly correlated fermionic systems, with accuracy competitive with or better than seniority zero for the Hubbard model and for the dissociation of the nitrogen molecule.
Zero-seniority methods have shown great promise for the description of strongly-correlated electronic systems. Other seniority sectors have been much less explored, and in particular the maximal seniority sector and zero seniority have the same underlying algebraic structure. We introduce a seniority eigenstate configuration interaction in which the wave function is constrained to have good fixed local seniority for each paired orbital, by which we mean we partition orbitals into a pairing set with seniority zero, and a spin set with seniority one. We show how to build the effective Hamiltonian for this ansatz, and demonstrate that high-seniority wave functions have unexpectedly excellent accuracy for strongly-correlated fermionic systems, with accuracy competitive with or better than seniority zero for the Hubbard model and for the dissociation of the nitrogen molecule.
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).
Background: The nuclear many-body system is a strongly correlated quantum system, posing serious challenges for perturbative approaches starting from uncorrelated reference states. The last decade has witnessed considerable progress in the accurate treatment of pairing correlations, one of the major components in medium-sized nuclei, reaching accuracies below the 1% level of the correlation energy. Purpose: Development of a quantum many-body method for pairing correlations that is (a) competitive in the 1% error range, and (b) can be systematically improved with a fast (exponential) convergence rate. Method: The present paper capitalizes upon ideas from Richardson-Gaudin integrability. The proposed method is a two-step approach. The first step consists of the optimization of a Richardson-Gaudin ground state as variational trial state. At the second step, the complete set of excited states on top of this Richardson-Gaudin ground state is used as an optimal basis for a Configuration Interaction method in an increasingly large effective Hilbert space. Results: The performance of the variational Richardson-Gaudin (varRG) and Richardson-Gaudin Configuration Interaction (RGCI) method is benchmarked against exact results using an effective GGG-matrix interaction for the Sn region. The varRG already reaches accuracies around the 1% level of the correlation energies, and the RGCI step sees an additional improvement scaling exponentially with the size of the effective Hilbert space. Conclusions: The Richardson-Gaudin models of integrability provide an optimized complete basis set for pairing correlations.