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.
Phase diagram of the two-dimensional Hubbard-Holstein model: enhancement of
s-wave pairing between charge and magnetic orders
We investigate the role of electron-electron and electron-phonon interactions in strongly correlated systems by performing unbiased quantum Monte Carlo simulations in the square lattice Hubbard-Holstein model at half-filling. We study the competition and interplay between antiferromagnetism (AFM) and charge-density wave (CDW), establishing its very rich phase diagram. In the region between AFM and CDW phases, we have found an enhancement of superconducting pairing correlations, favoring (nonlocal)
s-wave pairs. Our study sheds light over past inconsistencies in the literature, in particular the emergence of CDW in the pure Holstein model case.
By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping
t, the on-site electron-electron interaction
U, the phonon energy ${\ensuremath{\omega}}_{0}$, and the electron-phonon coupling
g. At half filling, the ground state is an antiferromagnetic insulator for $U\ensuremath{\gtrsim}2{g}^{2}/{\ensuremath{\omega}}_{0}$, while it is a charge-density-wave (or bipolaronic) insulator for $U\ensuremath{\lesssim}2{g}^{2}/{\ensuremath{\omega}}_{0}$. In addition to these phases, we find a superconducting phase that intrudes between them. For ${\ensuremath{\omega}}_{0}/t=1$, superconductivity emerges when both
U/t and $2{g}^{2}/t{\ensuremath{\omega}}_{0}$ are small; then, by increasing the value of the phonon energy ${\ensuremath{\omega}}_{0}$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case.
Abstract The electron–electron and electron–phonon interactions play an important role in correlated materials, being key features for spin, charge and pair correlations. Thus, here we investigate their effects in strongly correlated systems by performing unbiased quantum Monte Carlo simulations in the square lattice Hubbard-Holstein model at half-filling. We study the competition and interplay between antiferromagnetism (AFM) and charge-density wave (CDW), establishing its very rich phase diagram. In the region between AFM and CDW phases, we have found an enhancement of superconducting pairing correlations, favouring (nonlocal) s-wave pairs. Our study sheds light over past inconsistencies in the literature, in particular the emergence of CDW in the pure Holstein model case.
The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.
The Hubbard-Holstein model is one of the central models that describe the competition between electron-electron and electron-phonon interactions. In one dimension and at half-filling, the interplay between an electronic spin-density wave and a phonon-driven charge-density wave is considered to stabilize an intermediate Luther-Emery liquid. Here we show that, once the extended metallic regime disappears, a narrow bond-order-wave phase emerges. We use an exact directed-loop quantum Monte Carlo method for retarded interactions to simulate system sizes of several hundred sites, necessary to observe the weak signatures of this novel phase. Our results suggest a second-order quantum phase transition between the two charge orders that only turns first-order at strong coupling. Our findings are reminiscent of the extended Hubbard model, but they are driven by competing frequency dependencies of the same dynamically-screened Hubbard interaction.