In this thesis we consider the hybrid quantum Monte Carlo method for simulations of the Hubbard and Su-Schrieffer-Heeger model. In the first instance, we discuss the hybrid quantum Monte Carlo method for the Hubbard model on a square lattice. We point out potential ergodicity issues and provide a way to circumvent them by a complexification of the method. Furthermore, we compare the efficiency of the hybrid quantum Monte Carlo method with a well established determinantal quantum Monte Carlo method for simulations of the half-filled Hubbard model on square lattices. One reason why the hybrid quantum Monte Carlo method loses the comparison is that we do not observe the desired sub-quadratic scaling of the numerical effort. Afterwards we present a formulation of the hybrid quantum Monte Carlo method for the Su-Schrieffer-Heeger model in two dimensions. Electron-phonon models like this are in general very hard to simulate using other Monte Carlo methods in more than one dimensions. It turns out that the hybrid quantum Monte Carlo method is much better suited for this model . We achieve favorable scaling properties and provide a proof of concept. Subsequently, we use the hybrid quantum Monte Carlo method to investigate the Su-Schrieffer-Heeger model in detail at half-filling in two dimensions. We present numerical data for staggered valence bond order at small phonon frequencies and an antiferromagnetic order at high frequencies. Due to an O(4) symmetry the antiferromagnetic order is connected to a superconducting charge density wave. Considering the Su-Schrieffer-Heeger model without tight-binding hopping reveals an additional unconstrained Z_2 gauge theory. In this case, we find indications for π-fluxes and a possible Z_2 Dirac deconfined phase as well as for a columnar valence bond ordered state at low phonon energies. In our investigations of the several phase transitions we discuss the different possibilities for the underlying mechanisms and reveal first insights into a rich phase diagram.
In this thesis we consider the hybrid quantum Monte Carlo method for simulations of the Hubbard and Su-Schrieffer-Heeger model. In the first instance, we discuss the hybrid quantum Monte Carlo method for the Hubbard model on a square lattice. We point out potential ergodicity issues and provide a way to circumvent them by a complexification of the method. Furthermore, we compare the efficiency of the hybrid quantum Monte Carlo method with a well established determinantal quantum Monte Carlo method for simulations of the half-filled Hubbard model on square lattices. One reason why the hybrid quantum Monte Carlo method loses the comparison is that we do not observe the desired sub-quadratic scaling of the numerical effort. Afterwards we present a formulation of the hybrid quantum Monte Carlo method for the Su-Schrieffer-Heeger model in two dimensions. Electron-phonon models like this are in general very hard to simulate using other Monte Carlo methods in more than one dimensions. It turns out that the hybrid quantum Monte Carlo method is much better suited for this model . We achieve favorable scaling properties and provide a proof of concept. Subsequently, we use the hybrid quantum Monte Carlo method to investigate the Su-Schrieffer-Heeger model in detail at half-filling in two dimensions. We present numerical data for staggered valence bond order at small phonon frequencies and an antiferromagnetic order at high frequencies. Due to an O(4) symmetry the antiferromagnetic order is connected to a superconducting charge density wave. Considering the Su-Schrieffer-Heeger model without tight-binding hopping reveals an additional unconstrained Z_2 gauge theory. In this case, we find indications for π-fluxes and a possible Z_2 Dirac deconfined phase as well as for a columnar valence bond ordered state at low phonon energies. In our investigations of the several phase transitions we discuss the different possibilities for the underlying mechanisms and reveal first insights into a rich phase diagram.
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 unique feature of the hybrid quantum Monte Carlo (HQMC) method is the potential to simulate negative sign free lattice fermion models with subcubic scaling in system size. Here we will revisit the algorithm for various models. We will show that for the Hubbard model the HQMC suffers from ergodicity issues and unbounded forces in the effective action. Solutions to these issues can be found in terms of a complexification of the auxiliary fields. This implementation of the HQMC that does not attempt to regularize the fermionic matrix so as to circumvent the aforementioned singularities does not outperform single spin flip determinantal methods with cubic scaling. On the other hand we will argue that there is a set of models for which the HQMC is very efficient. This class is characterized by effective actions free of singularities. Using the Majorana representation, we show that models such as the Su-Schrieffer-Heeger Hamiltonian at half filling and on a bipartite lattice belong to this class. For this specific model subcubic scaling is achieved.
Monte Carlo methods have been strong tools in numerous fields. Their simple formulation is easy to implement to many problems and can be adapted from the classical Markov chain Monte Carlo (MCMC) method into any number of related techniques, including the quantum realm. This dissertation first presents three such cases in condensed matter physics. We begin by introducing electron-phonon interactions and building up the phonon models relevant to us in Chapter 2, then covering the groundwork for classical Monte Carlo methods along with how to adapt them for quantum systems in Chapter 3. Chapter 4 utilizes the MCMC formulation to study how “synthetic dimensions” within a classical model can be modeled and simulated. From there, Chapter 5 brings us into quantum systems in order to study the interplay present within a combined Holstein and Su-Schrieffer-Heeger model. Chapter 6 goes through the optical Su-Schrieffer-Heeger model on a triangular lattice to examine the effects of geometric frustration in conjunction with doping from half filling. Finally, Chapter 7 summarizes previous findings and proposes directions for future research, outlining topics of continuing interest.
Most nonperturbative numerical studies of electron-phonon interactions focus on model Hamiltonians where the electrons interact with a phonon branch via a single type of microscopic mechanism. Two commonly explored couplings in this context are the Holstein and Su-Schrieffer-Heeger (SSH) interactions, which describe phonons modulating the on-site energy and intersite electron hopping, respectively. Many materials, however, have multiple phonon branches that can each interact with electronic degrees of freedom in different ways. We present here a determinant quantum Monte Carlo study of the half-filled two-dimensional (bond) SSH-Holstein Hamiltonian, where electrons couple to different phonon branches via either the Holstein or SSH mechanism. We map the model's phase diagram and determine the nature of the transitions between charge-density wave, bond-order wave, and antiferromagnetic order.
Over the past several years, a new generation of quantum simulations has greatly expanded our understanding of charge density wave phase transitions in Hamiltonians with coupling between local phonon modes and the on-site charge density. A quite different, and interesting, case is one in which the phonons live on the bonds, and hence modulate the electron hopping. This situation, described by the Su-Schrieffer-Heeger (SSH) Hamiltonian, has so far only been studied with quantum Monte Carlo in one dimension. Here we present results for the 2D SSH model, show that a bond ordered wave (BOW) insulator is present in the ground state at half filling, and argue that a critical value of the electron-phonon coupling is required for its onset, in contradistinction with the 1D case where BOW exists for any nonzero coupling. We determine the precise nature of the bond ordering pattern, which has hitherto been controversial, and the critical transition temperature, which is associated with a spontaneous breaking of Z_{4} symmetry.
The Hubbard and Su-Schrieffer-Heeger (SSH) Hamiltonians are iconic models for understanding the qualitative effects of electron-electron and electron-phonon interactions, respectively. In the two-dimensional square lattice Hubbard model at half filling, the on-site Coulomb repulsion U between up and down electrons induces antiferromagnetic (AFM) order and a Mott insulating phase. On the other hand, for the SSH model, there is an AFM phase when the electron-phonon coupling $\ensuremath{\lambda}$ is less than a critical value ${\ensuremath{\lambda}}_{c}$ and a bond order wave when $\ensuremath{\lambda}>{\ensuremath{\lambda}}_{c}$. In this Letter, we perform numerical studies on the square lattice optical Su-Schrieffer-Heeger-Hubbard Hamiltonian, which combines both interactions. We use the determinant quantum Monte Carlo method which does not suffer from the fermionic sign problem at half filling. We map out the phase diagram and find that it exhibits a direct first-order transition between an antiferromagnetic phase and a bond-ordered wave as $\ensuremath{\lambda}$ increases. The AFM phase is characterized by two different regions. At smaller $\ensuremath{\lambda}$ the behavior is similar to that of the pure Hubbard model; the other region, while maintaining long-range AFM order, exhibits larger kinetic energies and double occupancy, i.e., larger quantum fluctuations, similar to the AFM phase found in the pure SSH model.