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:1505.01283 · Physical Review B

Transport through an Anderson impurity: Current ringing, nonlinear magnetization, and a direct comparison of continuous-time quantum Monte Carlo and hierarchical quantum master equations

R. Härtle, Guy Cohen, David R. Reichman, Andrew J. Millis

Published Aug 27, 2015• 78 citations
arXivHTMLPDF

Abstract

We give a detailed comparison of the hierarchical quantum master equation (HQME) method to a continuous-time quantum Monte Carlo (CT-QMC) approach, assessing the usability of these numerically exact schemes as impurity solvers in practical nonequilibrium calculations. We review the main characteristics of the methods and discuss the scaling of the associated numerical effort. We substantiate our discussion with explicit numerical results for the nonequilibrium transport properties of a single-site Anderson impurity. The numerical effort of the HQME scheme scales linearly with the simulation time but increases (at worst exponentially) with decreasing temperature. In contrast, CT-QMC is less restricted by temperature at short times, but in general the cost of going to longer times is also exponential. After establishing the numerical exactness of the HQME scheme, we use it to elucidate the influence of different ways to induce transport through the impurity on the initial dynamics, discuss the phenomenon of coherent current oscillations, known as current ringing, and explain the nonmonotonic temperature dependence of the steady-state magnetization as a result of competing broadening effects. We also elucidate the pronounced nonlinear magnetization dynamics, which appears on intermediate time scales in the presence of an asymmetric coupling to the electrodes.

DOI
10.1103/physrevb.92.085430

Related papers

Numerically exact long-time magnetization dynamics at the nonequilibrium Kondo crossover of the Anderson impurity model

Guy Cohen, Emanuel Gull, David R. Reichman +2 more

We investigate the dynamical and steady-state spin response of the nonequilibrium Anderson model to magnetic fields, bias voltage, and temperature using a numerically exact method combining a bold-line quantum Monte Carlo technique with a memory function formalism. We obtain converged results in a range of previously inaccessible regimes, in particular calculating the spin dynamics for a range of temperatures down to the crossover to the Kondo domain. We provide predictions for nonequilibrium phenomena, including nonmonotonic temperature dependence of observables at high bias voltage and oscillatory quench dynamics at high magnetic fields.

arXiv:1212.0546•Physical Review B•May 8, 2013•153 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

Exact real-time dynamics of single-impurity Anderson model from a single-spin hybridization-expansion

Patryk Kubiczek, Alexey Rubtsov, Alexander I. Lichtenstein

In this work we introduce a modified real-time continuous-time hybridization-expansion quantum Monte Carlo solver for a time-dependent single-orbital Anderson impurity model: CT-1/2-HYB-QMC. In the proposed method the diagrammatic expansion is performed only for one out of the two spin channels, while the resulting effective single-particle problem for the other spin is solved semi-analytically for each expansion diagram. CT-1/2-HYB-QMC alleviates the dynamical sign problem by reducing the order of sampled diagrams and makes it possible to reach twice as long time scales in comparison to the standard CT-HYB method. We illustrate the new solver by calculating an electric current through impurity in paramagnetic and spin-polarized cases.

arXiv:1904.12582•SciPost Physics•Aug 2, 2019•14 citations

Extending the hierarchical quantum master equation approach to low temperatures and realistic band structures

André Erpenbeck, Christopher Hertlein, C. Schinabeck +1 more

The hierarchical quantum master equation (HQME) approach is an accurate method to describe quantum transport in interacting nanosystems. It generalizes perturbative master equation approaches by including higher-order contributions as well as non-Markovian memory and allows for the systematic convergence to the numerically exact result. As the HQME method relies on a decomposition of the bath correlation function in terms of exponentials, however, its application to systems at low temperatures coupled to baths with complexer band structures has been a challenge. In this publication, we outline an extension of the HQME approach, which uses re-summation over poles and can be applied to calculate transient currents at a numerical cost that is independent of temperature and band structure of the baths. We demonstrate the performance of the extended HQME approach for noninteracting tight-binding model systems of increasing complexity as well as for the spinless Anderson-Holstein model.

arXiv:1805.10161•The Journal of Chemical Physics•Aug 9, 2018•46 citations

The Anderson impurity model out-of-equilibrium: Assessing the accuracy of simulation techniques with an exact current-occupation relation

Bijay Kumar Agarwalla, Dvira Segal

We study the interacting, symmetrically coupled single impurity Anderson model. By employing the nonequilibrium Green's function formalism, we reach an exact relationship between the steady-state charge current flowing through the impurity (dot) and its occupation. We argue that the steady-state current-occupation relation can be used to assess the consistency of simulation techniques and identify spurious transport phenomena. We test this relation in two different model variants: First, we study the Anderson-Holstein model in the strong electron-vibration coupling limit using the polaronic quantum master equation method. We find that the current-occupation relation is violated numerically in standard calculations, with simulations bringing up incorrect transport effects. Using a numerical procedure, we resolve the problem efficiently. Second, we simulate the Anderson model with electron-electron interaction on the dot using a deterministic numerically exact time-evolution scheme. Here, we observe that the current-occupation relation is satisfied in the steady-state limit-even before results converge to the exact limit.

arXiv:1704.00121•The Journal of Chemical Physics•Aug 4, 2017•6 citations

Transient dynamics of the Anderson impurity model out of equilibrium

Thomas L. Schmidt, Philipp Werner, Lothar Mühlbacher +1 more

We discuss the transient effects in the Anderson impurity model that occur when two fermionic continua with finite bandwidths are instantaneously coupled to a central level. We present results for the analytically solvable noninteracting resonant-level system first and then consistently extend them to the interacting case using the conventional perturbation theory and recently developed nonequilibrium Monte Carlo simulation schemes. The main goal is to gain an understanding of the full time-dependent nonlinear current-voltage characteristics and the population probability of the central level. We find that, contrary to the steady state, the transient dynamics of the system depends sensitively on the bandwidth of the electrode material.

arXiv:0808.0442•Physical Review B•Dec 10, 2008•122 citations
We give a detailed comparison of the hierarchical quantum master equation (HQME) method to a continuous-time quantum Monte Carlo (CT-QMC) approach, assessing the usability of these numerically exact schemes as impurity solvers in practical nonequilibrium calculations. We review the main characteristics of the methods and discuss the scaling of the associated numerical effort. We substantiate our discussion with explicit numerical results for the nonequilibrium transport properties of a single-site Anderson impurity. The numerical effort of the HQME scheme scales linearly with the simulation time but increases (at worst exponentially) with decreasing temperature. In contrast, CT-QMC is less restricted by temperature at short times, but in general the cost of going to longer times is also exponential. After establishing the numerical exactness of the HQME scheme, we use it to elucidate the influence of different ways to induce transport through the impurity on the initial dynamics, discuss the phenomenon of coherent current oscillations, known as current ringing, and explain the nonmonotonic temperature dependence of the steady-state magnetization as a result of competing broadening effects. We also elucidate the pronounced nonlinear magnetization dynamics, which appears on intermediate time scales in the presence of an asymmetric coupling to the electrodes.
We investigate the dynamical and steady-state spin response of the nonequilibrium Anderson model to magnetic fields, bias voltage, and temperature using a numerically exact method combining a bold-line quantum Monte Carlo technique with a memory function formalism. We obtain converged results in a range of previously inaccessible regimes, in particular calculating the spin dynamics for a range of temperatures down to the crossover to the Kondo domain. We provide predictions for nonequilibrium phenomena, including nonmonotonic temperature dependence of observables at high bias voltage and oscillatory quench dynamics at high magnetic fields.
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.
In this work we introduce a modified real-time continuous-time hybridization-expansion quantum Monte Carlo solver for a time-dependent single-orbital Anderson impurity model: CT-1/2-HYB-QMC. In the proposed method the diagrammatic expansion is performed only for one out of the two spin channels, while the resulting effective single-particle problem for the other spin is solved semi-analytically for each expansion diagram. CT-1/2-HYB-QMC alleviates the dynamical sign problem by reducing the order of sampled diagrams and makes it possible to reach twice as long time scales in comparison to the standard CT-HYB method. We illustrate the new solver by calculating an electric current through impurity in paramagnetic and spin-polarized cases.
The hierarchical quantum master equation (HQME) approach is an accurate method to describe quantum transport in interacting nanosystems. It generalizes perturbative master equation approaches by including higher-order contributions as well as non-Markovian memory and allows for the systematic convergence to the numerically exact result. As the HQME method relies on a decomposition of the bath correlation function in terms of exponentials, however, its application to systems at low temperatures coupled to baths with complexer band structures has been a challenge. In this publication, we outline an extension of the HQME approach, which uses re-summation over poles and can be applied to calculate transient currents at a numerical cost that is independent of temperature and band structure of the baths. We demonstrate the performance of the extended HQME approach for noninteracting tight-binding model systems of increasing complexity as well as for the spinless Anderson-Holstein model.
We study the interacting, symmetrically coupled single impurity Anderson model. By employing the nonequilibrium Green's function formalism, we reach an exact relationship between the steady-state charge current flowing through the impurity (dot) and its occupation. We argue that the steady-state current-occupation relation can be used to assess the consistency of simulation techniques and identify spurious transport phenomena. We test this relation in two different model variants: First, we study the Anderson-Holstein model in the strong electron-vibration coupling limit using the polaronic quantum master equation method. We find that the current-occupation relation is violated numerically in standard calculations, with simulations bringing up incorrect transport effects. Using a numerical procedure, we resolve the problem efficiently. Second, we simulate the Anderson model with electron-electron interaction on the dot using a deterministic numerically exact time-evolution scheme. Here, we observe that the current-occupation relation is satisfied in the steady-state limit-even before results converge to the exact limit.
We discuss the transient effects in the Anderson impurity model that occur when two fermionic continua with finite bandwidths are instantaneously coupled to a central level. We present results for the analytically solvable noninteracting resonant-level system first and then consistently extend them to the interacting case using the conventional perturbation theory and recently developed nonequilibrium Monte Carlo simulation schemes. The main goal is to gain an understanding of the full time-dependent nonlinear current-voltage characteristics and the population probability of the central level. We find that, contrary to the steady state, the transient dynamics of the system depends sensitively on the bandwidth of the electrode material.