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.
By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping
t, the on-site electron-electron interaction
U, the phonon energy ${\ensuremath{\omega}}_{0}$, and the electron-phonon coupling
g. At half filling, the ground state is an antiferromagnetic insulator for $U\ensuremath{\gtrsim}2{g}^{2}/{\ensuremath{\omega}}_{0}$, while it is a charge-density-wave (or bipolaronic) insulator for $U\ensuremath{\lesssim}2{g}^{2}/{\ensuremath{\omega}}_{0}$. In addition to these phases, we find a superconducting phase that intrudes between them. For ${\ensuremath{\omega}}_{0}/t=1$, superconductivity emerges when both
U/t and $2{g}^{2}/t{\ensuremath{\omega}}_{0}$ are small; then, by increasing the value of the phonon energy ${\ensuremath{\omega}}_{0}$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case.
The renormalization group technique is applied to one-dimensional electron-phonon Hubbard models at half filling and zero temperature. For the Holstein-Hubbard model, the results of one-loop calculations are congruent with the phase diagram obtained by quantum Monte Carlo simulations in the
(U,gph) plane for the phonon-mediated interaction
gph and the Coulomb interaction
U. The incursion of an intermediate phase between a fully gapped charge-density-wave state and a Mott antiferromagnet is supported along with the growth of its size with the molecular phonon frequency ${\ensuremath{\omega}}_{0}$. We find additional phases enfolding the base boundary of the intermediate phase. A Luttinger liquid line is found below some critical ${U}^{*}\ensuremath{\approx}{g}_{\mathrm{ph}}^{*}$, followed at larger $U\ensuremath{\sim}{g}_{\mathrm{ph}}$ by a narrow region of bond-order-wave ordering which is either charge or spin gapped depending on
U. For the Peierls-Hubbard model, the region of the
(U,gph) plane with a fully gapped Peierls-bond-order-wave state shows a growing domination over the Mott gapped antiferromagnet as the Debye frequency ${\ensuremath{\omega}}_{D}$ decreases. A power-law dependence ${g}_{\mathrm{ph}}\ensuremath{\sim}{U}^{2\ensuremath{\eta}}$ is found to map out the boundary between the two phases, whose exponent is in good agreement with the existing quantum Monte Carlo simulations performed when a finite nearest-neighbor repulsion term
V is added to the Hubbard interaction.
Phase diagram of the two-dimensional Hubbard-Holstein model: enhancement of
s-wave pairing between charge and magnetic orders
We investigate the role of electron-electron and electron-phonon interactions in strongly correlated systems by performing unbiased quantum Monte Carlo simulations in the square lattice Hubbard-Holstein model at half-filling. We study the competition and interplay between antiferromagnetism (AFM) and charge-density wave (CDW), establishing its very rich phase diagram. In the region between AFM and CDW phases, we have found an enhancement of superconducting pairing correlations, favoring (nonlocal)
s-wave pairs. Our study sheds light over past inconsistencies in the literature, in particular the emergence of CDW in the pure Holstein model case.
We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.
The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling
g to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction
U drives the fermions into a Mott antiferromagnet. At low
g and
U, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range $\ensuremath{\xi}$, interpolating between the Holstein and Fr\"ohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range $\ensuremath{\xi}$ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction
U will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the $U\ensuremath{\ne}0$ case, even for small values of the coupling range.