The repulsive two-dimensional Hubbard model is analysed within a functional renormalisation group (RG) approach. The flow equation for the effective interactions is evaluated on 1-loop level. The effective interactions diverge at a finite energy scale which is exponentially small for small bare interactions. To analyse the nature of the instabilities signalled by the diverging interactions, we compute the flow of the superconducting susceptibilities for various pairing symmetries and also charge and spin density susceptibilities. Depending on the choice of the model parameters (hopping amplitudes, interaction strength and band-filling) we find spin density wave instabilities or d-wave superconductivity as leading instability close to half-filling. This calculation establishes the existence of d-wave superconductivity in the Hubbard model beyond doubt.
We report on a numerical study of the two-dimensional Hubbard model and describe two new algorithms for the simulation of many-electron systems. These algorithms allow one to carry out simulations within the grand canonical ensemble at significantly lower temperatures than had previously been obtained and to calculate ground-state properties with fixed numbers of electrons. We present results for the two-dimensional Hubbard model with half- and quarter-filled bands. Our results support the existence of long-range antiferromagnetic order in the ground state at half-filling and its absence at quarter-filling. Results for the magnetic susceptibility and the momentum occupation along with an upper bound to the spin-wave spectrum are given. We find evidence for an attractive effective d-wave pairing interaction near half-filling but have not found evidence for a phase transition to a superconducting state.
In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the
d-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature
Tc, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the
d-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.