Calculating two-dimensional spectra with the mixed quantum-classical Ehrenfest method

We present a mixed quantum-classical simulation approach to calculate two-dimensional spectra of coupled two-level electronic model systems. We include the change in potential energy of the classical system due to transitions in the quantum system using the Ehrenfest method. We study how this feedback of the quantum system on the classical system influences the shape of two-dimensional spectra. We show that the feedback leads to the expected Stokes shift of the energy levels in the quantum system. This subsequently leads to changes in the population transfer between quantum sites, which in turn influence the intensities of the peaks in two-dimensional spectra. The obtained spectra are compared with spectra calculated using the Hierarchical Equations of Motion method which is exact. While the spectra match perfectly for short waiting times, clear differences are found for longer waiting times. This is attributed to a violation of detailed balance between the quantum states in the Ehrenfest method. The energy of the total quantum-classical system however does obey a Boltzmann distribution, when coupled to a stochastic heat bath.