Using numerical diagonalization of a $4\ifmmode\times\else\texttimes\fi{}4$ cluster, we calculate on-site s, extended-
s, and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ pairing correlation functions (PCF's) in an effective generalized Hubbard model for the cuprates, with nearest-neighbor correlated hopping and next-nearest-neighbor hopping ${t}^{\ensuremath{'}}.$ The vertex contributions to the PCF's are significantly enhanced, relative to the ${t\ensuremath{-}t}^{\ensuremath{'}}\ensuremath{-}U$ model. The behavior of the PCF's and their vertex contributions, and signatures of anomalous flux quantization, indicate superconductivity in the d-wave channel for moderate doping and in the s-wave channel for high doping and small U.
With the high-Tc cuprates in mind, properties of correlated d-wave superconducting (SC) states are studied for a Hubbard model on square lattices with a diagonal transfer (t-t’-U model), using a variational Monte Carlo method. We employ a simple wave function, which includes crucial parameters, in particular, a doublon-holon (D-H) binding factor important for correlated SC and normal states as doped Mott insulators. We first check that the range of dominant superconductivity is limited to a strongly correlated regime (U>W, U: onsite correlation strength, W: band width), and coincides with the effective range of the D-H binding factors. In this range of U/t and δ (doping rate), holons (in hole-doped cases) are classified into two types: doped holons and ones created as D-H pairs. Only the former holons participate in current, whereas the latter contribute to singlet-pair formation. Next, we show that the SC properties undergo a crossover at U=Uco∼W. For U<Uco, the behavior is appropriately interpreted with the conventional BCS mechanism, whereas for U>Uco (the regime of cuprates), a new idea is needed to understand the peculiar SC behavior.
We study the square-lattice extended Hubbard model with on-site
U and nearest-neighbor
V interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive
U, charge density wave (extended
s∗-wave superconductivity) at large positive (negative)
V, and
s-wave superconductivity at large negative
U with vanishing
V. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic
p-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by
37.5 holes, especially in the interaction parameter range relevant to the cuprate superconductors.
While both the hole and electron doped cuprates can exhibit ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity, the local distribution of the doped carriers is known to be significantly different with the doped holes going primarily on the O sites while the doped electrons go on the Cu sites. Here we report the results of density-matrix-renormalization-group calculations for a three-orbital model of a
CuO2 lattice. In addition to the asymmetric dependence of the intra-unit-cell occupation of the Cu and O for hole and electron doping, we find important differences in the longer range spin and charge correlations. As expected, the pair-field response has a ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-like structure for both the hole and electron doped systems.
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.