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:1206.2988 · Journal of Statistical Mechanics Theory and Experiment

Generalized Heine–Stieltjes and Van Vleck polynomials associated with two-level, integrable BCS models

Ian Marquette, Jon Links

Published Aug 31, 2012• 36 citations
arXivHTMLPDF

Abstract

We study the Bethe ansatz/ordinary differential equation (BA/ODE) correspondence for Bethe ansatz equations that belong to a certain class of coupled, nonlinear, algebraic equations. Through this approach we numerically obtain the generalized Heine–Stieltjes and Van Vleck polynomials in the degenerate, two-level limit for four cases of integrable Bardeen–Cooper–Schrieffer (BCS) pairing models. These are the s-wave pairing model, the p + ip-wave pairing model, the p + ip pairing model coupled to a bosonic molecular pair degree of freedom, and a newly introduced extended d + id-wave pairing model with additional interactions. The zeros of the generalized Heine–Stieltjes polynomials provide solutions of the corresponding Bethe ansatz equations. We compare the roots of the ground states with curves obtained from the solution of a singular integral equation approximation, which allows for a characterization of ground-state phases in these systems. Our techniques also permit the computation of the roots of the excited states. These results illustrate how the BA/ODE correspondence can be used to provide new numerical methods to study a variety of integrable systems.

DOI
10.1088/1742-5468/2012/08/p08019

Related papers

Heine-Stieltjes correspondence and the polynomial approach to the standard pairing problem

Xin Guan, Kristina D. Launey, Mingxia Xie +3 more

A new approach for solving the Bethe ansatz (Gaudin-Richardson) equations of the standard pairing problem is established based on the Heine-Stieltjes correspondence. For $k$ pairs of valence nucleons on $n$ different single-particle levels, it is found that solutions of the Bethe ansatz equations can be obtained from one $(k+1)\ifmmode\times\else\texttimes\fi{}(k+1)$ and one $(n\ensuremath{-}1)\ifmmode\times\else\texttimes\fi{}(k+1)$ matrices, which are associated with the extended Heine-Stieltjes and Van Vleck polynomials, respectively. Since the coefficients in these polynomials are free from divergence with variations in contrast to the original Bethe ansatz equations, the approach provides an efficient and systematic way to solve the problem, which by extension, can also be used to solve a large class of Gaudin-type quantum many-body problems, including an efficient angular momentum projection method for multiparticle systems.

arXiv:1106.5237•Physical Review C•Aug 28, 2012•79 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

The Heine–Stieltjes correspondence and the polynomial approach to the Gaudin–Richardson models

Feng Pan, Kristina D. Launey, Xin Guan

The Heine–Stieltjes correspondence is extended and applied to solve Bethe ansatz equations of the Gaudin–Richardson models, from which the extended Heine–Stieltjes polynomial approach to these models is proposed. As examples for the application of this approach, exact solutions of the standard two-site Bose–Hubbard model and the standard pairing model for nuclei are formulated from the corresponding polynomials.

Journal of Physics Conference Series•Dec 18, 2012•0 citations

Exact solution of the p + ip pairing Hamiltonian and a hierarchy of integrable models

Clare Dunning, Miguel Ibáñez-Berganza, Jon Links +2 more

Using the well-known trigonometric six-vertex solution of the Yang-Baxter equation we derive an integrable pairing Hamiltonian with anyonic degrees of freedom. The exact algebraic Bethe ansatz solution is obtained using standard techniques. From this model we obtain several limiting models, including the pairing Hamiltonian with p+ip-wave symmetry. An in-depth study of the p+ip model is then undertaken, including a mean-field analysis, analytic and numerical solution of the Bethe ansatz equations, and an investigation of the topological properties of the ground-state wavefunction. Our main result is that the ground-state phase diagram of the p+ip model consists of three phases. There is the known boundary line with gapless excitations that occurs for vanishing chemical potential, separating the topologically trivial strong pairing phase and the topologically non-trivial weak pairing phase. We argue that a second boundary line exists separating the weak pairing phase from a topologically trivial weak coupling BCS phase, which includes the Fermi sea in the limit of zero coupling. The ground state on this second boundary line is the Moore-Read state.

arXiv:1001.1591•Journal of Statistical Mechanics Theory and Experiment•Aug 26, 2010•57 citations

Gaudin models solver based on the correspondence between Bethe ansatz and ordinary differential equations

Alexandre Faribault, Omar El Araby, Christoph Sträter +1 more

We present a numerical approach which allows the solving of Bethe equations whose solutions define the eigenstates of Gaudin models. By focusing on a different set of variables, the canceling divergences which occur for certain values of the coupling strength no longer appear explicitly. The problem is thus reduced to a set of quadratic algebraic equations. The required inverse transformation can then be realized using only linear operations and a standard polynomial root-finding algorithm. The method is applied to Richardson's fermionic pairing model, the central spin model, and the generalized Dicke model.

arXiv:1103.0472•Physical Review B•Jun 13, 2011•100 citations

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
We study the Bethe ansatz/ordinary differential equation (BA/ODE) correspondence for Bethe ansatz equations that belong to a certain class of coupled, nonlinear, algebraic equations. Through this approach we numerically obtain the generalized Heine–Stieltjes and Van Vleck polynomials in the degenerate, two-level limit for four cases of integrable Bardeen–Cooper–Schrieffer (BCS) pairing models. These are the s-wave pairing model, the p + ip-wave pairing model, the p + ip pairing model coupled to a bosonic molecular pair degree of freedom, and a newly introduced extended d + id-wave pairing model with additional interactions. The zeros of the generalized Heine–Stieltjes polynomials provide solutions of the corresponding Bethe ansatz equations. We compare the roots of the ground states with curves obtained from the solution of a singular integral equation approximation, which allows for a characterization of ground-state phases in these systems. Our techniques also permit the computation of the roots of the excited states. These results illustrate how the BA/ODE correspondence can be used to provide new numerical methods to study a variety of integrable systems.
A new approach for solving the Bethe ansatz (Gaudin-Richardson) equations of the standard pairing problem is established based on the Heine-Stieltjes correspondence. For kkk pairs of valence nucleons on nnn different single-particle levels, it is found that solutions of the Bethe ansatz equations can be obtained from one $(k+1)\ifmmode\times\else\texttimes\fi{}(k+1)$ and one $(n\ensuremath{-}1)\ifmmode\times\else\texttimes\fi{}(k+1)$ matrices, which are associated with the extended Heine-Stieltjes and Van Vleck polynomials, respectively. Since the coefficients in these polynomials are free from divergence with variations in contrast to the original Bethe ansatz equations, the approach provides an efficient and systematic way to solve the problem, which by extension, can also be used to solve a large class of Gaudin-type quantum many-body problems, including an efficient angular momentum projection method for multiparticle systems.
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.
Using the well-known trigonometric six-vertex solution of the Yang-Baxter equation we derive an integrable pairing Hamiltonian with anyonic degrees of freedom. The exact algebraic Bethe ansatz solution is obtained using standard techniques. From this model we obtain several limiting models, including the pairing Hamiltonian with p+ip-wave symmetry. An in-depth study of the p+ip model is then undertaken, including a mean-field analysis, analytic and numerical solution of the Bethe ansatz equations, and an investigation of the topological properties of the ground-state wavefunction. Our main result is that the ground-state phase diagram of the p+ip model consists of three phases. There is the known boundary line with gapless excitations that occurs for vanishing chemical potential, separating the topologically trivial strong pairing phase and the topologically non-trivial weak pairing phase. We argue that a second boundary line exists separating the weak pairing phase from a topologically trivial weak coupling BCS phase, which includes the Fermi sea in the limit of zero coupling. The ground state on this second boundary line is the Moore-Read state.
We present a numerical approach which allows the solving of Bethe equations whose solutions define the eigenstates of Gaudin models. By focusing on a different set of variables, the canceling divergences which occur for certain values of the coupling strength no longer appear explicitly. The problem is thus reduced to a set of quadratic algebraic equations. The required inverse transformation can then be realized using only linear operations and a standard polynomial root-finding algorithm. The method is applied to Richardson's fermionic pairing model, the central spin model, and the generalized Dicke model.
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.