Robust fluctuating intertwined spin and charge stripes in Emery model

In the single-band Hubbard model, recent numerical exact methods have indicated that charge density wave competes closely or dominates over superconductivity. However, a full picture of charge and spin stripes in the Emery model is still lacking. In this talk, we will present determinant quantum Monte Carlo simulations that show a tendency towards intertwined charge stripes on oxygen orbitals lying on the antiferromagnetic domain walls for spin correlation on copper orbitals. This orbital character of intertwined charge stripe is consistent with the density matrix renormalization group ground state on 4-leg ladders, where the charge modulation is significantly stronger on $p_y$ oxygen orbitals than on copper orbitals, following the Yamada relation. This tendency of charge-spin-locking is already visible at temperatures as high as $T\sim J$, where $J$ is the magnetic exchange parameter. Both methods show that the tendency towards intertwined charge stripe is enhanced for a set of Emery model parameter with larger Hubbard interactions, possibly due to the enhanced local magnetic moments.