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:1009.4100 · Physical Review B

Strong-coupling approach to the Mott-Hubbard insulator on a Bethe lattice in dynamical mean-field theory

Daniel Ruhl, Florian Gebhard

Published Jan 19, 2011• 3 citations
arXivHTMLPDF

Abstract

We calculate the Hubbard bands for the half-filled Hubbard model on a Bethe lattice with an infinite coordination number up to and including third order in the inverse Hubbard interaction. We employ the Kato-Takahashi perturbation theory to solve the self-consistency equation of the dynamical mean-field theory analytically for the single-impurity Anderson model in multichain geometry. The weight of the secondary Hubbard sub-bands is of fourth order so that the two-chain geometry is sufficient for our study. Our results for the Mott-Hubbard gap and the density of states of the lower Hubbard band agree very well with those from numerical dynamical density-matrix renormalization group calculations, apart from a resonance contribution at the upper band edge, which cannot be reproduced in low-order perturbation theory.

DOI
10.1103/physrevb.83.035120

Related papers

Dynamical density-matrix renormalization group for the Mott–Hubbard insulator in high dimensions

Satoshi Nishimoto, Florian Gebhard, Eric Jeckelmann

We study the Hubbard model at half band-filling on a Bethe lattice with infinite coordination number in the paramagnetic insulating phase at zero temperature. We use the dynamical mean-field theory (DMFT) mapping to a single-impurity Anderson model with a bath whose properties have to be determined self-consistently. For a controlled and systematic implementation of the self-consistency scheme we use the fixed-energy (FE) approach to the DMFT. In FE-DMFT the onset and the width of the Hubbard bands are adjusted self-consistently but the energies of the bath levels are kept fixed relatively to both band edges during the calculation of self-consistent hybridization strengths between impurity and bath sites. Using the dynamical density-matrix renormalization group method (DDMRG) we calculate the density of states with a resolution ranging from 3% of the bare bandwidth $W=4t$ at high energies to 0.5% in the vicinity of the gap. The DDMRG resolution and accuracy for the density of states and the gap is superior to those obtained with other numerical methods in previous DMFT investigations. We find that the Mott gap closes at a critical coupling $U_{\rm c}/t=4.45 \pm 0.05$. At $U=4.5t$, we observe prominent shoulders near the onset of the Hubbard bands. They are the remainders of the quasi-particle resonance in the metallic phase which appears to split when the gap opens at $U_{\rm c}$.

arXiv:cond-mat/0406666•Journal of Physics Condensed Matter•Sep 21, 2004•61 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

Comparison of computer-algebra strong-coupling perturbation theory and dynamical mean-field theory for the Mott-Hubbard insulator in high dimensions

Martin Paech, W. Apel, Eva Kalinowski +1 more

We present a large-scale combinatorial-diagrammatic computation of high-order contributions to the strong-coupling Kato-Takahashi perturbation series for the Hubbard model in high dimensions. The ground-state energy of the Mott-insulating phase is determined exactly up to the 15th order in $1/U$. The perturbation expansion is extrapolated to infinite order and the critical behavior is determined using the Domb-Sykes method. We compare the perturbative results with two dynamical mean-field theory (DMFT) calculations using a quantum Monte Carlo method and a density-matrix renormalization group method as impurity solvers. The comparison demonstrates the excellent agreement and accuracy of both extrapolated strong-coupling perturbation theory and quantum Monte Carlo based DMFT, even close to the critical coupling where the Mott insulator becomes unstable.

arXiv:1410.6630•Physical Review B•Dec 29, 2014•0 citations

Fourth-Order Perturbation Expansion for Hubbard Model on a Two-Dimensional Square Lattice

Hiroaki Ikeda, Shogo Shinkai, Kōsaku Yamada

We investigate the Hubbard model on a two-dimensional square lattice by the perturbation expansion to the fourth order in the on-site Coulomb repulsion U. Numerically calculating all diagrams up to the fourth order in self-energy, we examine the convergence of perturbation series in the lattice system. We indicate that the coefficient of each order term rapidly decreases as in the impurity Anderson model for T > 0.1t in the half-filled case, but it holds in the doped case even at lower temperatures. Thus, we can expect that the convergence of perturbation expansion in U is very good in a wide parameter region also in the lattice system, except for T < 0.1t in the half-filled case. We next calculate the density of states in the fourth-order perturbation. In the half-filled case, the shape in a moderate correlation regime is quite different from the three peak structure in the second-order perturbation. Remarkable upper and lower Hubbard bands locate at w = +(-)U/2, and a pseudogap appears at the Fermi level w=0. This is considered as the precursor of the Mott-Hubbard antiferromagnetic structure. In the doped case, quasiparticles with very heavy mass are formed at the Fermi level. Thus, we conclude that the fourth-order perturbation theory overall well explain the asymptotic behaviors in a strong correlation regime.

arXiv:0810.1815•Journal of the Physical Society of Japan•May 26, 2008•5 citations

Finite-temperature dynamics of the Mott insulating Hubbard chain

Alberto Nocera, Fabian H. L. Eßler, Adrian E. Feiguin

We study the dynamical response of the half-filled one-dimensional Hubbard model for a range of interaction strengths $U$ and temperatures $T$ by a combination of numerical and analytical techniques. Using time-dependent density matrix renormalization group computations we find that the single-particle spectral function undergoes a crossover to a spin-incoherent Luttinger liquid regime at temperatures $T\ensuremath{\sim}J=4{t}^{2}/U$ for sufficiently large $U>4t$. At smaller values of $U$ and elevated temperatures the spectral function is found to exhibit two thermally broadened bands of excitations, reminiscent of what is found in the Hubbard-I approximation. The dynamical density-density response function is shown to exhibit a finite-temperature resonance at low frequencies inside the Mott gap, with a physical origin similar to the Villain mode in gapped quantum spin chains. We complement our numerical computations by developing an analytic strong-coupling approach to the low-temperature dynamics in the spin-incoherent regime.

arXiv:1710.06452•Physical review. B./Physical review. B•Jan 29, 2018•44 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
We calculate the Hubbard bands for the half-filled Hubbard model on a Bethe lattice with an infinite coordination number up to and including third order in the inverse Hubbard interaction. We employ the Kato-Takahashi perturbation theory to solve the self-consistency equation of the dynamical mean-field theory analytically for the single-impurity Anderson model in multichain geometry. The weight of the secondary Hubbard sub-bands is of fourth order so that the two-chain geometry is sufficient for our study. Our results for the Mott-Hubbard gap and the density of states of the lower Hubbard band agree very well with those from numerical dynamical density-matrix renormalization group calculations, apart from a resonance contribution at the upper band edge, which cannot be reproduced in low-order perturbation theory.
We study the Hubbard model at half band-filling on a Bethe lattice with infinite coordination number in the paramagnetic insulating phase at zero temperature. We use the dynamical mean-field theory (DMFT) mapping to a single-impurity Anderson model with a bath whose properties have to be determined self-consistently. For a controlled and systematic implementation of the self-consistency scheme we use the fixed-energy (FE) approach to the DMFT. In FE-DMFT the onset and the width of the Hubbard bands are adjusted self-consistently but the energies of the bath levels are kept fixed relatively to both band edges during the calculation of self-consistent hybridization strengths between impurity and bath sites. Using the dynamical density-matrix renormalization group method (DDMRG) we calculate the density of states with a resolution ranging from 3% of the bare bandwidth W=4tW=4tW=4t at high energies to 0.5% in the vicinity of the gap. The DDMRG resolution and accuracy for the density of states and the gap is superior to those obtained with other numerical methods in previous DMFT investigations. We find that the Mott gap closes at a critical coupling Uc/t=4.45±0.05U_{\rm c}/t=4.45 \pm 0.05Uc​/t=4.45±0.05. At U=4.5tU=4.5tU=4.5t, we observe prominent shoulders near the onset of the Hubbard bands. They are the remainders of the quasi-particle resonance in the metallic phase which appears to split when the gap opens at UcU_{\rm c}Uc​.
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.
We present a large-scale combinatorial-diagrammatic computation of high-order contributions to the strong-coupling Kato-Takahashi perturbation series for the Hubbard model in high dimensions. The ground-state energy of the Mott-insulating phase is determined exactly up to the 15th order in 1/U1/U1/U. The perturbation expansion is extrapolated to infinite order and the critical behavior is determined using the Domb-Sykes method. We compare the perturbative results with two dynamical mean-field theory (DMFT) calculations using a quantum Monte Carlo method and a density-matrix renormalization group method as impurity solvers. The comparison demonstrates the excellent agreement and accuracy of both extrapolated strong-coupling perturbation theory and quantum Monte Carlo based DMFT, even close to the critical coupling where the Mott insulator becomes unstable.
We investigate the Hubbard model on a two-dimensional square lattice by the perturbation expansion to the fourth order in the on-site Coulomb repulsion U. Numerically calculating all diagrams up to the fourth order in self-energy, we examine the convergence of perturbation series in the lattice system. We indicate that the coefficient of each order term rapidly decreases as in the impurity Anderson model for T > 0.1t in the half-filled case, but it holds in the doped case even at lower temperatures. Thus, we can expect that the convergence of perturbation expansion in U is very good in a wide parameter region also in the lattice system, except for T < 0.1t in the half-filled case. We next calculate the density of states in the fourth-order perturbation. In the half-filled case, the shape in a moderate correlation regime is quite different from the three peak structure in the second-order perturbation. Remarkable upper and lower Hubbard bands locate at w = +(-)U/2, and a pseudogap appears at the Fermi level w=0. This is considered as the precursor of the Mott-Hubbard antiferromagnetic structure. In the doped case, quasiparticles with very heavy mass are formed at the Fermi level. Thus, we conclude that the fourth-order perturbation theory overall well explain the asymptotic behaviors in a strong correlation regime.
We study the dynamical response of the half-filled one-dimensional Hubbard model for a range of interaction strengths UUU and temperatures TTT by a combination of numerical and analytical techniques. Using time-dependent density matrix renormalization group computations we find that the single-particle spectral function undergoes a crossover to a spin-incoherent Luttinger liquid regime at temperatures $T\ensuremath{\sim}J=4{t}^{2}/U$ for sufficiently large U>4tU>4tU>4t. At smaller values of UUU and elevated temperatures the spectral function is found to exhibit two thermally broadened bands of excitations, reminiscent of what is found in the Hubbard-I approximation. The dynamical density-density response function is shown to exhibit a finite-temperature resonance at low frequencies inside the Mott gap, with a physical origin similar to the Villain mode in gapped quantum spin chains. We complement our numerical computations by developing an analytic strong-coupling approach to the low-temperature dynamics in the spin-incoherent regime.
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.