Functional renormalization group study of the two-dimensional Su-Schrieffer-Heeger-Hubbard model

Abstract

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 (SDW)/charge-density-wave (CDW)/s-wave superconductivity (sSC) state at smaller electron-phonon coupling strength $\lambda$ and higher phonon frequency $\omega$, and the valence bond solid (VBS) state at larger $\lambda$ and lower $\omega$. After switching on a positive $U$, the VBS state is suppressed, while the SDW state is enhanced. At finite doping, the SSH phonon is found to favor 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 $\lambda$, we also find the incommensurate VBS state. The sSC and dSC here can be naturally related to the 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.