Abstract:Reaction–diffusion equations have been widely used to describe biological pattern formation. Nonuniform steady states of reaction–diffusion models correspond to stationary spatial patterns supported by these models. Frequently these steady states are not unique and correspond to various spatial patterns observed in biology. Traditionally, time-marching methods or steady state solvers based on Newton's method were used to compute such solutions. However, the solutions that these methods converge to highly depend on the initial conditions or guesses. In this paper, we present a systematic method to compute multiple nonuniform steady states for reaction–diffusion models and determine their dependence on model parameters. The method is based on homotopy continuation techniques and involves mesh refinement, which significantly reduces computational cost. The method generates one-parameter steady state bifurcation diagrams that may contain multiple unconnected components, as well as two-parameter solution maps that divide the parameter space into different regions according to the number of steady states. We applied the method to two classic reaction–diffusion models and compared our results with available theoretical analysis in the literature. The first is the Schnakenberg model which has been used to describe biological pattern formation due to diffusion-driven instability. The second is the Gray–Scott model which was proposed in the 1980s to describe autocatalytic glycolysis reactions. In each case, the method uncovers many, if not all, nonuniform steady states and their stabilities.

Rotating periodic patterns in reaction diffusion systems

The Existence of Positive Periodic Solutions for a Kind of Periodic Reaction-diffusion Equations

Target-Wave to Spiral-Wave Pattern Transition in A Discrete Belousov-Zhabotinsky Reaction Driven by Inactive Resin Beads

Pattern Selection in a Reaction-Diffusion Equation

Complex spatiotemporal phenomena in multi-dimensional reaction–diffusion system

Analysing transitions from a Turing instability to large periodic patterns in a reaction-diffusion system

A Selection Criterion for Patterns in Reaction-Diffusion Systems

Square Turing Patterns in Reaction-Diffusion Systems with Coupled Layers

Asymptotic Behavior Of Solutions For Systems Of Periodic Reaction-Diffusion Equations In Unbounded Domains

Solution landscape of reaction-diffusion systems reveals a nonlinear mechanism and spatial robustness of pattern formation

Spatial pattern formation in reaction–diffusion models: a computational approach

Spatiotemporal Dynamics of Reaction–Diffusion System and Its Application to Turing Pattern Formation in a Gray–Scott Model

Instability of Turing patterns in reaction-diffusion-ODE systems

Designing reaction-cross-diffusion systems with Turing and wave instabilities

Effects of square spatial periodic forcing on oscillatory hexagon patterns in coupled reaction-diffusion systems

Emergence of non-trivial solutions from trivial solutions in reaction–diffusion equations for pattern formation

Homoclinic snaking of contact defects in reaction-diffusion equations

Dynamics of periodic solutions in the reaction-diffusion glycolysis model: Mathematical mechanisms of Turing pattern formation

Two simple criterion to prove the existence of patterns in reaction-diffusion models of two components

Spatio-temporal dynamics of a reaction-diffusion system for a predator–prey model with hyperbolic mortality

Analysis of Spatiotemporal Patterns in a Single Species Reaction–diffusion Model with Spatiotemporal Delay