No Arabic abstract
We present numerically exact results from sign-problem free quantum Monte Carlo simulations for a spin-fermion model near an $O(3)$ symmetric antiferromagnetic (AFM) quantum critical point. We find a hierarchy of energy scales that emerges near the quantum critical point. At high energy scales, there is a broad regime characterized by Landau-damped order parameter dynamics with dynamical critical exponent $z=2$, while the fermionic excitations remain coherent. The quantum critical magnetic fluctuations are well described by Hertz-Millis theory, except for a $T^{-2}$ divergence of the static AFM susceptibility. This regime persists down to a lower energy scale, where the fermions become overdamped and concomitantly, a transition into a $d-$wave superconducting state occurs. These findings resemble earlier results for a spin-fermion model with easy-plane AFM fluctuations of an $O(2)$ SDW order parameter, despite noticeable differences in the perturbative structure of the two theories. In the $O(3)$ case, perturbative corrections to the spin-fermion vertex are expected to dominate at an additional energy scale, below which the $z=2$ behavior breaks down, leading to a novel $z=1$ fixed point with emergent local nesting at the hot spots [Schlief et al., PRX 7, 021010 (2017)]. Motivated by this prediction, we also consider a variant of the model where the hot spots are nearly locally nested. Within the available temperature range in our study ($Tge E_F/200$), we find substantial deviations from the $z=2$ Hertz-Millis behavior, but no evidence for the predicted $z=1$ criticality.
Metallic quantum critical phenomena are believed to play a key role in many strongly correlated materials, including high temperature superconductors. Theoretically, the problem of quantum criticality in the presence of a Fermi surface has proven to be highly challenging. However, it has recently been realized that many models used to describe such systems are amenable to numerically exact solution by quantum Monte Carlo (QMC) techniques, without suffering from the fermion sign problem. In this article, we review the status of the understanding of metallic quantum criticality, and the recent progress made by QMC simulations. We focus on the cases of spin density wave and Ising nematic criticality. We describe the results obtained so far, and their implications for superconductivity, non-Fermi liquid behavior, and transport in the vicinity of metallic quantum critical points. Some of the outstanding puzzles and future directions are highlighted.
We explore the Matsubara quasiparticle fraction and the pseudogap of the two-dimensional Hubbard model with the dynamical cluster quantum Monte Carlo method. The character of the quasiparticle fraction changes from non-Fermi liquid, to marginal Fermi liquid to Fermi liquid as a function of doping, indicating the presence of a quantum critical point separating non-Fermi liquid from Fermi liquid character. Marginal Fermi liquid character is found at low temperatures at a very narrow range of doping where the single-particle density of states is also symmetric. At higher doping the character of the quasiparticle fraction is seen to cross over from Fermi Liquid to Marginal Fermi liquid as the temperature increases.
The origin of the pseudogap behavior, found in many high-$T_c$ superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a gap-filling rather than gap-closing behavior, consistent with experimental observations. Our results provide the first unambiguous lattice model realization of a pseudogap state in a strongly correlated system, driven by superconducting fluctuations.
We study a two-dimensional model of an isolated narrow topological band at partial filling with local attractive interactions. Numerically exact quantum Monte Carlo calculations show that the ground state is a superconductor with a critical temperature that scales nearly linearly with the interaction strength. We also find a broad pseudogap regime at temperatures above the superconducting phase that exhibits strong pairing fluctuations and a tendency towards electronic phase separation; introducing a small nearest neighbor attraction suppresses superconductivity entirely and results in phase separation. We discuss the possible relevance of superconductivity in this unusual regime to the physics of flat band moir{e} materials, and as a route to designing higher temperature superconductors.
Fixed-node Greens function Monte Carlo calculations have been performed for very large 16x6 2D Hubbard lattices, large interaction strengths U=10,20, and 40, and many (15-20) densities between empty and half filling. The nodes were fixed by a simple Slater-Gutzwiller trial wavefunction. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of U. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape and boundary condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions.