Phonon Dispersion and the Competition between Pairing and Charge Order

The impact of Rashba spin-orbit coupling in charge-ordered systems

We study the impact of the Rashba spin-orbit coupling (RSOC) on the stability of charge-density wave (CDW) in systems with large electron-phonon coupling (EPC). Here, the EPC is considered in the framework of the Holstein model at the half-filled square lattice. We obtain the phase diagram of the Rashba-Holstein model using the Hartree-Fock mean-field theory, and identify the boundaries of the CDW and Rashba metal phases. We notice that the RSOC disfavors the CDW phase, driving the system to a correlated Rashba metal. Also, we employ a cluster perturbation theory (CPT) approach to investigate the phase diagram beyond the Hartree-Fock approximation. The quantum correlations captured by CPT indicate that the RSOC is even more detrimental to CDW than previously anticipated. That is, the Rashba metal region is observed to be expanded in comparison to the mean-field case. Additionally, we investigate pairing correlations, and the results further strengthen the identification of critical points.

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

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.

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

The origin of the pseudogap behavior, found in many high-T_c superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a “gap-filling" rather than “gap-closing" behavior, similar to the one observed in cuprate superconductors. Our results present direct evidence of the pseudogap state, driven by superconducting fluctuations.

The origin of pseudogap in high-T_c superconductors remains a big puzzle. Here, the authors report numerical evidence of pseudogap behavior employing Quantum Monte Carlo algorithm emerging from pairing fluctuations in a quantum-critical non-Fermi liquid, similar to the pseudogap phase observed in cuprate superconductors.

Many-Body Quantum State Diffusion for Non-Markovian Dynamics in Strongly Interacting Systems

Capturing non-Markovian dynamics of open quantum systems is generally a challenging problem, especially for strongly interacting many-body systems. In this Letter, we combine recently developed non-Markovian quantum state diffusion techniques with tensor network methods to address this challenge. As a first example, we explore a Hubbard-Holstein model with dissipative phonon modes, where this new approach allows us to quantitatively assess how correlations spread in the presence of non-Markovian dissipation in a 1D many-body system. We find regimes where correlation growth can be enhanced by these effects, offering new routes for dissipatively enhancing transport and correlation spreading, relevant for both solid state and cold atom experiments.

Antiferromagnetism Induced by Bond Su-Schrieffer-Heeger Electron-Phonon Coupling: A Quantum Monte Carlo Study

Antiferromagnetism (AFM) such as Néel ordering is often closely related to Coulomb interactions such as Hubbard repulsion in two-dimensional (2D) systems. Whether Néel AFM ordering in two dimensions can be dominantly induced by electron-phonon couplings (EPC) has not been completely understood. Here, by employing numerically exact sign-problem-free quantum Monte Carlo (QMC) simulations, we show that bond Su-Schrieffer-Heeger (SSH) phonons with frequency ω and EPC constant λ can induce AFM ordering for a wide range of phonon frequency ω>ω_{c}. For ω<ω_{c}, a valence-bond-solid (VBS) order appears and there is a direct quantum phase transition between VBS and AFM phases at ω_{c}. The phonon mechanism of the AFM ordering is related to the fact that SSH phonons directly couple to electron hopping whose second-order process can induce an effective AFM spin exchange. Our results shall shed new light on understanding AFM ordering in correlated quantum materials.

Reducing autocorrelation time in determinant quantum Monte Carlo using the Wang-Landau algorithm: Application to the Holstein model

When performing a Monte Carlo calculation, the running time should, in principle, be much longer than the autocorrelation time in order to get reliable results. Among different lattice fermion models, the Holstein model is notorious for its particularly long autocorrelation time. In this paper, we employ the Wang-Landau algorithm in the determinant quantum Monte Carlo to achieve the flat-histogram sampling in the "configuration weight space," which can greatly reduce the autocorrelation time by sacrificing some sampling efficiency. The proposal is checked in the Holstein model on both square and honeycomb lattices. Based on such a Wang-Landau assisted determinant quantum Monte Carlo method, some models with long autocorrelation times can now be simulated possibly.

Magnetism and Charge Order in the Honeycomb Lattice

Despite being relevant to better understand the properties of honeycomblike systems, as graphene-based compounds, the electron-phonon interaction is commonly disregarded in theoretical approaches. That is, the effects of phonon fields on interacting Dirac electrons is an open issue, in particular when investigating long-range ordering. Thus, here we perform unbiased quantum Monte Carlo simulations to examine the Hubbard-Holstein model (HHM) in the half-filled honeycomb lattice. By performing careful finite-size scaling analysis, we identify semimetal-to-insulator quantum critical points, and determine the behavior of the antiferromagnetic and charge-density wave phase transitions. We have, therefore, established the ground state phase diagram of the HHM for intermediate interaction strength, determining its behavior for different phonon frequencies. Our findings provide quantitative and qualitative descriptions of the model at intermediate coupling strengths, and may shed light on the emergence of many-body properties in honeycomblike systems.

Revealing fermionic quantum criticality from new Monte Carlo techniques

This review summarizes recent developments in the study of fermionic quantum criticality, focusing on new progress in numerical methodologies, especially quantum Monte Carlo methods, and insights that emerged from recently large-scale numerical simulations. Quantum critical phenomena in fermionic systems have attracted decades of extensive research efforts, partially lured by their exotic properties and potential technology applications, and partially awakened by the profound and universal fundamental principles that govern these quantum critical systems. Due to the complex and non-perturbative nature, these systems face the most difficult and challenging problems in the study of modern condensed matter physics, and many important fundamental problems remain open. Recently, new developments in model design and algorithm improvements enabled unbiased large-scale numerical solutions to be achieved in the close vicinity of these quantum critical points, which paves a new pathway towards achieving controlled conclusions through combined efforts of theoretical and numerical studies, as well as possible theoretical guidance for experiments in heavy-fermion compounds, Cu-based and Fe-based superconductors, ultra-cold fermionic atomic gas, twisted graphene layers, etc, where signatures of fermionic quantum criticality exist.

Charge-Density-Wave Transitions of Dirac Fermions Coupled to Phonons

The spontaneous generation of charge-density-wave order in a Dirac fermion system via the natural mechanism of electron-phonon coupling is studied in the framework of the Holstein model on the honeycomb lattice. Using two independent and unbiased quantum Monte Carlo methods, the phase diagram as a function of temperature and coupling strength is determined. It features a quantum critical point as well as a line of thermal critical points. Finite-size scaling appears consistent with fermionic Gross-Neveu-Ising universality for the quantum phase transition and bosonic Ising universality for the thermal phase transition. The critical temperature has a maximum at intermediate couplings. Our findings motivate experimental efforts to identify or engineer Dirac systems with sufficiently strong and tunable electron-phonon coupling.

Charge Order in the Holstein Model on a Honeycomb Lattice

The effect of electron-electron interactions on Dirac fermions, and the possibility of an intervening spin-liquid phase between the semimetal and antiferromagnetic (AF) regimes, has been a focus of intense quantum simulation effort over the last five years. We use determinant quantum Monte Carlo simulations to study the Holstein model on a honeycomb lattice and explore the role of electron- phonon interactions on Dirac fermions. We show that they give rise to charge-density-wave (CDW) order and present evidence that this occurs only above a finite critical interaction strength. We evaluate the temperature for the transition into the CDW which, unlike the AF transition, can occur at finite values owing to the discrete nature of the broken symmetry.