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

Richardson-Gaudin models and broken integrability

Pieter W. Claeys

Published Sep 11, 2018• 16 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.1809.04447

Related papers

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

Eigenvalue-based method and form-factor determinant representations for integrable XXZ Richardson-Gaudin models

Pieter W. Claeys, Stijn De Baerdemacker, Mario Van Raemdonck +1 more

We propose an extension of the numerical approach for integrable Richardson-Gaudin models based on a new set of eigenvalue-based variables [A. Faribault et al., Phys. Rev. B 83, 235124 (2011); O. El Araby et al., Phys. Rev. B 85, 115130 (2012)]. Starting solely from the Gaudin algebra, the approach is generalized towards the full class of XXZ Richardson-Gaudin models. This allows for a fast and robust numerical determination of the spectral properties of these models, avoiding the singularities usually arising at the so-called singular points. We also provide different determinant expressions for the normalization of the Bethe ansatz states and form factors of local spin operators, opening up possibilities for the study of larger systems, both integrable and nonintegrable. These expressions can be written in terms of the new set of variables and generalize the results previously obtained for rational Richardson-Gaudin models [A. Faribault and D. Schuricht, J. Phys. A 45, 485202 (2012)] and Dicke-Jaynes-Cummings-Gaudin models [H. Tschirhart and A. Faribault, J. Phys. A 47, 405204 (2014)]. Remarkably, these results are independent of the explicit parametrization of the Gaudin algebra, exposing a universality in the properties of Richardson-Gaudin integrable systems deeply linked to the underlying algebraic structure.

arXiv:1501.05827•Physical Review B•Apr 3, 2015•71 citations

Variational method for integrability-breaking Richardson-Gaudin models

Pieter W. Claeys, Jean-Sébastien Caux, Dimitri Van Neck +1 more

We present a variational method for approximating the ground state of spin models close to (Richardson-Gaudin) integrability. This is done by variationally optimizing eigenstates of integrable Richardson-Gaudin models, where the toolbox of integrability allows for an efficient evaluation and minimization of the energy functional. The method is shown to return exact results for integrable models and improve substantially on perturbation theory for models close to integrability. For large integrability-breaking interactions, it is shown how (avoided) level crossings necessitate the use of excited states of integrable Hamiltonians in order to accurately describe the ground states of general nonintegrable models.

arXiv:1707.06793•Physical review. B./Physical review. B•Oct 30, 2017•11 citations

Kaleidoscope Yang-Baxter Equation for Gaudin's Kaleidoscope models

Wen-Jie Qiu, Xi-Wen Guan, Yi-Cong Yu

Recently, researchers have proposed the Asymmetric Bethe ansatz method - a theoretical tool that extends the scope of Bethe ansatz-solvable models by "breaking" partial mirror symmetry via the introduction of a fully reflecting boundary. Within this framework, the integrability conditions which were originally put forward by Gaudin have been further generalized. In this work, building on Gaudin's generalized kaleidoscope model, we present a detailed investigation of the relationship between DN symmetry and its integrability. We demonstrate that the mathematical essence of integrability in this class of models is characterized by a newly proposed Kaleidoscope Yang-Baxter Equation. Furthermore, we show that the solvability of a model via the coordinate Bethe ansatz depends not only on the consistency relations satisfied by scattering matrices, but also on the model's boundary conditions and the symmetry of the subspace where solutions are sought. Through finite element method based numerical studies, we further confirm that Bethe ansatz integrability arises in a specific symmetry sector. Finally, by analyzing the algebraic structure of the Kaleidoscope Yang-Baxter Equation, we derive a series of novel quantum algebraic identities within the framework of quantum torus algebra.

nlin.SICross-listedarXiv:2601.13596•Jan 21, 2026•0 citations

Introduction to Part I

V. E. Korepin, N. M. Bogoliubov, Anatoli G. Izergin

A method of solution of a number of quantum field theory and statistical mechanics models in two space-time dimensions is presented in this Part. This method was first suggested by H. Bethe in 1931 and is traditionally called the Bethe Ansatz. Later on the method was developed by Hulthen, Yang and Yang, Lieb, Sutherland, Baxter, Gaudin and others (see, and). We begin the presentation with the coordinate Bethe Ansatz not only due to historical reasons but also because of its simplicity and clarity. The multi-particle scattering matrix appears to be equal to the product of two-particle matrices for integrable models. This property of two-particle reducibility is of primary importance when constructing the Bethe wave function. The important feature of integrable models is that there is no mass-shell multiple particle production. This property is closely connected to the existence of an infinite number of conservation laws in such models; this will be clear from Part II. Four main models, namely the one-dimensional Bose gas, the Heisenberg magnet, the massive Thirring model and the Hubbard model, are considered in Part I. Eigenfunctions of the Hamiltonians of these models are constructed. Imposing periodic boundary conditions leads to a system of equations for the permitted values of momenta. These are known as the Bethe equations. This system can also be derived from a certain variational principle, the corresponding action being called the Yang-Yang action.

Cambridge University Press eBooks•Aug 19, 1993•1 citation

The Roaming Bethe Roots: An Effective Bethe Ansatz Beyond Integrability

Wenlong Zhao, Yunfeng Jiang, Rui-Dong Zhu

Abstract We propose an effective Bethe ansatz (EBA) for solving quantum many-body systems near an integrable point. Our approach retains the functional form of the Bethe wave function while renormalizing the Bethe roots to account for integrability-breaking interactions. These effective roots are determined by minimizing physically motivated cost functions. The resulting off-shell Bethe states serve as approximate eigenstates of the non-integrable models. We assess the quality of the approximation using various physical observables, including the energy eigenvalue, state fidelity, and bipartite entanglement entropy. Our tests show that for models with weak integrability-breaking, the EBA provides a high-quality approximation to the exact eigenstates over a wide range of deformation parameters. In contrast, for models with strong integrability-breaking interactions, the efficacy of EBA degrades relatively quickly as the deformation parameter increases. The efficacy of the method thus offers a useful probe for characterizing the strength of integrability breaking. Within its regime of accuracy, it also provides a new representation of the eigenstates of nearly integrable models, enabling one to exploit the algebraic structure inherited from integrability.

arXiv:2604.04627•Chinese Physics Letters•Aug 28, 2026•0 citations
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.
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.
We propose an extension of the numerical approach for integrable Richardson-Gaudin models based on a new set of eigenvalue-based variables [A. Faribault et al., Phys. Rev. B 83, 235124 (2011); O. El Araby et al., Phys. Rev. B 85, 115130 (2012)]. Starting solely from the Gaudin algebra, the approach is generalized towards the full class of XXZ Richardson-Gaudin models. This allows for a fast and robust numerical determination of the spectral properties of these models, avoiding the singularities usually arising at the so-called singular points. We also provide different determinant expressions for the normalization of the Bethe ansatz states and form factors of local spin operators, opening up possibilities for the study of larger systems, both integrable and nonintegrable. These expressions can be written in terms of the new set of variables and generalize the results previously obtained for rational Richardson-Gaudin models [A. Faribault and D. Schuricht, J. Phys. A 45, 485202 (2012)] and Dicke-Jaynes-Cummings-Gaudin models [H. Tschirhart and A. Faribault, J. Phys. A 47, 405204 (2014)]. Remarkably, these results are independent of the explicit parametrization of the Gaudin algebra, exposing a universality in the properties of Richardson-Gaudin integrable systems deeply linked to the underlying algebraic structure.
We present a variational method for approximating the ground state of spin models close to (Richardson-Gaudin) integrability. This is done by variationally optimizing eigenstates of integrable Richardson-Gaudin models, where the toolbox of integrability allows for an efficient evaluation and minimization of the energy functional. The method is shown to return exact results for integrable models and improve substantially on perturbation theory for models close to integrability. For large integrability-breaking interactions, it is shown how (avoided) level crossings necessitate the use of excited states of integrable Hamiltonians in order to accurately describe the ground states of general nonintegrable models.
Recently, researchers have proposed the Asymmetric Bethe ansatz method - a theoretical tool that extends the scope of Bethe ansatz-solvable models by "breaking" partial mirror symmetry via the introduction of a fully reflecting boundary. Within this framework, the integrability conditions which were originally put forward by Gaudin have been further generalized. In this work, building on Gaudin's generalized kaleidoscope model, we present a detailed investigation of the relationship between DN symmetry and its integrability. We demonstrate that the mathematical essence of integrability in this class of models is characterized by a newly proposed Kaleidoscope Yang-Baxter Equation. Furthermore, we show that the solvability of a model via the coordinate Bethe ansatz depends not only on the consistency relations satisfied by scattering matrices, but also on the model's boundary conditions and the symmetry of the subspace where solutions are sought. Through finite element method based numerical studies, we further confirm that Bethe ansatz integrability arises in a specific symmetry sector. Finally, by analyzing the algebraic structure of the Kaleidoscope Yang-Baxter Equation, we derive a series of novel quantum algebraic identities within the framework of quantum torus algebra.
A method of solution of a number of quantum field theory and statistical mechanics models in two space-time dimensions is presented in this Part. This method was first suggested by H. Bethe in 1931 and is traditionally called the Bethe Ansatz. Later on the method was developed by Hulthen, Yang and Yang, Lieb, Sutherland, Baxter, Gaudin and others (see, and). We begin the presentation with the coordinate Bethe Ansatz not only due to historical reasons but also because of its simplicity and clarity. The multi-particle scattering matrix appears to be equal to the product of two-particle matrices for integrable models. This property of two-particle reducibility is of primary importance when constructing the Bethe wave function. The important feature of integrable models is that there is no mass-shell multiple particle production. This property is closely connected to the existence of an infinite number of conservation laws in such models; this will be clear from Part II. Four main models, namely the one-dimensional Bose gas, the Heisenberg magnet, the massive Thirring model and the Hubbard model, are considered in Part I. Eigenfunctions of the Hamiltonians of these models are constructed. Imposing periodic boundary conditions leads to a system of equations for the permitted values of momenta. These are known as the Bethe equations. This system can also be derived from a certain variational principle, the corresponding action being called the Yang-Yang action.
Abstract We propose an effective Bethe ansatz (EBA) for solving quantum many-body systems near an integrable point. Our approach retains the functional form of the Bethe wave function while renormalizing the Bethe roots to account for integrability-breaking interactions. These effective roots are determined by minimizing physically motivated cost functions. The resulting off-shell Bethe states serve as approximate eigenstates of the non-integrable models. We assess the quality of the approximation using various physical observables, including the energy eigenvalue, state fidelity, and bipartite entanglement entropy. Our tests show that for models with weak integrability-breaking, the EBA provides a high-quality approximation to the exact eigenstates over a wide range of deformation parameters. In contrast, for models with strong integrability-breaking interactions, the efficacy of EBA degrades relatively quickly as the deformation parameter increases. The efficacy of the method thus offers a useful probe for characterizing the strength of integrability breaking. Within its regime of accuracy, it also provides a new representation of the eigenstates of nearly integrable models, enabling one to exploit the algebraic structure inherited from integrability.