Exploring one-particle orbitals in large many-body localized systems

Strong disorder in interacting quantum systems can give rise to the phenomenon of many-body localization (MBL), which defies thermalization due to the formation of an extensive number of quasilocal integrals of motion. The one-particle operator content of these integrals of motion is related to the one-particle orbitals (OPOs) of the one-particle density matrix and shows a strong signature across the MBL transition as recently pointed out by Bera et al. [Phys. Rev. Lett. 115, 046603 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.046603; Ann. Phys. 529, 1600356 (2017)ANPYA20003-380410.1002/andp.201600356]. We study the properties of the OPOs of many-body eigenstates of an MBL system in one dimension. Using shift-and-invert MPS, a matrix product state method to target highly excited many-body eigenstates introduced previously [Phys. Rev. Lett. 118, 017201 (2017)PRLTAO0031-900710.1103/PhysRevLett.118.017201], we are able to obtain accurate results for large systems of sizes up to L=64. We find that the OPOs drawn from eigenstates at different energy densities have high overlap and their occupations are correlated with the energy of the eigenstates. Moreover, the standard deviation of the inverse participation ratio of these orbitals is maximal at the nose of the mobility edge. Also, the OPOs decay exponentially in real space, with a correlation length that increases at low disorder. In addition, we find that the probability distribution of the strength of the large-range coupling constants of the number operators generated by the OPOs approach a log-uniform distribution at strong disorder.