The existence and stability of spike equilibria in the one-dimensional Gray-Scott model: the pulse-splitting regime

Abstract The existence, stability, and pulse-splitting behavior of spike patterns in the one-dimensional Gray–Scott model on a finite domain is analyzed in the semi-strong spike-interaction regime. This regime is characterized by a localization of one of the components of the reaction near certain spike locations, while the other component exhibits a more global spatial variation across the domain. The method of matched asymptotic expansions is then used to construct k-spike equilibria in terms of a certain core problem. This core problem is studied numerically and asymptotically. For each integer k ≥ 1 , it is shown that there are two branches of k-spike equilibria that meet at a saddle-node bifurcation value. For small values of the diffusivity D of the second component, these saddle-node bifurcation points occur at approximately the same value. A combination of asymptotic and numerical methods is used to analyze the stability of these branches of k-spike equilibria with respect to both drift instabilities associated with the small eigenvalues and oscillatory instabilities of the spike profile. In this way, the key bifurcation and spectral conditions of Ei et al. [S. Ei, Y. Nishiura, K. Ueda, 2 n splitting or edge splitting? A manner of splitting in dissipative systems, Jpn. J. Ind. Appl. Math. 18 (2001) 181–205] believed to be essential for pulse-splitting behavior in a reaction–diffusion system are verified. By having verified these conditions, a simple analytical criterion for the occurrence of pulse-splitting is then formulated and confirmed with full numerical simulations of the Gray–Scott model. This criterion verifies a conjecture based on numerics and topological arguments reported in [A. Doelman, R.A. Gardner, T.J. Kaper, Stability analysis of singular patterns in the 1D Gray–Scott model: a matched asymptotics approach, Physica D 122 (1–4) (1998) 1–36]. The analytical results are compared with previously obtained results for pulse-splitting behavior.

[1]  Michael J. Ward,et al.  The Dynamics of Multispike Solutions to the One-Dimensional Gierer--Meinhardt Model , 2002, SIAM J. Appl. Math..

[2]  J. E. Pearson.  Complex Patterns in a Simple System , 1993, Science.

[3]  Cyrill B. Muratov,et al.  Static spike autosolitons in the Gray-Scott model , 2000 .

[4]  Guillermo H. Goldsztein,et al.  Dynamical Hysteresis without Static Hysteresis: Scaling Laws and Asymptotic Expansions , 1997, SIAM J. Appl. Math..

[5]  Yasumasa Nishiura,et al.  2n-splitting or edge-splitting? — A manner of splitting in dissipative systems — , 2001 .

[6]  U. Ascher,et al.  A collocation solver for mixed order systems of boundary value problems , 1979 .

[7]  Robert Gardner,et al.  A stability index analysis of 1-D patterns of the Gray-Scott model , 2002 .

[8]  M. Pino,et al.  THE GIERER & MEINHARDT SYSTEM: THE BREAKING OF HOMOCLINICS AND MULTI-BUMP GROUND STATES , 2001 .

[9]  Wentao Sun,et al.  The Slow Dynamics of Two-Spike Solutions for the Gray-Scott and Gierer-Meinhardt Systems: Competition and Oscillatory Instabilities , 2005, SIAM J. Appl. Dyn. Syst..

[10]  Global bifurcational approach to the onset of spatio-temporal chaos in reaction diffusion systems , 2001 .

[11]  Arjen Doelman,et al.  Pattern formation in the one-dimensional Gray - Scott model , 1997 .

[12]  Daishin Ueyama,et al.  Spatio-temporal chaos for the Gray—Scott model , 2001 .

[13]  Juncheng Wei,et al.  Hopf Bifurcations and Oscillatory Instabilities of Spike Solutions for the One-Dimensional Gierer-Meinhardt Model , 2003, J. Nonlinear Sci..

[14]  Robert D. Russell,et al.  Collocation Software for Boundary-Value ODEs , 1981, TOMS.

[15]  Cyrill B. Muratov,et al.  Traveling spike autosolitons in the Gray-Scott model , 2001 .

[16]  H. Swinney,et al.  Experimental observation of self-replicating spots in a reaction–diffusion system , 1994, Nature.

[17]  Michael J. Ward,et al.  Slow translational instabilities of spike patterns in the one-dimensional Gray-Scott model , 2006 .

[18]  Arjen Doelman,et al.  Slowly Modulated Two-Pulse Solutions in the Gray--Scott Model I: Asymptotic Construction and Stability , 2000, SIAM J. Appl. Math..

[19]  Kei-Ichi Ueda,et al.  Scattering and separators in dissipative systems. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  H. Meinhardt,et al.  A theory of biological pattern formation , 1972, Kybernetik.

[21]  Reynolds,et al.  Dynamics of self-replicating patterns in reaction diffusion systems. , 1994, Physical review letters.

[22]  Daishin Ueyama,et al.  Dynamics of self-replicating patterns in the one-dimensional Gray-Scott model , 1999 .

[23]  Daishin Ueyama,et al.  A skeleton structure of self-replicating dynamics , 1997 .

[24]  Shin-Ichiro Ei,et al.  The Motion of Weakly Interacting Pulses in Reaction-Diffusion Systems , 2002 .

[25]  Juncheng Wei,et al.  On ring-like solutions for the Gray–Scott model: existence, instability and self-replicating rings , 2005, European Journal of Applied Mathematics.

[26]  Cyrill B. Muratov,et al.  Stability of the Static Spike Autosolitons in the Gray--Scott Model , 2002, SIAM J. Appl. Math..

[27]  Michael J. Ward,et al.  The stability of spike solutions to the one-dimensional Gierer—Meinhardt model , 2001 .

[28]  Edmund J. Crampin,et al.  Reaction-Diffusion Models for Biological Pattern Formation , 2001 .

[29]  W. Burridge,et al.  “Excitability” , 1933 .

[30]  M. Georgiou.  Slow Passage Through Bifurcation and Limit Points. Asymptotic Theory and Applications in the Areas of Chemical and Laser Instabilities. , 1991 .

[31]  Robert Gardner,et al.  Stability analysis of singular patterns in the 1-D Gray-Scott model I: a matched asymptotics approach , 1998 .

[32]  Yasumasa Nishiura,et al.  Stability of singularly perturbed solutions to systems of reaction-diffusion equations , 1987 .

[33]  James Demmel,et al.  LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

[34]  Michael J. Ward,et al.  Numerical Challenges for Resolving Spike Dynamics for Two One‐Dimensional Reaction‐Diffusion Systems , 2003 .

[35]  Tsutomu Ikeda,et al.  Pattern Selection for Two Breathers , 1994, SIAM J. Appl. Math..

[36]  Michael J. Ward,et al.  Zigzag and Breakup Instabilities of Stripes and Rings in the Two-Dimensional Gray–Scott Model , 2006 .

[37]  Lee,et al.  Lamellar structures and self-replicating spots in a reaction-diffusion system. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[38]  Tasso J. Kaper,et al.  Axisymmetric ring solutions of the 2D Gray–Scott model and their destabilization into spots , 2004 .

[39]  Valery Petrov,et al.  Excitability, wave reflection, and wave splitting in a cubic autocatalysis reaction-diffusion system , 1994, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[40]  Stephen K. Scott,et al.  Autocatalytic reactions in the isothermal, continuous stirred tank reactor: Oscillations and instabilities in the system A + 2B → 3B; B → C , 1984 .

[41]  Masayasu Mimura,et al.  Layer oscillations in reaction-diffusion systems , 1989 .

[42]  John E. Pearson,et al.  Self-replicating spots in reaction-diffusion systems , 1997 .

[43]  Arjen Doelman,et al.  Slowly Modulated Two-Pulse Solutions in the Gray--Scott Model II: Geometric Theory, Bifurcations, and Splitting Dynamics , 2001, SIAM J. Appl. Math..