Arxiv:1805.09819V2 [Cond-Mat.Quant-Gas] 13 Jun 2018 Mained Elusive
Total Page:16
File Type:pdf, Size:1020Kb
Probing entanglement in a many-body-localized system Alexander Lukin,1 Matthew Rispoli,1 Robert Schittko,1 M. Eric Tai,1 Adam M. Kaufman,1, 2 Soonwon Choi,1 Vedika Khemani,1 Julian L´eonard,1 and Markus Greiner1, ∗ 1Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 2JILA, National Institute of Standards and Technology and University of Colorado, and Department of Physics, University of Colorado, Boulder, Colorado 80309, USA An interacting quantum system that is subject to disorder may cease to thermalize due to lo- calization of its constituents, thereby marking the breakdown of thermodynamics. The key to our understanding of this phenomenon lies in the system's entanglement, which is experimentally chal- lenging to measure. We realize such a many-body-localized system in a disordered Bose-Hubbard chain and characterize its entanglement properties through particle fluctuations and correlations. We observe that the particles become localized, suppressing transport and preventing the ther- malization of subsystems. Notably, we measure the development of non-local correlations, whose evolution is consistent with a logarithmic growth of entanglement entropy - the hallmark of many- body localization. Our work experimentally establishes many-body localization as a qualitatively distinct phenomenon from localization in non-interacting, disordered systems. INTRODUCTION A Number entanglement Configurational entanglement Isolated quantum many-body systems, undergoing uni- tary time evolution, maintain their initial global purity. However, the presence of interactions drives local ther- Subsystem A Subsystem B Subsystem A Subsystem B B malization: the coupling between any subsystem and its Thermal MBL Anderson remainder mimics the contact with a bath. This causes interacting number entanglement disordered the subsystem's degrees of freedom to be ultimately de- configurational entanglement scribed by a thermal ensemble, even if the full system is in a pure state [1{5]. A consequence of thermaliza- 0 AB 0 AB 0 AB tion is that local information about the initial state of 10 10 10 1 1 1 the subsystem gets scrambled and transferred into non- 10 10 10 Time 2 2 2 local correlations that are only accessible through global 10 10 10 3 3 3 observables [4{6]. 10 10 10 Disordered systems [7{16] can provide an exception to FIG. 1. Entanglement dynamics in non-equilibrium this paradigm of quantum thermalization. In such sys- quantum systems. (A) Subsystems A and B of an iso- tems, particles can localize and transport ceases, which lated system out of equilibrium entangle in two different ways: prevents thermalization. This phenomenon is called number entanglement stems from a superposition of states many-body localization (MBL) [6, 7, 17{22]. Experi- with different particle numbers in the subsystems and is gen- erated through particle motion across the boundary; config- mental studies have identified MBL through the persis- urational entanglement stems from a superposition of states tence of the initial density distribution [22{25] and two- with different particle arrangement within the subsystems and point correlation functions during transient dynamics requires both particle motion and interactions. (B) In the ab- [24]. However, while particle transport is frozen, the pres- sence of disorder, both types of entanglement rapidly spread ence of interactions gives rise to slow coherent many-body across the entire system due to delocalization of particles (left dynamics that generate non-local correlations, which are panel). The degree of entanglement and the timescales change inaccessible to local observables [26{28]. These dynamics drastically when applying disorder (central panel): particle localization spatially restricts number entanglement, yet in- are considered to be the hallmark of MBL and distinguish teractions allow configurational entanglement to form very it from its non-interacting counterpart, called Anderson slowly across the entire system. A disordered system with- localization [7{12]. Their observation, however, has re- out interactions shows only local number entanglement while arXiv:1805.09819v2 [cond-mat.quant-gas] 13 Jun 2018 mained elusive. the slow growth of configurational entanglement is completely We study these many-body dynamics by probing the absent (right panel). entanglement properties of an MBL system with fixed particle number [26{30]. We distinguish two types of entanglement that can exist between a subsystem and systems. Configurational entanglement implies that the its complement (Fig. 1A): Number entanglement implies configuration of the particles in one subsystem is corre- that the particle number in one subsystem is correlated lated with the configuration of the particles in the other. with the particle number in the other. It is generated It arises from a combination of particle motion and in- through tunneling across the boundary between the sub- teraction. The formation of particle and configurational 2 A B C Physical model Mott insulator 1xN initial state Initial state Quench Time evolution Pinning Readout J 45Er 45Er U 8Er 8Er 10211102 W atoms no atoms D E 1.5 1.4 t 100 = thermal prediction 1.3 ) ) 1 1 ( vN ( vN S S 1.0 1.2 0.7 site entropy site entropy 0.6 0.6 - - 1.1 1 = Thermal W 1.0J 0.4 n 0.5 P 0.5 n ( ) P Single MBL W 8.9J 0.2 Single 0.4 ( ) 0 1.0 0.3 012345 5 10 Site occupation W J ( ) 0.0 0.9 10 1 100 101 102 2 4 6 8 10 - Evolution time t Disorder strength W J () ( ) FIG. 2. Site-resolved measurement of thermalization breakdown. (A) One-dimensional Aubry-Andr´emodel with particle tunneling at rate J=~, on-site interaction energy U and quasi-periodic potential with amplitude W . (B) We prepare the initial state of eight unentangled atoms by projecting tailored optical potentials on a two-dimensional Mott insulator at 45Er lattice depth, where Er=h = 1:24 kHz is the recoil energy. (C) We create a non-equilibrium system by abruptly enabling tunneling dynamics. Following a variable evolution time, we project the many-body state back onto the number basis by increasing the lattice depth, and obtain the site-resolved atom number from a fluorescence image (SI). (D) We compute the (1) single-site von Neumann entropy SvN from the site-resolved atom number statistics (inset) after different evolution times (scaled with tunneling time τ = ~=J) in the presence of weak and strong disorder. (E) Probability p1 to retrieve the initial state (1) (inset) and SvN for different W , measured after 100τ evolution. The deviation from the thermal ensemble prediction for strong disorder signals the breakdown of thermalization in the system. All lines in (C-D) show the prediction of exact diagonalization calculations without any free parameters. Each data point is sampled from 197 disorder realizations (SI). Error bars denote the s.e.m. entanglement changes in the presence or absence of in- localized state. teractions and disorder in the system (Fig. 1B). In ther- mal systems without disorder, interacting particles de- localize and rapidly create both types of entanglement EXPERIMENTAL SYSTEM throughout the entire system. Contrarily, for Anderson localization, number entanglement builds up only locally In our experiments, we study MBL in the interacting at the boundary between the two subsystems. Here the Aubry-Andr´emodel for bosons in one dimension [31, 32], lack of interactions prevents the substantial formation of which is described by the Hamiltonian configurational entanglement. In MBL systems, number ^ X y entanglement builds up in a similarly local way as for H = − J a^i a^i+1 + h:c: Anderson localization. However, notably, the presence i (1) of interactions additionally enables the slow formation of U X X + n^ (^n − 1) + W h n^ , configurational entanglement throughout the entire sys- 2 i i i i i i tem. y In this work, we realize an MBL system and character- wherea ^i (^ai) is the creation (annihilation) operator for y ize these key properties: breakdown of quantum thermal- a boson on site i, andn ^i =a ^i a^i is the particle number op- ization, finite localization length of the particles, area-law erator on that site. The first term describes the tunneling scaling of the number entanglement, and slow growth of between neighboring lattice sites with the rate J=~, where the configurational entanglement that ultimately results ~ is the reduced Planck constant. The second term rep- in a volume-law scaling. Each property shows a contrast- resents the energy shift U when multiple particles occupy ing behavior when the system is prepared at weak disor- the same site. The last term introduces a site-resolved der in a thermalizing state. While the former three prop- potential offset, which is created with an incommensurate erties are also present for an Anderson localized state, lattice hi = cos (2πβi + φ) of period β ≈ 1:618 lattice the slowly growing configurational entanglement qualita- sites, phase φ, and amplitude W . In our experiment, we tively distinguishes our system from a non-interacting, achieve independent control over J, W , and φ (Fig. 2A). 3 Our experiments begin with a Mott-insulating state in low disorder depth (W = 1:0(1)J), the entropy grows the atomic limit with one 87Rb atom on each site of a over a few tunneling times and then reaches a station- two-dimensional optical lattice (Fig. 2B). The system is ary value (Fig. 2D). The stationary value is reduced for placed in the focus of a high-resolution imaging system deep disorder (W = 8:9(1)J) and remains constant over through which we project site-resolved optical potentials two orders of magnitude, up to several hundred tunneling [33]. We first isolate a single, one-dimensional chain from times. The lack of entropy increase indicates the absence the Mott insulator and then add the site-resolved poten- of heating in the system. The excellent agreement of the tial offsets Wi with the incommensurate lattice.