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:2405.00279 · The Journal of Physical Chemistry A

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

Paul A. Johnson

Published Jul 15, 2024• 5 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1021/acs.jpca.4c02857

Related papers

Transition density matrices of Richardson–Gaudin states

Paul A. Johnson, Hubert Fortin, Samuel Cloutier +1 more

Recently, ground state eigenvectors of the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian, Richardson-Gaudin (RG) states, have been employed as a wavefunction ansatz for strong correlation. This wavefunction physically represents a mean-field of pairs of electrons (geminals) with a constant pairing strength. To move beyond the mean-field, one must develop the wavefunction on the basis of all the RG states. This requires both practical expressions for transition density matrices and an idea of which states are most important in the expansion. In this contribution, we present expressions for the transition density matrix elements and calculate them numerically for half-filled picket-fence models (reduced BCS models with constant energy spacing). There are no Slater-Condon rules for RG states, though an analog of the aufbau principle proves to be useful in choosing which states are important.

arXiv:2012.10477•The Journal of Chemical Physics•Mar 28, 2021•14 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

Reduced density matrices of Richardson–Gaudin states in the Gaudin algebra basis

Charles‐Émile Fecteau, Hubert Fortin, Samuel Cloutier +1 more

Eigenvectors of the reduced Bardeen–Cooper–Schrieffer Hamiltonian have recently been employed as a variational wavefunction ansatz in quantum chemistry. This wavefunction is a mean-field of pairs of electrons (geminals). In this contribution, we report optimal expressions for their reduced density matrices in both the original physical basis and the basis of the Richardson–Gaudin pairs. Physical basis expressions were originally reported by Gorohovsky and Bettelheim [Phys. Rev. B 84, 224503 (2011)]. In each case, the expressions scale like O(N4), with the most expensive step being the solution of linear equations. Analytic gradients are also reported in the physical basis. These expressions are an important step toward practical mean-field methods to treat strongly correlated electrons.

arXiv:2008.12713•The Journal of Chemical Physics•Oct 28, 2020•34 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

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

Richardson-Gaudin integrable systems (and beyond) for strongly correlated quantum many-body systems

Stijn De Baerdemacker, Paul W. Ayers, Patrick Bultinck +4 more

Integrable systems take a special niche in quantum many-body systems. Due to the underlying algebraic structure, it is possible to capture all the quantum correlations exactly within a Bethe Ansatz wavefunction at a favourably scaling computational cost. Richardson [1] and Gaudin [2] have shown that the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian of superconductivity [3] is integrable and that the Bethe Ansatz state corresponds to a product wave function of generalized Cooper pairs [4]. In this formulation, the Bethe ansatz state is equivalent to a non-orthogonal geminal, which is widely used in quantum chemistry to describe electron pair correlations in atoms and molecules [5]. Inspired by Richardson-Gaudin integrable systems, it is possible to construct classes of non-orthogonal geminals, allowing for an accurate description of realistic and non-integrable strongly correlated atoms and molecules, while retaining the favorable computational scaling [6, 7]. In this contribution, I will discuss Richardson-Gaudin integrability and why it is useful for strongly correlated systems.

•Jan 1, 2013•0 citations
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.
Recently, ground state eigenvectors of the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian, Richardson-Gaudin (RG) states, have been employed as a wavefunction ansatz for strong correlation. This wavefunction physically represents a mean-field of pairs of electrons (geminals) with a constant pairing strength. To move beyond the mean-field, one must develop the wavefunction on the basis of all the RG states. This requires both practical expressions for transition density matrices and an idea of which states are most important in the expansion. In this contribution, we present expressions for the transition density matrix elements and calculate them numerically for half-filled picket-fence models (reduced BCS models with constant energy spacing). There are no Slater-Condon rules for RG states, though an analog of the aufbau principle proves to be useful in choosing which states are important.
Eigenvectors of the reduced Bardeen–Cooper–Schrieffer Hamiltonian have recently been employed as a variational wavefunction ansatz in quantum chemistry. This wavefunction is a mean-field of pairs of electrons (geminals). In this contribution, we report optimal expressions for their reduced density matrices in both the original physical basis and the basis of the Richardson–Gaudin pairs. Physical basis expressions were originally reported by Gorohovsky and Bettelheim [Phys. Rev. B 84, 224503 (2011)]. In each case, the expressions scale like O(N4), with the most expensive step being the solution of linear equations. Analytic gradients are also reported in the physical basis. These expressions are an important step toward practical mean-field methods to treat strongly correlated electrons.
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.
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.
Integrable systems take a special niche in quantum many-body systems. Due to the underlying algebraic structure, it is possible to capture all the quantum correlations exactly within a Bethe Ansatz wavefunction at a favourably scaling computational cost. Richardson [1] and Gaudin [2] have shown that the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian of superconductivity [3] is integrable and that the Bethe Ansatz state corresponds to a product wave function of generalized Cooper pairs [4]. In this formulation, the Bethe ansatz state is equivalent to a non-orthogonal geminal, which is widely used in quantum chemistry to describe electron pair correlations in atoms and molecules [5]. Inspired by Richardson-Gaudin integrable systems, it is possible to construct classes of non-orthogonal geminals, allowing for an accurate description of realistic and non-integrable strongly correlated atoms and molecules, while retaining the favorable computational scaling [6, 7]. In this contribution, I will discuss Richardson-Gaudin integrability and why it is useful for strongly correlated systems.