© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1106.5237 · Physical Review C

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

Xin Guan, Kristina D. Launey, Mingxia Xie, Lina Bao, Feng Pan, Jerry Paul Draayer

Published Aug 28, 2012• 79 citations

Abstract

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.

DOI: 10.1103/physrevc.86.024313

Related papers

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

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

Ian Marquette, Jon Links

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.

arXiv:1206.2988•Journal of Statistical Mechanics Theory and Experiment•Aug 31, 2012•36 citations

Solutions of nuclear pairing

A. B. Balantekin, Yamaç Pehlivan

We give the exact solution of orbit dependent nuclear pairing problem between two nondegenerate energy levels using the Bethe ansatz technique. Our solution reduces to previously solved cases in the appropriate limits including Richardson's treatment of reduced pairing in terms of rational Gaudin algebra operators.

arXiv:0710.3941•Physical Review C•Nov 14, 2007•18 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

New exact solutions of the standard pairing model for well-deformed nuclei

Feng Pan, Mingxia Xie, Xin Guan +2 more

A new step-by-step diagonalization procedure for evaluating exact solutions of the nuclear deformed mean-field plus pairing interaction model is proposed via a simple Bethe ansatz in each step from which the eigenvalues and corresponding eigenstates can be obtained progressively. This new approach draws upon an observation that the original one- plus two-body problem in a $k$-particle Hilbert subspace can be mapped onto a one-body grand hard-core boson picture that can be solved step by step with a simple Bethe ansatz known from earlier work. Based on this new procedure, it is further shown that the extended pairing model for deformed nuclei [Feng Pan, V. G. Gueorguiev, and J. P. Draayer, Phys. Rev. Lett. 92, 112503 (2004)] is similar to the standard pairing model with the first step approximation, in which only the lowest energy eigenstate of the standard pure pairing interaction part is taken into consideration. Our analysis shows that the standard pairing model with the first step approximation displays similar pair structures of the first few exact low-lying states of the model, which, therefore, provides a link between the two models.

arXiv:0904.2830•Physical Review C•Oct 9, 2009•21 citations

Numerical algorithm for the standard pairing problem based on the Heine–Stieltjes correspondence and the polynomial approach

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

Computer Physics Communications•Jun 4, 2014•36 citations

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.

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.

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.

A new step-by-step diagonalization procedure for evaluating exact solutions of the nuclear deformed mean-field plus pairing interaction model is proposed via a simple Bethe ansatz in each step from which the eigenvalues and corresponding eigenstates can be obtained progressively. This new approach draws upon an observation that the original one- plus two-body problem in a

k

-particle Hilbert subspace can be mapped onto a one-body grand hard-core boson picture that can be solved step by step with a simple Bethe ansatz known from earlier work. Based on this new procedure, it is further shown that the extended pairing model for deformed nuclei [Feng Pan, V. G. Gueorguiev, and J. P. Draayer, Phys. Rev. Lett. 92, 112503 (2004)] is similar to the standard pairing model with the first step approximation, in which only the lowest energy eigenstate of the standard pure pairing interaction part is taken into consideration. Our analysis shows that the standard pairing model with the first step approximation displays similar pair structures of the first few exact low-lying states of the model, which, therefore, provides a link between the two models.