We develop a stochastic Schrödinger equation (SSE) framework to simulate the real-time dynamics of Anderson-Holstein (AH) impurities coupled to a continuous Fermionic bath. The bath degrees of freedom are incorporated through fluctuating terms determined by exact system-bath correlations, which is derived in an ab initio manner. We show that such an SSE treatment provides a middle ground between numerically expansive microscopic simulations and macroscopic master equations. Computationally, the SSE model enables efficient numerical methods for propagating stochastic trajectories. We demonstrate that this approach not only naturally provides microscopically detailed information unavailable from reduced models but also captures effects beyond master equations, thus serving as a promising tool to study open quantum dynamics emerging in physics and chemistry.
We develop a stochastic Schrödinger equation (SSE) framework to simulate the real-time dynamics of Anderson-Holstein (AH) impurities coupled to a continuous Fermionic bath. The bath degrees of freedom are incorporated through fluctuating terms determined by exact system-bath correlations, which is derived in an ab initio manner. We show that such an SSE treatment provides a middle ground between numerically expansive microscopic simulations and macroscopic master equations. Computationally, the SSE model enables efficient numerical methods for propagating stochastic trajectories. We demonstrate that this approach not only naturally provides microscopically detailed information unavailable from reduced models but also captures effects beyond master equations, thus serving as a promising tool to study open quantum dynamics emerging in physics and chemistry.
Based on the stochastic unravelling of the reduced density operator in the Feynman path integral formalism for an open quantum system in touch with harmonic environments, a new non-Markovian stochastic Schrödinger equation (NMSSE) has been established that allows for the systematic perturbation expansion in the system-bath coupling to arbitrary order. This NMSSE can be transformed in a facile manner into the other two NMSSEs, i.e., non-Markovian quantum state diffusion and time-dependent wavepacket diffusion method. Benchmarked by numerically exact results, we have conducted a comparative study of the proposed method in its lowest order approximation, with perturbative quantum master equations in the symmetric spin-boson model and the realistic Fenna-Matthews-Olson complex. It is found that our method outperforms the second-order time-convolutionless quantum master equation in the whole parameter regime and even far better than the fourth-order in the slow bath and high temperature cases. Besides, the method is applicable on an equal footing for any kind of spectral density function and is expected to be a powerful tool to explore the quantum dynamics of large-scale systems, benefiting from the wavefunction framework and the time-local appearance within a single stochastic trajectory.
ABSTRACT Stochastic evolution underpins several approaches to the dynamics of open quantum systems, such as random modulation of Hamiltonian parameters, the stochastic Schrödinger equation (SSE), and the stochastic Liouville equation (SLE). These approaches replace the explicit system–environment coupling with an effective system‐only dynamics, where dissipative behavior emerges from ensemble averaging. Stochastic Hamiltonians, in particular, have long served as phenomenological tools in physical chemistry to include environmental effects without recourse to an explicit microscopic derivation. In this work, we aim at a self‐contained and accessible presentation of these approaches to further elaborate on their common roots in essential concepts of stochastic calculus and to delineate the conditions under which they are equivalent. We also discuss how different formulations naturally lead to different numerical time‐integration schemes, better suited for either classical simulation platforms, based on finite‐difference approximations, or quantum algorithms, that employ random unitary maps. Our analysis aims at providing a unified perspective and actionable recipes for classical and quantum implementations of stochastic evolution in the simulation of open quantum systems.
We investigate a simple surface hopping (SH) approach for modeling a single impurity level coupled to a single phonon and an electronic (metal) bath (i.e., the Anderson-Holstein model). The phonon degree of freedom is treated classically with motion along--and hops between--diabatic potential energy surfaces. The hopping rate is determined by the dynamics of the electronic bath (which are treated implicitly). For the case of one electronic bath, in the limit of small coupling to the bath, SH recovers phonon relaxation to thermal equilibrium and yields the correct impurity electron population (as compared with numerical renormalization group). For the case of out of equilibrium dynamics, SH current-voltage (I-V) curve is compared with the quantum master equation (QME) over a range of parameters, spanning the quantum region to the classical region. In the limit of large temperature, SH and QME agree. Furthermore, we can show that, in the limit of low temperature, the QME agrees with real-time path integral calculations. As such, the simple procedure described here should be useful in many other contexts.
The stochastic Schrödinger equation (SSE) provides a trajectory-level route to simulate the dynamics of open quantum systems with applications ranging from molecular processes to quantum technologies. We study a colored-noise extension of the SSE based on an Ornstein-Uhlenbeck (OU) noise drive and benchmark its ensemble-averaged dynamics against the standard white-noise SSE and against a fluctuating OU random Hamiltonian. When the environment exhibits a finite correlation time, averaging over pure-state trajectories yields master equations that are generally open-form and not of Lindblad type, yet remain positive by construction. By considering the differential of the OU process, we define an effective correlated noise, whose properties we analyze and use to formulate an SSE unraveling of its associated open-form quantum master equation. We show that the averaged dissipator separates into a Lindblad contribution stemming from the white-noise component and additional correlation terms arising from the fluctuations of the OU Hamiltonian. To obtain a practical closed description and physical intuition, we introduce a Redfield-inspired perturbative closure for these correlation terms, providing an effective master equation for the colored SSE. For a two-level system, the resulting dynamics exhibit long-lived coherences, nontrivial stationary (including oscillatory) states, and multi-timescale relaxation, rationalized through the components of a time-dependent Redfield tensor.
Stochastic Schrödinger equation (SSE) methods provide a powerful approach for simulating open quantum system dynamics by unraveling the influence functional into an ensemble average of pure-state trajectories. However, long-time simulations using SSE methods are hindered by the inherent growth of ensemble variance. An alternative strategy is to extract essential memory kernels or transfer tensors from short-time dynamics and extrapolate to arbitrarily long times, as in the transfer tensor method (TTM). Conventional TTM implementations, however, require dynamical maps constructed from an informationally complete set of N 2 initial states─where N is the system dimension─limiting their applicability to small-scale problems. Here, we show that the SSE naturally provides an overcomplete set of Kraus operators that captures the essential information needed to construct dynamical maps for TTM. This connection enables the extension of TTM to high-dimensional systems using results obtained from a single converged SSE simulation with a sufficient number of stochastic realizations. We validate our approach with two representative SSE methods applied to a symmetric spin-boson model and a 24-site Fenna–Matthews–Olson (FMO) complex. The results demonstrate that SSE-TTM offers an accurate and computationally efficient strategy for simulating long-time quantum dynamics in relatively large-scale systems.
We propose a method which combines the quantum-classical mapping approach to surface hopping (MASH) with the dissipative quantum dynamics of the Lindblad master equation. Like conventional surface-hopping methods, our approach is based on classical trajectories coupled to the dynamics of a quantum subsystem. However, instead of evolving the subsystem wavefunction according to the time-dependent Schrödinger equation, we use stochastic quantum trajectories derived from secular Redfield theory. This enables the simulation of open quantum systems coupled simultaneously to Markovian quantum baths and anharmonic non-Markovian classical degrees of freedom. Applications to the spin--boson model and to the cavity-enhanced fluorescence of an electronically nonadiabatic molecule show excellent agreement with fully quantum-mechanical benchmarks.