Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2601.13121

Charge Order in the half-filled bond-Holstein Model

Charles Jordan, George Issa, Ehsan Khatami, Richard T. Scalettar, Benjamin Cohen-Stead, Steven Johnston

Published Jan 19, 2026• 1 citation
arXivHTMLPDF

Abstract

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 $T_\mathrm{cdw}$ considerably higher than that of the canonical `site-Holstein' model. Using a finite-size scaling analysis of the charge structure factor $S_{\rm cdw}$, we obtain $T_\mathrm{cdw}$ 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\rightarrow 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 $T_\mathrm{cdw}$ but also provide insight into the high-temperature crossover between a metallic and bipolaron liquid regime.

DOI
10.48550/arxiv.2601.13121

Related papers

Competition Between Antiferromagnetic and Charge-Density-Wave Order in the Half-Filled Hubbard-Holstein Model

Elizabeth Nowadnick, Steven Johnston, Brian Moritz +2 more

We present a determinant quantum Monte Carlo study of the competition between instantaneous on-site Coulomb repulsion and retarded phonon-mediated attraction between electrons, as described by the two-dimensional Hubbard-Holstein model. At half filling, we find a strong competition between antiferromagnetism (AFM) and charge-density-wave (CDW) order. We demonstrate that a simple picture of AFM-CDW competition that incorporates the phonon-mediated attraction into an effective-U Hubbard model requires significant refinement. Specifically, retardation effects slow the onset of charge order so that CDW order remains absent even when the effective U is negative. This delay opens a window where neither AFM nor CDW order is well established and where there are signatures of a possible metallic phase.

arXiv:1209.0829•Physical Review Letters•Dec 10, 2012•82 citations

Superconductivity and charge density wave order in the two-dimensional Holstein model

Owen Bradley, G. G. Batrouni, Richard T. Scalettar

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 ${T}_{\mathrm{cdw}}$ in two dimensions. Quantum Monte Carlo calculations have determined ${T}_{\mathrm{cdw}}$ 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 ${T}_{\mathrm{sc}}$ 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.

arXiv:2011.11703•Physical review. B./Physical review. B•Jun 1, 2021•52 citations

Dominant charge density wave correlations in the Holstein model on the half-filled square lattice

Martin Hohenadler, G. G. Batrouni

We use an unbiased, continuous-time quantum Monte Carlo method to address the possibility of a zero-temperature phase without charge-density-wave (CDW) order in the Holstein and, by extension, the Holstein-Hubbard model on the half-filled square lattice. In particular, we present results spanning the whole range of phonon frequencies, allowing us to use the well understood adiabatic and antiadiabatic limits as reference points. For all parameters considered, our data suggest that CDW correlations are stronger than pairing correlations even at very low temperatures. These findings are compatible with a CDW ground state that is also suggested by theoretical arguments.

arXiv:1907.10866•Physical review. B./Physical review. B•Oct 9, 2019•38 citations

Intermediate bond-order-wave phase and nature of the order-to-order transition in the one-dimensional Hubbard-Holstein model

Weber, Manuel

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.

arXiv:2412.13263•Dec 17, 2024•0 citations

The Two-Dimensional Rashba-Holstein Model: A Quantum Monte Carlo Approach

Julián Faúndez, Rodrigo A. Fontenele, Sebastião dos Anjos Sousa-Júnior +2 more

In this work, we investigate the impact of Rashba spin-orbit coupling (RSOC) on the formation of charge-density wave (CDW) and superconducting (SC) phases in the Holstein model on a half-filled square lattice. Using unbiased finite-temperature Quantum Monte Carlo simulations, we go beyond mean-field approaches to determine the ground state order parameter as a function of RSOC and phonon frequency. Our results reveal that the Rashba metal is unstable due to particle-hole instabilities, favoring the emergence of a CDW phase for any RSOC value. In the limit of a pure Rashba hopping, the model exhibits a distinct behavior with the appearance of four Weyl cones at half-filling, where quantum phase transitions are expected to occur at strong interactions. Indeed, a quantum phase transition, belonging to the Gross-Neveu Ising universality class between a semi-metal and CDW emerges at finite phonon frequency dependent coupling $λ_c$. In the antiadiabatic limit we observe an enhanced symmetry in the infrared that unifies SC and CDW orders. These results advance our understanding of competing CDW and SC phases in systems with spin-orbit coupling, providing insights that may help clarify the behavior of related materials.

arXiv:2411.07119•Nov 11, 2024•0 citations

Phonon Dispersion and the Competition between Pairing and Charge Order

Natanael C. Costa, Thomas Blommel, Wei-ting Chiu +2 more

The Holstein model describes the interaction between fermions and a collection of local (dispersionless) phonon modes. In the dilute limit, the phonon degrees of freedom dress the fermions, giving rise to polaron and bipolaron formation. At higher densities, the phonons mediate collective superconducting (SC) and charge-density wave (CDW) phases. Quantum Monte Carlo (QMC) simulations have considered both these limits but have not yet focused on the physics of more general phonon spectra. Here we report QMC studies of the role of phonon dispersion on SC and CDW order in such models. We quantify the effect of finite phonon bandwidth and curvature on the critical temperature T_{cdw} for CDW order and also uncover several novel features of diagonal long-range order in the phase diagram, including a competition between charge patterns at momenta q=(π,π) and q=(0,π) which lends insight into the relationship between Fermi surface nesting and the wave vector at which charge order occurs. We also demonstrate SC order at half filling in situations where a nonzero bandwidth sufficiently suppresses T_{cdw}.

arXiv:1712.04072•Physical Review Letters•May 4, 2018•69 citations
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 $T_\mathrm{cdw}$ considerably higher than that of the canonical site-Holstein' model. Using a finite-size scaling analysis of the charge structure factor ScdwS_{\rm cdw}Scdw​, we obtain TcdwT_\mathrm{cdw}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→0t\rightarrow 0t→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 TcdwT_\mathrm{cdw}Tcdw​ but also provide insight into the high-temperature crossover between a metallic and bipolaron liquid regime.
We present a determinant quantum Monte Carlo study of the competition between instantaneous on-site Coulomb repulsion and retarded phonon-mediated attraction between electrons, as described by the two-dimensional Hubbard-Holstein model. At half filling, we find a strong competition between antiferromagnetism (AFM) and charge-density-wave (CDW) order. We demonstrate that a simple picture of AFM-CDW competition that incorporates the phonon-mediated attraction into an effective-U Hubbard model requires significant refinement. Specifically, retardation effects slow the onset of charge order so that CDW order remains absent even when the effective U is negative. This delay opens a window where neither AFM nor CDW order is well established and where there are signatures of a possible metallic phase.
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{T}_{\mathrm{cdw}}Tcdw​ in two dimensions. Quantum Monte Carlo calculations have determined Tcdw{T}_{\mathrm{cdw}}Tcdw​ for a variety of geometries, including square, honeycomb, and Lieb lattices. The superconducting transition, on the other hand, in d=2d=2d=2 is in the Kosterlitz-Thouless universality class and is much less well characterized. In this paper we determine Tsc{T}_{\mathrm{sc}}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 ttt is the near-neighbor hopping amplitude, consistent with previous rough estimates from simulations which extrapolated to only the temperatures we reach from considerably higher TTT. We also show evidence of a discontinuous evolution of the density as the CDW transition is approached at half filling.
We use an unbiased, continuous-time quantum Monte Carlo method to address the possibility of a zero-temperature phase without charge-density-wave (CDW) order in the Holstein and, by extension, the Holstein-Hubbard model on the half-filled square lattice. In particular, we present results spanning the whole range of phonon frequencies, allowing us to use the well understood adiabatic and antiadiabatic limits as reference points. For all parameters considered, our data suggest that CDW correlations are stronger than pairing correlations even at very low temperatures. These findings are compatible with a CDW ground state that is also suggested by theoretical arguments.
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.
In this work, we investigate the impact of Rashba spin-orbit coupling (RSOC) on the formation of charge-density wave (CDW) and superconducting (SC) phases in the Holstein model on a half-filled square lattice. Using unbiased finite-temperature Quantum Monte Carlo simulations, we go beyond mean-field approaches to determine the ground state order parameter as a function of RSOC and phonon frequency. Our results reveal that the Rashba metal is unstable due to particle-hole instabilities, favoring the emergence of a CDW phase for any RSOC value. In the limit of a pure Rashba hopping, the model exhibits a distinct behavior with the appearance of four Weyl cones at half-filling, where quantum phase transitions are expected to occur at strong interactions. Indeed, a quantum phase transition, belonging to the Gross-Neveu Ising universality class between a semi-metal and CDW emerges at finite phonon frequency dependent coupling λcλ_cλc​. In the antiadiabatic limit we observe an enhanced symmetry in the infrared that unifies SC and CDW orders. These results advance our understanding of competing CDW and SC phases in systems with spin-orbit coupling, providing insights that may help clarify the behavior of related materials.
The Holstein model describes the interaction between fermions and a collection of local (dispersionless) phonon modes. In the dilute limit, the phonon degrees of freedom dress the fermions, giving rise to polaron and bipolaron formation. At higher densities, the phonons mediate collective superconducting (SC) and charge-density wave (CDW) phases. Quantum Monte Carlo (QMC) simulations have considered both these limits but have not yet focused on the physics of more general phonon spectra. Here we report QMC studies of the role of phonon dispersion on SC and CDW order in such models. We quantify the effect of finite phonon bandwidth and curvature on the critical temperature T_{cdw} for CDW order and also uncover several novel features of diagonal long-range order in the phase diagram, including a competition between charge patterns at momenta q=(π,π) and q=(0,π) which lends insight into the relationship between Fermi surface nesting and the wave vector at which charge order occurs. We also demonstrate SC order at half filling in situations where a nonzero bandwidth sufficiently suppresses T_{cdw}.