Algebraic Bethe Ansatze and eigenvalue-based determinants for Dicke-Jaynes-Cummings-Gaudin quantum integrable models

In this work, we construct an alternative formulation to the traditional Algebraic Bethe ansatz for quantum integrable models derived from a generalised rational Gaudin algebra realised 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.

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

We present a numerical approach which allows the solving of Bethe equations whose solutions define the eigenstates of Gaudin models. By focusing on a new 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 Richardsons fermionic pairing model, the central spin model and generalized Dicke model.