Non-Markovian decoherence of a two-level system in a Lorentzian bosonic reservoir and a stochastic environment with finite correlation time

In this paper we investigate non-Markovian evolution of a two-level system (qubit) in a bosonic bath under influence of an external classical fluctuating environment. The interaction with the bath has the Lorentzian spectral density, and the fluctuating environment (stochastic field) is represented by a set of Ornstein-Uhlenbeck processes. Each of the subenvironments of the composite environment is able to induce non-Markovian dynamics of the two-level system. By means of the numerically exact method of hierarchical equations of motion, we study steady states of the two-level system, evolution of the reduced density matrix and the equilibrium emission spectra in dependence on the frequency cutoffs and the coupling strengths of the subenvironments. Additionally, we investigate the impact of the rotating wave approximation (RWA) for the interaction with the bath on accuracy of the results.