Rosenzweig-Porter Random Matrix Ensemble

• The Rosenzweig-Porter Random Matrix Ensemble is a one-parameter family characterized by ergodic, fractal, and localized phases controlled by the scaling parameter γ.
• Analytical techniques such as the supersymmetry formalism and functional Hubbard–Stratonovich transformation yield self-consistent equations for spectral density and transition order parameters.
• A rich (μ, γ) phase diagram illustrates nontrivial multifractality and ergodicity breaking, offering insights into localization phenomena in disordered quantum systems.

The Rosenzweig–Porter (RP) Random Matrix Ensemble is a one-parameter family of large random matrices that realizes distinct ergodic, fractal (non-ergodic extended), and localized phases of eigenstates as a function of the scaling parameter $\gamma$. The Lévy Rosenzweig–Porter ensemble, a fat-tailed extension, provides a framework for analytically tractable non-ergodic extended phases in systems where off-diagonal matrix elements are $\mu$-stable Lévy random variables. The model evidences sharp ergodic–fractal and fractal–localized transitions, described by nontrivial spectral and multifractal properties, and exhibits a rich $\mu$–$\gamma$ phase diagram. The analytical study of its spectral density relies on functional Hubbard–Stratonovich transformations combined with Efetov’s supersymmetry formalism, yielding a self-consistent order-parameter-driven description of the ergodic–fractal transition (Safonova et al., 2024).

1. Construction and Definition of the Lévy Rosenzweig–Porter Ensemble

The Lévy RP ensemble is defined by the Hamiltonian $H = H^{(1)} + H^{(2)}$, where

• $H^{(1)}$ is diagonal, with elements $H^{(1)}_{kk}$ drawn independently from a Gaussian distribution $$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$,
• $H^{(2)}$ has full off-diagonal structure, $H^{(2)}_{ij}$ ($\mu$0) independently drawn from a symmetric $\mu$1-stable Lévy law $\mu$2 with characteristic function $\mu$3, $\mu$4.

Off-diagonal elements are rescaled as

$\mu$5

so that their typical scale is $\mu$6, with $\mu$7 tuning the hopping-to-diagonal-disorder ratio. For $\mu$8, $\mu$9 reduces to a Gaussian and the model reduces to the classical (Gaussian) RP ensemble; for $\mu$0, the power-law tails render the second moment of off-diagonal elements divergent.

2. Analytical Framework: Supersymmetry and Functional HS Transformation

The non-Gaussian, power-law nature of the off-diagonal entries invalidates the standard (quadratic) Hubbard–Stratonovich transformation typically used for Gaussian RP models. Instead, the analytic approach proceeds as follows (Safonova et al., 2024):

• The partition function in the Efetov supervector representation is

$\mu$1

where $\mu$2 combines bosonic and Grassmann components.

• Averaging over $\mu$3, the cumulant generates a non-analytic $\mu$4 inter-site action term:

$\mu$5

• The required functional HS transformation linearizes the non-analytic action:

$\mu$6

and the large-$\mu$7 saddle point yields a functional self-consistency equation for $\mu$8.

3. Density of States and Scaling as an Order Parameter

The mean spectral density $\mu$9 is controlled by a nontrivial transcendental equation for $\gamma$0,

$\gamma$1

with $\gamma$2 the Euler beta function. Then,

$\gamma$3

For Gaussian $\gamma$4, the logarithmic term simplifies: $\gamma$5.

Rescaling to $\gamma$6 with $\gamma$7 shows that all $\gamma$8 and $\gamma$9 dependence enters through the single variable $H = H^{(1)} + H^{(2)}$0. In this way, the spectral density acts as a "spectral order parameter" distinguishing the ergodic and fractal regimes:

• $H = H^{(1)} + H^{(2)}$1 for $H = H^{(1)} + H^{(2)}$2 (ergodic phase),
• $H = H^{(1)} + H^{(2)}$3 for $H = H^{(1)} + H^{(2)}$4 (fractal phase).

The critical point $H = H^{(1)} + H^{(2)}$5 marks a nontrivial, $H = H^{(1)} + H^{(2)}$6-independent line; finite-size scaling for the order parameter yields a scaling exponent $H = H^{(1)} + H^{(2)}$7.

4. Phase Diagram and Multifractality

The Lévy-RP model displays a phase diagram governed jointly by $H = H^{(1)} + H^{(2)}$8 and $H = H^{(1)} + H^{(2)}$9, with three principal regimes (Biroli et al., 2020):

• For $H^{(1)}$0, eigenstates are fully delocalized (ergodic): level statistics are Wigner–Dyson and spectral density displays band spreading $H^{(1)}$1;
• For $H^{(1)}$2, eigenstates are non-ergodic extended (fractal or multifractal): they occupy $H^{(1)}$3 sites, with

$H^{(1)}$4

reflecting a nontrivial multifractal spectrum $H^{(1)}$5;

• For $H^{(1)}$6, eigenstates are Anderson localized.

At $H^{(1)}$7, a mobility edge appears, and at $H^{(1)}$8, $H^{(1)}$9 is a tricritical point. For $H^{(1)}_{kk}$0 (Lévy tails), genuine multifractality emerges; for $H^{(1)}_{kk}$1, fractal dimensions coincide and the classical Gaussian RP scaling $H^{(1)}_{kk}$2 is recovered.

5. Spectral, Dynamical, and Correlation Properties

Physical signatures of the three regimes are manifest in both spectral correlations and dynamics (Lunkin et al., 2024):

• Local density of states (LDOS) correlations $H^{(1)}_{kk}$3 exhibit universal scaling:

$H^{(1)}_{kk}$4

with $H^{(1)}_{kk}$5 the characteristic width of a local state.

• Return probability $H^{(1)}_{kk}$6 at long times displays a stretched-exponential decay governed by the Lévy index,

$H^{(1)}_{kk}$7

For $H^{(1)}_{kk}$8 (Gaussian RP), this becomes a simple exponential, $H^{(1)}_{kk}$9; for $$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$0, the stretched exponential reflects non-Fermi's Golden Rule (FGR) dynamics.

• Level compressibility and long-range spectral correlations at the Thouless scale are universal across Gaussian, Lévy, and Wishart RP variants.

6. Physical Interpretation and Applications

The Lévy RP ensemble encapsulates a minimal yet analytically tractable model for the breakdown of ergodicity via fractal or multifractal mini-bands, where non-ergodic extended states retain overlap on a sub-extensive number of sites, yielding Wigner–Dyson statistics on the mini-band scale but anomalous global features (Biroli et al., 2020). This ensemble provides a robust synthetic benchmark for non-ergodic extended behavior in many-body disordered systems, including analogs of "bad metallicity," sub-diffusive transport, and multifractal eigenstates observed in Hilbert space for interacting quantum systems.

Numerical studies confirm the predictions for $$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$1 and the scaling of LDOS correlations up to system sizes $$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$2 (Safonova et al., 2024). The established connection between the universal spectral order parameter and transitions in eigenfunction statistics (ergodicity breaking at $$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$3, localization at $$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$4) provides analytic tools for probing nontrivial intermediate regimes in both single-particle and many-body localization contexts.

7. Summary of Key Scaling Relations and Table of Phase Diagram

Below, the (μ, γ) phase diagram is summarized:

Regime	Condition	Eigenfunction Scaling	Level Statistics
Ergodic	$$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$5	$$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$6	Wigner–Dyson
Fractal/Multifractal	$$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$7	$$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$8	WD on mini-band, Poisson outside
Localized	$$P_1(H^{(1)}_{kk}) = (2\pi W^2)^{-1/2}\exp[-(H^{(1)}_{kk})^2/2W^2]$$9	$H^{(2)}$0	Poisson

The critical lines at $H^{(2)}$1 (ergodic–fractal transition) and $H^{(2)}$2 (fractal–localized transition) are validated analytically and numerically, with a universal finite-size scaling exponent $H^{(2)}$3 at the ergodic–fractal line. The order parameter $H^{(2)}$4 transitions smoothly across the phase boundary (Safonova et al., 2024).

The Lévy Rosenzweig–Porter ensemble thus stands as a paradigmatic, analytically tractable model for exploring ergodicity breaking, multifractality, and localization transitions in strongly disordered quantum systems.