Belousov-Zhabotinsky Reaction: Oscillating Chemistry & Spiral Waves

Three competing chemical processes

The classic BZ reaction oxidises malonic acid with bromate in acidic solution, using cerium ions (Ce³⁺/Ce⁴⁺) as a redox catalyst — the colour change reveals the oscillation, vivid red and blue with a ferroin indicator. The oscillation arises because HBrO₂ autocatalyses its own production (positive feedback) while the accumulated oxidised catalyst simultaneously regenerates bromide, which quenches HBrO₂ (negative feedback with a delay) — the quintessential recipe for chemical and biochemical oscillations. The system is held far from equilibrium by continuous consumption of malonic acid and bromate, so it does not violate the second law of thermodynamics — it is a transient trajectory sustained only while reactants remain.

The Oregonator: reducing chemistry to three ODEs

Richard Field, Endre Körös and Richard Noyes at the University of Oregon reduced the full mechanism to a handful of elementary steps and then to three coupled ODEs known as the Oregonator (1974):

dx/dt = (1/ε)·[x − x² − fz·(x−q)/(x+q)]
dz/dt = x − z

x ≈ [HBrO₂] (fast activator),  z ≈ [Ce⁴⁺] (slow catalyst)
ε ≈ 0.04 (fast-slow timescale ratio)

Because ε is small, x evolves on a fast timescale while z is slow — this fast-slow structure produces the classic relaxation-oscillator shape: slow drift interrupted by rapid jumps, giving the sharp colour-change pulses seen experimentally, not smooth sinusoidal waves. Typical periods are tens of seconds, consistent with the 30–60 second oscillations observed at room temperature.

Spiral waves and Turing patterns

In an unstirred thin layer, diffusion couples neighbouring regions and the Oregonator becomes a reaction-diffusion PDE, producing two spatial patterns: target patterns — concentric rings expanding from pacemaker sites like ripples from a dropped pebble — and spiral waves, self-sustaining rotating spirals that are robust attractors of the system. This is directly analogous to cardiac tissue: ventricular fibrillation is a broken wavefront organising into a rotating spiral, the biological equivalent of a BZ spiral, which is why understanding BZ dynamics informs defibrillation strategy. Separately, Alan Turing showed in 1952 that a homogeneous steady state can spontaneously break spatial symmetry when an inhibitor diffuses faster than an activator — the diffusion-driven instability behind stripes and spots on animal skins, later confirmed experimentally in mouse palate ridge formation.

Frequently asked questions

Try it live

Everything above runs in your browser — open Belousov-Zhabotinsky Reaction and watch spiral waves and concentric rings self-organize from a simple excitable-medium cellular automaton.