Phase Dynamics of Nearly Stationary Patterns in Activator-Inhibitor Systems

Aric Hagberg,Ehud Meron,Thierry Passot
DOI: https://doi.org/10.1103/PhysRevE.61.6471
2000-04-07
Abstract:The slow dynamics of nearly stationary patterns in a FitzHugh-Nagumo model are studied using a phase dynamics approach. A Cross-Newell phase equation describing slow and weak modulations of periodic stationary solutions is derived. The derivation applies to the bistable, excitable, and the Turing unstable regimes. In the bistable case stability thresholds are obtained for the Eckhaus and the zigzag instabilities and for the transition to traveling waves. Neutral stability curves demonstrate the destabilization of stationary planar patterns at low wavenumbers to zigzag and traveling modes. Numerical solutions of the model system support the theoretical findings.
Pattern Formation and Solitons
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the slow - dynamic behavior of periodic patterns close to the stable state in the activator - inhibitor system. Specifically, the author uses the phase - dynamics method to study the stability of these patterns and derives the Cross - Newell phase equation that describes the weak modulation of these periodic stable solutions. This equation is applicable to bistable, excitable media and Turing - unstable regions. ### Research Background and Problems 1. **Existing Research**: - Previous studies have mainly focused on the instabilities of local structures (such as pulses and spots) in excitable and bistable media. - Most of the research on periodic patterns has focused on the behavior near Turing bifurcations or the analysis near the codimension - two point of Turing and Hopf bifurcations. - Very few studies have explored the instabilities of periodic non - local stable patterns in excitable and bistable media, or periodic stable patterns far from Turing instabilities. 2. **Research Objectives**: - Derive a phase equation that describes the weak modulation of periodic stripe patterns in the FitzHugh - Nagumo (FHN) model. - Evaluate the stability thresholds of Eckhaus and zigzag instabilities and the transition from rest to traveling waves. - Verify the theoretical results through numerical simulations. ### Main Contributions - Propose a Cross - Newell phase equation applicable to periodic patterns with spatial - scale separation that occur far from Turing instabilities and in excitable and bistable media. - Through derivation and numerical simulation, determine the stability boundaries of periodic stripe patterns, including Eckhaus, zigzag instabilities and the transition to traveling waves. - Provide a theoretical basis for understanding complex spatial or spatio - temporal behaviors, such as the coupling between different wavenumber patterns. ### Mathematical Model and Formulas The FitzHugh - Nagumo (FHN) equations are as follows: \[ \frac{\partial u}{\partial t}=u - u^{3}-v+\nabla^{2}u, \] \[ \frac{\partial v}{\partial t}=\epsilon(u - a_{1}v - a_{0})+\delta\nabla^{2}v. \] Here, \(u\) is the activator, \(v\) is the inhibitor, and the parameters \(a_{0}\) and \(a_{1}\) can be chosen so that the FHN model represents an excitable medium, a bistable medium, or a system with Turing instability. ### Phase Equation Derivation By introducing a small parameter \(\lambda=\frac{1}{kL}\) and assuming that the modulation length scale \(L\) is much larger than the wavelength \(\frac{1}{k}\), the modulation solution can be expanded as an asymptotic series: \[ u(\theta, R, T)=u_{0}(\theta; k)+\lambda u_{1}(\theta, R, T)+\lambda^{2}u_{2}(\theta, R, T)+\cdots, \] \[ v(\theta, R, T)=v_{0}(\theta; k)+\lambda v_{1}(\theta, R, T)+\lambda^{2}v_{2}(\theta, R, T)+\cdots, \] where \(R = \lambda r\) and \(T=\lambda^{2}t\) are slow - space and time variables. The finally derived phase equation is: \[ \tau\frac{\partial\Theta}{\partial T}=-\nabla_{R}\cdot(kB), \] where: \[ \tau=\left<(\partial_{\theta}u_{0})^{2}\right>-\epsilon^{-1}\left<(\partial_{\theta}v_{0})^{2}\right>, \] \[ B = -\left<(\partial_{\theta}u_{0})^{2}\right>