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:cond-mat/0401305 · Physical Review B

Two-dimensional Hubbard-Holstein bipolaron

Alexandru Macridin, G. A. Sawatzky, Mark Jarrell

Published Jun 23, 2004• 110 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.69.245111

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

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

Laurent Proville, S. Aubry

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

•Jan 1, 1998•27 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

Comprehensive study of bond bipolaron superconductivity on the triangular lattice

Chao Zhang

We employ the diagrammatic Monte Carlo method with a lattice path-integral formulation for both electron and phonon degrees of freedom to investigate the formation and properties of bond polarons and bipolarons on a two-dimensional triangular lattice. In the deep adiabatic regime ($\ensuremath{\omega}/t<1.0$), single polarons remain light with a small effective mass, while bipolarons remain compact and lightweight, resulting in a high superfluid transition temperature ${T}_{c}$. We systematically study the dependence of ${T}_{c}$ on electron-phonon coupling strength and on-site interaction in the bond bipolaron model. Our results show that a moderate on-site repulsion enhances ${T}_{c}$ by stabilizing compact yet lightweight bipolarons, leading to high-${T}_{c}$ superconductivity in the dilute limit. Notably, the triangular lattice sustains relatively high ${T}_{c}$ across a wide range of phonon frequencies, outperforming square lattice geometry. This enhancement arises from the higher coordination number and the bond-centered nature of the electron-phonon coupling. These findings suggest that triangular lattice geometry offers a promising platform for realizing high-${T}_{c}$ bipolaronic superconductivity.

Physical review. B./Physical review. B•Nov 25, 2025•2 citations

Bipolarons in the extended Holstein Hubbard model

J. Bonča, S. A. Trugman

We numerically and analytically calculate the properties of the bipolaron in an extended Hubbard-Holstein model, which has a longer-range electron-phonon coupling like the Fr\"ohlich model. In the strong-coupling regime, the effective mass of the bipolaron in the extended model is much smaller than the Holstein bipolaron mass. In contrast to the Holstein bipolaron, the bipolaron in the extended model has a lower binding energy and remains bound with substantial binding energy even in the large-$U$ limit. In comparison with the Holstein model where only a singlet bipolaron is bound, in the extended Holstein model a triplet bipolaron can also form a bound state. We discuss the possibility of phase separation in the case of finite electron doping.

arXiv:cond-mat/0103457•Physical review. B, Condensed matter•Aug 7, 2001•108 citations

Bipolarons from long-range interactions: Singlet and triplet pairs in the screened Hubbard-Fröhlich model on the chain

J. P. Hague, P. E. Kornilovitch

We present details of a continuous-time quantum Monte Carlo algorithm for the screened Hubbard-Fr\"ohlich bipolaron. We simulate the bipolaron in one dimension with arbitrary interaction range in the presence of Coulomb repulsion, computing the effective mass, binding energy, total number of phonons associated with the bipolaron, mass isotope exponent, and bipolaron radius in a comprehensive survey of the parameter space. We discuss the role of the range of the electron-phonon interaction, demonstrating the evolution from Holstein to Fr\"ohlich bipolarons and we compare the properties of bipolarons with singlet and triplet pairing. Finally, we present simulations of the bipolaron dispersion. The bandwidth of the Fr\"ohlich bipolaron is found to be broad and the decrease in bandwidth as the two polarons bind into a bipolaron is found to be far less rapid than in the case of the Holstein interaction. The properties of bipolarons formed from long-range electron-phonon interactions, such as light strongly bound bipolarons and intersite pairing when Coulomb repulsion is large, are found to be robust against screening, with qualitative differences between Holstein and screened Fr\"ohlich bipolarons found even for interactions screened within a single lattice site.

arXiv:0904.3678•Physical Review B•Aug 3, 2009•21 citations
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.
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).
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
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 employ the diagrammatic Monte Carlo method with a lattice path-integral formulation for both electron and phonon degrees of freedom to investigate the formation and properties of bond polarons and bipolarons on a two-dimensional triangular lattice. In the deep adiabatic regime ($\ensuremath{\omega}/t<1.0$), single polarons remain light with a small effective mass, while bipolarons remain compact and lightweight, resulting in a high superfluid transition temperature Tc{T}_{c}Tc​. We systematically study the dependence of Tc{T}_{c}Tc​ on electron-phonon coupling strength and on-site interaction in the bond bipolaron model. Our results show that a moderate on-site repulsion enhances Tc{T}_{c}Tc​ by stabilizing compact yet lightweight bipolarons, leading to high-Tc{T}_{c}Tc​ superconductivity in the dilute limit. Notably, the triangular lattice sustains relatively high Tc{T}_{c}Tc​ across a wide range of phonon frequencies, outperforming square lattice geometry. This enhancement arises from the higher coordination number and the bond-centered nature of the electron-phonon coupling. These findings suggest that triangular lattice geometry offers a promising platform for realizing high-Tc{T}_{c}Tc​ bipolaronic superconductivity.
We numerically and analytically calculate the properties of the bipolaron in an extended Hubbard-Holstein model, which has a longer-range electron-phonon coupling like the Fr\"ohlich model. In the strong-coupling regime, the effective mass of the bipolaron in the extended model is much smaller than the Holstein bipolaron mass. In contrast to the Holstein bipolaron, the bipolaron in the extended model has a lower binding energy and remains bound with substantial binding energy even in the large-UUU limit. In comparison with the Holstein model where only a singlet bipolaron is bound, in the extended Holstein model a triplet bipolaron can also form a bound state. We discuss the possibility of phase separation in the case of finite electron doping.
We present details of a continuous-time quantum Monte Carlo algorithm for the screened Hubbard-Fr\"ohlich bipolaron. We simulate the bipolaron in one dimension with arbitrary interaction range in the presence of Coulomb repulsion, computing the effective mass, binding energy, total number of phonons associated with the bipolaron, mass isotope exponent, and bipolaron radius in a comprehensive survey of the parameter space. We discuss the role of the range of the electron-phonon interaction, demonstrating the evolution from Holstein to Fr\"ohlich bipolarons and we compare the properties of bipolarons with singlet and triplet pairing. Finally, we present simulations of the bipolaron dispersion. The bandwidth of the Fr\"ohlich bipolaron is found to be broad and the decrease in bandwidth as the two polarons bind into a bipolaron is found to be far less rapid than in the case of the Holstein interaction. The properties of bipolarons formed from long-range electron-phonon interactions, such as light strongly bound bipolarons and intersite pairing when Coulomb repulsion is large, are found to be robust against screening, with qualitative differences between Holstein and screened Fr\"ohlich bipolarons found even for interactions screened within a single lattice site.