Pairing, Charge, and Spin Correlations in the Three-Band Hubbard Model
  Using the Constrained Path Monte Carlo (CPMC) method, we simulated the two-dimensional, three-band Hubbard model to study pairing, charge, and spin correlations as a function of electron and hole doping and the Coulomb repulsion $V_{pd}$ between charges on neighboring Cu and O lattice sites. As a function of distance, both the $d_{x^2 - y^2}$-wave and extended s-wave pairing correlations decayed quickly. In the charge-transfer regime, increasing $V_{pd}$ decreased the long-range part of the correlation functions in both channels, while in the mixed-valent regime, it increased the long-range part of the s-wave behavior but decreased that of the d-wave behavior. Still the d-wave behavior dominated. At a given doping, increasing $V_{pd}$ increased the spin-spin correlations in the charge-transfer regime but decreased them in the mixed-valent regime. Also increasing $V_{pd}$ suppressed the charge-charge correlations between neighboring Cu and O sites. Electron and hole doping away from half-filling was accompanied by a rapid suppression of anti-ferromagnetic correlations. 
Z. B. Huang, H. Q. Lin, J. E. Gubernatis