Chaos and sensitivity to initial conditions

Level SpecializationDifficulty ★★★★★Concept⌖ Open in the map

What is it?

Deterministic systems whose nearby trajectories separate exponentially, ∣δ(t)∣≈∣δ0∣eλt|\delta(t)| \approx |\delta_0|e^{\lambda t} with λ>0\lambda > 0. Prediction is possible only up to a horizon of about 1λln⁡tolerance∣δ0∣\frac1\lambda\ln\frac{\text{tolerance}}{|\delta_0|}: better measurements buy time only logarithmically.

Why does it exist?

Laplace imagined that knowing the present exactly would reveal the future. Poincaré (three bodies) and Lorenz (weather, 1963) showed that even simple deterministic equations amplify any uncertainty so fast that long-term prediction is impossible in practice.

Intuition

Lorenz restarted a simulation from a printout rounded to 3 decimals instead of 6; after a few simulated weeks the two forecasts had nothing in common. Stretching and folding — like kneading dough — makes neighbouring points separate exponentially while staying in a bounded region.

Formulas

λ=lim⁡t→∞1tln⁡∣δ(t)∣∣δ0∣\lambda = \lim_{t\to\infty}\frac1t\ln\frac{|\delta(t)|}{|\delta_0|}
largest Lyapunov exponent
Tpred≈1λln⁡Δ∣δ0∣T_{\text{pred}} \approx \frac{1}{\lambda}\ln\frac{\Delta}{|\delta_0|}
predictability horizon

Why does it matter?

Weather services run ensembles (dozens of forecasts from perturbed initial states) and report probabilities instead of a single answer. Chaotic simulations are not reproducible bit-for-bit across compilers or GPUs: a different rounding order is a tiny perturbation that grows. And chaos sets hard limits on what any model — learned or physical — can predict.

Where it shows up in computing

  • Weather and climate modelling★★★★★fundamentalPhysics and simulation

    Chaos limits deterministic forecasts to about two weeks; ensemble forecasting is the response.

  • Floating point (IEEE 754)★★★★★frequentScientific computing and algorithms

    Rounding differences (operation order, FMA, GPU) are amplified exponentially in chaotic simulations.

  • N-body gravitational simulation★★★★★frequentPhysics and simulation

    Three or more gravitating bodies are generically chaotic; long-term planetary predictions have finite horizons.

Where is it used?

Computing topics reachable from here, through the chain of ideas that leads to them:

What depends on it

Exercises

1Computing

A system has λ=0.5\lambda = 0.5 per day and initial error 10−610^{-6}. When does the error reach 1? And if the initial error is 10−1210^{-12}?

Solution

t=ln⁡(106)/0.5≈27.6t = \ln(10^6)/0.5 \approx 27.6 days. With 10−1210^{-12}: 55.355.3 days. A million times better data only doubles the horizon.

↑ ↓ to navigate · ↵ · Esc