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:1310.4151 · Physical Review Letters

Green’s Functions from Real-Time Bold-Line Monte Carlo Calculations: Spectral Properties of the Nonequilibrium Anderson Impurity Model

Guy Cohen, Emanuel Gull, David R. Reichman, Andrew J. Millis

Published Apr 9, 2014• 108 citations
arXivHTMLPDF

Abstract

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real-time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.

DOI
10.1103/physrevlett.112.146802

Related papers

Green's functions from real-time bold-line Monte Carlo

Guy Cohen, David R. Reichman, Andrew J. Millis +1 more

We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation.

arXiv:1401.0526•Physical Review B•Mar 31, 2014•77 citations

Dynamics of Kondo voltage splitting after a quantum quench

Igor S. Krivenko, Joseph Kleinhenz, Guy Cohen +1 more

We analyze the time-dependent formation of the spectral function of an Anderson impurity model in the Kondo regime within a numerically exact real-time quantum Monte Carlo framework. At steady state, splitting of the Kondo peak occurs with nontrivial dependence on voltage and temperature, and with little effect on the location or intensity of high-energy features. Examining the transient development of the Kondo peak after a quench from an initially uncorrelated state reveals a two-stage process where the initial formation of a single central Kondo peak is followed by splitting. We analyze the time dependence of splitting in detail and demonstrate a strong dependence of its characteristic timescale on the voltage. We expect both the steady state and the transient phenomenon to be experimentally observable.

arXiv:1904.11527•Physical review. B./Physical review. B•Nov 15, 2019•51 citations

Auxiliary master equation approach within matrix product states: Spectral properties of the nonequilibrium Anderson impurity model

Antonius Dorda, Martin Ganahl, Hans Gerd Evertz +2 more

Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.

arXiv:1507.02982•Physical Review B•Sep 24, 2015•79 citations

Time Evolution of the Kondo Resonance in Response to a Quench

Hoa T. M. Nghiem, Theo A. Costi

We investigate the time evolution of the Kondo resonance in response to a quench by applying the time-dependent numerical renormalization group (TDNRG) approach to the Anderson impurity model in the strong correlation limit. For this purpose, we derive within the TDNRG approach a numerically tractable expression for the retarded two-time nonequilibrium Green function $G(t+{t}^{\ensuremath{'}},t)$, and its associated time-dependent spectral function, $A(\ensuremath{\omega},t)$, for times $t$ both before and after the quench. Quenches from both mixed valence and Kondo correlated initial states to Kondo correlated final states are considered. For both cases, we find that the Kondo resonance in the zero temperature spectral function, a preformed version of which is evident at very short times $t\ensuremath{\rightarrow}{0}^{+}$, only fully develops at very long times $t\ensuremath{\gtrsim}1/{T}_{K}$, where ${T}_{K}$ is the Kondo temperature of the final state. In contrast, the final state satellite peaks develop on a fast time scale $1/\mathrm{\ensuremath{\Gamma}}$ during the time interval $\ensuremath{-}1/\mathrm{\ensuremath{\Gamma}}\ensuremath{\lesssim}t\ensuremath{\lesssim}+1/\mathrm{\ensuremath{\Gamma}}$, where $\mathrm{\ensuremath{\Gamma}}$ is the hybridization strength. Initial and final state spectral functions are recovered in the limits $t\ensuremath{\rightarrow}\ensuremath{-}\ensuremath{\infty}$ and $t\ensuremath{\rightarrow}+\ensuremath{\infty}$, respectively. Our formulation of two-time nonequilibrium Green functions within the TDNRG approach provides a first step towards using this method as an impurity solver within nonequilibrium dynamical mean field theory.

arXiv:1701.07558•Physical Review Letters•Oct 13, 2017•39 citations

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

Self-energy method for time-dependent spectral functions of the Anderson impurity model within the time-dependent numerical renormalization group approach

Hoa T. M. Nghiem, Theo A. Costi

The self-energy method for quantum impurity models expresses the correlation part of the self-energy in terms of the ratio of two Green's functions and allows for a more accurate calculation of equilibrium spectral functions than is possible directly from the one-particle Green's function [Bulla et al., J. Phys.: Condens. Matter 10, 8365 (1998)], for example, within the numerical renormalization group method. In addition, the self-energy itself is a central quantity required in the dynamical mean field theory of strongly correlated lattice models. Here, we show how to generalize the self-energy method to the time-dependent situation for the prototype model of strong correlations, the Anderson impurity model. We use the equation-of-motion method to obtain closed expressions for the local Green's function in terms of a time-dependent correlation self-energy, with the latter being given as a ratio of a one-particle time-dependent Green's function and a higher-order correlation function. We benchmark this self-energy approach to time-dependent spectral functions against the direct approach within the time-dependent numerical renormalization group method. The self-energy approach improves the accuracy of time-dependent spectral function calculations, and the closed-form expressions for the Green's function allow for a clear picture of the time evolution of spectral features at the different characteristic time scales. The self-energy approach is of potential interest also for other quantum impurity solvers for real-time evolution, including time-dependent density matrix renormalization group and continuous-time quantum Monte Carlo techniques.

arXiv:2109.01580•Physical review. B./Physical review. B•Nov 12, 2021•3 citations
The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real-time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.
We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation.
We analyze the time-dependent formation of the spectral function of an Anderson impurity model in the Kondo regime within a numerically exact real-time quantum Monte Carlo framework. At steady state, splitting of the Kondo peak occurs with nontrivial dependence on voltage and temperature, and with little effect on the location or intensity of high-energy features. Examining the transient development of the Kondo peak after a quench from an initially uncorrelated state reveals a two-stage process where the initial formation of a single central Kondo peak is followed by splitting. We analyze the time dependence of splitting in detail and demonstrate a strong dependence of its characteristic timescale on the voltage. We expect both the steady state and the transient phenomenon to be experimentally observable.
Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.
We investigate the time evolution of the Kondo resonance in response to a quench by applying the time-dependent numerical renormalization group (TDNRG) approach to the Anderson impurity model in the strong correlation limit. For this purpose, we derive within the TDNRG approach a numerically tractable expression for the retarded two-time nonequilibrium Green function $G(t+{t}^{\ensuremath{'}},t)$, and its associated time-dependent spectral function, $A(\ensuremath{\omega},t)$, for times ttt both before and after the quench. Quenches from both mixed valence and Kondo correlated initial states to Kondo correlated final states are considered. For both cases, we find that the Kondo resonance in the zero temperature spectral function, a preformed version of which is evident at very short times $t\ensuremath{\rightarrow}{0}^{+}$, only fully develops at very long times $t\ensuremath{\gtrsim}1/{T}_{K}$, where TK{T}_{K}TK​ is the Kondo temperature of the final state. In contrast, the final state satellite peaks develop on a fast time scale $1/\mathrm{\ensuremath{\Gamma}}$ during the time interval $\ensuremath{-}1/\mathrm{\ensuremath{\Gamma}}\ensuremath{\lesssim}t\ensuremath{\lesssim}+1/\mathrm{\ensuremath{\Gamma}}$, where $\mathrm{\ensuremath{\Gamma}}$ is the hybridization strength. Initial and final state spectral functions are recovered in the limits $t\ensuremath{\rightarrow}\ensuremath{-}\ensuremath{\infty}$ and $t\ensuremath{\rightarrow}+\ensuremath{\infty}$, respectively. Our formulation of two-time nonequilibrium Green functions within the TDNRG approach provides a first step towards using this method as an impurity solver within nonequilibrium dynamical mean field theory.
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.
The self-energy method for quantum impurity models expresses the correlation part of the self-energy in terms of the ratio of two Green's functions and allows for a more accurate calculation of equilibrium spectral functions than is possible directly from the one-particle Green's function [Bulla et al., J. Phys.: Condens. Matter 10, 8365 (1998)], for example, within the numerical renormalization group method. In addition, the self-energy itself is a central quantity required in the dynamical mean field theory of strongly correlated lattice models. Here, we show how to generalize the self-energy method to the time-dependent situation for the prototype model of strong correlations, the Anderson impurity model. We use the equation-of-motion method to obtain closed expressions for the local Green's function in terms of a time-dependent correlation self-energy, with the latter being given as a ratio of a one-particle time-dependent Green's function and a higher-order correlation function. We benchmark this self-energy approach to time-dependent spectral functions against the direct approach within the time-dependent numerical renormalization group method. The self-energy approach improves the accuracy of time-dependent spectral function calculations, and the closed-form expressions for the Green's function allow for a clear picture of the time evolution of spectral features at the different characteristic time scales. The self-energy approach is of potential interest also for other quantum impurity solvers for real-time evolution, including time-dependent density matrix renormalization group and continuous-time quantum Monte Carlo techniques.