Superconductivity from repulsion: Variational results for the 2D Hubbard model in the limit of weak interaction

Abstract

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<4t. Both the gap parameter and the order parameter have a dome-like shape as a function of n with a maximum for n about 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 U_c. 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.

20 pages, 13 figures

References in corpus (13)

• Diagrammatic routes to nonlocal correlations beyond dynamical mean field theory
• Superconductivity in the repulsive Hubbard model: an asymptotically exact weak-coupling solution
• Unconventional Superconductivity
• Variational Monte Carlo Method Combined with Quantum-Number Projection and Multi-Variable Optimization
• Role of backflow correlations for the non-magnetic phase of the t-t' Hubbard model
• Crossover between BCS Superconductor and Doped Mott Insulator of d-wave Pairing State in Two-Dimensional Hubbard Model
• Competition among various charge-inhomogeneous states and d-wave superconducting state in Hubbard models on square lattices
• Energetics of superconductivity in the two dimensional Hubbard model
• Superconductivity and antiferromagnetism in the two-dimensional Hubbard model: a variational study
• Hidden Mott transition and large-$U$ superconductivity in the two-dimensional Hubbard model
• Multiloop functional renormalization group for the two-dimensional Hubbard model: Loop convergence of the response functions
• Crossover from weakly to strongly correlated regions in the two-dimensional Hubbard model -- Off-diagonal wave function Monte Carlo studies of Hubbard model II --
• Antiferromagnetism, Superconductivity and Phase Diagram in the Two-Dimensional Hubbard Model -- Off-Diagonal Wave Function Monte Carlo Studies of Hubbard Model III --

Cited by in corpus (2)

• Ground state properties of one-dimensional and two-dimensional Hubbard model from Gutzwiller conjugate gradient minimization theory
• Foundations of the Variational Discrete Action Theory