We introduce the numerical linked cluster (NLC) expansion as a controlled numerical tool for the study of the many-body localization (MBL) transition in a disordered system with continuous non-perturbative 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 find that our conclusion is the same at all examined energies.
PACS numbers:Introduction-The eigenstate thermalization hypothesis (ETH) is a powerful statement relating observables of the high energy eigenstates of a quantum many-body system with their thermal expectation values [1,2]. However, this principle can be violated in certain systems with strong enough disorder, where even the high energy eigenstates possess only local entanglement [3,4]. Anderson localization is a one-body example of this. An area of key interest is how far this localization persists in a many-body state in presence of interactions [5]. At what point are interactions strong enough that the localization is destroyed and the system obeys ETH? This is the problem of many-body localization (MBL) transition, which is a topic of active research both theoretically and experimentally .The surge of interest in many-body localized systems has motivated many numerical studies. Most studies have focused on exact diagonalization or Lanczos methods which are able to address both sides of the transition in small systems [32][33][34][35][36][37][38][39][40][41][42][43][44][45]. However, since much about this phase transition is still not well understood, extension of finite size results to the thermodynamic limit can prove difficult. We would like to examine this phase transition using expansion methods, which provide an alternate way of addressing the thermodynamic limit. While standard perturbative series expansions are very powerful [46][47][48], they suffer from small energy denominators in models with continuous non-perturbative disorder. Thus, we turn to the numerical linked cluster (NLC) expansion [49][50][51], which does not suffer from this problem of small energy denominators.In this paper, we provide evidence that for a prototypical model of MBL, approaching the critical disorder from the localized side, the localized phase actually becomes unstable only at increasingly long length scales inaccessible to most numerical techniques. This implies that finite si...