PHYSICAL REVIEW RESEARCH 2, 023405 (2020)
Strong coupling Bose polarons in a two-dimensional gas
L. A. Peña Ardila ,1,2,* G. E. Astrakharchik ,3 and S. Giorgini4 1Institut für Theoretische Physik, Leibniz Universität, 30167 Hannover, Germany 2Institut for Fysik og Astronomi, Aarhus Universitet, 8000 Aarhus C, Denmark 3Departament de Física, Universitat Politècnica de Catalunya, Campus Nord B4-B5, E-08034 Barcelona, Spain 4Dipartimento di Fisica, Università di Trento, CNR-INO BEC Center, I-38123 Povo, Trento, Italy
(Received 8 July 2019; revised manuscript received 16 January 2020; accepted 11 June 2020; published 26 June 2020)
We study the properties of Bose polarons in two dimensions using quantum Monte Carlo techniques. Results for the binding energy, the effective mass, and the quasiparticle residue are reported for a typical strength of interactions in the gas and for a wide range of impurity-gas coupling strengths. A lower and an upper branch of the quasiparticle exist. The lower branch corresponds to an attractive polaron and spans from the regime of weak coupling where the impurity acts as a small density perturbation of the surrounding medium to deep bound states which involve many particles from the bath and extend as far as the healing length. The upper branch corresponds to an excited state where due to repulsion a low-density bubble forms around the impurity but might be unstable against decay into many-body bound states. Interaction effects strongly affect the quasiparticle properties of the polaron. In particular, in the strongly correlated regime, the impurity features a vanishing quasiparticle residue, signaling the transition from an almost free quasiparticle to a bound state involving many atoms from the bath.
DOI: 10.1103/PhysRevResearch.2.023405
I. INTRODUCTION servation of Bose polarons in ultracold gases has triggered an intense research activity aiming at describing the crossover Impurities embedded in a quantum many-body environ- from weak- to strong-coupling regimes. In the former case, ment can lead to the formation of quasiparticles coined po- the so-called Bogoliubov-Fröhlich Hamiltonian describes ac- larons. The concept was first introduced by Landau and Pekar curately the ground state properties of the polaron [21–30]. in the solid-state context to describe an electron coupled to However, quantum fluctuations become relevant as interac- an ionic crystal [1]. Polarons are fundamental ingredients in tions are increased, making the description in terms of the many different transport phenomena across condensed-matter Fröhlich paradigm inadequate. By using the Gross-Pitaevskii physics. Electronic transport in polar crystals or semicon- equation, strongly interacting Bose polarons were predicted to ductors [2] as well as charge and spin transport in organic manifest exotic phenomena, such as self-localization [31–35], materials [3,4] can be understood in terms of polarons. Pair- but without experimental evidence so far. Recently, properties ing between polarons is relevant in the physics of high- of these strongly coupled impurities have also been addressed temperature superconductors [5], and polarons are candidates by techniques, such as theT matrix, diagrammatic, and for electronic transport in DNA and proteins [6]. Furthermore, variational approaches which go beyond the single-phonon polarons are used as probes of quantum many-body systems. excitation scheme of the Fröhlich model [36–42]. These For example, the low-energy excitations in a strongly corre- studies predict exotic out of equilibrium dynamics, nontriv- lated superfluid, such as 4He can be probed by 3He impurity ial quasiparticle splitting due to finite-temperature effects atoms [7]. [40,43–48], as well as important few-body effects [49,50]. Unprecedented control and versatility of ultracold gases [8] Furthermore, the regime of strong coupling should also fea- made it possible to experimentally observe dressed impurities ture the interchange of Bogoliubov modes between polarons named Fermi and Bose polarons depending on whether they via polaron-polaron interactions [51,52]. In the context of interact, respectively, with a degenerate Fermi gas [9–13] theoretical techniques suitable to investigate this latter regime, or a Bose-Einstein condensate (BEC). For Bose polarons, the quantum Monte Carlo (QMC) method is based on a micro- experiments have been carried out in three-dimensional (3D) scopic Hamiltonian and provides exact (within controllable [14–18] and one-dimensional (1D) [19,20] geometries. Ob- statistical errors) ground-state properties of the polaron for arbitrary coupling strengths [39,41,53,54]. Physically, the two-dimensional (2D) geometry is ap- *[email protected] pealing since the role of quantum fluctuations is enhanced whereas off-diagonal long-range order, responsible for BEC Published by the American Physical Society under the terms of the phenomena, still exists in the ground state. Polarons in 2D Creative Commons Attribution 4.0 International license. Further geometries have been extensively investigated in the context distribution of this work must maintain attribution to the author(s) of Fermi polarons [10,55–58] and exciton impurities coupled and the published article’s title, journal citation, and DOI. to semiconductors [59]. Bose polarons have been investigated
2643-1564/2020/2(2)/023405(8) 023405-1 Published by the American Physical Society ARDILA, ASTRAKHARCHIK, AND GIORGINI PHYSICAL REVIEW RESEARCH 2, 023405 (2020)
through the two-body potential VB, which depends on the dis- tance rij =|ri − r j| between a pair of bosons. Furthermore, 2∇2 − h¯ α is the kinetic energy of the impurity denoted by the 2mI coordinate vector rα, and VI is the boson-impurity potential depending on the distance riα =|rα − ri| between the impu- rity and the ith bath particle. Both interaction potentials VB and VI are short ranged and are parametrized by the scattering lengths aB and a, respectively. Within Bogoliubov theory, the Hamiltonian (1) can be written in second quantization as the sum of two terms H = H0 + Hint, where p2 H = + E + α†α (2) 0 2m B k k k I k is the unperturbed Hamiltonian of a free impurity moving with momentum p and a static host gas. The bath is described in FIG. 1. Average over many snapshots of the particle positions terms√ of noninteracting Bogoliubov excitations with energy 0 2 0 0 h¯2k2 around the impurity for three characteristic values of the impurity- k = ( ) + 2gBn , where = is the dispersion of k k k 2mB π 2/ bath coupling constant ln(kF a). The impurity is located in the center, = 4 h¯ mB free particles and gB / 2 is the 2D density-dependent and distances are in units of the healing length ξ. The local density ln(1 naB ) n(r) is estimated by summing particles over a square grid of size coupling constant of the Bose gas. The ground state of the L/100, where L is the size of the simulation box, and the color bar bath corresponds to the vacuum of excitations and has energy indicates the ratio n(r)/n over the bulk density n. At weak coupling, EB. The interaction Hamiltonian Hint is given by the sum of a impurities form polarons, i.e., almost free quasiparticles slightly mean-field shift and a term where the impurity is coupled to dressed by the medium (right and left upper panels). On the attractive the creation and annihilation operators of single excitations in branch as the coupling ln(kF a) > 0 is decreased, the impurity forms the Bose gas, a many-body bound state involving up to few tens of particles, and √ whose size is as large as the healing length (central upper panel). g n 0 = + √ iq·rα q α + α† . Hint gn e ( q −q ) (3) 2 L q q within the context of the Fröhlich model [60–62], but a π 2/ = 2 h¯ mr quantitatively precise description of 2D Bose polarons in the Here, g ln(1/na2 ) is the 2D effective coupling constant which mI mB strongly coupled regime is still lacking. contains the reduced mass mr = . It describes the scat- mI +mB Here, we use exact QMC methods [39,53,54] to study an tering processes between the impurity and the bath particles in impurity immersed in a 2D Bose superfluid and to compute terms of the scattering length a of the potential VI . The above the polaron energy, the effective mass, and the quasiparticle Hamiltonian H0 + Hint embodies the well-known Fröhlich residue for arbitrary coupling strength. Quantitatively signif- model which is expected to correctly describe the physics of icant deviations of the quasiparticle properties from pertur- Bose polarons in the weakly interacting limit where coupling bation theory are found already at weak-coupling strengths to multiple excitations of the bath can be neglected [49]. between the impurity and the bath. In the strongly interacting If E(p) is the energy of the impurity-bath system where the regime, the polaron loses the quasiparticle nature characteris- impurity has momentum p, the low-momentum expansion of tic of weak interactions: The wave-function residue vanishes the energy difference, indicating that the coherence is lost. In this regime, the impu- 2 rity is no longer free to move but, instead, is bound to a density − = μ + p +··· E(p) EB ∗ (4) perturbation which involves many particles from the bath (see 2mI Fig. 1). defines the binding energy μ of the impurity and its effective ∗ mass mI . By using perturbation theory one finds the follow- II. SYSTEM AND PERTUBATION THEORY ing results holding for mI = mB = m to lowest order in the coupling strength g of the interaction Hamiltonian H (see We consider an impurity of mass mI embedded in a 2D int Appendix A), Bose gas consisting of N atoms of mass mB at T = 0ina = N square box of size L with overall density n L2 . In the first μ 4 quantization formalism, the Hamiltonian of the system reads = , (5) μ0 ln(4π ) − 2ln(kF a) h¯2 N h¯2 N =− ∇2 + − ∇2 + . and H i VB(rij) α VI (riα ) 2mB 2mI m 1 [ln(4π ) − 2ln(k a )] i=1 i< j i=1 = 1 − F B . (6) ∗ π − 2 (1) m 2 [ln(4 ) 2ln(kF a)] h2k2 μ = ¯ F Here, the first two terms represent the kinetic and the inter- Here,√ we use 0 2m , involving the Fermi wave-vector action energies of the bosonic bath where particles interact kF = 4πn of a system having the same density n of the gas.
023405-2 STRONG COUPLING BOSE POLARONS IN A TWO- … PHYSICAL REVIEW RESEARCH 2, 023405 (2020)
the quasiparticle with μ>0. It is important to note that a two-body bound state with energy b exists for any value of the coupling constant ln (kF a). This is in contrast with the 3D polaron where the dimer state only appears as the s-wave scat- tering length turns positive on one side of the scattering res- onance [17]. In the weakly interacting regime | ln(kF a)|1, the QMC results are in good agreement with the prediction (5) of perturbation theory. Following the attractive branch, we note that the polaron energy is always much larger in absolute value than the dimer binding energy b. This is due to many- body effects which favor the formation of cluster states around the impurity involving many particles of the bath (see Fig. 1, central panel). In particular, in the vicinity of ln(kF a) ≈ 0, large fluctuations occur in our QMC simulations due to the formation of very deep many-body bound states which makes both the attractive and the repulsive branch of the polaron hard to follow further. The repulsive branch describes an excited state where the impurity repels the particles of the bath at FIG. 2. Polaron energy as a function of the coupling strength long distances but is unstable against cluster formation at short ln(kF a) for the attractive and repulsive branches (circles). The dashed distances. The state is well defined provided the typical size a line shows the dimer binding energy b, and the solid line shows the of bound states is small compared to the average interparticle perturbation result from Eq. (5). The coupling constant of the gas is distance k−1 but gets increasingly ill defined as the two length = . F g˜B 0 136. scales become comparable. This is exactly what we observe in our simulations where the excited state of the polaron is Moreover, the coupling strength of the impurity-bath interac- described using an appropriate choice of the wave function tion is expressed in terms of ln(kF a). These results were first used for importance sampling (see Appendix B for more derived in Ref. [62]. The same perturbation approach allows details). one to calculate the overlap Zp between the interacting and Furthermore, we study the mobility of the impurity by the noninteracting ground state of the impurity-bath system calculating its effective-mass m∗ as a function of the cou- with the impurity moving with momentum p. For the impurity pling strength ln kF a. The effective mass is determined by at rest (p = 0), one finds [62] (see Appendix A) computing the mean-square displacement of the impurity in imaginary time [39], [ln(4π ) − 2ln(kF aB )] Z0 = 1 − . (7) 2 [ln(4π ) − 2ln(k a)]2 m |rα (τ )| F = lim , (8) ∗ τ→∞ τ Note that the above results hold in the weak-coupling regime m 4D 2 | ln(kF a)|1 of the impurity-bath interaction. where D = h¯ /(2m) is the diffusion constant of a free parti- 2 2 cle and |rα (τ )| = |rα (τ ) − rα (0)| , being τ = it/h¯ the III. QMC RESULTS imaginary time of the QMC simulation. The effective mass is found by fitting the slope of the mean-square displacement In order to calculate the properties of the polaron for all for large values of τ. The residue Z0 of the polaron is obtained coupling strengths, we resort to QMC techniques. Details on from the one-body density matrix associated with the impu- the general method can be found in Refs. [39,53,63], whereas rity, an exhaustive discussion of the interatomic potentials used ψ (rα + r, r ,...,r ) in the simulations and on the trial wave function used for ρ(r) = T 1 N , (9) importance sampling are found in Appendix B. Simulations ψT (rα, r1,...,rN ) are performed for a gas of N identical particles and a single where ψ is the many-body guiding wave function of the impurity in a square box of size L with periodic boundary T QMC simulation. The above quantity is normalized to unity conditions. We choose the valueg ˜ = mBgB = 0.136 for the B h¯2 for r → 0, whereas its long-range limit gives the residue, dimensionless coupling constant in the bath. This value corre- | | sponds to ln(kF aB ) 45 and is typical for the experimental lim ρ(r) → Z0. (10) conditions of 2D Bose gases [64]. Furthermore, as in the per- r→∞ turbation theory study, we consider the case where impurity In Figs. 3 and 4, we show the results for the effective and particles in the bath have the same mass: mB = mI = m. mass and the quasiparticle residue, respectively. The calcu- We calculate the polaron energy μ from the direct calcula- lation of the residue Z0 is particularly sensitive to finite- tion of the ground-state energy of the bath with and without size effects which make the extrapolation to the thermody- the impurity μ = E(N, 1) − E(N), where N is the number of namic limit delicate. We have chosen different long-range particles in the bath. Results are shown in Fig. 2. In analogy asymptotic behaviors for the Jastrow terms entering the trial with the 2D Fermi polaron [10,58], we find two branches: One wave function (see Appendix B). In the bath, the long-range corresponds to the ground state of the attractive polaron with decay of boson-boson correlations is governed by phonons μ<0, and the second one corresponds to an excited state of as shown in Ref. [66]. For the impurity-boson correlations,
023405-3 ARDILA, ASTRAKHARCHIK, AND GIORGINI PHYSICAL REVIEW RESEARCH 2, 023405 (2020)
dressed by the particles from the bath. The perturbation caused by the impurity in the surrounding medium involves up to few tens of particles over a distance on the order of the healing length. We find a vanishing quasiparticle residue and a large effective mass which signal the transition to a many- body bound state (cluster state) without breaking of transla- tional symmetry (localization). A similar situation occurs for Fermi polarons [10,58] where Pauli exclusion principle only allows for the formation of a molecular state involving just one particle from the bath. A question which remains open also in the Fermi polaron case is whether the quasiparticle to bound-state transition is discontinuous or continuous. Our findings are also in contrast with Bose polarons in 3D and 1D. In fact, polarons in 3D remain well-defined quasiparticles up to the limit of resonant interactions [14,15], whereas, in 1D, they are never well-defined quasiparticles as the one-body density-matrix (9) decays to zero at long FIG. 3. Effective mass m∗ of the polaron as a function of the distances with algebraic law for any value of the coupling constant between the impurity and the bath. In this respect, coupling strength ln(kF a). The coupling constant of the Bose gas is 2D geometry is peculiar because the quasiparticle nature of g˜B = 0.136. polarons is rapidly suppressed by increasing the interaction strength. instead, we use the same functional form as from the Gross- Pitaevskii equation in the case of a static impurity. This choice of the trial wave function exhibits a fast convergence IV. EXPERIMENTAL IMPLEMENTATION of Z0 with increasing system sizes and allows us to keep In Ref. [64], a gas of 87Rb atoms in the hyperfine state finite-size effects under control. From Figs. 3 and 4 we |F = 1, m = 0 is confined in a 2D rectangular box with note that, even for the smallest reported values of the cou- dimensions Lx ≈ Ly ≈ 30 μm, at temperatures much below pling (| ln(kF a)|10), perturbation theory does not repro- / ∗ the Berezinskii-Kosterlitz-Thouless critical temperature. In duce the QMC results of m m and Z0. This is in contrast the transverse direction, a strong harmonic confinement is ap- with the results of the polaron energy reported in Fig. 2 plied with frequency ωz/(2π ) 4.6 kHz and by changing the and shows that higher-order terms, not accounted for by the number of trapped atoms the 2D density can be varied in the Fröhlich model, play an important role for these quantities n ≈ μ −2 | | range of 10–80 m . The 2D scattering length is given already at such large values of ln(kF a) [67]. In the regime = . − π / (3D) by aB 1 863 z exp ( 2 z aB ) (see Ref. [8]) in terms of of strong interactions, both the inverse effective mass and (3D) the quasiparticle residue become significantly smaller than the 3D√s-wave scattering length aB and the transverse length = / ω / (3D) the corresponding noninteracting values. Indeed, we find that z h¯ mB z. With the typical ratio of lengths z aB the polaron loses its quasiparticle nature as it gets more 30–50 reached in experiments the interaction strength is in the range ofg ˜B 0.10–0√.16, whereg ˜B is the dimensionless parameterg ˜ = mBgB = 8πa(3D)/ . Due to the exponential B h¯2 B z / (3D) dependence of aB on the ratio z aB , the 2D gas parameter 2 −30 −50 takes on very small values: naB 10 –10 . In our purely 2 = −40 2D simulations, we use the value naB 10 for the gas parameter of the bath, which corresponds to the effective 2D = 4π = . coupling strengthg ˜B / 2 0 136 close to the experi- ln(1 naB ) mental conditions of Ref. [64].
V. CONCLUSIONS We investigated the properties of Bose polarons in two dimensions. The polaron energy, effective mass, and quasi- particle residue have been calculated using QMC techniques for arbitrary coupling strength. We study the properties of the attractive and repulsive branch which correspond to the ground state and to a metastable state of the impurity. In the ground state, the polaron energy is much lower than the one of the two-body bound state, which, in 2D, is present for any
FIG. 4. Residue Z0 of the polaron as a function of the coupling value of the interaction strength. At stronger couplings, the strength ln(kF a)[65]. The coupling constant of the gas isg ˜B = impurity forms a bound state involving many particles from 0.136. the bath, which features a large effective mass and a vanishing
023405-4 STRONG COUPLING BOSE POLARONS IN A TWO- … PHYSICAL REVIEW RESEARCH 2, 023405 (2020) wave-function residue. A vanishing quasiparticle residue and energy μ0 and in terms of the wave-vector kF , is written as a large effective-mass signal the transition from a polaron to in Eq. (5)ofthemaintext. a many-body bound state without breaking of translational Effective mass. We assume that the energy of the bath symmetry. A similar behavior is found along the repulsive = + μ + p2 +··· with the impurity is written as E(p) EB 2m∗ , branch where a low-density bubble is formed around the holding at low-momenta p of the impurity. Thus, in terms impurity. However, this state rapidly becomes unstable against of the second-order energy correction, the polaron mass is cluster formation as the interaction strength is increased. Our renormalized as paper is important for the investigation of transport properties 2 E (2) − μ + gn in layered structures of ultracold atoms [68,69]aswellas 1 = 1 + 0 . ∗ lim (A4) layered solid-state materials [59]. m m p→0 p2 (2) The term E0 is given by ACKNOWLEDGMENTS ξ 2 1/2 (2) 2 n (k ) 1 This research was funded by the DFG Excellence Clus- E =−g , 0 L2 (kξ )2 + 2 (k) ter QuantumFrontiers. S.G. acknowledges funding from the k Provincia Autonoma di Trento. G.E.A. acknowledges funding 0 ξ 2 1/2 √ k = (k ) ξ = / whereweuse [ 2 ] with h¯ 2mgBn as the from the Spanish MINECO (Grant No. FIS2017-84114-C2-1- k (kξ ) +2 P). The Barcelona Supercomputing Center (The Spanish Na- healing length in the bath. In addition, we introduce the = (p−h¯k)2 + − p2 → tional Supercomputing Center-Centro Nacional de Supercom- quantity (k) 2m k 2m . At low-momenta p 0, putación) is acknowledged for the provided computational one can expand 1/ (k)as facilities (Grant No. RES-FI-2019-2-0033). · / 1 1 + h¯k p m h¯2k2 h¯2k2 (k) + k + k APPENDIX A: GROUND-STATE 2m 2m 2 PROPERTIES-PERTURBATION THEORY h¯ k2 p2/m2 + cos2 θ +··· , (A5) h¯2k2 + 3 The Fröhlich Hamiltonian in Eq. (3) of the main text reads 2m k √ where θ is the angle between p andh ¯k. After taking the sum 0 g n k † over k, we identify the first term in Eq. (A5) as the second- H = H0 + gn + √ exp(ik · rα ) (αk + α− ), 2 k order correction to the polaron energy μ whereas the second L k k (A1) term vanishes due to symmetry. Thus, one ends up with where H0 is the unperturbed Hamiltonian in Eq. (2), con- (2) − μ + E0 gn sisting of an impurity with momentum p and independent − 2 ξ 2 1/2 2 2 2 θ 2/ 2 Bogoliubov excitations with energy k. Hence, the unper- = g n (k ) h¯ k cos p m , 2 turbed ground state is represented by |0 =|p, 0k and cor- 2 ξ 2 + h¯ θk2 3 L (k ) 2 + k responds to the free impurity and the vacuum of Bogoliubov k 2m excitations. Relevant processes consist of states in which the and by performing the integration over momenta, one finds impurity is scattered by a single excitation. These states are 2 2 (2) 1 ln 1 naB p represented by |k =|p − hk¯ , 1k and correspond to unper- E − μ + gn =− . 0 2 2 (0) (p−h¯k)2 2 ln (1/na ) 2m turbed energies E = + EB + k. k 2mI Polaron energy. The Hamiltonian is split into H = H + The result (A4) is then given by 0 H . The energy expansion within perturbation theory is writ- int m 1 ln 1 na2 = (0) + (1) + (2) +··· (0) = p2 + = − B . ten as E0 E E E , where E EB ∗ 1 2 (A6) 0 0 0 0 2mI m 2 ln (1/na2 ) is the ground-state energy of the unperturbed system. The ∗ first- and second-order contributions to the energy are given Finally, in terms of the wave-vector kF , the ratio m/m is by written as in Eq. (6) of the main text. Quasiparticle residue. Within perturbation theory, one (1) = | | , E0 0 Hint 0 computes the correction to the ground state as | k|H |0 |2 k|H |0 E (2) = int . (A2) |p, 0 =|0 + int |k +··· , (A7) 0 (0) − (0) pert (0) − (0) k =0 E0 Ek k =0 E0 Ek
E (1) where we neglect terms of second order in Hint orthogonal By computing the matrix element 0 , one straightforwardly | obtains to 0 as well as higher-order contributions. The quasiparticle residue is defined as the square modulus of the overlap be- π 2 | (1) 4 h¯ n 1 tween the unperturbed state 0 and the normalized perturbed E = gn = , (A3) | , 0 m ln(1/na2 ) state Zp p 0 pert. Up to second-order contributions in Hint, one finds where we assumed equal masses for impurity and particles in | | | |2 = = μ 1 k Hint 0 the bath (mB mI m). The polaron energy is estimated Zp = = 1 − . p, 0|p, 0 (0) − (0) 2 using the lowest-order result (A3) which, in units of the pert pert k =0 E0 Ek
023405-5 ARDILA, ASTRAKHARCHIK, AND GIORGINI PHYSICAL REVIEW RESEARCH 2, 023405 (2020)
By carrying out the integral over momenta and by taking the The Jastrow term for boson-boson correlations is chosen of limit p → 0 similar to the case of the effective mass, one finds the following form: ⎧ ⎪I (k r), r < R , ln 1 na2 ⎨ 0 0 0 Z = − B . r , < , 0 1 2 (A8) fB(r) = A ln R0 r R (B2) ln (1/na2 ) ⎩⎪ a B −C + D , < / . B exp r r2 R r L 2 The above result, if written in terms of the wave-vector kF , = / / reduces to Eq. (7) of the main text. Here, A I0(k0R0 ) ln(R0 aB ) to ensure continuity of fB(r) at r = R0. Furthermore, the coefficients B–D are chosen such
that fB and its first derivative fB are continuous functions = / = APPENDIX B: TRIAL WAVE FUNCTIONS AND at the matching point R and fB(r L 2) 0, complying INTERATOMIC POTENTIALS with the periodic boundary conditions. The position R of the matching point is a parameter optimized by minimizing We describe the trial wave function which is used in QMC the energy in a variational calculation. As stated, the short- simulations as a guiding function for importance sampling and range part corresponds to the two-body scattering solution at to impose proper boundary conditions on the many-body state. zero energy, and the leading long-range part reproduces the In general, the trial wave function is written as a pair product phononic tail as predicted from hydrodynamic theory [66]. of Jastrow functions,