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.

Identification of the Gray-Scott Model Via Deterministic Learning

Fractional Gray–Scott model: Well-posedness, discretization, and simulations

Analysis and Simulation of a Nonlocal Gray–Scott Model

A high-precision numerical approach to solving space fractional Gray-Scott model

Study on Turing Patterns of Gray-Scott Model Via Amplitude Equation

A Parameter Dispersion Based Global Identification Method and Its Application

Numerical simulation of the time fractional Gray-Scott model on 2D space domains using radial basis functions

Early Warning of Tipping in a Chemical Model with Cross-Diffusion Via Spatiotemporal Pattern Formation and Transition

Composite Molecules and Decoupling in Reaction Diffusion Models

Stability analysis of the stochastic Grey‐Scott model using spectral method

Pattern Formation Induced by Internal Microscopic Fluctuations

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

Modeling Unknown Stochastic Dynamical System Subject to External Excitation

Spatiotemporal Chaos of Self-Replicating Spots in Reaction-Diffusion Systems

Identification of the Fitzhugh-Nagumo Model Dynamics Via Deterministic Learning

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

Pattern Selection of Three Components Gray-Scott Model

Complex spatiotemporal phenomena in multi-dimensional reaction–diffusion system

On a reversible Gray-Scott type system from energetic variational approach and its irreversible limit

Identification of Damped and Driven Sine-Gordon Equation Based on Deterministic Learning

Approximation Identification for the Stochastic Time-Delayed Dynamical System