We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the extended stochastic Liouville--von Neumann equation. Our approach generalizes previous work based on Caldeira-Leggett models and a partitioned initial density matrix. This provides a simple, yet exact, closed-form description for the evolution of open systems from equilibriated initial conditions. The applicability of this model and the potential for numerical implementations are also discussed.
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the extended stochastic Liouville--von Neumann equation. Our approach generalizes previous work based on Caldeira-Leggett models and a partitioned initial density matrix. This provides a simple, yet exact, closed-form description for the evolution of open systems from equilibriated initial conditions. The applicability of this model and the potential for numerical implementations are also discussed.
We present an application of the extended stochastic Liouville equation [G. M. G. McCaul et al., Phys. Rev. B 95, 125124 (2017)], which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time-dependent perturbation, as opposed to the usual assumption that the system and the bath are initially partitioned. Our method presents an exact technique capable of accounting for arbitrary parameter regimes of the model. Here, we present our first numerical implementation of the method in the simplest case of a Caldeira-Leggett representation of the bath Hamiltonian, and apply it to a spin-boson system driven from coupled equilibrium. We observe significant behaviors in both the transient dynamics and asymptotic states of the reduced density matrix not present in the usual approximation.
The exact evolution of a system coupled to a complex environment can be described by a stochastic mean-field evolution of the reduced system density. The formalism developed in Ref. [D.Lacroix, Phys. Rev. E77, 041126 (2008)] is illustrated in the Caldeira-Leggett model where a harmonic oscillator is coupled to a bath of harmonic oscillators. Similar exact reformulation could be used to extend mean-field transport theories in Many-body systems and incorporate two-body correlations beyond the mean-field one. The connection between open quantum system and closed many-body problem is discussed.
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.
Exact solutions of the Caldeira-Leggett Master equation for the reduced density matrix for a free particle and for a harmonic oscillator system coupled to a heat bath of oscillators are obtained for arbitrary initial conditions. The solutions prove that the Fourier transform of the density matrix at time t with respect to (x + x')/2, where x and x' are the initial and final coordinates, factorizes exactly into a part depending linearly on the initial density matrix and a part independent of it. The theorem yields the exact initial state dependence of the density operator at time t and its eventual diagonalization in the energy basis.
We present an application of the extended stochastic Liouville--von Neumann equation (ESLN) method introduced earlier [G. M. G. McCaul, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 95, 125124 (2017); G. M. G. McCaul, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 97, 224310 (2018)], which describes the dynamics of an exactly thermalized open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalized at finite temperature using an imaginary-time evolution procedure before the application of real-time evolution. This initializes the combined system in the correct canonical equilibrium state rather than being initially decoupled. Here we apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open-system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher-order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we first consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under nonequilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are nonzero, initially thermalizing the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than starting in a pure state.
We present a stochastic projection formalism for the description of quantum dynamics in bosonic or spin environments. The Schrödinger equation in the coherent state representation with respect to the environmental degrees of freedom can be reformulated by employing the Feshbach partitioning technique for open quantum systems based on the introduction of suitable non-Hermitian projection operators. In this picture the reduced state of the system can be obtained as a stochastic average over pure state trajectories, for any temperature of the bath. The corresponding non-Markovian stochastic Schrödinger equations include a memory integral over the past states. In the case of harmonic environments and linear coupling the approach gives a new form of the established non-Markovian quantum state diffusion stochastic Schrödinger equation without functional derivatives. Utilizing spin coherent states, the evolution equation for spin environments resembles the bosonic case with, however, a non-Gaussian average for the reduced density operator.