Systematic study of intertwined orders in the three-band Hubbard model

The three-band Hubbard model (also known as the Emery model) captures much of the physics of the copper oxide plane of cuprate high-temperature superconductors. We report a numerically exact determinant quantum Monte Carlo study of the three-band Hubbard model for various doping levels on square clusters and four-leg cylinders, focusing on the superconducting pair-field susceptibility, as well as the spin and charge correlation functions. The pair-field susceptibility indicates that d-wave pairing of holes on copper is dominant over an extended s-wave pairing of holes on oxygen and other copper-oxygen pairing channels. Our results shed light on the role of oxygen in multiple intertwined orders of the multi-orbital Hubbard model.