[論文レビュー] Dissecting the snake: the transition from localized patterns to isolated spikes in pattern formation systems

Nicolás Verschueren, Alan Champneys|Explore Bristol Research|Sep 20, 2018

Nonlinear Dynamics and Pattern Formation参考文献 43被引用数 4

ひとこと要約

本稿は、1次元反応拡散系における多重ピークを示す局所化されたパターンから孤立したスパイクへの遷移を調査し、空間固有値が複素数から実数に変わることで、無限個のマルチパルスホモクリニック軌道が同時に消失する、codimension-two点におけるBelyakov-Devaney遷移によってこの遷移が発生することを示している。解析により、遠方領域の減衰が振動的から単調的へと変わると、ホモクリニックスナキング構造が非連結な分岐へと崩壊することが明らかになった。これは一般化されたSchakenbergモデルおよび都市犯罪モデルにおける数値シミュレーションと整合的である。

ABSTRACT

An investigation is undertaken of coupled reaction-diffusion systems in one spatial dimension that are able to support, in different regions of their parameter space, either an isolated spike solution, or stable localized patterns with an arbitrary number of peaks. The distinction between the two cases is drawn through the behavior of the far field, where there is either an oscillatory or a monotonic decay. Several examples are studied, including the Lugiato-Lefever model, a generalized Schakenberg system that arises in cellular-level morphogensis and a continuum model of urban crime spread. In each, it is found that localized patterns connected via a so-called homoclinic snaking curve in parameter space transition into a single spike solution as a second parameter is varied, via a change in topology of the snake into a series of disconnected branches. In each case, the transition is caused by a so-called Belyakov-Devaney transition between complex and real spatial eigenvalues of the fair field of the primary pulse. A codimension-two problem is studied in detail where a non-transverse homoclinic orbit undergoes this transition. A Shilnikov-style analysis is undertaken which reveals the asymptotics of how the infinite family of folds of multi-pulse orbits are all destroyed at the same parameter value. The results are shown to be consistent with numerical experiments on two of the examples.

研究の動機と目的

• 1次元反応拡散系における多重ピークを示す局所化されたパターンがどのように孤立スパイク解に遷移するかのメカニズムを理解すること。
• 空間固有値の遷移（複素数から実数に）が、遠方領域の減衰挙動を振動的から単調的へと変化させることの役割を分析すること。
• codimension-two Belyakov-Devaney点における無限個のマルチパルスホモクリニック軌道の破壊を調査すること。
• ホモクリニックスナキング曲線の位相的変化とマルチパルス解の不安定化との関係を確立すること。
• 代表的な2つのモデル（一般化されたSchakenberg系と連続的都市犯罪モデル）における数値実験を通じて理論的予測の妥当性を検証すること。

提案手法

• 空間力学においてxを時間的変数として扱い、定常状態問題を4次元の可逆系として定式化する。
• 幾何学的特異摂動理論とShilnikovスタイルの解析を用いて、非横断的ホモクリニック軌道を含むcodimension-two分岐点近傍での挙動を研究する。
• 均一平衡状態の空間固有値が複素数から実数に変わる遷移を分析し、Belyakov-Devaney条件が臨界的メカニズムであることを同定する。
• マルチパルス軌道の折りたたみ位置の漸近的表現を導出し、ε → 0 のとき、それらが主パルスの折りたたみ点の近傍で指数関数的に集積することを示す。
• 数値継続法を用いて、ホモクリニックスナキング構造の崩壊が理論的予測と一致することを確認する。
• 細胞形態形成をモデル化した反応拡散系と都市犯罪拡散モデルの2つの物理的モデルを比較分析する。

実験結果

リサーチクエスチョン

• RQ1 遠方領域の減衰が振動的から単調的へと変わると、局所化パターン系におけるホモクリニックスナキング構造にどのような影響を与えるか？
• RQ2 Belyakov-Devaney遷移がマルチパルスホモクリニック軌道の破壊に果たす役割は何か？
• RQ3 なぜ無限個のマルチパルス解の折りたたみが、単一のパrameter値で同時に崩壊するのか？
• RQ4 系がcodimension-two分岐点を通過する際に、ホモクリニックスナキング曲線の位相的性質はどのように変化するか？
• RQ5 理論的漸近的予測と、現実のパターン形成系における数値計算との一致度はどの程度か？

主な発見

• 空間固有値が複素数から実数に変わることで発生するBelyakov-Devaney遷移が、ホモクリニックスナキング構造が非連結な分岐へと崩壊する原因となる。
• codimension-two点において、無限個のマルチパルスホモクリニック軌道のすべての折りたたみが同時に破壊され、ε → 0 のとき、それらのν値が主パルスのそれへと指数関数的に収束する。
• ν < 0 の場合のマルチパルス軌道の分離の漸近的レートはShilnikov解析によって予測され、一般化されたSchakenbergモデルおよび都市犯罪モデルの両方で数値的に確認された。
• この遷移は遠方領域の挙動の変化として特徴づけられる：振動的減衰（局所化されたパターン）から単調的減衰（孤立スパイク）に変わることで、グローバル分岐構造が変化する。
• 数値実験により、系がBelyakov-Devaney点を通過する際にホモクリニックスナキング曲線が非連結な孤立帯（isolas）に分裂することが確認され、理論的予測と整合的である。
• このメカニズムは頑健であり、対称性破壊項を含むSchnakenberg型モデルやLugiato-Lefever方程式など、多数の系で観察されており、広範な適用可能性を示唆している。