A large polaron is a quasiparticle that consists of a nearly free electron interacting with the phonons of a material, whose lattice parameters are much smaller than the polaron scale. The electron-phonon interaction also leads to an attractive interaction between electrons, which can allow two polarons to pair up and form a bipolaron. It has been shown that large bipolarons can form in theory due to strong one-electron--one-phonon coupling, but they have not been seen in real materials because the critical value of the required electron-phonon interaction is too large. Here, we investigate the effect of one-electron--two-phonon coupling on the large bipolaron problem. Starting from a generalization of the Fr\"ohlich Hamiltonian that includes both the standard one-electron--one-phonon interaction as well as an anharmonic one-electron--two-phonon interaction, we use the path-integral method to find a semianalytical upper bound for the bipolaron energy that is valid at all values of the Fr\"ohlich coupling strength $\ensuremath{\alpha}$. We find the bipolaron phase diagram and conditions for the bipolaron stability by comparing the bipolaron energy to the energy of two free polarons. The critical value of the Fr\"ohlich coupling strength ${\ensuremath{\alpha}}_{\text{crit}}$ is calculated as a function of the strength of the one-electron--two-phonon interaction. The results suggest that large bipolaron formation is more likely in materials with significant one-electron--two-phonon interaction as well as strong one-electron--one-phonon interaction, such as strontium titanate.
Read full abstract