We explore the properties of the bipolaron in a 1D Holstein-Hubbard model with dynamical quantum phonons. Using a recently developed variational method combined with analytical strong coupling calculations, we compute correlation functions, effective mass, bipolaron isotope effect and the phase diagram. The two site bipolaron has a significantly reduced mass and isotope effect compared to the on-site bipolaron, and is bound in the strong coupling regime up to twice the Hubbard U naively expected. The model can be described in this regime as an effective t-J-V model with nearest neighbor repulsion. These are the most accurate bipolaron calculations to date.We explore the properties of the bipolaron in a 1D Holstein-Hubbard model with dynamical quantum phonons. Using a recently developed variational method combined with analytical strong coupling calculations, we compute correlation functions, effective mass, bipolaron isotope effect and the phase diagram. The two site bipolaron has a significantly reduced mass and isotope effect compared to the on-site bipolaron, and is bound in the strong coupling regime up to twice the Hubbard U naively expected. The model can be described in this regime as an effective t-J-V model with nearest neighbor repulsion. These are the most accurate bipolaron calculations to date.

arXiv:cond-mat/0001325Jan 21, 20000 citations

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).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, 19980 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.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/0103457Physical review. B, Condensed matterAug 7, 2001108 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.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/0401305Physical Review BJun 23, 2004110 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.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 MatterOct 27, 19995 citations

Numerical study of the E⊗e Jahn-Teller polaron and bipolaron

Samir El Shawish, J. Bonča, Li-Chung Ku +1 more

The properties of the polaron and bipolaron are explored in the one-dimensional Jahn-Teller model with dynamical quantum phonons. The ground-state properties of the polaron and bipolaron are computed using a recently developed variational method. Dynamical properties of the ground state of a polaron are investigated by calculating the optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega}).$ Our numerical results suggest that the Jahn-Teller and Holstein polarons are similar. However, in the strong-coupling regime qualitative differences in $\ensuremath{\sigma}(\ensuremath{\omega})$ between the two models are found and discussed. The influence of the electron-phonon coupling and the electrostatic repulsion on the bipolaron binding energy, bipolaron masses, and correlation functions is investigated.The properties of the polaron and bipolaron are explored in the one-dimensional Jahn-Teller model with dynamical quantum phonons. The ground-state properties of the polaron and bipolaron are computed using a recently developed variational method. Dynamical properties of the ground state of a polaron are investigated by calculating the optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega}).$ Our numerical results suggest that the Jahn-Teller and Holstein polarons are similar. However, in the strong-coupling regime qualitative differences in $\ensuremath{\sigma}(\ensuremath{\omega})$ between the two models are found and discussed. The influence of the electron-phonon coupling and the electrostatic repulsion on the bipolaron binding energy, bipolaron masses, and correlation functions is investigated.

arXiv:cond-mat/0211099Physical review. B, Condensed matterJan 10, 200329 citations

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).

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.

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.