[Paper Review] Stochastic integrability and the KPZ equation

Herbert Spohn|arXiv (Cornell University)|Apr 12, 2012

Stochastic processes and financial applicationsEconomics, Econometrics and Finance23 references16 citations

TL;DR

This paper establishes the stochastic integrability of the KPZ equation through connections to quantum integrable systems, demonstrating that the exact solution of the asymmetric simple exclusion process (ASEP) in the continuum limit yields the KPZ equation. The key result is the derivation of a Fredholm determinant formula for the height fluctuations, which matches the Tracy-Widom distribution from random matrix theory, confirming universal scaling behavior in non-equilibrium statistical mechanics.

ABSTRACT

The emerging field of stochastic integrability is outlined.

Motivation & Objective

• To establish a rigorous connection between stochastic processes and integrable systems in non-equilibrium statistical mechanics.
• To demonstrate that the KPZ equation, a fundamental model of interface growth, is stochastically integrable through exact solutions derived from ASEP.
• To show that the universal fluctuations of the KPZ equation match those of the GUE Tracy-Widom distribution, confirming universality beyond microscopic details.
• To explore the role of duality, Bethe ansatz, and replica methods in deriving exact solutions for stochastic PDEs.
• To bridge quantum integrability (via Lieb-Liniger model) and stochastic integrability (via KPZ equation) through generating function identities.

Proposed method

• Use of duality in the symmetric simple exclusion process (SSEP) to decouple correlation functions and enable exact solution.
• Application of the Bethe ansatz to the asymmetric simple exclusion process (ASEP) to compute eigenvalues and eigenvectors of the generator.
• Employment of Tracy-Widom's Fredholm determinant formula for the ASEP step initial condition to describe particle position distributions.
• Continuum limit of ASEP with scaling $\varepsilon \to 0$, $\varepsilon^{-2}$ time, and $\sqrt{\varepsilon}$ asymmetry to derive the KPZ equation.
• Derivation of the generating function for the height at the origin via a Fredholm determinant involving the Airy kernel and projection onto $[0,\infty)$.
• Establishment of a connection between $\log Z_n(t)$ in a stochastic heat equation and the open quantum Toda chain, using Macdonald functions as eigenfunctions.

Experimental results

Research questions

• RQ1 Can the KPZ equation be shown to be stochastically integrable despite the absence of probability amplitudes and a normalized partition function?
• RQ2 How do exact solutions of the ASEP in the continuum limit yield the KPZ equation and its universal fluctuations?
• RQ3 What is the precise mathematical structure linking the KPZ equation to quantum integrable systems such as the Lieb-Liniger model?
• RQ4 How do initial conditions (e.g., sharp wedge vs. flat) affect the universal scaling statistics of height fluctuations?
• RQ5 In what way do replica methods and determinantal formulas provide rigorous access to the moments and distribution of the KPZ height function?

Key findings

• The generating function for the KPZ height at the origin is given by a Fredholm determinant $\det(1 - P_0 K_{s,t} P_0)$, with $K_{s,t}$ defined via the Airy function and a parameter-dependent kernel.
• For large $t$, the height $h(t)$ scales as $-t/24 + (t/2)^{1/3} \xi$, where $\xi$ follows the GUE Tracy-Widom distribution.
• The solution of the KPZ equation with sharp wedge initial data is linked to the $n$-particle attractive Lieb-Liniger Hamiltonian through $\mathbb{E}(Z(0,t)^n) = \langle 0| e^{-tH_n} |0\rangle$, confirming a deep quantum-stochastic duality.
• The moments $\mathbb{E}(Z(0,t)^n)$ grow as $\exp(n^3)$, making rigorous control difficult, but replica methods yield consistent results.
• Experimental verification in turbulent liquid crystal films confirms the GUE Tracy-Widom distribution for sharp wedge initial data with high precision.
• For flat initial conditions ($h(x,0) = 0$), the height fluctuations switch to GOE Tracy-Widom statistics, showing that initial condition features persist in the universal scaling limit.