The formation of charge-density-wave order in Dirac fermion systems via electron-phonon coupling represents a significant topic in condensed matter physics. In this work, we investigate this phenomenon within the Holstein model on the honeycomb lattice, with a specific focus on the effect of disorder. While the interplay between electron-electron interactions and disorder has long been a central theme in the field, recent attention has increasingly turned to the combined influence of disorder and electron-phonon coupling. Using determinant quantum Monte Carlo simulations, we concentrate on the phase transitions of charge-density-wave order on the honeycomb lattice. Disorder is introduced through the random hopping of electrons in the system, which can localize electrons via the Anderson effect. Our primary result is that disorder suppresses the charge-density-wave phase, and the interplay between disorder and electron-phonon interactions extends the phase area. We also determine the transition temperature (β_c) to the ordered phase as a function of the electron-phonon coupling. Additionally, we observed a suppression of electron kinetic energy and dc conductivity under disorder, highlighting the role of Anderson localization in the degradation of electronic transport. These findings offer significant theoretical insight into the stability and critical phenomena of correlated phases in disordered two-dimensional systems.
The interplay between electron-electron correlations and disorder has been a central theme of condensed matter physics over the past several decades, with particular interest in the possibility that interactions might cause delocalization of an Anderson insulator into a metallic state, and the disrupting effects of randomness on magnetic order and the Mott phase. Here we extend this physics to explore electron-phonon interactions and show, via exact quantum Monte Carlo simulations, that the suppression of the charge density wave correlations in the half-filled Holstein model by disorder can stabilize a superconducting phase. Our simulations thus capture qualitatively the suppression of charge ordered phases and emergent superconductivity recently seen experimentally.
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.
The effect of electron-phonon coupling on Dirac fermions was recently explored numerically on a honeycomb lattice, leading to precise quantitative values for the finite temperature and quantum critical points. In this paper, we use the unbiased determinant quantum Monte Carlo method to study the Holstein model on a half-filled staggered-flux square lattice and compare the results with the honeycomb lattice geometry, presenting results for a range of phonon frequencies $0.1\ensuremath{\leqslant}\ensuremath{\omega}\ensuremath{\leqslant}2.0$. We find that the interactions give rise to charge density wave order, but only above a finite coupling strength ${\ensuremath{\lambda}}_{\mathrm{crit}}$. The transition temperature is evaluated and presented in a ${T}_{c}\text{\ensuremath{-}}\ensuremath{\lambda}$ phase diagram. An accompanying mean-field theory calculation also predicts the existence of a quantum phase transition, but at a substantially smaller coupling strength.
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.
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.
The Holstein Hamiltonian describes fermions hopping on a lattice and interacting locally with dispersionless phonon degrees of freedom. In the low-density limit, dressed quasiparticles, polarons and bipolarons, propagate with an effective mass. At higher densities, pairs can condense into a low-temperature superconducting phase and, at or near commensurate filling on a bipartite lattice, to charge density wave (CDW) order. CDW formation breaks a discrete symmetry and hence occurs via a second-order (Ising) transition and therefore at a finite
Tcdw in two dimensions. Quantum Monte Carlo calculations have determined
Tcdw for a variety of geometries, including square, honeycomb, and Lieb lattices. The superconducting transition, on the other hand, in
d=2 is in the Kosterlitz-Thouless universality class and is much less well characterized. In this paper we determine
Tsc for the square lattice for several values of the density $\ensuremath{\rho}$ and phonon frequency ${\ensuremath{\omega}}_{0}$. We find that quasilong-range order sets in at ${T}_{\mathrm{sc}}\ensuremath{\lesssim}t/20$, where
t is the near-neighbor hopping amplitude, consistent with previous rough estimates from simulations which extrapolated to only the temperatures we reach from considerably higher
T. We also show evidence of a discontinuous evolution of the density as the CDW transition is approached at half filling.