We have calculated within the fluctuation exchange approximation the full momentum and real-frequency dependence of the quasiparticle self-energy, the gap function, and the dynamic spin susceptibility for the two-dimensional Hubbard model near half filling. The spin-density-wave (SDW) instability is reached only asymptotically in the limit of infinite U or zero temperature. Near the SDW instability we find approximately marginal Fermi-liquid behavior, a pseudogap in the density of states, and a singlet
dx2-
y2 pairing instability. Below
Tc, gaps open rapidly in the frequency spectra of the gap function, the quasiparticle scattering rate, and the density of states. The spectral density of the dynamic spin susceptibility agrees qualitatively with the neutron-scattering data on ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$
SrxCuO4. The quasiparticle scattering rate and spin-lattice relaxation rate ${\mathit{T}}_{1}^{\mathrm{\ensuremath{-}}1}$ decrease rapidly below
Tc. The photoemission intensities along the antinode and node of the gap are in good agreement with the angle-resolved photoemission spectroscopy data on Bi 2212. Including next-nearest-neighbor hopping we obtain for sufficiently large coupling (U/t\ensuremath{\gtrsim}7) a broad commensurate neutron-scattering peak like the observed one in
YBa2Cu3O6+x.
Using quantum Monte Carlo techniques we study the magnetic susceptibility of the Hubbard model on a two-dimensional square lattice for values of the coupling U/t ranging from 4 to 20, and as a function of temperature and doping. For U/t=4 the magnetic susceptibility has a maximum at half filling and decreases with doping, in qualitative agreement with recent random-phase-approximation calculations. However, for U/t=10 the susceptibility increases with hole doping near half filling resembling the experimental results obtained for ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$
SrxCuO4. At low electronic density the temperature dependence of the susceptibility is compatible with that of a Fermi liquid.
Quasiparticle bands of the two-dimensional Hubbard model are calculated using the Roth two-pole approximation to the one-particle Green's function. Excellent agreement is obtained with recent Monte Carlo calculations, including an anomalous volume of the Fermi surface near half-filling, which can possibly be explained in terms of a breakdown of Fermi liquid theory. The calculated bands are very flat around the (\ensuremath{\pi},0) points of the Brillouin zone in agreement with photoemission measurements of cuprate superconductors. With doping there is a shift in spectral weight from the upper band to the lower band. The Roth method is extended to deal with superconductivity within a four-pole approximation allowing electron-hole mixing. It is shown that triplet p-wave pairing never occurs. A self-consistent solution with singlet ${\mathit{d}}_{{\mathit{x}}^{2}\mathrm{\ensuremath{-}}{\mathit{y}}^{2}}$-wave pairing is found and optimal doping occurs when the van Hove singularity, corresponding to the flat band part, lies at the Fermi level. Nearest-neighbor antiferromagnetic correlations play an important role in flattening the bands near the Fermi level and in favoring superconductivity. However, the mechanism for superconductivity is a local one, in contrast to spin-fluctuation exchange models. For reasonable values of the hopping parameter the transition temperature
Tc is in the range 10--100 K. The optimum doping ${\mathrm{\ensuremath{\delta}}}_{\mathit{c}}$ lies between 0.14 and 0.25, depending on the ratio U/t. The gap equation has a BCS-like form and 2${\mathrm{\ensuremath{\Delta}}}_{\mathrm{max}}$/
kTc\ensuremath{\simeq}4.
We study the electronic thermal conductivity ${\ensuremath{\kappa}}_{\text{el}}$ and the thermal diffusion constant
DQ,el in the square lattice Hubbard model using the finite-temperature Lanczos method. We make use of the Nernst-Einstein relation for thermal transport and connect the strong nonmonotonous temperature dependence of ${\ensuremath{\kappa}}_{\text{el}}$ to that of
DQ,el and the electronic specific heat
cel. We present also the results for the Heisenberg model on a square lattice and ladder geometries. We study the effects of doping and consider the doped case also with the dynamical mean-field theory. We show that ${\ensuremath{\kappa}}_{\text{el}}$ is below the corresponding Mott-Ioffe-Regel value in almost all calculated regimes, while the mean free path is typically above or close to lattice spacing. We discuss the opposite effect of quasiparticle renormalization on charge and heat diffusion constants. We calculate the Lorenz ratio and show that it differs from the Sommerfeld value. We discuss our results in relation to experiments on cuprates. Additionally, we calculate the thermal conductivity of overdoped cuprates within the anisotropic marginal Fermi liquid phenomenological approach.
Monte Carlo results for a two-dimensional Hubbard model in the intermediate-coupling regime U=4t are compared with a diagrammatic spin-fluctuation approximation. Simulations on an 8\ifmmode\times\else\texttimes\fi{}8 lattice doped away from half filling were carried out down to temperatures of order 140 of the bandwidth. Results for the spin susceptibility \ensuremath{\chi}(q,i\ensuremath{\omega}
m), the electron self-energy \ensuremath{\Sigma}(p,i\ensuremath{\omega}
n), various pair-field susceptibilities, and the irreducible particle-particle scattering vertex \ensuremath{\Gamma}(p\ensuremath{'},i\ensuremath{\omega}${}_{\mathrm{n}\ensuremath{'}}$|p,i\ensuremath{\omega}
n) were obtained. A random-phase approximation for \ensuremath{\chi}(q,i\ensuremath{\omega}
m) with a renormalized Coulomb coupling U\ifmmode\bar\else\textasciimacron\fi{} is shown to provide a fit to the Monte Carlo data. A similar approximation for the Berk-Schrieffer spin-fluctuation interaction also provides a reasonable fit to the self-energy \ensuremath{\Sigma}(p,i\ensuremath{\omega}
n) in the region explored by the Monte Carlo data. However, a similar approximation for the irreducible particle-particle interaction failed to reproduce the Monte Carlo results. Higher-order vertex corrections were calculated, but significant discrepancies with Monte Carlo results for \ensuremath{\Gamma} remain.
The interplay between antiferromagnetism and superconductivity is studied by using the three-dimensional nearly half-filled Hubbard model with anisotropic transfer matrices
tz and ${t}_{\ensuremath{\perp}}.$ The phase diagrams are calculated for varying values of the ratio ${r}_{z}{=t}_{z}{/t}_{\ensuremath{\perp}}$ using the spin fluctuation theory within the fluctuation-exchange approximation. The antiferromagnetic phase around the half-filled electron density expands while the neighboring phase of the anisotropic ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity shrinks with increasing
rz. For small
rz Tc decreases slowly with increasing
rz. For moderate values of
rz we find a second order transition, with lowering temperature, from the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconducting phase to a phase where an incommensurate spin density wave (SDW) coexists with ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity. Resonance peaks discussed previously for two-dimensional (2D) superconductors are shown to survive in the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconducting phase of 3D systems. Soft components of the incommensurate SDW spin fluctuation mode grow as the coexistent phase is approached.