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

Bethe-Salpeter Equation at the Critical End Point of the Mott Transition

Erik G. C. P. van Loon, Friedrich Krien, A. A. Katanin

Published Sep 23, 2020• 43 citations
arXivHTMLPDF

Abstract

Strong repulsive interactions between electrons can lead to a Mott metal-insulator transition. The dynamical mean-field theory (DMFT) explains the critical end point and the hysteresis region usually in terms of single-particle concepts, such as the spectral function and the quasiparticle weight. In this Letter, we reconsider the critical end point of the metal-insulator transition on the DMFT's two-particle level. We show that the relevant eigenvalue and eigenvector of the nonlocal Bethe-Salpeter kernel in the charge channel provide a unified picture of the hysteresis region and of the critical end point of the Mott transition. In particular, they simultaneously explain the thermodynamics of the hysteresis region and the iterative stability of the DMFT equations. This analysis paves the way for a deeper understanding of phase transitions in correlated materials.

DOI
10.1103/physrevlett.125.136402

Related papers

Nonperturbative landscape of the Mott-Hubbard transition: Multiple divergence lines around the critical endpoint

Thomas Schäfer, Sergio Ciuchi, Markus Wallerberger +5 more

We analyze the highly nonperturbative regime surrounding the Mott-Hubbard metal-to-insulator transition (MIT) by means of dynamical mean field theory (DMFT) calculations at the two-particle level. By extending the results of Sch\"afer et al. [Phys. Rev. Lett. 110, 246405 (2013)] we show the existence of infinitely many lines in the phase diagram of the Hubbard model where the local Bethe-Salpeter equations, and the related irreducible vertex functions, become singular in the charge as well as the particle-particle channel. By comparing our numerical data for the Hubbard model with analytical calculations for exactly solvable systems of increasing complexity [disordered binary mixture (BM), Falicov-Kimball (FK), and atomic limit (AL)], we have (i) identified two different kinds of divergence lines; (ii) classified them in terms of the frequency structure of the associated singular eigenvectors; and (iii) investigated their relation to the emergence of multiple branches in the Luttinger-Ward functional. In this way, we could distinguish the situations where the multiple divergences simply reflect the emergence of an underlying, single energy scale ${\ensuremath{\nu}}^{*}$ below which perturbation theory is no longer applicable, from those where the breakdown of perturbation theory affects, not trivially, different energy regimes. Finally, we discuss the implications of our results on the theoretical understanding of the nonperturbative physics around the MIT and for future developments of many-body algorithms applicable in this regime.

arXiv:1606.03393•Physical review. B./Physical review. B•Dec 5, 2016•96 citations

Divergent Precursors of the Mott-Hubbard Transition at the Two-Particle Level

Thomas Schäfer, G. Rohringer, O. Gunnarsson +3 more

Identifying the fingerprints of the Mott-Hubbard metal-insulator transition may be quite elusive in correlated metallic systems if the analysis is limited to the single particle level. However, our dynamical mean-field calculations demonstrate that the situation changes completely if the frequency dependence of the two-particle vertex functions is considered: The first nonperturbative precursors of the Mott physics are unambiguously identified well inside the metallic regime by the divergence of the local Bethe-Salpeter equation in the charge channel. In the low-temperature limit this occurs for interaction values where incoherent high-energy features emerge in the spectral function, while at high temperatures it is traceable up to the atomic limit.

arXiv:1303.0246•Physical Review Letters•Jun 13, 2013•132 citations

Fate of Quasiparticle at Mott Transition and Interplay with Lifshitz Transition Studied by Correlator Projection Method

Kota Hanasaki, Masatoshi Imada

Filling-control metal-insulator transition on the two-dimensional Hubbard model is investigated by using the correlator projection method, which takes into account momentum dependence of the free energy beyond the dynamical mean-field theory. The phase diagram of metals and Mott insulators is analyzed. Lifshitz transitions occur simultaneously with metal-insulator transitions at large Coulomb repulsion. On the other hand, they are separated each other for lower Coulomb repulsion, where the phase sandwiched by the Lifshitz and metal-insulator transitions appears to show violation of the Luttinger sum rule. Through the metal-insulator transition, quasiparticles retain nonzero renormalization factor and finite quasi-particle weight in the both sides of the transition. This supports that the metal-insulator transition is caused not by the vanishing renormalization factor but by the relative shift of the Fermi level into the Mott gap away from the quasiparticle band, in sharp contrast with the original dynamical mean-field theory. Charge compressibility diverges at the critical end point of the first-order Lifshitz transition at finite temperatures. The origin of the divergence is ascribed to singular momentum dependence of the quasiparticle dispersion.

arXiv:cond-mat/0604395•Journal of the Physical Society of Japan•Jul 27, 2006•8 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

Cluster Dynamical Mean Field Theory of the Mott Transition

H. Park, Kristjan Haule, Gabriel Kotliar

We address the nature of the Mott transition in the Hubbard model at half-filling using cluster dynamical mean field theory (DMFT). We compare cluster-DMFT results with those of single-site DMFT. We show that inclusion of the short-range correlations on top of the on-site correlations does not change the order of the transition between the paramagnetic metal and the paramagnetic Mott insulator, which remains first order. However, the short range correlations reduce substantially the critical U and modify the shape of the transition lines. Moreover, they lead to very different physical properties of the metallic and insulating phases near the transition point. Approaching the transition from the metallic side, we find an anomalous metallic state with very low coherence scale. The insulating state is characterized by the narrow Mott gap with pronounced peaks at the gap edge.

arXiv:0803.1324•Physical Review Letters•Oct 30, 2008•306 citations

Dynamical mean field theory phase-space extension and critical properties of the finite temperature Mott transition

Hugo U. R. Strand, Andro Sabashvili, Mats Granath +2 more

We consider the finite temperature metal-insulator transition in the half-filled paramagnetic Hubbard model on the infinite dimensional Bethe lattice. A method for calculating the dynamical mean field theory fixpoint surface in the phase diagram is presented and shown to be free from the convergence problems of standard forward recursion. The fixpoint equation is then analyzed using dynamical systems methods. On the fixpoint surface the eigenspectra of its Jacobian is used to characterize the hysteresis boundaries of the first-order transition line and its second-order critical end point. The critical point is shown to be a cusp catastrophe in the parameter space, opening a pitchfork bifurcation along the first-order transition line, while the hysteresis boundaries are shown to be saddle-node bifurcations of two merging fixpoints. Using Landau theory the properties of the critical end point are determined and related to the critical eigenmode of the Jacobian. Our findings provide insights into basic properties of this intensively studied transition.

arXiv:1012.3829•Physical Review B•May 31, 2011•25 citations
Strong repulsive interactions between electrons can lead to a Mott metal-insulator transition. The dynamical mean-field theory (DMFT) explains the critical end point and the hysteresis region usually in terms of single-particle concepts, such as the spectral function and the quasiparticle weight. In this Letter, we reconsider the critical end point of the metal-insulator transition on the DMFT's two-particle level. We show that the relevant eigenvalue and eigenvector of the nonlocal Bethe-Salpeter kernel in the charge channel provide a unified picture of the hysteresis region and of the critical end point of the Mott transition. In particular, they simultaneously explain the thermodynamics of the hysteresis region and the iterative stability of the DMFT equations. This analysis paves the way for a deeper understanding of phase transitions in correlated materials.
We analyze the highly nonperturbative regime surrounding the Mott-Hubbard metal-to-insulator transition (MIT) by means of dynamical mean field theory (DMFT) calculations at the two-particle level. By extending the results of Sch\"afer et al. [Phys. Rev. Lett. 110, 246405 (2013)] we show the existence of infinitely many lines in the phase diagram of the Hubbard model where the local Bethe-Salpeter equations, and the related irreducible vertex functions, become singular in the charge as well as the particle-particle channel. By comparing our numerical data for the Hubbard model with analytical calculations for exactly solvable systems of increasing complexity [disordered binary mixture (BM), Falicov-Kimball (FK), and atomic limit (AL)], we have (i) identified two different kinds of divergence lines; (ii) classified them in terms of the frequency structure of the associated singular eigenvectors; and (iii) investigated their relation to the emergence of multiple branches in the Luttinger-Ward functional. In this way, we could distinguish the situations where the multiple divergences simply reflect the emergence of an underlying, single energy scale ${\ensuremath{\nu}}^{*}$ below which perturbation theory is no longer applicable, from those where the breakdown of perturbation theory affects, not trivially, different energy regimes. Finally, we discuss the implications of our results on the theoretical understanding of the nonperturbative physics around the MIT and for future developments of many-body algorithms applicable in this regime.
Identifying the fingerprints of the Mott-Hubbard metal-insulator transition may be quite elusive in correlated metallic systems if the analysis is limited to the single particle level. However, our dynamical mean-field calculations demonstrate that the situation changes completely if the frequency dependence of the two-particle vertex functions is considered: The first nonperturbative precursors of the Mott physics are unambiguously identified well inside the metallic regime by the divergence of the local Bethe-Salpeter equation in the charge channel. In the low-temperature limit this occurs for interaction values where incoherent high-energy features emerge in the spectral function, while at high temperatures it is traceable up to the atomic limit.
Filling-control metal-insulator transition on the two-dimensional Hubbard model is investigated by using the correlator projection method, which takes into account momentum dependence of the free energy beyond the dynamical mean-field theory. The phase diagram of metals and Mott insulators is analyzed. Lifshitz transitions occur simultaneously with metal-insulator transitions at large Coulomb repulsion. On the other hand, they are separated each other for lower Coulomb repulsion, where the phase sandwiched by the Lifshitz and metal-insulator transitions appears to show violation of the Luttinger sum rule. Through the metal-insulator transition, quasiparticles retain nonzero renormalization factor and finite quasi-particle weight in the both sides of the transition. This supports that the metal-insulator transition is caused not by the vanishing renormalization factor but by the relative shift of the Fermi level into the Mott gap away from the quasiparticle band, in sharp contrast with the original dynamical mean-field theory. Charge compressibility diverges at the critical end point of the first-order Lifshitz transition at finite temperatures. The origin of the divergence is ascribed to singular momentum dependence of the quasiparticle dispersion.
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.
We address the nature of the Mott transition in the Hubbard model at half-filling using cluster dynamical mean field theory (DMFT). We compare cluster-DMFT results with those of single-site DMFT. We show that inclusion of the short-range correlations on top of the on-site correlations does not change the order of the transition between the paramagnetic metal and the paramagnetic Mott insulator, which remains first order. However, the short range correlations reduce substantially the critical U and modify the shape of the transition lines. Moreover, they lead to very different physical properties of the metallic and insulating phases near the transition point. Approaching the transition from the metallic side, we find an anomalous metallic state with very low coherence scale. The insulating state is characterized by the narrow Mott gap with pronounced peaks at the gap edge.
We consider the finite temperature metal-insulator transition in the half-filled paramagnetic Hubbard model on the infinite dimensional Bethe lattice. A method for calculating the dynamical mean field theory fixpoint surface in the phase diagram is presented and shown to be free from the convergence problems of standard forward recursion. The fixpoint equation is then analyzed using dynamical systems methods. On the fixpoint surface the eigenspectra of its Jacobian is used to characterize the hysteresis boundaries of the first-order transition line and its second-order critical end point. The critical point is shown to be a cusp catastrophe in the parameter space, opening a pitchfork bifurcation along the first-order transition line, while the hysteresis boundaries are shown to be saddle-node bifurcations of two merging fixpoints. Using Landau theory the properties of the critical end point are determined and related to the critical eigenmode of the Jacobian. Our findings provide insights into basic properties of this intensively studied transition.