General conditions for Turing and wave instabilities in reaction -diffusion systems

The stability of a one-spike solution to a general class of reaction-diffusion (RD) system with both regular and anomalous diffusion is analyzed. The method of matched asymptotic expansions is used to construct a one-spike equilibrium solution and to derive a nonlocal eigenvalue problem (NLEP) that determines the stability of this solution on an O (1) time-scale. For a particular sub-class of the reaction kinetics, it is shown that the discrete spectrum of this NLEP is determined in terms of the roots of certain simple transcendental equations that involve two key parameters related to the choice of the nonlinear kinetics. From a rigorous analysis of these transcendental equations by using a winding number approach and explicit calculations, sufficient conditions are given to predict the occurrence of Hopf bifurcations of the one-spike solution. Our analysis determines explicitly the number of possible Hopf bifurcation points as well as providing analytical formulae for them. The analysis is implemented for the shadow limit of the RD system defined on a finite domain and for a one-spike solution of the RD system on the infinite line. The theory is illustrated for two specific RD systems. Finally, in parameter ranges for which the Hopf bifurcation is unique, it is shown that the effect of sub-diffusion is to delay the onset of the Hopf bifurcation.

The Turing (wave) instability is only possible in reaction-diffusion systems with more than one (two) components. Motivated by the fact that a time delay increases the dimension of a system, we investigate the presence of diffusion-driven instabilities in single-species reaction-diffusion systems with delay. The stability of arbitrary one-component systems with a single discrete delay, with distributed delay, or with a variable delay is systematically analyzed. We show that a wave instability can appear from an equilibrium of single-species reaction-diffusion systems with fluctuating or distributed delay, which is not possible in similar systems with constant discrete delay or without delay. More precisely, we show by basic analytic arguments and by numerical simulations that fast asymmetric delay fluctuations or asymmetrically distributed delays can lead to wave instabilities in these systems. Examples, for the resulting traveling waves are shown for a Fisher-KPP equation with distributed delay in the reaction term. In addition, we have studied diffusion-induced instabilities from homogeneous periodic orbits in the same systems with variable delay, where the homogeneous periodic orbits are attracting resonant periodic solutions of the system without diffusion, i.e., periodic orbits of the Hutchinson equation with time-varying delay. If diffusion is introduced, standing waves can emerge whose temporal period is equal to the period of the variable delay.

The role of a reaction-diffusion system in the initiation of skin organ primordia. I. The first wave of initiation

Reaction-Diffusion (RD) systems provide a computational framework that governs many pattern formation processes in nature. Current RD system design practices boil down to trial-and-error parameter search. We propose a differentiable optimization method for learning the RD system parameters to perform example-based texture synthesis on a 2D plane. We do this by representing the RD system as a variant of Neural Cellular Automata and using task-specific differentiable loss functions. RD systems generated by our method exhibit robust, non-trivial 'life-like' behavior.

The interplay between reaction kinetics and diffusion leads to a wide range of spatiotemporal behaviors in reaction–diffusion (RD) systems. This article presents a theoretical and computational study of a RD system inspired by experimental observations of the Rho-GEF-Myosin signaling network controlling cell contraction dynamics. The temporal reaction system dynamics range from periodicity to bistability. By employing a dual hybrid dynamical systems approach of numerical bifurcation analysis and Floquet theory, we characterize the spatiotemporal dynamics when stable limit-cycle oscillators or homogeneous time-periodic solutions undergo the process of Floquet–Turing-diffusion-driven-instability (FTDDI). FTDDI refers to the emergence of spatially nonuniform patterns when diffusion destabilizes an otherwise temporally stable limit cycle of the underlying reaction kinetics. For the temporal reaction system, numerical bifurcation allows us to identify regions in a two-parameter space defining the stability of the uniform steady states and regions where the system exhibits limit cycles, which are characterized by employing the Floquet theory. In the presence of spatial variations, Floquet theory classifies regions where diffusion destabilizes the limit cycle or maintains their homogeneous stability and the emerging spatiotemporal dynamics of the full system, far from equilibrium. In the bistable regime, diffusion differentiates regions into those that exhibit classical Turing diffusion-driven instability, leading to pattern formation; and those that exhibit FTDDI, leading to space-time periodic patterns, spatially inhomogeneous patterns, oscillatory pulses, and wave propagation, and those that remain unaffected by diffusion. These findings provide theoretical insights into understanding complex spatiotemporal dynamics relevant to biological, chemical, and ecological spatiotemporal systems.

We show that if a Turing instability occurs in a reaction–diffusion system with a nearly scalar diffusion matrix, then the parameters of the corresponding well-mixed system are necessarily such that the well-mixed system has at least two eigenvalues near zero. Conversely, if the corresponding well-mixed system is sufficiently close to a coalescence point of Hopf and saddle-node bifurcations (two eigenvalues are zero at such a point), and if the spatial domain is sufficiently large, then there exists a nearly scalar diffusion matrix such that a Turing instability occurs. These results imply that information on bifurcation loci from experiments in continuous-flow stirred tank reactor suffices to locate regions of parameter space where Turing instabilities are likely to occur; no knowledge of the reaction mechanism and rate constants is needed. In order to illustrate these results, we have analyzed a six-step model of the Belousov–Zhabotinskii reaction due to Showalter, Noyes, and Bar-Eli. In this model, the sixth step describes the reduction of the metal ion catalyst by the organic substrate. Depending on the rate constant of this step, we find that Turing instabilities can occur when the diffusion coefficient of either Ce4+ or of HOBr is very near, but slightly larger than that of Br−.

Time delays, modelling the process of intracellular gene expression, have been shown to have important impacts on the dynamics of pattern formation in reaction–diffusion systems. In particular, past work has shown that such time delays can shrink the Turing space, thereby inhibiting patterns from forming across large ranges of parameters. Such delays can also increase the time taken for pattern formation even when Turing instabilities occur. Here, we consider reaction–diffusion models incorporating fixed or distributed time delays, modelling the underlying stochastic nature of gene expression dynamics, and analyse these through a systematic linear instability analysis and numerical simulations for several sets of different reaction kinetics. We find that even complicated distribution kernels (skewed Gaussian probability density functions) have little impact on the reaction–diffusion dynamics compared to fixed delays with the same mean delay. We show that the location of the delay terms in the model can lead to changes in the size of the Turing space (increasing or decreasing) as the mean time delay, tau , is increased. We show that the time to pattern formation from a perturbation of the homogeneous steady state scales linearly with tau , and conjecture that this is a general impact of time delay on reaction–diffusion dynamics, independent of the form of the kinetics or location of the delayed terms. Finally, we show that while initial and boundary conditions can influence these dynamics, particularly the time-to-pattern, the effects of delay appear robust under variations of initial and boundary data. Overall, our results help clarify the role of gene expression time delays in reaction–diffusion patterning, and suggest clear directions for further work in studying more realistic models of pattern formation.

Global attractor of a coupled finite difference reaction diffusion system with delays

Reaction–diffusion processes across layered media arise in several scientific domains such as pattern-forming E. coli on agar substrates, epidermal–mesenchymal coupling in development, and symmetry-breaking in cell polarization. We develop a modeling framework for bilayer reaction–diffusion systems and relate it to a range of existing models. We derive conditions for diffusion-driven instability of a spatially homogeneous equilibrium analogous to the classical conditions for a Turing instability in the simplest nontrivial setting where one domain has a standard reaction–diffusion system, and the other permits only diffusion. Due to the transverse coupling between these two regions, standard techniques for computing eigenfunctions of the Laplacian cannot be applied, and so we propose an alternative method to compute the dispersion relation directly. We compare instability conditions with full numerical simulations to demonstrate impacts of the geometry and coupling parameters on patterning, and explore various experimentally relevant asymptotic regimes. In the regime where the first domain is suitably thin, we recover a simple modulation of the standard Turing conditions, and find that often the broad impact of the diffusion-only domain is to reduce the ability of the system to form patterns. We also demonstrate complex impacts of this coupling on pattern formation. For instance, we exhibit non-monotonicity of pattern-forming instabilities with respect to geometric and coupling parameters, and highlight an instability from a nontrivial interaction between kinetics in one domain and diffusion in the other. These results are valuable for informing design choices in applications such as synthetic engineering of Turing patterns, but also for understanding the role of stratified media in modulating pattern-forming processes in developmental biology and beyond.

Robust point control for a class of fractional-order reaction–diffusion systems via non-collocated point measurement

Turing bifurcation in a reaction–diffusion system with density-dependent dispersal

By discretising space into compartments and letting system dynamics be governed by the reaction–diffusion master equation, it is possible to derive and simulate a stochastic model of reaction and diffusion on an arbitrary domain. However, there are many implementation choices involved in this process, such as the choice of discretisation and method of derivation of the diffusive jump rates, and it is not clear a priori how these affect model predictions. To shed light on this issue, in this work we explore how a variety of discretisations and methods for derivation of the diffusive jump rates affect the outputs of stochastic simulations of reaction–diffusion models, in particular using Turing’s model of pattern formation as a key example. We consider both static and uniformly growing domains and demonstrate that, while only minor differences are observed for simple reaction–diffusion systems, there can be vast differences in model predictions for systems that include complicated reaction kinetics, such as Turing’s diffusion-driven instability model of pattern formation. Our work highlights that care must be taken in using the reaction–diffusion master equation framework to make predictions as to the dynamics of stochastic reaction–diffusion systems.

Conditional symmetries and exact solutions of nonlinear reaction–diffusion systems with non-constant diffusivities

We introduce a DNA-based reaction-diffusion (RD) system in which reaction and diffusion terms can be precisely and independently controlled. The effective diffusion coefficient of an individual reaction component, as we demonstrate on a traveling wave, can be reduced up to 2.7-fold using a self-assembled hydrodynamic drag. The intrinsic programmability of this RD system allows us to engineer, for the first time, orthogonal autocatalysts that counterpropagate with minimal interaction. Our results are in excellent quantitative agreement with predictions of the Fisher-Kolmogorov-Petrovskii-Piscunov model. These advances open the way for the rational engineering of pattern formation in pure chemical RD systems.

We investigate a reaction-diffusion system consisting of an activator and an inhibitor in a two-dimensional domain. There is a morphogen gradient in the domain. The production of the activator depends on the concentration of the morphogen. Mathematically, this leads to reaction-diffusion equations with explicitly space-dependent terms. It is well known that in the absence of an external morphogen, the system can produce either spots or stripes via the Turing bifurcation. We derive first-order expansions for the possible patterns in the presence of an external morphogen and show how both stripes and spots are affected. This work generalizes previous one-dimensional results to two dimensions. Specifically, we consider the quasi-one-dimensional case of a thin rectangular domain and the case of a square domain. We apply the results to a model of skeletal pattern formation in vertebrate limbs. In the framework of reaction-diffusion models, our results suggest a simple explanation for some recent experimental findings in the mouse limb which are much harder to explain in positional-information-type models.