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

arXiv:0810.2345 · Physical Review B

Diagrammatic Monte Carlo simulation of nonequilibrium systems

Philipp Werner, Takashi Oka, Andrew J. Millis

Published Jan 21, 2009• 289 citations
arXivHTMLPDF

Abstract

We generalize the recently developed diagrammatic Monte Carlo techniques for quantum impurity models from an imaginary time to a Keldysh formalism suitable for real-time and nonequilibrium calculations. Both weak-coupling and strong-coupling based methods are introduced, analyzed, and applied to the study of transport and relaxation dynamics in interacting quantum dots.

DOI
10.1103/physrevb.79.035320

Related papers

Real-time dynamics in quantum impurity models with diagrammatic Monte Carlo

Marco Schirò

We extend the recently developed real-time diagrammatic Monte Carlo method, in its hybridization expansion formulation, to the full Kadanoff-Baym-Keldysh contour. This allows us to study real-time dynamics in correlated impurity models starting from an arbitrary, even interacting, initial density matrix. As a proof of concept, we apply the algorithm to study the nonequilibrium dynamics after a local quantum quench in the Anderson impurity model. Being a completely general approach to real-time dynamics in quantum impurity models, it can be used as a solver for nonequilibrium dynamical mean field theory.

arXiv:0911.0184•Physical Review B•Feb 26, 2010•72 citations

Real-time diagrammatic Monte Carlo for nonequilibrium quantum transport

Marco Schirò, Michele Fabrizio

We propose a numerically exact approach to nonequilibrium real-time dynamics that is applicable to quantum impurity models coupled to biased noninteracting leads, such as those relevant to quantum transport in nanoscale devices. The method is based on a diagrammatic Monte Carlo sampling of the real-time perturbation theory along the Keldysh contour. We benchmark the method on a noninteracting resonant-level model and, as a first nontrivial application, we study zero-temperature nonequilibrium transport through a vibrating molecule.

arXiv:0808.0589•Physical Review B•Apr 13, 2009•170 citations

Diagrammatic Monte Carlo for dissipative quantum impurity models

Matthieu Vanhoecke, Marco Schirò

We develop a diagrammatic Monte Carlo method for the real-time dynamics of dissipative quantum impurity models. These are small open quantum systems with interaction and local Markovian dissipation, coupled to a large quantum bath. Our algorithm samples the hybridization expansion formulated on a single real-time contour, rather than on the double Keldysh one, as it naturally arises in the thermofield/vectorized representation of the Lindblad dynamics. We show that local Markovian dissipation generally helps the convergence of the diagrammatic Monte Carlo sampling by reducing the sign problem, thus allowing to reach longer timescales as compared to the conventional unitary case. We apply our method to an Anderson impurity model in the presence of local dephasing and discuss its effect on the charge and spin dynamics of the impurity.

arXiv:2311.17839•Physical review. B./Physical review. B•Mar 18, 2024•14 citations

Weak-coupling quantum Monte Carlo calculations on the Keldysh contour: Theory and application to the current-voltage characteristics of the Anderson model

Philipp Werner, Takashi Oka, Martin Eckstein +1 more

We present optimized implementations of the weak-coupling continuous-time Monte Carlo method defined for nonequilibrium problems on the Keldysh contour. We describe and compare two methods of preparing the system before beginning the real-time calculation: the ``interaction quench'' and the ``voltage quench,'' which are found to be suitable for large and small voltage biases, respectively. We also discuss technical optimizations which increase the efficiency of the real-time measurements. The methods allow the accurate simulation of transport through quantum dots over wider interaction ranges and longer times than have heretofore been possible. The current-voltage characteristics of the particle-hole symmetric Anderson-impurity model is presented for interactions $U$ up to ten times the intrinsic level width $\ensuremath{\Gamma}$. We compare the Monte Carlo results to fourth-order perturbation theory, finding that perturbation theory is accurate up to $U\ensuremath{\approx}4\ensuremath{\Gamma}$ or for a voltage bias $V\ensuremath{\gtrsim}4\ensuremath{\Gamma}$. The interplay of voltage and temperature and the Coulomb blockade conductance regime are studied.

arXiv:0911.0587•Physical Review B•Jan 14, 2010•136 citations

Diagrammatic Monte Carlo simulation of quantum impurity models - weak and strong coupling approach

Philipp Werner

An impurity model describes an atom or molecule embedded in some host or bath, with which it can exchange electrons. This exchange of electrons allows the impurity to make transitions between different quantum states, and leads to a non-trivial dynamics. Therefore, despite the zero dimensional nature (which makes impurity problems computationally much more tractable than fermionic lattice models), their numerical simulation remains a challenging task. Methods such as exact diagonalization or numerical RG, which explicitly treat a finite number of bath states, work well for single orbital models. However, because the number of bath states must be increased proportional to the number of orbitals, the computational effort grows exponentially with system size, and requires severe truncations of the bath already for two orbitals. Monte Carlo methods have the advantage that the bath is integrated out and thus the (infinite) size of the bath Hilbert space does not affect the simulation. While restricted to finite temperature, Monte Carlo methods are thus the method of choice for the solution of large multi-orbital or cluster impurity problems. Over the last few years, significant progress has been made (both in terms of efficiency and flexibility) with the development of diagrammatic Monte Carlo techniques. These notes provide an overview over two recently developed, complementary methods: (i) the weakcoupling approach, which scales favorably with system size and allows the efficient simulation of large impurity clusters, and (ii) the strong-coupling approach, which can handle impurity

•Jan 1, 2010•0 citations

Continuous-Time Monte Carlo Methods for Quantum Impu- rity Problems and Dynamical Mean Field Calculations

Philipp Werner, Andrew J. Millis

Andrew J. Millis Department of Physics, Columbia University, 538 West, 120th Street, New York, NY 10027, U.S.A. An atom or molecule embedded in an infinite host material defines a quantum impurity model (QIM). Exchange with the host causes the impurity to make transitions between different quantum states, leading to a nontrivial dynamics which may be described as a quantum field theory in zero space and one time dimension. QIMs exhibit many phenomena of intrinsic theoretical interest including quantum phase transitions (QPTs). The development over the past 15 years of dynamical mean field theory (DMFT) shows that they also serve as an auxiliary problem from whose solution important information relating to the physics of interacting lattice models can be obtained. Numerical simulation of these models remains a challenging task but the development of diagrammatic Monte Carlo techniques over the last few years has led to significant improvements in efficiency and flexibility. Highly precise results on few-orbital models can be computed and large Hilbert-space multiorbital models are now within reach. This chapter provides an overview of recent developments in quantum Monte Carlo methods and their application to multi-orbital systems. It presents both the weak-coupling approach, which scales favorably with system size and allows the efficient simulation of large impurity clusters, and the strong-coupling approach, which can handle QIMs with strong interactions.

CRC Press eBooks•Nov 2, 2010•2 citations
We generalize the recently developed diagrammatic Monte Carlo techniques for quantum impurity models from an imaginary time to a Keldysh formalism suitable for real-time and nonequilibrium calculations. Both weak-coupling and strong-coupling based methods are introduced, analyzed, and applied to the study of transport and relaxation dynamics in interacting quantum dots.
We extend the recently developed real-time diagrammatic Monte Carlo method, in its hybridization expansion formulation, to the full Kadanoff-Baym-Keldysh contour. This allows us to study real-time dynamics in correlated impurity models starting from an arbitrary, even interacting, initial density matrix. As a proof of concept, we apply the algorithm to study the nonequilibrium dynamics after a local quantum quench in the Anderson impurity model. Being a completely general approach to real-time dynamics in quantum impurity models, it can be used as a solver for nonequilibrium dynamical mean field theory.
We propose a numerically exact approach to nonequilibrium real-time dynamics that is applicable to quantum impurity models coupled to biased noninteracting leads, such as those relevant to quantum transport in nanoscale devices. The method is based on a diagrammatic Monte Carlo sampling of the real-time perturbation theory along the Keldysh contour. We benchmark the method on a noninteracting resonant-level model and, as a first nontrivial application, we study zero-temperature nonequilibrium transport through a vibrating molecule.
We develop a diagrammatic Monte Carlo method for the real-time dynamics of dissipative quantum impurity models. These are small open quantum systems with interaction and local Markovian dissipation, coupled to a large quantum bath. Our algorithm samples the hybridization expansion formulated on a single real-time contour, rather than on the double Keldysh one, as it naturally arises in the thermofield/vectorized representation of the Lindblad dynamics. We show that local Markovian dissipation generally helps the convergence of the diagrammatic Monte Carlo sampling by reducing the sign problem, thus allowing to reach longer timescales as compared to the conventional unitary case. We apply our method to an Anderson impurity model in the presence of local dephasing and discuss its effect on the charge and spin dynamics of the impurity.
We present optimized implementations of the weak-coupling continuous-time Monte Carlo method defined for nonequilibrium problems on the Keldysh contour. We describe and compare two methods of preparing the system before beginning the real-time calculation: the ``interaction quench'' and the ``voltage quench,'' which are found to be suitable for large and small voltage biases, respectively. We also discuss technical optimizations which increase the efficiency of the real-time measurements. The methods allow the accurate simulation of transport through quantum dots over wider interaction ranges and longer times than have heretofore been possible. The current-voltage characteristics of the particle-hole symmetric Anderson-impurity model is presented for interactions UUU up to ten times the intrinsic level width $\ensuremath{\Gamma}$. We compare the Monte Carlo results to fourth-order perturbation theory, finding that perturbation theory is accurate up to $U\ensuremath{\approx}4\ensuremath{\Gamma}$ or for a voltage bias $V\ensuremath{\gtrsim}4\ensuremath{\Gamma}$. The interplay of voltage and temperature and the Coulomb blockade conductance regime are studied.
An impurity model describes an atom or molecule embedded in some host or bath, with which it can exchange electrons. This exchange of electrons allows the impurity to make transitions between different quantum states, and leads to a non-trivial dynamics. Therefore, despite the zero dimensional nature (which makes impurity problems computationally much more tractable than fermionic lattice models), their numerical simulation remains a challenging task. Methods such as exact diagonalization or numerical RG, which explicitly treat a finite number of bath states, work well for single orbital models. However, because the number of bath states must be increased proportional to the number of orbitals, the computational effort grows exponentially with system size, and requires severe truncations of the bath already for two orbitals. Monte Carlo methods have the advantage that the bath is integrated out and thus the (infinite) size of the bath Hilbert space does not affect the simulation. While restricted to finite temperature, Monte Carlo methods are thus the method of choice for the solution of large multi-orbital or cluster impurity problems. Over the last few years, significant progress has been made (both in terms of efficiency and flexibility) with the development of diagrammatic Monte Carlo techniques. These notes provide an overview over two recently developed, complementary methods: (i) the weakcoupling approach, which scales favorably with system size and allows the efficient simulation of large impurity clusters, and (ii) the strong-coupling approach, which can handle impurity
Andrew J. Millis Department of Physics, Columbia University, 538 West, 120th Street, New York, NY 10027, U.S.A. An atom or molecule embedded in an infinite host material defines a quantum impurity model (QIM). Exchange with the host causes the impurity to make transitions between different quantum states, leading to a nontrivial dynamics which may be described as a quantum field theory in zero space and one time dimension. QIMs exhibit many phenomena of intrinsic theoretical interest including quantum phase transitions (QPTs). The development over the past 15 years of dynamical mean field theory (DMFT) shows that they also serve as an auxiliary problem from whose solution important information relating to the physics of interacting lattice models can be obtained. Numerical simulation of these models remains a challenging task but the development of diagrammatic Monte Carlo techniques over the last few years has led to significant improvements in efficiency and flexibility. Highly precise results on few-orbital models can be computed and large Hilbert-space multiorbital models are now within reach. This chapter provides an overview of recent developments in quantum Monte Carlo methods and their application to multi-orbital systems. It presents both the weak-coupling approach, which scales favorably with system size and allows the efficient simulation of large impurity clusters, and the strong-coupling approach, which can handle QIMs with strong interactions.