We study one-dimensional spinless fermions with random interactions, but without any on-site disorder. We find that random interactions generically stabilize a many-body localized phase, in spite of the completely extended single-particle degrees of freedom. In the large randomness limit, we construct ``bubble-neck'' eigenstates having a universal area-law entanglement entropy on average, with the number of volume-law states being exponentially suppressed. We argue that this statistical localization is beyond the phenomenological local-integrals-of-motion description of many-body localization. With exact diagonalization, we confirm the robustness of the many-body localized phase at finite randomness by investigating eigenstate properties such as level statistics, entanglement/participation entropies, and nonergodic quantum dynamics. At weak random interactions, the system develops a thermalization transition when the single-particle hopping becomes dominant.
In the presence of strong enough disorder one-dimensional systems of interacting spinless fermions at nonzero filling factor are known to be in a many-body localized phase. When represented in `Fock space,' the Hamiltonian of such a system looks like that of a single `particle' hopping on a Fock lattice in the presence of a random disordered potential. The coordination number of the Fock lattice increases linearly with the system size
L in one dimension. Thus in the thermodynamic limit $L\ensuremath{\rightarrow}\ensuremath{\infty}$, the disordered interacting problem in one dimension maps on to an Anderson model with infinite coordination number. Despite this, this system displays localization which appears counterintuitive. A close observation of the on-site disorder potentials on the Fock lattice reveals a large degree of correlation among them as they are derived from an exponentially smaller number of on-site disorder potentials in real space. This indicates that the correlations between the on-site disorder potentials on a Fock lattice has a strong effect on the localization properties of the corresponding many-body system. This intuition is also consistent with studies of quantum random energy model where the typical mid-spectrum states are ergodic and the on-site potentials in Fock space are completely uncorrelated. In this work we perform a systematic quantitative exploration of the nature of correlations of the Fock space potential required for localization. We study different functional variations of the disorder correlation in Fock lattice by analyzing the eigenspectrum obtained through exact diagonalization. Without changing the typical strength of the on-site disorder potential in Fock lattice we show that changing the correlation strength can induce thermalization or localization in systems. From among the various forms of correlations we study, we find that only the linear variation of correlations with Hamming distance in Fock space is able to drive a thermal-MBL phase transition where the transition is driven by the correlation strength. Systems with the other forms of correlations we study are found to be ergodic.
We study many-body localization (MBL) in a one-dimensional system of spinless fermions with a deterministic aperiodic potential in the presence of random interactions
Vij decaying as power-law ${V}_{ij}/{({r}_{ij})}^{\ensuremath{\alpha}}$ with distance
rij. We demonstrate that MBL survives even for $\ensuremath{\alpha}<1$ and is preceded by a broad nonergodic subdiffusive phase. Starting from parameters at which the short-range interacting system shows an infinite temperature MBL phase, turning on random power-law interactions results in many-body mobility edges in the spectrum with a larger fraction of ergodic delocalized states for smaller values of $\ensuremath{\alpha}$. Hence, the critical disorder
hcr, at which ergodic to nonergodic transition takes place, increases with the range of interactions. Time evolution of the density imbalance
I(t), which has power-law decay $I(t)\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\gamma}}$ in the intermediate to large time regime, shows that the critical disorder
hcI, above which the system becomes diffusionless (with $\ensuremath{\gamma}\ensuremath{\sim}0$) and transits into the MBL phase, is much larger than
hcr. In between
hcr and
hcI there is a broad nonergodic subdiffusive phase, which is characterized by the Poissonian statistics for the level spacing ratio, multifractal eigenfunctions, and a nonzero dynamical exponent $\ensuremath{\gamma}\ensuremath{\ll}1/2$. The system continues to be subdiffusive even on the ergodic side (
h<hcr) of the MBL transition, where the eigenstates near the mobility edges are multifractal. For
h<h0<hcr, the system is superdiffusive with $\ensuremath{\gamma}>1/2$. The rich phase diagram obtained here is unique to the random nature of long-range interactions. We explain this in terms of the enhanced correlations among local energies of the effective Anderson model induced by random power-law interactions.
Recent work by De Roeck et al. [Phys. Rev. B 95, 155129 (2017)] has argued that many-body localization is unstable in two and higher dimensions due to a thermalization avalanche triggered by rare regions of weak disorder. To examine these arguments, we construct several models of a finite ergodic bubble coupled to an Anderson insulator of noninteracting fermions. We first describe the ergodic region using a Gaussian orthogonal ensemble random matrix and perform an exact diagonalization study of small systems. The results are in excellent agreement with a refined theory of the thermalization avalanche that includes transient finite-size effects, lending strong support to the avalanche scenario. We then explore the limit of large system sizes by modeling the ergodic region via a Hubbard model with all-to-all random hopping: the combined system, consisting of the bubble and the insulator, can be reduced to an effective Anderson impurity problem. We find that the spectral function of a local operator in the ergodic region changes dramatically when coupling to a large number of localized fermionic states; this occurs even when the localized sites are weakly coupled to the bubble. In principle, for a given size of the ergodic region, this may arrest the avalanche. However, this back-action effect is suppressed and the avalanche can be recovered if the ergodic bubble is large enough. Thus, the main effect of the back-action is to renormalize the critical bubble size.
Strong long-range interaction leads to localization in a closed quantum system without disorders. Employing the exact diagonalization method, I numerically investigate thermalization and many-body localization in translational invariant quantum chains with finite Coulomb interactions. In the computational basis, excluding all trivial degeneracies, the interaction-induced localization is well demonstrated in aspects of level statistics, eigenstate expectation values, and Anderson localization on graphs constructed of many-body bases. The nature of localization for generic eigenstates is attributed to the quasidisorder from the power-law interactions, and the full localization in the Hilbert space is similar to that in the disorder case. However, due to real-space symmetries, the long-time dynamics is dominated by the degenerated eigenstates and eventually reaches homogeneity in real space. On the other hand, the entanglement entropy exhibits size dependence beyond the area law for the same reason, even deep in the localized state, indicating an incomplete localization in real space.
Many-body localized (MBL) systems lie outside the framework of statistical mechanics, as they fail to equilibrate under their own quantum dynamics. Even basic features of MBL systems, such as their stability to thermal inclusions and the nature of the dynamical transition to thermalizing behaviour, remain poorly understood. We study a simple central spin model to address these questions: a two-level system interacting with strength J with N ≫1 localized bits subject to random fields. On increasing J , the system transitions from an MBL to a delocalized phase on the vanishing scale J c ( N )∼1/ N , up to logarithmic corrections. In the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either ‘on’ (strongly entangled) or ‘off’ (weakly entangled) in eigenstates. The clusters of ‘on’ bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogeneous discontinuous transition out of the localized phase in single-site observables. We obtain these results by perturbative mapping to bond percolation on the hypercube at small J and by numerical exact diagonalization of the full many-body system. Our results support the arguments that the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions d that are high but finite. This article is part of the themed issue ‘Breakdown of ergodicity in quantum systems: from solids to synthetic matter’.
Understanding the nature of the transition from the delocalized to the many-body localized (MBL) phase is an important unresolved issue. To probe the nature of the MBL transition, we investigate the universal properties of single-particle excitations produced in highly excited many-body eigenstates of a disordered interacting quantum many-body system. In a class of one-dimensional spinless fermionic models with random disorder, we study the finite-size scaling of the ratio of typical to average values of the single-particle local density of states and the scattering rates across the MBL transition. Our results indicate that the MBL transition in this class of one-dimensional models of spinless fermions is continuous in nature. For various ranges of interactions in the system, the critical exponent $\ensuremath{\nu}$ with which the correlation length $\ensuremath{\xi}$ diverges at the transition point ${W}_{c}, \ensuremath{\xi}\ensuremath{\sim}|W\ensuremath{-}{W}_{c}{|}^{\ensuremath{-}\ensuremath{\nu}}$, satisfies the Chayes-Chayes-Fisher-Spencer (CCFS) bound $\ensuremath{\nu}\ensuremath{\ge}2/d$ where
d is the physical dimension of the system. We also discuss why the critical exponent obtained from finite-size scaling of the conventional diagnostic of many-body localization, the level spacing ratio, strongly violates the CCFS bound while the single-particle density of states and scattering rates are consistent with the CCFS criterion.