We present in this article results of Quantum Monte Carlo simulations for the three-band Hubbard model obtained both with the grand canonical approach as well as with the T =0 formalism. A direct comparison of Quantum Monte Carlo results with experimental values for the charge-transfer gap leads to a set of parameters for the three-band Hubbard model that consistently describes salient features of the high- T c superconductors. We find long-range antiferromagnetic order in the undoped case as well as an incommensurate structure for a doping concentration of δ≈0.20. Moreover the magnetic susceptibility as a function of doping δ is in a qualitative agreement with experiments. Concerning dynamical properties, we have evaluated the spectral density [Formula: see text]. The dispersion of the state with the lowest excitation energy is in very good agreement with experimental photoemission and inverse photoemission experiments for a doping concentration of δ=0.25. In the superconducting regime, we find that the interaction vertex for the pairing correlation function in the extended s-channel shows a maximum for δ~0.20, resembling the dependence of T c on doping. This feature is not present for other symmetries of the order parameter or in other parameter regions.
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.
The concentration dependences of the band structure, spectral weight, density of states, and Fermi surface in the paramagnetic state are studied in the Hubbard model within cluster pertubation theory with 2 × 2 clusters. Representation of the Hubbard X operators makes it possible to control conservation of the spectral weight in constructing cluster perturbation theory. The calculated value of the ground-state energy is in good agreement with the results obtained using nonperturbative methods such as the quantum Monte Carlo method, exact diagonalization of a 4 × 4 cluster, and the variational Monte Carlo method. It is shown that in the case of hole doping, the states in the band gap (in-gap states) lie near the top of the lower Hubbard band for large values of U and near the bottom of the upper band for small U . The concentration dependence of the Fermi surface strongly depends on hopping to second ( t ′) and third ( t ″) neighbors. For parameter values typical of HTSC cuprates, the existence of three concentration regions with different Fermi surfaces is demonstrated. It is shown that broadening of the spectral electron density with an energy resolution typical of contemporary ARPES leads to a pattern of arcs with a length depending on the concentration. Only an order-of-magnitude decrease in the linewidth makes it possible to obtain the true Fermi surface from the spectral density. The kinks associated with strong electron correlations are detected in the dispersion relation below the Fermi level.
We compute high-resolution angle-resolved photoemission spectroscopy of the Hubbard model using the unbiased determinant quantum Monte Carlo algorithm, revealing an asymmetry between electron and hole doping. Electron doping exhibits more coherent quasiparticles and stronger antiferromagnetic correlations compared to hole doping. At low doping, a nodal-antinodal dichotomy on the Fermi surface is observed, similar to cuprate experiments. The dichotomy reflects the momentum dependence of the Mott gap, as manifested in both the spectral function and the self-energy. For hole doping, we observe a transition towards the pseudogap, without signature of pocket formation. The simulated nuclear magnetic resonance pseudogap temperatures do not necessarily agree with the temperature determined by spectroscopy. These findings collectively suggest the pseudogap is a smooth crossover driven by strong correlations.
We present results of a systematic quantum Monte Carlo study for the single-band Hubbard model. Thereby we evaluated single-particle spectra (PES and IPES), two-particle spectra (spin and density correlation functions), and the dynamical correlation function of suitably defined diagnostic operators, all as a function of temperature and hole doping. The results allow us to identify different physical regimes. Near half-filling we find an anomalous ``Hubbard-I phase,'' where the band structure is, up to some minor modifications, consistent with the Hubbard-I predictions. At lower temperatures, where the spin response becomes sharp, additional dispersionless ``bands'' emerge due to the dressing of electrons/holes with spin excitations. We present a simple phenomenological fit that reproduces the band structure of the insulator quantitatively. The Fermi surface volume in the low-doping phase, as derived from the single-particle spectral function, is not consistent with the Luttinger theorem, but qualitatively in agreement with the predictions of the Hubbard-I approximation. The anomalous phase extends up to a hole concentration of $\ensuremath{\approx}15%,$ i.e., the underdoped region in the phase diagram of high-
Tc superconductors. We also investigate the nature of the magnetic ordering transition in the single-particle spectra. We show that the transition to a spin-density wave-like band structure is not accomplished by the formation of any resolvable ``precursor bands,'' but rather by a (spectroscopically invisible) band of spin-3/2 quasiparticles. We discuss implications for the ``remnant Fermi surface'' in insulating cuprate compounds and the shadow bands in the doped materials.
A detailed numerical study of the three-band Hubbard (3BH) model in the underdoped region is carried out, both in the hole- as well as in the electron-doped case, by means of the variational cluster approach. The phase diagram and the low-energy single-particle spectrum are found to be very similar to recent results for the single-band Hubbard model with next-nearest-neighbor hoppings t '. In particular, we obtain a mixed antiferromagnetic+superconducting phase at low doping with a first-order transition to a pure superconducting phase accompanied by phase separation. Despite the fact that the physically relevant and generally accepted parameters ( t pp / t pd , t pd /Δ pd ) of our 3BH-model are not in the small regime, required for a 'downward' projection onto a single-band model, in the single-particle spectrum a clear Zhang–Rice (ZR) singlet band (with a projection on d- and p-orbitals consistent with the phase factor of the ZR singlet wave function) can be seen. It displays an incoherent and a coherent part, in which holes enter upon doping around (π/2, π/2). The latter is very similar to the coherent quasi-particle band crossing the Fermi surface in the t – t '– U single-band model, with generally accepted values of the parameters. Doped electrons go instead into the upper Hubbard band, first filling the regions of the Brillouin zone around (π, 0). This fact can be related to the enhanced robustness of the antiferromagnetic phase as a function of electron doping compared to hole doping.