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.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of
U/t from 2 to 8, for a range of densities including half-filling and
n=0.3,0.5,0.6,0.75, and
0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
This thesis investigates the interplay between electron-electron and electron-phonon interactions in low-dimensional lattice models, specifically within the Hubbard model incorporating Su-Schrieffer-Heeger (SSH) coupling with phonons. SSH interactions, which modulate the electron hopping amplitude through lattice displacements, have recently attracted renewed interest due to their potential relevance in strongly correlated materials. The main objective of this work is to determine whether the inclusion of SSH couplings in the Hubbard model can stabilize superconducting or spin-gapped metallic states, as well as other exotic phases of matter. While the half-filled SSH-Hubbard model has been extensively studied in one dimension and only partially examined in two dimensions, its phase diagram in the doped regime remains largely unexplored. This thesis aims to fill this gap, focusing on the emergence of correlated phases away from half-filling in both one and two spatial dimensions. We use the variational Monte Carlo method with a trial wave function defined in a hybrid electron-phonon Hilbert space. In addition to Jastrow factors that account for correlation effects, the wave function includes an antisymmetric component that captures backflow-inspired electron-phonon correlations. This approach enables us to accurately model both electron-electron and electron-phonon interactions. Moreover, since the variational method avoids the sign problem, we can simulate the system in the doped regime, an area typically inaccessible to exact quantum Monte Carlo methods (QMC), which are limited to half-filling. In one dimension, we observe the stabilization of a Luther-Emery liquid, a metallic phase with gapless charge excitations and a finite spin gap, upon doping a Peierls insulating phase. This finding resolves a long-standing debate regarding its realization in SSH chains. Benchmarking with density-matrix renormalization group results confirms the reliability of our approach. In two dimensions, we investigate the SSH-Hubbard model on the square lattice, obtaining results that are consistent with previous exact QMC studies at half-filling. Upon doping, we observe that SSH phonons can induce s-wave superconductivity in the weakly correlated regime. However, we suggest that the enhancement of d-wave superconductivity by SSH coupling in the strongly correlated regime may significantly depend on the specific implementation of the SSH interaction. In one particular case, we find that d-wave superconductivity is slightly suppressed by SSH coupling. These results are further confirmed by a Hartree-Fock analysis performed in the anti-adiabatic limit. Overall, this thesis highlights the crucial role of lattice effects in correlated electron systems, shedding light on the distinct contributions of SSH-type couplings in one and two dimensions. Our findings provide new insights into the microscopic mechanisms behind exotic metallic phases and set the stage for further studies of phonon-mediated phenomena in correlated systems.
A systematically improvable wave function is proposed for the numerical solution of strongly correlated systems. With a stochastic optimization method, based on the auxiliary field quantum Monte Carlo technique, an effective temperature Teff is defined, probing the distance of the ground state properties of the model in the thermodynamic limit from the ones of the proposed correlated mean-field ansatz. In this way their uncertainties from the unbiased zero temperature limit may be estimated by simple and stable extrapolations well before the so called sign problem gets prohibitive. At finite Teff the convergence of the energy to the thermodynamic limit is indeed shown to be possible in the Hubbard model already for relatively small square lattices with linear dimension L ~10, thanks to appropriate averages over several twisted boundary conditions. Within the estimated energy accuracy of the proposed variational ansatz, two clear phases are identified, as the energy is lowered by spontaneously breaking some symmetries satisfied by the Hubbard Hamiltonian: a) a stripe phase where both spin and translation symmetries are broken, and b) a strong coupling d-wave superconducting phase when the particle number is not conserved and global U(1) symmetry is broken. On the other hand the symmetric phase is stable in a wide region at large doping and small coupling.
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.