Arxiv:2102.13616V2 [Cond-Mat.Quant-Gas] 30 Jul 2021 That Preserve Stability of the Underlying Problem1
Total Page:16
File Type:pdf, Size:1020Kb
Self-stabilized Bose polarons Richard Schmidt1, 2 and Tilman Enss3 1Max-Planck-Institute of Quantum Optics, Hans-Kopfermann-Straße 1, 85748 Garching, Germany 2Munich Center for Quantum Science and Technology, Schellingstraße 4, 80799 Munich, Germany 3Institut f¨urTheoretische Physik, Universit¨atHeidelberg, 69120 Heidelberg, Germany (Dated: August 2, 2021) 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¨ohlich 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. I. INTRODUCTION At weak interactions between the impurity and the BEC, the impurity deforms the BEC only slightly and Impurities in a Bose-Einstein condensate (BEC) ex- physics is well described within the Fr¨ohlich model in hibit a multitude of fundamental physical phenomena: terms of long-wavelength phonon excitations [1, 14{24]. the formation of quasiparticles [1], Efimov bound states For strong local attraction, instead, the Fr¨ohlich model is [2, 3], synthetic Lamb shifts [4, 5], Casimir interactions incomplete and needs to be amended by quadratic terms induced by a fluctuating medium [6, 7], and quantum that absorb and re-emit phonons [25]. Importantly, these criticality [8]. Current experiments with ultracold atomic term are also required to correctly capture the formation gases are investigating several of these effects reaching of bound states between the impurity and bath atoms far into the strong-coupling regime [4, 8{13]. For un- which is a crucial ingredient for the physics of polorans derstanding experimental observations it is thus vital at strong coupling [26]. Variational wave functions based to develop a theoretical model that applies to impurity on a single phonon excitation [27, 28] are able to describe systems at strong coupling and that can address both the single occupation of such a bound state. However, ground state and non-equilibrium phenomena. the bosons that make up the BEC tend to bunch, and at strong coupling it is energetically favorable to occupy the bound state multiple times leading to a gain of several energy times the binding energy. This process, recently observed with Rydberg atoms immersed in a BEC [29], is described Bose-Einstein condensate by a coherent state variational ansatz [25, 26, 30, 31] that allows for an arbitrary number of excitations and strong local deformations of the condensate [32{38]. x Generally, the application of the variational principle effective potential for the determination of the ground state is only viable bound state if the Hamiltonian is bounded from below. Under this condition a stable solution can be found, and theoreti- impurity potential cal approaches should rely only on such approximations arXiv:2102.13616v2 [cond-mat.quant-gas] 30 Jul 2021 that preserve stability of the underlying problem1. For FIG. 1. Illustration of the self-stabilized Bose polaron. The the strongly interacting Bose polaron, the resulting the- strong-coupling Bose polaron mimics a microtrap (red dashed oretical challenge can be understood in the simple toy line) with bound state (red bar) within the Bose-Einstein con- densate. If several bosons (blue dots) occupy the bound state, boson repulsion results in a shallower effective potential seen by the remaining bosons (blue solid line) that no longer ad- 1 In this work we restrict ourselves to models of cold dilute gases mits a bound state. The dressing cloud itself thus stabilizes that disregard deeply bound states, transitions to solid or liquid the Bose polaron. phases, as well as large cluster formation. 2 model illustrated in Fig. 1. Here a static attractive im- and finite-range attractive impurity potentials employing purity is represented by a local potential well of finite a coherent variational state that allows for large dress- range around the impurity. This well acts as a microtrap ing clouds and strong local deformations of the BEC. In within the BEC. At strong coupling beyond a scattering Sec. II we introduce the stable Bose polaron model and resonance [39], the potential well (red dashed line) is deep discuss its solution within Gross-Pitaevskii theory. In enough to admit a bound state with energy "i = −"B < 0 Section III we minimize the resulting GP energy func- (red bar). In absence of boson repulsion, in the many- tional and obtain the condensate wave function around body ground state all bosons would occupy the bound the impurity. We find that even if the bare impurity state and the ground-state energy E0 = −N"B is indeed potential admits a bound state, the emerging effective unbounded from below in the thermodynamic limit: the potential does not, thus providing a simple mechanism whole BEC is collapsed onto the impurity. for the self-stabilization of Bose polarons. Section IV Naturally, a local boson repulsion counteracts this pro- presents our results for the polaron energy, the size of cess by balancing the impurity attraction [40] and thus the polaron dressing cloud and the Tan contact across providing a lower bound to the ground-state energy. the resonance. The universality of the polaron energy is This mechanism can be understood in terms of a single- discussed in Sec. V, and we show that the energy at uni- site Bose-Hubbard model [41, 42] with local repulsion tarity depends universally on the effective range, as long U as Efimov states can be neglected. Finally, in Sec. VI we 2 ni(ni − 1) on the impurity site ri, that competes with a local attractive potential energy −"Bni for occupation compare variational approaches to the Bose polaron and number ni. For U > 0 the ground state has a finite occu- discuss which Hamiltonians and energy functionals can 2"B provide rigorous bounds on the ground-state energy. pation ni ' U , which is nonperturbative in the strength of the interaction: the size of this polaron dressing cloud grows for weaker repulsion. It is crucial to adequately capture this repulsive effect in the theoretical description II. MODEL of Bose polarons. Previous variational approaches have included the bo- We consider a single impurity particle immersed in an son repulsion only at the level of the Bogoliubov approx- interacting Bose gas. The combined system is described imation. Here the interaction between the Bogoliubov by the Hamiltonian quasiparticles is neglected, and thus bosons in the bound p^2 X state fail to generate the compensating pressure required H = + V (x^ − x^) to ensure the stability of the ground state. As a conse- 2m IB B;i I i quence this approach falsely predicts instead a dynami- 2 X p^B;i X cal instability in presence of boson repulsion toward infi- + + V (x^ − x^ ): (1) 2m BB B;i B;j nite occupation of the bound state in the strong coupling i B i<j regime [30, 31, 36]. This shows the crucial importance of including the boson repulsion beyond the Bogoliubov Here, p^ and x^ denote the momentum and position of approximation. the impurity of mass mI, and p^B;i and x^B;i character- In this work, we present a stable variational approach ize the bosons i = 1;:::;N of mass mB. The boson to the Bose polaron problem at strong coupling. Our interaction VBB(x) is assumed to be repulsive and the approach applies to arbitrary dimension and impurity- impurity-boson interaction VIB(x) is attractive. boson scattering lengths and it provides a basis for the The coupled system of impurity and bosons is conve- study of dynamical properties of Bose polarons. In niently analyzed in the reference frame comoving with this way our work complements and extends previous the impurity. This is achieved by a canonical trans- approaches using quantum Monte Carlo to determine formation introduced by Lee, Low, and Pines [53] and ground-state properties [43], studies of the role of bo- elaborated on by Girardeau [15]. The transformation P son repulsion in one-dimensional systems [44{48], or a S = exp(ix^ · i p^B;i) leads to the LLP Hamiltonian recently developed nonlocal Gross-Pitaevskii (GP) the- P 2 ory for nonequilibrium dynamics [40]. Moreover, previ- −1 (p0 − i p^B;i) X HLLP = S HS = + VIB(x^B;i) ous works raised the question to which extent polaron 2mI properties are universal in the short-range limit [49, 50]. i 2 In the following we study this question across the Fesh- X p^B;i X + + V (x^ − x^ ); (2) bach resonance, and in particular in the regime where a 2m BB B;i B;j i B i<j bound state is supported by the impurity-bath potential, extending previous work that considered how the prop- in which the impurity operators have been eliminated, erties of the Bose polaron depend on the range of the and p0 denotes the conserved total momentum of the interactions, both for the finite-range boson repulsion in system.