The two-photon quantum Rabi model is analyzed using the Z 4 -symmetry and su ( 1 , 1 ) -algebra. We derive the G -function whose zeros give the exact eigenvalues of the Hamiltonian. The singularity structure of the G -function allows to draw conclusions about the distribution of these eigenvalues along the real axis and we derive the spectral collapse phenomenon at critical coupling g c , found numerically before. The spectrum at g c consists of a discrete and a continuous part: the ground state is always separated from the continuum by a finite excitation gap, ruling out a quantum phase transition in the usual sense. For large qubit splitting, also other low lying states split off from the continuum. However, perturbation theory predicts the vanishing of the gap to all orders, demonstrating its non-perturbative nature. We corroborate this result with a variational calculation for the ground state.

The quantum Rabi–Stark model, where the linear dipole coupling and the nonlinear Stark-like coupling are present on an equal footing, are studied within the Bogoliubov operators approach. Transcendental functions responsible for the exact solutions are derived in a compact way, much simpler than previous ones obtained in the Bargmann representation. The zeros of transcendental functions reproduce completely the regular spectra. In terms of the explicit pole structure of these functions, two kinds of exceptional eigenvalues are obtained and distinguished in a transparent manner. Very interestingly, a first-order quantum phase transition indicated by level crossing of the ground state and the first excited state is induced by the positive nonlinear Stark-like coupling, which is however absent in any previous isotropic quantum Rabi models. When the absolute value of the nonlinear coupling strength is equal to twice the cavity frequency, this model can be reduced to an effective quantum harmonic oscillator, and solutions are then obtained analytically. The spectra collapse phenomenon is observed at a critical coupling, while below this critical coupling, infinite discrete spectra accumulate into a finite energy from below.

Dynamics of the sub-Ohmic spin-boson model is examined using three numerical approaches, namely the Dirac-Frenkel time-dependent variation with the Davydov D${}_{1}$ ansatz, the adaptive time-dependent density matrix renormalization group method within the representation of orthogonal polynomials, and a perturbative approach based on a unitary transformation. In order to probe the validity regimes of the three approaches, we study the dynamics of a qubit coupled to a bosonic bath with and without a local field. Comparison of the up-state population evolution shows that the three approaches are in agreement in the weak-coupling regime but exhibit marked differences when the coupling strength is large. The Davydov D${}_{1}$ ansatz and the time-dependent density matrix renormalization group can both be reliably employed in the weak-coupling regime, while the former is also valid in the strong-coupling regime as judged by how faithfully the trial state follows the Schr\"odinger equation. We further explore the bipartite entanglement dynamics between two qubits coupled with individual bosonic baths which reveals entanglement sudden death and revival.

The Dirac-Frenkel time-dependent variation is employed to probe the dynamics of the zero temperature sub-Ohmic spin-boson model with strong friction utilizing the Davydov D1 ansatz. It is shown that initial conditions of the phonon bath have considerable influence on the dynamics. Counterintuitively, even in the very strong coupling regime, quantum coherence features still manage to survive under the polarized bath initial condition, while such features are absent under the factorized bath initial condition. In addition, a coherent-incoherent transition is found at a critical coupling strength α ≈ 0.1 for s = 0.25 under the factorized bath initial condition. We quantify how faithfully our ansatz follows the Schrödinger equation, finding that the time-dependent variational approach is robust for strong dissipation and deep sub-Ohmic baths (s ≪ 1).

A generalized quantum Rabi Hamiltonian with both one- and two-photon terms has emerged in the circuit quantum electrodynamics system for a decade. The usual parity symmetry is broken naturally in the simultaneous presence of both couplings, which complicates analytical treatments, even in the rotating-wave approximations. In this paper, we propose an adiabatic approximation to this generic model by using Bogoliubov operators and obtain a very concise analytical solution for both eigenvalues and eigenstates. Although the adiabatic approximation is only exact in the vanishing limit of the qubit frequency, the results for some physical observables nevertheless agree well with the numerical ones in a wide parameter regime. In the rotating-wave approximations, we also derive an analytical eigensolution. Two dominant Rabi frequencies are found in the Rabi oscillations of this generalized model. We also apply the present analytical theory to the vacuum Rabi splitting. It is found that some different phenomena emerge just because of the presence of the additional two-photon coupling term.