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/0009367 · Physical review. B, Condensed matter

Quantum Monte Carlo calculation of the finite temperature Mott-Hubbard transition

Jaewook Joo, V. S. Oudovenko

Published Oct 19, 2001• 41 citations
arXivHTMLPDF

Abstract

We present clear numerical evidence for the coexistence of metallic and insulating dynamical mean-field theory (DMFT) solutions in a half-filled single-band Hubbard model at finite temperature. The quantum Monte Carlo (QMC) method is used to solve the DMFT equation. We discuss important technical aspects of the DMFT-QMC method, which should be properly dealt with in order to obtain reliable DMFT solutions at low temperatures. Among these technical aspects are the critical slowing down near the boundaries of the coexistence region, the convergence criteria for the DMFT iterations, and the correct interpolation of the discretized Green's functions.

DOI
10.1103/physrevb.64.193102

Related papers

Coexistence of solutions in dynamical mean-field theory of the Mott transition

Werner Krauth

In this paper, I discuss the finite-temperature metal-insulator transition of the paramagnetic Hubbard model within dynamical mean-field theory. I show that coexisting solutions, the hallmark of such a transition, can be obtained in a consistent way both from quantum Monte Carlo (QMC) simulations and from the exact diagonalization method. I pay special attention to discretization errors within the QMC method. These errors explain why it is difficult to obtain the solutions from the QMC method close to the boundaries of the coexistence region.

arXiv:cond-mat/9908221•Physical review. B, Condensed matter•Sep 15, 2000•14 citations

Dynamical mean-field study of the Mott transition in the half-filled Hubbard model on a triangular lattice

Karan Aryanpour, Warren E. Pickett, Richard T. Scalettar

We employ dynamical mean-field theory (DMFT) with a quantum Monte Carlo atomic solver to investigate the finite-temperature Mott transition in the Hubbard model with nearest-neighbor hopping on a triangular lattice at half-filling. We estimate the value of the critical interaction to be ${U}_{c}=12.0\ifmmode\pm\else\textpm\fi{}0.5$ in units of the hopping amplitude $t$ through the evolution of the magnetic moment, spectral function, internal energy, and specific heat as the interaction $U$ and temperature $T$ are varied. This work also presents a comparison between DMFT and finite-size determinant quantum Monte Carlo calculations and a discussion of the advantages and limitations of both methods.

arXiv:cond-mat/0604609•Physical Review B•Aug 30, 2006•54 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

Mott-Hubbard metal-insulator transition at noninteger filling

Krzysztof Byczuk, Walter Hofstetter, D. Vollhardt

Correlated electrons in a binary alloy ${A}_{x}{B}_{1\ensuremath{-}x}$ are investigated within the Hubbard model and dynamical mean-field theory (DMFT). The random energies ${\ensuremath{\epsilon}}_{i}$ have a bimodal probability distribution and an energy separation $\ensuremath{\Delta}.$ We solve the DMFT equations by the numerical renormalization-group method at zero temperature and calculate the spectral density as a function of disorder strength $\ensuremath{\Delta}$ and interaction U at different fillings. For filling factors $\ensuremath{\nu}=x$ or $1+x$ the lower or upper alloy subband is half filled and the system becomes a Mott insulator at strong interactions, with a correlation gap at the Fermi level. At the metal-insulator transition hysteresis is observed. We also analyze the effective theory in the $\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\Delta}}\ensuremath{\infty}$ limit and find good agreement between analytical and numerical results for the critical interaction ${U}_{c}$ at which the metal-insulator transition occurs.

arXiv:cond-mat/0309067•Physical Review B•Jan 29, 2004•47 citations

The numerical calculation of Green's function in Hubbard model

Song Shu

From the Hubbard model,we calculated the Green's function of the system with the dynamical mean-field theory and the quantum Monte Carlo simulation.We also had a thorough discussion on the numerical results of the Green's function in paramagnetic and antiferromagnetic orders and gave the variations of Green's function at different temperatures and different Hubbard energies.From the analysis of these results,we concluded that with the decrease of the temperature or the increase of the Hubbard energy,the system would have Mott transition from metal to insulator.

Journal of Hubei University•Jan 1, 2012•0 citations

Quasi-continuous-time impurity solver for the dynamical mean-field theory with linear scaling in the inverse temperature

D. Rost, Fakher F. Assaad, Nils Blümer

We present an algorithm for solving the self-consistency equations of the dynamical mean-field theory (DMFT) with high precision and efficiency at low temperatures. In each DMFT iteration, the impurity problem is mapped to an auxiliary Hamiltonian, for which the Green function is computed by combining determinantal quantum Monte Carlo (BSS-QMC) calculations with a multigrid extrapolation procedure. The method is numerically exact, i.e., yields results which are free of significant Trotter errors, but retains the BSS advantage, compared to direct QMC impurity solvers, of linear (instead of cubic) scaling with the inverse temperature. The new algorithm is applied to the half-filled Hubbard model close to the Mott transition; detailed comparisons with exact diagonalization, Hirsch-Fye QMC, and continuous-time QMC are provided.

arXiv:1303.2004•Physical Review E•May 29, 2013•16 citations
We present clear numerical evidence for the coexistence of metallic and insulating dynamical mean-field theory (DMFT) solutions in a half-filled single-band Hubbard model at finite temperature. The quantum Monte Carlo (QMC) method is used to solve the DMFT equation. We discuss important technical aspects of the DMFT-QMC method, which should be properly dealt with in order to obtain reliable DMFT solutions at low temperatures. Among these technical aspects are the critical slowing down near the boundaries of the coexistence region, the convergence criteria for the DMFT iterations, and the correct interpolation of the discretized Green's functions.
In this paper, I discuss the finite-temperature metal-insulator transition of the paramagnetic Hubbard model within dynamical mean-field theory. I show that coexisting solutions, the hallmark of such a transition, can be obtained in a consistent way both from quantum Monte Carlo (QMC) simulations and from the exact diagonalization method. I pay special attention to discretization errors within the QMC method. These errors explain why it is difficult to obtain the solutions from the QMC method close to the boundaries of the coexistence region.
We employ dynamical mean-field theory (DMFT) with a quantum Monte Carlo atomic solver to investigate the finite-temperature Mott transition in the Hubbard model with nearest-neighbor hopping on a triangular lattice at half-filling. We estimate the value of the critical interaction to be ${U}_{c}=12.0\ifmmode\pm\else\textpm\fi{}0.5$ in units of the hopping amplitude ttt through the evolution of the magnetic moment, spectral function, internal energy, and specific heat as the interaction UUU and temperature TTT are varied. This work also presents a comparison between DMFT and finite-size determinant quantum Monte Carlo calculations and a discussion of the advantages and limitations of both methods.
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.
Correlated electrons in a binary alloy ${A}_{x}{B}_{1\ensuremath{-}x}$ are investigated within the Hubbard model and dynamical mean-field theory (DMFT). The random energies ${\ensuremath{\epsilon}}_{i}$ have a bimodal probability distribution and an energy separation $\ensuremath{\Delta}.$ We solve the DMFT equations by the numerical renormalization-group method at zero temperature and calculate the spectral density as a function of disorder strength $\ensuremath{\Delta}$ and interaction U at different fillings. For filling factors $\ensuremath{\nu}=x$ or 1+x1+x1+x the lower or upper alloy subband is half filled and the system becomes a Mott insulator at strong interactions, with a correlation gap at the Fermi level. At the metal-insulator transition hysteresis is observed. We also analyze the effective theory in the $\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\Delta}}\ensuremath{\infty}$ limit and find good agreement between analytical and numerical results for the critical interaction Uc{U}_{c}Uc​ at which the metal-insulator transition occurs.
From the Hubbard model,we calculated the Green's function of the system with the dynamical mean-field theory and the quantum Monte Carlo simulation.We also had a thorough discussion on the numerical results of the Green's function in paramagnetic and antiferromagnetic orders and gave the variations of Green's function at different temperatures and different Hubbard energies.From the analysis of these results,we concluded that with the decrease of the temperature or the increase of the Hubbard energy,the system would have Mott transition from metal to insulator.
We present an algorithm for solving the self-consistency equations of the dynamical mean-field theory (DMFT) with high precision and efficiency at low temperatures. In each DMFT iteration, the impurity problem is mapped to an auxiliary Hamiltonian, for which the Green function is computed by combining determinantal quantum Monte Carlo (BSS-QMC) calculations with a multigrid extrapolation procedure. The method is numerically exact, i.e., yields results which are free of significant Trotter errors, but retains the BSS advantage, compared to direct QMC impurity solvers, of linear (instead of cubic) scaling with the inverse temperature. The new algorithm is applied to the half-filled Hubbard model close to the Mott transition; detailed comparisons with exact diagonalization, Hirsch-Fye QMC, and continuous-time QMC are provided.