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.

Pattern Formation in a Reaction-Diffusion-ode Model with Hysteresis in Spatially Heterogeneous Environments

Existence and Stability of Patterns in a Reaction-Diffusion-ode System with Hysteresis in Non-Uniform Media

Hysteresis-driven Pattern Formation in Reaction-Diffusion-ode Systems

Pattern Formation in a Diffusion-Ode Model with Hysteresis

Dynamics of a Reaction-diffusion-ODE System in a Heterogeneous Media

Stable patterns with jump discontinuity in systems with Turing instability and hysteresis

Pattern Selection in a Reaction-Diffusion Equation

Instability of Turing patterns in reaction-diffusion-ODE systems

Systems of reaction-diffusion equations with spatially distributed hysteresis

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

Exact Solutions for Pattern Formation in a Reaction-Diffusion System

Existence of Traveling Wave Solutions to Reaction-Diffusion-ode Systems with Hysteresis

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

Dynamics of a Reaction–diffusion–ode System with Quiescence

Stability And Hopf Bifurcation In A Reaction-Diffusion Model With Chemotaxis And Nonlocal Delay Effect

Pattern formation in reaction-diffusion systems with piece-wise kinetic modulation: an example study of heterogeneous kinetics

Pattern formations and instability waves for a Reaction–Diffusion system

Pattern Formation in a Cross-Diffusion System

Bifurcation Analysis of a Diffusion-Ode Model with Turing Instability and Hysteresis

Spatiotemporal Pattern Formation of A Diffusive Bimolecular Model with Autocatalysis and Saturation Law

Spatial Patterns of a Reaction-Diffusion Population System with Cross-Diffusion and Habitat Complexity