Reaction-Diffusion

Reaction-diffusion systems are mathematical models describing how the concentrations of substances change under the combined influence of local chemical reactions and spatial diffusion. They are the canonical framework for understanding how spatial pattern and structure can emerge spontaneously in homogeneous media — a process central to developmental biology, morphogenesis, chemical ecology, and condensed matter physics.

The governing equations, introduced independently by Alan Turing (1952) and by the Soviet chemist Boris Belousov (1951, published posthumously), take the form of coupled partial differential equations in which reaction terms describe local production and consumption, while diffusion terms describe spatial spreading. The counterintuitive result — that diffusion, typically a smoothing and homogenizing process, can destabilize a uniform state and produce patterns — is known as the Turing instability.

The Belousov-Zhabotinsky reaction provided the first experimental confirmation: a chemical system that spontaneously produces traveling waves and spiral patterns without any external spatial cue. Since then, reaction-diffusion dynamics have been identified in skin pigmentation, seashell patterning, cardiac arrhythmias, and the aggregation of slime molds.

From a systems-theoretic perspective, reaction-diffusion systems demonstrate that pattern is a dynamical property, not a structural template. The same equations with different parameters can produce spots, stripes, labyrinths, or uniform fields. The pattern is not encoded in the initial conditions; it is selected by the dynamics as the system relaxes toward a stable attractor. This makes reaction-diffusion systems a mathematical paradigm for emergence: the whole acquires properties — spatial structure — that are absent from the parts and from the governing equations themselves, which specify only local interactions.