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.
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.
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.
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.
Electron-phonon (e-ph) interactions arise in many strongly correlated quantum materials from the modulation of the nearest-neighbor hopping integrals, as in the celebrated Su-Schrieffer-Heeger (SSH) model. Nevertheless, relatively few non-perturbative studies of correlated SSH models have been conducted in dimensions greater than one, and those that have been done have primarily focused on bond models, where generalized displacements independently modulate each hopping integral. We conducted a sign-problem free determinant quantum Monte Carlo study of the optical SSH-Hubbard model on a two-dimensional square lattice, where site-centered phonon modes simultaneously modulate pairs of nearest-neighbor hopping integrals. We report the model's low-temperature phase diagram in the challenging adiabatic regime (Ω/EF∼1/8). It exhibits insulating antiferromagnetic Mott and bond-order-wave (BOW) phases with a narrow region of coexistence between them. We also find that a critical e-ph coupling is required to stabilize the BOW phase in the small U limit. Lastly, in stark contrast to recent findings for the model's bond variant, we find no evidence for a long-range antiferromagnetism in the pure (U/t=0) optical SSH model.
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.
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.
This thesis investigates the interplay between electron-electron and electron-phonon interactions in low-dimensional lattice models, specifically within the Hubbard model incorporating Su-Schrieffer-Heeger (SSH) coupling with phonons. SSH interactions, which modulate the electron hopping amplitude through lattice displacements, have recently attracted renewed interest due to their potential relevance in strongly correlated materials. The main objective of this work is to determine whether the inclusion of SSH couplings in the Hubbard model can stabilize superconducting or spin-gapped metallic states, as well as other exotic phases of matter. While the half-filled SSH-Hubbard model has been extensively studied in one dimension and only partially examined in two dimensions, its phase diagram in the doped regime remains largely unexplored. This thesis aims to fill this gap, focusing on the emergence of correlated phases away from half-filling in both one and two spatial dimensions. We use the variational Monte Carlo method with a trial wave function defined in a hybrid electron-phonon Hilbert space. In addition to Jastrow factors that account for correlation effects, the wave function includes an antisymmetric component that captures backflow-inspired electron-phonon correlations. This approach enables us to accurately model both electron-electron and electron-phonon interactions. Moreover, since the variational method avoids the sign problem, we can simulate the system in the doped regime, an area typically inaccessible to exact quantum Monte Carlo methods (QMC), which are limited to half-filling. In one dimension, we observe the stabilization of a Luther-Emery liquid, a metallic phase with gapless charge excitations and a finite spin gap, upon doping a Peierls insulating phase. This finding resolves a long-standing debate regarding its realization in SSH chains. Benchmarking with density-matrix renormalization group results confirms the reliability of our approach. In two dimensions, we investigate the SSH-Hubbard model on the square lattice, obtaining results that are consistent with previous exact QMC studies at half-filling. Upon doping, we observe that SSH phonons can induce s-wave superconductivity in the weakly correlated regime. However, we suggest that the enhancement of d-wave superconductivity by SSH coupling in the strongly correlated regime may significantly depend on the specific implementation of the SSH interaction. In one particular case, we find that d-wave superconductivity is slightly suppressed by SSH coupling. These results are further confirmed by a Hartree-Fock analysis performed in the anti-adiabatic limit. Overall, this thesis highlights the crucial role of lattice effects in correlated electron systems, shedding light on the distinct contributions of SSH-type couplings in one and two dimensions. Our findings provide new insights into the microscopic mechanisms behind exotic metallic phases and set the stage for further studies of phonon-mediated phenomena in correlated systems.