An analytical theory, based on the perturbative treatment of the disorder and extended into a self-consistent set of equations for the dynamical density correlations, is developed and applied to the prototype one-dimensional model of many-body localization. Results show a qualitative agreement with the numerically obtained dynamical structure factor in the whole range of frequencies and wave vectors, as well as across the transition to nonergodic behavior. The theory reveals the singular nature of the one-dimensional problem, whereby on the ergodic side the dynamics is subdiffusive with dynamical conductivity $\ensuremath{\sigma}(\ensuremath{\omega})\ensuremath{\propto}{|\ensuremath{\omega}|}^{\ensuremath{\alpha}}$, i.e., with vanishing dc limit ${\ensuremath{\sigma}}_{0}=0$ and $\ensuremath{\alpha}<1$ varying with disorder, while we get $\ensuremath{\alpha}>1$ in the localized phase.
An analytical theory, based on the perturbative treatment of the disorder and extended into a self-consistent set of equations for the dynamical density correlations, is developed and applied to the prototype one-dimensional model of many-body localization. Results show a qualitative agreement with numerically obtained dynamical structure factor in the whole range of frequencies and wavevectors, as well as across the transition to the nonergodic behavior. The theory reveals the singular nature of the one-dimensional problem, whereby on the ergodic side the dynamics is subdiffusive with a dynamical exponent varying with the disorder.
We study the infinite temperature dynamics of a prototypical one-dimensional system expected to exhibit many-body localization. Using numerically exact methods, we establish the dynamical phase diagram of this system based on the statistics of its eigenvalues and its dynamical behavior. We show that the nonergodic phase is reentrant as a function of the interaction strength, illustrating that localization can be reinforced by sufficiently strong interactions even at infinite temperature. Surprisingly, within the accessible time range, the ergodic phase shows subdiffusive behavior, suggesting that the diffusion coefficient vanishes throughout much of the phase diagram in the thermodynamic limit. Our findings strongly suggest that Wigner-Dyson statistics of eigenvalue spacings may appear in a class of ergodic but subdiffusive systems.
We study the $t\text{\ensuremath{-}}V$ disordered spinless fermionic chain in the strong-coupling regime, $t/V\ensuremath{\rightarrow}0$. Strong interactions highly hinder the dynamics of the model, fragmenting its Hilbert space into exponentially many blocks in system size. Macroscopically, these blocks can be characterized by the number of new degrees of freedom, which we refer to as movers. We focus on two limiting cases: blocks with only one mover and ones with a finite density of movers. The former many-particle block can be exactly mapped to a single-particle Anderson model with correlated disorder in one dimension. As a result, these eigenstates are always localized for any finite amount of disorder. The blocks with a finite density of movers, on the other side, show a many-body localized (MBL) transition that is tuned by the disorder strength. Moreover, we provide numerical evidence that its ergodic phase is diffusive at weak disorder. Approaching the MBL transition, we observe subdiffusive dynamics at finite timescales and find indications that this might be only a transient behavior before crossing over to diffusion.
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.
Within one-dimensional disordered models of interacting fermions, we perform a numerical study of several dynamical density correlations, which can serve as hallmarks of the transition to the many-body localized state. The results confirm that density-wave correlations exhibit quite an abrupt change with increasing disorder, with a nonvanishing long-time value characteristic for the nonergodic phase. In addition, our results reveal a logarithmic variation of correlations in time in a wide time window, which we can bring in connection with the anomalous behavior of the dynamical conductivity near the transition. Our results support the view that the transition to many-body localization can be characterized by universal dynamical exponents.