We present the complete ground-state phase diagram of the Holstein model in two and three dimension considering the phonon variables to be classical. We first establish the overall structure of the phase diagram by using exact-diagonalisation–based Monte Carlo (ED-MC) on small lattices and then use a new "travelling cluster" approximation (TCA) for annealing the phonon degrees of freedom on large lattices. The phases that emerge include a Fermi liquid (FL), with no lattice distortions, an insulating polaron liquid (PL) at strong coupling, and a charge-ordered insulating (COI) phase around half-filling. The COI phase is separated from the Fermi liquid by a regime of phase coexistence whose width grows with increasing electron-phonon coupling. We provide results on the electronic density of states, the COI order parameter, and the spatial organisation of polaronic states, for arbitrary density and electron-phonon coupling. The results highlight the crucial role of spatial correlations in this strong-coupling problem.
We present the complete ground-state phase diagram of the Holstein model in two and three dimension considering the phonon variables to be classical. We first establish the overall structure of the phase diagram by using exact-diagonalisation–based Monte Carlo (ED-MC) on small lattices and then use a new "travelling cluster" approximation (TCA) for annealing the phonon degrees of freedom on large lattices. The phases that emerge include a Fermi liquid (FL), with no lattice distortions, an insulating polaron liquid (PL) at strong coupling, and a charge-ordered insulating (COI) phase around half-filling. The COI phase is separated from the Fermi liquid by a regime of phase coexistence whose width grows with increasing electron-phonon coupling. We provide results on the electronic density of states, the COI order parameter, and the spatial organisation of polaronic states, for arbitrary density and electron-phonon coupling. The results highlight the crucial role of spatial correlations in this strong-coupling problem.
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.
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.
We provide a comprehensive solution to the lattice dynamics problem in the two dimensional Holstein model at finite electron density and finite temperature. We work in the physically relevant adiabatic regime and vary the electron-phonon interaction from the weak coupling perturbative window to the strong coupling polaronic regime. We explore three typical electron densities, dilute - where spatial correlations between polarons is weak, intermediate - where correlations are significant, and half-filling - where there is long range checkerboard order at low temperature. We use two methods both of which exploit the "slowness" of the phonons to handle the problem. These are (i)~a standard random phase approximation (RPA), adapted to capture small quantum fluctuations on Monte Carlo generated classical thermal backgrounds, and (ii)~a Langevin dynamics scheme, with a simplified ``thermal noise'', that can address large amplitude dynamical fluctuations. The Langevin scheme, as we argue in the paper, is the superior method in the strong coupling part of the phase diagram, where lattice distortions are large. It reveals a non trivial multi-peak momentum resolved spectrum with a high energy part, on the scale of the bare phonon frequency Ω, and a low energy peak at ω≪Ω. Below the polaronic threshold, the high energy dispersion changes only modestly with temperature T, while the broadening, arising from mode coupling, increases linearly with T at low temperature. The low energy peak shows up at strong coupling and finite temperature and arises from the slow tunneling of polarons. The tunneling events become spatially correlated as electron density increases towards half-filling, and the weight becomes strongly momentum and temperature dependent. We suggest the analytic basis of these results.
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 electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.
We investigate the magnetotransport properties of the one-band Hubbard-Holstein model at half-filling in three dimensions (3D) using exact diagonalization based semi-classical Monte Carlo simulations with phonons treated in the adiabatic limit. The low-temperature electronic correlation U vs electron-phonon coupling V phase diagram reveals two insulating phases--antiferromagnetic (AF-I) and charge-ordered (CO-I)--separated by a first-order transition, with no metallic phase observed at their intersection, indicating robustness of these phases in 3D. For U∼bandwidth, the V vs temperature T phase diagram exhibits multiple phases including AF-I, CO-I, Mott-Hubbard insulator, bipolaronic insulator, and two bipolaronic metallic states. Several first-order transitions occur near V∼3.75. Above the ordering temperature, the density of states shows universal behavior dominated by electronic contribution, while susceptibility and DOS analyses reveal pseudogap features. Magnetic and transport properties along the phase boundary highlight strong proximity effects between competing phases, suggesting routes for tuning correlated materials and emergent electronic states.