Level statistics across the many-body localization transition

Level statistics of systems that undergo many-body localization transition are studied. An analysis of the gap ratio statistics from the perspective of inter- and intrasample randomness allows us to pin point differences between transitions in random and quasirandom disorder, showing the effects due to Griffiths rare events for the former case. It is argued that the transition in the case of random disorder exhibits universal features that are identified by constructing an appropriate model of intermediate spectral statistics which is a generalization of the family of short-range plasma models. The considered weighted short-range plasma model yields a very good agreement both for level spacing distribution including its exponential tail and the number variance up to tens of level spacings outperforming previously proposed models. In particular, our model grasps the critical level statistics which arise at disorder strength for which the intersample fluctuations are the strongest. Going beyond the paradigmatic examples of many-body localization in spin systems, we show that the considered model also grasps the level statistics of disordered Bose- and Fermi-Hubbard models. The remaining deviations for long-range spectral correlations are discussed and attributed mainly to the intricacies of level unfolding.