Thermalization of open quantum systems using the multiple-Davydov-D2 variational approach

Abstract

Numerical implementation of an explicit phonon bath requires a large number of oscillator modes in order to maintain oscillators at the initial temperature when modeling energy relaxation processes. An additional thermalization algorithm may be useful in controlling the local temperature. In this paper we extend our previously proposed thermalization algorithm [M. Jakučionis and D. Abramavičius, Phys. Rev. A 103, 032202 (2021) ] to be used with the numerically exact multiple-Davydov-D2 trial wave function for simulation of relaxation dynamics and spectroscopic signals of open quantum systems using the time-dependent Dirac-Frenkel variational principle. By applying it to the molecular aggregate model, we demonstrate how the thermalization approach significantly reduces the numerical cost of simulations by decreasing the number of oscillators needed to explicitly simulate the aggregate's environment fluctuations while maintaining correspondence to the exact population relaxation dynamics. Additionally, we show how the thermalization can be used to find the equilibrium state of the excited molecular aggregate, which is necessary for simulation of the fluorescence and other spectroscopic signals. The thermalization algorithm we present offers the possibility to investigate larger system-bath models than was previously possible using the multiple-Davydov-D2 trial wave function and local heating effects in molecular complexes.

References in corpus (9)

• The quantum Rabi model: solution and dynamics
• The equilibrium states of open quantum systems in the strong coupling regime
• Polaron dynamics with a multitude of Davydov D${}_2$ trial states
• Excitons in Molecular Aggregates with Lévy Disorder: Anomalous Localization and Exchange Broadening of Optical Spectra
• Schrödinger-cat states in Landau-Zener-Stückelberg-Majorana interferometry: a multiple Davydov Ansatz approach
• Modeling Molecular J and H Aggregates using Multiple-Davydov D2 Ansatz
• Simulation of ab initio optical absorption spectrum of β-carotene with fully resolved S_{0} and S_{2} vibrational normal modes
• Temperature Controlled Open Quantum System Dynamics using Time-dependent Variational Method
• Inspecting molecular aggregate quadratic vibronic coupling effects using squeezed coherent states