Abstract
We present a general qubit-boson interaction Hamiltonian that describes the Jaynes-Cummings model and its extensions as a single Hamiltonian class. Our model includes non-linear processes for both the free qubit and boson field as well as non-linear, multi-boson excitation exchange between them. It shows an underlying algebra with supersymmetric quantum mechanics features allowing an operator based diagonalization that simplifies the calculations of observables. As a practical example, we show the evolution of the population inversion and the boson quadratures for an initial state consisting of the qubit in the ground state interacting with a coherent field for a selection of cases covering the standard Jaynes-Cummings model and some of its extensions including Stark shift, Kerr-like, intensity dependent coupling, multi-boson exchange and algebraic deformations.