We thoroughly analyze the divergences of the irreducible vertex functions occurring in the charge channel of the half-filled Hubbard model in close proximity to the Mott metal-insulator transition (MIT). In particular, by systematically performing dynamical mean-field theory (DMFT) calculations on the two-particle level, we determine the location and the number of the vertex divergences across the whole coexistence region adjacent to the first-order metal-to-insulator transition. We find that the lines in the parameter space, along which the vertex divergences occur, display a qualitatively different shape in the coexisting metallic and insulating phase, which is also associated to an abrupt jump of the number of divergences across the MIT. Physically, the systematically larger number of divergences on the insulating side of the transition reflects the sudden suppression of local charge fluctuation at the MIT. Further, a systematic analysis of the results demonstrates that the number of divergence lines increases as a function of the inverse temperature $\ensuremath{\beta}={({k}_{\mathrm{B}}T)}^{\ensuremath{-}1}$ by approaching the Mott transition in the zero-temperature limit. This makes it possible to identify the zero-temperature MIT as an accumulation point of an infinite number of vertex divergence lines, unveiling the highly nonperturbative nature of the underlying transition.
We thoroughly analyze the divergences of the irreducible vertex functions occurring in the charge channel of the half-filled Hubbard model in close proximity to the Mott metal-insulator transition (MIT). In particular, by systematically performing dynamical mean-field theory (DMFT) calculations on the two-particle level, we determine the location and the number of the vertex divergences across the whole coexistence region adjacent to the first-order metal-to-insulator transition. We find that the lines in the parameter space, along which the vertex divergences occur, display a qualitatively different shape in the coexisting metallic and insulating phase, which is also associated to an abrupt jump of the number of divergences across the MIT. Physically, the systematically larger number of divergences on the insulating side of the transition reflects the sudden suppression of local charge fluctuation at the MIT. Further, a systematic analysis of the results demonstrates that the number of divergence lines increases as a function of the inverse temperature \ensuremathβ=(kBT)\ensuremath−1 by approaching the Mott transition in the zero-temperature limit. This makes it possible to identify the zero-temperature MIT as an accumulation point of an infinite number of vertex divergence lines, unveiling the highly nonperturbative nature of the underlying transition.
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.
Here we analyze the divergences of the irreducible vertex function in dynamical mean field theory, which may indicate either a nonphysical breakdown of the perturbation theory or a response to some physical phenomenon. To investigate this question, we construct a quasiparticle vertex from the diverging irreducible vertex functions. This vertex describes the scattering between quasiparticles and quasiholes in a Fermi liquid. We show that the quasparticle vertex does not diverge in the charge channel, wherein the irreducible vertex does diverge; and we show that the quasiparticle vertex does diverge in the spin channel, wherein the irreducible vertex does not diverge. This divergence occurs at the Mott transition, wherein the Fermi-liquid theory breaks down. Both the half filled Hubbard and Anderson lattices are investigated. In general, our results support that the divergences of the irreducible vertex function do not indicate a nonphysical failure of the perturbation theory. Instead, the divergences are the mathematical consequence of inverting a matrix (the local charge susceptibility) which accumulates increasingly negative diagonal elements as the Hubbard interaction suppresses charge fluctuations. Indeed, we find that the first divergences of the irreducible vertex in both Hubbard and Anderson lattices occurs near the maximum magnitude of the (negative) vertex-connected part of the charge susceptibility.
In this work, we analyze in detail the occurrence of divergences in the irreducible vertex functions for one of the fundamental models of many-body physics: the Anderson impurity model (AIM). These divergences, a surprising hallmark of the breakdown of many-electron perturbation theory, have been recently observed in several contexts, including the dynamical mean-field solution of the Hubbard model. The numerical calculations for the AIM presented in this work, as well as their comparison with the corresponding results for the Hubbard model, allow us to clarify several open questions about the properties of vertex divergences in a particularly interesting context, the correlated metallic regime at low temperatures. Specifically, our analysis (i) rules out explicitly the transition to a Mott-insulating phase, but not the more general suppression of charge fluctuations (proposed in [O. Gunnarsson et al., Phys. Rev. B 93, 245102 (2016)]), as a necessary condition for the occurrence of vertex divergences, (ii) clarifies their relation with the underlying Kondo physics, and, eventually, (iii) individuates which divergences might also appear on the real-frequency axis in the limit of zero temperature, through the discovered scaling properties of the singular eigenvectors.
At the Mott transition, electron-electron interaction changes a metal, in which electrons are itinerant, to an insulator, in which electrons are localized. This phenomenon is central to quantum materials. Here we contribute to its understanding by studying the two-dimensional Hubbard model at finite temperature with plaquette cellular dynamical mean-field theory. We provide an exhaustive thermodynamic description of the correlation-driven Mott transition of the half-filled model by calculating pressure, charge compressibility, entropy, kinetic energy, potential energy, and free energy across the first-order Mott transition and its high-temperature crossover (Widom line). The entropy is extracted from the Gibbs-Duhem relation and shows complex behavior near the transition, marked by discontinuous jumps at the first-order boundary, singular behavior at the Mott endpoint, and inflections marking sharp variations in the supercritical region. The free energy allows us to identify the thermodynamic phase boundary, to discuss phases stability and metastability, and to touch upon nucleation and spinodal decomposition mechanisms for the transition. We complement this thermodynamic description of the Mott transition by an information-theoretic description. We achieve this by calculating the local entropy, which is a measure of entanglement, and the single-site total mutual information, which quantifies quantum and classical correlations. These information-theoretic measures exhibit characteristic behaviors that allow us to identify the first-order coexistence regions, the Mott critical endpoint and the crossovers along the Widom line in the supercritical region.
For doped two-dimensional Mott insulators in their normal state, the challenge is to understand the evolution from a conventional metal at high doping to a strongly correlated metal near the Mott insulator at zero doping. To this end, we solve the cellular dynamical mean-field equations for the two-dimensional Hubbard model using a plaquette as the reference quantum impurity model and continuous-time quantum Monte Carlo method as impurity solver. The normal-state phase diagram as a function of interaction strength U, temperature T, and filling n shows that, upon increasing n toward the Mott insulator, there is a surface of first-order transition between two metals at nonzero doping. That surface ends at a finite temperature critical line originating at the half-filled Mott critical point. Associated with this transition, there is a maximum in scattering rate as well as thermodynamic signatures. These findings suggest a new scenario for the normal-state phase diagram of the high temperature superconductors. The criticality surmised in these systems can originate not from a T = 0 quantum critical point, nor from the proximity of a long-range ordered phase, but from a low temperature transition between two types of metals at finite doping. The influence of Mott physics therefore extends well beyond half-filling.