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.
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.
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.
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.
It was recently shown that, for central spin-1/2 and central spin-1, the XX central spin model is integrable in the presence of a magnetic field oriented perpendicular to the XY plane in which the coupling exists. In the spin-1/2 case, it was also shown, through an appropriate limit of the non-skew symmetric XXZ Richardson-Gaudin models, that it remained integrable even when the magnetic field is tilted to contain an in-plane component. Although the model has not yet been shown to explicitly belong to a known class of Richardson-Gaudin models, we show, in this work, that the spin-1 case also remains integrable in a tilted magnetic field. We do so by writing explicitly the complete set of conserved charges, then showing that these operators obey polynomial relations. It is finally demonstrated numerically that dark states, for which the central spin is completely unentangled with the bath, can emerge at strong enough coupling just as they do in the central spin-1/2 model in an arbitrarily oriented magnetic field.
Central spin models describe a variety of quantum systems in which a spin-1 2 qubit interacts with a bath of surrounding spins, as realized in quantum dots and defect centers in diamond.We show that the fully anisotropic central spin Hamiltonian with (XX) Heisenberg interactions is integrable.Building on the class of integrable Richardson-Gaudin models, we derive an extensive set of conserved quantities and obtain the exact eigenstates using the Bethe ansatz.These states divide into two exponentially large classes: bright states, where the qubit is entangled with the bath, and dark states, where it is not.We discuss how dark states limit qubit-assisted spin bath polarization and provide a robust long-lived quantum memory for qubit states.
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.
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.