Mobile impurities in a Bose-Einstein condensate form quasiparticles called polarons. Here, we show that two such polarons can bind to form a bound bipolaron state. Its emergence is caused by an induced nonlocal interaction mediated by density oscillations in the condensate, and we derive using field theory an effective Schrödinger equation describing this for an arbitrarily strong impurity-boson interaction. We furthermore compare with quantum Monte Carlo simulations finding remarkable agreement, which underlines the predictive power of the developed theory. It is found that bipolaron formation typically requires strong impurity interactions beyond the validity of more commonly used weak-coupling approaches that lead to local Yukawa-type interactions. We predict that the bipolarons are observable in present experiments, and we describe a procedure to probe their properties.
Mobile impurities in a Bose-Einstein condensate form quasiparticles called polarons. Here, we show that two such polarons can bind to form a bound bipolaron state. Its emergence is caused by an induced nonlocal interaction mediated by density oscillations in the condensate, and we derive using field theory an effective Schrödinger equation describing this for an arbitrarily strong impurity-boson interaction. We furthermore compare with quantum Monte Carlo simulations finding remarkable agreement, which underlines the predictive power of the developed theory. It is found that bipolaron formation typically requires strong impurity interactions beyond the validity of more commonly used weak-coupling approaches that lead to local Yukawa-type interactions. We predict that the bipolarons are observable in present experiments, and we describe a procedure to probe their properties.
Bose polarons, quasiparticles 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 nonperturbative theory and complementing it with exact numerical simulations. We develop an analytic approach for weak boson-boson interactions and arbitrarily strong impurity-boson couplings. Our approach 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 the Born-Oppenheimer approximation, leading to an 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 the leading-order Born-Huang corrections, we calculate bipolaron binding energies for impurity-boson mass ratios as low as 3 and find excellent agreement with QMC results.
Quasiparticles and their interactions are a key part of our understanding of quantum many-body systems. Quantum simulation experiments with cold atoms have in recent years advanced our understanding of isolated quasiparticles, but so far they have provided limited information regarding their interactions and possible bound states. Here, we show how exploring mobile impurities immersed in a Bose-Einstein condensate (BEC) in a two-dimensional lattice can address this problem. First, the spectral properties of individual impurities are examined, and in addition to the attractive and repulsive polarons known from continuum gases, we identify a new kind of quasiparticle stable for repulsive boson-impurity interactions. The spatial properties of polarons are calculated showing that there is an increased density of bosons at the site of the impurity both for repulsive and attractive interactions. We then derive an effective Schrödinger equation describing two polarons interacting via the exchange of density oscillations in the BEC, which takes into account strong impurity-boson two-body correlations. Using this, we show that the attractive nature of the effective interaction between two polarons combined with the two-dimensionality of the lattice leads to the formation of bound states - i.e. bipolarons. The wave functions of the bipolarons are examined showing that the ground state is symmetric under particle exchange and therefore relevant for bosonic impurities, whereas the first excited state is doubly degenerate and odd under particle exchange making it relevant for fermionic impurities. Our results show that quantum gas microscopy in optical lattices is a promising platform to explore the spatial properties of polarons as well as to finally observe the elusive bipolarons.
The variational Feynman formalism for the polaron, extended to an all-coupling treatment of bipolarons, is applied for two impurity atoms in a Bose-Einstein condensate. This shows that if the polaronic coupling strength is large enough, the impurities will form a bound state (the bipolaron). As a function of the mutual repulsion between the impurities, two types of bipolaron are distinguished: a tightly bound bipolaron at weak repulsion and a dumbbell bipolaron at strong repulsion. Apart from the binding energy, the evolution of the bipolaron radius and its effective mass are also examined as a function of the strength of the repulsive interaction between the impurities and of the polaronic coupling strength. We then apply the strong-coupling formalism to multiple-impurity atoms in a condensate, which leads to the prediction of multipolaron formation in the strong-coupling regime. The results of the two formalisms are compared for two impurities in a condensate, which results in a general qualitative agreement and a quantitative agreement at strong coupling. Typically, the system of impurity atoms in a Bose-Einstein condensate is expected to exhibit the polaronic weak-coupling regime. However, the polaronic coupling strength is, in principle, tunable with a Feshbach resonance.
Impurities immersed in a Bose-Einstein condensate (BEC) can interact indirectly through the exchange of Bogoliubov excitations. These impurities, which form dressed quasiparticles known as Bose polarons due to their interaction with the BEC, can pair up to form a bound state called a bipolaron, via an induced interaction. Previous studies on induced interactions have focused primarily on cases with an isotropic excitation spectrum. In this work, we investigate the properties of induced interactions and bipolarons mediated by anisotropic Bogoliubov excitations using field theory. Taking a BEC with spin-orbit coupling as an example, we show that the induced interaction becomes anisotropic. Notably, a double-minima feature appears in the induced interaction in momentum space due to the exchange of roton excitations. Additionally, we calculate the binding energy and wave functions of these bipolarons induced by anisotropic interactions. Unlike previously studied bipolarons formed through the exchange of isotropic phonon excitations, we identify an alternative type of bipolaron whose wave functions feature a double-peak structure under strong impurity-boson interactions. Our work extends the theory of induced interactions from isotropic to anisotropic systems, and it reveals the novel features in both the induced interactions and bipolarons arising from BEC with an unconventional excitation spectrum.
Cold atom experiments show that a mobile impurity particle immersed in a weakly interacting Bose-Einstein condensate forms a well-defined quasiparticle (Bose polaron) for weak to moderate impurity-boson interaction strengths, whereas a significant line broadening is consistently observed for strong interactions. Motivated by this, we introduce a phenomenological theory based on the assumption that the most relevant states are characterized by the impurity correlated with at most one boson, since they have the largest overlap with the uncorrelated states to which the most common experimental probes couple. These experimentally relevant states can, however, decay to lower energy states characterized by correlations involving multiple bosons, and we model this using a minimal variational wave function combined with a complex impurity-boson interaction strength. We first motivate this approach by comparing to a more elaborate theory that includes correlations with up to two bosons. Our phenomenological model is shown to recover the main results of two recent experiments probing both the spectral and the nonequilibrium properties of the Bose polaron. Our work offers an intuitive framework for analyzing experimental data and highlights the importance of understanding the complicated problem of the Bose polaron decay in a many-body setting.
The mobile impurity in a Bose-Einstein condensate (BEC) is a paradigmatic many-body problem. For weak interaction between the impurity and the BEC, the impurity deforms the BEC only slightly and it is well described within the Fröhlich model and the Bogoliubov approximation. For strong local attraction this standard approach, however, fails to balance the local attraction with the weak repulsion between the BEC particles and predicts an instability where an infinite number of bosons is attracted toward the impurity. Here we present a solution of the Bose polaron problem beyond the Bogoliubov approximation which includes the local repulsion between bosons and thereby stabilizes the Bose polaron even near and beyond the scattering resonance. We show that the Bose polaron energy remains bounded from below across the resonance and the size of the polaron dressing cloud stays finite. Our results demonstrate how the dressing cloud replaces the attractive impurity potential with an effective many-body potential that excludes binding. We find that at resonance, including the effects of boson repulsion, the polaron energy depends universally on the effective range. Moreover, while the impurity contact is strongly peaked at positive scattering length, it remains always finite. Our solution highlights how Bose polarons are self-stabilized by repulsion, providing a mechanism to understand quench dynamics and nonequilibrium time evolution at strong coupling.