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:2008.12713 · The Journal of Chemical Physics

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

Charles‐Émile Fecteau, Hubert Fortin, Samuel Cloutier, Paul A. Johnson

Published Oct 28, 2020• 34 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1063/5.0027393

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

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

Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach

Diego R. Alcoba, Ofelia B. Oña, Luis Laín +6 more

This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.

The Journal of Chemical Physics•Jun 10, 2021•12 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

Connections between Richardson–Gaudin states, perfect-pairing, and pair coupled-cluster theory

Paul A. Johnson, Charles‐Émile Fecteau, Samuel Nadeau +2 more

Slater determinants underpin most electronic structure methods, but orbital-based approaches often struggle to describe strong correlation efficiently. Geminal-based theories, by contrast, naturally capture static correlation in bond-breaking and multi-reference problems, though at the expense of implementation complexity and limited treatment of dynamic effects. In this work, we examine the interplay between orbital and geminal frameworks, focusing on perfect-pairing (PP) wavefunctions and their relation to pair coupled-cluster doubles (pCCD) and Richardson-Gaudin states. We show that PP arises as an eigenvector of a simplified, reduced Bardeen-Cooper-Schrieffer Hamiltonian expressed in bonding/antibonding orbital pairs, with the complementary eigenvectors enabling a systematic treatment of weak correlation. Second-order Epstein-Nesbet perturbation theory on top of PP is found to yield energies nearly equivalent to pCCD. These results clarify the role of pair-based ansätze and open avenues for hybrid approaches that combine the strengths of orbital- and geminal-based methods.

arXiv:2510.06144•The Journal of Chemical Physics•Dec 17, 2025•4 citations
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.
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.
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.
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.
This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.
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.
Slater determinants underpin most electronic structure methods, but orbital-based approaches often struggle to describe strong correlation efficiently. Geminal-based theories, by contrast, naturally capture static correlation in bond-breaking and multi-reference problems, though at the expense of implementation complexity and limited treatment of dynamic effects. In this work, we examine the interplay between orbital and geminal frameworks, focusing on perfect-pairing (PP) wavefunctions and their relation to pair coupled-cluster doubles (pCCD) and Richardson-Gaudin states. We show that PP arises as an eigenvector of a simplified, reduced Bardeen-Cooper-Schrieffer Hamiltonian expressed in bonding/antibonding orbital pairs, with the complementary eigenvectors enabling a systematic treatment of weak correlation. Second-order Epstein-Nesbet perturbation theory on top of PP is found to yield energies nearly equivalent to pCCD. These results clarify the role of pair-based ansätze and open avenues for hybrid approaches that combine the strengths of orbital- and geminal-based methods.