Pattern selection of three components Gray-Scott model

Analytical stability conditions and weakly nonlinear analysis (2019)

Key insight

Pattern formation in reaction–diffusion systems is governed by a balance between:

  • Reaction kinetics → drives local growth or decay
  • Diffusion → redistributes species and can destabilize uniform states

In this system:

  • A homogeneous steady state can be stable without diffusion
  • The same state becomes unstable with diffusion, leading to spatial patterns (Turing instability)

Near the instability threshold, the type of pattern (stripes vs spots) is not arbitrary —
it is determined by nonlinear interactions captured by amplitude equations.

System

Three-component Gray–Scott reaction–diffusion model:

  • One activator + two inhibitors
  • Coupled nonlinear reaction terms with diffusion
  • Spatially extended system exhibiting diffusion-driven instability

This framework captures how microscopic reaction dynamics give rise to macroscopic spatial structure.

What I did

  • Derived analytical conditions for Turing instability using linear stability analysis
  • Computed the dispersion relation to identify unstable wavelength bands
  • Applied Routh–Hurwitz criteria to determine stability of the homogeneous steady state
  • Performed weakly nonlinear analysis to derive amplitude equations near the bifurcation

Methods & tools

  • Model: three-component Gray–Scott reaction–diffusion system
  • Linear analysis:
    • steady-state solutions and Jacobian eigenvalues
    • Routh–Hurwitz stability conditions
    • dispersion relation for diffusion-driven instability
  • Nonlinear analysis:
    • multiple-scale expansion
    • Fredholm solvability condition
    • derivation of amplitude equations
  • Outputs: instability criteria, dispersion curves, amplitude equations, pattern classification

Key results

1. Diffusion induces instability

  • A stable homogeneous state becomes unstable when diffusion is introduced
  • This leads to stationary spatial patterns (Turing instability)

2. Diffusion controls pattern scale

  • The dispersion relation defines a band of unstable wave numbers
  • Changing diffusion coefficients shifts the system across the instability threshold

3. Nonlinear effects determine pattern type

  • Linear theory predicts whether patterns form
  • Nonlinear analysis determines which pattern emerges

4. Pattern selection mechanism

Amplitude equations classify stable regimes:

  • Stripe (labyrinth) patterns
  • Hexagonal (spot) patterns
  • Mixed states depending on parameter regime

Takeaway

Pattern formation is not random —
it emerges from a precise interplay between diffusion, reaction kinetics, and nonlinear interactions.

This work shows how analytical methods can predict not only when patterns form, but which structures are selected.

Publication

H. Rahim, N. Iqbal, C. Cong, and Z. Ding
Pattern selection of three components Gray–Scott model
Journal of Physics: Conference Series (2019)
DOI · PDF

← Back to Projects