We obtain a complete and numerically exact in the weak-coupling limit $(U\ensuremath{\rightarrow}0)$ ground-state phase diagram of the repulsive fermionic Hubbard model on the square lattice for filling factors
0<n<2 and next-nearest-neighbor hopping amplitudes $0\ensuremath{\le}{t}^{\ensuremath{'}}\ensuremath{\le}0.5$. Phases are distinguished by the symmetry and the number of nodes of the superfluid order parameter. The phase diagram is richer than may be expected and typically features states with a high---higher than that of the fundamental mode of the corresponding irreducible representation---number of nodes. The effective coupling strength in the Cooper channel $\ensuremath{\lambda}$, which determines the critical temperature
Tc of the superfluid transition, is calculated in the whole parameter space and regions with high values of $\ensuremath{\lambda}$ are identified. It is shown that besides the expected increase of $\ensuremath{\lambda}$ near the Van Hove singularity line, joining the ferromagnetic and antiferromagnetic points, another region with high values of $\ensuremath{\lambda}$ can be found at quarter filling and ${t}^{\ensuremath{'}}=0.5$ due to the presence of a line of nesting at ${t}^{\ensuremath{'}}\ensuremath{\ge}0.5$. The results can serve as benchmarks for controlled nonperturbative methods and guide the ongoing search for high-
Tc superconductivity in the Hubbard model.
The phase diagram of the Hubbard model is studied at weak coupling in two and three spatial dimensions. It is shown that the N\'eel temperature and the order parameter in
d=3 are smaller than the Hartree-Fock predictions by a factor of
q=0.2599. For
d=2 we show that the self-consistent (sc) perturbation series bears no relevance to the behavior of the exact solution of the Hubbard model in the symmetry-broken phase. We also investigate an anisotropic model and show that the coupling between planes is essential for the validity of mean-field-type order parameters.
Motivated by the first experimental realization of the antiferromagnetic phase transition in the three-dimensional Fermi-Hubbard model, we present a theoretical study of the model's 3D superconducting phase. By formulating a 3D extension of the Feshbach mechanism, we provide a unified microscopic picture at strong coupling in which pairing is driven by near-resonant bound states of dopants. These long-lived bound states acquire a qualitatively different internal structure in three dimensions compared to their two-dimensional counterparts, giving rise to a distinct superconducting state, namely a time-reversal symmetry breaking
dx2−y2+idz2 pairing state. We also provide an estimate of the corresponding critical temperature, and provide evidence that the resulting superconducting phase hosts gapless Weyl points. Our results represent a new milestone for the field of quantum simulation, challenging experiments and large-scale numerics alike to test our predictions.
Extensive numerical studies have demonstrated that the two-dimensional single-band Hubbard model contains much of the key physics in cuprate high-temperature superconductors. However, there is no definitive proof that the Hubbard model truly possesses a superconducting ground state or, if it does, of how it depends on model parameters. To answer these longstanding questions, we study an extension of the Hubbard model including an infinite-range
d-wave pair field term, which precipitates a superconducting state in the
d-wave channel. Using exact diagonalization on 16-site square clusters, we study the evolution of the ground state as a function of the strength of the pairing term. This is achieved by monitoring the fidelity metric of the ground state, as well as determining the ratio between the two largest eigenvalues of the
d-wave pair/spin/charge-density matrices. The calculations show a
d-wave superconducting ground state in doped clusters bracketed by a strong antiferromagnetic state at half filling controlled by the Coulomb repulsion
U and a weak short-range checkerboard charge ordered state at larger hole doping controlled by the next-nearest-neighbor hopping ${t}^{\ensuremath{'}}$. We also demonstrate that negative ${t}^{\ensuremath{'}}$ plays an important role in facilitating
d-wave superconductivity.
The two-dimensional Hubbard model is studied for small values of the interaction strength (
U of the order of the hopping amplitude
t), using a variational ansatz well suited for this regime. The wave function, a refined Gutzwiller ansatz, has a BCS mean-field state with
d-wave symmetry as its reference state. Superconducting order is found for densities
n<1 (but not for
n=1). This resolves a discrepancy between results obtained with the functional renormalization group, which do predict superconducting order for small values of
U, and numerical simulations, which did not find superconductivity for $U\ensuremath{\lesssim}4t$. Both the gap parameter and the order parameter have a domelike shape as a function of
n with a maximum for $n\ensuremath{\approx}0.8$. Expectation values for the energy, the particle number, and the superconducting order parameter are calculated using a linked-cluster expansion up to second order in
U. In this way large systems (millions of sites) can be readily treated and well converged results are obtained. A big size is indeed required to see that the gap becomes very small at half filling and probably tends to zero in the thermodynamic limit, whereas away from half filling a finite asymptotic limit is reached. For a lattice of a given size the order parameter becomes finite only above a minimal coupling strength
Uc. This threshold value decreases steadily with increasing system size, which indicates that superconductivity exists for arbitrarily small
U for an infinite system. For moderately large systems the size dependence is quite irregular, due to variations in level spacings at the Fermi energy. The fluctuations die out if the gap parameter spans several level spacings.
In superconducting phases, repulsive interactions between the electrons usually lead to higher-angular momentum Cooper pairing. Here, the authors employ the weak-coupling renormalization group approach to study unconventional superconducting phases emerging in the extended, repulsive Hubbard model on paradigmatic two-dimensional lattices. Among their findings are extended regions of chiral topological superconductivity on the square, triangular, and honeycomb lattices. They provide detailed superconducting phase diagrams as well as the coupling strengths that quantify the corresponding critical temperatures depending on the band-structure parameters, band filling, and interaction parameters. The importance of the resolution and accuracy of the numerical integration for the method is discussed and a way to higher numerical efficiency provided.