Arxiv:2107.06931V1 [Cond-Mat.Quant-Gas] 14 Jul 2021

Disorder in order: localization in a randomless cold atom system

Félix Rose1, 2, 3 and Richard Schmidt1, 2 1Max Planck Institute of Quantum Optics, Hans-Kopfermann-Str. 1, 85748 Garching, Germany 2Munich Center for Quantum Science and Technology (MCQST), Schellingstraße 4, 80799 Munich, Germany 3Technische Universität München, Physik Department, James-Franck-Strasse 1, 85748 Garching, Germany (Dated: July 14, 2021) We present a mapping between the Edwards model of disorder describing the motion of a single particle subject to randomly-positioned static scatterers and the Bose polaron problem of a light quantum impurity interacting with a Bose-Einstein condensate (BEC) of heavy atoms. The mapping offers an experimental setting to investigate the physics of Anderson localization where, by exploiting the quantum nature of the BEC, the time evolution of the quantum impurity emulates the disorder- averaged dynamics of the Edwards model. Valid in any space dimension, the mapping can be extended to include interacting particles, arbitrary disorder or confinement, and can be generalized to study many-body localization. Moreover, the corresponding exactly-solvable disorder model offers means to benchmark variational approaches used to study polaron physics. Here, we illustrate the mapping by focusing on the case of an impurity interacting with a one-dimensional BEC through a contact interaction. While a simple wave function based on the expansion in the number of bath excitations misses the localization physics entirely, a coherent state Ansatz combined with a canonical transformation captures the physics of disorder and Anderson localization.

Introduction. Coupling a particle to the collective an impurity immersed in a disorder-free Bose-Einstein excitations of a system with many degrees of freedom condensate (BEC). This mapping, illustrated in Fig.1, can radically alter the particle’s properties. While this can experimentally be realized using a mass-imbalanced paradigm has been first proposed by Landau and Pekar mixture of light impurities immersed in a bath of heavy to describe how the interaction between electrons and bosons and provides a new theoretical tool to include lattice phonons gives rise to quasiparticles named po- disorder effects in many-body approaches. larons [1], it has been extended to give insight into nu- Indeed, both the theoretical description of disorder and merous systems [2,3], including 3He–4He mixtures [4], Bose polarons face challenges of different origins. While semiconductors [5] and high-temperature superconduc- the rich physics of the Bose polaron problem arises from tors [6]. Although described by simple models such as the hard to capture many-body effects, performing the dis- Fröhlich Hamiltonian and in spite of intensive efforts [7– order average for even single-particle models is challeng- 22], a comprehensive solution of the polaron problem still ing. As such, a solution of the polaron problem can give escapes theory. The advent of ultracold atoms allows the insight into the corresponding disorder model, and con- realization of polaron models with high tunability and versely the disorder model can serve as an exactly solv- brought means to probe such models, as evidenced by able benchmark for polaron theories. To demonstrate the recent observation of the Bose polaron spectral func- this connection we show how a variational method ap- tion in impurity-boson mixtures [23–27]. plied to the study of impurity models can reproduce the In this letter, we bring together the seemingly discon- exact short-time solution of the corresponding disorder nected fields of polarons and disorder. Following An- model and discuss implications of the mapping for stud- derson’s realization that quantum interference can hinder the diffusion of a particle to the point that it be- comes localized [28], the interplay of disorder and quantum physics has been extensively studied, revealing in- tricate phenomena such as magnetoresistance [29, 30], coherent backscattering [31–33] and many-body localization [34–37]. While Anderson localization was first ob- arXiv:2107.06931v1 [cond-mat.quant-gas] 14 Jul 2021 served in wave systems [38, 39], the development of cold atom physics allowed to achieve the localization of matter waves [40, 41] and realize the quantum kicked rotor model [42, 43] which can be mapped onto the Anderson model of disorder [44]. Here, we explore an alternative way to study the Figure 1. The disorder-averaged evolution —represented by physics of disorder by establishing a mapping between the overline— of a particle subject to a random scattering po- the disorder-averaged motion of a single particle evolv- tential (left) is mapped onto the time evolution of an impurity ing through a random disorder potential, and that of immersed in a homogeneous, disorder-free BEC (right). 2 ies of polaron and disorder physics. Polaron model Disorder model The mapping. Many aspects of disorder physics such Impurity Particle as Anderson localization are universal. They may depend Heavy boson Static scatterer on dimension and symmetries but not on specific details Interspecies interaction Scattering potential of the Hamiltonian, leaving freedom as to which model N-boson state |{r }i Disorder configuration {r } of disorder one studies. We consider the Edwards model i i BEC state Sampling of disorder [45, 46] that describes the evolution of a single particle Quantum measurement hOˆ(t)i Disorder average hoˆ(t)i through a medium of N randomly-positioned static scatterers given by the Hamiltonian Table I. Correspondence between the heavy Bose polaron and Edwards model of disorder. pˆ2 N ˆ {ri} X HEd = + V (ˆr),V (ˆr) = v(ˆr − ri). (1) 2mI i=1 denotes the state where the N bosons have well-defined Here ˆr and pˆ are the position and momentum operators positions ri. The combined state of the system reads of the particle and v(ˆr − ri) is the potential created by a scatterer at site r . In one dimension and with contact |Ψi = |ψi ⊗ |BECi = |ψi ⊗ |{r }i. (6) i ˆ i interactions, this model is known as the random Kronig– {ri} Penney model (RKPM) [47, 48]. The randomness comes For infinitely massive bosons, mB/mI = ∞, the time from the positions ri of the scatterers, distributed e.g. evolution of each state |ψi ⊗ |{ri}i contributing to the uniformly in a volume Ω. Denoting hoˆ(t)i{ri} the expectation value of an observable oˆ at time t for a given superposition (4) can be determined exactly. The boson realization of disorder {r }, its disorder average is kinetic energy drops out and the total Hamiltonian com- i mutes with the bosonic position operators. Hence, the bosons remain in the state |{ri}i. Physically, the heavy hoˆ(t)i = hoˆ(t)i , (2) ˆ {ri} bosons’ positions are not affected by the interaction with {ri} the impurity. On the other hand, the impurity views −N d d the N localized bosons as scatterers at positions ri and where = Ω d r1 ··· d rN is the normalized in- {ri} evolves through the Edwards Hamiltonian (1) such that tegral over´ all possible´ scatterer positions r . i the system evolves into We now show that for any observable oˆ the disorder average (2) can be computed by means of a disorder-free h −iHˆ {ri}t i polaron model. In this model one considers a system |Ψ(t)i = e Ed |ψi ⊗ |{ri}i; (7) ˆ{r } where a mobile impurity of position ˆr and momentum pˆ i is immersed in a bosonic bath, described by the Hamil- i.e., the system evolves as a superposition over all possible tonian disorder realizations.

2 Hence, the expectation value of any observable of the ˆ pˆ X ˆ† ˆ 0 ˆ† ˆ ˆ H = + ωkbkbk + v(ˆr − r )br0 br0 (3) impurity O =o ˆ ⊗ 1 with respect to the state |Ψ(t)i, 2mI ˆ 0 k r iHˆ {ri}t −iHˆ {ri}t hOˆ(t)i = hψ|e Ed oeˆ Ed |ψi = hoˆ(t)i, (8) where ω = k2/2m describes the dispersion relation of k B ˆ{ri} bosons at momentum k annihilated (created) by the op- ˆ(†) erator bk , and mI and mB are the masses of the impurity realizes the disorder average hoˆ(t)i (2). All states in the and the bosons, respectively. Both species interact with superposition contribute equally to the measurement of the density-density interaction v(r). Oˆ, thus carrying out an average over all possible {ri} To establish the mapping, we quench the system by conﬁgurations. In other words, the fact that the BEC is preparing the impurity in a given wavepacket |ψi and a quantum superposition with equal weight of orthogonal the bosons in a Bose-Einstein condensate (BEC) of N states |{ri}i enables disorder averaging in an analog way non-interacting particles, using quantum systems such as ultracold atoms. The correspondence between the two models is summa- (ˆb† )N rized in TableI. We emphasize that hOˆ(t)i and hoˆ(t)i have √k=0 |BECi = |0i = |{ri}i (4) very diﬀerent meanings. While hOˆ(t)i represents a many- N! ˆ{ri} body measurement of the impurity evolving through the where interaction with a bath of heavy bosons, hoˆ(t)i corre- sponds to the measurement of the corresponding observ- ˆ† ˆ† br1 ··· brN able in the single-particle Edwards model averaged over |{ri}i = √ |0i (5) N! many classical realizations of disorder. 3

Generalizations. Most assumptions made for simplic- ity in the proof can be relaxed, as long as the crucial Chevy 20 Coherent state ingredient that the |{ri}i are eigenstates of the Hamilto- Exact nian remains valid. In particular, the proof applies to any Free ( ) dimension. Moreover, it is possible to include arbitrary 15 confining potentials or interactions between the bosons, 푔 = 0 and even prepare the boson bath in a mixed state. In Δ푟/휎 all these cases, the corresponding distribution of disor- 10 der would not be uniform but set by the boson state. Crucially, the mapping holds also for multiple particles. That way, it is possible to simulate transport prop- 5 erties of realistic metals or to consider interactions between the impurities in order to investigate the interplay 0 of many-body and disorder effects. The latter scenario 0 5 10 15 20 25 30 푡/푡휎 has recently been considered in a study of many-body lo- Figure 2. Time evolution of the width ∆r of the impurity calization, where the thermalization properties of mass- wavepacket. We compare the exact solution (red) to vari- imbalanced mixtures was investigated [49]. ational results based on a coherent state (yellow) or Chevy Ansatz (blue). The initial width of the wavepacket σ sets a The mapping holds when the bosons have a mass ex- 2 characteristic timescale tσ = 2mIσ /~ and the interaction and ceedingly large compared to the impurity. When this is −1 boson density n are g/Ω = 1.5~tσ , σn = 0.4. Free diffusion not the case, the feedback of the interaction with the obtained for g = 0 is shown for reference (black). impurity as well as eventual interactions among bosons will induce dynamics of the bath, defining a characteristic timescale τφ diverging with mB/mI. In the disorder exp(−|r|/2ξ) at large |r| with ξ a localization length. model, this would correspond to dynamical disorder and This contrasts with the free case (g = 0) where dif- cause dephasing. While this may hinder the observation fusion takes place. A simple observable that distin- of phenomena associated with disorder, one could con- guishes the two regimes is the width of the wavepacket p 2 versely take advantage of the presence of τφ to investigate ∆r = h(r − hri) i. For a localized wavepacket, ∆r re- the physics of dephasing in a controlled manner. Further- mains finite at all times, while in the free case ∆r grows more, the paradigm of using a controlled bath as an “aux- linearly with time. iliary disorder” can be reversed. Should one be interested We first compare the exact solution of the disorder in the bosons’ properties, the dynamics of the impurity model to the Chevy Ansatz, which includes for the time- can give information about the time correlations of the dependent wavefunction |Ψp(t)i at most one bosonic ex- bosons, akin to diffusion wave spectroscopy [46, 50]. citation [52, 53], i.e. Benchmarking variational solutions to Bose polarons. p So far, we have discussed how the exact mapping can |Ψ (t)i = α0(t)|pi ⊗ |BECi provide new approaches to disorder theory by enabling X ˆ† ˆ + αq(t)|p + qi ⊗ b−qb0|BECi (9) the realization of disorder in a controlled manner. How- q6=0 ever, one can also turn the correspondence around and use disorder theory to gain insight into polaron formation where αq(t) are time-dependent variational parame- by providing an exactly solvable limit. As an example, we ters [51, 54]. use the mapping to test variational solutions of the Bose The time-resolved spreading of the wavepacket is com- polaron problem. While widely used, the quality of such pared to the exact solution in Fig.2. At short times, we approximations is hard to gauge, given their nonpertur- observe good agreement with the exact solution. This is bative nature as there is no small control parameter. to be expected as, in variational methods, the deviation Specifically, we consider the one-dimensional Bose po- from the exact solution arises from the iterated projection laron, where a single impurity interacts with the bosonic of the Schrödinger equation and the associated error did bath through a contact interaction, v(r) = gδ(r), cho- not built up at early times. This is also reflected in the sen to be repulsive (g > 0) so that there are no boson- evolution of the spatially resolved wavepacket shown in impurity bound state. We turn our attention to the Fig.3, where the whole density profile is well described spread of a particle prepared in a Gaussian wavepacket by the Chevy Ansatz at short times. Since the Chevy ψ(r, t = 0) ∝ exp(−r2/2σ2). The corresponding disorder Ansatz specifically relies on including only a small num- problem can be solved exactly [48, 51] and can thus serve ber of excitations, it is most accurate at short times at as a means to evaluate variational approximations. which the polaron cloud is still developing. However, the Since all states are localized in presence of disor- Chevy Ansatz breaks down at larger times, and, rather der in one dimension, at any finite interaction g > 0 than the localization clearly observed in the exact solu- the wavepacket localizes at t → ∞, i.e., ψ(r, t) → tion, it predicts either diffusive or sub-diffusive behavior. 4

a) b) c)

휎 휎 휎 2 푡/푡 = 2 푡/푡 = 10 푡/푡 = 20 −1 10 |휓(퐫, 푡)| −2 10

−3 10 −5 0 5 −5 0 5 −5 0 5 Figure 3. Density proﬁle |ψ퐫/휎(r, t)|2 as function of r of the impurity퐫/휎 interacting with a BEC for increasing퐫/휎 times t. We compare the exact solution (red) to the variational result based on a coherent state (yellow) and Chevy Ansatz (blue). We use σn = 0.4, −1 g/Ω = 1.5~tσ . The free diﬀusion (g = 0) is shown for reference (black).

This prompts the search for more involved variational Application to the Anderson transition. A long- solutions. We consider an Ansatz based on the Lee– standing question in disorder physics is that of the An- Low–Pines (LLP) transformation [55] to transform into derson transition in three dimensions (3D). Unlike in one the comoving frame of the impurity. In the new frame dimension, where all states are localized, in 3D a mobility the impurity momentum p represents the conserved total edge separates low-energy localized and high-energy ex- momentum of the complete system. Impurity operators tended states [28]. The observation of the mobility edge are hence eliminated and the bosons evolve according to using the expansion of atomic matter waves through dis- the p-dependent Hamiltonian ordered speckle potentials has sparked recent theoretical [71, 72] and experimental interest [41, 73, 74]. Dis- 2 p − P kˆb† ˆb crepancies between theory and experiments remain, how- ˆ BEC k k k X ˆ† ˆ g X ˆ† ˆ Hp = + ωkbkbk + bkbk0 ever, as to the position of the mobility edge [72, 75]. On 2mI Ω k kk0 the theory side, the problem is difficult because it requires (10) going beyond the perturbative weak disorder regime and which now contains an transformation-induced interac- treating fully the speckle potential, while for experiments tion between bosons. Ω is the system volume. it is difficult to prepare a sufficiently narrow energy dis- The Hamiltonian (10) has been studied using various tribution of the atomic cloud. Ansätze [56–63]. We approximate the boson wavefunc- The polaron-disorder correspondence can help to ob- tion |Φp(t)i by a product of coherent states, tain further insight. A specific candidate is 6Li impu- 133 rities embedded in Cs with mB/mI = 22.1 [76, 77]. p −iφp(t) P βp(t)ˆb† −h.c. |Φ (t)i = e e k k k |0i (11) In particular, the Li-Cs interaction can be tuned by a Feshbach resonance at 889 G at which the Cs-Cs inter- p p with variational parameters φ (t) and βk (t). By consid- action is small (aCsCs ' 190a0). The cloud of impurities ering this wavefunction, we neglect the possibility of cor- expands through the bosonic medium as through a dis- relations induced between the infinitely massive bosons. ordered Edwards potential. Theory-wise, one does not The results for the wavepacket are shown in Figs.2 have to deal with the speckle potential, and can use the and3. Again, the exact solution is well reproduced at arsenal developed to tackle strongly interacting problems short times, unsurprisingly, since for a state close to the to study the dynamics of the mixture. One expects two BEC, Eq. (11) reduces to Eq. (9). However, the coher- issues when trying to study the localization transition in ent Ansatz directly addresses one of the shortcomings of this setting. A first limit is set by the aforementioned the Chevy Ansatz by allowing the description of a large decoherence due to the finite boson mass. A second con- number of bosons excited to small momenta and thus straint is given by the lifetime of the system, set e.g. by captures localization, although ξ is underestimated. three-body loss. This may guide the choice of the system These examples demonstrate how the exact mapping and, for instance, molecular mixtures that feature sup- can help to benchmark variational solutions by provid- pressed three-body loss might prove better candidates ing an exact reference solution of the polaron problem. despite having lower mass ratios than Li-Cs [78]. Indeed both the Chevy and coherent Ansätze can be Conclusion and outlook. We have shown that a improved systematically, respectively, by allowing more generic model of disorder, the Edwards model, is mapped bosonic excitations [10–13] or by considering Gaussian to the problem of a light impurity coupled to a BEC, states [64] with the mapping allowing to quantify the im- showing a deep connection between two seemingly very provement. This equally applies to other schemes such different problems. Our exact mapping has several impli- as non-Gaussian states or quantum Monte-Carlo [65–70]. cations. Experimentally, it offers an alternative setup to 5 investigate the Anderson transition. On the theory side, [5] M. E. Gershenson, V. Podzorov, and A. F. Morpurgo, it can help to develop reliable approximation schemes to Colloquium: Electronic transport in single-crystal or- describe the Bose polaron problem, while conversely of- ganic transistors, Rev. Mod. Phys. 78, 973 (2006). fering an efficient way to bypass numerically expensive [6] E. Dagotto, Correlated electrons in high-temperature su- perconductors, Rev. Mod. Phys. 66, 763 (1994). disorder averaging. However, the mapping applies to a [7] J. T. Devreese and A. S. Alexandrov, Fröhlich polaron much larger class of models. For instance, in the case of and bipolaron: recent developments, Rep. Prog. Phys. a mixture with both interactions and disorder, the map- 72, 066501 (2009). ping can be used to gauge away interactions, reducing it [8] A. G. Volosniev, H.-W. Hammer, and N. T. Zinner, Real- to a single particle problem in presence of two sources time dynamics of an impurity in an ideal Bose gas in a of disorder. Another possibility is the realization of more trap, Phys. Rev. A 92, 023623 (2015). complicated models of disorder, e.g. with long-range cor- [9] M. A. García-March, A. S. Dehkharghani, and N. T. Zin- ner, Entanglement of an impurity in a few-body one- relations, by tuning the boson Hamiltonian. dimensional ideal Bose system, J. Phys. B 49, 075303 A promising avenue of research could originate from (2016). the field theoretical treatment of the Bose polaron prob- [10] J. Levinsen, M. M. Parish, and G. M. Bruun, Impurity in lem. For instance, the Chevy Ansatz is equivalent to a a Bose-Einstein Condensate and the Efimov Effect, Phys. non-selfconsistent T -matrix approach [52, 53]. It would Rev. Lett. 115, 125302 (2015). [11] S. M. Yoshida, S. Endo, J. Levinsen, and M. M. Parish, be interesting to see how such methods compare to the di- Universality of an Impurity in a Bose-Einstein Conden- agrammatic treatment of disorder, and whether they can sate, Physical Review X 8, 011024 (2018). be used to shed new light on either side of the mapping. [12] Z.-Y. Shi, S. M. Yoshida, M. M. Parish, and J. Levinsen, Further intriguing applications would be the study of the Impurity-Induced Multibody Resonances in a Bose Gas, interplay of disorder and interactions, especially in lower Phys Rev. Lett. 121, 243401 (2018). dimensions where many-body methods such as the use of [13] J. Levinsen, L. A. Peña Ardila, S. M. Yoshida, and variational matrix product states could help understand M. M. Parish, Quantum behavior of a heavy impurity strongly coupled to a Bose gas, arXiv:2102.06368 [cond- thermalization and the spread of entanglement [79–82]. mat.quant-gas] (2021). A recent proposal of a different mapping between dis- [14] A. S. Dehkharghani, A. G. Volosniev, and N. T. Zin- ordered systems and interacting semi-metals brings to- ner, Coalescence of Two Impurities in a Trapped One- gether disorder and many-body physics in a solid state dimensional Bose Gas, Phys. Rev. Lett. 121, 080405 setting [83]. While our work focusses on cold atom ap- (2018). [15] J. Jager, R. Barnett, M. Will, and M. Fleischhauer, plications, the mapping presented in this work can be ex- Strong-coupling Bose polarons in one dimension: Con- tended to other condensed matter systems. Indeed, the densate deformation and modified Bogoliubov phonons, crucial ingredient behind the mapping is the absence of Phys. Rev. Research 2, 033142 (2020). boson kinetic energy. Hence, solid state realizations using [16] D. Dzsotjan, R. Schmidt, and M. Fleischhauer, Dynami- as a bath a system with flat bands —such as encountered cal Variational Approach to Bose Polarons at Finite Tem- within the quantum Hall effect— are also conceivable. peratures, Phys. Rev. Lett. 124, 223401 (2020). [17] N.-E. Guenther, R. Schmidt, G. M. Bruun, V. Gu- Acknowledgments. The authors thank Cheng Chin, rarie, and P. Massignan, Mobile impurity in a Bose- Arthur Christianen, Eugene Demler, Falko Pientka, Einstein condensate and the orthogonality catastrophe, Lode Pollet, Tao Shi, Binh Tran, Matteo Zacanti and Phys. Rev. A 103, 013317 (2021). Wilhelm Zwerger for fruitful discussions. They are grate- [18] P. Massignan, N. Yegovtsev, and V. Gurarie, Universal ful to Jean-Marc Luck for his insight regarding the Aspects of a Strongly Interacting Impurity in a Dilute RKPM and pointing out Ref. [48]. We acknowledge sup- Bose Condensate, Phys. Rev. Lett. 126, 123403 (2021). port from the Deutsche Forschungsgemeinschaft (DFG, [19] F. Isaule, I. Morera, P. Massignan, and B. Juliá- Díaz, Renormalization group study of Bose polarons, German Research Foundation) under Germany’s Excel- arXiv:2105.10801 [cond-mat.quant-gas] (2021). lence Strategy–EXC–2111–390814868. [20] R. Schmidt and T. Enss, Self-stabilized Bose polarons, arXiv:2102.13616 [cond-mat.quant-gas] (2021). [21] C. Franchini, M. Reticcioli, M. Setvin, and U. Diebold, Polarons in materials, Nat. Rev. Mater. (2021). [22] K. Seetharam, Y. Shchadilova, F. Grusdt, M. B. [1] L. Landau and S. Pekar, Effective mass of a polaron, J. Zvonarev, and E. Demler, Dynamical quantum Exp. Theor. Phys. 18, 419 (1948). Cherenkov transition of fast impurities in quantum [2] C. B. Duke and G. D. Mahan, Phonon-Broadened Impu- liquids, arXiv:2101.00030 [cond-mat.quant-gas] (2021). rity Spectra. I. Density of States, Phys. Rev. 139, A1965 [23] M.-G. Hu, M. J. Van de Graaff, D. Kedar, J. P. Cor- (1965). son, E. A. Cornell, and D. S. Jin, Bose Polarons in [3] A. S. Alexandrov and N. F. Mott, Polarons and Bipo- the Strongly Interacting Regime, Phys. Rev. Lett. 117, larons (World Scientific, Singapore, 1995). 055301 (2016). [4] G. Baym and C. Pethick, Landau Fermi-Liquid The- [24] N. B. Jørgensen, L. Wacker, K. T. Skalmstang, M. M. ory: Concepts and Applications (Wiley-VCH, New York, Parish, J. Levinsen, R. S. Christensen, G. M. Bruun, 1991). and J. J. Arlt, Observation of Attractive and Repulsive 6

Polarons in a Bose-Einstein Condensate, Phys. Rev. Lett. [43] C. Hainaut, A. Rançon, J.-F. Clément, I. Manai, 117, 055302 (2016). P. Szriftgiser, D. Delande, J. C. Garreau, and R. Chicire- [25] Z. Z. Yan, Y. Ni, C. Robens, and M. W. Zwierlein, anu, Experimental realization of an ideal Floquet disor- Bose polarons near quantum criticality, Science 368, 190 dered system, New J. Phys. 21, 035008 (2019). (2020). [44] S. Fishman, D. R. Grempel, and R. E. Prange, Chaos, [26] L. A. Peña Ardila, N. B. Jørgensen, T. Pohl, S. Giorgini, Quantum Recurrences, and Anderson Localization, Phys. G. M. Bruun, and J. J. Arlt, Analyzing a Bose polaron Rev. Lett. 49, 509 (1982). across resonant interactions, Phys. Rev. A 99, 063607 [45] S. F. Edwards, A New Method for the Evaluation of Elec- (2019). tric Conductivity in Metals, Phil. Mag. 3, 1020 (1958). [27] M. G. Skou, T. G. Skov, N. B. Jørgensen, K. K. Nielsen, [46] E. Akkermans and G. Montambaux, Mesoscopic Physics A. Camacho-Guardian, T. Pohl, G. M. Bruun, and J. J. of Electrons and Photons (Cambridge University Press, Arlt, Non-equilibrium quantum dynamics and formation Cambridge, 2007). of the Bose polaron, Nat. Phys. 17, 2005.00424 (2021). [47] M. Lax and J. C. Phillips, One-Dimensional Impurity [28] P. W. Anderson, Absence of Diffusion in Certain Random Bands, Phys. Rev. 110, 41 (1958). Lattices, Phys. Rev. 109, 1492 (1958). [48] T. Nieuwenhuizen, Exact electronic spectra and inverse [29] D. Y. Sharvin and Y. V. Sharvin, Magnetic-flux quan- localization lengths in one-dimensional random systems: tization in a cylindrical film of a normal metal, J. Exp. I. Random alloy, liquid metal and liquid alloy, Physica A Theor. Phys. Lett. 34, 272 (1981). 120, 468 (1983). [30] B. Pannetier, J. Chaussy, R. Rammal, and P. Gandit, [49] T. Grover and M. P. A. Fisher, Quantum disentangled Magnetic Flux Quantization in the Weak-Localization liquids, J. Stat. Mech: Theory Exp. 2014, 10010 (2014). Regime of a Nonsuperconducting Metal, Phys. Rev. Lett. [50] B. J. Berne and R. Pecora, Dynamic Light Scattering with 53, 718 (1984). Applications to Chemistry, Biology and Physics (John [31] Y. Kuga and A. Ishimaru, Retroreflectance from a dense Wiley, New York, 1976). distribution of spherical particles, J. Opt. Soc. Am. A 1, [51] See supplementary online material. 831 (1984). [52] S. P. Rath and R. Schmidt, Field-theoretical study of the [32] P.-E. Wolf and G. Maret, Weak Localization and Co- Bose polaron, Phys. Rev. A 88, 053632 (2013). herent Backscattering of Photons in Disordered Media, [53] W. Li and S. Das Sarma, Variational study of polarons Phys. Rev. Lett. 55, 2696 (1985). in Bose-Einstein condensates, Phys. Rev. A 90, 013618 [33] M. P. V. Albada and A. Lagendijk, Observation of Weak (2014). Localization of Light in a Random Medium, Phys. Rev. [54] P. Kramer and M. Saraceno, Geometry of the Time- Lett. 55, 2692 (1985). Dependent Variational Principle in Quantum Mechanics [34] L. Fleishman and P. W. Anderson, Interactions and the (Lecture Notes in Physics) (Springer, Berlin, Heidelberg, Anderson transition, Phys. Rev. B 21, 2366 (1980). 1981). [35] D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Metal [55] T. D. Lee, F. E. Low, and D. Pines, The Motion of Slow insulator transition in a weakly interacting many-electron Electrons in a Polar Crystal, Phys. Rev. 90, 297 (1953). system with localized single-particle states, Ann. Phys. [56] A. Shashi, F. Grusdt, D. A. Abanin, and E. Demler, 321, 1126 (2006). Radio-frequency spectroscopy of polarons in ultracold [36] V. Oganesyan and D. A. Huse, Localization of interacting Bose gases, Phys. Rev. A 89, 053617 (2014). fermions at high temperature, Phys. Rev. B 75, 155111 [57] F. Grusdt, A. Shashi, D. Abanin, and E. Demler, Bloch (2007). oscillations of bosonic lattice polarons, Phys. Rev. A 90, [37] R. Nandkishore and D. A. Huse, Many-Body Localization 063610 (2014). and Thermalization in Quantum Statistical Mechanics, [58] Y. E. Shchadilova, R. Schmidt, F. Grusdt, and E. Demler, Annu. Rev. Condens. Matter Phys. 6, 15 (2015). Quantum Dynamics of Ultracold Bose Polarons, Phys. [38] R. Dalichaouch, J. P. Armstrong, S. Schultz, P. M. Platz- Rev. Lett. 117, 113002 (2016). man, and S. L. McCall, Microwave localization by two- [59] B. Kain and H. Y. Ling, Hartree-Fock treatment of Fermi dimensional random scattering, Nature 354, 53 (1991). polarons using the Lee-Low-Pine transformation, Phys. [39] L. Ye, G. Cody, M. Zhou, P. Sheng, and A. N. Norris, Rev. A 96, 033627 (2017). Observation of bending wave localization and quasi mo- [60] E. Nakano, H. Yabu, and K. Iida, Bose-Einstein- bility edge in two dimensions, Phys. Rev. Lett. 69, 3080 condensate polaron in harmonic trap potentials in the (1992). weak-coupling regime: Lee-Low-Pines-type approach, [40] J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, Phys. Rev. A 95, 023626 (2017). P. Lugan, D. Clément, L. Sanchez-Palencia, P. Bouyer, [61] E. Yakaboylu, B. Midya, A. Deuchert, N. Leopold, and and A. Aspect, Direct observation of Anderson localiza- M. Lemeshko, Theory of the rotating polaron: Spectrum tion of matter waves in a controlled disorder, Nature 453, and self-localization, Phys. Rev. B 98, 224506 (2018). 891 (2008). [62] M. Drescher, M. Salmhofer, and T. Enss, Real-space [41] F. Jendrzejewski, A. Bernard, K. Müller, P. Cheinet, dynamics of attractive and repulsive polarons in Bose- V. Josse, M. Piraud, L. Pezzé, L. Sanchez-Palencia, Einstein condensates, Phys. Rev. A 99, 023601 (2019). A. Aspect, and P. Bouyer, Three-dimensional localiza- [63] M. Drescher, M. Salmhofer, and T. Enss, Theory of a tion of ultracold atoms in an optical disordered potential, resonantly interacting impurity in a Bose-Einstein con- Nat. Phys. 8, 398 (2012). densate, Phys. Rev. Res. 2, 032011 (2020). [42] A. Eckardt and E. Anisimovas, High-frequency approx- [64] T. Shi, E. Demler, and J. Ignacio Cirac, Variational study imation for periodically driven quantum systems from of fermionic and bosonic systems with non-Gaussian a Floquet-space perspective, New J. Phys. 17, 093039 states: Theory and applications, Ann. Phys. 390, 245 (2015). (2018). 7

[65] L. A. P. Ardila and S. Giorgini, Impurity in a Bose- [83] S. Sun and S. Syzranov, Equivalence of interact- Einstein condensate: Study of the attractive and repul- ing semimetals and low-density many-body systems sive branch using quantum Monte Carlo methods, Phys. to single-particle systems with quenched disorder, Rev. A 92, 033612 (2015). arXiv:2104.02720 [cond-mat.mes-hall] (2021). [66] L. A. P. n. Ardila and S. Giorgini, Bose polaron problem: [84] A.D. McLachlan, A variational solution of the time- Effect of mass imbalance on binding energy, Phys. Rev. dependent Schrodinger equation, Molecular Physics 8, A 94, 063640 (2016). 39 (1964). [67] F. Grusdt, G. E. Astrakharchik, and E. Demler, Bose [85] A. G. Basile and V. Elser, Equations of motion for su- polarons in ultracold atoms in one dimension: beyond perfluids, Phys. Rev. E 51, 5688 (1995). the Fröhlich paradigm, New J. Phys. 19, 103035 (2017). [86] J. Haegeman, J. I. Cirac, T. J. Osborne, I. Pižorn, H. Ver- [68] L. Parisi and S. Giorgini, Quantum Monte Carlo study of schelde, and F. Verstraete, Time-Dependent Variational the Bose-polaron problem in a one-dimensional gas with Principle for Quantum Lattices, Phys. Rev. Lett. 107, contact interactions, Phys. Rev. A 95, 023619 (2017). 070601 (2011). [69] E. Akaturk and B. Tanatar, Two-dimensional Bose polaron using diffusion Monte Carlo method, Int. J. Mod. Phys. B 33, 1950238 (2019). [70] J. Vlietinck, W. Casteels, K. V. Houcke, J. Tempere, J. Ryckebusch, and J. T. Devreese, Diagrammatic Monte Carlo study of the acoustic and the Bose–Einstein condensate polaron, New Journal of Physics 17, 033023 (2015). [71] D. Delande and G. Orso, Mobility Edge for Cold Atoms in Laser Speckle Potentials, Phys. Rev. Lett. 113, 060601 (2014). [72] M. Pasek, G. Orso, and D. Delande, Anderson Localiza- tion of Ultracold Atoms: Where is the Mobility Edge?, Phys. Rev. Lett. 118, 170403 (2017). [73] S. S. Kondov, W. R. McGehee, J. J. Zirbel, and B. De- Marco, Three-Dimensional Anderson Localization of Ul- tracold Matter, Science 334, 66 (2011). [74] G. Semeghini, M. Landini, P. Castilho, S. Roy, G. Spag- nolli, A. Trenkwalder, M. Fattori, M. Inguscio, and G. Modugno, Measurement of the mobility edge for 3D Anderson localization, Nat. Phys. 11, 554 (2015). [75] N. Cherroret, T. Scoquart, and D. Delande, Coherent multiple scattering of out-of-equilibrium interacting Bose gases, Ann. Phys. 2021, 168543 (2021). [76] D. S. Petrov and F. Werner, Three-body recombination in heteronuclear mixtures at finite temperature, Phys. Rev. A 92, 022704 (2015). [77] S. Häfner, J. Ulmanis, E. D. Kuhnle, Y. Wang, C. H. Greene, and M. Weidemüller, Role of the intraspecies scattering length in the Efimov scenario with large mass difference, Phys. Rev. A 95, 062708 (2017). [78] A. Ciamei, S. Finelli, M. Inguscio, A. Trenkwalder, and M. Zaccanti, Ultracold 6Li-53Cr Fermi mixtures with tun- able interactions, Bulletin of the American Physical So- ciety (2021). [79] J. H. Bardarson, F. Pollmann, and J. E. Moore, Un- bounded Growth of Entanglement in Models of Many- Body Localization, Phys. Rev. Lett. 109, 017202 (2012). [80] M. Serbyn, Z. Papić, and D. A. Abanin, Universal Slow Growth of Entanglement in Interacting Strongly Disor- dered Systems, Phys. Rev. Lett. 110, 260601 (2013). [81] K. Agarwal, E. Altman, E. Demler, S. Gopalakrishnan, D. A. Huse, and M. Knap, Rare-region effects and dynamics near the many-body localization transition, Ann. Phys. (Berl.) 529, 1600326 (2017). [82] M. Schreiber, S. S. Hodgman, P. Bordia, H. P. Lüschen, M. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Observation of many-body localization of interacting fermions in a quasirandom optical lattice, Science 349, 842 (2015). 8

Supplementary Material The spectrum is fixed by the boundary conditions; for instance periodic boundary conditions impose M = 1 or In this Supplementary Material, we provide the exact equivalently Tr M = 2, as det M = 1. solution to the random Kronig-Penney model (RKPM) and explain the time-dependent variational approach Once the spectrum is determined, one finds for each used to study the corresponding polaron problem. eigenvalue the amplitudes and phases Ai, φi defining the eigenstate. Using the continuity of ψ(x) together with Eq. (S4) one relates Ai+1 and φi+1 to Ai and φi. Thus EXACT SOLUTION OF THE RKPM by recursion AN and φN can be expressed as a function of φ0 and A0. Finally, a suitable value of φ0 is deter- The random Kronig-Penney model (RKPM) is defined mined numerically such that the boundary condition is by the one-dimensional Hamiltonian fulfilled, while the amplitudes are fixed by the normalization. As an illustration, we show the first eigenstates N ∂2 X thus obtained for a given realization of disorder in Fig.4. Hˆ = − x + g δ(x − x ), (S1) RKPM 2m i i=1 For attractive interactions (g < 0), not considered in this letter, this method remains valid. However, in describing a single particle in a box of size L with mass this case one also needs to consider eigenstates of neg- m and position x ∈ [0,L] interacting with fixed random ative energy E = −κ2/2m. On each subinterval, the scatterers with a contact potential of strength g, with corresponding wavefunctions then take the form ψ(x) = set to 1. We consider a given realization of disorder ~ A e−κx + B eκx, physically representing states that are defined by the positions of the N scatterers x , which i i i exponentially localized about the scatterers. The spec- we chose to be ennumerated as x < x and note for i i+1 trum and wavefunctions can again be determined by the convenience x = 0, x = L. 0 N+1 transfer matrix method. We seek to solve the Schrödinger equation

N 1 X Eψ(x) = − ψ00(x) + g δ(x − x )ψ(x). (S2) 2m i i=1

The general solution is found using a standard transfer TIME-DEPENDENT VARIATIONAL METHOD matrix method, see e.g. Ref. [48]. We consider only the case of a repulsive scattering potential, g > 0, for which E is positive, as can been seen by multiplying Eq. (S2) by In this section we briefly review how to determine the ψ(x)∗ and integrating over the box. On any subinterval time evolution of a ket |ψi under a Hamiltonian Hˆ us- ]xi, xi+1[ the eigenfunctions ψ(x) take the form ing the variational method; see Refs. [54, 64, 84–86] for a more in-depth discussion. We first provide a general dis- ψ(x) = Ai sin(kx + φi) (S3) cussion before applying the formalism to the variational states (9) and (11) introduced in the main text. with Ai, φi an amplitude and a phase and k ≥ 0 defined by E = k2/2m. While ψ is continuous, there is a jump We seek for the best approximation of the Schrödinger 0 ˆ in ψ at each scatterer xi, equation (i∂t − H)|ψi = 0 with |ψi being restrained to a variational manifold M. Let us assume that M = ψ0(x+) − ψ0(x−) = 2mgψ(x ) (S4) i i i {|ψ(zi)i, zi ∈ C} is parameterized by complex numbers z , with |ψ(z )i an holomorphic funtion of the z -s. Under as can been seen by integrating Eq. (S2) over [x−, x+]. i i i i i that assumption, the following three approaches yield the Constructing the vector Ψ(x) = (ψ0(x), kψ(x))T in- same equations of motions for the variational parameters dexed by the coordinate x, the conditions (S3) and (S4) z . The first is to minimize at all times the norm of can be rewritten as i the ket (i∂t − Hˆ )|ψi with respect to the z˙i. The second − + + − approach is to extremize the Lagrangian Ψ(xi+1) = R[k(xi+1 − xi)]Ψ(xi ), Ψ(xi ) = T Ψ(xi ), (S5) where i L = [hψ|(∂t|ψi) − (∂thψ|)|ψi] − hψ|Hˆ |ψi (S8) ! ! 2 cos(θ) − sin(θ) 1 2mg/k R(θ) = ,T = . (S6) sin(θ) cos(θ) 0 1 with respect to the variational parameters (zi, z˙i). The Hence, Ψ(L) = MΨ(0) with third approach is to project at every time step the Schrödinger equation onto the tangent space to M. We M = R[k(L − xN )]TR[k(xN − xN−1)] ··· TR[kx1]. (S7) use in practice the last method which we present below. 9

4 −1/2 (푥)/퐿 푛

휓 0 a) b) c)

푔/푔0 = 1 푔/푔0 = 10 푔/푔0 = 300 −4 0 1 0 1 0 1 푥/퐿 푥/퐿 푥/퐿 Figure 4. First three eigenstates ψn=1,2,3(x) of HˆRPKM for a given realization of disorder for three diﬀerent interaction −1 strengths g, expressed in units of g0 = (2mL) . The ground state wavefunction and the ﬁrst two excited states are shown in black, red and blue, respectively. We use periodic boundary conditions in a box of size L. The positions of the scatterers xi are represented by grey vertical lines. For weak disorder (a) the wavefunctions closely resemble the solutions of a free particle in a box while at larger interactions (b,c) the eigenstates become localized over distances smaller than the box size.

Gram matrix formulation more involved, and we present here the detailed deriva- tion. Indeed, when applying the Gram formalism, one The tangent space to M at |ψi is spanned by the vec- encounters two issues for the coherent Ansatz (11). First, ∗ tors it is not holomorphic (as both βk and βk appear), and, second, there is no straight-forward equation of motion ∂ for the phase, which represents a gauge degree of free- |∂iψi = |ψ(z)i. (S9) ∂zi dom. Hence the projected Schrödinger equation is satisﬁed if We tackle both issues at once by rather considering the and only if for all i variational state P ˆ† βkbk h∂iψ|(i∂t − Hˆ )|ψi = 0. (S10) |βi = N e k |0i, (S15)

We now deﬁne the energy of the state E and the Gram parameterized by the complex numbers N and βk, with matrix G using the overlaps of the vectors |∂iψi, N introduced for normalization. The Ansätze (11) and (S15) are completely equivalent provided N = ∂ P 2 ˆ ˆ −iφ− |βk| /2 E = hψ|H|ψi,Ej = E = h∂jψ|H|ψi, (S11) e k . ∂z∗ j For clarity we drop in this Section the t and p depen- Gij = h∂iψ|∂jψi. (S12) dencies. We rewrite Eq. (S15) by deﬁning the vectors ˆ bˆ†·β (β)k = βk and (bˆ)k = bk, such that |βi = N e |0i. Eq. (S10) can then be rewritten as † P 2 We also note S = (β) · β = k |βk| . X We use the Gram matrix formulation where from the iGijz˙j = Ei. (S13) j gradients of |βi,

In the specific case where G is invertible the equivalent ∂ 1 |∂ βi = |βi = |βi, (S16) form is N ∂N N ∂ X ˆ† iz˙ = [G−1] E . (S14) |∂kβi = |βi = bk|βi, (S17) i ij j ∂βk j we deduce the Gram matrix [Eq. (S12)] Otherwise there is some indeterminacy in the equations of motion, i.e. it is possible to find two different solutions † ! 0 0 1 N β S z˙i and z˙i fulfilling Eq. (S14) provided Gij(z ˙j − z˙j) = 0. G = e . (S18) N ∗β |N |2(δ + β · β†)

Application to coherent states In the above expression, the ﬁrst line and column stand for the N direction, while the rest of the matrix corre- While determining the equations of motion using the sponds to the momentum modes labelled by k. In the Gram formalism for the Chevy Ansatz (9) is straightfor- k-k sector, δ is understood as the identity matrix and † † ∗ ward, the case of the coherent state Ansatz is slightly β · β as the matrix with elements (β · β )k,k0 = βkβk0 . 10

p − PB[β] g X G is invertible with the inverse + + ω − k · β + β . k k m k Ω q q  1  1 + S − β† (S23) −1  N ∗  −S G =  1 1  e . (S19) − β δ Applying Eq. (S14) we obtain N |N |2 p − PB[β] g X iβ˙ = + ω − k · β + β , The energy [Eq. (S11)] is given by k k k m k Ω q q (S24) E = |N |2eSE , (S20) 0 p2 − P [β]2 iN˙ = N B . (S25) 2m with From these equations of motion we deduce that |N | 2 2 X (p − PB[β]) g X as well as S are constant in time, as expected since E = ( + ω )|β |2 + + β 0 k k k 2m Ω q the Hamiltonian conserves the number of particles and k q S = hβ|Nˆ|βi/hβ|βi = N. Hence the norm of |βi, (S21) hβ|βi = |N |2eS, is conserved. We can rewrite N = exp(−S/2) exp(−iφ) as P 2 and PB[β] = k k|βk| is the total momentum of the bosons. The gradients of E read p2 − P [β]2 φ˙ = B . (S26) 2mI ∂ S EN = E = N e E0, (S22) Using this new notation, (S15) is exactly equivalent ∂N ∗ to (11). The equations of motion thus obtained are the ∂ 2 S Ek = ∗ E = |N | e βkE0 same as those present in the literature, see e.g. Ref. [58]. ∂βk