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.
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.
In this work, we generalize the numerical approach to Gaudin models developed earlier by us [Faribault, El Araby, Str\"ater, and Gritsev, Phys. Rev. B 83, 235124 (2011)] to degenerate systems, showing that their treatment is surprisingly convenient from a numerical point of view. In fact, high degeneracies not only reduce the number of relevant states in the Hilbert space by a non-negligible fraction, they also allow us to write the relevant equations in the form of sparse matrix equations. Moreover, we introduce an inversion method based on a basis of barycentric polynomials that leads to a more stable and efficient root extraction, which most importantly avoids the necessity of working with arbitrary precision. As an example, we show the results of our procedure applied to the Richardson model on a square lattice.
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.
Abstract In this work we construct the eigenstates of the most general spin-1/2 Richardson–Gaudin model integrable in an external magnetic field. This includes the possibility for fully anisotropic XYZ coupling such that the S i x S j x , S i y S j y and S i z S j z terms all have distinct coupling strengths. While insuring that integrability is maintained in the presence of an external field excludes the elliptic XYZ model which is only integrable at zero field, this work still covers a wide class of fully anisotropic (XYZ) models associated with non skew-symmetric r -matrices. The eigenstates, as constructed here, do not require any usable Bethe ansatz and therefore: no proper pseudo-vacuum, Bethe roots, or generalised spin raising (Gaudin) operators have to be defined. Indeed, the eigenstates are generically built only through the conserved charges which define the model of interest and the specification of the set of eigenvalues defining the particular eigenstate. Since these eigenvalues are, in general, solutions to a simple set of quadratic equations, the proposed approach is simpler to implement than any Bethe ansatz and, moreover, it remains completely identical independently of the symmetries of the model. Indeed, the construction removes any distinction between XYZ models and XXZ/XXX models and, generically, that between models with or without U (1) symmetry so that any difficulties associated with the use of a Bethe ansatz in any of these cases are avoided.
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.
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.