Well-defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé-resummed master equations with numerically exact results in the context of the spin-boson model demonstrates that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques.
Well-defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé-resummed master equations with numerically exact results in the context of the spin-boson model demonstrates that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques.
In this thesis, I explore dissipative quantum dynamics of several prototypical model systems via various approaches, ranging from approximate to numerically exact schemes. In particular, in the realm of the approximate I explore the accuracy of Padé–resummed master equations and the fewest switches surface hopping (FSSH) algorithm for the spin–boson model, and non-crossing approximations (NCA) for the Anderson–Holstein model. Next, I develop new and exact Monte Carlo approaches and test them on the spin–boson model. I propose well–defined criteria for assessing the accuracy of Padé-resummed quantum master equations, which correctly demarcate the regions of parameter space where the Padé approximation is reliable. I continue the investigation of spin–boson dynamics by benchmark comparisons of the semiclassical FSSH algorithm to exact dynamics over a wide range of parameters. Despite small deviations from golden-rule scaling in the Marcus regime, standard surface hopping algorithm is found to be accurate over a large portion of parameter space. The inclusion of decoherence corrections via the augmented FSSH algorithm improves the accuracy of dynamical behavior compared to exact simulations, but the effects are generally not dramatic for the cases I consider. Next, I introduce new methods for numerically exact real-time simulation based on real-time diagrammatic Quantum Monte Carlo (dQMC) and the inchworm algorithm. These methods optimally recycle Monte Carlo information from earlier times to greatly suppress the dynamical sign problem. In the context of the spin–boson model, I formulate the inchworm expansion in two distinct ways: the first with respect to an expansion in the system–bath coupling and the second as an expansion in the diabatic coupling. In addition, a cumulant version of the inchworm Monte Carlo method is motivated by the latter expansion, which allows for further suppression of the growth of the sign error. I provide a comprehensive comparison of the performance of the inchworm Monte Carlo algorithms to other exact methodologies as well as a discussion of the relative advantages and disadvantages of each. Finally, I investigate the dynamical interplay between the electron–electron interaction and the electron–phonon coupling within the Anderson–Holstein model via two complementary NCAs: the first is constructed around the weak-coupling limit and the second around the polaron limit. The influence of phonons on spectral and transport properties is explored in equilibrium, for non-equilibrium steady state and for transient dynamics after a quench. I find the two NCAs disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The complementary frameworks used here pave the way for numerically exact methods based on inchworm dQMC algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.
Generalized master equations provide a concise formalism for studying reduced population dynamics. Usually, these master equations require a perturbative expansion of the memory kernels governing the dynamics; in order to prevent divergences, these expansions must be resummed. Resummation techniques of perturbation series are ubiquitous in physics, but they have not been readily studied for the time-dependent memory kernels used in generalized master equations. In this paper, we present a comparison of different resummation techniques for such memory kernels up to fourth order. We study specifically the spin-boson Hamiltonian as a model system bath Hamiltonian, treating the diabatic coupling between the two states as a perturbation. A novel derivation of the fourth-order memory kernel for the spin-boson problem is presented; then, the second- and fourth-order kernels are evaluated numerically for a variety of spin-boson parameter regimes. We find that resumming the kernels through fourth order using a Padé approximant results in divergent populations in the strong electronic coupling regime due to a singularity introduced by the nature of the resummation, and thus recommend a non-divergent exponential resummation (the "Landau-Zener resummation" of previous work). The inclusion of fourth-order effects in a Landau-Zener-resummed kernel is shown to improve both the dephasing rate and the obedience of detailed balance over simpler prescriptions like the non-interacting blip approximation, showing a relatively quick convergence on the exact answer. The results suggest that including higher-order contributions to the memory kernel of a generalized master equation and performing an appropriate resummation can provide a numerically-exact solution to system-bath dynamics for a general spectral density, opening the way to a new class of methods for treating system-bath dynamics.
The reduced dynamics of an open quantum system obtained from an underlying microscopic Hamiltonian can in general only approximately be described by a time-local master equation. The quality of that approximation depends primarily on the coupling strength and the structure of the environment. Various such master equations have been proposed with different aims. Choosing the most suitable one for a specific system is not straightforward. By focusing on the accuracy of the reduced dynamics we provide a thorough assessment for a selection of methods (Redfield equation, quantum optical master equation, coarse-grained master equation, a related dynamical map approach, and a partial-secular approximation). We consider secondary, here, whether or not an approach guarantees positivity. We use two qubits coupled to a Lorentzian environment in a spin-boson-like fashion modeling a generic situation with various system and bath time scales. We see that, independent of the initial state, the simple Redfield equation with time-dependent coefficients is significantly more accurate than all other methods under consideration. We emphasize that positivity violation in the Redfield equation formalism becomes relevant only in a regime where any of the perturbative master equations considered here are rendered invalid anyway. This implies that the loss of positivity should in fact be welcomed as an important feature: it indicates the breakdown of the weak-coupling assumption. In addition we present the various approaches in a self-contained way and use the behavior of their errors to provide further insight into the range of validity of each method.
The urgent need for reliable simulation tools to match the extreme accuracy needed to control tailored quantum devices highlights the importance of understanding open quantum systems and their modeling. To this end, we compare here the commonly used Redfield and Lindblad master equations against numerically exact results in the case of one and two resonant qubits transversely coupled at a single point to a Drude-cut ohmic bath. All the relevant parameters are varied over a broad range, which allows us to give detailed predictions about the validity and physically meaningful applicability of the weak-coupling approaches. We characterize the accuracy of the approximate approaches by comparing the maximum difference of their system evolution superoperators with numerically exact results. After optimizing the parameters of the approximate models to minimize the difference, we also explore if and to what extent the weak-coupling equations can be applied at least as phenomenological models. Optimization may lead to an accurate reproduction of experimental data, but yet our results are important to estimate the reliability of the extracted parameter values such as the bath temperature. Our findings set general guidelines for the range of validity of the usual Born-Markov master equations and indicate that they fail to accurately describe the physics in a surprisingly broad range of parameters, in particular, at low temperatures. Since quantum-technological devices operate there, their accurate modeling calls for a careful choice of methods.
The Nakajima-Zwanzig generalized master equation provides a formally exact framework to simulate quantum dynamics in condensed phases. Yet, the exact memory kernel is hard to obtain and calculations based on perturbative expansions are often employed. By using the spin-boson model as an example, we assess the convergence of high order memory kernels in the Nakajima-Zwanzig generalized master equation. The exact memory kernels are calculated by combining the hierarchical equation of motion approach and the Dyson expansion of the exact memory kernel. High order expansions of the memory kernels are obtained by extending our previous work to calculate perturbative expansions of open system quantum dynamics [M. Xu et al., J. Chem. Phys. 146, 064102 (2017)]. It is found that the high order expansions do not necessarily converge in certain parameter regimes where the exact kernel show a long memory time, especially in cases of slow bath, weak system-bath coupling, and low temperature. Effectiveness of the Padé and Landau-Zener resummation approaches is tested, and the convergence of higher order rate constants beyond Fermi's golden rule is investigated.
A generalized master equation for the spin-boson model is proposed that does not require a perturbative treatment of the interstate or system-bath couplings. The nonperturbative formulation is based on the assumption of a decoupling of system and bath density operators, and results in a generalized memory function that reduces in the limit of weak system-bath coupling to the standard result of the perturbative non-interacting-blip approximation. Numerical studies in comparison to exact path-integral calculations demonstrate that the nonperturbative master equation represents a clear improvement to perturbation theory as long as the overall coupling is still small enough to justify the underlying decoupling assumption. Finally, possible generalizations of the method to the description of a dissipative N-level system and the limitations of the decoupling ansatz are discussed in some detail.