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

Physical review. B, Condensed matter

Peierls instability in the two-dimensional half-filled Hubbard model

Sanyee Tang, J. E. Hirsch

Published Jun 1, 1988• 85 citations
Publisher

Abstract

We have studied the tendency of the two-dimensional half-filled Hubbard model towards dimerization as a function of the on-site repulsion U. We have found that the phonon with (\ensuremath{\pi},\ensuremath{\pi}) wave vector is the most favorable vibration mode when U is less than a critical value ${U}_{c}$, where ${U}_{c}$ is a function of the dimerization. When U is greater than ${U}_{c}$, the (\ensuremath{\pi},0) phonon becomes preferred. In leading order, the energy gain under dimerization appears to be proportional to ${\ensuremath{\delta}}^{2}$ln\ensuremath{\delta} if the (\ensuremath{\pi},\ensuremath{\pi}) phonon is favorable, so that the system displays a Peierls instability in this case. For large U the energy gain to leading order appears to be proportional to ${\ensuremath{\delta}}^{2}$ for both the (\ensuremath{\pi},\ensuremath{\pi}) and (\ensuremath{\pi},0) phonons, and the latter one is preferred; in that case the system is stable unless the electron-phonon coupling exceeds a critical value. Unlike the effect in the one-dimensional case, the on-site repulsion U seems to have essentially no effect on the tendency towards dimerization for small U, and suppresses it for large U. This difference may be related to the fact that there exists long-range antiferromagnetic order in the two-dimensional half-filled Hubbard model when U is present.

DOI
10.1103/physrevb.37.9546

Related papers

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

F. Hébert, Bo Xiao, V. G. Rousseau +2 more

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.

arXiv:1811.07811•Physical review. B./Physical review. B•Feb 6, 2019•17 citations

Finite-size studies of phases and dimerization in one-dimensional extended Peierls-Hubbard models

V. Waas, H. Büttner, Johannes Voit

We study extended Hubbard and Peierls-Hubbard models with up to N=10 sites and half-filled bands with modified boundary conditions. For the extended Hubbard model, the finite-size dependence of the transition from the charge-density-wave regime (CDW) to the spin-density-wave (SDW) regime and of the condensation transition is critically examined, clarifying some previous results. For the Peierls-Hubbard model and N\ensuremath{\le}10 with modified periodic boundary conditions, we always see a decrease of dimerization with U. Extrapolating these results to infinite N, we find that dimerization is enhanced by U for small electron-phonon coupling \ensuremath{\alpha}\ifmmode \tilde{}\else \~{}\fi{}0.75, whereas it decreases with U for 0.75\ensuremath{\le}1. By extending this model by a nearest-neighbor repulsion V, the amplitude of dimerization reaches a maximum at the same value U\ensuremath{\le}2V, where the CDW-to-SDW transition occurs in the extended Hubbard model.

Physical review. B, Condensed matter•May 1, 1990•37 citations

Coexistence of the bond-order wave and antiferromagnetism in a two-dimensional half-filled Peierls-Hubbard model

Qingshan Yuan, Thilo Kopp

The two-dimensional Peierls-Hubbard model is studied at half filling within both Hartree-Fock and Kotliar-Ruckenstein slave-boson theory. The interplay between two types of long-range order, bond-order wave (BOW) and antiferromagnetism (AFM), is analyzed for two representative dimerization patterns, corresponding both to the same wave vector $(\ensuremath{\pi},\ensuremath{\pi}).$ For each pattern, the Peierls dimerization (and associated BOW) is weakened and finally suppressed with increasing Hubbard on-site interaction, and correspondingly AFM is gradually enhanced. In particular, a coexistence regime with both BOW and AFM order is established in the parameter space of electron-lattice and Hubbard interactions.

arXiv:cond-mat/0109435•Physical review. B, Condensed matter•Feb 4, 2002•17 citations

Computational Studies of Correlated Electrons in One-Dimension

Pinaki Sengupta

"We have studied the extended Hubbard and the Peierls-extended Hubbard models in one dimension using primarily computational methods. We discuss the derivation of the models and the computational methods used to study them, focusing primarily on our extensions over existing numerical methods. In particular, we present the construction of ""operator-loop"" updates for spin-1/2 fermions and the inclusion of electron-phonon interactions in the Stochastic Series Expansion Quantum Monte Carlo method and the use of Quantum Parallel Tempering to improve the quality of the numerical data. For the extended Hubbard model at half-filling, we have confirmed the existence of a recently reported bond-order wave phase for the ground state in a small region of parameter space between the well-known charge- and spin-density wave phases. We have studied in detail the nature of the associated quantum phase transitions between all the phases. Upon doping away from half-filling, the doubly-degenerate, spin and charge gapped long-range-ordered ground states of the model are shown to evolve to a Luther-Emery liquid with no charge gap but a finite spin gap. For the Peierls-extended Hubbard model, which includes both electron-phonon and electron-electron interactions, we have studied the dimerization transition of the half-filled system in the presence of fully quantum mechanical phonons. We show that the transition to a dimerized state occurs only when the electron-phonon coupling exceeds a finite critical value that depends on the bare phonon frequency. This is in contrast to previous studies in similar models with adiabatic phonons that had predicted a dimerized state for arbitrarily small values of the electron-phonon coupling."

•Jan 1, 2001•0 citations

Effect of Quantum Fluctuations on the Peierls Instability: A Monte Carlo Study

J. E. Hirsch, Eduardo Fradkin

The effect of quantum fluctuations of the phonon field on the ground-state character of one-dimensional electron-phonon systems is studied with a numerical simulation technique. Two different electron-phonon models are studied, both for spinless and spin-\textonehalf{} electrons, in the half-filled-band case. In both models, it is found that in the spinless case quantum fluctuations destroy the long-range dimerization order for small electron-phonon coupling constant $\ensuremath{\lambda}$. For the spin-\textonehalf{} case, the results are consistent with the Peierls transition occurring at $\ensuremath{\lambda}=0$.

Physical Review Letters•Aug 9, 1982•115 citations

Incremental expansions for Hubbard-Peierls systems

J. Malek, K. Kladko, Sergej Flach A

The ground state energies of infinite half-filled Hubbard-Peierls chains are investigated by combining an incremental expansion with exact diagonalization of finite chain segments. The ground-state energy of equidistant infinite Hubbard (Heisenberg) chains is calculated with a relative error of less than 3×10−3 for all values of U using diagonalations of 12-site (20-site) chain segments. For dimerized chains the dimerization order parameter d as a function of the on-site repulsion interaction U has a maximum at nonzero values of U, if the electron-phonon coupling g is lower than a critical value g c . The critical value g c is found with high accuracy to be g c =0.69. For smaller values of g the position of the maximum of d(U) is approximately 3t, and rapidly tends to zero as g approaches g c from below. We show how our method can be applied to calculate breathers for the problem of phonon dynamics in Hubbard-Peierls systems.

arXiv:cond-mat/9710112•Journal of Experimental and Theoretical Physics Letters•Jun 1, 1998•5 citations
We have studied the tendency of the two-dimensional half-filled Hubbard model towards dimerization as a function of the on-site repulsion U. We have found that the phonon with (\ensuremath{\pi},\ensuremath{\pi}) wave vector is the most favorable vibration mode when U is less than a critical value Uc{U}_{c}Uc​, where Uc{U}_{c}Uc​ is a function of the dimerization. When U is greater than Uc{U}_{c}Uc​, the (\ensuremath{\pi},0) phonon becomes preferred. In leading order, the energy gain under dimerization appears to be proportional to \ensuremathδ2{\ensuremath{\delta}}^{2}\ensuremathδ2ln\ensuremath{\delta} if the (\ensuremath{\pi},\ensuremath{\pi}) phonon is favorable, so that the system displays a Peierls instability in this case. For large U the energy gain to leading order appears to be proportional to \ensuremathδ2{\ensuremath{\delta}}^{2}\ensuremathδ2 for both the (\ensuremath{\pi},\ensuremath{\pi}) and (\ensuremath{\pi},0) phonons, and the latter one is preferred; in that case the system is stable unless the electron-phonon coupling exceeds a critical value. Unlike the effect in the one-dimensional case, the on-site repulsion U seems to have essentially no effect on the tendency towards dimerization for small U, and suppresses it for large U. This difference may be related to the fact that there exists long-range antiferromagnetic order in the two-dimensional half-filled Hubbard model when U is present.
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{\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 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{\ne}0$ case, even for small values of the coupling range.
We study extended Hubbard and Peierls-Hubbard models with up to N=10 sites and half-filled bands with modified boundary conditions. For the extended Hubbard model, the finite-size dependence of the transition from the charge-density-wave regime (CDW) to the spin-density-wave (SDW) regime and of the condensation transition is critically examined, clarifying some previous results. For the Peierls-Hubbard model and N\ensuremath{\le}10 with modified periodic boundary conditions, we always see a decrease of dimerization with U. Extrapolating these results to infinite N, we find that dimerization is enhanced by U for small electron-phonon coupling \ensuremath{\alpha}\ifmmode \tilde{}\else \~{}\fi{}0.75, whereas it decreases with U for 0.75\ensuremath{\le}1. By extending this model by a nearest-neighbor repulsion V, the amplitude of dimerization reaches a maximum at the same value U\ensuremath{\le}2V, where the CDW-to-SDW transition occurs in the extended Hubbard model.
The two-dimensional Peierls-Hubbard model is studied at half filling within both Hartree-Fock and Kotliar-Ruckenstein slave-boson theory. The interplay between two types of long-range order, bond-order wave (BOW) and antiferromagnetism (AFM), is analyzed for two representative dimerization patterns, corresponding both to the same wave vector $(\ensuremath{\pi},\ensuremath{\pi}).$ For each pattern, the Peierls dimerization (and associated BOW) is weakened and finally suppressed with increasing Hubbard on-site interaction, and correspondingly AFM is gradually enhanced. In particular, a coexistence regime with both BOW and AFM order is established in the parameter space of electron-lattice and Hubbard interactions.
"We have studied the extended Hubbard and the Peierls-extended Hubbard models in one dimension using primarily computational methods. We discuss the derivation of the models and the computational methods used to study them, focusing primarily on our extensions over existing numerical methods. In particular, we present the construction of ""operator-loop"" updates for spin-1/2 fermions and the inclusion of electron-phonon interactions in the Stochastic Series Expansion Quantum Monte Carlo method and the use of Quantum Parallel Tempering to improve the quality of the numerical data. For the extended Hubbard model at half-filling, we have confirmed the existence of a recently reported bond-order wave phase for the ground state in a small region of parameter space between the well-known charge- and spin-density wave phases. We have studied in detail the nature of the associated quantum phase transitions between all the phases. Upon doping away from half-filling, the doubly-degenerate, spin and charge gapped long-range-ordered ground states of the model are shown to evolve to a Luther-Emery liquid with no charge gap but a finite spin gap. For the Peierls-extended Hubbard model, which includes both electron-phonon and electron-electron interactions, we have studied the dimerization transition of the half-filled system in the presence of fully quantum mechanical phonons. We show that the transition to a dimerized state occurs only when the electron-phonon coupling exceeds a finite critical value that depends on the bare phonon frequency. This is in contrast to previous studies in similar models with adiabatic phonons that had predicted a dimerized state for arbitrarily small values of the electron-phonon coupling."
The effect of quantum fluctuations of the phonon field on the ground-state character of one-dimensional electron-phonon systems is studied with a numerical simulation technique. Two different electron-phonon models are studied, both for spinless and spin-\textonehalf{} electrons, in the half-filled-band case. In both models, it is found that in the spinless case quantum fluctuations destroy the long-range dimerization order for small electron-phonon coupling constant $\ensuremath{\lambda}$. For the spin-\textonehalf{} case, the results are consistent with the Peierls transition occurring at $\ensuremath{\lambda}=0$.
The ground state energies of infinite half-filled Hubbard-Peierls chains are investigated by combining an incremental expansion with exact diagonalization of finite chain segments. The ground-state energy of equidistant infinite Hubbard (Heisenberg) chains is calculated with a relative error of less than 3×10−3 for all values of U using diagonalations of 12-site (20-site) chain segments. For dimerized chains the dimerization order parameter d as a function of the on-site repulsion interaction U has a maximum at nonzero values of U, if the electron-phonon coupling g is lower than a critical value g c . The critical value g c is found with high accuracy to be g c =0.69. For smaller values of g the position of the maximum of d(U) is approximately 3t, and rapidly tends to zero as g approaches g c from below. We show how our method can be applied to calculate breathers for the problem of phonon dynamics in Hubbard-Peierls systems.