Moment Lyapunov Exponents in Random Systems

• Moment Lyapunov exponents are quantitative measures that capture the exponential growth rate of moment norms in random dynamical processes.
• They are computed through spectral analysis of twisted transfer operators and cumulant expansion, providing detailed insights into stability and fluctuations.
• Applications span disordered spin chains, stochastic PDEs, and statistical mechanics, illustrating practical implications in stability and large deviation phenomena.

A moment Lyapunov exponent quantifies the exponential growth rate of the expectations of powers of norms of random dynamical processes, most commonly matrix products or stochastic flows. Formally, for i.i.d. random matrices $\{M_i\}$ in $\mathrm{SL}(2, \mathbb{R})$, and product $\Pi_n = M_n \cdots M_1$, the moment Lyapunov exponent is defined as $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$, extracting the large-$n$ growth of the $q$th moment independent of the initial direction $\mathbf{x}_0$. The function $\Lambda(q)$, often called the generalized or moment Lyapunov exponent, characterizes the full set of moment growth rates and encodes the cumulants of $\ln \|\Pi_n\|$ through its Taylor expansion about $q=0$. This concept generalizes to random dynamical systems on manifolds, products of transfer operators, and stochastic PDEs. Moment Lyapunov exponents provide a systematic hierarchy beyond the typical (quenched) Lyapunov exponent, with applications spanning localization, large deviations, intermittency, and stability in stochastic systems (Texier, 2019, Baxendale, 21 Jul 2025, Charbonneau et al., 2017, Bao et al., 2011, Borodin et al., 2012, Vladimirov, 2012, Duarte et al., 8 Sep 2025).

1. Formal Definitions and Properties

For i.i.d. random matrices $\mathrm{SL}(2, \mathbb{R})$0 in $\mathrm{SL}(2, \mathbb{R})$1, the product is $\mathrm{SL}(2, \mathbb{R})$2. The $\mathrm{SL}(2, \mathbb{R})$3th moment Lyapunov exponent is

$\mathrm{SL}(2, \mathbb{R})$4

independent of $\mathrm{SL}(2, \mathbb{R})$5 by strong irreducibility and subadditivity (Texier, 2019). Its expansion at $\mathrm{SL}(2, \mathbb{R})$6,

$\mathrm{SL}(2, \mathbb{R})$7

produces cumulants $\mathrm{SL}(2, \mathbb{R})$8 of $\mathrm{SL}(2, \mathbb{R})$9. For Markov processes, random SDE flows, or SPDEs, analogous asymptotics are

$\Pi_n = M_n \cdots M_1$0

for the derivative cocycle of a stochastic flow $\Pi_n = M_n \cdots M_1$1 (Baxendale, 21 Jul 2025), or

$\Pi_n = M_n \cdots M_1$2

for the $\Pi_n = M_n \cdots M_1$3th norm moment of SPDE solutions (Bao et al., 2011).

$\Pi_n = M_n \cdots M_1$4 is convex, $\Pi_n = M_n \cdots M_1$5, and completely characterizes the large deviations of empirical exponents via the Legendre transform $\Pi_n = M_n \cdots M_1$6 (Baxendale, 21 Jul 2025, Charbonneau et al., 2017). In random media, the Lyapunov exponent is related to integrated densities of states and exhibits subtle regularity (e.g., Hölder continuity demands strong moment constraints on the potential) (Duarte et al., 8 Sep 2025).

2. Spectral Characterizations via Transfer Operators

Central to the computation of $\Pi_n = M_n \cdots M_1$7 in the matrix product context is the spectral analysis of a family of twisted transfer operators. For $\Pi_n = M_n \cdots M_1$8 acting as a Möbius transformation, one introduces a transfer operator $\Pi_n = M_n \cdots M_1$9 acting on functions $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$0 by

$$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$1

where $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$2 is the projective Jacobian, and averages over $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$3 yield the operator $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$4 (Texier, 2019). The principal eigenvalue $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$5 of $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$6 obeys

$$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$7

with decay constraints on $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$8. The moment Lyapunov exponent is then $$\Lambda(q) = \lim_{n \to \infty} \frac{1}{n} \ln \mathbb{E}[\|\Pi_n \mathbf{x}_0\|^{q}]$$9. In the context of SDEs on manifolds, this spectral characterization translates to the twisted generator $n$0 acting on the sphere bundle, whose principal eigenvalue gives $n$1 (Baxendale, 21 Jul 2025).

For homogeneous Gaussian random fields in combinatorial/statistical mechanics, the MLE emerges as the spectral radius of a transfer operator $n$2 on a Hilbert space, reducing partition function asymptotics to eigenvalue computations (Vladimirov, 2012).

3. Calculation: Algorithms and Explicit Formulas

Several analytical and numerical approaches are developed for concrete computation:

• Random Matrix Products: For $n$3 products (including the Schrödinger equation with random potential), explicit spectral problems, often reducible in Fourier space, enable calculation of $n$4 and its derivatives as functionals of invariant densities (Texier, 2019).
• Cycle-Expansion (Ruelle Zeta Function): For one-dimensional matrix products, the Ruelle dynamical zeta function yields a systematic cycle-expansion in terms of prime and pseudocycles, producing cumulants $n$5 of finite-time Lyapunov exponents via analytic derivatives (Charbonneau et al., 2017).
• Gaussian Random Fields: In elastic monomer–dimer models, the growth rate of product moments reduces to the spectral radius of a compact transfer operator, tractable using matrix truncation or recursion in the orthonormal Hermite basis, and alternatively as solutions to “pantograph” functional-differential equations (Vladimirov, 2012).
• Stochastic PDEs with Markovian Switching: The $n$6-th moment Lyapunov exponent is given by

$n$7

where $n$8 is the Donsker–Varadhan rate functional for Markov chain occupation measures, and $n$9 encodes drift and noise coefficients (Bao et al., 2011).

• Parabolic Anderson Model: The $q$0-th moment Lyapunov exponent can be derived via multidimensional contour integrals or a variational principle involving the Bethe ansatz, producing explicit formulas for all $q$1 (Borodin et al., 2012).

4. Large Deviations and Fluctuations

The full function $q$2 controls the large deviations and fluctuations of empirical Lyapunov exponents:

• By the Gärtner–Ellis theorem, the rate function for the LDP of the empirical exponent $q$3 is $q$4 (Baxendale, 21 Jul 2025, Charbonneau et al., 2017).
• The second derivative $q$5 (variance of $q$6) determines the central limit theorem and moderate deviation regime.
• In high-dimensional or non-compact settings, growth conditions, hypoellipticity, and positivity of transition laws are necessary to ensure spectral existence, differentiability, and validity of fluctuation theorems (Baxendale, 21 Jul 2025).

In random Schrödinger cocycles, the modulus of continuity of the Lyapunov exponent as a function of parameters depends sharply on moment conditions; sufficient exponential moments ensure Hölder regularity, while weaker moments only imply weaker continuity (Duarte et al., 8 Sep 2025).

5. Applications: Physics and Probability

Moment Lyapunov exponents are integral in:

• Disordered Spin Chains and Ising Models: Systematic computation of cumulants and scaling exponents via cycle-expansion, with direct comparison to Monte Carlo simulations (Charbonneau et al., 2017).
• Stochastic Stability Analysis of PDEs: Assessing $q$7-th moment stability/instability of hybrid stochastic heat equations reveals phenomena such as stabilization of unstable deterministic systems via rapid regime switching (Bao et al., 2011).
• Statistical Mechanics Partition Functions: In monomer–dimer models, the MLE encapsulates the free energy per site, providing a probabilistic route to combinatorial enumeration in the thermodynamic limit (Vladimirov, 2012).
• Parabolic Anderson Model: Explicit formulas for all moment Lyapunov exponents demonstrate intermittency and growth exponents of solutions to stochastic equations with broad implications in random media (Borodin et al., 2012).
• Localization and Integrated Density of States: Regularity of Lyapunov exponents, and by duality, the density of states, is determined by sharp moment/probability constraints on the random environment (Duarte et al., 8 Sep 2025).

6. Exact Results and Limiting Cases

Exact solutions are available in several illustrative cases:

• For products of $q$8 matrices close to identity, the spectral problem simplifies to a second-order differential operator with explicit Gaussian solutions for the generalized Lyapunov exponent (Texier, 2019).
• In discrete-space parabolic Anderson models, all orders of moment exponents are computable via contour integrals and saddle points, with full characterization of intermittency (Borodin et al., 2012).
• For one-dimensional AR(1) Gaussian fields (elastic monomer–dimer), the largest eigenvalue of a transfer operator or solution to a pantograph equation directly yields the MLE (Vladimirov, 2012).
• In SDEs of Ornstein–Uhlenbeck or pitchfork type, explicit functional forms for $q$9 are obtainable, delineating regimes of finite/infinite values and singular growth (Baxendale, 21 Jul 2025).

7. Regularity, Obstructions, and Open Problems

Continuity and regularity of the Lyapunov exponent as a function of parameters like energy in random Schrödinger equations hinge crucially on suitable exponential moment conditions. Obstructions occur in the absence of such moments: while positivity and spectral gaps guarantee existence, Hölder continuity can fail catastrophically, confirmed by explicit constructions (Duarte et al., 8 Sep 2025). This connects to deep aspects of spectral theory and ergodic theory. The study of moment Lyapunov exponents thus interfaces with transfer operator theory, spectral analysis, large deviations, statistical mechanics, and stability of stochastic systems, with ongoing research targeting universality, finer fluctuation results, and extension to broader classes of random systems.