Abstract. This article presents a review of recent developments on various aspects of the quantum Rabi model. Particular emphasis is given on the exact analytic solution obtained in terms of confluent Heun functions. The analytic solutions for various generalisations of the quantum Rabi model are also discussed. Results are also reviewed on the level statistics and the dynamics of the quantum Rabi model. The article concludes with an introductory overview of several experimental realisations of the quantum Rabi model. An outlook towards future developments is also given.

We develop a method to find analytical solutions for the eigenstates of the quantum Rabi model. These include symmetric, anti-symmetric and asymmetric analytic solutions given in terms of the confluent Heun functions. Both regular and exceptional solutions are given in a unified form. In addition, the analytic conditions for determining the energy spectrum are obtained. Our results show that conditions proposed by Braak [Phys. Rev. Lett. 107, 100401 (2011)] are a type of sufficiency condition for determining the regular solutions. The well-known Judd isolated exact solutions appear naturally as truncations of the confluent Heun functions.

We introduce a novel concept of the pseudo-parity-time (pseudo-PT) symmetry in periodically modulated optical systems with balanced gain and loss. We demonstrate that whether or not the original system is PT symmetric, we can manipulate the property of the PT symmetry by applying a periodic modulation in such a way that the effective system derived by the high-frequency Floquet method is PT symmetric. If the original system is non-PT symmetric, the PT symmetry in the effective system will lead to quasistationary propagation that can be associated with the pseudo-PT symmetry. Our results provide a promising approach for manipulating the PT symmetry of realistic systems.

We investigate the energy spectrum for hybrid mechanical systems described by non-parity-symmetric quantum Rabi models. A set of analytical solutions in terms of the confluent Heun functions and their analytical energy spectrum is obtained. The analytical energy spectrum includes regular and exceptional parts, which are both confirmed by direct numerical simulation. The regular part is determined by the zeros of the Wronskian for a pair of analytical solutions. The exceptional part is relevant to the isolated exact solutions and its energy eigenvalues are obtained by analyzing the truncation conditions for the confluent Heun functions. By analyzing the energy eigenvalues for exceptional points, we obtain the analytical conditions for the energy-level crossings, which correspond to two-fold energy degeneracy.