Separation of variables for integrable spin-boson models

We formulate the functional Bethe ansatz for bosonic (infinite dimensional) representations of the Yang-Baxter algebra. The main deviation from the standard approach consists in a half infinite 'Sklyanin lattice' made of the eigenvalues of the operator zeros of the Bethe annihilation operator. By a separation of variables, functional TQ equations are obtained for this half infinite lattice. They provide valuable information about the spectrum of a given Hamiltonian model. We apply this procedure to integrable spin-boson models subject to both twisted and open boundary conditions. In the case of general twisted and certain open boundary conditions polynomial solutions to these TQ equations are found and we compute the spectrum of both the full transfer matrix and its quasi-classical limit. For generic open boundaries we present a two-parameter family of Bethe equations, derived from TQ equations that are compatible with polynomial solutions for Q. A connection of these parameters to the boundary fields is still missing.