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.
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.
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.
Over the past several decades, computational approaches to studying strongly-interacting systems have become increasingly varied and sophisticated. This book provides a comprehensive introduction to state-of-the-art quantum Monte Carlo techniques relevant for applications in correlated systems. Providing a clear overview of variational wave functions, and featuring a detailed presentation of stochastic samplings including Markov chains and Langevin dynamics, which are developed into a discussion of Monte Carlo methods. The variational technique is described, from foundations to a detailed description of its algorithms. Further topics discussed include optimisation techniques, real-time dynamics and projection methods, including Green's function, reptation and auxiliary-field Monte Carlo, from basic definitions to advanced algorithms for efficient codes, and the book concludes with recent developments on the continuum space. Quantum Monte Carlo Approaches for Correlated Systems provides an extensive reference for students and researchers working in condensed matter theory or those interested in advanced numerical methods for electronic simulation.
This work deals with the development and application of novel quantum Monte Carlo methods to simulate fermion-boson models. Our developments are based on the path-integral formalism, where the bosonic degrees of freedom are integrated out exactly to obtain a retarded fermionic interaction. We give an overview of three methods that can be used to simulate retarded interactions. In particular, we develop a novel quantum Monte Carlo method with global directed-loop updates that solves the autocorrelation problem of previous approaches and scales linearly with system size. We demonstrate its efficiency for the Peierls transition in the Holstein model and discuss extensions to other fermion-boson models as well as spin-boson models. Furthermore, we show how with the help of generating functionals bosonic observables can be recovered directly from the Monte Carlo configurations. This includes estimators for the boson propagator, the fidelity susceptibility, and the specific heat of the Holstein model. The algorithmic developments of this work allow us to study the specific heat of the spinless Holstein model covering its entire parameter range. Its key features are explained from the single-particle spectral functions of electrons and phonons. In the adiabatic limit, the spectral properties are calculated exactly as a function of temperature using a classical Monte Carlo method and compared to results for the Su-Schrieffer-Heeger model.
Phase transitions have been an active area of research in statistical mechanics for almost a century and have recently been integrated into quantum mechanics. Many phenomena such as superconductivity and unconventional magnetism are understood to arise from exotic quantum phases and at points describing quantum phase transitions. A detailed understanding of these phase transitions requires numerical simulations of models which benchmark realistic models against theoretical frameworks. The topic of this thesis is the implementation of Quantum Monte Carlo simulation, which is a powerful technique to understand quantum condensed matter, in interesting models to illustrate novel phenomena in magnetic systems. The novel features of condensed matter systems described in this thesis consist of emergent symmetries at critical points, interesting dynamical features of such systems and the drastic effects of defects in spin systems used in the field of adiabatic quantum computing. Emergent symmetries are shown by condensed matter systems especially at critical points and are features which cannot be shown by individual or a small number of spins. Examples of this in one and two dimensions are presented in an early chapter of this thesis. In addition to this, spin systems can show excitations which have an interesting spatial structure as a consequence of restricted dynamics which only allow the excitations to spread in a particular region. This is presented in the context of a simple model in the following chapter along with numerical support. The following chapter contains a description of adiabatic quantum computing along with a particular model which we study. The phase transition and the effects on the performance of adiabatic quantum computing are studied in this context.
We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and we report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but it incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.