Bond Bipolarons: Sign-free Monte Carlo Approach

Polarons originating from phonon displacement modulated hopping have relatively light masses and, thus, are of significant current interest as candidates for bipolaron mechanism of high-temperature superconductivity [Phys. Rev. Lett. {\bf 121}, 247001 (2018)]. We observe that the bond model, when the dominant coupling comes from atomic vibrations on lattice bonds, can be solved by efficient sign-free Monte Carlo methods based on the path-integral formulation of the particle sector in combination with either the (real-space) diagrammatic or Fock-path-integral representation of the phonon sector. We introduce the corresponding algorithms and provide illustrative results for bipolarons in two dimensions. The results suggest that the route towards high-temperature superconductivity (if any) in the multiparametric space of the model lies between the Scylla of large size of moderately light bipolarons and Charybdis of large mass of compact bipolarons. As a result, on-site repulsion is helping $s$-wave superconductivity in sharp contrast with existing expectations.