We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from ``standard'' staggered spin order to more ``exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (
xc), where long-range antiferromagnetic order disappears, are presented. It was found that
xc is very small in agreement with experiments for the high-
Tc superconductors. For example, in the t-J model, ${\mathit{x}}_{\mathit{c}}$0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, ${\mathrm{\ensuremath{\sigma}}}_{1}$(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$
SrxCuO4, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-
Tc superconductors.
Numerical results for the density of states to remove or add electrons to a finite Hubbard cluster are presented for both positive and negative on-site interactions. The behavior of the filling versus the chemical potential is discussed. At half filling there is a gap in the density of states. When the repulsive-U Hubbard model is doped spectral weight is transferred into the gap, and the chemical potential is found to move across the gap as one goes from electron to hole doping. We discuss the relationship of these results to experiment. For the negative-U Hubbard model, the chemical potential shifts continuously and lies in the middle of the superconducting gap.
We investigate the Hubbard model on a two-dimensional cluster of 18 sites using Lanczos algorithm for different strengths of Coulomb repulsion in the low doping region. Energies and correlation functions are given for half-filling and two holes. In the strong coupling regime we compare our results with the $t\ensuremath{-}J$ model. We confirm the tendency of holes to repel, although a significant value of the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductive order parameter is measured. These data represent a benchmark for future quantum simulations.
Numerical results obtained from Monte Carlo simulation in the ground state and at finite temperatures as well as the exact diagonalization and the transfer matrix method are reported. Efficiency of the Monte Carlo method in the ground state is examined. Spin, charge and superconducting correlations are investigated for the Hubbard and the t - J model in one and two dimensions. The momentum distribution in the one-dimensional Hubbard model shows fermi-liquid-like behavior at least in two orders of magnitude smaller energy scale than the band width. The short-range incommensurate spin correlation is observed and analyzed in the doped Hubbard and the t - J model both in one and two dimensions. In the ground state, the superconducting correlation shows the absence of system size dependence in two dimensions at the filling smaller than 0.8 and in one dimension at any filling. The binding energy of two fermions and the spin and charge distortion around an itinerant fermion are also discussed in the t - J and the coupled spin-fermion model.
Exact identities among the nth extended s-wave-pairing amplitudes in the Hubbard model are presented. They provide more constraints on the possible symmetries of the superconducting gap in the Hubbard model than that between the on-site and the first extended s-wave pairings pointed out by Zhang. For the t-J model, we also find an exact identity between the on-site and the first extended s-wave-pairing amplitudes, which is essentially consistent with that in the Hubbard model. Besides the s-wave pairings, we discuss the constraints on the other pairing amplitudes.
We present simple lower bounds on the ground-state energy of the two-dimensional (2D) Hubbard and t-J models for arbitrary values of band filling and coupling constant. For the Hubbard model we derive two types of bounds, both based on decomposing the model Hamiltonian into a sum of sub-Hamiltonians. For a decomposition into local cluster sub-Hamiltonians, we perform a generalized Legendre transform on previously derived bounds for the grand-canonical potential. For a decomposition into spin-up and spin-down parts, previous results on the spinless Falicov-Kimball model may be used to obtain bounds for the Hubbard model, generalizing a result of Langer and Mattis to arbitrary filling. For the 2D t-J model we have only considered the decomposition into clusters. The 1D Hubbard model is used as a test case. The bounds may be improved by diagonalizing the Hamiltonian for larger clusters.