Raymond E. Goldstein,David J. Muraki,Dean M. Petrich
Abstract:In the bistable regime of the FitzHugh-Nagumo model of reaction-diffusion systems, spatially homogeneous patterns may be nonlinearly unstable to the formation of compact ``localized states." The formation of space-filling patterns from instabilities of such structures is studied in the context of a nonlocal contour dynamics model for the evolution of boundaries between high and low concentrations of the activator. An earlier heuristic derivation {[}D.M. Petrich and R.E. Goldstein, Phys. Rev. Lett. {\bf 72}, 1120 (1994){]} is made more systematic by an asymptotic analysis appropriate to the limits of fast inhibition, sharp activator interfaces and small asymmetry in the bistable minima. The resulting contour dynamics is temporally local, with the normal component of the velocity involving a local contribution linear in the interface curvature and a nonlocal component having the form of a screened Biot-Savart interaction. The amplitude of the nonlocal interaction is set by the activator-inhibitor coupling and controls the ``lateral inhibition" responsible for the destabilization of localized structures such as spots and stripes, and the repulsion of nearby interfaces in the later stages of those instabilities. The phenomenology of pattern formation exhibited by the contour dynamics is consistent with that seen by Lee, McCormick, Ouyang, and Swinney in experiments on the iodide-ferrocyanide-sulfite reaction in a gel reactor. Extensive numerical studies of the underlying partial differential equations are presented and compared in detail with the contour dynamics. The similarity of these phenomena (and their mathematical description) with those observed in amphiphilic monolayers, Type-I superconductors in the intermediate state, and magnetic fluids in Hele-Shaw geometry are emphasized.
What problem does this paper attempt to address?
This paper aims to study how space - filling patterns (such as labyrinth - like patterns) develop from finite - amplitude perturbations of a linearly stable homogeneous state in reaction - diffusion systems. Specifically, the authors focus on the situation in the bistable region of the FitzHugh - Nagumo model where the spatially homogeneous pattern may be nonlinearly unstable to form compact "local states". The paper studies the formation process of the space - filling patterns of these structures through a non - local contour dynamics model, which describes the evolution of the boundary between high - and low - concentration activators.
### Main problems and methods
1. **Research background**:
- Recent experimental studies have shown that in reaction - diffusion systems, the generation mechanism of space - filling patterns is significantly different from the classical Turing bifurcation. In the Turing bifurcation, periodic patterns are generated throughout the space through the linear instability of the homogeneous state. In this study, however, the space - filling "labyrinth - like" patterns develop from finite - amplitude perturbations of a linearly stable homogeneous state.
- It has been experimentally observed that these labyrinth - like patterns are composed of patches with different chemical compositions, which are separated by relatively sharp interfaces or fronts. The pattern formation process involves the movement of these fronts, and it is usually observed that these fronts repel each other, thus preventing self - crossing and topological changes.
2. **Model and hypothesis**:
- Based on the fast - inhibition limit in the FitzHugh - Nagumo model, the authors propose a non - local contour dynamics model. This model combines the Young - Laplace force related to the interface curvature and the screened Biot - Savart coupling between the distal segments of the interface.
- The non - local contribution reflects the "lateral inhibition" phenomenon, that is, the inhibition effect generated on a scale much larger than the thickness of the activator front, which leads to the destabilization of local structures such as points and stripes and the repulsion of nearby interfaces in the later stage.
3. **Theoretical analysis**:
- Through asymptotic analysis, the authors derive the detailed form of the contour dynamics. The results show that the normal component of the interface velocity contains a local contribution that is linear with the interface curvature and a non - local contribution in the form of a screened Biot - Savart interaction.
- The paper also explores three basic processes: stripe stabilization and repulsion, domain localization, and interface proliferation and transverse - front instability, which together explain the complex pattern - formation phenomenon.
4. **Numerical simulation**:
- The authors conduct extensive numerical simulations to verify the consistency between the contour - dynamics model and the original reaction - diffusion partial - differential equation. The simulation results show that by changing the control parameters, the entire process from the instability of the compact structure, the repulsion of the chemical front, to the relaxation of the branched structure back to the compact state can be observed.
### Conclusion
Through theoretical analysis and numerical simulation, the paper reveals the mechanism of complex pattern formation in reaction - diffusion systems, especially how the space - filling labyrinth - like patterns develop from finite - amplitude perturbations of a linearly stable homogeneous state. These findings not only help to understand pattern formation in chemical reactions but also have important connections with similar phenomena in other systems (such as amphiphilic - molecule monolayers, the intermediate state of type - I superconductors, and the Hele - Shaw geometry of magnetic fluids).