We implement a quantum Monte Carlo calculation for a repulsive Hubbard model with nearest- and next-nearest-neighbor hopping interactions on clusters up to 12\ifmmode\times\else\texttimes\fi{}12. A parameter region where the Fermi level lies close to the van Hove singularity at the saddle points in the bulk band structure is investigated. A pairing tendency in the ${\mathit{d}}_{{\mathit{x}}^{2}\mathrm{\ensuremath{-}}{\mathit{y}}^{2}}$ symmetry channel, but no other channel, is found. Estimates of the effective pairing interaction show that it is close to the value required for a 40 K superconductor. Finite-size scaling compares with the attractive Hubbard model. \textcopyright{}1996 The American Physical Society.
We implement a quantum Monte Carlo calculation for a repulsive Hubbard model with nearest- and next-nearest-neighbor hopping interactions on clusters up to 12\ifmmode\times\else\texttimes\fi{}12. A parameter region where the Fermi level lies close to the van Hove singularity at the saddle points in the bulk band structure is investigated. A pairing tendency in the dx2\ensuremath−y2 symmetry channel, but no other channel, is found. Estimates of the effective pairing interaction show that it is close to the value required for a 40 K superconductor. Finite-size scaling compares with the attractive Hubbard model. \textcopyright{}1996 The American Physical Society.
Quantum Monte Carlo is used to calculate various pairing correlations of the 2D Hubbard model possessing band features experimentally observed in the cuprates. In the hole-doped case, where the Fermi level lies close to the van Hove singularities around (0,π), the dx2−y2 pairing correlation is selectively enhanced, while in the electron-doped case, where the singularities are far below the Fermi level and the Fermi surface runs through (±π/2,±π/2), both dx2−y2 and dxy correlations are enhanced with the latter having a 2×2 structure. The two pairing symmetries can mix to result in a nodeless gap.
A quantum critical point is found in the phase diagram of the two-dimensional Hubbard model [Vidhyadhiraja et al., Phys. Rev. Lett. 102, 206407 (2009)]. It is due to the vanishing of the critical temperature associated with a phase-separation transition, and it separates the non-Fermi-liquid region from the Fermi liquid. Near the quantum critical point, the pairing is enhanced since the real part of the bare d-wave pairing susceptibility exhibits an algebraic divergence with decreasing temperature, replacing the logarithmic divergence found in a Fermi liquid [Yang et al., Phys. Rev. Lett. 106, 047004 (2011)]. In this paper, we explore the single-particle and transport properties near the quantum critical point using high-quality estimates of the self-energy obtained by direct analytic continuation of the self-energy from the continuous-time quantum Monte Carlo method. We focus mainly on a van Hove singularity coming from the relatively flat dispersion that crosses the Fermi level near the quantum critical filling. The flat part of the dispersion orthogonal to the antinodal direction remains pinned near the Fermi level for a range of doping that increases when we include a negative next-near-neighbor hopping ${t}^{\ensuremath{'}}$ in the model. For comparison, we calculate the bare d-wave pairing susceptibility for noninteracting models with the usual two-dimensional tight-binding dispersion and a hypothetical quartic dispersion. We find that neither model yields a van Hove singularity that completely describes the critical algebraic behavior of the bare d-wave pairing susceptibility found in the numerical data. The resistivity, thermal conductivity, thermopower, and the Wiedemann-Franz law are examined in the Fermi liquid, marginal Fermi liquid, and pseudogap doping regions. A negative next-near-neighbor hopping ${t}^{\ensuremath{'}}$ increases the doping region with marginal Fermi liquid character. Both T and negative ${t}^{\ensuremath{'}}$ are relevant variables for the quantum critical point, and both the transport and the displacement of the van Hove singularity with filling suggest that they are qualitatively similar in their effect.
Using large-scale dynamical cluster quantum Monte Carlo simulations, we study the Lifshitz transition of the two-dimensional Hubbard model with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. At ${t}^{\ensuremath{'}}\ensuremath{\le}0$, we identify a line of Lifshitz transition points associated with a change in the Fermi surface topology at zero temperature. In the overdoped region, the Fermi surface is complete and electron-like; across the Lifshitz transition, the Fermi surface becomes hole-like and develops a pseudogap. At (or very close to) the Lifshitz transition points, a van Hove singularity in the density of states crosses the Fermi level. The van Hove singularity occurs at finite doping due to correlation effects and becomes more singular when ${t}^{\ensuremath{'}}$ becomes more negative. The resulting temperature dependence on the bare d-wave pairing susceptibility close to the Lifshitz points is significantly different from that found in the traditional van Hove scenarios. Such unambiguous numerical observation of the Lifshitz transition at ${t}^{\ensuremath{'}}\ensuremath{\le}0$ extends our understanding of the quantum critical region in the phase diagram and shines lights on future investigations of the nature of the quantum critical point in the two-dimensional Hubbard model.
The search for superconductivity with unconventional pairing symmetry has been a central focus in the study of strongly correlated electron systems. In this work, we report a numerically exact study of the spin-triplet pairing in a two-dimensional Hubbard model with repulsive interactions, employing the determinant quantum Monte Carlo method. The model includes next-nearest-neighbor and third-nearest-neighbor hopping terms and maintains spin balance. By tuning the fermion filling close to a type-II Van Hove singularity (VHS) in the model, we numerically investigate the ordering tendencies of several possible pairing channels with different symmetries. Our numerical results provide clear evidence of the spin-triplet p-wave pairing instability approaching low temperatures, as revealed by the vertex contribution to the pairing susceptibility. This signature becomes increasingly pronounced as the interaction strength increases in the weak to intermediate regime. We further find that, near the type-II VHS, the dominant spin-spin correlations in the system are ferromagnetic, suggesting its close relation to the spin-triplet pairing instability. Our findings offer a reliable approach to realize spin-triplet p-wave superfluidity in the repulsive Hubbard model from an unbiased numerical perspective.
A combination of analytical approaches and quantum Monte Carlo simulations is used to study both magnetic and pairing correlations for a version of the Hubbard model that includes second-neighbor hopping ${\mathit{t}}^{\ensuremath{'}}$=-0.35t as a model for high-temperature superconductors. Magnetic properties are analyzed using the two-particle self-consistent approach. The maximum in magnetic susceptibility as a function of doping appears both at finite ${\mathit{t}}^{\ensuremath{'}}$ and at ${\mathit{t}}^{\ensuremath{'}}$=0 but for two totally different physical reasons. When ${\mathit{t}}^{\ensuremath{'}}$=0, it is induced by antiferromagnetic correlations while at ${\mathit{t}}^{\ensuremath{'}}$=-0.35t it is a band structure effect amplified by interactions. Finally, pairing fluctuations are compared with T-matrix results to disentangle the effects of van Hove singularity and of nesting on superconducting correlations. The addition of antiferromagnetic fluctuations increases slightly the d-wave superconducting correlations despite the presence of a van Hove singularity which tends to decrease them in the repulsive model. Some aspects of the phase diagram and some subtleties of finite-size scaling in Monte Carlo simulations, such as inverted finite-size dependence, are also discussed.
Quasiparticle bands of the two-dimensional Hubbard model are calculated using the Roth two-pole approximation to the one-particle Green's function. Excellent agreement is obtained with recent Monte Carlo calculations, including an anomalous volume of the Fermi surface near half-filling, which can possibly be explained in terms of a breakdown of Fermi liquid theory. The calculated bands are very flat around the (\ensuremath{\pi},0) points of the Brillouin zone in agreement with photoemission measurements of cuprate superconductors. With doping there is a shift in spectral weight from the upper band to the lower band. The Roth method is extended to deal with superconductivity within a four-pole approximation allowing electron-hole mixing. It is shown that triplet p-wave pairing never occurs. A self-consistent solution with singlet ${\mathit{d}}_{{\mathit{x}}^{2}\mathrm{\ensuremath{-}}{\mathit{y}}^{2}}$-wave pairing is found and optimal doping occurs when the van Hove singularity, corresponding to the flat band part, lies at the Fermi level. Nearest-neighbor antiferromagnetic correlations play an important role in flattening the bands near the Fermi level and in favoring superconductivity. However, the mechanism for superconductivity is a local one, in contrast to spin-fluctuation exchange models. For reasonable values of the hopping parameter the transition temperature Tc is in the range 10--100 K. The optimum doping ${\mathrm{\ensuremath{\delta}}}_{\mathit{c}}$ lies between 0.14 and 0.25, depending on the ratio U/t. The gap equation has a BCS-like form and 2${\mathrm{\ensuremath{\Delta}}}_{\mathrm{max}}$/kTc\ensuremath{\simeq}4.