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

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, Eric Jeckelmann

Published Dec 29, 2014• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.90.245147

Related papers

Strong-coupling solver for the quantum impurity model

Xi Dai, Kristjan Haule, Gabriel Kotliar

We propose a fast impurity solver for the general quantum impurity model based on the perturbation theory around the atomic limit, which can be used in combination with the local density approximation (LDA) and the dynamical mean-field theory (DMFT). We benchmark the solver in the two-band Hubbard model within DMFT against quantum Monte Carlo (QMC) and numerical renormalization-group (NRG) results. We find that the solver works very well in the paramagnetic Mott insulator phase. We also apply this impurity solver to the DMFT study of the antiferromagnetic phase transition in the unfrustrated Bethe lattice. The Neel temperature obtained by the fast impurity solver agrees very well with the QMC results in the large Hubbard U limit. The method is a promising tool to be used in combination with the $\mathrm{LDA}+\mathrm{DMFT}$ to study Mott insulators starting from first principles.

arXiv:cond-mat/0412505•Physical Review B•Jul 8, 2005•26 citations

Mott insulator: Tenth-order perturbation theory extended to infinite order using a quantum Monte Carlo scheme

Nils Blümer, Eva Kalinowski

We present a method based on the combination of analytical and numerical techniques within the framework of the dynamical mean-field theory. Building upon numerically exact results obtained in an improved quantum Monte Carlo scheme, tenth-order strong-coupling perturbation theory for the Hubbard model on the Bethe lattice is extrapolated to infinite order. We obtain continuous estimates of energy $E$ and double occupancy $D$ with unprecedented precision $\mathcal{O}({10}^{\ensuremath{-}5})$ for the Mott insulator above its stability edge ${U}_{c1}\ensuremath{\approx}4.78$ as well as critical exponents. The relevance for recent experiments on Cr-doped ${\mathrm{V}}_{2}{\mathrm{O}}_{3}$ is pointed out.

arXiv:cond-mat/0404568•Physical Review B•May 4, 2005•21 citations

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

Daniel Ruhl, Florian Gebhard

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.

arXiv:1009.4100•Physical Review B•Jan 19, 2011•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

Fast impurity solver for dynamical mean field theory based on second order perturbation around the atomic limit

Zhuang Jia-Ning, Liu Qing-Mei, Fang Zhong +1 more

This paper proposes an impurity solver for the dynamical mean field theory (DMFT) study of the Mott insulators, which is based on the second order perturbation of the hybridization function. After careful benchmarking with quantum Monte Carlo results on the anti-ferromagnetic phase of the Hubbard model, it concludes that this impurity solver can capture the main physical features in the strong coupling regime and can be a very useful tool for the LDA (local density approximation) + DMFT studies of the Mott insulators with long range order.

arXiv:0911.1422•Chinese Physics B•Aug 1, 2010•8 citations

Strong-coupling expansion for the momentum distribution of the Bose-Hubbard model with benchmarking against exact numerical results

James K. Freericks, H. R. Krishnamurthy, Yasuyuki Kato +2 more

A strong-coupling expansion for the Green's functions, self-energies, and correlation functions of the Bose-Hubbard model is developed. We illustrate the general formalism, which includes all possible (normal-phase) inhomogeneous effects in the formalism, such as disorder or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator--to--superfluid transition along with a generalization of the random-phase-approximation-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions. The accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott-phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling-theory results are benchmarked against numerically exact quantum Monte Carlo simulations in two and three dimensions and against density-matrix renormalization-group calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.

arXiv:0902.3435•Physical Review A•May 20, 2009•84 citations
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 propose a fast impurity solver for the general quantum impurity model based on the perturbation theory around the atomic limit, which can be used in combination with the local density approximation (LDA) and the dynamical mean-field theory (DMFT). We benchmark the solver in the two-band Hubbard model within DMFT against quantum Monte Carlo (QMC) and numerical renormalization-group (NRG) results. We find that the solver works very well in the paramagnetic Mott insulator phase. We also apply this impurity solver to the DMFT study of the antiferromagnetic phase transition in the unfrustrated Bethe lattice. The Neel temperature obtained by the fast impurity solver agrees very well with the QMC results in the large Hubbard U limit. The method is a promising tool to be used in combination with the LDA+DMFT\mathrm{LDA}+\mathrm{DMFT}LDA+DMFT to study Mott insulators starting from first principles.
We present a method based on the combination of analytical and numerical techniques within the framework of the dynamical mean-field theory. Building upon numerically exact results obtained in an improved quantum Monte Carlo scheme, tenth-order strong-coupling perturbation theory for the Hubbard model on the Bethe lattice is extrapolated to infinite order. We obtain continuous estimates of energy EEE and double occupancy DDD with unprecedented precision $\mathcal{O}({10}^{\ensuremath{-}5})$ for the Mott insulator above its stability edge ${U}_{c1}\ensuremath{\approx}4.78$ as well as critical exponents. The relevance for recent experiments on Cr-doped V2O3{\mathrm{V}}_{2}{\mathrm{O}}_{3}V2​O3​ is pointed out.
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 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 paper proposes an impurity solver for the dynamical mean field theory (DMFT) study of the Mott insulators, which is based on the second order perturbation of the hybridization function. After careful benchmarking with quantum Monte Carlo results on the anti-ferromagnetic phase of the Hubbard model, it concludes that this impurity solver can capture the main physical features in the strong coupling regime and can be a very useful tool for the LDA (local density approximation) + DMFT studies of the Mott insulators with long range order.
A strong-coupling expansion for the Green's functions, self-energies, and correlation functions of the Bose-Hubbard model is developed. We illustrate the general formalism, which includes all possible (normal-phase) inhomogeneous effects in the formalism, such as disorder or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator--to--superfluid transition along with a generalization of the random-phase-approximation-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions. The accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott-phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling-theory results are benchmarked against numerically exact quantum Monte Carlo simulations in two and three dimensions and against density-matrix renormalization-group calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.