[論文レビュー] Analysis of biological integrate-and-fire oscillators

Marat Akhmet|arXiv (Cornell University)|Nov 25, 2010

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

ひとこと要約

本稿では、連続的な相互結合を伴う非同一なパルス結合インテグレート・アンド・ファイア（I&F）オシレーターに対する画期的な地図ベースの解析を提案し、特定の条件下では非同一なオシレーター間でも同期が発生することを示している。主な貢献は、ペスキンの第二予想を解決するための顕著な進展であり、同期領域の同定と、初期条件の集合が拡大するに従い同期時間が増加することの示唆を含む。

ABSTRACT

We consider discontinuous dynamics of integrate-and-fire models, which consist of pulse-coupled biological oscillators. A thoroughly constructed map is in the basis of the analysis. Synchronization of non-identical oscillators is investigated. Significant advances for the solution of second Peskin's conjecture have been made. Examples with numerical simulations are given to validate the theoretical results. Perspectives are discussed.

研究の動機と目的

• 非同一でパルス結合されたインテグレート・アンド・ファイア（I&F）オシレーターにおける同期に関する未解決のペスキンの第二予想に対処すること。
• 非発火間隔におけるオシレーター間の連続的結合を考慮する動的システムの枠組みを構築し、生物学的現実性を高めること。
• オシレーター間の状態遷移を追跡する特殊な地図（Poincaré地図ではない）を構築し、同期および周期的運動の解析を可能にすること。
• 非同一なオシレーターが初期条件が同一でない場合でも同期する条件を同定すること。
• 不連続な動的システムにおける周期的・準周期的・ほぼ周期的運動を研究するための理論的・数値的基盤を提供すること。

提案手法

• 1回の発火イベント後の1つのオシレーターの状態を、別のオシレーターの状態に写像するプロトタイプ地図 $ L $ が定義され、その後の写像では役割が入れ替わる。
• 地図 $ L $ は微分方程式 $ u' = f(u) $ の解を用いて構築され、$ s(v) $ は閾値に到達するまでの時間を表し、$ \bar{L}(v) = u(s, 0, 0) $ は別のオシレーターの状態を追跡する。
• 解析は動的システムの単調性および関数 $ g $ の性質に依存しており、特に $ g $ が $ \mathcal{K} $ 型であり、条件 (A2) を満たす場合に注目する。
• すべてのオシレーターが相互に結合されており、各発火は、発火していないすべてのオシレーターの状態にジャンプを引き起こし、ジャンプの大きさは $ \epsilon + \epsilon_i $ である。
• パrameter $ \mu_i, \xi_i, \epsilon_i $ を用いた摂動が導入され、非同一なダイナミクスを維持しながらも、特定の条件下で同期を保つことが可能になる。
• 理論的結果は、2オシレーターおよび多オシレーター系の数値シミュレーションにより検証され、同期時間および同期領域の大きさに注目する。

実験結果

リサーチクエスチョン

• RQ1 非同一でパルス結合されたインテグレート・アンド・ファイア（I&F）オシレーターが、初期条件が同一でない場合でも、どのような条件下で同期するか？
• RQ2 初期条件領域の大きさが拡大するに従い、同期時間はどのように変化するか、特にその拡大に伴う変化は？
• RQ3 非同一系において、発火イベントとオシレーター間の状態遷移の相互作用を捉える地図ベースのアプローチを構築可能か？
• RQ4 結合強度、閾値、ダイナミクスにおける摂動が、同期の維持または破壊に果たす役割は何か？
• RQ5 提案された地図フレームワークは、不連続な動的システムにおける周期的・ほぼ周期的・準周期的運動の解析に拡張可能か？

主な発見

• 関数 $ g $ が $ \mathcal{K} $ 型であり、条件 (A2) が成り立つ限り、$ x_2(t_0+) \in S_m $ のとき、同期は有限時間間隔 $[t_0 + \frac{m}{2}\tilde{T}, t_0 + mT]$ 内で達成される。
• 関数 $ g $ が $ \mathcal{K} $ 型ではあるが、(A2) が成立しない場合、システムは同期しない。これは同期のための必要条件であることを示唆する。
• 多オシレーター系では、発火していないすべてのオシレーターが $ S_0 $ にある場合、同期は $[t_0, t_0 + T]$ 内で達成されるが、一部のオシレーターが $ S_{k_i} $ 集合に属する場合は、より長い時間間隔で同期する。
• 同期に要する時間が、同期が必要な初期条件領域が拡大するに従い有界に拡大するため、非同期の初期値が存在することが示唆される。
• 非同期初期値の測度がゼロであるとは証明されていない（同一の場合とは異なり）が、本稿ではそのような領域を同定し、特徴づけている。
• 本手法は、時間変動する閾値 $ 1 + \phi(t,x) $、状態依存のジャンプ、不連続な右辺を伴うシステムへ一般化可能であり、地図 $ L $ が所定の性質を満たす限り有効である。