
Wavefunction extreme values statistics in Anderson localization
P. R. N. Falc˜ao∗
Instituto de F´ısica, Universidade Federal de Alagoas, 57072-900 Macei´o, AL, Brazil and
Institute of Theoretical Physics, Jagiellonian University in Krakow, ul. Lojasiewicza 11, 30-348 Krak´ow, Poland
M. L. Lyra†
Instituto de F´ısica, Universidade Federal de Alagoas, 57072-900 Macei´o, AL, Brazil
We consider a disordered one-dimensional tight-binding model with power-law decaying hopping
amplitudes to disclose wavefunction maximum distributions related to the Anderson localization
phenomenon. Deeply in the regime of extended states, the wavefunction intensities follow the
Porter-Thomas distribution while their maxima assume the Gumbel distribution. At the critical
point, distinct scaling laws govern the regimes of small and large wavefunction intensities with a
multifractal singularity spectrum. The distribution of maxima deviates from the usual Gumbel
form and some characteristic finite-size scaling exponents are reported. Well within the localization
regime, the wavefunction intensity distribution is shown to develop a sequence of pre-power-law,
power-law, exponential and anomalous localized regimes. Their values are strongly correlated,
which significantly affects the emerging extreme values distribution.
I. INTRODUCTION
The Anderson localization transition is a key physical
phenomenon related to a drastic change in the spatial
distribution of Hamiltonian eigenfunctions promoted by
disorder [1–3]. In general lines, the eigenfunctions be-
come exponentially localized in the regime of strong dis-
order while remaining spatially extended for weak disor-
der. Anderson transition is a quite general phenomenon
affecting the transport of electronic, acoustic, magnetic
and optical waves [4–8].
The critical behavior of the Anderson transition has
been extensively studied over the past decades and shown
to strongly depend on the system’s dimensionality, sym-
metries and the short or long-range character of the un-
derlying couplings [3]. In particular, all eigenstates be-
come exponentially localized in one-dimensional systems
with short-range couplings for any finite amount of un-
correlated disorder. In contrast, one-dimensional systems
with power-law decaying couplings support an Anderson
transition and have been frequently used as a prototype
model to investigate its universality classes and critical
behavior [9–24].
Field theoretical renormalization group and random
matrix (RM) theories yield several relevant aspects of the
Anderson localization transition. In particular, random
matrix theory unveiled universal characteristics of the
eigenvalues statistics in the localized and extended phases
[25]. Exponentially localized states are uncorrelated and
randomly distributed along the energy band, resulting
in a Poissonian probability distribution function (PDF)
of the level spacements. On the other hand, the level
repulsion typical of spatially extended states leads to a
new spacement PDF that depends of the nature of the
∗Electronic address: pedro.falcao@fis.ufal.br
†Electronic address: marcelo@fis.ufal.br
random matrix ensemble (Gaussian orthogonal, unitary
or symplectic). RM theory has also been used to explore
the statistical properties of the extreme eigenvalues [26].
Extreme events play a fundamental role in the study
of disordered physical systems [27,28]. These are events
with magnitude much larger than the average value
of a physical stochastic process. In condensed matter
physics, they have a key impact on transport phenomena
where the largest energy barriers and maximally localized
modes hinder the transmission of physical excitations
[29,30]. The extreme values statistics of uncorrelated
and identically distributed (IID) random variables is well
understood [31]. Unbounded random variables having a
PDF with a faster than power-law tail have a Gumbel dis-
tribution of the extreme values in long sub-sequences. A
Fr´echet extreme values PDF is in order for random vari-
ables with power-law decaying PDF tails. In the case of
bounded stochastic variables, the extreme values are dis-
tributed according to an asymptotic Weilbull PDF. For
sequences of strongly correlated random variables, very
little is known regarding the extreme values statistics (for
a recent review see [28]). In the context of Anderson lo-
calization, although RM theory discloses the PDF of the
extreme values for the correlated eigenvalues of Gaussian
ensembles [26], studies of the extreme values statistics of
the own eigenfunction intensities are missing. Although
the structureless nature of Gaussian extended eigenfunc-
tions with an exponentially decaying Porter-Thomas in-
tensity distribution[32] allows us to anticipate that the
maximum distribution of extended states shall fall in the
Gumbel class, multifractal correlations present in critical
states as well as the structured aspect of exponentially
localized states point to new distributions of the wave-
function intensity maxima that are still unexplored.
In this work, we address the above relevant open ques-
tion aiming to unveil how the extended or localized na-
ture of single-particle eigenstates in a disordered sys-
tem is reflected in the extreme values statistics of their
eigenfunctions. We will consider the prototype one-
arXiv:2210.06365v1 [cond-mat.dis-nn] 12 Oct 2022