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

Steady-State Nonequilibrium Density of States of Driven Strongly Correlated Lattice Models in Infinite Dimensions

A. V. Joura, James K. Freericks, Thomas Pruschke

Published Nov 3, 2008• 80 citations
arXivHTMLPDF

Abstract

An exact formalism for calculating the retarded and advanced Green's functions of strongly correlated lattice models in a uniform electric field is derived within dynamical mean-field theory. To illustrate the method, we solve for the nonequilibrium density of states of the Hubbard model in both the metallic and Mott-insulating phases at half-filling (with an arbitrary strength electric field) by employing the approximate numerical renormalization group as the impurity solver. This general approach can be applied to any strongly correlated lattice model in the limit of large dimensions.

DOI
10.1103/physrevlett.101.196401

Related papers

Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions

Antoine Georges, Gabriel Kotliar, Werner Krauth +1 more

We review the dynamical mean-field theory of strongly correlated electron systems which is based on a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. This mapping is exact for models of correlated electrons in the limit of large lattice coordination (or infinite spatial dimensions). It extends the standard mean-field construction from classical statistical mechanics to quantum problems. We discuss the physical ideas underlying this theory and its mathematical derivation. Various analytic and numerical techniques that have been developed recently in order to analyze and solve the dynamical mean-field equations are reviewed and compared to each other. The method can be used for the determination of phase diagrams (by comparing the stability of various types of long-range order), and the calculation of thermodynamic properties, one-particle Green's functions, and response functions. We review in detail the recent progress in understanding the Hubbard model and the Mott metal-insulator transition within this approach, including some comparison to experiments on three-dimensional transition-metal oxides. We present an overview of the rapidly developing field of applications of this method to other systems. The present limitations of the approach, and possible extensions of the formalism are finally discussed. Computer programs for the numerical implementation of this method are also provided with this article.

Reviews of Modern Physics•Jan 1, 1996•6,900 citations

The Local Impurity Self Consistent Approximation (LISA) to Strongly Correlated Fermion Systems and the Limit of Infinite Dimensions

Antoine Georges, Gabriel Kotliar, Werner Krauth +1 more

We review the dynamical mean-field theory of strongly correlated fermion systems which is based on a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. This mapping is exact in the limit of large lattice coordination (or infinite spatial dimensions). This method can be used for the determination of phase diagrams and the calculation of thermodynamic properties, one-particle Green's functions, and response functions, using analytic and numerical techniques which are described. We review the recent progress in understanding the Hubbard model and the Mott metal-insulator transition within this approach, including some comparison to experiments on three-dimensional transition-metal oxides. Some applications of this method to other models are also reviewed. The present limitations of the approach, and possible extensions of the formalism are finally discussed. Computer programs for the numerical implementation of this method are provided with this article.

arXiv:cond-mat/9510091•Oct 17, 1995•3 citations

Mott-Hubbard Insulator in Infinite Dimensions

Eva Kalinowski, Florian Gebhard

We calculate the one-particle density of states for the Mott-Hubbard insulating phase of the Hubbard model on a Bethe lattice in the limit of infinite coordination number. We employ the Kato-Takahashi perturbation theory around the strong-coupling limit to derive the Green function. We show that the Green function for the lower Hubbard band can be expressed in terms of polynomials in the bare hole-hopping operator. We check our technique against the exact solution of the Falicov-Kimball model and give explicit results up to and including second order in the inverse Hubbard interaction. Our results provide a stringent test for analytical and numerical investigations of the Mott-Hubbard insulator and the Mott-Hubbard transition within the dynamical mean-field theory. We find that the Hubbard-III approximation is not satisfactory beyond lowest order, but the local-moment approach provides a very good description of the Mott-Hubbard insulator at strong coupling.

arXiv:cond-mat/0109535•Sep 28, 2001•0 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

Transport through a correlated interface: Auxiliary master equation approach

Irakli Titvinidze, Antonius Dorda, Wolfgang von der Linden +1 more

We present improvements of a recently introduced numerical method [E. Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)] to compute steady-state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a nonequilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the nonequilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, ${N}_{b}$ uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of ${N}_{b}$ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semiquantitative description of the Kondo regime. To illustrate the approach, we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The nonequilibrium situation is driven by a bias voltage applied to the leads. For this system, we investigate the spectral function and the steady-state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular, we investigate the nonequilibrium behavior of quasiparticle excitations within the Mott gap of the correlated layer. We find for low-bias voltage Kondo-type behavior in the vicinity of the insulating phase. In particular, we observe a splitting of the Kondo resonance as a function of the bias voltage.

arXiv:1508.02953•Physical Review B•Dec 21, 2015•24 citations

Quenching Bloch oscillations in a strongly correlated material: Nonequilibrium dynamical mean-field theory

James K. Freericks

Dynamical mean-field theory is generalized to solve the nonequilibrium Keldysh boundary problem: a system is started in equilibrium at a temperature $T=0.1$, a uniform electric field is turned on at $t=0$, and the system is monitored as it approaches the steady state. The focus here is on the Bloch oscillations of the current and how they decay after their initial appearance near $t=0$. The system is evolved out to the longest time allowed by our computational resources---in most cases we are unable to reach the steady state. The strongly correlated material is modeled by the spinless Falicov-Kimball model at half-filling on a hypercubic lattice in $d=\ensuremath{\infty}$ dimensions, which has a Mott-like metal-insulator transition at $U=\sqrt{2}$. The computational algorithm employed is highly efficient, parallelizes well, and scales to thousands of processors. For strong fields, we find beats develop with a period of $2\ensuremath{\pi}∕U$, while for strong interactions, the Bloch oscillations are sharply damped and become quite irregular in time.

arXiv:0711.0999•Physical Review B•Feb 12, 2008•115 citations
An exact formalism for calculating the retarded and advanced Green's functions of strongly correlated lattice models in a uniform electric field is derived within dynamical mean-field theory. To illustrate the method, we solve for the nonequilibrium density of states of the Hubbard model in both the metallic and Mott-insulating phases at half-filling (with an arbitrary strength electric field) by employing the approximate numerical renormalization group as the impurity solver. This general approach can be applied to any strongly correlated lattice model in the limit of large dimensions.
We review the dynamical mean-field theory of strongly correlated electron systems which is based on a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. This mapping is exact for models of correlated electrons in the limit of large lattice coordination (or infinite spatial dimensions). It extends the standard mean-field construction from classical statistical mechanics to quantum problems. We discuss the physical ideas underlying this theory and its mathematical derivation. Various analytic and numerical techniques that have been developed recently in order to analyze and solve the dynamical mean-field equations are reviewed and compared to each other. The method can be used for the determination of phase diagrams (by comparing the stability of various types of long-range order), and the calculation of thermodynamic properties, one-particle Green's functions, and response functions. We review in detail the recent progress in understanding the Hubbard model and the Mott metal-insulator transition within this approach, including some comparison to experiments on three-dimensional transition-metal oxides. We present an overview of the rapidly developing field of applications of this method to other systems. The present limitations of the approach, and possible extensions of the formalism are finally discussed. Computer programs for the numerical implementation of this method are also provided with this article.
We review the dynamical mean-field theory of strongly correlated fermion systems which is based on a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. This mapping is exact in the limit of large lattice coordination (or infinite spatial dimensions). This method can be used for the determination of phase diagrams and the calculation of thermodynamic properties, one-particle Green's functions, and response functions, using analytic and numerical techniques which are described. We review the recent progress in understanding the Hubbard model and the Mott metal-insulator transition within this approach, including some comparison to experiments on three-dimensional transition-metal oxides. Some applications of this method to other models are also reviewed. The present limitations of the approach, and possible extensions of the formalism are finally discussed. Computer programs for the numerical implementation of this method are provided with this article.
We calculate the one-particle density of states for the Mott-Hubbard insulating phase of the Hubbard model on a Bethe lattice in the limit of infinite coordination number. We employ the Kato-Takahashi perturbation theory around the strong-coupling limit to derive the Green function. We show that the Green function for the lower Hubbard band can be expressed in terms of polynomials in the bare hole-hopping operator. We check our technique against the exact solution of the Falicov-Kimball model and give explicit results up to and including second order in the inverse Hubbard interaction. Our results provide a stringent test for analytical and numerical investigations of the Mott-Hubbard insulator and the Mott-Hubbard transition within the dynamical mean-field theory. We find that the Hubbard-III approximation is not satisfactory beyond lowest order, but the local-moment approach provides a very good description of the Mott-Hubbard insulator at strong coupling.
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.
We present improvements of a recently introduced numerical method [E. Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)] to compute steady-state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a nonequilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the nonequilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, Nb{N}_{b}Nb​ uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of Nb{N}_{b}Nb​ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semiquantitative description of the Kondo regime. To illustrate the approach, we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The nonequilibrium situation is driven by a bias voltage applied to the leads. For this system, we investigate the spectral function and the steady-state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular, we investigate the nonequilibrium behavior of quasiparticle excitations within the Mott gap of the correlated layer. We find for low-bias voltage Kondo-type behavior in the vicinity of the insulating phase. In particular, we observe a splitting of the Kondo resonance as a function of the bias voltage.
Dynamical mean-field theory is generalized to solve the nonequilibrium Keldysh boundary problem: a system is started in equilibrium at a temperature T=0.1T=0.1T=0.1, a uniform electric field is turned on at t=0t=0t=0, and the system is monitored as it approaches the steady state. The focus here is on the Bloch oscillations of the current and how they decay after their initial appearance near t=0t=0t=0. The system is evolved out to the longest time allowed by our computational resources---in most cases we are unable to reach the steady state. The strongly correlated material is modeled by the spinless Falicov-Kimball model at half-filling on a hypercubic lattice in $d=\ensuremath{\infty}$ dimensions, which has a Mott-like metal-insulator transition at U=2U=\sqrt{2}U=2​. The computational algorithm employed is highly efficient, parallelizes well, and scales to thousands of processors. For strong fields, we find beats develop with a period of $2\ensuremath{\pi}∕U$, while for strong interactions, the Bloch oscillations are sharply damped and become quite irregular in time.