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

Imaginary-Time Formulation of Steady-State Nonequilibrium: Application to Strongly Correlated Transport

Jong E. Han, Ryan Heary

Published Dec 7, 2007• 93 citations
arXivHTMLPDF

Abstract

We extend the imaginary-time formulation of the equilibrium quantum many-body theory to steady-state nonequilibrium with an application to strongly correlated transport. By introducing the Matsubara voltage, we maintain the finite chemical potential shifts in the Fermi-Dirac function, in agreement with the Keldysh formulation. The formulation is applied to strongly correlated transport in the Kondo regime using the quantum Monte Carlo method.

DOI
10.1103/physrevlett.99.236808

Related papers

Imaginary-time formulation of steady-state nonequilibrium in quantum dot models

Jong E. Han

We examine the recently proposed imaginary-time formulation for strongly correlated steady-state nonequilibrium for its range of validity and discuss significant improvements in the analytic continuation of the Matsubara voltage as well as the fermionic Matsubara frequency. The discretization error in the conventional Hirsch-Fye algorithm has been compensated in the Fourier transformation with reliable small frequency behavior of self-energy. Here we give detailed discussions for generalized spectral representation ansatz by including high-order vertex corrections and its numerical analytic continuation procedures. The differential conductance calculations agree accurately with available data from other nonequilibrium transport theories. It is verified that, at finite source-drain voltage, the Kondo resonance is destroyed at bias comparable to the Kondo temperature. Calculated scaling coefficients of differential conductance for zero-bias anomaly fall within the range of experimental estimates.

arXiv:1001.4989•Physical Review B•Jun 4, 2010•23 citations

Continuous-time quantum Monte Carlo and maximum entropy approach to an imaginary-time formulation of strongly correlated steady-state transport

Andreas Dirks, Philipp Werner, Mark Jarrell +1 more

Recently, Han and Heary [Phys. Rev. Lett. 99, 236808 (2007)] proposed an approach to steady-state quantum transport through mesoscopic structures, which maps the nonequilibrium problem onto a family of auxiliary quantum impurity systems subject to imaginary voltages. We employ continuous-time quantum Monte-Carlo solvers to calculate accurate imaginary time data for the auxiliary models. The spectral function is obtained from a maximum entropy analytical continuation in both Matsubara frequency and complexified voltage. To enable the analytical continuation we construct a kernel which is compatible with the analytical structure of the theory. While it remains a formidable task to extract reliable spectral functions from this unbiased procedure, particularly for large voltages, our results indicate that the method in principle yields results in agreement with those obtained by other methods.

arXiv:1002.4081•Physical Review E•Aug 5, 2010•30 citations

CT-QMC and Maximum Entropy Approach to a Scattering-States Formulation of Strongly Correlated Steady-State Transport

Andreas Dirks, Philipp Werner, Mark Jarrell +1 more

In a recent work [17] Han and Heary introduced a formalism approaching steady-state quantum transport through mesoscopic structures which maps the non-equilibrium problem onto a family of auxiliary equilibrium quantum impurity systems by introducing imaginary voltages. We apply continuous-time quantum Monte-Carlo solvers to obtain precise and unbiased imaginary-time data for these auxiliary models. Physical observables are obtained by an analytical continuation in both Matsubara frequency and complexified voltage using a Maximum Entropy Method (MEM). The MEM is introduced by means of a kernel function compatible with the analytical structure of the theory. While it remains a yet challenging task to obtain reliable spectral functions, this unbiased procedure seems to indicate that the formalism yields results which are compatible with those of other methods.

Journal of Physics Conference Series•Jan 1, 2011•1 citation

Imaginary-time quantum many-body theory out of equilibrium: Formal equivalence to Keldysh real-time theory and calculation of static properties

Jong E. Han, Andreas Dirks, Thomas Pruschke

We discuss the formal relationship between the real-time Keldysh and imaginary-time theory for nonequilibrium in quantum dot systems. The latter can be reformulated using the recently proposed Matsubara-voltage approach. We establish general conditions for correct analytic continuation procedure on physical observables, and apply the technique to the calculation of static quantities in steady-state nonequilibrium for a quantum dot subject to a finite bias voltage and external magnetic field. Limitations of the Matsubara voltage approach are also pointed out.

arXiv:1205.1816•Physical Review B•Oct 17, 2012•16 citations

Nonequilibrium steady state for strongly correlated many-body systems: Variational cluster approach

Michael Knap, Wolfgang von der Linden, Enrico Arrigoni

A numerical approach is presented that allows to compute nonequilibrium steady-state properties of strongly correlated quantum many-body systems. The method is imbedded in the Keldysh Green's function formalism and is based upon the idea of the variational cluster approach as far as the treatment of strong correlations is concerned. It appears that the variational aspect is crucial as it allows for a suitable optimization of a ``reference'' system to the nonequilibrium target state. The approach is neither perturbative in the many-body interaction nor in the field, that drives the system out of equilibrium, and it allows to study strong perturbations and nonlinear responses of systems in which also the correlated region is spatially extended. We apply the presented approach to nonlinear transport across a strongly correlated quantum wire described by the fermionic Hubbard model. We illustrate how the method bridges to cluster dynamical mean-field theory upon coupling two baths containing an increasing number of uncorrelated sites.

arXiv:1104.3838•Physical Review B•Sep 27, 2011•32 citations

Many-body wavefunctions for quantum impurities out of equilibrium. I. The nonequilibrium Kondo model

Adrian B. Culver, Natan Andrei

We present here the details of a method [A. B. Culver and N. Andrei, Many-body wavefunctions for quantum impurities out of equilibrium, Phys. Rev. B 103, L201103 (2021)] for calculating the time-dependent many-body wavefunction that follows a local quench. We apply the method to the voltage-driven nonequilibrium Kondo model to find the exact time-evolving wavefunction following a quench where the dot is suddenly attached to the leads at $t=0$. The method, which does not use Bethe ansatz, also works in other quantum impurity models and may be of wider applicability. We show that the long-time limit (with the system size taken to infinity first) of the time-evolving wavefunction of the Kondo model is a current-carrying nonequilibrium steady state that satisfies the Lippmann-Schwinger equation. We show that the electric current in the time-evolving wavefunction is given by a series expression that can be expanded either in weak coupling or in strong coupling, converging to all orders in the steady-state limit in either case. The series agrees to leading order with known results in the well-studied regime of weak antiferromagnetic coupling and also reveals a universal regime of strong ferromagnetic coupling with Kondo temperature ${T}_{K}^{(F)}=D{e}^{\ensuremath{-}\frac{3{\ensuremath{\pi}}^{2}}{8}\ensuremath{\rho}|J|}$ ($J<0, \ensuremath{\rho}|J|\ensuremath{\rightarrow}\ensuremath{\infty}$). In this regime, the differential conductance $dI/dV$ reaches the unitarity limit $2{e}^{2}/h$ asymptotically at large voltage or temperature.

arXiv:1912.00281•Physical review. B./Physical review. B•May 4, 2021•5 citations
We extend the imaginary-time formulation of the equilibrium quantum many-body theory to steady-state nonequilibrium with an application to strongly correlated transport. By introducing the Matsubara voltage, we maintain the finite chemical potential shifts in the Fermi-Dirac function, in agreement with the Keldysh formulation. The formulation is applied to strongly correlated transport in the Kondo regime using the quantum Monte Carlo method.
We examine the recently proposed imaginary-time formulation for strongly correlated steady-state nonequilibrium for its range of validity and discuss significant improvements in the analytic continuation of the Matsubara voltage as well as the fermionic Matsubara frequency. The discretization error in the conventional Hirsch-Fye algorithm has been compensated in the Fourier transformation with reliable small frequency behavior of self-energy. Here we give detailed discussions for generalized spectral representation ansatz by including high-order vertex corrections and its numerical analytic continuation procedures. The differential conductance calculations agree accurately with available data from other nonequilibrium transport theories. It is verified that, at finite source-drain voltage, the Kondo resonance is destroyed at bias comparable to the Kondo temperature. Calculated scaling coefficients of differential conductance for zero-bias anomaly fall within the range of experimental estimates.
Recently, Han and Heary [Phys. Rev. Lett. 99, 236808 (2007)] proposed an approach to steady-state quantum transport through mesoscopic structures, which maps the nonequilibrium problem onto a family of auxiliary quantum impurity systems subject to imaginary voltages. We employ continuous-time quantum Monte-Carlo solvers to calculate accurate imaginary time data for the auxiliary models. The spectral function is obtained from a maximum entropy analytical continuation in both Matsubara frequency and complexified voltage. To enable the analytical continuation we construct a kernel which is compatible with the analytical structure of the theory. While it remains a formidable task to extract reliable spectral functions from this unbiased procedure, particularly for large voltages, our results indicate that the method in principle yields results in agreement with those obtained by other methods.
In a recent work [17] Han and Heary introduced a formalism approaching steady-state quantum transport through mesoscopic structures which maps the non-equilibrium problem onto a family of auxiliary equilibrium quantum impurity systems by introducing imaginary voltages. We apply continuous-time quantum Monte-Carlo solvers to obtain precise and unbiased imaginary-time data for these auxiliary models. Physical observables are obtained by an analytical continuation in both Matsubara frequency and complexified voltage using a Maximum Entropy Method (MEM). The MEM is introduced by means of a kernel function compatible with the analytical structure of the theory. While it remains a yet challenging task to obtain reliable spectral functions, this unbiased procedure seems to indicate that the formalism yields results which are compatible with those of other methods.
We discuss the formal relationship between the real-time Keldysh and imaginary-time theory for nonequilibrium in quantum dot systems. The latter can be reformulated using the recently proposed Matsubara-voltage approach. We establish general conditions for correct analytic continuation procedure on physical observables, and apply the technique to the calculation of static quantities in steady-state nonequilibrium for a quantum dot subject to a finite bias voltage and external magnetic field. Limitations of the Matsubara voltage approach are also pointed out.
A numerical approach is presented that allows to compute nonequilibrium steady-state properties of strongly correlated quantum many-body systems. The method is imbedded in the Keldysh Green's function formalism and is based upon the idea of the variational cluster approach as far as the treatment of strong correlations is concerned. It appears that the variational aspect is crucial as it allows for a suitable optimization of a ``reference'' system to the nonequilibrium target state. The approach is neither perturbative in the many-body interaction nor in the field, that drives the system out of equilibrium, and it allows to study strong perturbations and nonlinear responses of systems in which also the correlated region is spatially extended. We apply the presented approach to nonlinear transport across a strongly correlated quantum wire described by the fermionic Hubbard model. We illustrate how the method bridges to cluster dynamical mean-field theory upon coupling two baths containing an increasing number of uncorrelated sites.
We present here the details of a method [A. B. Culver and N. Andrei, Many-body wavefunctions for quantum impurities out of equilibrium, Phys. Rev. B 103, L201103 (2021)] for calculating the time-dependent many-body wavefunction that follows a local quench. We apply the method to the voltage-driven nonequilibrium Kondo model to find the exact time-evolving wavefunction following a quench where the dot is suddenly attached to the leads at t=0t=0t=0. The method, which does not use Bethe ansatz, also works in other quantum impurity models and may be of wider applicability. We show that the long-time limit (with the system size taken to infinity first) of the time-evolving wavefunction of the Kondo model is a current-carrying nonequilibrium steady state that satisfies the Lippmann-Schwinger equation. We show that the electric current in the time-evolving wavefunction is given by a series expression that can be expanded either in weak coupling or in strong coupling, converging to all orders in the steady-state limit in either case. The series agrees to leading order with known results in the well-studied regime of weak antiferromagnetic coupling and also reveals a universal regime of strong ferromagnetic coupling with Kondo temperature ${T}_{K}^{(F)}=D{e}^{\ensuremath{-}\frac{3{\ensuremath{\pi}}^{2}}{8}\ensuremath{\rho}|J|}$ ($J<0, \ensuremath{\rho}|J|\ensuremath{\rightarrow}\ensuremath{\infty}$). In this regime, the differential conductance dI/dVdI/dVdI/dV reaches the unitarity limit 2e2/h2{e}^{2}/h2e2/h asymptotically at large voltage or temperature.