---
abstract: 'Rabi oscillations is a key phenomenon among the variety of quantum optical effects that manifests itself in the periodic oscillations of a two-level system between the ground and excited states when interacting with electromagnetic field. Commonly, the rate of these oscillations scales proportionally with the magnitude of the electric field probed by the two-level system. Here, we investigate the interaction of light with a two-level quantum emitter possessing permanent dipole moments. The semi-classical approach to this problem predicts slowing down and even full suppression of Rabi oscillations due to asymmetry in diagonal components of the dipole moment operator of the two-level system. We consider behavior of the system in the fully quantized picture and establish the analytical condition of Rabi oscillations collapse. These results for the first time emphasize the behavior of two-level systems with permanent dipole moments in the few photon regime, and suggest observation of novel quantum optical effects.'
author:
- 'Denis G. Baranov'
- 'Mihail I. Petrov'
- 'Alexander E. Krasnok'
bibliography:
- 'tls.bib'
title: 'Decoupling light and matter: permanent dipole moment induced collapse of Rabi oscillations'
---

Introduction
============

Theory of a two-level system (TLS) interacting with electromagnetic field is of prime importance for a wide spectrum of applied problems, including laser science [@Sargent; @Scully], fluorescent spectroscopy [@Lakowicz2008], nano-imaging [@SandoghdarN2000; @Acuna2012; @Beams2013], design of single photon sources [@Arcari2014; @Lodahl2015; @Rao2007; @Hoang2016] and efficient light emitting devices [@Lodahl2015; @LEDs]. It also plays the central role in the quantum information theory in the context of coherent qubits control [@Chiorescu2003; @Neeley2008; @Wallraff2004; @Lupascu2007]. Due to its very general formalism, it is equally applicable to the wide range of different electronic systems: electric and magnetic atomic and molecular transitions, quantum dots and quantum wells, superconducting Josephson junctions, defect centers in nanocrystals and others.

In the theoretical description of light-matter interaction, it is often assumed that dynamics of the system is governed by the non-diagonal matrix element of the dipole moment operator. However, certain systems possess non-zero permanent dipole moments (PDM). An example of such system is a polar molecule, an atom polarized by static electric field [@Potassium] or an asymmetric quantum dot [@Savenko]. Magnetic dipole atomic transitions, e.g., in rare-earth ions of Eu$^{3+}$ [@Zia2012; @Carminati2014; @Novotny2015] may also have non-zero permanent magnetic dipoles, in contrast to the case of electric dipole transitions, where diagonal elements of the electric dipole moment operator are *always zero* for atomic eigenstates [@Landau], what follows from the parity of the wave function.

Quantum systems with electric PDM have been widely investigated in the context of multi-photon processes [@Meath84; @Meath85; @Band90; @Meath92; @AIP]. It was shown that presence of permanent dipole moment modifies multi-photon absorption rates. Emission spectrum features of quantum systems possessing permanent dipole moments were studied in Refs. [@Hoffmann; @Savenko], where it was shown that such a system can radiate at Rabi frequency and serve as an emitter in the THz range. More recently, it was demonstrated that a two-level system with permanent dipoles can be inversely populated in the steady state if acted by two monochromatic fields [@Macovei].

![\[fig1\]A schematic of the system under study. A two-level system modeling a quantum emitter interacts with electromagnetic field of a cavity resulting in Rabi oscillations of the two-level system population inversion. Inset on the right shows the internal structure of the two-level system with permanent dipole moments.](concept.pdf){width="45.00000%"}

Here, we investigate in more details the Rabi oscillations of a TLS with PDM. One of the main effects in such systems is that PDM enables multi-photon Rabi oscillations and strongly affects the Rabi frequency for single photon process. Since the pioneering paper of Meth et. al. [@Meath84], it was shown that the Rabi oscillations in principle may be suppressed by tuning the magnitude of the emitter PDM. However, the treatment of this problem has been mostly limited to the semiclassical approach. In this work, we address the fully quantum picture of interaction of light with a permanent dipole TLS. By means of numerical simulations and theoretical analysis, we find new eigenstates of the PDM TLS–single mode cavity Hamiltonian and observe the collapse of Rabi oscillations in the few photon regime. Our findings demonstrate previously unexplored regime of light-matter interaction and provide novel tools for coherent control of quantum emitters possessing PDM.

Semiclassical description
=========================

We start from a brief overview of the semiclassical description of the dynamics of a TLS with PDM in the presence of light field. The system under study is schematically depicted in Fig. 1. It consists of a generalized TLS interacting with electromagnetic field. The TLS can make transition between the ground $\left| g \right\rangle $ and excited $\left| e \right\rangle $ states, respectively, separated by the energy $\hbar \omega_0$.

To begin with, we consider the simplest scenario of a TLS driven by a classical monochromatic electromagnetic wave ${\bf{E}}\left( t \right) = {{\bf{E}}_0}\cos \left( {{\omega _0}t} \right)$. In what follows we will use the notation ${\bf{E}}\left( t \right)$ for description of the field component interacting with the TLS, implying the electric field. However, this formalism is equally applicable to any other TLS interacting with either electric or magnetic oscillating field and therefore this does not limit the ubiquity of our results. Despite simplicity of this model, it captures important signatures in the dynamics of a TLS with PDM. The Hamiltonian of the system has the standard form $\hat H = {\hat H_{{\rm{TLS}}}} + {\hat H_{{\mathop{\rm int}} }}$, where the TLS part is $\hat H = \hbar {\omega _0}{{\hat \sigma }^\dag }\hat \sigma$ and the interaction part is ${\hat H_{{\mathop{\rm int}} }} =  - {\bf{E}}\left( t \right){\bf{\hat d}}$. The dipole moment operator ${\bf{\hat d}}$ has the form $${\bf{\hat d}} = {{\bf{d}}_{ge}}\hat \sigma  + {\bf{d}}_{eg}{\hat \sigma ^\dag } + {{\bf{d}}_{gg}}\hat \sigma {\hat \sigma ^\dag } + {{\bf{d}}_{ee}}{\hat \sigma ^\dag }\hat \sigma,
\label{eq1}$$ Here $\hat \sigma = \left| g \right\rangle \left\langle e \right|$ is the lowering operator, ${{\bf{d}}_{eg}}={\bf{d}}_{ge}^{*}$ is the transition dipole moment, and ${{\bf{d}}_{ee}}$ and ${{\bf{d}}_{gg}}$ are the permanent dipole moments which are equal to zero in the common Jaynes-Cummings model [@Scully]. Without loss of generality we will assume the that ${\bf{d}}_{eg}$ along with ${{\bf{d}}_{ee}}$ and ${\bf d}_{gg}$ are all real-valued.

With the use of expression (\[eq1\]) for the dipole moment operator we write the Hamiltonian in the symmetric form: $$\begin{gathered}
\hat H =  - \frac{{\hbar {\omega _0}}}{2}\hat \sigma {{\hat \sigma }^\dag } + \frac{{\hbar {\omega _0}}}{2}{{\hat \sigma }^\dag }\hat \sigma  - \hbar {\Omega _R}\cos \left( {{\omega _0}t} \right){{\hat \sigma }^\dag }\hat \sigma \\
 - {{\bf{d}}_{gg}}{{\bf{E}}_0}\cos \left( {{\omega _0}t} \right)\hat \sigma {{\hat \sigma }^\dag } - {{\bf{d}}_{ee}}{{\bf{E}}_0}\cos \left( {{\omega _0}t} \right){{\hat \sigma }^\dag }\hat \sigma 
\end{gathered}$$ where $\hbar {\Omega _R} = {{\bf{E}}_0}{{\bf{d}}_{eg}}$ is the Rabi frequency. In the absence of permanent dipoles, the population inversion of the TLS oscillates in time with frequency $\Omega_R$.

![(a, b) Population inversion of a TLS with non-zero permanent dipoles (solid) and an equivalent TLS without permanent dipoles (dashed) in a monochromatic electric field at different values of $\Omega_p$. Significant delay becomes visible in the permanent dipoles TLS dynamics. (c) Modified Rabi frequency of the TLS with permanent dipole moments for different $\Omega_p/\Omega_R$ as a function of applied electric field $E_0$. Dashed line shows unchanged Rabi frequency of a TLS without permanent dipoles. Inset: temporal dynamics of the population inversion at the collapse point corresponding to the first zero of the first order Bessel function $J_1$, $\Omega_{p}\approx 3.8\omega_0$.[]{data-label="fig2"}](2.pdf){width="1\columnwidth"}

Representing the TLS wave function in the general form $\left| {\Psi \left( t \right)} \right\rangle  = {c_e}\left( t \right)\left| e \right\rangle  + {c_g}\left( t \right)\left| g \right\rangle $ and substituting into the Schroedinger equation, we arrive at the following equations of motion describing temporal evolution of the TLS in the interaction representation: $$\begin{gathered}
 {{\dot c}_g}(t) = i \Omega_{R}\cos(\omega t)\exp\left(-i\left(\omega_0 t - \dfrac{\Omega_{p}}{\omega} \sin(\omega t)\right)\right)c_{e}(t), \\
 {{\dot c}_e}(t) = i \Omega_{R}\cos(\omega t)\exp\left(i\left(\omega_0 t - \dfrac{\Omega_{p}}{\omega} \sin(\omega t)\right)\right)c_{g}(t).
\end{gathered}
\label{eq3}$$ Here we have introduced the quantity $\hbar {\Omega _{{\text{p}}}} = {{\bf{E}}_0}\left( {{{\bf{d}}_{ee}} - {{\bf{d}}_{gg}}} \right)$, which will be referred to as the *permanent* coupling constant. It shows the strength of asymmetry between the dipole moment elements in ground and excited states.

Figures \[fig2\](a),(b) show the TLS population inversion $W(t)={\left| {c_e\left( t \right)} \right|^2} - {\left| {c_g\left( t \right)} \right|^2}$ as a function of time obtained from the numerical solution of Eq. \[eq3\] for different ratios $\Omega_{p}/\omega_0$. For comparison we also plot the population inversion of a pure non-diagonal system with ${{{\bf{d}}_{ee}} - {{\bf{d}}_{gg}}}=0$. The Rabi frequency is fixed at value $\Omega_R=0.01\omega_0$. The TLS is assumed to be initially in the ground state $\left| g \right\rangle $. It is seen that, in the regime of weak permanent coupling, ${\Omega _{p}} \ll {\omega_0}$, the TLS dynamics is unaffected by permanent dipoles \[Fig. 2(a)\]. However, in the strong permanent coupling regime, ${\Omega _{p}} \sim \omega_0$, the effect of permanent dipoles on the dynamics of the TLS becomes clearly visible, Fig. 2(b). Specifically, this effect consists in a larger period of Rabi oscillations cycle.

The slowing of Rabi oscillations observed in the numerical simulations can be described by an elegant analytical expression derived in Refs. [@Meath2000; @Raikh; @Hoffmann]. The resulting modified Rabi frequency is given by $${\tilde \Omega _R} = 2{\omega _0}\frac{{{\Omega _R}}}{{{\Omega _{\text{p}}}}}{J_1}\left( {\frac{{{\Omega _{\text{p}}}}}{{{\omega _0}}}} \right),
\label{eq5}$$ which yields $\Omega_R$ in the limit ${\Omega _{p}} \to 0$. The analytical estimation for the modified Rabi frequency $\tilde \Omega_R$ for emitters with different ratios of $\Omega_{\text{p}}/\Omega_R$ is presented in Fig. 2(c) as a function of applied oscillating electric field $E_0$. It clearly demonstrates the nonlinear dependence of the modified Rabi frequency shift with increasing $E_0$.

It immediately follows from Eq. (\[eq5\]) that the strength of the PDM effect is dictated by the relation between between the permanent coupling and the transition frequency. It also shows that for observation of the predicted modification of inversion oscillations one does not need to achieve the ultrastrong coupling characterized by $\Omega_R \sim \omega_0$ [@Sirtori2010USC; @Liberato2014]. Instead, it is sufficient to achieve only permanent ultrastrong coupling with $\Omega_{\text{p}} \sim \omega_0$ which is easier for quantum systems with large permanent dipoles $\left| {{{\bf{d}}_{ee}} - {{\bf{d}}_{gg}}} \right| \gg \left| {{{\bf{d}}_{eg}}} \right|$.

Equation \[eq5\] also suggests that the modified Rabi frequency may reach zero what would indicate *collapse of Rabi oscillations*. The value of incident field $E_{\text{coll}}$ at which Rabi oscillations collapse occurs can be simply estimated by equating the argument of Bessel function to its first zero: $${E_{{\text{coll}}}} = {x_1}\frac{{\hbar {\omega _0}}}{{\left| {{{\mathbf{d}}_{ee}} - {{\mathbf{d}}_{gg}}} \right|}}
\label{eq6}$$ with $x_1 \approx 3.8$ being the first zero of the $J_1$ Bessel function. The point of Rabi oscillations collapse can be found in Fig. \[fig2\](c) for the TLS with $\Omega_{\text{p}}/\Omega_R=10$. The corresponding temporal dynamics of the population inversion at the point of Rabi oscillations collapse $\hbar {\Omega _{{{p}}}} = {E_{{\rm{coll}}}}\left| {{{\bf{d}}_{ee}} - {{\bf{d}}_{gg}}} \right|$ is shown in the inset in Fig. \[fig2\](c). Surprisingly, inversion of the TLS experiences only small deviations from its initial value $-1$ despite the considerable classical Rabi frequency $\Omega_R/\omega_0 \approx 0.38$. At this point we observe collapse of Rabi oscillations. With further increase of $\Omega_R$ and $\Omega_{p}$ Rabi oscillations recover, in consistence with expression (\[eq5\]).

The additional term in Hamiltonian (2) containing $\Omega_{p}$ describes fast modulation of the transition frequency $\omega_0$. From this point of view, the system can be treated as a parametric oscillator. Given that the external electromagnetic field is tuned to the TLS transition frequency $\omega_0$, this modulation induced by the permanent dipoles introduces quickly variable detuning between the TLS and external field at certain time intervals. A similar Rabi oscillations slowing has been reported recently for a spin-1/2 system in a slowly modulated longitudinal magnetic field [@Raikh].

The semiclassical approach presented above is commonly used for description of TLS dynamics with PDM as it draws a simple yet very illustrative picture. However, it does not allow one to understand the interaction of such TLS with single photons. To the best of our knowledge, only two works have treated the problem of PDM TLS interaction with light within the fully quantum approach [@Macovei; @Savenko]. Nevertheless, the intriguing effect of Rabi oscillations collapse in the quantum regime has not been addressed.

![image](Perm_RWA.pdf){width=".9\textwidth"}

Quantum description
===================

We start our consideration with the standard fully quantized Hamiltonian of a TLS interacting with a single electromagnetic mode: $$\hat H = \hbar {\omega _0}{\hat \sigma ^\dag }\hat \sigma  + \hbar \omega {\hat a^\dag }\hat a - {\mathbf{\hat d\hat E}}\left( {{{\mathbf{r}}_d}} \right)
\label{eq7}$$ where $\omega$ is the resonator frequency, $\hat a$ is the annihilation operator of the resonator mode, ${\mathbf{\hat d}}$ is the dipole moment operator given by Eq. (1), ${\bf{\hat E}}\left( {{{\bf{r}}_d}} \right) = {\cal E}_{0}\left( {{{\hat a}^\dag } + \hat a} \right)$ is the electric field operator at the TLS position $\mathbf{r}_{d}$, and ${\cal E }_{0}=\sqrt{2\pi\hbar \omega_0/V}$ is the vacuum field amplitude with $V$ being the mode volume.

In order to draw an analytical picture of the interaction, we find the energy spectrum of the full Hamiltonian (\[eq7\]). Following the semiclassical case, where the system is written in the interaction picture, we rewrite the Hamiltonian using the unitary transformation: $${\widehat{H}}'={\widehat{D}}(z){\widehat{H}}{\widehat{D}}^{-1}(z), \qquad z=\dfrac{iz_0}{{\omega}}e^{-\frac{i{\omega}t}2}\sin({\omega}t),$$ where $\hat{D}(z)=\exp(z {\widehat{a}}^{+} - z^{*}{\widehat{a}} )$ is the displacement operator [@Scully], and $z_0=-\mathcal{E}_{0}({\bf d_{ee}} + {\bf d_{gg}})/\hbar$. The new Hamiltonian $\hat{H}'$ can be separated into the following three components: $$\begin{gathered}
{\widehat{H}}'={\widehat{H}}_{0}+ {\widehat{H}}_{JC}+{\widehat{H}}_{p},\\
{\widehat{H}}_{0}=\dfrac{\hbar{\omega}_{0}}2{{\widehat{\sigma}}}_{z}+\hbar {\omega}{\widehat{b}}^{+}{\widehat{b}},\\ 
{\widehat{H}}_{JC}=\hbar g({\widehat{b}}^{+}+{\widehat{b}})({{\widehat{\sigma}}}_{+}+{{\widehat{\sigma}}}_{-}),\\
{\widehat{H}}_{p}=\dfrac{\hbar g_{p}}{2}({\widehat{b}}^{+}+{\widehat{b}}){{\widehat{\sigma}}}_{z},\\
\end{gathered}
\label{eq9}$$ where ${{\widehat{\sigma}}}_z={{\widehat{\sigma}}}_+{{\widehat{\sigma}}}_--{{\widehat{\sigma}}}_-{{\widehat{\sigma}}}_+$, $g=-\mathcal{E}_{0}{\bf d_{eg}}/\hbar$ is the Jaynes-Cummings (JC) coupling constant, $g_{\rm p}=-\mathcal{E}_{0}({\bf d_{gg}} - {\bf d_{ee}})/\hbar$ is the coupling constant related to the PDM contribution, and ${\widehat{b}}={{\widehat{a}}}-z$ is the new cavity annihilation operator. In Eq. (\[eq9\]) the component ${\widehat{H}}_{0}$ corresponds to the energy of non-interacting free atom and field, the second ${\widehat{H}}_{JC}$ corresponds to the JC part, and the last component ${\widehat{H}}_p$ is the energy operator related to the permanent dipole contribution. We will also use the notation ${\widehat{H}}_{JC}^{RWA}=\hbar g({\widehat{b}}^{+}{{\widehat{\sigma}}}_{-}+{{\widehat{\sigma}}}_{+}{\widehat{b}})$ for the JC term written in the rotating wave approximation (RWA).

Now we have two terms in the full Hamiltonian ${\widehat{H}}'$ that modify the interaction-free atom-field system: the JC part ${\widehat{H}}_{JC}$ and the PDM part ${\widehat{H}}_{\rm p}$. This interaction leads to Rabi splitting, which manifests itself in Rabi oscillations of the TLS population inversion. The frequency of these oscillations is directly related to the magnitude of this splitting. Thus, the collapse of the Rabi oscillations in the quantum picture is accompanied by zero Rabi splitting, i.e., degeneracy of the eigen states of the system.

A natural way to approach the eigen states problem seems to treat the PDM part of the Hamiltonian as a perturbation over the classical JC Hamiltonian ${\widehat{H}}_{0}+{\widehat{H}}_{JC}$, which solution can be obtained either in RWA approximation \[35\] or beyond that \[36\]. However, the effect of Rabi oscillations collapse is expected to appear in the regime of strong permanent coupling, thus, such perturbative approach would not be valid. We will start with considering the JC part as a perturbation over the permanently polarized atomic system described by the Hamiltonian ${\widehat{H}}_{0-p}={\widehat{H}}_{0}+{\widehat{H}}_{p}$. The case of diagonal part as a perturbation will be also considered for illustrative purposes in the end of this section. In the following, to get a simple and clear picture of the effect of permanent contribution, we will focus on the resonant case, which implies $\omega=\omega_{0}$.

Strong permanent coupling: the collapse of Rabi splitting
---------------------------------------------------------

We start with finding the eigen states of ${\widehat{H}}_{0-p}$ Hamiltonian, which is already diagonal in terms of the atomic states as contains the ${\widehat{\sigma}}_{z}$ atomic operator only. To diagonalize the photonic part we introduce the generalized (or extended) [@DeOliveira1990; @Philbin2014] coherent states following the recent works on exact solvability of the Rabi model reported in Ref. [@Braak; @Chen2012]. Using simple Bogoluybov’s transformation ${\widehat{A}}={\widehat{b}}+g_{p}/2$ and ${\widehat{B}}={\widehat{b}}-g_{p}/2$ one can get the diagonal form of the Hamiltonian $$\begin{gathered}
{\widehat{H}}_{0-p}^{'}=\dfrac{\hbar {\omega}_{0}}2 +\hbar {\omega}_{0}\left({\widehat{A}}^{+}{\widehat{A}} -\dfrac{g^{2}}{4} \right)|e \rangle \langle e|\\ 
+\hbar {\omega}_{0}\left({\widehat{B}}^{+}{\widehat{B}} -\dfrac{g^{2}}{4} \right)|g \rangle \langle g|,
\end{gathered}
\label{H0p}$$ with the eigen functions $$\label{eigH0p}
|\psi_{n,A}\rangle=|n\rangle_{A}|e\rangle \qquad \mbox{and }\quad |\psi_{n,B}\rangle=|n\rangle_{B}|g\rangle.$$ Here $|n\rangle_{A}$ and $|n\rangle_{B}$ are the eigen states (number of photon states) for ${\widehat{A}}^{+}{\widehat{A}}$ and ${\widehat{B}}^{+}{\widehat{B}}$ operators respectively, and can be considered also as generalized coherent states constructed from the displaced Fock’s states: $$|n\rangle_{A}={\widehat{D}}({g_{p}}/2)|n\rangle \qquad |n\rangle_{B}={\widehat{D}}(-{g_{p}}/2)|n\rangle.$$ One can see from the form of the Hamiltonian that the energy levels corresponding to the $|\psi_{n-1,A}\rangle $ and $|\psi_{n,B}\rangle $ states are doubly degenerate and equal to $E_{n}=\hbar {\omega}_{0}n+g_{p}^{2}/4$, and the energy is quadratic with respect to the permanent coupling constant. Using the degenerate perturbation method in the basis we can find the energy splitting between the levels, which is defined by the matrix elements $\langle\psi_{n-1,A}|{\widehat{H}}_{JC}^{RWA}|\psi_{n,B}\rangle$. This Rabi splitting is related to the displacement operator matrix elements in the Fock’s basis: $$\begin{gathered}
\tilde{\Omega}_{R, n}=2\langle\psi_{n-1,A}|{\widehat{H}}_{JC}^{RWA}|\psi_{n,B}\rangle=\\ 
g\left(\sqrt{n+1}\langle n+1|{\widehat{D}}(-g_{p}/\omega_{0})|n+1\rangle-\right.\\ 
\left.\dfrac{g_{p}}{2\omega_{0}}\langle n+1|{\widehat{D}}(-g_{p}/\omega_{0})|n\rangle\right).\end{gathered}$$

The matrix elements of the displacement operator can be directly computed [@Tanas1992] that gives the final expression containing Laguerre polynomial $L_{n}$: $$\begin{gathered}
\label{QuantumSplit}
\tilde{\Omega}_{R, n}=
\dfrac{\Omega_{R, n}}{2}\exp\left(-\dfrac{g_{p}^{2}}{2\omega_{0}^{2}}\right)\left[L_{n+1}\left(\dfrac{g_{p}^{2}}{\omega_{0}^{2}}\right)\right. \\
\left.+L_{n}\left(\dfrac{g_{p}^{2}}{\omega_{0}^{2}}\right)\right], 
\end{gathered}$$ where $\Omega_{R, n}=2g\sqrt{n+1}$ is Rabi splitting in the JC model.

The dependence of Rabi splitting on the permanent coupling constant $g_{p}$ is illustrated in Fig. \[Permanent\_JC\](a). The calculations were performed with the JC part written in the RWA. The nonzero value of the coupling constant $g$ gives finite splitting of the eigenstates at $g_{p}=0$. Increase of the permanent coupling constant $g_p$ results in oscillatory behavior of the energy levels, clearly seen in Fig. \[Permanent\_JC\](a). At the intersection of two spectral branches the Rabi splitting vanishes, which is the trace of the Rabi oscillations collapse. This occurs for a discrete series of values of $g_{p}$, which are determined as the roots of the expression . This expression is the analogue of the formula for the fully quantized case. According to the properties of the Laguerre polynomials [@Arfken1972], the number of Rabi collapse points is the same as the order of the energy band (see Fig. \[Permanent\_JC\]).

The permanent dipole term ${\widehat{H}}_{p}$ is not energy conserving since it contains the terms ${{\widehat{a}}}{{\widehat{\sigma}}}_{-}$ and ${{\widehat{a}}^{+}}{{\widehat{\sigma}}}_{+}$, which rotate at high frequency. Thus, it is not fully consistent to consider the permanent dipole Hamiltonian without keeping the non-RWA terms in the JC term. This is the well-known limitation of the RWA model [@Irish; @Braak]. In order to investigate the energy spectrum of the system beyond the RWA, we calculate the eigenvalues of the Hamiltonian (\[eq9\]) numerically, Fig. \[Permanent\_JC\](b). The non-RWA terms lead to the interaction of the neighboring energy terms and the degeneracy in the points of Rabi collapse in the RWA model \[see dashed line in Fig. \[Permanent\_JC\] (b)\] is lifted and non-zero splitting is observed.

![Simulated population inversion dynamics of a TLS with non-zero permanent dipoles (solid) and an equivalent TLS without permanent dipoles (dashed) interacting with the three-photon Fock state of the cavity for different ratios $g_{p}/\omega_0$. Panel (b) corresponds to the intersection of the two spectral branches with vanishing Rabi splitting leading to collapse of the Rabi oscillations.[]{data-label="fig5"}](5.pdf){width="1\columnwidth"}

The predicted collapse of Rabi oscillations in the quantum regime is directly illustrated via numerical modeling of temporal dynamics of the system described by Hamiltonian . Dynamical behavior of the system initially placed in the Fock state $\left| {{\Psi _0}} \right\rangle  = \left| {g,n} \right\rangle$ is presented in Fig. \[fig5\]. The Schroedinger equations of motion in the non-RWA regime are solved for different ratios $g_p/\omega_0$ and the fixed cavity photons number $n=3$. The JC coupling constant is fixed at the value $g=0.15 \omega_0$ being consistent with Fig. \[Permanent\_JC\](b). Similarly to the results of the semiclassical approach, frequency of the Rabi oscillations decreases with increasing $g_p$, Fig. \[fig5\](a). However, in contrast to the semiclassical picture, not only the oscillations period is affected, but the amplitude of Rabi oscillations also changes \[see Fig. \[fig2\](b) for comparison\]. Finally, at the intersection point of the two JC branches, shown in the inset of Fig. \[Permanent\_JC\] (b) with circle, we observe collapse of Rabi oscillations, Fig. \[fig5\](b).

![image](JC_Perm.pdf){width=".9\textwidth"}

Strong non-diagonal coupling.
-----------------------------

In order to accomplish the picture of the eigen states of the system under study, we also consider the case of strong non-diagonal coupling in more details, assuming that the diagonal part can be treated as a perturbation. Similar to the previous section, we begin with the JC Hamiltonian in the simplified RWA form, where the terms ${{\widehat{a}}}{{\widehat{\sigma}}}_{-}$ and ${{\widehat{a}}^{+}}{{\widehat{\sigma}}}_{+}$. The Eigenstates of the non-perturbed Hamiltonian $H_{JC}^{RWA}$ are the well-known dressed light-matter states [@Scully], which in the case of resonant excitation ${\omega}={\omega}_{0}$ have simple form [@Meystre]: $$\begin{gathered}
  \left| {\psi _n^ + } \right\rangle  = \frac{1}
{{\sqrt 2 }}\left| {g,n+1} \right\rangle  + \frac{1}
{{\sqrt 2 }}\left| {e,n } \right\rangle , \hfill \\
  \left| {\psi _{n }^ - } \right\rangle  = \frac{1}
{{\sqrt 2 }}\left| {e,n } \right\rangle  - \frac{1}
{{\sqrt 2 }}\left| {g,n+1} \right\rangle , \hfill \\ 
\end{gathered}$$ with the corresponding energy spectrum of the excited states given by $$\begin{gathered}
  E_n^ \pm  = \left(n+\dfrac12\right)\hbar {\omega _0} \pm g\sqrt {n+1} ,n \ge 0 \hfill \\
\end{gathered}$$ The ground state has energy $ -{\hbar {\omega _0}}/2 $. The energy spectrum of the non-perturbed JC part is shown in Fig. \[JC\_Permanent\](a) and (b) with dashed red lines. Each line corresponds to a state with a fixed number of photons $n$. Without interaction, $g=0$, the states with the same number of quanta are degenerate, but the energy splitting between them $E_{n}^{+}-E_{n}^{-}$ increases linearly with increase of $g$ (red dashed line). Thus, in the JC-RWA model only neighboring states $\left| {\psi _n^ + } \right\rangle $ and $\left| {\psi _n^ - } \right\rangle $ interact. However, after introducing the permanent component ${\widehat{H}}_{p}$, which contains terms ${\widehat{b}}^{+}$ and ${\widehat{b}}^{-}$ changing the number of photons by 1, the neighboring polaritonic states $\left| {\psi _{n-1}^ + } \right\rangle $ and $\left| {\psi _{n}^ - } \right\rangle $ start to interact resulting in an anticrossing behavior (blue solid lines in Fig. \[JC\_Permanent\](a)). The magnitude of the energy gap equals to approximately $(\sqrt{n+1}+\sqrt{n+2})g_{p}/2$, which can be found from the perturbation theory by taking the $\left| {\psi _{n-1}^ + } \right\rangle $ and $\left| {\psi _{n}^ - } \right\rangle $ as non-perturbated solutions.

![ The two-dimensional plots of the eigen energies as the function of Permanent coupling constant $g_{p}$ and JC coupling constant $g$ with JC part written in the RWA (a) and non-RWA (b) from. The four polaritonic branches with $n=2$ and $n=3$ are shown. []{data-label="Permanent_JC_2D"}](2D_energy.pdf){width="1\columnwidth"}

Although the RWA allows us obtaining simple analytical results in the vicinity of energy terms crossings, the permanent Hamiltonian component contains non-RWA terms. In order to stay at the same level of generality we have numerically computed the energy of the eigen states considering JC counterpart in non-RWA approximation, and the results are shown in the Fig. \[JC\_Permanent\] (b). Similar to the RWA case the gap opening for neighboring states still exists, but due to the higher order terms ${\widehat{b}}^{+} {\widehat{\sigma}}_{+}$ and ${\widehat{b}} {\widehat{\sigma}}_{-}$ in JC Hamiltonian, the states differing in number of quanta $n$ by two also start to interact, which was reported previously in Ref. [@Irish].

Finally, it is illustrative to summarize the analyzed cases in a 2D plot showing the eigen energy levels of the full Hamiltonian as a function of two coupling constants $g$ and $g_{p}$. This plot is show in Fig. \[Permanent\_JC\_2D\] for two cases of Jaynes-Cummings part written in RWA and non-RWA form. We have discussed in more details the cases of weak permanent and weak non-diagonal coupling, which correspond to cross section of the Fig. \[Permanent\_JC\_2D\] along $x$ and $y$ axis. One can also see the ripple-like behavior of the 2D surface in RWA approximation, which indicates that the collapse of Rabi splitting occurs not only for weak coupling regimes. However, this structure of the eigen states is smeared out in the non-RWA regime, shown in Fig. \[Permanent\_JC\_2D\] (b).

Outlook and Conslusion
======================

So far we have treated interaction of a TLS with a closed lossless cavity. At the same time, the Rabi oscillations slowing demonstrated here may have important effect on spontaneous decay of a TLS with permanent dipole. Indeed, in the weak coupling regime the spontaneous decay rate of a quantum emitter is given by $\Gamma  = 4\Omega _R^2/{\gamma _a}$ where $\gamma_a$ is the cavity mode decay rate [@Scully]. In the particular case of a TLS with permanent dipoles Rabi frequency gets smaller with increasing $\Omega_p$. Thus we may expect that decay rate $\Gamma$ will also decrease for a TLS with permanent dipoles in comparison to a pure non-diagonal quantum emitter.

Another consequence of the predicted modification of Rabi oscillations may be found in the resonance florescence spectra of quantum emitters with permanent dipoles. In the simplest case, it is represented by the Mollow triplet [@Scully] which consists of the central transition line and the two sidebands with frequencies $\omega_0 \pm \Omega_R$. Therefore, modification of $\Omega_R$ will shift the sidebands of the Mollow triplet.

To conclude, we have investigated the energy spectrum and temporal dynamics of a two-level quantum emitter with permanent dipole moments interacting with electromagnetic field. The presence of a permanent dipole moment has a crucial role on the Rabi oscillations of the quantum emitter. When the two-level system is driven by a classical field, Rabi oscillations frequency gradually decreases until oscillations completely collapse. We established the condition of Rabi oscillations collapse in semiclassical and quantum regime and discussed the implications of this phenomenon on other quantum optical effects.

Authors acknowledge fruitful discussion with Prof. O. V. Kibis. D.G.B. acknowledges support from RFBR project No 16-32-00444. M.I.P acknowledges support from RFBR project No 16-32-60167.
