Strongly coupling bipolarons and polaron interactions in 1D Bose gases

Bose polarons, quasi-particles composed of mobile impurities surrounded by cold Bose gas, can experience strong interactions mediated by the many-body environment and form bipolaron bound states. Here we present a detailed study of heavy polarons in a one-dimensional Bose gas by formulating a non-perturbative theory and complementing it with exact numerical simulations. We developed a semi analytic approach for weak boson-boson interactions and arbitrarily strong impurity-boson couplings. It is based on a mean-field theory that accounts for deformations of the superfluid by the impurities and in this way minimizes quantum fluctuations. The mean-field equations are solved exactly in Born-Oppenheimer (BO) approximation leading to a semi analytic expression for the interaction potential of heavy polarons which is found to be in excellent agreement with quantum Monte-Carlo (QMC) results. In the strong-coupling limit the potential substantially deviates from the exponential form valid for weak coupling and has a linear shape at short distances. Taking into account lowest-order BO corrections we calculate bipolaron binding energies and find excellent agreement with QMC results for impurity-boson mass ratios as low as 3.