Efficient choice of colored noise in the stochastic dynamics of open quantum systems

The stochastic Liouville-von Neumann (SLN) equation describes the dynamics of an open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. The interaction with the environment is represented by complex colored noises which drive the system, and whose correlation functions are set by the properties of the environment. We present a number of schemes capable of generating colored noises of this kind that are built on a noise amplitude reduction procedure [Imai et al., Chem. Phys. 446, 134 (2015)CMPHC20301-010410.1016/j.chemphys.2014.11.014], including two analytically optimized schemes. In doing so, we pay close attention to the properties of the correlation functions in Fourier space, which we derive in full. For some schemes the method of Wiener filtering for deconvolutions leads to the realization that weakening causality in one of the noise correlation functions improves numerical convergence considerably, allowing us to introduce a well-controlled method for doing so. We compare the ability of these schemes, along with an alternative optimized scheme [Schmitz and Stockburger, Eur. Phys. J.: Spec. Top. 227, 1929 (2019)1951-635510.1140/epjst/e2018-800094-y], to reduce the growth in the mean and variance of the trace of the reduced density matrix, and their ability to extend the region in which the dynamics is stable and well converged for a range of temperatures. By numerically optimizing an additional noise scaling freedom, we identify the scheme which performs best for the parameters used, improving convergence by orders of magnitude and increasing the time accessible by simulation.