Impurities in a trapped one-dimensional Bose gas of arbitrary interaction strength: Localization-delocalization transition and absence of self-localization
Dennis Breu, Eric Vidal Marcos, Martin Will, Michael Fleischhauer
We discuss impurities in a one-dimensional Bose gas with arbitrary boson-boson and boson-impurity interactions. To fully account for quantum effects, we employ numerical simulations based on the density-matrix renormalization group (DMRG) and—in the regime of strong boson-boson interactions—the mapping to weakly interacting fermions. A mean-field description of the Bose polaron based on coupled Gross–Pitaevskii Schrödinger equations predicts the existence of a self-localized polaron. We here show that such a solution does not exist and is an artifact of the underlying decoupling approximation. To this end we consider a mobile impurity in a box potential. Our work demonstrates that correlations between the impurity position and the bosons are important even in the limit where mean-field approaches are expected to work well. Furthermore, we derive analytical approximations for the energy of a single polaron formed by a heavy impurity for arbitrary interaction strengths and large but finite boson-boson couplings which accurately reproduce DMRG results. This demonstrates that the polaron problem of a heavy impurity in a one-dimensional Bose gas can be accurately approximated by a proper mean-field description plus a linearized treatment of quantum fluctuations for arbitrary boson-boson and impurity-boson couplings. Finally, we determine the polaron-polaron interaction potential V ( r ) in the Born-Oppenheimer approximation for small and intermediate distances r , which in the Tonks gas limit is oscillatory due to Friedel oscillations in the Bose gas.
We discuss impurities in a one-dimensional Bose gas with arbitrary boson-boson and boson-impurity interactions. To fully account for quantum effects, we employ numerical simulations based on the density-matrix renormalization group (DMRG) and—in the regime of strong boson-boson interactions—the mapping to weakly interacting fermions. A mean-field description of the Bose polaron based on coupled Gross–Pitaevskii Schrödinger equations predicts the existence of a self-localized polaron. We here show that such a solution does not exist and is an artifact of the underlying decoupling approximation. To this end we consider a mobile impurity in a box potential. Our work demonstrates that correlations between the impurity position and the bosons are important even in the limit where mean-field approaches are expected to work well. Furthermore, we derive analytical approximations for the energy of a single polaron formed by a heavy impurity for arbitrary interaction strengths and large but finite boson-boson couplings which accurately reproduce DMRG results. This demonstrates that the polaron problem of a heavy impurity in a one-dimensional Bose gas can be accurately approximated by a proper mean-field description plus a linearized treatment of quantum fluctuations for arbitrary boson-boson and impurity-boson couplings. Finally, we determine the polaron-polaron interaction potential V ( r ) in the Born-Oppenheimer approximation for small and intermediate distances r , which in the Tonks gas limit is oscillatory due to Friedel oscillations in the Bose gas.
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.
We study a weakly interacting one-dimensional Bose gas with two impurities coupled locally to the boson density. We derive analytical results for the induced interaction between the impurities at arbitrary positive coupling and separation r. At $r\ensuremath{\lesssim}\ensuremath{\xi}$, where $\ensuremath{\xi}$ denotes the healing length of the Bose gas, the interaction is well described by the mean-field contribution. Its form changes as the coupling is increased, approaching a linear function of r at short distances in the regime of strong coupling. The mean-field contribution decays exponentially at arbitrary coupling for $r\ensuremath{\gg}\ensuremath{\xi}$. At such long distances, however, the effect of quantum fluctuations becomes important, giving rise to a long-ranged quantum contribution to the induced interaction. At longest distances it behaves as 1/r3, while at strong coupling we find an intermediate distance regime with a slower decay, 1/r. The quantum contribution in the crossover regime is also calculated. The induced interaction between impurities (i.e., polarons) is attractive and leads to the formation of their bound state, known as bipolaron. We discuss its binding energy.
We present a theoretical study based upon quantum Monte Carlo methods of the Bose polaron in one-dimensional systems with contact interactions. In this instance of the problem of a single impurity immersed in a quantum bath, the medium is a Lieb-Liniger gas of bosons ranging from the weakly interacting to the Tonks-Girardeau regime, whereas the impurity is coupled to the bath via a different contact potential, producing both repulsive and attractive interactions. Both the case of a mobile impurity, having the same mass as the particles in the medium, and the case of a static impurity with infinite mass are considered. We make use of numerical techniques that allow us to calculate the ground-state energy of the impurity, its effective mass, and the contact parameter between the impurity and the bath. These quantities are investigated as a function of the strength of interactions between the impurity and the bath and within the bath. In particular, we find that the effective mass rapidly increases to very large values when the impurity gets strongly coupled to an otherwise weakly repulsive bath. This heavy impurity hardly moves within the medium, thereby realizing the ``self-localization'' regime of the Landau-Pekar polaron. Furthermore, we compare our results with predictions of perturbation theory valid for weak interactions and with exact solutions available when the bosons in the medium behave as impenetrable particles.
One-dimensional Bose gases present an interesting setting to study the physics of Bose polarons, as density fluctuations play an enhanced role due to reduced dimensionality. Theoretical descriptions of this system have predominantly relied on contact pseudopotentials to model the impurity-bath interaction, leading to unphysical results in the strongly coupled limit. In this work, we analytically solve the Gross-Pitaevskii equation, using a square-well potential instead of a zero-range potential, for the ground-state wave function of a static impurity. We compute perturbative corrections arising from infinitesimally slow impurity motion. The polaron energy and effective mass remain finite in the strongly coupled regime, in contrast to the divergent behavior obtained using a contact potential. In this limit, we characterize the polaron properties in terms of the dimensionless ratio w ̃ ≡ w / ξ between the interaction range w of the impurity-bath potential and the coherence length ξ of the Bose gas. The effective mass exhibits a 1 / w ̃ scaling. The energy of the attractive polaron scales as − 1 / w ̃ 3 , whereas the repulsive polaron features subleading corrections to the dark soliton energy at the order w ̃ 3 .
We investigate self-localization of a polaron in a homogeneous Bose-Einstein condensate in one dimension. This effect, where an impurity is trapped by the deformation that it causes in the surrounding Bose gas, has been first predicted by mean-field calculations, but has not been seen in experiments. We study the system in one dimension, where, according to the mean-field approximation, the self-localization effect is particularly robust and present for arbitrarily weak impurity-boson interactions. We address the question whether self-localization is a real effect by developing a variational method which incorporates impurity-boson correlations nonperturbatively and solving the resulting inhomogeneous correlated polaron equations. We find that correlations inhibit self-localization except for very strongly repulsive or attractive impurity-boson interactions. Our prediction for the critical interaction strength for self-localization agrees with a sharp drop of the inverse effective mass found in quantum Monte Carlo simulations of polarons in one dimension. Published by the American Physical Society 2024
A few years ago, flow equations were introduced as a technique for calculating the ground-state energies of cold Bose gases with and without impurities[1,2]. In this paper, we extend this approach to compute observables other than the energy. As an example, we calculate the densities, and phase fluctuations of one-dimensional Bose gases with one and two impurities. For a single mobile impurity, we use flow equations to validate the mean-field results obtained upon the Lee-Low-Pines transformation. We show that the mean-field approximation is accurate for all values of the boson-impurity interaction strength as long as the phase coherence length is much larger than the healing length of the condensate. For two static impurities, we calculate impurity-impurity interactions induced by the Bose gas. We find that leading order perturbation theory fails when boson-impurity interactions are stronger than boson-boson interactions. The mean-field approximation reproduces the flow equation results for all values of the boson-impurity interaction strength as long as boson-boson interactions are weak.