Hydrodynamic Limit Theorem

Convergence of the empirical density of an interacting particle system (Markov dynamics) to a deterministic PDE under macroscopic scaling, exemplified by exclusion processes.
Hydrodynamic Limit Theorem

Prerequisites

Empirical density field

Let (ηt)t0(\eta_t)_{t\ge 0} be a conservative interacting particle system on a discrete torus TNd=(Z/NZ)d\mathbb{T}_N^d=(\mathbb{Z}/N\mathbb{Z})^d with Markov generator LNL_N (described by a / ).

A standard macroscopic observable is the empirical measure (density field)

πtN(du)=1NdxTNdηt(x)δx/N(du), \pi_t^N(du) ={} \frac{1}{N^d}\sum_{x\in\mathbb{T}_N^d}\eta_t(x)\,\delta_{x/N}(du),

viewed as a random measure on the continuum torus Td\mathbb{T}^d. Equivalently, for a smooth test function GG,

πtN,G=1NdxTNdηt(x)G(x/N). \langle \pi_t^N, G\rangle ={} \frac{1}{N^d}\sum_{x\in\mathbb{T}_N^d}\eta_t(x)\,G(x/N).

Theorem (hydrodynamic limit, typical form)

Fix a macroscopic time scale and a macroscopic initial profile ρ0:Td[0,1]\rho_0:\mathbb{T}^d\to[0,1]. Assume the initial distributions of η0\eta_0 are “associated” to ρ0\rho_0 in the sense that π0N\pi_0^N converges (in probability) to ρ0(u)du\rho_0(u)\,du.

Then, under an appropriate scaling of time (diffusive or hyperbolic, depending on the model), the path (πtN)t[0,T](\pi_t^N)_{t\in[0,T]} converges in probability to a deterministic trajectory ρ(t,u)du\rho(t,u)\,du, where ρ\rho solves a macroscopic PDE of conservation/diffusion type.

Example: symmetric simple exclusion process (SSEP)

For SSEP, the natural scaling is diffusive: observe the process at times tN2tN^2. The hydrodynamic limit states that ρ(t,u)\rho(t,u) solves the heat equation on Td\mathbb{T}^d:

tρ(t,u)=Δρ(t,u),ρ(0,u)=ρ0(u). \partial_t \rho(t,u) = \Delta \rho(t,u), \qquad \rho(0,u)=\rho_0(u).

Example: asymmetric exclusion (ASEP, 1D)

For ASEP with nonzero drift, the natural scaling is often hyperbolic: observe at times tNtN. In one dimension, the limiting density typically solves a scalar conservation law (inviscid Burgers type):

tρ(t,u)+u(ρ(t,u)(1ρ(t,u)))=0,ρ(0,u)=ρ0(u), \partial_t \rho(t,u) + \partial_u\big(\rho(t,u)(1-\rho(t,u))\big)=0, \qquad \rho(0,u)=\rho_0(u),

under appropriate entropy-solution interpretation.

What “hydrodynamic limit” encodes

  • Law of large numbers at the macroscopic scale: fluctuations vanish after scaling, yielding deterministic PDE behavior.
  • Local equilibrium principle: microscopic configurations near a macroscopic point behave approximately like equilibrium Gibbs measures at the local density.
  • Connection to dynamical large deviations: hydrodynamic limits are often paired with path-space LDPs; the time-averaged DV theory ( ) is one cornerstone in this direction.

Typical proof inputs (conceptual)

Many rigorous proofs rely on:

  • martingale formulations of πtN,G\langle \pi_t^N, G\rangle derived from the generator,
  • replacement lemmas (local functions replaced by functions of the empirical density),
  • and entropy/relative-entropy methods comparing the law of ηt\eta_t to local equilibrium references.