We use thermalization indicators and numerical linked cluster expansions to probe the onset of many-body localization in a disordered one-dimensional hard-core boson model in the thermodynamic limit. We show that after equilibration following a quench from a delocalized state, the momentum distribution indicates a freezing of one-particle correlations at higher values than in thermal equilibrium. The position of the delocalization to localization transition, identified by the breakdown of thermalization with increasing disorder strength, is found to be consistent with the value from the level statistics obtained via full exact diagonalization of finite chains. Our results strongly support the existence of a many-body localized phase in the thermodynamic limit.We use thermalization indicators and numerical linked cluster expansions to probe the onset of many-body localization in a disordered one-dimensional hard-core boson model in the thermodynamic limit. We show that after equilibration following a quench from a delocalized state, the momentum distribution indicates a freezing of one-particle correlations at higher values than in thermal equilibrium. The position of the delocalization to localization transition, identified by the breakdown of thermalization with increasing disorder strength, is found to be consistent with the value from the level statistics obtained via full exact diagonalization of finite chains. Our results strongly support the existence of a many-body localized phase in the thermodynamic limit.

Nov 3, 20140 citations

Many-body localization and delocalization dynamics in the thermodynamic limit

Jonas Richter, Arijeet Pal

Disordered quantum systems undergoing a many-body localization (MBL) transition fail to reach thermal equilibrium under their own dynamics. Distinguishing between asymptotically localized or delocalized dynamics based on numerical results is, however, nontrivial due to finite-size effects. Numerical linked cluster expansions (NLCE) provide a means to tackle quantum systems directly in the thermodynamic limit but are challenging for models without translational invariance. Here we demonstrate that NLCE provides a powerful tool to explore MBL by simulating quench dynamics in disordered spin- $1/2$ two-leg ladders and Fermi-Hubbard chains. Combining NLCE with an efficient real-time evolution of pure states, we obtain converged results for the decay of the imbalance on long time scales and show that, especially for intermediate disorder below the putative MBL transition, NLCE outperforms direct simulations of finite systems with open or periodic boundaries. Furthermore, while spin is delocalized even in strongly disordered Hubbard chains with frozen charge, we unveil that an additional tilted potential leads to a drastic slowdown of the spin imbalance and nonergodic behavior on accessible times. Our work sheds light on MBL in systems beyond the well-studied disordered Heisenberg chain and emphasizes the usefulness of NLCE for this purpose.Disordered quantum systems undergoing a many-body localization (MBL) transition fail to reach thermal equilibrium under their own dynamics. Distinguishing between asymptotically localized or delocalized dynamics based on numerical results is, however, nontrivial due to finite-size effects. Numerical linked cluster expansions (NLCE) provide a means to tackle quantum systems directly in the thermodynamic limit but are challenging for models without translational invariance. Here we demonstrate that NLCE provides a powerful tool to explore MBL by simulating quench dynamics in disordered spin-$1/2$ two-leg ladders and Fermi-Hubbard chains. Combining NLCE with an efficient real-time evolution of pure states, we obtain converged results for the decay of the imbalance on long time scales and show that, especially for intermediate disorder below the putative MBL transition, NLCE outperforms direct simulations of finite systems with open or periodic boundaries. Furthermore, while spin is delocalized even in strongly disordered Hubbard chains with frozen charge, we unveil that an additional tilted potential leads to a drastic slowdown of the spin imbalance and nonergodic behavior on accessible times. Our work sheds light on MBL in systems beyond the well-studied disordered Heisenberg chain and emphasizes the usefulness of NLCE for this purpose.

arXiv:2202.10498Physical review. B./Physical review. BJun 24, 20225 citations

Quantum Quenches and Relaxation Dynamics in the Thermodynamic Limit

Krishnanand Mallayya, Marcos Rigol

We implement numerical linked cluster expansions (NLCEs) to study dynamics of lattice systems following quantum quenches, and focus on a hard-core boson model in one-dimensional lattices. We find that, in the nonintegrable regime and within the accessible times, local observables exhibit exponential relaxation. We determine the relaxation rate as one departs from the integrable point and show that it scales quadratically with the strength of the integrability breaking perturbation. We compare the NLCE results with those from exact diagonalization calculations on finite chains with periodic boundary conditions, and show that NLCEs are far more accurate.We implement numerical linked cluster expansions (NLCEs) to study dynamics of lattice systems following quantum quenches, and focus on a hard-core boson model in one-dimensional lattices. We find that, in the nonintegrable regime and within the accessible times, local observables exhibit exponential relaxation. We determine the relaxation rate as one departs from the integrable point and show that it scales quadratically with the strength of the integrability breaking perturbation. We compare the NLCE results with those from exact diagonalization calculations on finite chains with periodic boundary conditions, and show that NLCEs are far more accurate.

arXiv:1711.01284Physical Review LettersFeb 16, 201854 citations

Early Breakdown of Area-Law Entanglement at the Many-Body Delocalization Transition

Trithep Devakul, Rajiv R. P. Singh

We introduce the numerical linked cluster expansion as a controlled numerical tool for the study of the many-body localization transition in a disordered system with continuous nonperturbative disorder. Our approach works directly in the thermodynamic limit, in any spatial dimension, and does not rely on any finite size scaling procedure. We study the onset of many-body delocalization through the breakdown of area-law entanglement in a generic many-body eigenstate. By looking for initial signs of an instability of the localized phase, we obtain a value for the critical disorder, which we believe should be a lower bound for the true value, that is higher than current best estimates from finite size studies. This implies that most current methods tend to overestimate the extent of the localized phase due to finite size effects making the localized phase appear stable at small length scales. We also study the mobility edge in these systems as a function of energy density, and we find that our conclusion is the same at all examined energies.We introduce the numerical linked cluster expansion as a controlled numerical tool for the study of the many-body localization transition in a disordered system with continuous nonperturbative disorder. Our approach works directly in the thermodynamic limit, in any spatial dimension, and does not rely on any finite size scaling procedure. We study the onset of many-body delocalization through the breakdown of area-law entanglement in a generic many-body eigenstate. By looking for initial signs of an instability of the localized phase, we obtain a value for the critical disorder, which we believe should be a lower bound for the true value, that is higher than current best estimates from finite size studies. This implies that most current methods tend to overestimate the extent of the localized phase due to finite size effects making the localized phase appear stable at small length scales. We also study the mobility edge in these systems as a function of energy density, and we find that our conclusion is the same at all examined energies.

arXiv:1508.04813Physical Review LettersOct 28, 2015158 citations

Expansion of one-dimensional lattice hard-core bosons at finite temperature

Wei Xu, Marcos Rigol

We develop an exact approach to study the quench dynamics of hard-core bosons initially in thermal equilibrium in one-dimensional lattices. This approach is used to study the sudden expansion of thermal states after confining potentials are switched off. We find that a dynamical fermionization of the momentum distribution occurs at all temperatures. This phenomenon is studied for low initial site occupations, for which the expansion of the cloud is self-similar. In this regime, the occupation of the natural orbitals allows one to distinguish hard-core bosons from noninteracting fermions. We also study the free expansion of initial Mott insulating domains at finite temperature and show that the emergence of off-diagonal one-body correlations is suppressed gradually with increasing temperature. Surprisingly, the melting of the Mott domain is accompanied by an effective cooling of the system. We explain this phenomenon analytically using an equilibrium description based on an emergent local Hamiltonian.We develop an exact approach to study the quench dynamics of hard-core bosons initially in thermal equilibrium in one-dimensional lattices. This approach is used to study the sudden expansion of thermal states after confining potentials are switched off. We find that a dynamical fermionization of the momentum distribution occurs at all temperatures. This phenomenon is studied for low initial site occupations, for which the expansion of the cloud is self-similar. In this regime, the occupation of the natural orbitals allows one to distinguish hard-core bosons from noninteracting fermions. We also study the free expansion of initial Mott insulating domains at finite temperature and show that the emergence of off-diagonal one-body correlations is suppressed gradually with increasing temperature. Surprisingly, the melting of the Mott domain is accompanied by an effective cooling of the system. We explain this phenomenon analytically using an equilibrium description based on an emergent local Hamiltonian.

arXiv:1612.08988Physical Review AMar 16, 201734 citations

Slow Many-Body Delocalization beyond One Dimension

Elmer V. H. Doggen, I. V. Gornyi, A. D. Mirlin +1 more

We study the delocalization dynamics of interacting disordered hard-core bosons for quasi-1D and 2D geometries, with system sizes and timescales comparable to state-of-the-art experiments. The results are strikingly similar to the 1D case, with slow, subdiffusive dynamics featuring power-law decay. From the freezing of this decay we infer the critical disorder W_{c}(L,d) as a function of length L and width d. In the quasi-1D case W_{c} has a finite large-L limit at fixed d, which increases strongly with d. In the 2D case W_{c}(L,L) grows with L. The results are consistent with the avalanche picture of the many-body localization transition.We study the delocalization dynamics of interacting disordered hard-core bosons for quasi-1D and 2D geometries, with system sizes and timescales comparable to state-of-the-art experiments. The results are strikingly similar to the 1D case, with slow, subdiffusive dynamics featuring power-law decay. From the freezing of this decay we infer the critical disorder W_{c}(L,d) as a function of length L and width d. In the quasi-1D case W_{c} has a finite large-L limit at fixed d, which increases strongly with d. In the 2D case W_{c}(L,L) grows with L. The results are consistent with the avalanche picture of the many-body localization transition.

arXiv:2002.07635Physical Review LettersOct 5, 202057 citations

1/2

two-leg ladders and Fermi-Hubbard chains. Combining NLCE with an efficient real-time evolution of pure states, we obtain converged results for the decay of the imbalance on long time scales and show that, especially for intermediate disorder below the putative MBL transition, NLCE outperforms direct simulations of finite systems with open or periodic boundaries. Furthermore, while spin is delocalized even in strongly disordered Hubbard chains with frozen charge, we unveil that an additional tilted potential leads to a drastic slowdown of the spin imbalance and nonergodic behavior on accessible times. Our work sheds light on MBL in systems beyond the well-studied disordered Heisenberg chain and emphasizes the usefulness of NLCE for this purpose.