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/0312390 · Physical Review B

Dynamical mean-field theory of the Hubbard-Holstein model at half filling: Zero temperature metal-insulator and insulator-insulator transitions

Gun Sang Jeon, Tae-Ho Park, Jung Hoon Han, Hyun‐Woo Lee, Han-Yong Choi

Published Sep 29, 2004• 49 citations
arXivHTMLPDF

Abstract

We study the Hubbard-Holstein model which includes both the electron-electron and electron-phonon interactions, characterized by $U$ and $g$, respectively. The model is solved with $U$ and $g$ on an equal footing in the infinite dimensions by employing the dynamical mean-field theory in combination with Wilson's numerical renormalization group. A zero temperature phase diagram of symmetry unbroken states at half filling is mapped out which exhibits the interplay between the two kinds of interactions and combines the two separately studied interaction-driven metal-insulator transitions of the Hubbard and Holstein models within a single frame. The ground state is metallic when both $U$ and $g$ are small, but is insulating when $U$ or $g$ is large, referred to as, respectively, Mott-Hubbard insulator (MHI) and bipolaron insulator (BPI). As the phase boundary between the metallic and MHI (BPI) states is approached from the metallic region, the quasiparticle weight $z$ goes to 0 continuously (discontinuously). Moreover, the two insulating states are distinct and cannot be adiabatically connected, and there is a first order phase transition between them.

DOI
10.1103/physrevb.70.125114

Related papers

First- and second-order phase transitions in the Holstein-Hubbard model

Winfried Koller, D. Meyer, Yūichi Ono +1 more

We investigate metal-insulator transitions in the half-filled Holstein-Hubbard model as a function of the on-site electron-electron interaction U and the electron-phonon coupling g . We use several different numerical methods to calculate the phase diagram, the results of which are in excellent agreement. When the electron-electron interaction U is dominant, the transition is to a Mott insulator; when the electron-phonon interaction dominates, the transition is to a localized bipolaronic state. In the former case, the transition is always found to be second order. This is in contrast to the transition to the bipolaronic state, which is clearly first order for larger values of U . We also present results for the quasiparticle weight and the double occupancy as functions of U and g .

arXiv:cond-mat/0312367•Europhysics Letters (EPL)•May 1, 2004•87 citations

First- and Second Order Phase Transitions in the Holstein-Hubbard Mode

W. Koller, D. Meyer, Ono, Y. +1 more

We investigate metal-insulator transitions in the half-filled Holstein-Hubbard model as a function of the on-site electron-electron interaction U and the electron-phonon coupling g. We use several different numerical methods to calculate the phase diagram, the results of which are in excellent agreement. When the electron-electron interaction U is dominant the transition is to a Mott-insulator; when the electron-phonon interaction dominates, the transition is to a localised bipolaronic state. In the former case, the transition is always found to be second order. This is in contrast to the transition to the bipolaronic state, which is clearly first order for larger values of U. We also present results for the quasiparticle weight and the double-occupancy as function of U and g.

•May 1, 2004•0 citations

Isotope effects in the Hubbard-Holstein model within dynamical mean-field theory

P. Paci, Massimo Capone, E. Cappelluti +2 more

We study the isotope effects arising from the coupling of correlated electrons with dispersionless phonons by considering the Hubbard-Holstein model at half-filling within the dynamical mean-field theory. In particular we calculate the isotope effects on the quasiparticle spectral weight $Z$, the renormalized phonon frequency, and the static charge and spin susceptibilities. In the weakly correlated regime $U∕t\ensuremath{\lesssim}1.5$, where $U$ is the Hubbard repulsion and $t$ is the bare electron half-bandwidth, the physical properties are qualitatively similar to those characterizing the Holstein model in the absence of Coulomb repulsion, where the bipolaronic binding takes place at large electron-phonon coupling and it is reflected in divergent isotope responses. On the contrary in the strongly correlated regime $U∕t\ensuremath{\gtrsim}1.5$, where the bipolaronic metal-insulator transition becomes of first order, the isotope effects are bounded, suggesting that the first-order transition is likely driven by an electronic mechanism, rather then by a lattice instability. These results point out how the isotope responses are extremely sensitive to phase boundaries and they may be used to characterize the competition between the electron-phonon coupling and the Hubbard repulsion.

arXiv:cond-mat/0609630•Physical Review B•Nov 9, 2006•18 citations

Dynamical mean field theory of polarons and bipolarons in the half-filled Holstein model

Massimo Capone, P. Carta, Sergio Ciuchi

The evolution of the properties of a finite density electronic system as a function of the electron-phonon coupling is investigated in the Holstein model using the dynamical mean-field theory (DMFT) that becomes exact in infinite dimensions. We compare the spinless fermion case, in which only isolated polarons can be formed, with the spinful model in which the polarons can bind and form bipolarons. In the latter case, the bipolaronic binding occurs within DMFT as a metal-insulator transition. In the adiabatic regime in which the phonon energy is small with respect to the electron hopping we compare numerically exact DMFT results with an analytical scheme inspired by the Born-Oppenheimer procedure. Within the latter approach, a truncation of the phononic Hilbert space leads to a mapping of the original model onto an Anderson spin-fermion model. In the anti-adiabatic regime (where the phonon energy exceeds the electronic scales) the standard treatment based on Lang-Firsov canonical transformation allows one to map the original model on to an attractive Hubbard model in the spinful case. The separate analysis of the two regimes supports the numerical evidence that the presence of well-defined polaronic lattice displacements is not necessarily associated to a metal-insulator transition, which is instead due to pairing between the carriers. The finite-dimensionality effects neglected in DMFT may lead to a finite conductivity in the bipolaronic state which is, however, not always associated with polaronic distortions. At the polaron crossover the Born-Oppenheimer approximation is shown to break down due to the entanglement of the electron-phonon state.

arXiv:cond-mat/0509542•Physical Review B•Jul 10, 2006•38 citations

Metallicity in the half-filled Holstein-Hubbard model

Holger Fehske, Georg Hager, Eric Jeckelmann

We re-examine the Peierls insulator to Mott insulator transition scenario in the one-dimensional Holstein-Hubbard model where, at half-filling, electron-phonon and electron-electron interactions compete for establishing charge- and spin-density-wave states, respectively. By means of large-scale density-matrix renormalization group calculations we determine the spin, single-particle and two-particle excitation gaps and prove —in the course of a careful finite-size scaling analysis— recent claims for an intervening metallic phase in the weak-coupling regime. We show that for large phonon frequencies the metallic region is even more extended than previously expected, and subdivided into ordinary Luttinger-liquid and bipolaronic-liquid phases.

arXiv:0808.1675•Europhysics Letters (EPL)•Nov 24, 2008•80 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 study the Hubbard-Holstein model which includes both the electron-electron and electron-phonon interactions, characterized by UUU and ggg, respectively. The model is solved with UUU and ggg on an equal footing in the infinite dimensions by employing the dynamical mean-field theory in combination with Wilson's numerical renormalization group. A zero temperature phase diagram of symmetry unbroken states at half filling is mapped out which exhibits the interplay between the two kinds of interactions and combines the two separately studied interaction-driven metal-insulator transitions of the Hubbard and Holstein models within a single frame. The ground state is metallic when both UUU and ggg are small, but is insulating when UUU or ggg is large, referred to as, respectively, Mott-Hubbard insulator (MHI) and bipolaron insulator (BPI). As the phase boundary between the metallic and MHI (BPI) states is approached from the metallic region, the quasiparticle weight zzz goes to 0 continuously (discontinuously). Moreover, the two insulating states are distinct and cannot be adiabatically connected, and there is a first order phase transition between them.
We investigate metal-insulator transitions in the half-filled Holstein-Hubbard model as a function of the on-site electron-electron interaction U and the electron-phonon coupling g . We use several different numerical methods to calculate the phase diagram, the results of which are in excellent agreement. When the electron-electron interaction U is dominant, the transition is to a Mott insulator; when the electron-phonon interaction dominates, the transition is to a localized bipolaronic state. In the former case, the transition is always found to be second order. This is in contrast to the transition to the bipolaronic state, which is clearly first order for larger values of U . We also present results for the quasiparticle weight and the double occupancy as functions of U and g .
We investigate metal-insulator transitions in the half-filled Holstein-Hubbard model as a function of the on-site electron-electron interaction U and the electron-phonon coupling g. We use several different numerical methods to calculate the phase diagram, the results of which are in excellent agreement. When the electron-electron interaction U is dominant the transition is to a Mott-insulator; when the electron-phonon interaction dominates, the transition is to a localised bipolaronic state. In the former case, the transition is always found to be second order. This is in contrast to the transition to the bipolaronic state, which is clearly first order for larger values of U. We also present results for the quasiparticle weight and the double-occupancy as function of U and g.
We study the isotope effects arising from the coupling of correlated electrons with dispersionless phonons by considering the Hubbard-Holstein model at half-filling within the dynamical mean-field theory. In particular we calculate the isotope effects on the quasiparticle spectral weight ZZZ, the renormalized phonon frequency, and the static charge and spin susceptibilities. In the weakly correlated regime $U∕t\ensuremath{\lesssim}1.5$, where UUU is the Hubbard repulsion and ttt is the bare electron half-bandwidth, the physical properties are qualitatively similar to those characterizing the Holstein model in the absence of Coulomb repulsion, where the bipolaronic binding takes place at large electron-phonon coupling and it is reflected in divergent isotope responses. On the contrary in the strongly correlated regime $U∕t\ensuremath{\gtrsim}1.5$, where the bipolaronic metal-insulator transition becomes of first order, the isotope effects are bounded, suggesting that the first-order transition is likely driven by an electronic mechanism, rather then by a lattice instability. These results point out how the isotope responses are extremely sensitive to phase boundaries and they may be used to characterize the competition between the electron-phonon coupling and the Hubbard repulsion.
The evolution of the properties of a finite density electronic system as a function of the electron-phonon coupling is investigated in the Holstein model using the dynamical mean-field theory (DMFT) that becomes exact in infinite dimensions. We compare the spinless fermion case, in which only isolated polarons can be formed, with the spinful model in which the polarons can bind and form bipolarons. In the latter case, the bipolaronic binding occurs within DMFT as a metal-insulator transition. In the adiabatic regime in which the phonon energy is small with respect to the electron hopping we compare numerically exact DMFT results with an analytical scheme inspired by the Born-Oppenheimer procedure. Within the latter approach, a truncation of the phononic Hilbert space leads to a mapping of the original model onto an Anderson spin-fermion model. In the anti-adiabatic regime (where the phonon energy exceeds the electronic scales) the standard treatment based on Lang-Firsov canonical transformation allows one to map the original model on to an attractive Hubbard model in the spinful case. The separate analysis of the two regimes supports the numerical evidence that the presence of well-defined polaronic lattice displacements is not necessarily associated to a metal-insulator transition, which is instead due to pairing between the carriers. The finite-dimensionality effects neglected in DMFT may lead to a finite conductivity in the bipolaronic state which is, however, not always associated with polaronic distortions. At the polaron crossover the Born-Oppenheimer approximation is shown to break down due to the entanglement of the electron-phonon state.
We re-examine the Peierls insulator to Mott insulator transition scenario in the one-dimensional Holstein-Hubbard model where, at half-filling, electron-phonon and electron-electron interactions compete for establishing charge- and spin-density-wave states, respectively. By means of large-scale density-matrix renormalization group calculations we determine the spin, single-particle and two-particle excitation gaps and prove —in the course of a careful finite-size scaling analysis— recent claims for an intervening metallic phase in the weak-coupling regime. We show that for large phonon frequencies the metallic region is even more extended than previously expected, and subdivided into ordinary Luttinger-liquid and bipolaronic-liquid phases.
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.