We investigate the doped two-dimensional Hubbard model at finite temperature using controlled diagrammatic Monte Carlo calculations allowing for the computation of spectral properties in the infinite-size limit and, crucially, with arbitrary momentum resolution. We show that three distinct regimes are found as a function of doping and interaction strength, corresponding to a weakly correlated metal with properties close to those of the non-interacting system, a correlated metal with strong interaction effects including a reshaping of the Fermi surface, and a pseudogap regime at low doping in which quasiparticle excitations are selectively destroyed near the antinodal regions of momentum space. We study the physical mechanism leading to the pseudogap and show that it forms both at weak coupling when the magnetic correlation length is large and at strong coupling when it is shorter. In both cases, we show that spin-fluctuation theory can be modified in order to account for the behavior of the non-local component of the self-energy. We discuss the fate of the pseudogap as temperature goes to zero and show that, remarkably, this regime extrapolates precisely to the ordered stripe phase found by ground-state methods. This handshake between finite temperature and ground-state results significantly advances the elaboration of a comprehensive picture of the physics of the doped Hubbard model.
The relationship between the pseudogap and underlying ground-state phases has not yet been rigorously established. We investigated the doped two-dimensional Hubbard model at finite temperature using controlled diagrammatic Monte Carlo calculations, allowing for the computation of spectral properties in the infinite-size limit and with arbitrary momentum resolution. We found three distinct regimes as a function of doping and interaction strength: a weakly correlated metal, a correlated metal with strong interaction effects, and a pseudogap regime at low doping. We show that the pseudogap forms both at weak coupling, when the magnetic correlation length is large, and at strong coupling, when it is shorter. As the temperature goes to zero, the pseudogap regime extrapolates precisely to the ordered stripe phase found by ground-state methods.
Exact diagrammatic quantum Monte Carlo (DiagMC) results for the nearest-neighbor Hubbard model motivate a closer study of the pseudogap. Using the improved two-particle self-consistent approach (TPSC+), we analyze Energy (EDC) and Momentum Distribution Curves (MDC) simultaneously. We show that Fermi-liquid terminology breaks down in the pseudogap regime, requiring a distinction between true and false Fermi surfaces and zero-energy quasiparticle (ZEQ) lines. A false Fermi surface has a momentum-space spectral maximum but a frequency-space depression at zero energy, while a false ZEQ line violates the standard quasiparticle condition
∂Σ′(k,ω)/∂ω∣ω=0<0.The pseudogap is driven by critical thermal spin fluctuations, which occur in two dimensions because of the Mermin-Wagner theorem. Commensurate fluctuations first open an antinodal pseudogap, leaving true Fermi arcs. With decreasing temperature, the Fermi surface evolves into hole- and electron-like false Fermi surfaces. Incommensurate fluctuations generate hot spots near the diagonal, where the pseudogap persists to the quantum critical point (QCP), whereas at
kAN it disappears before the QCP doping. There, the spectrum has two precursor antiferromagnetic (AFM) bands, both in the unoccupied (
ω>0) region. We benchmark TPSC+ against DiagMC results for the Matsubara spectral proxy
−Im[G(k,iπT)]/π. TPSC+ underestimates pseudogap suppression in the strong-interaction regime, reproducing DiagMC behavior at lower temperatures or doping. This proxy is equivalent to the spectral function with thermal broadening
η=πT, which obscures features when
πT is not the smallest energy scale. Using
A(k,0), we find that hole-like Fermi surfaces emerge at any interaction strength at low temperature, even at weak coupling.
Long-wavelength spin fluctuations prohibit antiferromagnetic long-range order at finite temperature in two dimensions. Nevertheless, the correlation length starts to grow rapidly at a crossover temperature, leading to critical slowing down and to a renormalized-classical regime over a wide range of temperature, between a fraction of the mean-field transition temperature and the zero-temperature ordered state. This leads to a single-particle pseudogap of the kind observed in electron-doped cuprates. Very few theoretical methods can claim an accurate description of this regime. The challenge is that in this regime Fermi-liquid quasiparticles are already destroyed and new quasiparticles of the ordered state are not fully formed yet. Here, we study this problem for the two-dimensional Hubbard model by first generalizing the two-particle self-consistent approach. Using a special algorithm, spin fluctuations are treated in the thermodynamic limit even for large correlation lengths. The effects of Kanamori-Br\"uckner screening, of classical and of quantum fluctuations are taken into account. Results are presented at half filling for the one-band Hubbard model with nearest-neighbor hopping. They agree well with available benchmark diagrammatic quantum Monte Carlo at high temperature where the pseudogap opens up. In addition to temperature-dependent spectral properties, we find quantum corrections to the zero-temperature renormalized mean-field antiferromagnetic gap. Finally, analytic continuation of the Matsubara results for spectral functions show that the pseudogap opens up at significantly higher temperature than was previously identified based on the Matsubara data only.
We compute high-resolution angle-resolved photoemission spectroscopy of the Hubbard model using the unbiased determinant quantum Monte Carlo algorithm, revealing an asymmetry between electron and hole doping. Electron doping exhibits more coherent quasiparticles and stronger antiferromagnetic correlations compared to hole doping. At low doping, a nodal-antinodal dichotomy on the Fermi surface is observed, similar to cuprate experiments. The dichotomy reflects the momentum dependence of the Mott gap, as manifested in both the spectral function and the self-energy. For hole doping, we observe a transition towards the pseudogap, without signature of pocket formation. The simulated nuclear magnetic resonance pseudogap temperatures do not necessarily agree with the temperature determined by spectroscopy. These findings collectively suggest the pseudogap is a smooth crossover driven by strong correlations.
We study pseudogap behavior in a metal near a spin density wave (SDW) instability due to thermal magnetic fluctuations. We consider the $t\ensuremath{-}{t}^{\ensuremath{'}}$ Hubbard model on a square lattice at a finite doping, at intermediate coupling strength, and analyze the thermal evolution of the electron spectral function between a SDW ordered state at low temperatures and a normal Fermi liquid at high temperatures. We argue that for proper description of the pseudogap one needs to sum up infinite series of diagrams for both the fermionic self-energy and the SDW order parameter in the SDW state or the magnetic correlation length in the paramagnetic state. We use the eikonal approach to sum up an infinite series of diagrammatic contributions from thermal fluctuations. Earlier studies found that in the SDW state the spectral function ${A}_{\mathbit{k}}(\ensuremath{\omega})$ of a hot fermion at a finite
T is exponentially small below the energy scale $\mathrm{\ensuremath{\Delta}}(T)$, which scales with SDW order and vanishes at the ordering temperature
TN, and has a hump at a larger frequency ${\mathrm{\ensuremath{\Delta}}}_{\mathrm{PG}}$, comparable to the zero-temperature SDW gap $\mathrm{\ensuremath{\Delta}}(T=0)$. We argue that the hump, which we associate with the pseudogap, survives in some
T range above
TN. We show that this range is split by regions of strong and weak pseudogap behavior. In the first region, ${\mathrm{\ensuremath{\Delta}}}_{\mathrm{PG}}$ is weakly temperature dependent, despite the fact that it comes from thermal fluctuations. Such a behavior has been seen in numerical studies of the Hubbard model. We show that to obtain it one needs to go beyond the one-loop approximation and sum up the infinite series of diagrams. In the second regime, ${\mathrm{\ensuremath{\Delta}}}_{\mathrm{PG}}$ decreases with increasing
T and eventually vanishes. We further argue that a magnetic pseudogap at a finite
T emerges only if the ground state is magnetically ordered. We present the phase diagram and apply the results to high-
Tc cuprates.
We introduce a method for summing Feynman's perturbation series based on diagrammatic Monte Carlo that significantly improves its convergence properties. This allows us to investigate in a controllable manner the pseudogap regime of the Hubbard model and to study the nodal/antinodal dichotomy at low doping and intermediate coupling. Marked differences from the weak-coupling scenario are manifest, such as a higher degree of incoherence at the antinodes than at the ``hot spots''. Our results show that the pseudogap and reduction of quasiparticle coherence at the antinode is due to antiferromagnetic spin correlations centered around the commensurate $(\ensuremath{\pi},\ensuremath{\pi})$ wave vector. In contrast, the dominant source of scattering at the node is associated with incommensurate momentum transfer. Umklapp scattering is found to play a key role in the nodal/antinodal dichotomy.