We analyze the transition temperature ${T}_{c}$ of bipolaronic superconductivity in a bond Su-Schrieffer-Heeger (bond-SSH) model---also known as a bond Peierls model---where bond phonons modulate the electron hoppings. Using a semiclassical instanton approximation justifiable in the adiabatic limit of slow phonons, we find that the bipolaron mass is only weakly enhanced in contrast to the typically large mass enhancement found in models where the phonons couple to electron density. Strikingly, in the ``atomic'' limit, we find a bond bipolaron can freely slide within a continuously degenerate energy manifold, explaining the origin of the weak mass enhancement in the strong coupling limit. A gas of these bipolarons can undergo a superfluid transition at a critical temperature for which we obtain an upper bound. This bound is exponentially larger than that in the Holstein model. Our study provides an analytical understanding of the mechanism behind the high-${T}_{c}$ bipolaronic superconductivity numerically observed in [Phys. Rev. X 13, 011010 (2023)].
We analyze the transition temperature Tc of bipolaronic superconductivity in a bond Su-Schrieffer-Heeger (bond-SSH) model---also known as a bond Peierls model---where bond phonons modulate the electron hoppings. Using a semiclassical instanton approximation justifiable in the adiabatic limit of slow phonons, we find that the bipolaron mass is only weakly enhanced in contrast to the typically large mass enhancement found in models where the phonons couple to electron density. Strikingly, in the ``atomic'' limit, we find a bond bipolaron can freely slide within a continuously degenerate energy manifold, explaining the origin of the weak mass enhancement in the strong coupling limit. A gas of these bipolarons can undergo a superfluid transition at a critical temperature for which we obtain an upper bound. This bound is exponentially larger than that in the Holstein model. Our study provides an analytical understanding of the mechanism behind the high-Tc bipolaronic superconductivity numerically observed in [Phys. Rev. X 13, 011010 (2023)].
Using numerically exact diagrammatic Monte Carlo simulations in the two-electron (single-bipolaron) sector, we explore the impact of long-range Coulomb repulsion on the dilute-limit Berezinskii-Kosterlitz-Thouless transition temperature T c of bipolarons on a two-dimensional square lattice. We study the bond Su-Schrieffer-Heeger model, in which bond phonons modulate the electron hopping. In the absence of long-range repulsion, this model was shown to support small, light bipolarons with a comparatively high transition temperature [C. Zhang, J. Sous, D. R. Reichman, M. Berciu, A. J. Millis, N. V. Prokof’ev, and B. V. Svistunov, Bipolaronic high-temperature superconductivity, ]. Here, we find that long-range Coulomb repulsion suppresses the optimal T c but leaves it appreciable over a broad parameter window, including the adiabatic regime ω / t = 0.5 at a representative Coulomb strength V = U / 10 (with U the on-site repulsion). Our results provide controlled single-bipolaron inputs for dilute-limit T c estimates in the presence of long-range repulsion.
Light-mass bipolarons in off-diagonally coupled electron-phonon systems provide a potential route to bipolaronic high-Tc superconductivity. While there has been numerical progress in the physically relevant limit of slow phonons, more insights are needed to fully understand to what extent this mechanism survives at finite densities and in different regimes. We address these questions using advanced tensor-network methods applied to the Su-Schrieffer-Heeger model. Studying both the spectral properties of isolated bipolarons as well as the pairing correlations at small finite densities, we find evidence that the conventional picture of bipolarons as molecular bound states which undergo Bose-Einstein condensation may need to be reconsidered. Instead, our findings suggest correlation-driven formation of a fragmented condensate with spatially separated polaron pairs, stabilized by strong repulsive electron-electron interactions at moderate values of the electron-phonon coupling. These spatially modulated charge clouds exhibit pairing with a typical length-scale set by the Fermi momentum.
Characterizing bipolaron binding, and understanding how it depends on electron-phonon interaction, is crucial to unraveling the nature of emergent many-body states in strongly interacting electron-phonon systems. So far, most studies of bipolarons have been limited to the Holstein model, in which the coupling constant is momentum-independent. The paradigmatic example of momentum-dependent electron-phonon interaction comes from the system in which phonon distortions modify electron hopping, the SSH model. Already individual polarons in the SSH model are richer than the Holstein model counterparts, and feature a phase transition into the finite momentum ground state with increasing electron-phonon interaction. In this paper, we use a variational approach to study bipolarons in the one-dimensional SSH model and discuss their ground state, dispersion, and excitation spectra. We explore the full parameter range of the system, including the adiabatic regime of slow phonons, which was inaccessible to previous theoretical studies. In agreement with earlier studies, we find that in the anti-adiabatic strongly interacting regime, bipolarons have low effective mass. By contrast, in the adiabatic case, we find that increasing electron-phonon interactions results in an exponential increase of the bipolaron mass. We establish the existence of multiple branches of bound excited states of SSH bipolaron and discuss the signatures of these bound states in dynamics. We show that in the anti-adiabatic regime, response functions obey a parity selection rule, that imposes symmetry constraints on the excitation spectra and provides a clear signature of SSH bipolarons.
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.
Theoretical studies suggest that Su-Schrieffer-Heeger-like electron-phonon (e-ph) interactions can mediate high-temperature bipolaronic superconductivity that is robust against repulsive electron-electron interactions. Here we present a comparative analysis of the pairing and competing charge/bond correlations in the two-dimensional Holstein and optical Su-Schrieffer-Heeger (SSH) models using numerically exact determinant quantum Monte Carlo. We find that the SSH interactions support light bipolarons and strong superconducting correlations out to relatively large values of the e-ph coupling λ and densities near half-filling, while the Holstein interaction does not due to the formation of heavy bipolarons and competing charge-density-wave order. We further find that the Holstein and SSH models have comparable pairing correlations in the weak coupling limit for carrier concentrations ⟨n⟩≪1, where competing orders and polaronic effects are absent. These results support the proposal that SSH (bi)polarons can support superconductivity to larger values of λ in comparison to the Holstein polaron, but that the resulting Tc gains are small in the weak coupling limit. We also find that the SSH model's pairing correlations are suppressed after including a weak on-site Hubbard repulsion. These results have important implications for identifying and engineering bipolaronic superconductivity.
Su-Schrieffer-Heeger (SSH) electron-phonon ($e\ensuremath{-}p\phantom{\rule{0}{0ex}}h$) interactions have been theorized to play critical roles in several novel states of matter, ranging from nontrivial topological states to high-temperature bipolaronic superconductivity. This work compares the superconducting and competing charge and bond correlations of the two-dimensional Holstein and optical (SSH) models. The authors find that near half-filling, light SSH (bi)polarons support superconductivity to larger values of $e\ensuremath{-}p\phantom{\rule{0}{0ex}}h$ coupling compared to the Holstein polaron. These results are essential for identifying and engineering bipolaronic superconductivity in quantum materials.