We present the results of a quantum Monte Carlo study of the extended
s and the ${d}_{{x}^{2}{\ensuremath{-}y}^{2}}$ pairing correlation functions for the two-dimensional Hubbard model, computed with the constrained-path Monte Carlo method. For small lattice sizes and weak interactions, we find that the ${d}_{{x}^{2}{\ensuremath{-}y}^{2}}$ pairing correlations are stronger than the extended
s pairing correlations and are positive when the pair separation exceeds several lattice constants. As the system size or the interaction strength increases, the magnitude of the long-range part of both correlation functions vanishes.
Superconducting properties of coupled Hubbard planes are studied by quantum Monte Carlo simulations. It is found that the interlayer coupling supports the pairing mechanism but it does not lead to off-diagonal long-range order in the purely electronic system. In addition, we have performed Monte Carlo studies based on a recently proposed variational wave function. Yielding identical results, the variational wave function overcomes the sign problem and allows us to analyze ground-state properties of strongly correlated many-body problems for system sizes and complexities not accessible otherwise.
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 this study we have computed the pair correlation functions in the two-dimensional Hubbard model using a quantum Monte Carlo method. We employ a new diagonalization algorithm in the quantum Monte Carlo method which is free from the negative sign problem. We show that the d-wave pairing correlation function is indeed enhanced slightly for the positive on-site Coulomb interaction U when doping away from the half-filling. When the system size becomes large, the pair correlation function P d increases for U > 0 compared to the non-interacting case, while P d is suppressed for U > 0 when the system size is small. The enhancement ratio P d [ U ]/ P d [ U = 0] will give a criterion on the existence of superconductivity. The ratio P d [ U ]/ P d [ U = 0] increases almost linearly ∝ L when the system size L × L is increased. This increase is a good indication of the existence of a superconducting phase in the two-dimensional Hubbard model. There is, however, no enhancement of pair correlation functions in the half-filled case, which indicates the absence of superconductivity without hole doping.
We investigate the ground-state properties of the two-dimensional Hubbard model, based on the off-diagonal wave function variational Monte Carlo method. We use an optimized wave function that is improved from an initial one-body wave function by multiplying by multiple correlation operators that are given by the form of exp(
−S)-type where
S is a suitable operator. We examine the inhomogeneous ground state at near 1/8 doping in the strongly correlated region where the on-site Coulomb interaction is larger than the bandwidth. The
d-wave superconductivity with a spatially oscillating gap function can coexist with the charge ordering in the ground state without magnetic ordering, and superconducting condensation energy increases by this coexistence. We also show the
d-wave pair correlation function as a function of lattice sites. The correlation function indicates that the long-range superconducting order indeed exists, and also that the strong electron correlation suppresses the pair correlation function.
We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the
d-wave and extended
s-wave pairing symmetries. The pairing structure factor shows a strong tendency to
d-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended
s-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the
d-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the
d-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.