Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.
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.
The two-dimensional Hubbard model near half-filling is studied using a configuration-interaction method. Superconducting order-parameter and spin-spin correlation functions are calculated as functions of distance to probe for superconductivity and antiferromagnetism. Consistent with our previous results, these more extensive calculations point toward d-wave superconductivity in the ground state of the two-dimensional repulsive Hubbard model. The d-wave correlations tend to increase with interaction strength, size of the state space used in the calculation, and doping away from half-filling. The superconducting correlations on the Fermi surface are interpreted as d-wave pairing between spin density-wave quasiparticles.
We use a variational Monte Carlo technique to study the ground state of the two-dimensional Hubbard and t-J models on a square lattice. We use a trial wave function that allows a continuous description of the paramagnetic, antiferromagnetic, and superconducting phases, as well as the coexistence of these phases, with no a priori constraint on double occupancy. The phase diagram of both models is given for intermediate coupling (U=10t). We show that the apparent discrepancies between the two models, that appeared in previous variational Monte Carlo studies, are not present when this sufficiently good variational wave function is used. The two models are in good qualitative agreement. We find that except at half-filling a pure antiferromagnetic phase does not exist but is always in coexistence with a d-wave superconducting phase, and is followed by a pure d-wave superconducting phase at doping up to \ensuremath{\delta}\ensuremath{\sim}0.3. The staggered magnetization and the superconducting gap are measured, and a comparison with other analytical or numerical results is made.
Superconducting correlation of the Hubbard model on a square lattice is reinvestigated by quantum Monte Carlo method. Results from more than one order of magnitude larger Monte Carlo steps than before are presented to discuss the convergence of the data to the ground state. By examining the dependence on the choice of initial trial states as well as on the Trotter number in detail, the superconducting correlation in the ground state is reliably obtained up to 12×12 lattice away from the half-filling and up to 16×16 lattice at the half-filling, which substantially improve former results near the half-filling. Various channels of pairing in the ground state show no system size dependence and remain quite short ranged at any filling including the region close to the half-filling. Numerical data show further evidence that neither extended s - nor d -wave superconducting state is possible in the two-dimensional Hubbard model.
Magnetic properties of the two-dimensional Hubbard-model are inferred from results of Monte Carlo simulations. Lattice sizes up to 8\ifmmode\times\else\texttimes\fi{}8 and temperatures down to
T=20t (
t=hopping) were studied. The half-filled system is found to exhibit antiferromagnetic long-range order for all values of the Coulomb repulsion
U. The low-temperature magnetic properties are found to be well described by spin-wave theory with renormalized local moment and spin-wave velocity. Numerical evidence presented suggests that when doped the system loses the long-range order immediately away from half filling.