---
title: p-th Moment Lyapunov Exponent in Random Systems
url: https://www.emergentmind.com/topics/p-th-moment-lyapunov-exponent
type: topic
---

# p-th Moment Lyapunov Exponent in Random Systems

The \(p\)-th moment Lyapunov exponent is an asymptotic growth rate for moments of multiplicative random evolutions. Across stochastic PDEs, random media, random matrix products, and derived tangent flows, it appears in limits of the form
\[
\lim_{t\to\infty} t^{-a}\log \mathbb{E}[X_t^p]
\quad\text{or}\quad
\lim_{N\to\infty}\frac{1}{N}\log \mathbb{E}[e^{pS_N}],
\]
with the normalization \(a\) determined by the model. In the literature it plays several roles simultaneously: an annealed growth exponent, a scaled cumulant-generating function, a principal eigenvalue of a twisted operator, and, via Legendre transform, an input to large-deviation theory for finite-time Lyapunov exponents [1211.7125] [1910.09271] [1707.00708] [2507.16092].

## 1. Definitions and normalization conventions

The notion is not tied to a single normalization. In discrete and continuum parabolic Anderson models with white-in-time forcing, the natural scaling is linear in time. In long-range temporal-covariance models and fractional parabolic equations, the logarithm of the \(p\)-th moment grows superlinearly, so the correct normalization is \(t^{-\rho}\) or \(t^{-\kappa}\). In random matrix products, the corresponding object is usually the scaled cumulant-generating function of the finite-sample Lyapunov exponent. Some authors use the unnormalized quantity \(L(p)\), while others divide by \(p\) [1211.7125] [1509.08120] [1707.00708] [1604.03493].

| Setting | Definition | Time/length scale |
|---|---|---|
| Discrete-space PAM | \(\gamma_p(\beta;\nu)=\lim_{t\to\infty}\frac{1}{t}\log \mathbb{E}[u(t,\lfloor \nu t\rfloor)^p]\) | \(t\) |
| SHE with delta data | \(\lambda_p=\lim_{s\to\infty}\frac{1}{s}\log \mathbb{E}[\mathcal{Z}(s,0)^p]\) | \(s\) |
| Fractional-time Gaussian PAM | \(\lambda(p;\lambda)=\lim_{t\to\infty} t^{-\rho}\log \mathbb{E}[u_\lambda(t,x)^p]\) | \(t^\rho\) |
| Fractional PAM | \(\lambda_{p,\rho}=\lim_{t\to\infty} t^{-\kappa}\log \|u^\rho(t,x)\|_p\) | \(t^\kappa\) |
| Random matrices | \(L(p)=\lim_{N\to\infty}\frac{1}{N}\log \mathbb{E}\big[\|\prod_{i=1}^N T_i\|^p\big]\) | \(N\) |
| Hybrid stochastic heat equation | \(\lambda_p(u_0)=\limsup_{t\to\infty}\frac{1}{t}\log \mathbb{E}|u(t)|^p\) | \(t\) |

These conventions are not interchangeable without renormalization. The hybrid stochastic heat equation paper explicitly states that the \(p\)-th moment exponent is defined with no \(1/p\) normalization, while the dynamic random environment PAM defines \(\lambda_p(\kappa)\) with a prefactor \(1/(pt)\) [1111.1229] [1208.0330]. In random matrix theory, the paper on cycle expansions notes that some authors define \(\Lambda_p\equiv L(p)\), whereas others use \(\frac{1}{p}L(p)\) as a normalized exponent [1707.00708].

## 2. Exact formulas in integrable stochastic PDEs

In the \((1+1)\)-dimensional parabolic Anderson model with nearest-neighbor jumps and space-time white noise, Borodin and Corwin obtain exact contour-integral formulas for moments and then extract moment Lyapunov exponents by steepest descent and residue analysis. For the nearest-neighbor model with \(p+q=2\), they prove an explicit two-contour formula for \(\mathbb{E}[Z_\beta(t,n_1)Z_\beta(t,n_2)]\). In the symmetric case \(p=q=1\) and velocity \(\nu=0\), the first moment exponent is \(\gamma_1(\beta;0)=0\), while the second moment exponent is
\[
\gamma_2(\beta;0)=2\big(\sqrt{4+\beta^4}-2\big),
\]
with small-\(\beta\) asymptotics \(\gamma_2(\beta;0)\sim \beta^4/2\) and large-\(\beta\) asymptotics \(\gamma_2(\beta;0)\sim 2\beta^2-4\). In the totally asymmetric right-jump-only model, they derive all \(k\)-point moments by nested contour integrals and compute
\[
\gamma_k(1;\nu)=H_k(z_k^0),\qquad
H_k(z)=\frac{k(k-3)}{2}+kz-\nu\log\prod_{i=0}^{k-1}(z+i),
\]
where \(z_k^0\in(0,\infty)\) solves \(H_k'(z)=0\) [1211.7125].

The same paper compares the discrete one-sided model with the continuum stochastic heat equation and records the continuum value
\[
\gamma_k^{\mathrm{SHE}}=\frac{k^3-k}{24}
\]
for integer \(k\) at \(\beta=1\) and \(x=0\) [1211.7125]. A later paper establishes the continuum result for all real \(p>0\): for the one-dimensional SHE with delta initial data,
\[
\lambda_p=\lim_{s\to\infty}\frac{1}{s}\log \mathbb{E}[\mathcal{Z}(s,0)^p]=\frac{p^3-p}{24}.
\]
Equivalently,
\[
\lim_{t\to\infty}\frac{1}{t}\log \mathbb{E}\Big[\big(e^{t/12}\mathcal{Z}(2t,0)\big)^p\Big]=\frac{p^3}{12}.
\]
The same asymptotics imply the one-point upper-tail large deviation principle for the centered Hopf–Cole height \(V_t=\mathcal{H}(2t,0)+t/12\), with speed \(t\) and rate function \(\Phi_+(y)=\frac{4}{3}y^{3/2}\) [1910.09271].

A central point in these integrable models is that the exponent is not merely asymptotic but explicitly computable. In the discrete PAM this occurs through contour deformation, residues, and saddle-point analysis; in the SHE it occurs through refined asymptotics for fractional moments, so that the exact cubic law in \(p\) persists beyond integer replica indices [1211.7125] [1910.09271].

## 3. Fractional and colored-noise generalizations

For the parabolic Anderson model driven by Gaussian noise fractional in time and with spatial covariance \(\gamma(cx)=c^{-\alpha}\gamma(x)\), the correct time normalization is
\[
\rho=\frac{4-\alpha-2\alpha_0}{2-\alpha},
\]
equivalently \(\rho=\frac{4H-\alpha}{2-\alpha}\) in terms of the Hurst index \(H\). Under the stated scaling and technical assumptions, Xia Chen’s integer-moment formula extends to all real \(p\ge 2\):
\[
\lambda(p;\lambda)
=
\lim_{t\to\infty} t^{-\rho}\log \mathbb{E}[u_\lambda(t,x)^p]
=
p\left(\frac{p-1}{2}\right)^{\frac{2}{2-\alpha}}\mathcal{E}(\lambda),
\]
with
\[
\mathcal{E}(\lambda)=\lambda^{\frac{2}{2-\alpha}}\mathcal{E}(1).
\]
The lower bound in fact holds for all \(p>1\), but the exact upper bound in the note is proved for \(p\ge 2\) by hypercontractivity of the Ornstein–Uhlenbeck semigroup and the comparison
\[
\|u_{\frac{p-1}{q-1}\lambda}(t,x)\|_{L^q(\Omega)}
\le
\|u_\lambda(t,x)\|_{L^p(\Omega)},
\qquad q\ge p>1
\]
[1509.08120].

A different fractional generalization replaces the Laplacian by \(-(-\Delta)^{\alpha/2}\) and allows Gaussian noise colored in both time and space. In that setting the unified Stratonovich/Skorohod family \(u^\rho\) satisfies
\[
\lambda_{p,\rho}
=
\lim_{t\to\infty} t^{-\kappa}\log \|u^\rho(t,x)\|_p,
\qquad
\kappa=\frac{2\alpha-\beta-\alpha\beta_0}{\alpha-\beta},
\]
and the exact asymptotic constant is
\[
\lim_{t\to\infty} t^{-\kappa}\log \|u^\rho(t,x)\|_p
=
(p-\rho)^{\frac{\alpha}{\alpha-\beta}}\,M(\alpha,\beta_0,d,\gamma).
\]
Hence
\[
\lambda_{p,\mathrm{Strat}}=p^{\frac{\alpha}{\alpha-\beta}}M(\alpha,\beta_0,d,\gamma),
\qquad
\lambda_{p,\mathrm{Skor}}=(p-1)^{\frac{\alpha}{\alpha-\beta}}M(\alpha,\beta_0,d,\gamma).
\]
The difference reflects the renormalization of diagonal terms in the Skorohod interpretation. The basic integrability threshold for the \(\alpha\)-stable path functional is \(\alpha\beta_0+\beta<\alpha\) [1604.03493].

These two lines of work show that the label “\(p\)-th moment Lyapunov exponent” covers both linear-time and superlinear-time growth laws. The shared structure is the extraction of a deterministic rate from \(\log\)-moments, while the scaling and the \(p\)-dependence encode the underlying covariance singularities, spatial homogeneity, and interpretation of the noise [1509.08120] [1604.03493].

## 4. Analytical mechanisms for computing the exponent

In integrable stochastic PDEs, the starting point is often a Feynman–Kac or directed-polymer representation. For the discrete PAM, moments reduce to expectations over multiple random walks with pairwise local-time interactions, yielding a “quantum delta Bose gas” structure. For \(k=2\) in the nearest-neighbor case, and for all \(k\) in the one-sided case, the moment equations close on ordered domains and are solved explicitly by coordinate Bethe ansatz and nested contour integrals. The large-\(t\) exponent then emerges from contour deformation, simple-pole residues, and steepest descent; in the one-sided model, the dominant term comes from the “ground state” partition \(\lambda=(k)\) [1211.7125].

For the SHE with delta initial data, the Laplace transform of \(U=e^{t/12}\mathcal{Z}(2t,0)\) admits a Fredholm determinant representation. Fractional moments are accessed through the identity
\[
\mathbb{E}[U^{n-1+\alpha}]
=
\frac{(-1)^n}{\Gamma(1-\alpha)}
\int_0^\infty s^{-\alpha}\frac{d^n}{ds^n}\mathbb{E}[e^{-sU}]\,ds,
\]
and the dominant contribution is the \(L=1\) trace term in the Fredholm expansion. Higher wedge terms are shown to be exponentially smaller, which isolates the leading growth \(e^{p^3 t/12}\) and hence \(\lambda_p=(p^3-p)/24\) [1910.09271].

In non-integrable Gaussian-noise PAMs, explicit formulas are replaced by variational principles and semigroup methods. The extension from integer to real \(p\ge 2\) in the fractional-time model relies on Mehler’s formula and hypercontractivity for the Ornstein–Uhlenbeck semigroup [1509.08120]. In the fractional Laplacian model, the exponent is identified through a variational inequality and a Feynman–Kac type large deviation principle for additive functionals of symmetric \(\alpha\)-stable processes, producing the variational constant \(M(\alpha,\beta_0,d,\gamma)\) [1604.03493].

The same broad theme persists in dynamical-systems formulations. In cycle expansions, the exponent is encoded by the zero of a Ruelle dynamical zeta function. In transfer-operator approaches for random matrices, it is the logarithm of the largest eigenvalue of an averaged operator. In random dynamical systems on manifolds, it is the principal eigenvalue of a twisted generator acting on a weighted function space [1707.00708] [1907.08512] [2507.16092].

## 5. Random matrices, maps, and generalized moment exponents

For products of i.i.d. random matrices, the generalized Lyapunov exponent is the scaled cumulant-generating function
\[
L(p)=\lim_{N\to\infty}\frac{1}{N}\log \mathbb{E}\Big[\Big\|\prod_{i=1}^{N}T_i\Big\|^p\Big]
=
\lim_{N\to\infty}\frac{1}{N}\log \mathbb{E}[e^{pN\lambda_N}],
\]
and its derivatives at \(p=0\) generate the usual quenched Lyapunov exponent, the Lyapunov susceptibility, and higher cumulant densities. The cycle-expansion formalism expresses the relevant zeta function as a sum over pseudocycles,
\[
\zeta^{-1}(z,n)=\sum_{\Gamma\in P^\star}\mathcal{M}(\Gamma)\,p(\Gamma)\,z^{\ell(\Gamma)}\,e^{n\ln \rho(\Gamma)},
\]
so that \(\lambda=L'(0)\) and \(\chi_{\rm L}=L''(0)\) are computable as formally exact cycle sums. The method is particularly effective for susceptibilities and higher moments in disordered transfer-matrix models [1707.00708].

A complementary operator-theoretic construction for \(\mathrm{SL}(2,\mathbb{R})\) products defines
\[
\widetilde{\Lambda}(q)=\lim_{n\to\infty}\frac{1}{n}\ln \mathbb{E}\|\Pi_n x_0\|^q
\]
and reduces its computation to the largest eigenvalue of an averaged transfer operator. In the solvable continuum-limit \(AN\) subgroup case, the generalized moment exponent is obtained in closed form:
\[
\widetilde{\Lambda}(q)=-qw+\frac{q^2}{8}.
\]
The same framework yields exact single-integral formulas for the variance rates of \(\ln\|\Pi_n\|\) and \(\ln|\psi(x)|\) in one-dimensional Schrödinger transfer-matrix problems with Lévy noise [1907.08512].

For deterministic one-dimensional maps, the paper on generalized Lyapunov exponents shifts attention from the linear clock \(n\) to a dynamical instability sequence \(a_n=n^\alpha L(n)\). It defines
\[
\Lambda_\alpha
=
\left\langle
\frac{1}{n^\alpha L(n)}\sum_{k=0}^{n-1}\ln|T'(x_k)|
\right\rangle,
\]
and states that the deterministic framework suggests the natural extension
\[
\Lambda(p)\equiv \lim_{n\to\infty}\frac{1}{a_n}\ln\langle e^{pS_n}\rangle
=
\lim_{n\to\infty}\frac{1}{a_n}\ln\left\langle\left|\frac{\Delta x_n}{\Delta x_0}\right|^p\right\rangle.
\]
This accommodates super-exponential and sub-exponential instability classes. The cited examples include the infinite Bernoulli scheme with Lyapunov pair \((n\ln n,\infty)\), the ant-lion map with \((n^2,\gamma(A)/2)\), Thaler maps with \((n/\log n,\Lambda_1)\) or \((n^\alpha,\Lambda_\alpha)\), and the log-Weibull map with \((\ln n\ln\ln n,A_{LW})\) [1412.6867].

## 6. Spectral characterization on manifolds

For random dynamical systems generated by SDEs on a manifold \(M\), the derived process on the tangent bundle leads to a canonical finite-time Lyapunov observable. If \(v_t\) is the linearized flow and \(\theta_t=v_t/\|v_t\|\) is the induced diffusion on the unit sphere bundle \(SM\), then
\[
\log\|v_t\|-\log\|v_0\|
=
A_t
=
\int_0^t q_0(\theta_s)\,ds+\sum_{j=1}^m\int_0^t q_j(\theta_s)\circ dW_s^j.
\]
The \(p\)-th moment Lyapunov exponent is
\[
\Lambda(p)=\lim_{t\to\infty}\frac{1}{t}\log E^\theta[e^{pA_t}]
=
\lim_{t\to\infty}\frac{1}{t}\log \mathbb{E}[\|v_t\|^p].
\]
Its twisted generator is
\[
L_p=L+pY+pQ+\frac{p^2}{2}R,
\qquad
Y=\sum_{j=1}^m q_jX_j,\quad R=\sum_{j=1}^m q_j^2,
\]
or, on \(SM\),
\[
\mathcal{L}_p
=
L_{SM}
+
p\sum_{j=1}^m q_j(\theta)\widetilde{X}_j
+
pQ(\theta)
+
\frac{p^2}{2}\sum_{j=1}^m q_j(\theta)^2
\]
[2507.16092].

On compact manifolds, the principal eigenvalue of \(\mathcal{L}_p\) equals \(\Lambda(p)\). The non-compact case requires a weighted Banach space
\[
B_V=\{g:\|g\|_V=\sup_x |g(x)|/V(x)<\infty\},
\]
where \(V\) satisfies the growth condition (G). Under the paper’s hypotheses (H), (P), and (G), the twisted semigroup
\[
T_t^p g(x)=E^x[g(x_t)e^{pA_t}]
\]
is bounded and compact on \(B_V\), possesses a positive eigenfunction \(\phi_p\in B_V\), and satisfies
\[
\Lambda(p)=\log\kappa_p,
\qquad
L_p\phi_p=\Lambda(p)\phi_p.
\]
The map \(p\mapsto \Lambda(p)\) is convex and real analytic on intervals where the hypotheses hold [2507.16092].

This spectral identification feeds directly into finite-time fluctuation theory. When \(\Lambda\) is finite and essentially smooth, the Gärtner–Ellis theorem gives the large deviation rate function
\[
I(a)=\sup_{p\in\mathbb{R}}\{pa-\Lambda(p)\}.
\]
At \(p=0\), one has \(\Lambda'(0)=\mu\), the almost-sure top Lyapunov exponent, and
\[
\sqrt{t}\,\big(X_t-\mu\big)\Rightarrow \mathcal{N}(0,\Lambda''(0)),
\qquad
X_t=\frac{1}{t}\log\|v_t\|.
\]
If \(\Lambda''(0)>0\), moderate deviations occur with speed \(b_t^2/t\) and quadratic rate \(s\mapsto s^2/(2\Lambda''(0))\). The Ornstein–Uhlenbeck quadratic-functional example in the paper shows why the weighted-space formulation is essential on non-compact state spaces: without it, the eigenvalue problem may admit multiple positive solutions that do not characterize \(\Lambda(p)\) [2507.16092].

## 7. Intermittency, stability, and unresolved directions

A recurrent qualitative signature is intermittency. In the discrete PAM with delta initial data, the paper proves
\[
\widetilde{\gamma}_1(\beta;\nu)<\gamma_1(\beta;\nu)<\frac{\gamma_2(\beta;\nu)}{2}<\frac{\gamma_3(\beta;\nu)}{3}<\cdots,
\]
and equivalently that \(p\mapsto \gamma_p(\beta;\nu)/p\) is strictly increasing. The one-sided model and the symmetric nearest-neighbor model both exhibit this strict ordering, meaning that large moments are dominated by rare high peaks. The fractional PAM with colored noise yields the same qualitative picture through the superlinear factor \((p-\rho)^{\alpha/(\alpha-\beta)}\), so that \(p\mapsto \lambda_{p,\rho}/p\) is strictly increasing when \(\beta>0\) [1211.7125] [1604.03493].

In hybrid stochastic heat equations with Markovian switching, the \(p\)-th moment exponent is an exact variational quantity. For deterministic \(u_0\neq 0\),
\[
\lambda_p(u_0)
=
-p\lambda_{n_0}
+
\sup_{\mu\in M}
\left\{
\sum_{i=1}^N
\left[
pa_i+\frac{p(p-1)}{2}\beta_i^2
\right]\mu(i)-I(\mu)
\right\},
\]
and in the \(m\)-dimensional noise case one replaces \(\beta_i^2\) by \(\sum_{j=1}^m\beta_j(i)^2\). The paper also gives the almost-sure exponent
\[
\lambda_{\mathrm{as}}(u_0)
=
-
\left[
\lambda_{n_0}
-
\sum_{i=1}^N \pi_i\left(a_i-\frac{1}{2}\beta_i^2\right)
\right].
\]
It follows that \(\lambda_p(u_0)\) can exceed \(p\lambda_{\mathrm{as}}(u_0)\), and the paper provides an explicit two-regime example in which the solution is almost surely exponentially stable while the second moment grows exponentially [1111.1229].

The dynamic random environment PAM adds a different perspective. There the annealed exponent is defined, when it exists, by
\[
\lambda_p(\kappa)=\lim_{t\to\infty}\frac{1}{pt}\log \mathbb{E}[u(0,t)^p],\qquad p\in\mathbb{N},
\]
while the main proved results concern the quenched exponent \(\lambda_0(\kappa)\). The paper formulates the conjecture
\[
\lim_{\kappa\to\infty}\big[\lambda_p(\kappa)-\lambda_0(\kappa)\big]=0
\qquad\text{for all }p\in\mathbb{N},
\]
motivated by the heuristic that large diffusion lets the random walk spend most of its time in the largest favorable space-time clumps, thereby suppressing the difference between quenched and annealed exponential growth rates [1208.0330].

Taken together, these results show that the \(p\)-th moment Lyapunov exponent is not a single formula but a structural object. Its normalization reflects the temporal scaling of the model, its explicit form ranges from closed contour-integral saddle values to variational constants and principal eigenvalues, and its probabilistic content spans intermittency, moment stability, susceptibility, and large deviations [1211.7125] [1910.09271] [1707.00708] [2507.16092].

Source: https://www.emergentmind.com/topics/p-th-moment-lyapunov-exponent