Localization transition in random Lévy matrices: multifractality of eigenvectors in the localized phase and at criticality

For random Lévy matrices of size $N\times N$, where matrix elements are drawn with some heavy-tailed distribution $P\left(H_{ij}\right)\propto N^{-1}|H_{ij}{|}^{-1-\mu }$ with $0<\mu <2$ (infinite variance), there exists an extensive number of finite eigenvalues E  =  O(1), while the maximal eigenvalue grows as $E_{\text{max}}\sim N^{\frac{1}{\mu }}$. Here we study the localization properties of the corresponding eigenvectors via some strong disorder perturbative expansion that remains consistent within the localized phase and that yields their inverse participation ratios (IPR) Y _q as a function of the continuous parameter $0<q<+\infty$. In the region $0<\mu <1$, we find that all eigenvectors are localized but display some multifractality: the IPR are finite above some threshold q  >  q _c but diverge in the region 0  <  q  <  q _c near the origin. In the region $1<\mu <2$, only the sub-extensive fraction $N^{\frac{3}{2+\mu }}$ of the biggest eigenvalues corresponding to the region $|E|⩾N^{\frac{(\mu -1)}{\mu (2+\mu )}}$ remains localized, while the extensive number of other states of smaller energy are delocalized. For the extensive number of finite eigenvalues E  =  O(1), the localization/delocalization transition thus takes place at the critical value ${\mu }_{\text{c}}=1$ corresponding to Cauchy matrices: the IPR Y _q of the corresponding critical eigenstates follow the strong-multifractality spectrum characterized by the generalized fractal dimensions $D^{\text{criti}}\left(q\right)=\frac{1-2q}{1-q}\theta (0⩽q⩽\frac{1}{2})$, which has been found previously in various other Localization problems in spaces of effective infinite dimensionality. (paper: disordered systems, classical and quantum)