Numerical time propagation of quantum systems in radiation fields

Abstract

Atoms, molecules or excitonic quasiparticles, for which excitations are induced by external radiation fields and energy is dissipated through radiative decay, are examples of driven open quantum systems. We explain the use of commutator-free exponential time propagators for the numerical solution of the associated Schrödinger or master equations with a time-dependent Hamilton operator. These time propagators are based on the Magnus series but avoid the computation of commutators, which makes them suitable for the efficient propagation of systems with a large number of degrees of freedom. We present an optimized fourth-order propagator and demonstrate its efficiency in comparison to the direct Runge-Kutta computation. As an illustrative example we consider the parametrically driven dissipative Dicke model, for which we calculate the periodic steady state and the optical emission spectrum. © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.