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

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

Stijn De Baerdemacker, Paul W. Ayers, Patrick Bultinck, Paul A. Johnson, Peter A. Limacher, Dimitri Van Neck, Mario Van Raemdonck

Published Jan 1, 2013• 0 citations
OpenAlex

Abstract

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.

Related papers

Richardson-Gaudin models and broken integrability

Pieter W. Claeys

This thesis presents an introduction to the class of Richardson-Gaudin integrable models, with special focus on the Bethe ansatz wave function, and investigates ways of applying the properties of Richardson-Gaudin models both in and out of integrability. A framework is outlined for the numerical and theoretical treatment of these systems, exposing a duality allowing the Bethe equations to be solved numerically. This is extended to the calculation of inner products and correlation functions. Using this framework, the influence of particle exchange on the Bethe ansatz is discussed, after which it is shown how the Bethe ansatz is able to accurately model wave functions of non-integrable models in two different settings. First, a variational approach is outlined for stationary models where integrability-breaking perturbations are explicitly introduced. Second, an alternative way of breaking integrability is through the introduction of dynamics and periodic driving, where it is shown how integrability can be used to model the resulting Floquet many-body resonances. Throughout this work, it is shown how the clear-cut structure and relatively large freedom in Richardson-Gaudin models makes them ideal for an investigation of the general principles of integrability, as well as being a perfect testing ground for the development of new quantum many-body techniques beyond integrability.

arXiv:1809.04447•Sep 11, 2018•16 citations

Superconducting correlations in metallic nanoparticles: Exact solution of the BCS model by the algebraic Bethe ansatz

Huan-Qiang Zhou, Jon Links, Ross H. McKenzie +1 more

Superconducting pairing of electrons in nanoscale metallic particles with discrete energy levels and a fixed number of electrons is described by the reduced Bardeen, Cooper, and Schrieffer model Hamiltonian. We show that this model is integrable by the algebraic Bethe ansatz. The eigenstates, spectrum, conserved operators, integrals of motion, and norms of wave functions are obtained. Furthermore, the quantum inverse problem is solved, meaning that form factors and correlation functions can be explicitly evaluated. Closed form expressions are given for the form factors and correlation functions that describe superconducting pairing.

arXiv:cond-mat/0106390•Physical review. B, Condensed matter•Jan 11, 2002•98 citations

Colloquium : Exactly solvable Richardson-Gaudin models for many-body quantum systems

J. Dukelsky, Stuart Pittel, Germán Sierra

The use of exactly solvable Richardson-Gaudin models to describe the physics of systems with strong pair correlations is reviewed. The article begins with a brief discussion of Richardson's early work, which demonstrated the exact solvability of the pure pairing model, and then shows how that work has evolved recently into a much richer class of exactly solvable models. The Richardson solution leads naturally to an exact analogy between these quantum models and classical electrostatic problems in two dimensions. This analogy is then used to demonstrate formally how BCS theory emerges as the large-$N$ limit of the pure pairing Hamiltonian. Several applications to problems of relevance to condensed-matter physics, nuclear physics, and the physics of confined systems are considered. Some of the interesting effects that are discussed in the context of these exactly solvable models include (i) the crossover from superconductivity to a fluctuation-dominated regime in small metallic grains; (ii) the role of the nucleon Pauli principle in suppressing the effects of high-spin bosons in interacting boson models of nuclei, and (iii) the possibility of fragmentation in confined boson systems. Interesting insight is also provided into the origin of the superconducting phase transition both in two-dimensional electronic systems and in atomic nuclei, based on the electrostatic image of the corresponding exactly solvable quantum pairing models.

arXiv:nucl-th/0405011•Reviews of Modern Physics•Aug 6, 2004•407 citations

THE EXACTLY SOLVABLE RICHARDSON MODEL IN THE BCS-to-BEC CROSSOVER

J. Dukelsky, Gerardo Ortíz

The exact Richardson solution of the reduced BCS Hamiltonian is reviewed with special emphasis on the structure of the many-body wavefunction. A clear definition of the Cooper pairs in the superconducting or Fermi superfluid medium emerges. The structure of the Cooper pair wavefunction is general enough to describe large overlapping pair resonances in the weak coupling limit (BCS) as well as a Bose-Einstein condensate (BEC) of quasi-bound molecules in the strong coupling limit.

International Journal of Modern Physics E•Feb 20, 2006•3 citations

Non-skew-symmetric classical r-matrices, algebraic Bethe ansatz, and Bardeen–Cooper–Schrieffer–type integrable systems

T. Skrypnyk

We construct quantum integrable systems associated with non-skew-symmetric gl(2)-valued classical r-matrices. We find a new explicit multiparametric family of such the non-skew-symmetric classical r-matrices. We consider two classes of examples of the corresponding integrable systems, namely generalized Gaudin systems with and without an external magnetic field. In the case of arbitrary r-matrices diagonal in a standard gl(2)-basis, we calculate the spectrum of the corresponding quantum integrable systems using the algebraic Bethe ansatz. We apply these results to a construction of integrable fermionic models and obtain a wide class of integrable Bardeen–Cooper–Schrieffer (BCS)-type fermionic Hamiltonians containing the pairing and electrostatic interaction terms. We also consider special cases when the corresponding integrable Hamiltonians contain only pairing interaction term and are exact analogs of the “reduced BCS Hamiltonian” of Richardson.

Journal of Mathematical Physics•Mar 1, 2009•49 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
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.
This thesis presents an introduction to the class of Richardson-Gaudin integrable models, with special focus on the Bethe ansatz wave function, and investigates ways of applying the properties of Richardson-Gaudin models both in and out of integrability. A framework is outlined for the numerical and theoretical treatment of these systems, exposing a duality allowing the Bethe equations to be solved numerically. This is extended to the calculation of inner products and correlation functions. Using this framework, the influence of particle exchange on the Bethe ansatz is discussed, after which it is shown how the Bethe ansatz is able to accurately model wave functions of non-integrable models in two different settings. First, a variational approach is outlined for stationary models where integrability-breaking perturbations are explicitly introduced. Second, an alternative way of breaking integrability is through the introduction of dynamics and periodic driving, where it is shown how integrability can be used to model the resulting Floquet many-body resonances. Throughout this work, it is shown how the clear-cut structure and relatively large freedom in Richardson-Gaudin models makes them ideal for an investigation of the general principles of integrability, as well as being a perfect testing ground for the development of new quantum many-body techniques beyond integrability.
Superconducting pairing of electrons in nanoscale metallic particles with discrete energy levels and a fixed number of electrons is described by the reduced Bardeen, Cooper, and Schrieffer model Hamiltonian. We show that this model is integrable by the algebraic Bethe ansatz. The eigenstates, spectrum, conserved operators, integrals of motion, and norms of wave functions are obtained. Furthermore, the quantum inverse problem is solved, meaning that form factors and correlation functions can be explicitly evaluated. Closed form expressions are given for the form factors and correlation functions that describe superconducting pairing.
The use of exactly solvable Richardson-Gaudin models to describe the physics of systems with strong pair correlations is reviewed. The article begins with a brief discussion of Richardson's early work, which demonstrated the exact solvability of the pure pairing model, and then shows how that work has evolved recently into a much richer class of exactly solvable models. The Richardson solution leads naturally to an exact analogy between these quantum models and classical electrostatic problems in two dimensions. This analogy is then used to demonstrate formally how BCS theory emerges as the large-NNN limit of the pure pairing Hamiltonian. Several applications to problems of relevance to condensed-matter physics, nuclear physics, and the physics of confined systems are considered. Some of the interesting effects that are discussed in the context of these exactly solvable models include (i) the crossover from superconductivity to a fluctuation-dominated regime in small metallic grains; (ii) the role of the nucleon Pauli principle in suppressing the effects of high-spin bosons in interacting boson models of nuclei, and (iii) the possibility of fragmentation in confined boson systems. Interesting insight is also provided into the origin of the superconducting phase transition both in two-dimensional electronic systems and in atomic nuclei, based on the electrostatic image of the corresponding exactly solvable quantum pairing models.
We construct quantum integrable systems associated with non-skew-symmetric gl(2)-valued classical r-matrices. We find a new explicit multiparametric family of such the non-skew-symmetric classical r-matrices. We consider two classes of examples of the corresponding integrable systems, namely generalized Gaudin systems with and without an external magnetic field. In the case of arbitrary r-matrices diagonal in a standard gl(2)-basis, we calculate the spectrum of the corresponding quantum integrable systems using the algebraic Bethe ansatz. We apply these results to a construction of integrable fermionic models and obtain a wide class of integrable Bardeen–Cooper–Schrieffer (BCS)-type fermionic Hamiltonians containing the pairing and electrostatic interaction terms. We also consider special cases when the corresponding integrable Hamiltonians contain only pairing interaction term and are exact analogs of the “reduced BCS Hamiltonian” of Richardson.
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.