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:1811.07811 · Physical review. B./Physical review. B

One-dimensional Hubbard-Holstein model with finite-range electron-phonon coupling

F. Hébert, Bo Xiao, V. G. Rousseau, Richard T. Scalettar, G. G. Batrouni

Published Feb 6, 2019• 17 citations
arXivHTMLPDF

Abstract

The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling $g$ to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction $U$ drives the fermions into a Mott antiferromagnet. At low $g$ and $U$, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range $\ensuremath{\xi}$, interpolating between the Holstein and Fr\"ohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range $\ensuremath{\xi}$ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction $U$ will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the $U\ensuremath{\ne}0$ case, even for small values of the coupling range.

DOI
10.1103/physrevb.99.075108

Related papers

Superconductivity, charge-density waves, antiferromagnetism, and phase separation in the Hubbard-Holstein model

Seher Karakuzu, Luca F. Tocchio, Sandro Sorella +1 more

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.

arXiv:1709.00278•Physical review. B./Physical review. B•Nov 27, 2017•54 citations

Insulating and metallic phases in the one-dimensional Hubbard-Su-Schrieffer-Heeger model: Insights from a backflow-inspired variational wave function

Davide Piccioni, Francesco Ferrari, Michele Fabrizio +1 more

The interplay between electron-electron and electron-phonon interactions is studied in a one-dimensional lattice model by means of a variational Monte Carlo method based on generalized Jastrow-Slater wave functions. Here, the fermionic part is constructed by a pair-product state, which explicitly depends on the phonon configuration, thus including the electron-phonon coupling in a backflow-inspired way. We report the results for the Hubbard model in the presence of the Su-Schrieffer-Heeger coupling to optical phonons, both at half filling and upon hole doping. At half filling, the ground state is either a translationally invariant Mott insulator, with gapless spin excitations, or a Peierls insulator, which breaks translations and has fully gapped excitations. Away from half filling, the charge gap closes in both Mott and Peierls insulators, turning the former into a conventional Luttinger liquid (gapless in all excitation channels). In the latter case, instead, a finite spin gap remains at small doping. Even though consistent with the general theory of interacting electrons in one dimension, the existence of such a phase (with gapless charge but gapped spin excitations) has never been demonstrated in a model with repulsive interaction and with only two Fermi points. Since the spin-gapped metal represents the one-dimensional counterpart of a superconductor, our results furnish evidence that a true off-diagonal long-range order may exist in the two-dimensional case.

arXiv:2407.03046•Physical review. B./Physical review. B•Jan 10, 2025•0 citations

Nature of ground states in one-dimensional electron-phonon Hubbard models at half filling

H. Bakrim, C. Bourbonnais

The renormalization group technique is applied to one-dimensional electron-phonon Hubbard models at half filling and zero temperature. For the Holstein-Hubbard model, the results of one-loop calculations are congruent with the phase diagram obtained by quantum Monte Carlo simulations in the $(U,{g}_{\mathrm{ph}})$ plane for the phonon-mediated interaction ${g}_{\mathrm{ph}}$ and the Coulomb interaction $U$. The incursion of an intermediate phase between a fully gapped charge-density-wave state and a Mott antiferromagnet is supported along with the growth of its size with the molecular phonon frequency ${\ensuremath{\omega}}_{0}$. We find additional phases enfolding the base boundary of the intermediate phase. A Luttinger liquid line is found below some critical ${U}^{*}\ensuremath{\approx}{g}_{\mathrm{ph}}^{*}$, followed at larger $U\ensuremath{\sim}{g}_{\mathrm{ph}}$ by a narrow region of bond-order-wave ordering which is either charge or spin gapped depending on $U$. For the Peierls-Hubbard model, the region of the $(U,{g}_{\mathrm{ph}})$ plane with a fully gapped Peierls-bond-order-wave state shows a growing domination over the Mott gapped antiferromagnet as the Debye frequency ${\ensuremath{\omega}}_{D}$ decreases. A power-law dependence ${g}_{\mathrm{ph}}\ensuremath{\sim}{U}^{2\ensuremath{\eta}}$ is found to map out the boundary between the two phases, whose exponent is in good agreement with the existing quantum Monte Carlo simulations performed when a finite nearest-neighbor repulsion term $V$ is added to the Hubbard interaction.

arXiv:1502.01981•Physical Review B•Feb 19, 2015•25 citations

Finite-size effects in Luther-Emery phases of Holstein and Hubbard models

Jonas Greitemann, Stephan Heßelmann, Stefan Weßel +2 more

The one-dimensional Holstein model and its generalizations have been studied extensively to understand the effects of electron-phonon interaction. The half-filled case is of particular interest, as it describes a transition from a metallic phase with a spin gap due to attractive backscattering to a Peierls insulator with charge-density-wave order. Our quantum Monte Carlo results support the existence of a metallic phase with dominant power-law charge correlations, as described by the Luther-Emery fixed point. We demonstrate that for Holstein and also for purely fermionic models the spin gap significantly complicates finite-size numerical studies, and explains inconsistent previous results for Luttinger parameters and phase boundaries. On the other hand, no such complications arise in spinless models. The correct low-energy theory of the spinful Holstein model is argued to be that of singlet bipolarons with a repulsive, mutual interaction. This picture naturally explains the existence of a metallic phase, but also implies that gapless Luttinger liquid theory is not applicable.

arXiv:1510.07454•Physical Review B•Dec 28, 2015•25 citations

Interplay of Site and Bond Electron-Phonon Coupling in One Dimension

Martin Hohenadler

The interplay of bond and charge correlations is studied in a one-dimensional model with both Holstein and Su-Schrieffer-Heeger (SSH) couplings to quantum phonons. The problem is solved exactly by quantum Monte Carlo simulations. If one of the couplings dominates, the ground state is a Peierls insulator with long-range bond or charge order. At weak coupling, the results suggest a spin-gapped and repulsive metallic phase arising from the competing order parameters and lattice fluctuations. Such a phase is absent from the pure SSH model even for quantum phonons. At strong coupling, evidence for a continuous transition between the two Peierls states is presented.

arXiv:1610.02426•Physical Review Letters•Nov 11, 2016•18 citations

Phonon spectral function of the one-dimensional Holstein-Hubbard model

Manuel Weber, Fakher F. Assaad, Martin Hohenadler

We use the continuous-time interaction expansion (CT-INT) quantum Monte Carlo method to calculate the phonon spectral function of the one-dimensional Holstein-Hubbard model at half-filling. Our results are consistent with a soft-mode Peierls transition in the adiabatic regime, and the existence of a central peak related to long-range order in the Peierls phase. We explain a previously observed feature at small momenta in terms of a hybridization of charge and phonon excitations. Tuning the system from a Peierls to a metallic phase with a nonzero Hubbard interaction suppresses the central peak, but a significant renormalization of the phonon dispersion remains. In contrast, the dispersion is only weakly modified in the Mott phase. We discuss finite-size effects, the relation to the dynamic charge structure factor, as well as additional sum rules and their implications. Finally, we reveal the existence of a discrete symmetry in a continuum field theory of the Holstein model, which is spontaneously broken in the Peierls phase.

arXiv:1504.07446•Physical Review B•Jun 29, 2015•22 citations
The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling ggg to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction UUU drives the fermions into a Mott antiferromagnet. At low ggg and UUU, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range \ensuremathξ\ensuremath{\xi}\ensuremathξ, interpolating between the Holstein and Fr"ohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range \ensuremathξ\ensuremath{\xi}\ensuremathξ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction UUU will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the U\ensuremath≠0U\ensuremath{\ne}0U\ensuremath=0 case, even for small values of the coupling range.
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 ttt, the on-site electron-electron interaction UUU, the phonon energy ${\ensuremath{\omega}}_{0}$, and the electron-phonon coupling ggg. 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/tU/tU/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.
The interplay between electron-electron and electron-phonon interactions is studied in a one-dimensional lattice model by means of a variational Monte Carlo method based on generalized Jastrow-Slater wave functions. Here, the fermionic part is constructed by a pair-product state, which explicitly depends on the phonon configuration, thus including the electron-phonon coupling in a backflow-inspired way. We report the results for the Hubbard model in the presence of the Su-Schrieffer-Heeger coupling to optical phonons, both at half filling and upon hole doping. At half filling, the ground state is either a translationally invariant Mott insulator, with gapless spin excitations, or a Peierls insulator, which breaks translations and has fully gapped excitations. Away from half filling, the charge gap closes in both Mott and Peierls insulators, turning the former into a conventional Luttinger liquid (gapless in all excitation channels). In the latter case, instead, a finite spin gap remains at small doping. Even though consistent with the general theory of interacting electrons in one dimension, the existence of such a phase (with gapless charge but gapped spin excitations) has never been demonstrated in a model with repulsive interaction and with only two Fermi points. Since the spin-gapped metal represents the one-dimensional counterpart of a superconductor, our results furnish evidence that a true off-diagonal long-range order may exist in the two-dimensional case.
The renormalization group technique is applied to one-dimensional electron-phonon Hubbard models at half filling and zero temperature. For the Holstein-Hubbard model, the results of one-loop calculations are congruent with the phase diagram obtained by quantum Monte Carlo simulations in the (U,gph)(U,{g}_{\mathrm{ph}})(U,gph​) plane for the phonon-mediated interaction gph{g}_{\mathrm{ph}}gph​ and the Coulomb interaction UUU. The incursion of an intermediate phase between a fully gapped charge-density-wave state and a Mott antiferromagnet is supported along with the growth of its size with the molecular phonon frequency ${\ensuremath{\omega}}_{0}$. We find additional phases enfolding the base boundary of the intermediate phase. A Luttinger liquid line is found below some critical ${U}^{*}\ensuremath{\approx}{g}_{\mathrm{ph}}^{*}$, followed at larger $U\ensuremath{\sim}{g}_{\mathrm{ph}}$ by a narrow region of bond-order-wave ordering which is either charge or spin gapped depending on UUU. For the Peierls-Hubbard model, the region of the (U,gph)(U,{g}_{\mathrm{ph}})(U,gph​) plane with a fully gapped Peierls-bond-order-wave state shows a growing domination over the Mott gapped antiferromagnet as the Debye frequency ${\ensuremath{\omega}}_{D}$ decreases. A power-law dependence ${g}_{\mathrm{ph}}\ensuremath{\sim}{U}^{2\ensuremath{\eta}}$ is found to map out the boundary between the two phases, whose exponent is in good agreement with the existing quantum Monte Carlo simulations performed when a finite nearest-neighbor repulsion term VVV is added to the Hubbard interaction.
The one-dimensional Holstein model and its generalizations have been studied extensively to understand the effects of electron-phonon interaction. The half-filled case is of particular interest, as it describes a transition from a metallic phase with a spin gap due to attractive backscattering to a Peierls insulator with charge-density-wave order. Our quantum Monte Carlo results support the existence of a metallic phase with dominant power-law charge correlations, as described by the Luther-Emery fixed point. We demonstrate that for Holstein and also for purely fermionic models the spin gap significantly complicates finite-size numerical studies, and explains inconsistent previous results for Luttinger parameters and phase boundaries. On the other hand, no such complications arise in spinless models. The correct low-energy theory of the spinful Holstein model is argued to be that of singlet bipolarons with a repulsive, mutual interaction. This picture naturally explains the existence of a metallic phase, but also implies that gapless Luttinger liquid theory is not applicable.
The interplay of bond and charge correlations is studied in a one-dimensional model with both Holstein and Su-Schrieffer-Heeger (SSH) couplings to quantum phonons. The problem is solved exactly by quantum Monte Carlo simulations. If one of the couplings dominates, the ground state is a Peierls insulator with long-range bond or charge order. At weak coupling, the results suggest a spin-gapped and repulsive metallic phase arising from the competing order parameters and lattice fluctuations. Such a phase is absent from the pure SSH model even for quantum phonons. At strong coupling, evidence for a continuous transition between the two Peierls states is presented.
We use the continuous-time interaction expansion (CT-INT) quantum Monte Carlo method to calculate the phonon spectral function of the one-dimensional Holstein-Hubbard model at half-filling. Our results are consistent with a soft-mode Peierls transition in the adiabatic regime, and the existence of a central peak related to long-range order in the Peierls phase. We explain a previously observed feature at small momenta in terms of a hybridization of charge and phonon excitations. Tuning the system from a Peierls to a metallic phase with a nonzero Hubbard interaction suppresses the central peak, but a significant renormalization of the phonon dispersion remains. In contrast, the dispersion is only weakly modified in the Mott phase. We discuss finite-size effects, the relation to the dynamic charge structure factor, as well as additional sum rules and their implications. Finally, we reveal the existence of a discrete symmetry in a continuum field theory of the Holstein model, which is spontaneously broken in the Peierls phase.