Starting from integrable su (2) (quasi-)spin Richardson–Gaudin (RG) XXZ models we derive several properties of integrable spin models coupled to a bosonic mode. We focus on the Dicke–Jaynes–Cummings–Gaudin models and the two-channel ( p + i p )-wave pairing Hamiltonian. The pseudo-deformation of the underlying su (2) algebra is here introduced as a way to obtain these models in the contraction limit of different RG models. This allows for the construction of the full set of conserved charges, the Bethe ansatz state, and the resulting RG equations. For these models an alternative and simpler set of quadratic equations can be found in terms of the eigenvalues of the conserved charges. Furthermore, the recently proposed eigenvalue-based determinant expressions for the overlaps and form factors of local operators are extended to these models, linking the results previously presented for the Dicke–Jaynes–Cummings–Gaudin models with the general results for RG XXZ models.
Starting from integrable su (2) (quasi-)spin Richardson–Gaudin (RG) XXZ models we derive several properties of integrable spin models coupled to a bosonic mode. We focus on the Dicke–Jaynes–Cummings–Gaudin models and the two-channel ( p + i p )-wave pairing Hamiltonian. The pseudo-deformation of the underlying su (2) algebra is here introduced as a way to obtain these models in the contraction limit of different RG models. This allows for the construction of the full set of conserved charges, the Bethe ansatz state, and the resulting RG equations. For these models an alternative and simpler set of quadratic equations can be found in terms of the eigenvalues of the conserved charges. Furthermore, the recently proposed eigenvalue-based determinant expressions for the overlaps and form factors of local operators are extended to these models, linking the results previously presented for the Dicke–Jaynes–Cummings–Gaudin models with the general results for RG XXZ models.
In this work, we construct an alternative formulation to the traditional algebraic Bethe ansätze for quantum integrable models derived from a generalized rational Gaudin algebra realized in terms of a collection of spins 1/2 coupled to a single bosonic mode. The ensemble of resulting models which we call Dicke–Jaynes–Cummings–Gaudin models are particularly relevant for the description of light–matter interaction in the context of quantum optics. Having two distinct ways to write any eigenstate of these models we then combine them in order to write overlaps and form factors of local operators in terms of partition functions with domain wall boundary conditions. We also demonstrate that they can all be written in terms of determinants of matrices whose entries only depend on the eigenvalues of the conserved charges. Since these eigenvalues obey a much simpler set of quadratic Bethe equations, the resulting expressions could then offer important simplifications for the numerical treatment of these models.
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.
The Dicke model is derived in the contraction limit of a pseudo-deformation of the quasispin algebra in the su (2)-based Richardson-Gaudin models. Likewise, the integrability of the Dicke model is established by constructing the full set of conserved charges, the form of the Bethe Ansatz state, and the associated Richardson-Gaudin equations. Thanks to the formulation in terms of the pseudo-deformation, the connection from the su (2)-based Richardson-Gaudin model towards the Dicke model can be performed adiabatically.
In this work, we construct an alternative formulation to the traditional Algebraic Bethe ansatz for rational Gaudin models realised in terms of a collection of spins 1/2 coupled to a single bosonic mode. In doing so, we obtain two distinct ways to write any eigenstate of these models which can then be combined to write overlaps and form factors of local operators in terms of partition functions with domain wall boundary conditions. We also demonstrate that they can therefore all be written in terms of determinants of matrices whose entries only depend on the eigenvalues of the conserved charges. Since these eigenvalues obey a much simpler set of quadratic Bethe equations, the resulting expressions could then offer important simplifications for the numerical treatment of these models.
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.