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:cond-mat/0508060 · Physical Review B

Quantum simulation of many-body effects in steady-state nonequilibrium: Electron-phonon coupling in quantum dots

Jong E. Han

Published Mar 15, 2006• 36 citations
arXivHTMLPDF

Abstract

We develope a method of mapping quantum nonequilibrium steady-state to an effective equilibrium system and present an algorithm to calculate electron-transport using an equilibrium technique. A systematic implementation of boundary conditions in steady-state nonequilibrium is made in the statistical operator $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y}$ constructed from scattering state operators. We explicitly demonstrate the equivalence of this method to nonequilibrium Green function techniques for a noninteracting quantum dot model. In electron-phonon coupled quantum dot systems, we formulate an algorithm to construct the statistical bias operator $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y}$ and perform a full many-body calculation with the quantum Monte Carlo technique. The results coherently demonstrate various transport behaviors such as phonon dephasing, $I\ensuremath{-}V$ staircase, and phonon-assisted tunneling phenomena. This formulation makes the existing computational quantum many-body techniques applicable to quantum steady-state nonequilibrium problems, which will complement the theories based on the diagrammatic approach.

DOI
10.1103/physrevb.73.125319

Related papers

Mapping of strongly correlated steady-state nonequilibrium system to an effective equilibrium

Jong E. Han

By mapping steady-state nonequilibrium to an effective equilibrium, we formulate nonequilibrium problems within an equilibrium picture. The Hamiltonian in the open system is rewritten in terms of scattering states with appropriate boundary condition. We first study the analytic properties of many-body scattering states, impose the boundary-condition operator in a statistical operator and prove that this mapping is equivalent to the linear-response theory in the low-bias limit. In an example of infinite-$U$ Anderson impurity model, we approximately solve the scattering state creation operators, based on which we derive the bias operator $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y}$ to construct the nonequilibrium ensemble in the form of the Boltzmann factor ${e}^{\ensuremath{-}\ensuremath{\beta}(\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}\ensuremath{-}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y})}$. The resulting effective Hamiltonian is solved by noncrossing approximation. We obtain the $I\text{\ensuremath{-}}V$ features of Kondo anomaly conductance at zero bias, inelastic transport via the charge excitation on the quantum dot, and significant inelastic current background over a wide range of bias.

arXiv:cond-mat/0604583•Physical Review B•Mar 30, 2007•24 citations

Quantum transport through a quantum dot: Combining the scattering-states numerical renormalization group with nonequilibrium Green functions

Frithjof B. Anders, Sebastian Schmitt

Scattering states fulfill the correct boundary conditions of a current carrying open quantum system. Discretizing the energy continuum of these states allows for employing Wilson's numerical renormalization group approach without violating the boundary conditions by using a finite size system. We evolve the analytically known steady-state density operator for a non-interacting quantum-system at finite bias to the full interacting problem by the time-dependent numerical renormalization group after switching on the local charging energy. Using a newly developed algorithm for steady-state nonequilibrium Green functions, we can calculate the current I as function of bias voltage V for arbitrary temperature and magnetic field. A comparison with second-order and GW Kadanoff-Baym-Keldysh results shows excellent agreement for weak interaction strength U .

Journal of Physics Conference Series•Apr 1, 2010•6 citations

Quantum Monte Carlo Method in the Steady State

André Erpenbeck, Emanuel Gull, Guy Cohen

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.

arXiv:2207.07547•Physical Review Letters•May 2, 2023•41 citations

In and Out-of-Equilibrium Ab Initio Theory of Electrons and Phonons

Gianluca Stefanucci, Robert van Leeuwen, Enrico Perfetto

In this work, we lay down the ab initio many-body quantum theory of electrons and phonons in equilibrium as well as in steady-state or time-varying settings.Our focus is on the harmonic approximation, but the developed tools can readily incorporate anharmonic effects.We begin by showing the necessity of determining the ab initio Hamiltonian in a self-consistent manner to ensure the existence of an equilibrium state.We then identify the correct partitioning into a "noninteracting" and an "interacting" part to carry out diagrammatic expansions in terms of dressed propagators and screened interactions.The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, showcasing the emergence of the coupling between electrons and coherent phonons through the time-local Ehrenfest diagram.The Hedin equations have limited practical utility for real-time simulations of systems driven out of equilibrium by external fields.To overcome this limitation, we leverage the versatility of the diagrammatic approach to generate a closed system of differential equations for the dressed propagators and nuclear displacements.These are the Kadanoff-Baym equations for electrons and phonons.The formalism naturally merges with the theory of conserving approximations, which guarantee the satisfaction of the continuity equation and the conservation of total energy during time evolution.As an example, we show that the popular Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving, provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electronphonon coupling.We also derive the formal solution of the Kadanoff-Baym equations for nonequilibrium steady states, which is useful for studies in photovoltaics and optoelectronics.Interestingly, the expansion of the phononic Green's function around the quasiphonon energies points to a correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.

Physical Review X•Sep 11, 2023•32 citations

Effective Equilibrium Description of Nonequilibrium Quantum Transport II: Perturbation Theory for Interacting Models

Prasenjit Dutt, Jens Koch, Jong E. Han +1 more

The theoretical description of strongly correlated quantum systems out of equilibrium presents several challenges and a number of open questions persist. In this paper we focus on nonlinear electronic transport through an interacting quantum dot maintained at finite bias using a concept introduced by Hershfield [Phys. Rev. Lett. 70, 2134 (1993)] whereby one can express such nonequilibrium quantum impurity models in terms of the system's Lippmann-Schwinger operators. These scattering operators allow one to reformulate the nonequilibrium problem as an effective equilibrium problem associated with a modified Hamiltonian. In this paper we provide a pedagogical analysis of the core concepts of the effective equilibrium theory. First, we demonstrate the equivalence between observables computed using the Schwinger-Keldysh framework and the effective equilibrium approach, and relate the Green's functions in the two theoretical frameworks. Second, we expound some applications of this method in the context of interacting quantum impurity models. We introduce a novel framework to treat effects of interactions perturbatively while capturing the entire dependence on the bias voltage. For the sake of concreteness, we employ the Anderson model as a prototype for this scheme. Working at the particle-hole symmetric point, we investigate the fate of the Abrikosov-Suhl resonance as a function of bias voltage and magnetic field.

arXiv:1101.1526•Jan 7, 2011•0 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
We develope a method of mapping quantum nonequilibrium steady-state to an effective equilibrium system and present an algorithm to calculate electron-transport using an equilibrium technique. A systematic implementation of boundary conditions in steady-state nonequilibrium is made in the statistical operator Y\ifmmode^\elseˆ\fi\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y}Y\ifmmode^\elseˆ\fi​ constructed from scattering state operators. We explicitly demonstrate the equivalence of this method to nonequilibrium Green function techniques for a noninteracting quantum dot model. In electron-phonon coupled quantum dot systems, we formulate an algorithm to construct the statistical bias operator Y\ifmmode^\elseˆ\fi\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y}Y\ifmmode^\elseˆ\fi​ and perform a full many-body calculation with the quantum Monte Carlo technique. The results coherently demonstrate various transport behaviors such as phonon dephasing, I\ensuremath−VI\ensuremath{-}VI\ensuremath−V staircase, and phonon-assisted tunneling phenomena. This formulation makes the existing computational quantum many-body techniques applicable to quantum steady-state nonequilibrium problems, which will complement the theories based on the diagrammatic approach.
By mapping steady-state nonequilibrium to an effective equilibrium, we formulate nonequilibrium problems within an equilibrium picture. The Hamiltonian in the open system is rewritten in terms of scattering states with appropriate boundary condition. We first study the analytic properties of many-body scattering states, impose the boundary-condition operator in a statistical operator and prove that this mapping is equivalent to the linear-response theory in the low-bias limit. In an example of infinite-UUU Anderson impurity model, we approximately solve the scattering state creation operators, based on which we derive the bias operator $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y}$ to construct the nonequilibrium ensemble in the form of the Boltzmann factor ${e}^{\ensuremath{-}\ensuremath{\beta}(\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}\ensuremath{-}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{Y})}$. The resulting effective Hamiltonian is solved by noncrossing approximation. We obtain the $I\text{\ensuremath{-}}V$ features of Kondo anomaly conductance at zero bias, inelastic transport via the charge excitation on the quantum dot, and significant inelastic current background over a wide range of bias.
Scattering states fulfill the correct boundary conditions of a current carrying open quantum system. Discretizing the energy continuum of these states allows for employing Wilson's numerical renormalization group approach without violating the boundary conditions by using a finite size system. We evolve the analytically known steady-state density operator for a non-interacting quantum-system at finite bias to the full interacting problem by the time-dependent numerical renormalization group after switching on the local charging energy. Using a newly developed algorithm for steady-state nonequilibrium Green functions, we can calculate the current I as function of bias voltage V for arbitrary temperature and magnetic field. A comparison with second-order and GW Kadanoff-Baym-Keldysh results shows excellent agreement for weak interaction strength U .
We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.
In this work, we lay down the ab initio many-body quantum theory of electrons and phonons in equilibrium as well as in steady-state or time-varying settings.Our focus is on the harmonic approximation, but the developed tools can readily incorporate anharmonic effects.We begin by showing the necessity of determining the ab initio Hamiltonian in a self-consistent manner to ensure the existence of an equilibrium state.We then identify the correct partitioning into a "noninteracting" and an "interacting" part to carry out diagrammatic expansions in terms of dressed propagators and screened interactions.The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, showcasing the emergence of the coupling between electrons and coherent phonons through the time-local Ehrenfest diagram.The Hedin equations have limited practical utility for real-time simulations of systems driven out of equilibrium by external fields.To overcome this limitation, we leverage the versatility of the diagrammatic approach to generate a closed system of differential equations for the dressed propagators and nuclear displacements.These are the Kadanoff-Baym equations for electrons and phonons.The formalism naturally merges with the theory of conserving approximations, which guarantee the satisfaction of the continuity equation and the conservation of total energy during time evolution.As an example, we show that the popular Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving, provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electronphonon coupling.We also derive the formal solution of the Kadanoff-Baym equations for nonequilibrium steady states, which is useful for studies in photovoltaics and optoelectronics.Interestingly, the expansion of the phononic Green's function around the quasiphonon energies points to a correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.
The theoretical description of strongly correlated quantum systems out of equilibrium presents several challenges and a number of open questions persist. In this paper we focus on nonlinear electronic transport through an interacting quantum dot maintained at finite bias using a concept introduced by Hershfield [Phys. Rev. Lett. 70, 2134 (1993)] whereby one can express such nonequilibrium quantum impurity models in terms of the system's Lippmann-Schwinger operators. These scattering operators allow one to reformulate the nonequilibrium problem as an effective equilibrium problem associated with a modified Hamiltonian. In this paper we provide a pedagogical analysis of the core concepts of the effective equilibrium theory. First, we demonstrate the equivalence between observables computed using the Schwinger-Keldysh framework and the effective equilibrium approach, and relate the Green's functions in the two theoretical frameworks. Second, we expound some applications of this method in the context of interacting quantum impurity models. We introduce a novel framework to treat effects of interactions perturbatively while capturing the entire dependence on the bias voltage. For the sake of concreteness, we employ the Anderson model as a prototype for this scheme. Working at the particle-hole symmetric point, we investigate the fate of the Abrikosov-Suhl resonance as a function of bias voltage and magnetic field.
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.