Charge Order in the Holstein Model on a Honeycomb Lattice

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.

A non-perturbative study of the interplay between electron-phonon interaction and Coulomb interaction in undoped graphene

In condensed-matter systems, electrons are subjected to two different interactions under certain conditions. Even if both interactions are weak, it is difficult to perform perturbative calculations due to the complexity caused by the interplay of two interactions. When one or two interactions are strong, ordinary perturbation theory may become invalid. Here we consider undoped graphene as an example and provide a non-perturbative quantum-field-theoretic analysis of the interplay of electron-phonon interaction and Coulomb interaction. We treat these two interactions on an equal footing and derive the exact Dyson-Schwinger (DS) integral equation of the full Dirac-fermion propagator. This equation depends on several complicated correlation functions and thus is difficult to handle. Fortunately, we find that these correlation functions obey a number of exact identities, which allows us to prove that the DS equation of full fermion propagator is self-closed. After solving this self-closed equation, we obtain the renormalized fermion velocity and show that its energy (momentum) dependence of renormalized fermion velocity is dominantly determined by the electron-phonon (Coulomb) interaction. In particular, the renormalized velocity exhibits a logarithmic momentum dependence and a non-monotonic energy dependence.

Chiral Ising Gross-Neveu Criticality of a Single Dirac Cone: A Quantum Monte Carlo Study

We perform large-scale quantum Monte Carlo simulations of SLAC fermions on a two-dimensional square lattice at half filling with a single Dirac cone with N=2 spinor components and repulsive on-site interactions. Despite the presence of a sign problem, we accurately identify the critical interaction strength U_{c}=7.28±0.02 in units of the hopping amplitude, for a continuous quantum phase transition between a paramagnetic Dirac semimetal and a ferromagnetic insulator. Using finite-size scaling, we extract the critical exponents for the corresponding N=2 chiral Ising Gross-Neveu universality class: the inverse correlation length exponent ν^{-1}=1.19±0.03, the order parameter anomalous dimension η_{ϕ}=0.31±0.01, and the fermion anomalous dimension η_{ψ}=0.136±0.005.

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.

Real-time evolution of static electron-phonon models in time-dependent electric fields

We present an exact Monte Carlo method to simulate the nonequilibrium dynamics of electron-phonon models in the adiabatic limit of zero phonon frequency. The classical nature of the phonons allows us to sample the equilibrium phonon distribution and efficiently evolve the electronic subsystem in a time-dependent electromagnetic field for each phonon configuration. We demonstrate that our approach is particularly useful for charge-density-wave systems experiencing pulsed electric fields, as they appear in pump-probe experiments. For the half-filled Holstein model in one and two dimensions, we calculate the out-of-equilibrium response of the current and the energy after a pulse is applied as well as the photoemission spectrum before and after the pump. Finite-size effects are under control for chains of 162 sites (in one dimension) or 16×16 square lattices (in two dimensions).

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.

Quantum Monte Carlo Simulations of the 2D Su-Schrieffer-Heeger Model

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.

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.