“…a (a † ) is the annihilation (creation) operator of the EM field, σ z = |e e| − |g g| and σ + = |e g| (σ − = σ † + ) are atomic operators, with |g and |e denoting the ground and excited atomic states, respectively. Although largely studied over the last decades, up to now its exact analytical solution is lacking and only numerical [2,3,4,5] and approximate analytical solutions are available [6,7,8], despite the conjecture by Reik and Doucha [9] that an exact solution of RH in terms of known functions is possible. The most used analytical approach to RH is to make the rotate wave approximation (RWA), where the antirotating term g a † σ + + aσ − is neglected , since in the weak coupling regime g/ω 1, small detuning |∆| ω (∆ = ω 0 − ω), and weak field amplitude its contribution to the evolution of the system is quite small [10,11].…”