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:1912.05692 · SciPost Physics Core

Ground-state energy of a Richardson-Gaudin integrable BCS model

Yibing Shen, Phillip S. Isaac, Jon Links

Published Feb 4, 2020• 7 citations
arXivHTMLPDF

Abstract

We investigate the ground-state energy of a Richardson-Gaudin integrable BCS model, generalizing the closed and open p+ip models. The Hamiltonian supports a family of mutually commuting conserved operators satisfying quadratic relations. From the eigenvalues of the conserved operators we derive, in the continuum limit, an integral equation for which a solution corresponding to the ground state is established. The energy expression from this solution agrees with the BCS mean-field result.

DOI
10.21468/scipostphyscore.2.1.001

Related papers

Ground-state energies of the open and closed p + ip-pairing models from the Bethe Ansatz

Yibing Shen, Phillip S. Isaac, Jon Links

Using the exact Bethe Ansatz solution, we investigate methods for calculating the ground-state energy for the p + i p -pairing Hamiltonian. We first consider the Hamiltonian isolated from its environment (closed model) through two forms of Bethe Ansatz solutions, which generally have complex-valued Bethe roots. A continuum limit approximation, leading to an integral equation, is applied to compute the ground-state energy. We discuss the evolution of the root distribution curve with respect to a range of parameters, and the limitations of this method. We then consider an alternative approach that transforms the Bethe Ansatz equations to an equivalent form, but in terms of the real-valued conserved operator eigenvalues. An integral equation is established for the transformed solution. This equation is shown to admit an exact solution associated with the ground state. Next we discuss results for a recently derived Bethe Ansatz solution of the open model. With the aforementioned alternative approach based on real-valued roots, combined with mean-field analysis, we are able to establish an integral equation with an exact solution that corresponds to the ground-state for this case.

arXiv:1807.00428•Nuclear Physics B•Aug 29, 2018•10 citations

Bethe Ansatz analysis of the open and closed p + ip-pairing BCS models

Yibing Shen

The p+ip-wave pairing BCS model has attracted much attention due to its potential applications inexperimental physics, and the developing theory of topological superconductors. In this thesis, weinvestigate the p + ip BCS model as an integrable model with exact Bethe Ansatz solutions. Theground-state energies of this model will be studied with various methods.  First a brief introduction on quantum integrability is given, followed by an introduction to theQuantum Inverse Scattering Method, algebraic Bethe Ansatz and BCS theory of superconductivity. Wethen introduce a two-dimensional BCS pairing Hamiltonian with px +ipy symmetry. In the frameworkof the Quantum Inverse Scattering Method, the p+ip Hamiltonian can be constructed from a specificsolution of the non-standard classical Yang-Baxter equation. This solution produces a family ofconserved operators satisfying quadratic relations. Remarks on topological superconductivity for thep+ip-wave pairing model are provided.  We move on to investigate methods for calculating the ground-state energy for the p+ip Hamil-tonian. We first consider the Hamiltonian isolated from its environment (closed model) through twoforms of Bethe Ansatz solution, which generally have complex-valued Bethe roots. A continuum limitapproximation, leading to an integral equation, is applied to compute the ground-state energy. Theevolution of the root distribution curve with respect to a range of parameters is presented, followed bya discussion of the limitations of this method. Despite the limitations, the ground-state energy andthe ‘gap’ equation derived from the integral approximation are consistent with mean-field results fromprevious study.  We then consider the Hamiltonian interacting with its environment (open model) and begin witha mean-field analysis. We discuss the relation between the Bethe roots and the conserved operatoreigenvalues. The Bethe Ansatz solution of the open model is studied numerically, revealing complicatedbehaviour of the Bethe roots without a clear indication for an integral approximation.  Next we consider an alternative approach based on the quadratic relations among the real-valuedconserved operator eigenvalues. An integral equation is first established for the closed model. Thisequation is shown to admit an exact solution associated with the ground state.  Finally we extend the aforementioned alternative approach to the open model. Combined withresults from mean-field analysis, we are able to establish an integral equation with an exact solutionthat corresponds to the ground state for the open case.

•Sep 27, 2019•0 citations

Quadratic operator relations and Bethe equations for spin-1/2 Richardson–Gaudin models

Claude Dimo, Alexandre Faribault

Abstract In this work we demonstrate how one can, in a generic approach, derive a set of N simple quadratic Bethe equations for integrable Richardson–Gaudin (RG) models built out of N spins-1/2. These equations depend only on the N eigenvalues of the various conserved charges so that any solution of these equations defines, indirectly through the corresponding set of eigenvalues, one particular eigenstate. The proposed construction covers the full class of integrable RG models of the XYZ (including the subclasses of XXZ and XXX models) type realised in terms of spins-1/2, coupled with one another through , , terms, including, as well, magnetic field-like terms linear in the Pauli matrices. The approach exclusively requires integrability, defined here only by the requirement that N conserved charges R i (with ) such that exist. The result is therefore valid, and equally simple, for models with or without symmetry, with or without a properly defined pseudo-vacuum as well as for models with non-skew symmetric couplings.

arXiv:1805.03427•Journal of Physics A Mathematical and Theoretical•Jun 14, 2018•16 citations

Integrable spin- 1 2 Richardson–Gaudin XYZ models in an arbitrary magnetic field

Pieter W. Claeys, Claude Dimo, Stijn De Baerdemacker +1 more

Abstract We establish the most general class of spin- integrable Richardson–Gaudin models including an arbitrary magnetic field, returning a fully anisotropic (XYZ) model. The restriction to spin- relaxes the usual integrability constraints, allowing for a general solution where the couplings between spins lack the usual antisymmetric properties of Richardson–Gaudin models. The full set of conserved charges are constructed explicitly and shown to satisfy a set of quadratic equations, allowing for the numerical treatment of a fully anisotropic central spin in an external magnetic field. While this approach does not provide expressions for the exact eigenstates, it allows their eigenvalues to be obtained, and expectation values of local observables can then be calculated from the Hellmann–Feynman theorem.

arXiv:1810.06059•Journal of Physics A Mathematical and Theoretical•Jan 15, 2019•29 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

Energy of N Cooper pair by analytically solving Richardson-Gaudin equations

Michel Crouzeix, Monique Combescot

This Letter provides the solution to a yet unsolved basic problem of Solid State Physics: the ground state energy of an arbitrary number of Cooper pairs interacting via the Bardeen-Cooper-Schrieffer potential. We here break a 50 year old math problem by analytically solving Richardson-Gaudin equations which give the exact energy of these $N$ pairs via $N$ parameters coupled through $N$ non-linear equations. Our result fully supports the standard BCS result obtained for a pair number equal to half the number of states feeling the potential. More importantly, it shows that the interaction part of the $N$-pair energy depends on $N$ as $N(N-1)$ only from N=1 to the dense regime, a result which evidences that Cooper pairs interact via Pauli blocking only.

arXiv:1111.5801•Nov 24, 2011•0 citations
We investigate the ground-state energy of a Richardson-Gaudin integrable BCS model, generalizing the closed and open p+ip models. The Hamiltonian supports a family of mutually commuting conserved operators satisfying quadratic relations. From the eigenvalues of the conserved operators we derive, in the continuum limit, an integral equation for which a solution corresponding to the ground state is established. The energy expression from this solution agrees with the BCS mean-field result.
Using the exact Bethe Ansatz solution, we investigate methods for calculating the ground-state energy for the p + i p -pairing Hamiltonian. We first consider the Hamiltonian isolated from its environment (closed model) through two forms of Bethe Ansatz solutions, which generally have complex-valued Bethe roots. A continuum limit approximation, leading to an integral equation, is applied to compute the ground-state energy. We discuss the evolution of the root distribution curve with respect to a range of parameters, and the limitations of this method. We then consider an alternative approach that transforms the Bethe Ansatz equations to an equivalent form, but in terms of the real-valued conserved operator eigenvalues. An integral equation is established for the transformed solution. This equation is shown to admit an exact solution associated with the ground state. Next we discuss results for a recently derived Bethe Ansatz solution of the open model. With the aforementioned alternative approach based on real-valued roots, combined with mean-field analysis, we are able to establish an integral equation with an exact solution that corresponds to the ground-state for this case.
The p+ip-wave pairing BCS model has attracted much attention due to its potential applications inexperimental physics, and the developing theory of topological superconductors. In this thesis, weinvestigate the p + ip BCS model as an integrable model with exact Bethe Ansatz solutions. Theground-state energies of this model will be studied with various methods.  First a brief introduction on quantum integrability is given, followed by an introduction to theQuantum Inverse Scattering Method, algebraic Bethe Ansatz and BCS theory of superconductivity. Wethen introduce a two-dimensional BCS pairing Hamiltonian with px +ipy symmetry. In the frameworkof the Quantum Inverse Scattering Method, the p+ip Hamiltonian can be constructed from a specificsolution of the non-standard classical Yang-Baxter equation. This solution produces a family ofconserved operators satisfying quadratic relations. Remarks on topological superconductivity for thep+ip-wave pairing model are provided.  We move on to investigate methods for calculating the ground-state energy for the p+ip Hamil-tonian. We first consider the Hamiltonian isolated from its environment (closed model) through twoforms of Bethe Ansatz solution, which generally have complex-valued Bethe roots. A continuum limitapproximation, leading to an integral equation, is applied to compute the ground-state energy. Theevolution of the root distribution curve with respect to a range of parameters is presented, followed bya discussion of the limitations of this method. Despite the limitations, the ground-state energy andthe ‘gap’ equation derived from the integral approximation are consistent with mean-field results fromprevious study.  We then consider the Hamiltonian interacting with its environment (open model) and begin witha mean-field analysis. We discuss the relation between the Bethe roots and the conserved operatoreigenvalues. The Bethe Ansatz solution of the open model is studied numerically, revealing complicatedbehaviour of the Bethe roots without a clear indication for an integral approximation.  Next we consider an alternative approach based on the quadratic relations among the real-valuedconserved operator eigenvalues. An integral equation is first established for the closed model. Thisequation is shown to admit an exact solution associated with the ground state.  Finally we extend the aforementioned alternative approach to the open model. Combined withresults from mean-field analysis, we are able to establish an integral equation with an exact solutionthat corresponds to the ground state for the open case.
Abstract In this work we demonstrate how one can, in a generic approach, derive a set of N simple quadratic Bethe equations for integrable Richardson–Gaudin (RG) models built out of N spins-1/2. These equations depend only on the N eigenvalues of the various conserved charges so that any solution of these equations defines, indirectly through the corresponding set of eigenvalues, one particular eigenstate. The proposed construction covers the full class of integrable RG models of the XYZ (including the subclasses of XXZ and XXX models) type realised in terms of spins-1/2, coupled with one another through , , terms, including, as well, magnetic field-like terms linear in the Pauli matrices. The approach exclusively requires integrability, defined here only by the requirement that N conserved charges R i (with ) such that exist. The result is therefore valid, and equally simple, for models with or without symmetry, with or without a properly defined pseudo-vacuum as well as for models with non-skew symmetric couplings.
Abstract We establish the most general class of spin- integrable Richardson–Gaudin models including an arbitrary magnetic field, returning a fully anisotropic (XYZ) model. The restriction to spin- relaxes the usual integrability constraints, allowing for a general solution where the couplings between spins lack the usual antisymmetric properties of Richardson–Gaudin models. The full set of conserved charges are constructed explicitly and shown to satisfy a set of quadratic equations, allowing for the numerical treatment of a fully anisotropic central spin in an external magnetic field. While this approach does not provide expressions for the exact eigenstates, it allows their eigenvalues to be obtained, and expectation values of local observables can then be calculated from the Hellmann–Feynman theorem.
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.
This Letter provides the solution to a yet unsolved basic problem of Solid State Physics: the ground state energy of an arbitrary number of Cooper pairs interacting via the Bardeen-Cooper-Schrieffer potential. We here break a 50 year old math problem by analytically solving Richardson-Gaudin equations which give the exact energy of these NNN pairs via NNN parameters coupled through NNN non-linear equations. Our result fully supports the standard BCS result obtained for a pair number equal to half the number of states feeling the potential. More importantly, it shows that the interaction part of the NNN-pair energy depends on NNN as N(N−1)N(N-1)N(N−1) only from N=1 to the dense regime, a result which evidences that Cooper pairs interact via Pauli blocking only.