We study the half-filled Haldane-Holstein model, where a paradigmatic Chern insulator is coupled to fully dynamical phonons, and provide an unbiased characterization of how retarded electron-phonon interactions destabilize Chern topology. Using determinant quantum Monte Carlo, we find that increasing the coupling drives an abrupt, first-order transition from a Chern insulator to a staggered charge-density wave that acts as a dynamical sublattice (Semenoff) mass. The transition is simultaneously signaled by a nearly quantized many-body Bott index and a real-space local Chern marker constructed from the interacting Green's function, both of which collapse as the charge order parameter becomes extensive. Spectral and open-boundary calculations reveal concomitant gap closing and the loss of boundary spectral weight at the critical coupling. Despite the generic phase problem induced by broken time-reversal symmetry, we show that it remains mild in the low-frequency regime and that the average phase factor sharply tracks the Chern insulator-charge density wave boundary. Our results establish a concrete route by which electron-phonon coupling can trigger a discontinuous collapse of Chern topology and provide experimentally relevant signatures for correlated topological platforms.
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.
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 quantum phase transition between a repulsive Luttinger liquid and an insulating Peierls state is studied in the framework of the one-dimensional spinless Holstein model. We focus on the adiabatic regime but include the full quantum dynamics of the phonons. Using continuous-time quantum Monte Carlo simulations, we track, in particular, the dynamic charge structure factor and the single-particle spectrum across the transition. With increasing electron-phonon coupling, the dynamic charge structure factor reveals the emergence of a charge gap and a clear signature of phonon softening at the zone boundary. The single-particle spectral function evolves continuously across the transition. Hybridization of the charge and phonon modes of the Luttinger liquid description leads to two modes, one of which corresponds to the coherent polaron band. This band acquires a gap upon entering the Peierls phase, whereas the other mode constitutes the incoherent, high-energy spectrum with backfolded shadow bands. Coherent polaronic motion is a direct consequence of quantum lattice fluctuations. In the strong-coupling regime, the spectrum is described by the static, mean-field limit. Importantly, whereas finite electron density in general leads to screening of polaron effects, the latter reappear at half-filling due to charge ordering and lattice dimerization.
Existing quantum Monte Carlo studies have investigated the properties of fermions on a Lieb (
CuO2) lattice interacting with on-site or near-neighbor electron-electron coupling. Attention has focused on the interplay of such interactions with the macroscopic degeneracy of local zero-energy modes, from which Bloch states can be formed to produce a flat band in which energy is independent of momentum. The resulting high density of states, in combination with the Stoner criterion, suggests that there should be pronounced instabilities to ordered phases. Indeed, a theorem by Lieb rigorously establishes the existence of ferrimagnetic order. Here we study the charge density wave phases induced by electron-phonon coupling on the Lieb lattice, as opposed to previous work on electron-electron interactions. Our key result is the demonstration of charge density wave phases at one third and two thirds fillings, characterized by long-range density-density correlations between doubly occupied sites on the minority or majority sublattices, and an accompanying gap. We also compute the transition temperature to the ordered phase as a function of the electron-phonon coupling.
We use determinant quantum Monte Carlo to study the half-filled `bond-Holstein' model on a square lattice. We find that the model exhibits a charge-density-wave (CDW) phase transition with a critical temperature
Tcdw considerably higher than that of the canonical `site-Holstein' model. Using a finite-size scaling analysis of the charge structure factor
Scdw, we obtain
Tcdw to greater than one percent accuracy. At the same time, local observables also show clear signatures consistent with the transition temperatures inferred from our scaling analysis. We attribute the enhanced CDW tendencies to a phonon-mediated nearest-neighbor electron repulsion that is directly proportional to the dimensionless electron-phonon coupling
λ in the atomic (
t→0) limit. This behavior contrasts with the site-Holstein case, where the same limit yields only an on-site attraction. We supplement our analysis with results from several unsupervised machine learning methods, which not only confirm our estimates of
Tcdw but also provide insight into the high-temperature crossover between a metallic and bipolaron liquid regime.
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.