An Asymptotic Analysis of Spike Self-Replication and Spike Nucleation of Reaction-Diffusion Patterns on Growing 1-D Domains

TL;DR

This paper analyzes the conditions under which spike self-replication and nucleation occur in reaction-diffusion systems on growing 1-D domains, providing asymptotic criteria, phase diagrams, and numerical validation for different models.

Contribution

It establishes precise asymptotic conditions for spike bifurcations in reaction-diffusion systems on growing domains, validated by numerical simulations.

Findings

Spike self-replication occurs via saddle-node bifurcation.

Spike nucleation occurs through boundary-value problem bifurcation.

Analytical thresholds accurately predict bifurcation points in simulations.

Abstract

In the asymptotic limit of a large diffusivity ratio, certain two-component reaction-diffusion (RD) systems can admit localized spike solutions on a 1-D finite domain in a far-from-equilibrium nonlinear regime. It is known that two distinct bifurcation mechanisms can occur which generate spike patterns of increased spatial complexity as the domain half-length L slowly increases; so-called spike nucleation and spike self-replication. Self-replication is found to occur via the passage beyond a saddle-node bifurcation point that can be predicted through linearization around the inner spike profile. In contrast, spike nucleation occurs through slow passage beyond the saddle-node of a nonlinear boundary-value problem defined in the outer region away from the core of a spike. Here, by treating L as a static parameter, precise conditions are established within the semi-strong interaction asymptotic regime to determine which occurs, conditions that are confirmed by numerical simulation and continuation. For the Schnakenberg and Brusselator RD models, phase diagrams in parameter space are derived that predict whether spike self-replication or spike nucleation will occur first as L is increased, or whether no such instability will occur. For the Gierer-Meinhardt model with a non-trivial activator background, spike nucleation is shown to be the only possible spike-generating mechanism. From time-dependent PDE numerical results on an exponentially slowly growing domain, it is shown that the analytical thresholds derived from the asymptotic theory accurately predict critical values of L where either spike self-replication or spike-nucleation will occur. The global bifurcation mechanism for transitions to patterns of increased spatial complexity is further elucidated by superimposing time-dependent PDE simulation results on the numerically computed solution branches of spike equilibria.