[論文レビュー] Singularly perturbed phase response curves

Pierre Sacré, Alessio Franci|arXiv (Cornell University)|Jun 5, 2015

Nonlinear Dynamics and Pattern Formation参考文献 19被引用数 3

ひとこと要約

本稿では、特異摂動理論に基づく幾何学的手法を導入し、特異極限における緩和発振子の有限位相応答曲線（PRC）を予測する。これは無限小PRCアプローチの限界を克服するものであり、等位相線と発振子の位相空間のグローバル構造を分析することで、時間スケールの分離が大きくなるにつれて真のPRCに収束する半解析的公式を提供する。FitzHugh-Nagumoモデルに対してインパルスおよび方形パルスのシナリオで数値的に検証された。

ABSTRACT

In this paper we propose a novel geometric method, based on singular perturbations, to approximate isochrones of relaxation oscillators and predict the qualitative shape of their (finite) phase response curve. This approach complements the infinitesimal phase response curve approach to relaxation oscillators and overcomes its limitations near the singular limit. We illustrate the power of the methodology by deriving semi-analytic formula for the (finite) phase response curve of generic planar relaxation oscillators to impulses and square pulses of finite duration and verify its goodness numerically on the FitzHugh-Nagumo model.

研究の動機と目的

• 緩和発振子における有限振幅入力の特異極限において、無限小位相応答曲線（PRC）アプローチが失敗することを解消すること。
• 等位相線と特異摂動理論を用いて、有限PRCの定性的な形状を予測する幾何学的フレームワークを構築すること。
• 時間スケールの分離が大きくなるにつれて精度が向上するが、無限小アプローチとは異なり、半解析的手法を提供すること。
• 線形近似を超えて、高速・低速系における非線形的で有限入力効果を捉えるために、位相応答曲線解析を拡張すること。
• バースターや他の複雑な発振子における有限PRCの幾何学的理論の基盤を築くこと。

提案手法

• 本手法は、特異摂動理論を用いて、平面上の緩和発振子の高速・低速ダイナミクスを分析し、特異極限（ε → 0）における層状動的および低次元動的へ分解する。
• 特異周期軌道 γ⁰ は、臨界多様体 S⁰ 上の遅いセグメントと臨界繊維に沿った速いジャンプの和集合として構築される。
• 等位相線は、特異軌道の構造と S⁰ 上の遅い流れを用いて、漸近的位相写像を追跡することで幾何学的に近似される。
• 有限PRCは、摂動が特異軌道に沿って位相をどの程度ずらすかを、摂動の位置が遅いセグメントと速いセグメントのどちらに近いかに基づいて予測する。
• 摂動が同じブランチにとどまるか、反対側のブランチにジャンプを引き起こすかを区別する、区分的法則が導出される。
• 本手法はFitzHugh-Nagumoモデルに適用され、インパルスおよび有限持続時間の方形パルスに対するPRCの半解析的表現が得られる。

実験結果

リサーチクエスチョン

• RQ1 特異摂動理論の範囲を超えて、特異極限における緩和発振子の有限位相応答曲線（PRC）をどのように予測できるか。
• RQ2 等位相線や特異周期軌道といった位相空間の幾何的特徴が、有限PRCの形状にどのように寄与するか。
• RQ3 時間スケールの分離が大きくなるにつれて、なぜ無限小PRCが有限振幅入力効果を捉えられなくなるのか。
• RQ4 周期軌道に沿った摂動の位置が位相シフトに与える影響は何か。また、これらは区分的解析的法則で捉えられるか。
• RQ5 提案された幾何学的手法は、高速・低速系におけるインパルスおよび有限持続時間パルスの両方について、PRCを正確に予測できるか。

主な発見

• ε が小さい場合、数値シミュレーションと非常に良好に一致する。特にインパルスおよび方形パルスに対して顕著である。
• ε が小さい場合、予測されたPRCは、位相進みと遅れの大きさと位置を捉えており、ε → 0 に伴い誤差が減少する。
• 方形パルスに対しては、PRCに二つの明確な領域が存在することが正しく特定された：上側の折りたたみ付近での摂動により位相進みが生じる領域と、下側の折りたたみ付近での遅れが生じる領域。
• インパルスに対するPRCは、上側ブランチの開始を促進するため、上側折りたたみ領域付近でのみ非ゼロであり、位相進みの大きさは単調にゼロへ減少する。
• ε が大きい場合でも、予測は定性的に正確であり、PRCの主要な特徴を捉えているが、上側折りたたみ付近での小さな位相シフトはやや低く見積もられている。
• 本手法により、有限パルスが上側ブランチの終了を遅らせる現象が明らかになった。これはインパルスベースのPRCでは捉えられていないため、入力持続時間の重要性が強調される。