Andy Tanjaroon Ly

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 $\lambda$ 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 $\langle n\rangle \ll 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 $\lambda$ in comparison to the Holstein polaron, but that the resulting $T_{\mathrm{c}}$ 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.

We compare the acoustic Su–Schrieffer–Heeger (SSH) model with two of its optical variants where the phonons are defined on either on the sites or bonds of the system. First, we discuss how to make fair comparisons between these models in any dimension by ensuring their dimensionless coupling $\lambda$ and relevant phonon energies are the same. We then use determinant quantum Monte Carlo to perform nonperturbative and sign-problem-free simulations of all three models on one-dimensional chains at and away from half-filling. By comparing the results obtained from each model, we demonstrate that the optical and acoustic models produce near identical results within error bars for suitably chosen phonon energies and $\lambda$ at half-filling. In contrast, the bond model has quantitatively different behavior due to its coupling to the $\mathbf{q}=0$ phonon mode. These differences also manifest in the total length of the chain, which shrinks for the bond model but not for the acoustic and optical models when $\lambda \ne 0$. Our results have important implications for quantum Monte Carlo modeling of SSH-like interactions, where these models are sometimes regarded as being interchangeable.

Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon (e-ph) coupling, and it is unclear how this interaction influences stripe correlations. We study static and fluctuating stripe orders in the doped single-band Hubbard-Holstein model using zero temperature variational Monte Carlo and finite temperature determinant quantum Monte Carlo. We find that the lattice couples more strongly with the charge component of the stripes, leading to an enhancement or suppression of stripe correlations, depending on model parameters like the next-nearest-neighbor hopping $${t}^{{\prime} }$$t′ or phonon energy Ω. Our results help elucidate how thee-ph interaction can tip the delicate balance between stripe and superconducting correlations in the Hubbard-Holstein model with implications for our understanding of the high-T_ccuprates.