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:1304.6919 · physica status solidi (b)

Iterative path integral summation for nonequilibrium quantum transport

S. Weiss, R. Hützen, D. Becker, Jens Eckel, Reinhold Egger, Michael Thorwart

Published Sep 10, 2013• 26 citations
arXivHTMLPDF

Abstract

We have developed a numerically exact approach to compute real‐time path integral expressions for quantum transport problems out of equilibrium. The scheme is based on a deterministic iterative summation of the path integral (ISPI) for the generating function of nonequilibrium observables of interest. We apply the scheme to compute the charge current or dynamical quantities of small strongly correlated quantum systems as a quantum dot or molecules that are tunnel coupled to leads. Self‐energies due to the leads, being nonlocal in time, are fully taken into account within a finite memory time, thereby including non‐Markovian effects. Numerical results are extrapolated first to vanishing (Trotter) time discretization and, second, to infinite memory time. The method is applied to nonequilibrium transport through a single‐impurity Anderson dot in the first place. We benchmark our results in various regimes of the rich parameter space. In the respective regime of validity, iterative summation of real‐time path integrals (ISPI) results are shown to match those of other state‐of‐the art methods. Especially, we have chosen the mixed valence regime of the Anderson model to compare ISPI to time dependent density matrix renormalization group (tDMRG) and functional RG calculations. Secondly, we determine the nonequilibrium current through a molecular junction in presence of a vibrational mode. We have found an exact mapping of the single impurity Anderson–Holstein model to an effective spin‐1 problem. In analytically tractable regimes, as the adiabatic phonon or weak molecule‐lead coupling regime, we reproduce known perturbative results. Studying the crossover regime between those limits shows that the Franck–Condon blockade persists in the quantum limit. At low temperature, the Franck–Condon steps in the characteristics are smeared due to nonequilibrium conditions. The third system under investigation here is the magnetic Anderson model which consists of a spinful single‐orbital quantum dot with an incorporated quantum mechanical spin–1/2 magnetic impurity. Coulomb interaction together with the exchange coupling of the magnetic impurity with the electron spins strongly influence the dynamics. We investigate the nonequilibrium tunneling current through the system as a function of exchange and Coulomb interaction as well as the real‐time impurity polarization. From the real‐time evolution of physical observables, we are able to determine characteristics of the time‐dependent nonequilibrium current and the relaxation dynamics of the impurity. These examples illustrate that the ISPI technique is particularly well suited for the quantum regime, when all time and energy scales are of the same order of magnitude.

DOI
10.1002/pssb.201349187

Related papers

Iterative real-time path integral approach to nonequilibrium quantum transport

S. Weiss, Jens Eckel, Michael Thorwart +1 more

We have developed a numerical approach to compute real-time path integral expressions for quantum transport problems out of equilibrium. The scheme is based on a deterministic iterative summation of the path integral for the generating function of the nonequilibrium current. Self-energies due to the leads, being nonlocal in time, are fully taken into account within a finite memory time, thereby including non-Markovian effects, and numerical results are extrapolated both to vanishing (Trotter) time discretization and to infinite memory time. This extrapolation scheme converges except at very low temperatures, and the results are then numerically exact. The method is applied to nonequilibrium transport through an Anderson dot.

arXiv:0802.3374•Physical Review B•May 15, 2008•211 citations

Non-equilibrium quantum dynamics of the magnetic Anderson model

D Becker, S Weiss, Michael Thorwart +1 more

We study the non-equilibrium dynamics of a spinful single-orbital quantum dot with an incorporated quantum mechanical spin-1/2 magnetic impurity. Due to the spin degeneracy, double occupancy is allowed, and Coulomb interaction together with the exchange coupling of the magnetic impurity influence the dynamics. By extending the iterative summation of real-time path integrals (ISPI) to this coupled system, we monitor the time-dependent non-equilibrium current and the impurity spin polarization to determine features of the time-dependent non-equilibrium dynamics. We particulary focus on the deep quantum regime, where all time and energy scales are of the same order of magnitude and no small parameter is available. We observe a significant influence of the non-equilibrium decay of the impurity spin polarization both in the presence and in the absence of Coulomb interaction. The exponential relaxation is faster for larger bias voltages, electron–impurity interactions and temperatures. We show that the exact relaxation rate deviates from the corresponding perturbative result. In addition, we study in detail the impurity's back action on the charge current and find a reduction of the stationary current for increasing coupling to the impurity. Moreover, our approach allows us to systematically distinguish mean-field Coulomb and impurity effects from the influence of quantum fluctuations and flip-flop scattering, respectively. In fact, we find a local maximum of the current for a finite Coulomb interaction due to the presence of the impurity.

arXiv:1205.2462•New Journal of Physics•Jul 26, 2012•10 citations

Comparative study of theoretical methods for non-equilibrium quantum transport

Jens Eckel, Fabian Heidrich‐Meisner, Severin G. Jakobs +3 more

We present a detailed comparison of three different methods designed to tackle non-equilibrium quantum transport, namely the functional renormalization group (fRG), the time-dependent density matrix renormalization group (tDMRG) and the iterative summation of real-time path integrals (ISPI). For the non-equilibrium single-impurity Anderson model (including a Zeeman term at the impurity site), we demonstrate that the three methods are in quantitative agreement over a wide range of parameters at the particle–hole symmetric point as well as in the mixed-valence regime. We further compare these techniques with two quantum Monte Carlo approaches and the time-dependent numerical renormalization group method.

arXiv:1001.3773•New Journal of Physics•Apr 28, 2010•126 citations

Numerically exact path-integral simulation of nonequilibrium quantum transport and dissipation

Dvira Segal, Andrew J. Millis, David R. Reichman

We develop an iterative, numerically exact approach for the treatment of nonequilibrium quantum transport and dissipation problems that avoids the real-time sign problem associated with standard Monte Carlo techniques. The method requires a well-defined decorrelation time of the nonlocal influence functional for proper convergence to the exact limit. Since finite decorrelation times may arise either from temperature or from a voltage drop at zero temperature, the approach is well suited for the description of the real-time dynamics of single-molecule devices and quantum dots driven to a steady state via interaction with two or more electron leads. We numerically investigate two nontrivial models: the evolution of the nonequilibrium population of a two-level system coupled to two electronic reservoirs, and quantum transport in the nonequilibrium Anderson model. For the latter case, two distinct formulations are described. Results are compared to those obtained by other techniques.

arXiv:1008.5200•Physical Review B•Nov 18, 2010•160 citations

Real-time effective-action approach to the Anderson quantum dot

Dénes Sexty, Thomas Gasenzer, Jan Martin Pawlowski

The nonequilibrium time evolution of an Anderson quantum dot is investigated. The quantum dot is coupled between two leads forming a chemical-potential gradient. We use Kadanoff-Baym dynamic equations within a nonperturbative resummation of the $s$-channel bubble chains. The effect of the resummation leads to the introduction of a frequency-dependent four-point vertex. The tunneling to the leads is taken into account exactly. The method allows the determination of the transient as well as stationary transport through the quantum dot, and results are compared with different schemes discussed in the literature (functional renormalization group, iterative real-time summation of the path integral, time-dependent matrix renormalization group, and quantum Monte Carlo methods).

arXiv:1012.4293•Physical Review B•Apr 22, 2011•7 citations

Open Wilson chain numerical renormalization group approach to steady-state nonequilibrium quantum transport

Anand Manaparambil, Frithjof B. Anders

The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as $T_K/D \approx 10^{-8}$ in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.

arXiv:2606.14635•Physical review. B./Physical review. B•Oct 23, 2026•0 citations
We have developed a numerically exact approach to compute real‐time path integral expressions for quantum transport problems out of equilibrium. The scheme is based on a deterministic iterative summation of the path integral (ISPI) for the generating function of nonequilibrium observables of interest. We apply the scheme to compute the charge current or dynamical quantities of small strongly correlated quantum systems as a quantum dot or molecules that are tunnel coupled to leads. Self‐energies due to the leads, being nonlocal in time, are fully taken into account within a finite memory time, thereby including non‐Markovian effects. Numerical results are extrapolated first to vanishing (Trotter) time discretization and, second, to infinite memory time. The method is applied to nonequilibrium transport through a single‐impurity Anderson dot in the first place. We benchmark our results in various regimes of the rich parameter space. In the respective regime of validity, iterative summation of real‐time path integrals (ISPI) results are shown to match those of other state‐of‐the art methods. Especially, we have chosen the mixed valence regime of the Anderson model to compare ISPI to time dependent density matrix renormalization group (tDMRG) and functional RG calculations. Secondly, we determine the nonequilibrium current through a molecular junction in presence of a vibrational mode. We have found an exact mapping of the single impurity Anderson–Holstein model to an effective spin‐1 problem. In analytically tractable regimes, as the adiabatic phonon or weak molecule‐lead coupling regime, we reproduce known perturbative results. Studying the crossover regime between those limits shows that the Franck–Condon blockade persists in the quantum limit. At low temperature, the Franck–Condon steps in the characteristics are smeared due to nonequilibrium conditions. The third system under investigation here is the magnetic Anderson model which consists of a spinful single‐orbital quantum dot with an incorporated quantum mechanical spin–1/2 magnetic impurity. Coulomb interaction together with the exchange coupling of the magnetic impurity with the electron spins strongly influence the dynamics. We investigate the nonequilibrium tunneling current through the system as a function of exchange and Coulomb interaction as well as the real‐time impurity polarization. From the real‐time evolution of physical observables, we are able to determine characteristics of the time‐dependent nonequilibrium current and the relaxation dynamics of the impurity. These examples illustrate that the ISPI technique is particularly well suited for the quantum regime, when all time and energy scales are of the same order of magnitude.
We have developed a numerical approach to compute real-time path integral expressions for quantum transport problems out of equilibrium. The scheme is based on a deterministic iterative summation of the path integral for the generating function of the nonequilibrium current. Self-energies due to the leads, being nonlocal in time, are fully taken into account within a finite memory time, thereby including non-Markovian effects, and numerical results are extrapolated both to vanishing (Trotter) time discretization and to infinite memory time. This extrapolation scheme converges except at very low temperatures, and the results are then numerically exact. The method is applied to nonequilibrium transport through an Anderson dot.
We study the non-equilibrium dynamics of a spinful single-orbital quantum dot with an incorporated quantum mechanical spin-1/2 magnetic impurity. Due to the spin degeneracy, double occupancy is allowed, and Coulomb interaction together with the exchange coupling of the magnetic impurity influence the dynamics. By extending the iterative summation of real-time path integrals (ISPI) to this coupled system, we monitor the time-dependent non-equilibrium current and the impurity spin polarization to determine features of the time-dependent non-equilibrium dynamics. We particulary focus on the deep quantum regime, where all time and energy scales are of the same order of magnitude and no small parameter is available. We observe a significant influence of the non-equilibrium decay of the impurity spin polarization both in the presence and in the absence of Coulomb interaction. The exponential relaxation is faster for larger bias voltages, electron–impurity interactions and temperatures. We show that the exact relaxation rate deviates from the corresponding perturbative result. In addition, we study in detail the impurity's back action on the charge current and find a reduction of the stationary current for increasing coupling to the impurity. Moreover, our approach allows us to systematically distinguish mean-field Coulomb and impurity effects from the influence of quantum fluctuations and flip-flop scattering, respectively. In fact, we find a local maximum of the current for a finite Coulomb interaction due to the presence of the impurity.
We present a detailed comparison of three different methods designed to tackle non-equilibrium quantum transport, namely the functional renormalization group (fRG), the time-dependent density matrix renormalization group (tDMRG) and the iterative summation of real-time path integrals (ISPI). For the non-equilibrium single-impurity Anderson model (including a Zeeman term at the impurity site), we demonstrate that the three methods are in quantitative agreement over a wide range of parameters at the particle–hole symmetric point as well as in the mixed-valence regime. We further compare these techniques with two quantum Monte Carlo approaches and the time-dependent numerical renormalization group method.
We develop an iterative, numerically exact approach for the treatment of nonequilibrium quantum transport and dissipation problems that avoids the real-time sign problem associated with standard Monte Carlo techniques. The method requires a well-defined decorrelation time of the nonlocal influence functional for proper convergence to the exact limit. Since finite decorrelation times may arise either from temperature or from a voltage drop at zero temperature, the approach is well suited for the description of the real-time dynamics of single-molecule devices and quantum dots driven to a steady state via interaction with two or more electron leads. We numerically investigate two nontrivial models: the evolution of the nonequilibrium population of a two-level system coupled to two electronic reservoirs, and quantum transport in the nonequilibrium Anderson model. For the latter case, two distinct formulations are described. Results are compared to those obtained by other techniques.
The nonequilibrium time evolution of an Anderson quantum dot is investigated. The quantum dot is coupled between two leads forming a chemical-potential gradient. We use Kadanoff-Baym dynamic equations within a nonperturbative resummation of the sss-channel bubble chains. The effect of the resummation leads to the introduction of a frequency-dependent four-point vertex. The tunneling to the leads is taken into account exactly. The method allows the determination of the transient as well as stationary transport through the quantum dot, and results are compared with different schemes discussed in the literature (functional renormalization group, iterative real-time summation of the path integral, time-dependent matrix renormalization group, and quantum Monte Carlo methods).
The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as TK/D≈10−8T_K/D \approx 10^{-8}TK​/D≈10−8 in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.