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 Chemical Theory and Computation

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

Paul A. Johnson, A. Eugene DePrince

Published Nov 13, 2023• 18 citations
Publisher

Abstract

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.

DOI
10.1021/acs.jctc.3c00807

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

arXiv:2307.13690•Jul 25, 2023•0 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

Reduced density matrices/static correlation functions of Richardson–Gaudin states without rapidities

Alexandre Faribault, Claude Dimo, Jean-David Moisset +1 more

Seniority-zero geminal wavefunctions are known to capture bond-breaking correlation. Among this class of wavefunctions, Richardson-Gaudin states stand out as they are eigenvectors of a model Hamiltonian. This provides a clear physical picture, clean expressions for reduced density matrix (RDM) elements, and systematic improvement (with a complete set of eigenvectors). Known expressions for the RDM elements require the computation of rapidities, which are obtained by first solving for the so-called eigenvalue based variables (EBV) and then root-finding a Lagrange interpolation polynomial. In this paper, we obtain expressions for the RDM elements directly in terms of the EBV. The final expressions can be computed at the same cost as the rapidity expressions. Therefore, except, in particular, circumstances, it is entirely unnecessary to compute rapidities at all. The RDM elements require numerically inverting a matrix, and while this is usually undesirable, we demonstrate that it is stable, except when there is degeneracy in the single-particle energies. In such cases, a different construction would be required.

arXiv:2209.00623•The Journal of Chemical Physics•Dec 1, 2022•22 citations

Beyond a Richardson–Gaudin Mean-Field: Slater–Condon Rules and Perturbation Theory

Paul A. Johnson

Richardson-Gaudin states provide a basis of the Hilbert space for strongly correlated electrons. In this study, optimal expressions for the transition density matrix elements between Richardson-Gaudin states are obtained with a cost comparable with the corresponding reduced density matrix elements. Analogues of the Slater-Condon rules are identified based on the number of near-zero singular values of the RG state overlap matrix. Finally, a perturbative approach is shown to be close in quality to a configuration interaction of Richardson-Gaudin states while being feasible to compute.

arXiv:2405.00279•The Journal of Physical Chemistry A•Jul 15, 2024•5 citations
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.
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.
Seniority-zero geminal wavefunctions are known to capture bond-breaking correlation. Among this class of wavefunctions, Richardson-Gaudin states stand out as they are eigenvectors of a model Hamiltonian. This provides a clear physical picture, clean expressions for reduced density matrix (RDM) elements, and systematic improvement (with a complete set of eigenvectors). Known expressions for the RDM elements require the computation of rapidities, which are obtained by first solving for the so-called eigenvalue based variables (EBV) and then root-finding a Lagrange interpolation polynomial. In this paper, we obtain expressions for the RDM elements directly in terms of the EBV. The final expressions can be computed at the same cost as the rapidity expressions. Therefore, except, in particular, circumstances, it is entirely unnecessary to compute rapidities at all. The RDM elements require numerically inverting a matrix, and while this is usually undesirable, we demonstrate that it is stable, except when there is degeneracy in the single-particle energies. In such cases, a different construction would be required.
Richardson-Gaudin states provide a basis of the Hilbert space for strongly correlated electrons. In this study, optimal expressions for the transition density matrix elements between Richardson-Gaudin states are obtained with a cost comparable with the corresponding reduced density matrix elements. Analogues of the Slater-Condon rules are identified based on the number of near-zero singular values of the RG state overlap matrix. Finally, a perturbative approach is shown to be close in quality to a configuration interaction of Richardson-Gaudin states while being feasible to compute.