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

Small Bipolarons in the 2-dimensional Holstein-Hubbard Model I The Adiabatic Limit

Laurent Proville, S. Aubry

Published Jan 1, 1998• 27 citations
OpenAlex

Abstract

Abstract. The spatially localized bound states of two electrons in the adiabatic two-dimensional Holstein-Hubbard model on a square lattice are investigated both numerically and analytically. The interplay between the electron-phonon coupling g, which tends to form bipolarons and the repulsive Hubbard interaction υ ≥ 0, which tends to break them, generates many different ground-states. There are four domains in the g,υ phase diagram delimited by first order transition lines. Except for the domain at weak electronphonon coupling (small g) where the electrons remain free, the electrons form bipolarons which can 1) be mostly located on a single site (small υ, large g); 2) be an anisotropic pair of polarons lying on two neighboring sites in the magnetic singlet state (large υ, large g); or 3) be a ”quadrisinglet state ” which is the superposition of 4 electronic singlets with a common central site. This quadrisinglet bipolaron is the most stable in a small central domain in between the three other phases. The pinning modes and the Peierls-Nabarro barrier of each of these bipolarons are calculated and the barrier is found to be strongly depressed in the region of stability of the quadrisinglet bipolaron. PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) – 71.38.+i Polarons and electron phonon interactions – 74.20.Mn Nonconventional mechanisms (spin fluctuations, polarons and bipolarons, resonating valence bond model, anyon mechanism, marginal Fermi liquid, Luttinger liquid, etc.) – 74.25.Jb Electronic structure

Related papers

Mobile Bipolarons in the Adiabatic Holstein-Hubbard Model in 1 and 2 dimensions

Laurent Proville, S. Aubry

The bound states of two electrons in the adiabatic Holstein-Hubbard model are studied numerically in one and two dimensions from the anticontinuous limit. This model involves a competition between a local electron-phonon coupling (with a classical lattice) which tends to form pairs of electrons and the repulsive Hubbard interaction $U \geq 0$ which tends to break them. In 1D, the ground-state always consists in a pair of localized polarons in a singlet state. They are located at the same site for U=0. Increasing U, there is a first order transition at which the bipolaron becomes a spin singlet pair of two polarons bounded by a magnetic interaction. The pinning mode of the bipolaron soften in the vicinity of this transition leading to a higher mobility of the bipolaron which is tested numerically. In 2D, and for any $U$, the electron-phonon coupling needs to be large enough in order to form small polarons or bipolarons instead of extended electrons. We calculate the phase diagram of the bipolaron involving first order transitions lines with a triple point. A pair of polarons can form three types of bipolarons: a) on a single site at small $U$, b) a spin singlet state on two nearest neighbor sites for larger $U$ as in 1D and c) a new intermediate state obtained as the resonant combination of four 2-sites singlet states sharing a central site, called quadrisinglet. The breathing and pinning internal modes of bipolarons in 2D generally only weakly soften and thus, they are practically not mobile. On the opposite, in the vicinity of the triple point involving the quadrisinglet, both modes exhibit a significant softening. However, it was not sufficient for allowing the existence of a classical mobile bipolaron (at least in that model).

•Jan 1, 1998•0 citations

Two-dimensional Hubbard-Holstein bipolaron

Alexandru Macridin, G. A. Sawatzky, Mark Jarrell

We present a diagrammatic Monte Carlo study of the properties of the Hubbard-Holstein bipolaron on a two-dimensional square lattice. With a small Coulomb repulsion $U$ and with increasing electron-phonon interaction, and when reaching a value about two times smaller than the one corresponding to the transition of light polaron to heavy polaron, the system suffers a sharp transition from a state formed by two weakly bound light polarons to a heavy, strongly bound on-site bipolaron. Aside from this rather conventional bipolaron a new bipolaron state is found for large $U$ at intermediate and large electron-phonon coupling, corresponding to two polarons bound on nearest-neighbor sites. We discuss both the properties of the different bipolaron states and the transition from one state to another. We present a phase diagram in parameter space defined by the electron-phonon coupling and $U$. Our numerical method does not use any artificial approximation and can be easily modified to other bipolaron models with longer range electron-phonon and∕or electron-electron interaction.

arXiv:cond-mat/0401305•Physical Review B•Jun 23, 2004•110 citations

Formation of bipolarons in the Holstein Hubbard model in one, two and three dimensions

Shreekantha Sil

In this paper we have studied the properties of a two-fermion system interacting with the phonon field in the framework of the Holstein Hubbard model in one, two and three dimensions. We have chosen the modified Lang-Firsov variational wave function with on-site and nearest-neighbour lattice distortion for the phonon subsystem to obtain an effective electronic Hamiltonian. This effective electronic Hamiltonian is then solved exactly to obtain the binding energy and effective mass of the bipolaron. We investigated the stability of the bipolaronic phase by comparing the ground-state energy of a bipolaron and two free polarons. We observed a critical repulsive on-site Hubbard interaction U c below which the bipolaronic phases are stable for a fixed electron-phonon coupling g . In the absence of on-site repulsion, bipolaronic phases are stable over the entire range of electron-phonon coupling for one dimension, whereas there is a critical electron-phonon coupling g c for formation of a stable bipolaron in two and three dimensions.

Journal of Physics Condensed Matter•Oct 27, 1999•5 citations

High- Tc superconductivity with polarons and bipolarons: an approach from the insulating states

S. Aubry

We review a series of exact, analytical and numerical results obtained on the adiabatic Holstein–Hubbard model, many of which are new and non-trivial. We study next the role of the quantum lattice fluctuations that were initially neglected. The possibility of having high- T c bipolaronic superconductors is analysed on the basis of these results. We suggest that models that involve only electron–phonon coupling are very unlikely to produce bipolaronic superconductivity, which is prevented by the spatial ordering of the bipolarons associated with a lattice instability. This is due to a very large effective mass of the bipolarons that can be related to the large Peierls–Nabarro energy barrier required to move these bipolarons through the lattice in the adiabatic limit. We conjecture that high- T c superconductivity originates specifically from the exceptionally well-balanced competition between electron–phonon coupling and electron–electron repulsion. In the restricted region in which, within the mean field approach, the energy of a bipolaron is close to those of two polarons, a new type of electron pairing occurs by formation of pairs of polarons in the spin singlet state. Such a polaron pair, called a spin resonant bipolaron (SRB), is not the standard bipolaron (both could exist in the adiabatic case). Its Peierls–Nabarro energy barrier can be sharply depressed, almost to zero. As a result of quantum lattice fluctuations, its tunnelling energy is sharply enhanced. Then, superconductivity could persist for a relatively large electron–phonon coupling, with an unusually large critical temperature as a Bose condensate of SRBs before becoming, at larger coupling and in any case, an insulating spin–Peierls polaronic phase.

Cambridge University Press eBooks•Sep 7, 1995•7 citations

Polarons and bipolarons in strongly interacting electron-phonon systems

Gerhard Wellein, Heinrich Röder, Holger Fehske

The Holstein Hubbard and Holstein t-J models are studied for a wide range of phonon frequencies, electron-electron, and electron-phonon interaction strengths on finite lattices with up to ten sites by means of direct Lanczos diagonalization. Previously the necessary truncation of the phononic Hilbert space caused serious limitations to either very small systems (four or even two sites) or to weak electron-phonon coupling, in particular in the adiabatic regime. Using parallel computers we were able to investigate the transition from ``large'' to ``small'' polarons in detail. By resolving the low-lying eigenstates of the Hamiltonian and by calculating the spectral function, we can identify a polaron band in the strong-coupling case, whose dispersion deviates from the free-particle dispersion at low and intermediate phonon frequencies. For two electrons (holes) we establish the existence of bipolaronic states and discuss the formation of a bipolaron band. For the two-dimensional Holstein t-J model, we demonstrate that the formation of hole polarons is favored by strong Coulomb correlations. Analyzing hole-hole correlation functions, we find that hole binding is enhanced as a dynamical effect of the electron-phonon interaction. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9512143•Physical review. B, Condensed matter•Apr 15, 1996•163 citations

Bipolaronic charge density waves and polaronic spin density waves in the Holstein-Hubbard model

S. Aubry

At large electron-phonon coupling, the phonon of the Holstein-Hubbard model can be (generally) treated classically. Exact results can be obtained for this model at any band filling and any finite dimension where the Hubbard term with coefficient U is treated rigorously. The existence of chaotic bipolaronic states is proven when U is smaller than a certain value Uc and for an electron-phonon coupling k sufficiently large. The existence of magnetic polaronic states is also proven for positive U and k sufficiently large. For 0<U<Uc, there exists an overlap region with chaotic mixed polaronic-bipolaronic states. There is a region for U=Uc/2 and large k of special interest, where the bipolarons in their ground-state become singlet pairs of polarons bounded by a magnetic resonance. Then, the Peierls-Nabarro barrier which usually pins polarons and bipolarons to the lattice depresses to almost zero and the validity of the classical lattice approximation breaks down so that the initial exact results do not apply in that specific region. The formation of a superconducting state with possibly a high critical temperature, is then conjectured at a large coupling k in a narrow domain in-between the two insulating domains corresponding to Bipolaronic Charge Density Waves and to magnetic Polaronic Spin Density Waves respectively.

Journal de Physique IV (Proceedings)•Jul 1, 1993•7 citations
Abstract. The spatially localized bound states of two electrons in the adiabatic two-dimensional Holstein-Hubbard model on a square lattice are investigated both numerically and analytically. The interplay between the electron-phonon coupling g, which tends to form bipolarons and the repulsive Hubbard interaction υ ≥ 0, which tends to break them, generates many different ground-states. There are four domains in the g,υ phase diagram delimited by first order transition lines. Except for the domain at weak electronphonon coupling (small g) where the electrons remain free, the electrons form bipolarons which can 1) be mostly located on a single site (small υ, large g); 2) be an anisotropic pair of polarons lying on two neighboring sites in the magnetic singlet state (large υ, large g); or 3) be a ”quadrisinglet state ” which is the superposition of 4 electronic singlets with a common central site. This quadrisinglet bipolaron is the most stable in a small central domain in between the three other phases. The pinning modes and the Peierls-Nabarro barrier of each of these bipolarons are calculated and the barrier is found to be strongly depressed in the region of stability of the quadrisinglet bipolaron. PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) – 71.38.+i Polarons and electron phonon interactions – 74.20.Mn Nonconventional mechanisms (spin fluctuations, polarons and bipolarons, resonating valence bond model, anyon mechanism, marginal Fermi liquid, Luttinger liquid, etc.) – 74.25.Jb Electronic structure
The bound states of two electrons in the adiabatic Holstein-Hubbard model are studied numerically in one and two dimensions from the anticontinuous limit. This model involves a competition between a local electron-phonon coupling (with a classical lattice) which tends to form pairs of electrons and the repulsive Hubbard interaction U≥0U \geq 0U≥0 which tends to break them. In 1D, the ground-state always consists in a pair of localized polarons in a singlet state. They are located at the same site for U=0. Increasing U, there is a first order transition at which the bipolaron becomes a spin singlet pair of two polarons bounded by a magnetic interaction. The pinning mode of the bipolaron soften in the vicinity of this transition leading to a higher mobility of the bipolaron which is tested numerically. In 2D, and for any UUU, the electron-phonon coupling needs to be large enough in order to form small polarons or bipolarons instead of extended electrons. We calculate the phase diagram of the bipolaron involving first order transitions lines with a triple point. A pair of polarons can form three types of bipolarons: a) on a single site at small UUU, b) a spin singlet state on two nearest neighbor sites for larger UUU as in 1D and c) a new intermediate state obtained as the resonant combination of four 2-sites singlet states sharing a central site, called quadrisinglet. The breathing and pinning internal modes of bipolarons in 2D generally only weakly soften and thus, they are practically not mobile. On the opposite, in the vicinity of the triple point involving the quadrisinglet, both modes exhibit a significant softening. However, it was not sufficient for allowing the existence of a classical mobile bipolaron (at least in that model).
We present a diagrammatic Monte Carlo study of the properties of the Hubbard-Holstein bipolaron on a two-dimensional square lattice. With a small Coulomb repulsion UUU and with increasing electron-phonon interaction, and when reaching a value about two times smaller than the one corresponding to the transition of light polaron to heavy polaron, the system suffers a sharp transition from a state formed by two weakly bound light polarons to a heavy, strongly bound on-site bipolaron. Aside from this rather conventional bipolaron a new bipolaron state is found for large UUU at intermediate and large electron-phonon coupling, corresponding to two polarons bound on nearest-neighbor sites. We discuss both the properties of the different bipolaron states and the transition from one state to another. We present a phase diagram in parameter space defined by the electron-phonon coupling and UUU. Our numerical method does not use any artificial approximation and can be easily modified to other bipolaron models with longer range electron-phonon and∕or electron-electron interaction.
In this paper we have studied the properties of a two-fermion system interacting with the phonon field in the framework of the Holstein Hubbard model in one, two and three dimensions. We have chosen the modified Lang-Firsov variational wave function with on-site and nearest-neighbour lattice distortion for the phonon subsystem to obtain an effective electronic Hamiltonian. This effective electronic Hamiltonian is then solved exactly to obtain the binding energy and effective mass of the bipolaron. We investigated the stability of the bipolaronic phase by comparing the ground-state energy of a bipolaron and two free polarons. We observed a critical repulsive on-site Hubbard interaction U c below which the bipolaronic phases are stable for a fixed electron-phonon coupling g . In the absence of on-site repulsion, bipolaronic phases are stable over the entire range of electron-phonon coupling for one dimension, whereas there is a critical electron-phonon coupling g c for formation of a stable bipolaron in two and three dimensions.
We review a series of exact, analytical and numerical results obtained on the adiabatic Holstein–Hubbard model, many of which are new and non-trivial. We study next the role of the quantum lattice fluctuations that were initially neglected. The possibility of having high- T c bipolaronic superconductors is analysed on the basis of these results. We suggest that models that involve only electron–phonon coupling are very unlikely to produce bipolaronic superconductivity, which is prevented by the spatial ordering of the bipolarons associated with a lattice instability. This is due to a very large effective mass of the bipolarons that can be related to the large Peierls–Nabarro energy barrier required to move these bipolarons through the lattice in the adiabatic limit. We conjecture that high- T c superconductivity originates specifically from the exceptionally well-balanced competition between electron–phonon coupling and electron–electron repulsion. In the restricted region in which, within the mean field approach, the energy of a bipolaron is close to those of two polarons, a new type of electron pairing occurs by formation of pairs of polarons in the spin singlet state. Such a polaron pair, called a spin resonant bipolaron (SRB), is not the standard bipolaron (both could exist in the adiabatic case). Its Peierls–Nabarro energy barrier can be sharply depressed, almost to zero. As a result of quantum lattice fluctuations, its tunnelling energy is sharply enhanced. Then, superconductivity could persist for a relatively large electron–phonon coupling, with an unusually large critical temperature as a Bose condensate of SRBs before becoming, at larger coupling and in any case, an insulating spin–Peierls polaronic phase.
The Holstein Hubbard and Holstein t-J models are studied for a wide range of phonon frequencies, electron-electron, and electron-phonon interaction strengths on finite lattices with up to ten sites by means of direct Lanczos diagonalization. Previously the necessary truncation of the phononic Hilbert space caused serious limitations to either very small systems (four or even two sites) or to weak electron-phonon coupling, in particular in the adiabatic regime. Using parallel computers we were able to investigate the transition from ``large'' to ``small'' polarons in detail. By resolving the low-lying eigenstates of the Hamiltonian and by calculating the spectral function, we can identify a polaron band in the strong-coupling case, whose dispersion deviates from the free-particle dispersion at low and intermediate phonon frequencies. For two electrons (holes) we establish the existence of bipolaronic states and discuss the formation of a bipolaron band. For the two-dimensional Holstein t-J model, we demonstrate that the formation of hole polarons is favored by strong Coulomb correlations. Analyzing hole-hole correlation functions, we find that hole binding is enhanced as a dynamical effect of the electron-phonon interaction. \textcopyright{} 1996 The American Physical Society.
At large electron-phonon coupling, the phonon of the Holstein-Hubbard model can be (generally) treated classically. Exact results can be obtained for this model at any band filling and any finite dimension where the Hubbard term with coefficient U is treated rigorously. The existence of chaotic bipolaronic states is proven when U is smaller than a certain value Uc and for an electron-phonon coupling k sufficiently large. The existence of magnetic polaronic states is also proven for positive U and k sufficiently large. For 0<U<Uc, there exists an overlap region with chaotic mixed polaronic-bipolaronic states. There is a region for U=Uc/2 and large k of special interest, where the bipolarons in their ground-state become singlet pairs of polarons bounded by a magnetic resonance. Then, the Peierls-Nabarro barrier which usually pins polarons and bipolarons to the lattice depresses to almost zero and the validity of the classical lattice approximation breaks down so that the initial exact results do not apply in that specific region. The formation of a superconducting state with possibly a high critical temperature, is then conjectured at a large coupling k in a narrow domain in-between the two insulating domains corresponding to Bipolaronic Charge Density Waves and to magnetic Polaronic Spin Density Waves respectively.