Variational and auxiliary field Monte Carlo for the Hubbard and Hubbard-Holstein models: an accurate finite-size scaling and a "sign problem" solution

Understanding the properties of strongly correlated materials is a tough work since most of the times analytical tools cannot be applied due to the complicated nature of the problem. Here we apply quantum Monte Carlo techniques in order to tackle with one of the most important and unknown phenomena of many materials: the mechanism of superconductivity. We study the Hubbard-Holstein model which is one of the most simple theoretical models including strong correlation and electron-phonon coupling. By using twist-averaging boundary conditions we eliminate finite-size errors in the Monte Carlo simulations. With this useful tool we first solve the Hubbard-Holstein model by variational Monte Carlo by reporting the phase diagram of the model for different phonon frequencies at half filling. Then we investigate the phase separation away from half filling. Finally we try to attack the problem with an essentially unbiased method based on the auxiliary field quantum Monte Carlo technique and accelerated Langevin dynamics. By curing the sign problem via Cauchy integration in the complex plane where the auxiliary fields have been defined, we report the effect of the electron-phonon coupling on some observables such as magnetization, spin and charge structure factors. We show preliminary results that are already meaningful to understand the nature of the transition between magnetism and charge order in the model at half filling.