Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Delving Into Dissipative Quantum Dynamics: From Approximate to Numerically Exact Approaches

Hsing-Ta Chen

Published Jan 1, 2016• 0 citations
Publisher

Abstract

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.

DOI
10.7916/d8m32wbn

Related papers

Inchworm Monte Carlo for exact non-adiabatic dynamics. I. Theory and algorithms

Hsing-Ta Chen, Guy Cohen, David R. Reichman

In this paper, we provide a detailed description of the inchworm Monte Carlo formalism for the exact study of real-time non-adiabatic dynamics. This method optimally recycles Monte Carlo information from earlier times to greatly suppress the dynamical sign problem. Using the example of the spin-boson model, we 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. The latter approach motivates the development of a cumulant version of the inchworm Monte Carlo method, which has the benefit of improved scaling. This paper deals completely with methodology, while Paper II provides 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.

arXiv:1610.09566•The Journal of Chemical Physics•Feb 1, 2017•82 citations

Inchworm Monte Carlo for exact non-adiabatic dynamics. II. Benchmarks and comparison with established methods

Hsing-Ta Chen, Guy Cohen, David R. Reichman

In this second paper of a two part series, we present extensive benchmark results for two different inchworm Monte Carlo expansions for the spin-boson model. Our results are compared to previously developed numerically exact approaches for this problem. A detailed discussion of convergence and error propagation is presented. Our results and analysis allow for an understanding of the benefits and drawbacks of inchworm Monte Carlo compared to other approaches for exact real-time non-adiabatic quantum dynamics.

arXiv:1610.09402•The Journal of Chemical Physics•Feb 1, 2017•66 citations

Nonequilibrium Dynamics of Strongly Correlated Systems

Joseph Kleinhenz

In this thesis, we describe methods for solving quantum impurity models out of equilibrium, and consider several applications of these methods. Quantum impurity models describe systems in which a small subsystem (the impurity) is coupled to a large external environment (the bath). These models are useful in describing a wide variety of physical systems, including quantum dots, absorbed adatoms, and molecular junctions. They are also important for the study of lattice models, such as the Hubbard model, where they appear as auxiliary problems within the framework of dynamical mean field theory (DMFT). Many interesting phenomena in these systems, such as transport through a quantum dot, or optical pumping of a Mott insulator, occur away from equilibrium. Theoretical description of these phenomena requires methods that are capable of handling both real time evolution and strong correlation effects presenting a formidable challenge. In this thesis, we describe several methods for the solution of non-equilibrium quantum impurity models which meet this challenge. These methods are based on an expansion in the coupling between the impurity and the bath known as the hybridization expansion. Chapter 1 gives an introduction to the problems addressed in this thesis. Chapter 2 provides a brief overview of Keldysh formalism for nonequilibrium problems. Chapter 3 describes the bare hybridization expansion for the Anderson impurity model as well as its bold extensions. In particular we derive the non-crossing approximation (NCA), the one-crossing approximation (OCA), the inchworm expansion, and a bold expansion for the impurity Green's function. Chapter 4 describes the implementation of a quantum Monte Carlo algorithm based on the inchworm expansion. The rest of the thesis describes applications of these algorithms to problems of physical interest. Chapters 5 and 6 apply the methods previously described to understanding the Kondo effect in quantum dot systems. The first chapter examines the behavior of the Kondo effect under an applied voltage bias, a classic non-equilibrium strongly-correlated problem. The second chapter examines different ways of measuring the cloud of electrons which screen the impurity spin in the Kondo regime. Chapters 7 and 8 apply the same methods to the Hubbard model using the Dynamical mean-field theory (DMFT) mapping. The first chapter studies the Hubbard model in equilibrium using the numerically exact inchworm quantum Monte Carlo method. The second chapter studies the possibility of dynamic switching between the metal and Mott-insulator states of the Hubbard model using the one-crossing approximation (OCA). Finally chapter 9 gives some overall conclusions and suggests future directions for the development of hybridization expansion methods.

•Jan 1, 2021•0 citations

Non-Equilibrium Dynamical Mean Field Theory with Quantum Monte Carlo

Qiaoyuan Dong

This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.

•Jan 1, 2017•0 citations

Anderson-Holstein model in two flavors of the noncrossing approximation

Hsing-Ta Chen, Guy Cohen, Andrew J. Millis +1 more

The dynamical interplay between electron-electron interactions and electron-phonon coupling is investigated within the Anderson-Holstein model, a minimal model for open quantum systems that embody these effects. The influence of phonons on spectral and transport properties is explored in equilibrium, for nonequilibrium steady state and for transient dynamics after a quench. Both the particle-hole symmetric and the more generic particle-hole asymmetric cases are studied. The treatment is based on two complementary noncrossing approximations, the first of which is constructed around the weak-coupling limit and the second around the polaron limit. In general, the two methods disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The frameworks used here can form the starting point for numerically exact methods based on bold-line continuous-time quantum Monte Carlo algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.

arXiv:1601.05755•Physical review. B./Physical review. B•May 23, 2016•50 citations

Modèles d'impuretés quantiques dissipatifs

Matthieu Vanhoecke

This thesis focuses on the study of dissipative quantum impurity models. Quantum impurity models describe a localized quantum system, known as the impurity, that interacts with a large quantum bath, enabling the exploration of essential phenomena such as the Kondo effect in magnetic impurities and the transport properties in quantum dots. An intriguing extension of these models is the exploration of the interplay between strong correlations and dissipation. Dissipation, often described by the Lindblad master equation in its Markovian limit, introduces challenges and opportunities as it interacts with coherent many-body effects. In this thesis, we want to investigate the non-equilibrium dynamics of quantum impurity systems influenced by both strong correlations and dissipation across four main parts. We begin by developing numerical methods, specifically a diagrammatic Monte Carlo technique, to analyze the real-time dynamics of dissipative quantum impurity models. Our findings demonstrate that local Markovian dissipation enhances the convergence of Monte Carlo sampling by constraining the process to a reduced subspace of diagrams. This effectively mitigates the sign problem and facilitates the exploration of longer time scales compared to conventional unitary approaches. In the second part, we derive a self-consistent dynamical map using the Non-Crossing Approximation to study the long-term behavior of dissipative quantum impurity models. This method effectively resums diagrams into a self-consistent Dyson equation, capturing essential non-Markovian environmental effects, including the Kondo effect. In a third part, we examine the relationship between Kondo physics and Markovian dissipation by investigating an interacting quantum dot coupled to a metallic bath under continuous charge monitoring. We identify a crossover from Kondo screening to the Quantum Zeno effect, effectively described by a non-Hermitian Kondo model in the weak dissipation regime. Finally, the last part explores the effects of particle losses on the Kondo effect and the Kondo-Zeno crossover. Preliminary results indicate that while single-particle losses reduce many-body effects, two-body losses introduce intriguing phenomena, such as spin decay and hints of the Kondo-Zeno crossover, along with an unusual spectral function exhibiting a low-frequency lineshape.

•Nov 25, 2024•0 citations
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.
In this paper, we provide a detailed description of the inchworm Monte Carlo formalism for the exact study of real-time non-adiabatic dynamics. This method optimally recycles Monte Carlo information from earlier times to greatly suppress the dynamical sign problem. Using the example of the spin-boson model, we 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. The latter approach motivates the development of a cumulant version of the inchworm Monte Carlo method, which has the benefit of improved scaling. This paper deals completely with methodology, while Paper II provides 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.
In this second paper of a two part series, we present extensive benchmark results for two different inchworm Monte Carlo expansions for the spin-boson model. Our results are compared to previously developed numerically exact approaches for this problem. A detailed discussion of convergence and error propagation is presented. Our results and analysis allow for an understanding of the benefits and drawbacks of inchworm Monte Carlo compared to other approaches for exact real-time non-adiabatic quantum dynamics.
In this thesis, we describe methods for solving quantum impurity models out of equilibrium, and consider several applications of these methods. Quantum impurity models describe systems in which a small subsystem (the impurity) is coupled to a large external environment (the bath). These models are useful in describing a wide variety of physical systems, including quantum dots, absorbed adatoms, and molecular junctions. They are also important for the study of lattice models, such as the Hubbard model, where they appear as auxiliary problems within the framework of dynamical mean field theory (DMFT). Many interesting phenomena in these systems, such as transport through a quantum dot, or optical pumping of a Mott insulator, occur away from equilibrium. Theoretical description of these phenomena requires methods that are capable of handling both real time evolution and strong correlation effects presenting a formidable challenge. In this thesis, we describe several methods for the solution of non-equilibrium quantum impurity models which meet this challenge. These methods are based on an expansion in the coupling between the impurity and the bath known as the hybridization expansion. Chapter 1 gives an introduction to the problems addressed in this thesis. Chapter 2 provides a brief overview of Keldysh formalism for nonequilibrium problems. Chapter 3 describes the bare hybridization expansion for the Anderson impurity model as well as its bold extensions. In particular we derive the non-crossing approximation (NCA), the one-crossing approximation (OCA), the inchworm expansion, and a bold expansion for the impurity Green's function. Chapter 4 describes the implementation of a quantum Monte Carlo algorithm based on the inchworm expansion. The rest of the thesis describes applications of these algorithms to problems of physical interest. Chapters 5 and 6 apply the methods previously described to understanding the Kondo effect in quantum dot systems. The first chapter examines the behavior of the Kondo effect under an applied voltage bias, a classic non-equilibrium strongly-correlated problem. The second chapter examines different ways of measuring the cloud of electrons which screen the impurity spin in the Kondo regime. Chapters 7 and 8 apply the same methods to the Hubbard model using the Dynamical mean-field theory (DMFT) mapping. The first chapter studies the Hubbard model in equilibrium using the numerically exact inchworm quantum Monte Carlo method. The second chapter studies the possibility of dynamic switching between the metal and Mott-insulator states of the Hubbard model using the one-crossing approximation (OCA). Finally chapter 9 gives some overall conclusions and suggests future directions for the development of hybridization expansion methods.
This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.
The dynamical interplay between electron-electron interactions and electron-phonon coupling is investigated within the Anderson-Holstein model, a minimal model for open quantum systems that embody these effects. The influence of phonons on spectral and transport properties is explored in equilibrium, for nonequilibrium steady state and for transient dynamics after a quench. Both the particle-hole symmetric and the more generic particle-hole asymmetric cases are studied. The treatment is based on two complementary noncrossing approximations, the first of which is constructed around the weak-coupling limit and the second around the polaron limit. In general, the two methods disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The frameworks used here can form the starting point for numerically exact methods based on bold-line continuous-time quantum Monte Carlo algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.
This thesis focuses on the study of dissipative quantum impurity models. Quantum impurity models describe a localized quantum system, known as the impurity, that interacts with a large quantum bath, enabling the exploration of essential phenomena such as the Kondo effect in magnetic impurities and the transport properties in quantum dots. An intriguing extension of these models is the exploration of the interplay between strong correlations and dissipation. Dissipation, often described by the Lindblad master equation in its Markovian limit, introduces challenges and opportunities as it interacts with coherent many-body effects. In this thesis, we want to investigate the non-equilibrium dynamics of quantum impurity systems influenced by both strong correlations and dissipation across four main parts. We begin by developing numerical methods, specifically a diagrammatic Monte Carlo technique, to analyze the real-time dynamics of dissipative quantum impurity models. Our findings demonstrate that local Markovian dissipation enhances the convergence of Monte Carlo sampling by constraining the process to a reduced subspace of diagrams. This effectively mitigates the sign problem and facilitates the exploration of longer time scales compared to conventional unitary approaches. In the second part, we derive a self-consistent dynamical map using the Non-Crossing Approximation to study the long-term behavior of dissipative quantum impurity models. This method effectively resums diagrams into a self-consistent Dyson equation, capturing essential non-Markovian environmental effects, including the Kondo effect. In a third part, we examine the relationship between Kondo physics and Markovian dissipation by investigating an interacting quantum dot coupled to a metallic bath under continuous charge monitoring. We identify a crossover from Kondo screening to the Quantum Zeno effect, effectively described by a non-Hermitian Kondo model in the weak dissipation regime. Finally, the last part explores the effects of particle losses on the Kondo effect and the Kondo-Zeno crossover. Preliminary results indicate that while single-particle losses reduce many-body effects, two-body losses introduce intriguing phenomena, such as spin decay and hints of the Kondo-Zeno crossover, along with an unusual spectral function exhibiting a low-frequency lineshape.