Arxiv:2102.13616V2 [Cond-Mat.Quant-Gas] 30 Jul 2021 That Preserve Stability of the Underlying Problem1

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ürTheoretische Physik, UniversitätHeidelberg, 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ö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.

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öhlich 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öhlich 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 condensate. 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

1 X † † + VBB(q)â 0 aˆ aˆkaˆk0 . (3) Following this strategy we focus on a variational ap- 2V k +q k−q kk0q proach to the ground state of the full Hamiltonian (2). Our ansatz is based on a coherent state. For large parti- ˆ† † Here V is the system volume, dp anda ˆp are the cre- cle number N the results for the coherent state become ation operators of impurity and bosons, respectively, and indistinguishable from those obtained from the product VBB(q) and VIB(q) are the Fourier transforms of the re- form [54] spective interactions in Eq. (1). Next, the ladder operators of bosons are shifted using a simple canonical coher-√ Ψ(x1,..., xN ) = φ(x1) ··· φ(xN ) (6) † † ent state transformation leading toa ˆp → aˆp + δp,0 N0. This is then followed by the crucial Bogoliubov approx- where the condensate wave function φ(x) is normalized R d 2 3 imation: All terms beyond quadratic order bosonic oper- to the condensate particle number d x |φ(x)| = N0 . ators are neglected, leading to the truncated Hamiltonian Both φ(x) and the ground-state energy are found by minimizing the resulting variational Gross-Pitaevskii (GP) X p2 X p2 energy functional Hˆ 0 = dˆ† dˆ + aˆ† aˆ 2m p p 2m p p p I p B R d ¯ 2 √ p0 − d x φ(−i∇)φ N0 N0 X ˆ† ˆ † E[φ] = + V (0) + V (q)d 0 d 0 (â +a ˆ ) 2m V IB V IB k +q k −q q I k0q Z h |∇φ|2 g i + ddx +V (x) |φ(x)|2+ |φ(x)|4−µ |φ(x)|2 . 1 X † † IB + V (q)dˆ dˆ 0 aˆ aˆ (4) 2mred 2 V IB k0+q k k−q k kk0q (7) 2 gBBN0 gBBN0 X † † † Here µ is the chemical potential and we assume weak bo- + + 2â aˆq +a ˆ aˆ +a ˆqaˆ−q , 2V 2V q q −q son repulsion represented in three dimensions by a con- q6=0 tact interaction of strength g = 4πaBB/mB with boson where, for simplicity, we have chosen the example of a scattering length aBB > 0. Note that normal ordering of boson contact interaction of strength gBB. As a result the impurity kinetic term in Eq. (2) contributes to the −1 −1 −1 of the Bogoliubov approximation the bosonic part of the boson kinetic term with reduced mass mred = mB +mI model can be diagonalized using the standard Bogoliubov [21, 30, 40]. rotation The energy functional Eq. (7) exhibits two impor- tant limiting cases: (i) for an infinitely heavy impu- ˆ ∗ † ˆ† ∗ † bp = upaˆp + v−paˆ−p , bp = upaˆp + v−paˆ−p. (5) rity mI → ∞, the normal-ordered kinetic recoil term in the first line of (7) vanishes and the standard GP en- While the Bogoliubov approximation thus allows to ob- ergy functional for bosons in a static external potential tain a simple dispersion relation for bosonic quasipar- is recovered; (ii) for a Bose polaron at rest (p0 = 0), ticles, it, however, captures the effect of repulsion only and a radially symmetric impurity potential VIB(|x|), within the Bogoliubov coefficients up and vp. They give the wave function φ(x) is spherically symmetric and the rise, e.g., to the modified quasiparticle dispersion ωp = p p(p + 2n0gBB) where, notably, the non-deformed, homogenous boson density n0 appears. It turns out that this approximate account of boson repulsion is in- 2 The instability becomes physical for a non-interacting BEC [30, sufficient to self-stabilize the Bose polaron. Indeed, the 38]. 3 neglect of terms beyond quadratic order is responsible For weak boson interactions N ≈ N0. 4 recoil term again vanishes —it only re-appears beyond and an exponentially decaying potential, the product ansatz when boson correlations are included [40, 55]. Vexpon(r) = −V0 exp[−r/R], (10) We find the condensate wave function φ(x) by minimizing the GP functional (7) in the thermodynamic limit both of depth V0 > 0 and range R. subject to the boundary conditions |φ(x → 0)| < ∞, The low-energy scattering properties of the impurity √ and boson are characterized by the impurity-boson scat- |φ(|x| → ∞)| = n0 in terms of the condensate den- tering length a and the effective range r , which deter- sity n0 far away from the impurity. For a radially sym- IB eff mine the leading terms of the effective range expansion of metric impurity potential VIB(|x|) at rest (p0 = 0), the ground-state wave function is spherically symmetric and the momentum-dependent scattering phase shift in three real. In the solution of the GP functional the energy is dimensions. They are found by numerically solving the universally expressed in units of the BEC bulk chemi- Schrödingerequation for the scattering of a boson with the impurity in the center-of-mass system via the poten- cal potential µ = gn0 and the distance from the impurity r = |x| is measured in units of the modified heal- tial VIB(r). Equivalently, one may solve the first-order √ p nonlinear variable phase equation [56], which yields ing length ξ = 1/ 2mredµ = 1/ 8π(mred/mB)aBBn0, which involves the reduced mass and is therefore larger √ a0(r) = 2m V (r)[r − a(r)]2, than the usual bulk healing length ξ = 1/ 2m µ of the red IB 0 B r BEC without impurity. 0 2 reff(r) = −4mredVIB(r)r − 1 We express the energy functional for the three di- a(r) mensional case in terms of the radial function u(r) = r r (r) √ · − 1 + eff . (11) rφ(r)/ n that is scaled in relation to the unperturbed 0 √ 3a(r) r BEC expectation value n0 as Here a(r) and reff(r) obey the boundary conditions E[u] Z ∞ h du2 V (r) a(0) = 0, reff(0) = 0, and account for the phase shift = 4πn dr ξ2 − 1 + IB u(r)2 accumulated by the scattering wave function (for the gen- µ 0 dr µ 0 eralization of the variable phase equation to singular po- (u(r)2 − r2)2 i + . (8) tentials see [57]). Correspondingly, the differential equa- 2r2 tions (11) are integrated from r = 0 ... ∞ which yields the scattering length aIB = a(r → ∞) and effective range Here the boundary conditions for φ(x) translate to r = r (r → ∞). u(0) = 0 and u(r → ∞) = r. Numerically, E[u] is mini- eff eff mized by global optimization on r ∈ [0,L] with L = 10 ξ and an r-grid spacing ∆r = 0.05 ξ much smaller than the III. EFFECTIVE POTENTIAL potential range. We thus obtain the scaling solution of the GP equation First, we present results for the condensate profile (GPE) in units of µ and ξ for a given dimensionless poten- ˜ around the impurity that we find by minimizing the tial shape V (x = r/ξ) = VIB(r)/µ. This scaling solution Gross-Pitaevskii energy functional (7)–(8). For at- is universal for arbitrary values of the condensate density tractive impurity-boson interaction, the wave function n0, boson scattering length aBB and mass ratio mred/mB φ(r) —which, in the comoving frame, directly yields as long as GP theory applies [40, 42], with these parame- the impurity-boson density-density correlation function ters entering only indirectly via µ and ξ. For comparison g(2)(r) = |φ(r)|2/n — is enhanced near the impurity, with experiment, and in order to visualize the effect of IB 0 boson repulsion, the universal solution can be rescaled to as shown in Fig. 2(a). In this figure we have chosen a obtain the specific solution for desired values of the bo- potential VIB(r) (dashed line in Fig. 2(b)) that is suffi- ciently deep to support a two-body bound state at energy son scattering length aBB and the mass ratio mred/mB 2 2/3 −εB < 0. Correspondingly, the potential is character- in density units of energy En = ~ n0 /2mred and length ized by a positive impurity-boson scattering length, in −1/3 n0 . Fig. 2(b) chosen as aIB ≈ 4ξ > 0. We thus realize a scenario as described in the introduction, which is mimicked by a microtrap or a single-site B. Impurity potential Bose-Hubbard model. While at the two-body level it suggests a dynamical instability where the occupation of For the impurity potential VIB(r) we consider two dif- the bound state would grow without bounds, many-body ferent functional forms. This allows us to study the uni- effects come to the rescue. Indeed, due to the repulsion versality of the Bose polaron by analyzing how quasipar- between bosons, each additional boson trying to partic- ticle properties depend on the potential shape and range. ipate in the formation of the Bose polaron is subject to Specifically, we consider an attractive Gaussian potential, the effective potential

2 2 Vgauss(r) = −V0 exp[−(r/R) ], (9) Veﬀ(r) = VIB(r) + g |φ(r)| − µ, (12) 5

4 0 (a) 0 (a) n condensate density / 2

| 3 unperturbed BEC ) 50 r ( 3

| 2 0 = n 100 ) r ( 1 / ) 2 B E ( I g 0 150 scaling solution (b) 0 200 / aIB 1 0

/ 2 (b) ) r ( 3 50 V 2 eff. potential Veff(r) = VIB(r) + g| (r)| 4 100 bare potential V (r) = V exp[ (r/R)2] IB 0 n 5 E

/ 150

0.0 0.5 1.0 1.5 2.0 2.5 3.0 E 8 (m /m )a n1/3 = 1.0 1/3 200 red B BB 0 r n0 1/3 8 (mred/mB)aBBn0 = 0.5 250 FIG. 2. (a) Impurity-boson density-density correlation func- mean field 300 (2) 1/3 tion gIB (r) as determined from the condensate wave function 1/aIB n0 φ(r) as function of the distance r from the impurity (blue (c) 200 solid). For an attractive impurity potential the wave func- 4 3 tion is enhanced near the impurity as compared to that of 150 / 1 √ 0 n 2 an unperturbed BEC φ0(r) = n0 (orange dashed). (b) The f f d e u r

bare impurity potential (orange dashed) of Gaussian shape o l 100 (V /µ = 5.1651, R/ξ = 0.81892) has an effective range r = ξ c 0 eff N 0 and positive scattering length aIB = 4ξ, and correspondingly 2 1 0 1 2 1/a n1/3 admits a bound state. An extra bosonic test particle is, how- 50 IB 0 ever, subject to the effective potential (blue solid) that is weakened by the repulsion from the polaron cloud; while still 0 2.0 1.5 1.0 0.5 0.0 0.5 1.0 1.5 2.0 attractive the effective potential is characterized by an effec- 1/3 1/aIB n0 tive, negative scattering length aIB,eff = −0.1ξ (and renor- malizedr ˜eff = −207ξ) which thus no longer supports a bound state. Hence, additional bosonic quantum fluctuations lead FIG. 3. Bose polaron energy across an impurity-boson Fesh- bach resonance. (a) Scaling solution for the polaron energy only to a weak, additional dressing of the impurity particle. 3 E/µn0ξ as function of the impurity-Bose interaction ξ/aIB for a Gaussian potential of fixed range R = ξ. (b) Polaron 2 2/3 energy E in density units En = ~ n0 /2mred in dependence that results from by the competition of the bare attrac- 1/3 tive potential V (r) (dashed line in Fig. 2(b)) and the on the impurity-Bose interaction 1/aIBn0 for different BEC IB gas parameters. The polaron binding energy is larger for weak repulsion created by the already existing polaron cloud. 1/3 boson repulsion 8π(mred/mB)aBBn0 = 0.5 (orange dashed) In essence, this effect can be understood as arising in 1/3 an effective density-functional theory (DFT) for the BEC compared to stronger repulsion 8π(mred/mB)aBBn0 = 1.0 particles in the presence of the impurity: the bosons (blue solid); at weak coupling it approaches the mean-field result (13) (green dotted). (c) Particle number N within already attracted to the impurity screen the attractive cloud the polaron dressing cloud. Inset: the effective range reff for potential and make it shallower, as shown by the blue 1/3 fixed potential range Rn0 = 1 is smallest on the repulsive line in Fig. 2(b). Moreover, we find that they create side. a small repulsive barrier at intermediate distance. As a consequence, the effective potential no longer admits a bound state and it is correspondingly characterized BEC. The energy is shown in Fig. 3(a) as a function of by a negative effective impurity-boson scattering length the dimensionless impurity-boson interaction ξ/aIB. It aIB,eff ≈ −0.1ξ. Hence, the single-particle excitation is always negative for an attractive impurity potential. spectrum for each additional boson is bounded from be- In Fig. 3(b) we present the energy for a mobile impu- low: the dynamical instability is replaced by Bose po- rity of arbitrary mass in units of En which depends on larons self-stabilized by their dressing cloud. the BEC density n0. Specifically, we show results for two BEC gas parameters, which for equal mass of boson and impurity, i.e., mred/mB = 1/2, correspond to values IV. BOSE POLARON ENERGY AND CONTACT 3 −3 −4 3 n0aBB = (4π) = 5.0 × 10 (n0ξ = 1) at stronger 3 −3 −4 boson repulsion and n0aBB = (8π) = 0.63 × 10 3 The value of the Gross-Pitaevskii functional (7) eval- (n0ξ = 2.8) at weaker boson repulsion. We find that, uated at the ground-state wave function determines the in absolute terms, the polaron binding energy is larger polaron energy E = E[φ] relative to the homogeneous for weaker boson repulsion where the BEC can be more 6 strongly deformed and thus acquire more attractive po- 5000 8 (m /m )a n1/3 = 1.0 tential energy. red B BB 0 8 (m /m )a n1/3 = 0.5 1/3 4000 red B BB 0 At weak impurity-bath attraction 1/aIBn −1 the 0 mean field polaron energy approaches the mean-field value

3000 3 /

2πaIBn0 1/3 1 Emf = = 4πaIBn0 En (13) 0 mred n C 2000 irrespective of aBB. Within mean-field theory the polaron energy diverges to −∞ at unitarity 1/aIB = 0. 1000 Variational approaches based on the trunctated Hamil- tonian Eq. (4) predict that the inclusion of Bogoliubov 0 corrections is not sufficient to heal this instability, but 2.0 1.5 1.0 0.5 0.0 0.5 1.0 1.5 2.0 1/3 rather results in a shift of the instability to the repul- 1/aIB n0 sive side of the Feshbach resonance [30]. In contrast, we find that going beyond the Bogoliubov approximation by FIG. 4. Tan contact C of the Bose polaron for different working with the full Hamiltonian (1) allows the boson boson repulsion (BEC gas parameter) across the impurity- −1/3 repulsion to stabilize the polaron at a finite ground-state boson Feshbach resonance. The contact Cn0 obtained energy that smoothly crosses over from the attractive to from Eq. (15) reaches its maximum on the repulsive side of the repulsive side of the Feshbach resonance. the resonance and increases for weaker boson repulsion. At The deformation of the BEC is also reflected in the weak coupling it approaches the mean-field result (16) (green number of bosons participating in the formation of the dotted). polaron dressing cloud

Z ∞ in an ideal BEC (aBB = 0) [38] 2 2 Ncloud = 4πn0 dr[u(r) − r ]. (14) 2 2 0 Cmf = 16π n0aIB, (16) As shown in Fig. 3(c), for a smaller gas parameter (orange which agrees with the mean-field result prediction. dashed line) the impurity attracts a larger polaron cloud because the bosons are less repulsive. Naturally, this larger dressing cloud corresponds to the larger polaron V. UNIVERSALITY binding energy found in Fig. 3(b). Our results in Fig. 3 are shown for a constant range Finally, we test the notion of universality of the Bose R across the Feshbach resonance as applicable to exper- polaron by studying different shapes and ranges R of the iments where the microscopic range of interactions can impurity potential. Generally, we work in the regime typically not be tuned synchronously with the scattering where the potential range is much larger than the boson length. Thus the effective range reff varies in dependence scattering length, R aBB, where local Gross-Pitaevskii on aIB: The inset of Fig. 3(c) shows the effective range theory has been shown to be applicable [35]. To test uni- that is obtained from the scattering phase shift using versality in this regime, we specifically compare the pre- 1/3 Eq. (11). At constant potential range Rn0 = 1, the dictions following from the Gaussian potential Vgauss(r) effective range reff reaches a minimum on the repulsive in Eq. (9), and the exponential potential Vexpon(r) (10) side of the resonance (aIB > 0) and grows towards both for various ranges R. ± weak-coupling limits aIB → 0 . We find that when tuning the depth of the potentials to The variation of the polaron energy with aIB defines obtain equal scattering length aIB = −∞ at equal range the impurity contact parameter [38, 40, 42, 51, 58] which R, the polaron energies are very different. Instead, if R (2) is tuned to yield the same effective range reff for both characterizes the impurity-boson correlations gIB (r) at short distances outside the impurity potential: potentials, as shown in Fig. 5(a), remarkable agreement between the polaron energies is found. Indeed, as shown 8πm ∂E ∂(E/E ) in Fig. 5(b), for both potentials, the polaron energy ap- C = red = 4πn1/3 n . (15) 2 0 1/3 ∂(−1/aIB) proximately follows a linear scaling law for reff & 0.2ξ, ~ ∂(−1/aIBn0 ) E(1/aIB = 0) 1/3 The contact is shown in Fig. 4: we find that it reaches = −5.2(2) − 9.0(1) reffn0 . (17) a maximal value on the repulsive side of the resonance. En Similar to the energy, we find that the contact is larger This shows that in the Bose polaron problem the for smaller boson repulsion aBB (blue solid line) where momentum-dependent scattering phase shift is probed in (2) the BEC is more strongly deformed and thus gIB (r) is the regime where the effective range expansion is valid. enhanced (see Fig. 2). At weak attractive interaction the Remarkably, the polaron energy remains sizeable even contact approaches the ground-state value of an impurity for short effective range. Indeed, our result complements 7

(a) 3 1/3 Gauss side, at 1/aIBn = 1: universality can at best hold in R = 0.69 reff 0 2 the vicinity of the ﬁrst bound state, and only for a limited 3 / 1 0

n interval of reff values. Furthermore, since the Efimov ef- expon R R = 0.28 r 1 eff fect modifies the polaron energy spectrum [2, 3, 49], universality can also depend on the three-body parameter 0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 [50, 59–61]. 1/3 reff n0 0 (b) Gauss 1/3 10 E /En 5.0 9.2 reffn0 expon 1/3 VI. DISCUSSION E /En 5.4 8.9 reffn n 20 0 E /

E 30 Gaussian potential The variational principle provides a powerful tool to 40 exponential potential compute both ground-state and dynamical properties of 50 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 quantum many-body systems. However, in order to make r n1/3 eff 0 full use of its predictive power it is essential to under- (c) 0 (d) 0 1/aIB = 1 1/aIB = + 1 stand the limitations of approximations applied to the 2 50 Hamiltonian to be analyzed. In this regard the strong-

n 4 n 100 E E

/ / coupling Bose polaron is a case in point. The full Bose E 6 E 150 polaron Hamiltonian (1) —and equivalently Eq. (3)— is 8 200 bounded from below for repulsive boson interaction, and 10 250 0 1 2 3 4 0 1 2 3 4 hence variational wave functions give rigorous bounds on 1/3 1/3 reff n0 reff n0 the ground-state energy. For finite Bose repulsion, the polaron energy remains finite for any Bose-impurity scat- FIG. 5. Universality of the Bose polaron. (a) In order to tering length, including resonant interactions, and the yield the same effective range reff, the Gaussian and exponen- ground state represents a strong-coupling Bose polaron tial potentials need to be tuned to different potential ranges that is self-stabilized by its own dressing cloud. R; here shown for fixed unitary scattering length 1/aIB = 0. In contrast, when applying the Bogoliubov approxima- (b) Polaron energy E as function of the effective range r at eff tion to the full model (1) by truncating terms of higher- unitarity 1/aIB = 0 for two different potential shapes. The 1/3 than quadratic order in the boson operators, the result- boson repulsion is set to 8π(mred/mB)aBBn0 = 1. The po- ˆ 0 laron energy coincides for both potential shapes and increases ing, truncated Hamiltonian H in Eq. (4) becomes un- bound from below. Crucially, this results in an instabil- linearly with reff. For short range it approaches the value for reff = 2aBB from [42] (star). (c) The polaron energy at ity of the Bose polaron problem that is solely an artefact 1/3 of this approximation. negative scattering length 1/aIBn0 = −1 is universal for r n1/3 1 for the parameters chosen. (d) Polaron energy at Quite remarkably, the instability of the truncated eff 0 & 0 1/3 Hamiltonian Hˆ becomes, however, only evident when positive scattering length 1/aIBn0 = 1. For the Gaussian potential there is no unique polaron energy for intermediate considering wavefunctions that account for more than values of reff (see text). two phonon excitations from the homogenous BEC. For instance, for a Chevy-type wavefunction [25, 27, 28] that itself is truncated at the single excitation level, the terms beyond quadratic order in the repulsive interactions have a recent study which found a power-law scaling of the vanishing expectation value. Thus, incidentally, the unitary polaron energy at ranges much shorter than the 1/3 Chevy ansatz yields the same prediction when applied healing length, E/En ∼ (reff/ξ) for aBB . reff ξ to both the full and the truncated model. Thus due to in the case of a square-well potential [42]. In the region its tremendous simplicity the Chevy ansatz becomes im- for reff = 2aBB this approach yields a polaron energy of mune to the instability of approximate Hamiltonian Hˆ 0. E ≈ −5 En as indicated by the star in Fig. 5(a), which However, while being a well-defined approach, the sim- is consistent with the linear extrapolation of our result. ple Chevy ansatz misses the fact that at the weak Boson These findings can be further extended to even shorter repulsion (as typically present in cold gases) the polaron ranges reff aBB by using a nonlocal generalization of cloud—even within the full model Eq. (1)—can contain Gross-Pitaevskii theory [40]. an exceedingly large number of bosons (Fig. 3(b)). Such Also on the attractive side of the resonance, see a large local deformation of the BEC is naturally cap- 1/3 Fig. 5(c) for 1/aIBn0 = −1, we find a universal effec- tured by the coherent state (6). Crucially, while account- tive range dependence of the polaron energy, albeit not ing for an arbitrary number of boson excitations, when a linear one. Note that to cover the domain of effective applied to the full Hamiltonian (1), it still leads to a ranges shown in Fig. 5(d), for our parameters the Gaus- bounded energy functional (7). Its solution shows that sian potential, unlike the exponential, must be tuned over the smaller the boson repulsion and the wider the impu- such depths that it supports more than one bound state. rity potential, the larger the polaron cloud becomes. The Consequently, there is no longer a unique mapping from coherent state approach is complementary to the Chevy reff to E. This is illustrated in Fig. 5(d) for the repulsive ansatz including its extensions to multi-boson excitations 8

[3, 49, 62], and it becomes particularly accurate for soft boson correlations, three-body and higher-order corre- impurity potentials where it is justified to ignore bosonic lations give rise to the Efimov effect and three-body re- correlations. Remarkably, the case of extremely soft po- combination. The Efimov effect can occur either between tentials is realized with Rydberg excitation immersed in one impurity and two bosons [3, 49, 62] or between two BECs. In this case it was predicted that up to hundreds impurities and one boson [2, 6]. These few-body effects of atoms can be bound to the single impurity leading to can be captured by Gaussian variational wave functions the creation of Rydberg polarons [63]. Since for Ryd- [55]. Alternatively, extensions of the Chevy ansatz to berg impurities the range of interactions R dramatically two or more independent bosonic excitations [3] can be exceeds the interparticle distance, our local GP theory employed. The latter approach was applied in the anal- applies and provides a so far missing explanation as to ysis of the truncated model (4) and universal scaling de- why the experimental observation of Rydberg polarons pending on the three-body parameter was found [49, 62]. [29] is described exceptionally well by a coherent state It remains an interesting open question how universal approach [26]. three-body physics carries over to the many-body case Recently also first steps to the understanding of the in a dense bosonic medium when the full Hamiltonian complementary, intermediate regime of short-range im- Eq. (1) is considered. purity potentials —with yet large dressing clouds— has Finally, in ultracold atomic gases the boson repulsion been achieved by using an extension to nonlocal Gross- originates from an attractive van-der-Waals potential be- Pitaevskii theory [40]. In conjunction with our present tween atoms. This results in the existence of deeply result, these combined new approaches resolve a funda- bound molecular states into which atoms can decay in mental shortcoming of the Bogoliubov approximation: three-body recombination. These deeply bound states while the formulation of the interacting Bose gas in terms are neither accounted for in variational approaches nor of Bogoliubov quasiparticles is exact, the additional ap- in quantum Monte Carlo [43]. Up to now most experi- proximation to neglect the residual interaction between ments probe the metastable Bose polaron state of matter phonons is not. In particular, the quadratic Bogoli- on transient time scales where the deeply bound states — ubov mean-field Hamiltonian is unbounded from below arising both from the fundamentally attractive Bose-Bose for the strong-coupling Bose polaron, most obviously in and Bose-impurity potentials— can be ignored, and thus the regime where a two-body bound state appears on the the stable variational approach discussed in our work is repulsive side of the resonance [30, 31]. well applicable. However, as one starts to explore longer time scales or the build-up of more complex correlated As discussed above, the Chevy-type ansatz applied to states of impurities, understanding the impact of the fun- the truncated the Bogoliubov Hamiltonian [28] still yields damental dissipative nature arising from deeply bound a finite polaron energy since it is of such low order in bo- states becomes essential and requires the development of son excitations that it is not sensitive to the truncated new theoretical approaches to quantum impurity prob- part of the Hamiltonian. However, when the coherent lems. state ansatz or higher-order excitation extensions of the Chevy ansatz are applied to the truncated Hamiltonian (4) they can lead to divergencies in the ground-state en- ACKNOWLEDGMENTS ergy. In the case of the coherent state ansatz [30, 36] the divergence is due to the large occupation of excita- This work is supported by the Deutsche Forschungs- tions in the vicinity of the impurity not counteracted by gemeinschaft (DFG, German Research Foundation), boson repulsion. Instead the local extension to the trun- project-ID 273811115 (SFB1225 ISOQUANT) and under cated Bogoliubov approach studied in the present work Germany’s Excellence Strategy EXC2181/1-390900948 (see also [42, 50, 51]) as well as its nonlocal counterpart (the Heidelberg STRUCTURES Excellence Cluster). [40] provides a stable starting point for strong-coupling R. S. is supported by the Deutsche Forschungsgemein- Bose polaron dynamics, and we showed how it can find schaft (DFG, German Research Foundation) under Ger- an effective description in terms of a renormalization of many’s Excellence Strategy – EXC-2111 – 390814868 the impurity-boson potential (Fig. 2). (Excellence Cluster ‘Munich Center for Quantum Science Beyond our treatment of the two-particle impurity- and Technology’).

[1] J. T. Devreese and A. S. Alexandrov, Fröhlich polaron Rev. Lett. 115, 125302 (2015). and bipolaron: recent developments, Rep. Prog. Phys. [4] T. Rentrop, A. Trautmann, F. A. Olivares, F. Jendrze- 72, 066501 (2009). jewski, A. Komnik, and M. K. Oberthaler, Observation of [2] N. T. Zinner, Efimov states of heavy impurities in a Bose- the phononic Lamb shift with a synthetic vacuum, Phys. Einstein condensate, Europhys. Lett. 101, 60009 (2013). Rev. X 6, 041041 (2016). [3] J. Levinsen, M. M. Parish, and G. M. Bruun, Impurity in [5] R. Schmidt and M. Lemeshko, Rotation of quantum im- a Bose-Einstein condensate and the Efimov effect, Phys. purities in the presence of a many-body environment, 9

Phys. Rev. Lett. 114, 203001 (2015). T. Pohl, Critical slowdown of non-equilibrium polaron [6] P. Naidon, Two Impurities in a Bose-Einstein Conden- dynamics, New J. Phys. 21, 043014 (2019). sate: From Yukawa to Efimov Attracted Polarons, J. [25] S. P. Rath and R. Schmidt, Field-theoretical study of the Phys. Soc. Japan 87, 043002 (2018). Bose polaron, Phys. Rev. A 88, 053632 (2013). [7] A. Camacho-Guardian, L. A. PeñaArdila, T. Pohl, and [26] R. Schmidt, J. D. Whalen, R. Ding, F. Camargo, G. M. Bruun, Bipolarons in a Bose-Einstein condensate, G. Woehl, S. Yoshida, J. Burgdörfer,F. B. Dunning, Phys. Rev. Lett. 121, 013401 (2018). E. Demler, H. R. Sadeghpour, and T. C. Killian, Theory [8] Z. Z. Yan, Y. Ni, C. Robens, and M. W. Zwierlein, of excitation of Rydberg polarons in an atomic quantum Bose polarons near quantum criticality, Science 368, 190 gas, Phys. Rev. A 97, 022707 (2018). (2020). [27] F. Chevy, Universal phase diagram of a strongly interact- [9] J. Catani, G. Lamporesi, D. Naik, M. Gring, M. Inguscio, ing Fermi gas with unbalanced spin populations, Phys. F. Minardi, A. Kantian, and T. Giamarchi, Quantum Rev. A 74, 063628 (2006). dynamics of impurities in a one-dimensional Bose gas, [28] W. Li and S. Das Sarma, Variational study of polarons Phys. Rev. A 85, 023623 (2012). in Bose-Einstein condensates, Phys. Rev. A 90, 013618 [10] R. Scelle, T. Rentrop, A. Trautmann, T. Schuster, and (2014). M. K. Oberthaler, Motional coherence of fermions im- [29] F. Camargo, R. Schmidt, J. D. Whalen, R. Ding, mersed in a Bose gas, Phys. Rev. Lett. 111, 070401 G. Woehl, S. Yoshida, J. Burgdörfer,F. B. Dunning, (2013). H. R. Sadeghpour, E. Demler, and T. C. Killian, Cre- [11] M.-G. Hu, M. J. Van de Graaff, D. Kedar, J. P. Corson, ation of Rydberg Polarons in a Bose Gas, Phys. Rev. E. A. Cornell, and D. S. Jin, Bose polarons in the strongly Lett. 120, 083401 (2018). interacting regime, Phys. Rev. Lett. 117, 055301 (2016). [30] Y. E. Shchadilova, R. Schmidt, F. Grusdt, and E. Dem- [12] N. B. Jørgensen, L. Wacker, K. T. Skalmstang, M. M. ler, Quantum dynamics of ultracold Bose polarons, Phys. Parish, J. Levinsen, R. S. Christensen, G. M. Bruun, Rev. Lett. 117, 113002 (2016). and J. J. Arlt, Observation of attractive and repulsive [31] F. Grusdt, R. Schmidt, Y. E. Shchadilova, and E. Demler, polarons in a Bose-Einstein condensate, Phys. Rev. Lett. Strong-coupling Bose polarons in a Bose-Einstein con- 117, 055302 (2016). densate, Phys. Rev. A 96, 013607 (2017). [13] M. G. Skou, T. G. Skov, N. B. Jørgensen, K. K. Nielsen, [32] P. Massignan, C. J. Pethick, and H. Smith, Static proper- A. Camacho-Guardian, T. Pohl, G. M. Bruun, and J. J. ties of positive ions in atomic bose-einstein condensates, Arlt, Non-equilibrium quantum dynamics and formation Phys. Rev. A 71, 023606 (2005). of the Bose polaron, Nature Phys. 17, 731 (2021). [33] G. E. Astrakharchik and B. A. Malomed, Quantum ver- [14] H. Fröhlich, Electrons in lattice fields, Advances in sus mean-field collapse in a many-body system, Phys. Physics 3, 325 (1954). Rev. A 92, 043632 (2015). [15] M. Girardeau, Motion of an Impurity Particle in a Boson [34] Y. Ashida, R. Schmidt, L. Tarruell, and E. Demler, Superfluid, Phys. Fluids 4, 279 (1961). Many-body interferometry of magnetic polaron dynam- [16] G. E. Astrakharchik and L. P. Pitaevskii, Motion of ics, Phys. Rev. B 97, 060302 (2018). a heavy impurity through a Bose-Einstein condensate, [35] K. Chen, N. V. Prokof’ev, and B. V. Svistunov, Trapping Phys. Rev. A 70, 013608 (2004). collapse: Infinite number of repulsive bosons trapped by [17] F. M. Cucchietti and E. Timmermans, Strong-coupling a generic short-range potential, Phys. Rev. A 98, 041602 polarons in dilute gas Bose-Einstein condensates, Phys. (2018). Rev. Lett. 96, 210401 (2006). [36] M. Drescher, M. Salmhofer, and T. Enss, Real-space [18] M. Bruderer, W. Bao, and D. Jaksch, Self-trapping of dynamics of attractive and repulsive polarons in Bose- impurities in Bose-Einstein condensates: Strong attrac- Einstein condensates, Phys. Rev. A 99, 023601 (2019). tive and repulsive coupling, Europhys. Lett. 82, 30004 [37] O. Hryhorchak, G. Panochko, and V. Pastukhov, Im- (2008). purity in a three-dimensional unitary Bose gas, Phys. [19] J. Tempere, W. Casteels, M. K. Oberthaler, S. Knoop, Lett. A 384, 126934 (2020). E. Timmermans, and J. T. Devreese, Feynman path- [38] M. Drescher, M. Salmhofer, and T. Enss, Quench Dy- integral treatment of the BEC-impurity polaron, Phys. namics of the Ideal Bose Polaron at Zero and Nonzero Rev. B 80, 184504 (2009). Temperatures, Phys. Rev. A 103, 033317 (2021). [20] W. Casteels, T. Van Cauteren, J. Tempere, and J. T. De- [39] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Fesh- vreese, Strong coupling treatment of the polaronic system bach resonances in ultracold gases, Rev. Mod. Phys. 82, consisting of an impurity in a condensate, Laser Phys. 21, 1225 (2010). 1480 (2011). [40] M. Drescher, M. Salmhofer, and T. Enss, Theory of a [21] A. Shashi, F. Grusdt, D. A. Abanin, and E. Demler, resonantly interacting impurity in a Bose-Einstein con- Radio-frequency spectroscopy of polarons in ultracold densate, Phys. Rev. Research 2, 032011(R) (2020). Bose gases, Phys. Rev. A 89, 053617 (2014). [41] I. Bloch, J. Dalibard, and W. Zwerger, Many-body [22] F. Grusdt, Y. E. Shchadilova, A. N. Rubtsov, and physics with ultracold gases, Rev. Mod. Phys. 80, 885 E. Demler, Renormalization group approach to the (2008). Fröhlich polaron model: application to impurity-BEC [42] P. Massignan, N. Yegovtsev, and V. Gurarie, Universal problem, Sci. Rep. 5, 12124 (2015). aspects of a strongly interacting impurity in a dilute Bose [23] J. Vlietinck, W. Casteels, K. Van Houcke, J. Tempere, condensate, Phys. Rev. Lett. 126, 123403 (2021). J. Ryckebusch, and J. T. Devreese, Diagrammatic Monte [43] L. A. PeñaArdila and S. Giorgini, Impurity in a Bose- Carlo study of the acoustic and the Bose–Einstein con- Einstein condensate: study of the attractive and repul- densate polaron, New J. Phys. 17, 033023 (2015). sive branch using quantum Monte Carlo methods, Phys. [24] K. K. Nielsen, L. A. PeñaArdila, G. M. Bruun, and Rev. A 92, 033612 (2015). 10

[44] A. S. Dehkharghani, A. G. Volosniev, and N. T. Zinner, Einstein condensate, Comm. Phys. 4, 1 (2021). Quantum impurity in a one-dimensional trapped bose [53] T. D. Lee, F. E. Low, and D. Pines, The motion of slow gas, Phys. Rev. A 92, 031601 (2015). electrons in a polar crystal, Phys. Rev. 90, 297 (1953). [45] S. I. Mistakidis, A. G. Volosniev, N. T. Zinner, and [54] L. P. Pitaevskii and S. Stringari, Bose-Einstein Conden- P. Schmelcher, Effective approach to impurity dynam- sation, Vol. 116 (Oxford University Press, 2003). ics in one-dimensional trapped bose gases, Phys. Rev. A [55] A. Christianen, I. Cirac, and R. Schmidt, in prep. (2021). 100, 013619 (2019). [56] F. Calogero, Variable Phase Approach to Potential Scat- [46] S. I. Mistakidis, F. Grusdt, G. M. Koutentakis, and tering (Academic Press, New York, 1967). P. Schmelcher, Dissipative correlated dynamics of a mov- [57] T. Enss, B. Tran, M. Rautenberg, M. Gerken, E. Lippi, ing impurity immersed in a Bose-Einstein condensate, M. Drescher, B. Zhu, M. Weidemüller,and M. Salmhofer, New J. Phys. 21, 103026 (2019). Scattering of two heavy Fermi polarons: resonances and [47] J. Jager, R. Barnett, M. Will, and M. Fleischhauer, quasibound states, Phys. Rev. A 102, 063321 (2020). Strong-coupling bose polarons in one dimension: Con- [58] S. Tan, Energetics of a strongly correlated Fermi gas, densate deformation and modified bogoliubov phonons, Ann. Phys. (N.Y.) 323, 2952 (2008). Phys. Rev. Res. 2, 033142 (2020). [59] Y. Wang, J. Wang, J. P. D’Incao, and C. H. Greene, [48] S. I. Mistakidis, G. M. Koutentakis, G. C. Katsimiga, Universal three-body parameter in heteronuclear atomic T. Busch, and P. Schmelcher, Many-body quantum dy- systems, Phys. Rev. Lett. 109, 243201 (2012). namics and induced correlations of Bose polarons, New [60] R. Schmidt, S. P. Rath, and W. Zwerger, Efimov physics J. Phys. 22, 043007 (2020). beyond universality, The European Physical Journal B [49] S. M. Yoshida, S. Endo, J. Levinsen, and M. M. Parish, 85, 386 (2012). Universality of an impurity in a Bose-Einstein conden- [61] C. Langmack, R. Schmidt, and W. Zwerger, Efimov sate, Phys. Rev. X 8, 011024 (2018). states near a Feshbach resonance and the limits of van der [50] B. Field, J. Levinsen, and M. M. Parish, Fate of the Bose Waals universality at finite background scattering length, polaron at finite temperature, Phys. Rev. A 101, 013623 Phys. Rev. A 97, 033623 (2018). (2020). [62] Z.-Y. Shi, S. M. Yoshida, M. M. Parish, and J. Levinsen, [51] N.-E. Guenther, R. Schmidt, G. M. Bruun, V. Gu- Impurity-induced multibody resonances in a Bose gas, rarie, and P. Massignan, Mobile impurity in a Bose- Phys. Rev. Lett. 121, 243401 (2018). Einstein condensate and the orthogonality catastrophe, [63] R. Schmidt, H. R. Sadeghpour, and E. Demler, Meso- Phys. Rev. A 103, 013317 (2021). scopic Rydberg Impurity in an Atomic Quantum Gas, [52] G. E. Astrakharchik, L. A. PeñaArdila, R. Schmidt, Phys. Rev. Lett. 116, 105302 (2016). K. Jachymski, and A. Negretti, Ionic polaron in a Bose-