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.
Most nonperturbative numerical studies of electron-phonon interactions focus on model Hamiltonians where the electrons interact with a phonon branch via a single type of microscopic mechanism. Two commonly explored couplings in this context are the Holstein and Su-Schrieffer-Heeger (SSH) interactions, which describe phonons modulating the on-site energy and intersite electron hopping, respectively. Many materials, however, have multiple phonon branches that can each interact with electronic degrees of freedom in different ways. We present here a determinant quantum Monte Carlo study of the half-filled two-dimensional (bond) SSH-Holstein Hamiltonian, where electrons couple to different phonon branches via either the Holstein or SSH mechanism. We map the model's phase diagram and determine the nature of the transitions between charge-density wave, bond-order wave, and antiferromagnetic order.
Electron-phonon (e-ph) interactions play a crucial role in determining many properties of materials. In this context, the Su-Schrieffer-Heeger (SSH) model, where atomic motion modulates the electronic hopping, has gained significant attention due to its potential for strong electron pairing in relation to high-Tc superconductivity. Previous studies of the SSH models have addressed many aspects of this problem, but have focused heavily on either dilute or half-filled models with dispersionless (Einstein) phonons. Here, we study the effects of dispersive optical phonons on the lightly doped one-dimensional optical Hubbard-SSH model using the density matrix renormalization group. We observe a significant enhancement in singlet binding driven by phonon dispersion; however, by calculating various correlation functions, we find that the enhanced binding does not translate to increased superconducting correlations but rather robust bond correlations in the studied parameter regime. Nevertheless, the significant impact of phonon dispersion on these correlations highlights the need to go beyond the Einstein phonon limit while modeling realistic quantum materials.
The Su-Schrieffer-Heeger (SSH) model, with bond phonons modulating electron tunneling, is a paradigmatic electron-phonon model that hosts an antiferromagnetic order to bond order transition at half-filling. In the presence of a repulsive Hubbard interaction, the antiferromagnetic phase is enhanced, but the phase transition remains first order. Here, we explore the physics of the SSH model with an attractive Hubbard interaction, which hosts an interesting interplay among charge order,
s-wave pairing, and bond order. Using the numerically exact determinant quantum Monte Carlo method, we show that both charge order, present at weak electron-phonon coupling, and bond order, at large coupling, give way to
s-wave pairing when the system is doped. Furthermore, we demonstrate that the SSH electron-phonon interaction competes with the attractive Hubbard interaction and reduces the
s-wave pairing correlation.
We study the Hubbard model on the square lattice coupled in addition to the optical Su-Schrieffer-Heeger (SSH) phonons, using the singular-mode functional renormalization group method. At half filling and in the absence of the Hubbard interaction
U, we find the degenerate spin-density-wave--charge-density-wave--
s-wave-superconductivity state at smaller electron-phonon coupling strength $\ensuremath{\lambda}$ and higher phonon frequency $\ensuremath{\omega}$, and the valence bond solid (VBS) state at larger $\ensuremath{\lambda}$ and lower $\ensuremath{\omega}$. After switching on a positive
U, the VBS state is suppressed, while the spin-density wave (SDW) state is enhanced. At finite doping, the SSH phonon is found to favor
s-wave superconductivity (sSC) at
U=0. With increasing positive
U, we find
d-wave superconductivity (dSC) and incommensurate SDW states. In a narrow window of moderate
U and $\ensuremath{\lambda}$, we also find the incommensurate VBS state. The sSC and dSC here can be naturally related to the charge-density wave (CDW) and SDW fluctuations, both of which can be triggered by the SSH phonons. In contrast, the repulsive interaction
U enhances SDW but suppresses CDW fluctuations. Our results at half filling are consistent with quantum Monte Carlo (QMC) and provide insights at finite doping where QMC may suffer from the minus sign problem.
Abstract The electron–lattice interaction gives rise to a rich set of phenomena in quantum materials. Microscopically, this interaction often arises from the modulation of orbital overlaps; however, many theoretical studies neglect such couplings. Here, we present an exact diagonalization and determinant quantum Monte Carlo study of a three-orbital Su–Schrieffer–Heeger (SSH) model, on a two-dimensional Lieb lattice and in the negative charge transfer regime. At half-filling (one hole/unit cell), we observe a bipolaron insulating phase with a bond-disproportionate lattice. This phase is robust against moderate hole doping but is suppressed at large hole concentrations, leading to a metallic polaron-liquid-like state with fluctuating patches of local distortions. We also find an s -wave superconducting state at large hole doping that primarily appears on the oxygen sublattice. Our work provides a non-perturbative view of SSH-type couplings in two dimensions with implications for materials where such couplings are dominant.