Delocalization and quantum diffusion of random band matrices in high dimensions I: Self-energy renormalization
Fan Yang, Horng-Tzer Yau, Jun Yin
TL;DR
The authors address delocalization and quantum diffusion for high-dimensional random band matrices by developing a rigorous, high-order T-expansion with self-energies. Their method translates resolvent behavior into a diffusive kernel through a renormalized diffusion matrix $\Theta$ and systematic self-energy corrections, enabling precise local laws and diffusion-type asymptotics for $\mathbb{E}|G_{xy}(z)|^2$. The framework, built from detailed graph operations and a doubly connected structure, yields delocalization results for bulk eigenvectors in $d\ge 8$ under a mild band-width condition and confirms a quadratic-in-$p$ diffusion form for the Fourier transform of the Green’s function squared. The approach is robust to Gaussian assumptions and is expected to extend to non-Gaussian band matrices, offering a powerful toolkit for non-mean-field disordered systems and connecting spectral properties to quantum diffusion on finite lattices.
Abstract
We consider Hermitian random band matrices $H=(h_{xy})$ on the $d$-dimensional lattice $(\mathbb Z/L\mathbb Z)^d$. The entries $h_{xy}$ are independent (up to Hermitian conditions) centered complex Gaussian random variables with variances $s_{xy}=\mathbb E|h_{xy}|^2$. The variance matrix $S=(s_{xy})$ has a banded structure so that $s_{xy}$ is negligible if $|x-y|$ exceeds the band width $W$. In dimensions $d\ge 8$, we prove that, as long as $W\ge L^ε$ for a small constant $ε>0$, with high probability most bulk eigenvectors of $H$ are delocalized in the sense that their localization lengths are comparable to $L$. Denote by $G(z)=(H-z)^{-1}$ the Green's function of the band matrix. For ${\mathrm Im}\, z\gg W^2/L^2$, we also prove a widely used criterion in physics for quantum diffusion of this model, namely, the leading term in the Fourier transform of $\mathbb E|G_{xy}(z)|^2$ with respect to $x-y$ is of the form $({\mathrm Im}\, z + a(p))^{-1}$ for some $a(p)$ quadratic in $p$, where $p$ is the Fourier variable. Our method is based on an expansion of $T_{xy}=|m|^2 \sum_αs_{xα}|G_{αy}|^2$ and it requires a self-energy renormalization up to error $W^{-K}$ for any large constant $K$ independent of $W$ and $L$. We expect that this method can be extended to non-Gaussian band matrices.
