Abstract
Clarifying the interplay of interactions and disorder is fundamental to the understanding of many quantum systems, including superfluid helium in porous media1, granular and thin-film superconductors2,3,4,5 and light propagating in disordered media6,7,8. One central aspect for bosonic systems9,10,11 is the competition between disorder, which tends to localize particles, and weak repulsive interactions, which instead have a delocalizing effect. As the required degree of independent control of disorder and of interactions is not easily achievable in most available physical systems, a systematic experimental investigation of this competition has not so far been possible. Here we use a degenerate Bose gas with tunable repulsive interactions in a quasiperiodic lattice potential to study this interplay in detail. We characterize the entire delocalization crossover through the study of the average local shape of the wavefunction, the spatial correlations and the phase coherence. Three different regimes are identified and compared with theoretical expectations12,13,14,15,16,17: an exponentially localized Anderson glass and the formation of locally coherent fragments as well as a coherent, extended state. These results provide insight into the role of weak repulsive interactions in disordered bosonic systems, but our approach should also enable investigations of disordered systems with interactions in the strongly correlated regime18,19,20.
Similar content being viewed by others
Main
The interplay of disorder and interactions lies at the heart of the behaviour of many physical systems. Notable examples are the transitions to insulators observed in superconductors and metals2,3,4,5, quantum Hall physics21, electrical conduction in DNA (ref. 22), or light propagation in nonlinear disordered media7,8. An important step towards their full comprehension is understanding disordered bosonic systems at zero temperature, where a competition between disorder and weak repulsive interactions is expected. Indeed, whereas disorder tends to localize non-interacting particles, giving rise to Anderson localization23, weak repulsive interactions can counteract this localization to minimize the energy. Eventually, interactions can screen the disorder and bring the system towards a coherent, extended ground state, that is, a Bose–Einstein condensate (BEC). In many years of research, mainly theoretical predictions have been made about the properties of the complex phases expected to appear as a result of this competition9,10,11,12,13,14,15,16,17. A systematic experimental study has not so far been possible, because on the one hand interactions in condensed-matter systems are strong but difficult to control1, and on the other hand in photonic systems only nonlinearities corresponding to attractive interactions7,8 have been explored in experiments. Instead, ultracold atoms in disordered optical potentials are a promising system for such investigations15,24,25, and have already enabled the observation of Anderson localization for bosons in the regime of negligible interactions26,27. Using one of these systems in a disordered lattice, we characterize the whole crossover from the regime of disorder-induced localization to that of Bose–Einstein condensation by tuning repulsive interactions in a controlled manner. The simultaneous measurement of localization properties, spatial correlations and phase coherence properties, and the comparison with the predictions of a theoretical model, enable us to identify the different regimes of this delocalization crossover.
The system used consists of a three-dimensional degenerate Bose gas of 39K in a one-dimensional quasiperiodic potential, which is generated by perturbing a strong primary optical lattice of periodicity d=π/k1 with a weak secondary lattice of incommensurate periodicity π/k2 (k=2π/λ, where λ is the wavelength of the light generating the lattice). The corresponding Hamiltonian is characterized by the site-to-site tunnelling energy J of the primary lattice, which is kept fixed in the experiment, and the disorder strength Δ. The interatomic interactions can be controlled by changing the atomic s-wave scattering length a by means of a Feshbach resonance28, which in turn determines the mean interaction energy per particle Eint (see Methods). In the case of non-interacting atoms, such a system is a realization of the Aubry–André model29, which shows an Anderson-like localization transition for a finite value of the disorder Δ/J=2. Above the transition, the non-interacting eigenstates of the potential are exponentially localized owing to the quasiperiodic perturbation of the lattice on-site energies and the energy spectrum is split into ‘minibands’13,30. The localization properties in this case have been studied experimentally in detail in ref. 27, where it was seen that several low-lying eigenstates, separated on average by d/(β−1)≈4.4d, where β=k2/k1, are typically populated in the experiment. Adding weak interactions, the different regimes that appear as a result of the interplay of disorder and interactions can be explored. For very weak repulsive interactions, the occupation of several eigenstates in the lowest miniband is favoured (Fig. 1a). This regime, in which several exponentially localized states coexist without phase coherence, is often identified with an Anderson glass11,15 (AG). As Eint is increased, coherent fragments, which extend over more than one well of the quasiperiodic potential, are expected to form (Fig. 1b). In this case, global phase coherence would not yet be restored, and the local shape of the states might be modified. Some authors have called this regime a ‘fragmented BEC’ (fBEC; ref. 12). Finally, for large enough Eint a single, extended phase-coherent state is expected to be formed, that is, a macroscopic BEC (Fig. 1c).
The system is prepared by first loading an interacting condensate adiabatically from the ground state of a harmonic trap into the quasiperiodic lattice. The interaction energy is then slowly changed to its final value Eint, and the confining potential is reduced. This process is adiabatic for most of the parameter range explored until Eint becomes sufficiently low for the system to enter the fully localized regime. Here, several independent low-lying excited states are populated even when it would be energetically favourable to populate just the ground state. This loss of adiabaticity is seen experimentally as a transfer of energy into the radial direction (see Supplementary Information).
The system is characterized in detail by analysing its momentum distribution, which is recovered by taking an image of the condensate after a long ballistic expansion without interactions (see Methods). From the momentum distribution and derived Fourier transforms, of which we show examples in Fig. 2, we extract the local shape of the wavefunction, spatial correlations and phase coherence properties for different values of Δ/J and Eint/J. The system can be approximately described as the superposition of states with the same envelope separated by 4.4d. First, the mean extension of individual states can be quantified by measuring the root-mean-squared width of the momentum distribution (Fig. 2a,b). A smaller (larger) width indicates a more extended (localized) state. Next, the mean local shape of the wavefunction on a length scale of 4.4d is extracted from the Fourier transform of the square root of the momentum distribution. From a fit to a generalized exponential function, the localization exponent α is recovered (see also Methods), as shown in Fig. 2c,d. The measured momentum width and exponent are shown in Fig. 3. We find that for very small Eint the states are exponentially localized, because α≈1, and the momentum width is large, consistent with the AG regime. Increasing Eint, the width decreases and the exponent increases to α≈2. Repulsive interactions therefore delocalize the system as expected, or, alternatively, the localization crossover is shifted to higher values of the disorder strength Δ/J when interactions are introduced into the system. The position of the delocalization crossover is in good agreement with the expectations of a simple screening argument14. The increasing interaction energy serves to smooth over the disordering potential in the occupied sites, providing a flatter energetic landscape on which more extended states can form. The centre of the crossover is therefore expected to occur when Eint is comparable to the standard deviation of energies in the lowest miniband of the non-interacting spectrum, 0.05Δ (white line in Fig. 3, see also Supplementary Information).
The correlation properties of neighbouring states can be extracted from a Fourier transform of the momentum distribution itself, which gives the spatially averaged first-order correlation function g(x) (see Methods). In Fig. 4a,b, g(x) at x=4.4d is shown for both the experiment and a ground-state theory that we have developed, with generally good agreement. In the localized regime, the correlation is exactly zero in the theory, because no neighbouring states are occupied. In contrast, the correlation is finite in the experiment owing to the occupation of neighbouring localized states arising from the non-adiabatic loading, but is small because the states are independent. As Eint is increased, the correlation features a crossover towards larger values, signalling that coherence is progressively established locally over distances of at least 4.4d. The shape of the crossover in the experiment is again in qualitative agreement with the expectation of the screening argument above (see the white line in Fig. 4a).
Finally, information about the phase coherence of neighbouring states can be obtained by measuring the phase φ of the interference pattern in the momentum distribution for repeated runs of the experiment with the same parameters (see Methods for details). If the states are not phase locked, φ changes almost randomly at each repetition of the experimental sequence. In Fig. 4d we show the standard deviation of φ, estimated from a large number of repetitions of the experiment, for fixed Δ/J=12. We see a slight decrease of the phase fluctuations with increasing Eint, which nevertheless remain relatively large in the crossover region where the correlation increases (Fig. 4c). The fluctuations finally drop to the background value only when Eint is comparable to the full width of the lowest miniband of the non-interacting spectrum, 0.17Δ. These observations confirm that in the localized regime the states are totally independent, which together with the localization properties (Fig. 3) indicates that the system can indeed be described as an AG (refs 11, 15). The system crosses a large region of only partial coherence while becoming progressively less localized as Eint is increased. This is consistent with the formation of locally coherent fragments expected for an fBEC. An analogous fragmentation behaviour was reported in ref. 31. Finally, the features of a single extended, coherent state are seen, that is, a BEC.
In the mean-field theory, boundaries between the different regimes expected for the system can be defined (see Methods). In particular, the crossover from the AG phase to an fBEC (white line in Fig. 4b) occurs when g(4.4d) starts to increase. Similarly, the orange line in Fig. 4b shows where the fragments are locked together in phase to form a single macroscopic condensate for very large interactions. The generally good agreement between the experimental observables and theory indicates that our system is well described by the mean-field theory for most of the parameter space explored experimentally.
In conclusion, we have provided the first experimental characterization of the localization, correlation and coherence properties of the various regimes owing to the competition of disorder and weak repulsive interactions in a bosonic system. Other aspects of the delocalization crossover worth further study are, for example, the detailed properties of the ground state of the AG regime, which was not possible to study in the present set-up, and the presence of a superfluid–insulator transition at the BEC–fBEC boundary analogous to the one observed in superconductors2. Regarding the latter, in transport experiments analogous to the ones described in ref. 27 we have been able to verify that the AG and fBEC regimes are not inconsistent with being insulating, as is the case in the regime of vanishing Eint (see Supplementary Information). Finally, it would be appealing to use the present system and the correlation analysis introduced here to explore the regime of strong correlations, Eint≫J, which could be reached by using a quasi-one-dimensional system with strong radial confinement. There, another elusive insulating phase due to the cooperation of disorder and interactions, the so-called Bose-glass phase, is expected to appear, although there is debate on the exact shape of the phase diagram11,15,18,19,32.
Methods
Condensate with tunable interactions.
A 39K condensate of about N=20,000 atoms with an s-wave scattering length of 250 a0, where a0=52.9 pm is the Bohr radius, is prepared in a harmonic optical trap. The condensate is loaded into the quasiperiodic potential while the optical trap is decompressed in about 250 ms to reduce the harmonic confinement, and a gravity-compensating magnetic-field gradient is added. At the same time, the scattering length a is changed by means of a broad Feshbach resonance to values ranging from a≤0.1 a0 to about a=300 a0 (ref. 28).
Quasiperiodic potential.
The quasiperiodic potential is created by two vertically oriented laser beams in standing-wave configuration. The primary lattice is generated by a Nd:YAG laser with a wavelength of λ1=1,064.4 nm and has a strength of s1=V1/ER,1=10.5 (corresponding to J/h=79 Hz), as measured in units of the recoil energy ER,1=h2/(2M λ12). The secondary lattice is generated by a Ti:sapphire laser of wavelength λ2=866.6 nm, the strength being adjustable up to s2=V2/ER,2=1.7. Both beams are focused onto the condensate with a beam waist of about 150 μm. The lattice lasers give a harmonic confinement of in the radial direction. In the vertical (axial) direction, a weak confinement of 5 Hz is given by a weak optical trap as well as by a curvature from the gravity-compensating magnetic field.
Energy scales.
In the tight-binding limit, the hopping energy J and disorder strength Δ can be estimated as and Δ=0.5s2β2(1.0264 exp{−2.3624/s10.59})ER,1 (ref. 30). The experimental uncertainty on Δ/J is around 15%. We estimate that around 30 lattice sites, corresponding to about seven localized states, are populated during the loading of the lattice. We then define a mean interaction energy per particle , where g=4πℏ2a/m, m is the atomic mass and ϕ(r) is a Gaussian approximation to the on-site Wannier function. We include coupling in the radial directions of our system, with the consequence that the interaction energy is nonlinear in the scattering length. Though this definition of the energy is strictly valid only in the localized regime, comparison with a numerical simulation of our experimental procedure has shown that it is a good approximation for all values of the scattering length up to an error of 30%. Note that the potential energy from the residual harmonic confinement is approximately 3×10−3 J over a distance 4.4d.
Momentum-distribution analysis.
The images of the momentum distribution are taken by absorption imaging with a CCD (charge-coupled device) camera after 36.5 ms ballistic expansion. At the time of release, the scattering length is set to below 1 a0 in less than 1 ms and kept there until the Feshbach magnetic field is switched off 10 ms before taking the image—at this point, the system has expanded a sufficient amount to minimize the effect of interactions. For such a free expansion, the image is approximately the in-trap momentum distribution . The acquired images are integrated along the radial direction to obtain a profile. In momentum space, the width of the central peak is calculated by taking the root-mean-square width within the first Brillouin zone. Owing to the quasiperiodic lattice potential, for a sufficiently homogeneous system the in-trap wavefunction can be decomposed into copies of a single state with real and non-negative envelope ξ(x)∼exp(−|x/L|α), spaced by 4.4d. Therefore, in momentum space, , where is an interference term, and ξ(x) can be extracted from a Fourier transform of (see also Supplementary Information). We fit to the sum of two generalized exponential functions A exp(−|(x−xc)/L|α)·[1+B cos(k1(x−xc)+δ)], where xc denotes the centre of each of these functions, spaced by 4.4d. From this fit, the exponent α is recovered. In addition, from the Wiener–Khinchin theorem, the momentum distribution can be expressed in terms of the first-order correlation function as . By taking the Fourier transform of the momentum distribution itself, we can therefore recover the spatially averaged correlation function . With the same fitting function as above, we evaluate the spatially averaged correlation between two states 4.4 lattice sites apart, A2/A1. Experimentally, the correlation function saturates at a value around 0.5 owing to the finite momentum resolution. The fluctuations in phase between neighbouring states are seen as a fluctuation of the phase φ of the interference pattern of the momentum distribution, which is directly extracted from a fit (see also Supplementary Information). The two-dimensional graphs in Figs 3 and 4 were generated by linearly interpolating a total of 130 averaged datapoints at nine different values of disorder, changing the interactions. Typical experimental scatter and statistical errors are seen in Fig. 4c.
Theory of the ground state.
The theoretical calculations presented in the paper rely on a mean-field approach. This is an effective one-dimensional model, which also partially includes the radial-to-axial coupling, and is known to provide an accurate description in the two limiting cases of Anderson localization and a BEC. The boundaries between the different regimes shown in Fig. 4b are obtained by analysing the correlation function g(x) and the density distribution. In the theory we define the AG phase as the one in which the correlation g(4.4d) is zero. To enter the fBEC phase, we require g(4.4d)>0, which implies that coherent fragments composed of adjacent localized states can start to form. For increasing Eint the extension of the fragments increases, until most of the system remains in a single component, which corresponds to a macroscopic BEC. To define the boundary between the fBEC and the BEC regimes, we first identify as fragments the parts of the system separated by low-density regions for which an applied relative phase twist does not affect the energy of the system. When one single macroscopic fragment forms, we assume the system to be in the BEC regime. A more detailed description of the theoretical methods can be found in the Supplementary Information.
References
Reppy, J. D. Superfluid helium in porous media. J. Low Temp. Phys. 87, 205–245 (1992).
Dubi, Y., Meir, Y. & Avishai, Y. Nature of the superconductor–insulator transition in disordered superconductors. Nature 449, 876–880 (2007).
Beloborodov, I. S., Lopatin, A. V., Vinokur, V. M. & Efetov, K. B. Granular electronic systems. Rev. Mod. Phys. 79, 469–518 (2007).
Goldman, A. M. & Marković, N. Superconductor–insulator transitions in the two-dimensional limit. Phys. Today 51, 39–44 (1998).
Phillips, P. & Dalidovich, D. The elusive Bose metal. Science 302, 243–247 (2003).
Pertsch, T. et al. Nonlinearity and disorder in fibre arrays. Phys. Rev. Lett. 93, 053901 (2004).
Schwartz, T., Bartal, G., Fishman, S. & Segev, M. Transport and Anderson localization in disordered two-dimensional photonic lattices. Nature 446, 52–55 (2007).
Lahini, Y. et al. Anderson localization and nonlinearity in one-dimensional disordered photonic lattices. Phys. Rev. Lett. 100, 013906 (2008).
Giamarchi, T. & Schulz, H. J. Anderson localization and interactions in one-dimensional metals. Phys. Rev. B 37, 325–340 (1988).
Fisher, D. S. & Fisher, M. P. A. Onset of superfluidity in random media. Phys. Rev. Lett. 61, 1847–1850 (1988).
Scalettar, R. T., Batrouni, G. G. & Zimanyi, G. T. Localization in interacting, disordered, Bose systems. Phys. Rev. Lett. 66, 3144–3147 (1991).
Lugan, P. et al. Ultracold Bose gases in 1D disorder: From Lifshits glass to Bose–Einstein condensate. Phys. Rev. Lett. 98, 170403 (2007).
Roux, G. et al. Quasiperiodic Bose–Hubbard model and localization in one-dimensional cold atomic gases. Phys. Rev. A 78, 023628 (2008).
Lee, D. K. K. & Gunn, J. M. F. Bosons in a random potential: Condensation and screening in a dense limit. J. Phys. Condens. Matter 2, 7753–7768 (1990).
Damski, B., Zakrzewski, J., Santos, L., Zoller, P. & Lewenstein, M. Atomic Bose and Anderson glasses in optical lattices. Phys. Rev. Lett. 91, 080403 (2003).
Schulte, T. et al. Routes towards Anderson-like localization of Bose–Einstein condensates in disordered optical lattices. Phys. Rev. Lett. 95, 170411 (2005).
Falco, G. M., Nattermann, T. & Pokrovsky, V. L. Localized states and interaction-induced delocalization in Bose gases with quenched disorder. Europhys. Lett. 85, 30002 (2009).
Fisher, M. P. A., Weichman, P. B., Grinstein, G. & Fisher, D. S. Boson localization and the superfluid–insulator transition. Phys. Rev. B 40, 546–570 (1989).
Fallani, L., Lye, J. E., Guarrera, V., Fort, C. & Inguscio, M. Ultracold atoms in a disordered crystal of light: Towards a Bose glass. Phys. Rev. Lett. 98, 130404 (2007).
Pasienski, M., McKay, D., White, M. & DeMarco, B. Disordered insulator in an optical lattice. Preprint at http://arxiv.org/abs/0908.1182 (2009).
Martin, J. et al. The nature of localization in graphene under quantum Hall conditions. Nature Phys. 5, 669–674 (2009).
Endres, R. G., Cox, D. L. & Singh, R. R. P. Colloquium: The quest for high-conductance DNA. Rev. Mod. Phys. 76, 195–214 (2004).
Anderson, P. W. Absence of diffusion in certain random lattices. Phys. Rev. 109, 1492–1505 (1958).
Fallani, L., Fort, C. & Inguscio, M. Bose–Einstein condensates in disordered potentials. Adv. At. Mol. Opt. Phys. 56, 119–160 (2008).
Sanchez-Palencia, L. & Lewenstein, M. Disordered quantum gases under control. Nature Phys. 6, 87–95 (2010).
Billy, J. et al. Direct observation of Anderson localization of matter waves in a controlled disorder. Nature 453, 891–894 (2008).
Roati, G. et al. Anderson localization of a non-interacting Bose–Einstein condensate. Nature 453, 895–898 (2008).
Roati, G. et al. 39K Bose–Einstein condensate with tunable interactions. Phys. Rev. Lett. 99, 010403 (2007).
Aubry, S. & André, G. Analycity breaking and Anderson localization in incommensurate lattices. Ann. Israel Phys. Soc. 3, 133–140 (1980).
Modugno, M. Exponential localization in one-dimensional quasi-periodic optical lattices. New J. Phys. 11, 033023 (2009).
Chen, Y. P. et al. Phase coherence and superfluid–insulator transition in a disordered Bose–Einstein condensate. Phys. Rev. A 77, 033632 (2008).
Pollet, L., Prokof’ev, N. V., Svistunov, B. V. & Troyer, M. Absence of a direct superfluid to Mott insulator transition in disordered Bose systems. Phys. Rev. Lett. 103, 140402 (2009).
Acknowledgements
We thank B. Altshuler and M. Larcher for discussions, S. Müller for experimental contributions and all the colleagues of the Quantum Gases group at LENS. This work has been supported by the EU (MEIF-CT-2004-009939), by the ESA (SAI 20578/07/NL/UJ), by the ERC through the Starting Grant QUPOL and by the ESF and INO-CNR through the EuroQUASAR program.
Author information
Authors and Affiliations
Contributions
B.D., G.R. and M.Z. carried out the experiments; B.D. and C.D. analysed data; M.F. contributed to the experiment; M.M. carried out the theoretical calculations; B.D. wrote the manuscript; M.I. and G.M. codirected the project; all authors participated in discussions on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Information (PDF 1047 kb)
Rights and permissions
About this article
Cite this article
Deissler, B., Zaccanti, M., Roati, G. et al. Delocalization of a disordered bosonic system by repulsive interactions. Nature Phys 6, 354–358 (2010). https://doi.org/10.1038/nphys1635
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys1635
This article is cited by
-
Many-body dynamical delocalization in a kicked one-dimensional ultracold gas
Nature Physics (2022)
-
Universality and quantum criticality in quasiperiodic spin chains
Nature Communications (2020)
-
Precise programmable quantum simulations with optical lattices
npj Quantum Information (2020)
-
Localisation of weakly interacting bosons in two dimensions: disorder vs lattice geometry effects
Scientific Reports (2019)
-
Bogoliubov Theory for a Superfluid Bose Gas Flowing in a Random Potential: Stability and Critical Velocity
Journal of Low Temperature Physics (2018)