We study the three-band Hubbard model for the cuprate superconductors using Quantum Monte Carlo simulations and the fluctuation exchange approximation including a self-consistent many-body renormalization of the single-particle propagator. The energy dispersion of the low-lying one-particle excitations is in very good agreement with spectroscopic measurements for a hole doping concentration δ = 0.25 of the CuO2 planes. At the same doping, we examine the spin susceptibility yielding a non-ordered phase with short-range incommensurate fluctuations near the antiferromagnetic wave vector. In the superconducting regime, the interaction vertex for the pairing correlation function in the extended s-wave channel shows a maximum for δ = 0.20, resembling the dependence of transition temperature Tc on doping.
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.
We present recent theoretical results on superconductivity in correlated-electron systems, especially in the two-dimensional Hubbard model and the three-band d-p model. The mechanism of superconductivity in high-temperature superconductors has been extensively studied on the basis of various electronic models and also electron-phonon models. In this study, we investigate the properties of superconductivity in correlated-electron systems by using numerical methods such as the variational Monte Carlo method and the quantum Monte Carlomethod. The Hubbard model is one of basic models for strongly correlated electron systems, and is regarded as the model of cuprate high temperature superconductors. The d-p model is more realistic model for cuprates. The superconducting condensation energy obtained by adopting the Gutzwiller ansatz is in reasonable agreement with the condensation energy estimated for YBa2Cu3O7. We show the phase diagram of the ground state using this method. We have further investigated the stability of striped and checkerboard states in the under-doped region. Holes doped in a half-filled square lattice lead to an incommensurate spin and charge density wave. The relationship of the hole density x and incommensurability δ, δ~x, is satisfied in the lower doping region, as indicated by the variationalMonte Carlocalculations for the two-dimensional Hubbard model. A checkerboard-like charge-density modulation with a roughly period has also been observed by scanning tunneling microscopy experiments in Bi2212 and Na-CCOC compounds. We have performed a variational Monte Carlo simulation on a two-dimensional t-t′-t″- U Hubbard model with a Bi-2212 type band structure and found that the period checkerboard spin modulation, that is characterized by multi Q vectors, is indeed stabilized. We have further performed an investigation by using a quantumMonte Carlomethod, which is a numerical method that can be used to simulate the behavior of correlated electron systems. We present a new algorithm of the quantum Monte Carlo diagonalization that is a method for the evaluation of expectation value without the negative sign problem. We compute pair correlation functions and show that pair correlation is indeed enhanced with hole doping.
Quantum Monte Carlo (QMC) and maximum-entropy techniques are used to study the spectral function $A(\mathbf{p},\ensuremath{\omega})$ of the one-band Hubbard model with strong coupling including a next-nearest-neighbor electronic hopping with amplitude ${t}^{\ensuremath{'}}/t=\ensuremath{-}0.35$. These values of parameters are chosen to improve the comparison of the Hubbard model with angle-resolved photoemission (ARPES) data for
Sr2CuO2Cl2. A narrow quasiparticle (qp) band is observed in the QMC analysis at the temperature of the simulation
T=t/3, both at and away from half-filling. Such a narrow band produces a large accumulation of weight in the density of states at the top of the valence band. As the electronic density
〈n〉 decreases further away from half-filling, the chemical potential travels through this energy window with a large number of states, and by $〈n〉\ensuremath{\sim}0.70$ it has crossed it entirely. The region near momentum $(0,\ensuremath{\pi})$ and $(\ensuremath{\pi},0)$ in the spectral function is more sensitive to doping than momenta along the diagonal from
(0,0) to $(\ensuremath{\pi},\ensuremath{\pi})$. The evolution with hole density of the quasiparticle dispersion contains some of the features observed in recent ARPES data in the underdoped regime. For sufficiently large hole densities the ``flat'' bands at $(\ensuremath{\pi},0)$ cross the Fermi energy, a prediction that could be tested with ARPES techniques applied to overdoped cuprates. The population of the qp band introduces a hidden density in the system which produces interesting consequences when the quasiparticles are assumed to interact through antiferromagnetic fluctuations and studied with the BCS gap-equation formalism. In particular, a region of extended
s-wave character is found to compete with the
d wave in the overdoped regime, i.e., when the chemical potential has almost entirely crossed the qp band as
〈n〉 is reduced. The present study also shows that previous ``real-space'' pairing theories for the cuprates, such as the antiferromagnetic Van Hove scenario, originally constructed based on information gathered at half-filling, do not change their predictions if hole dispersions resembling noninteracting electrons with renormalized parameters are used.
To gain a deeper understanding of the high-temperature superconductivity in doped cuprates, the constrained-path quantum Monte Carlo method was used to explore the superconducting pairing correlation function and effective pairing interaction in the three-band Hubbard model. Our main observation was that in the low doping region, compared to other symmetries, [Formula: see text]-wave superconducting pairing dominates. The Coulomb interaction [Formula: see text] suppresses [Formula: see text]-wave pairing correlation, while the effective pairing interaction is both enhanced by the increasing [Formula: see text] and the nearest attraction. These results are used to understand the superconducting state in copper oxide high-temperature superconductors.