Data-driven analysis of metastability in a stochastic bistable system

技術要約：確率的双安定系におけるメタ安定性のデータ駆動型解析

問題提起
複雑なシステムにおいて、複数の競合する状態が共存する現象であるメタ安定性（metastability）は、気候科学、生態学、工学などの分野で普遍的に見られる現象である。これらのシステムを理解する上での重要な課題は、メタ安定状態の安定性を定量化し、「ノイズ誘起型ティッピング（noise-induced tipping）」、すなわちノイズによって駆動される状態間の確率的な遷移を予測することである。平均脱出時間や脱出経路（インスタントン）を記述するために、大偏差原理（LDT）のような厳密な数学的枠組みが存在するが、これらは通常、基礎となる動力学方程式や位相空間の幾何学的構造（例：準ポテンシャル）に関する完全な知識を必要とする。観測データの部分的な時系列のみが利用可能な実世界の複雑なシステムにおいては、これらの幾何学依存的な手法は大きな限界に直面する。著者らは、システムの支配方程式や盆地境界（basin boundary）の幾何学に関する事前知識なしに、遅い状態間遷移と高速なサドル脱出プロセスを特定することに焦点を当てた、純粋にデータ駆動型のメタ安定性解析手法の必要性に取り組んでいる。

手法
著者らは、非線形な状態変数の進化から、無限次元空間における観測量の線形進化へと解析の対象を転換するコープマン作用素（Koopman operator）形式に基づくフレームワークを提案している。本研究では、平衡状態（$\mu=0$）および非平衡状態（$\mu \neq 0$）の両方を表現するように調整可能な、典型的な二次元確率的双安定系（回転する二重井戸ポテンシャル）を利用している。

コアとなる手法は以下の通りである：

1. データ生成： オイラー・丸山法（Euler-Maruyama scheme）を用いた確率微分方程式の数値シミュレーションを行い、様々なノイズ強度（$\alpha$）および回転パラメータ（$\mu$）の下でのアンサンブル時系列を生成する。
2. コープマン近似： **拡張動的モード分解（EDMD）**アルゴリズムを適用する。時系列データからコープマン作用素とその生成子を近似するために、多項式関数の辞書を使用する。
3. スペクトル解析： コープマン作用素の固有値および固有関数（モード）の特定を行う。本研究では以下に焦点を当てる：
   • 第一劣支配モード（$\phi_2$）： 最も遅い減衰率に関連し、メタ安定状態間の遷移を支配する。
   • より深い劣支配モード： より速い減衰率に関連し、井戸内のダイナミクスやサドルからアトラクタへの緩和を捉える。
4. 検証： データ駆動による推定値を、理論的な予測（具体的にはクラマースの法則およびその非可逆系への拡張、ならびにフリードリン・ウェントヴェル理論）と比較する。

主な貢献および結果

• 幾何学に依存しない盆地の特定： 第一劣支配コープマンモード（$\phi_2$）が、競合する吸引領域（basin of attraction）に対応する二つのほぼ不変な集合へと位相空間を分割することが示された。$\phi_2$ のゼロレベルセット（等位面）は、ノイズが存在する場合や系の数式が未知の場合であっても、決定論的な盆地境界を正確に再構成する。これは平衡および非平衡（回転）の両方のケースで成立する。
• 平均脱出時間の推定： 著者らは、平均脱出時間を推定するための二つのデータ駆動型アプローチを提示している：
  1. 時系列トラッキング： $\phi_2(x(t))$ のスカラー時系列を監視することにより、モードがゼロを横切る際に遷移を特定する。これらの横切り時間の統計は、指数関数的なスケーリングとサブ指数関数的な前因子（prefactor）の両方を含め、LDTの予測とよく一致する。
  2. スペクトル推定： コープマン生成子の第二固有値の絶対値の逆数（$1/|\lambda_2|$）は、アトラクタからの平均脱出時間に等しいことが示される。この関係は可逆系および非可逆系の両方で成立し、近年の理論的結果を裏付けている。
• 脱出プロセスの分解： コープマン・スペクトルのより深いモードを分析することで、以下の区別に成功している：
  • 井戸間の混合（Inter-well mixing）： $\phi_2$ によって支配される。
  • 井戸内の変動（Intra-well variability）： 安定不動点付近の複素固有値に関連する。
  • サドル脱出（Saddle escape）： 特定のモードが、サドル（エッジ状態）からアトラクタへの急速な緩和を捉えることを特定した。このモードの減衰率は、標準的なLDT近似ではしばしば無視される量である、サドル付近での滞在時間の平均値を近似する。
• 非平衡条件の影響： 回転成分（$\mu$）がメタ安定性にどのように影響するかを定量化している。$\mu$ を増加させるとシステムの不安定性が強まり、遷移の頻度が高まる。脱出時間のノイズ強度に対する理論的なスケーリングは依然として有効であるが、高い $\mu$ では横切り点の幾何学的な広がりが増大する。これは、インスタントン軌道がより広い領域にわたって盆地境界に接するため、サドルから離れた場所での脱出が発生するようになるためである。

意義および主張
本論文は、純粋にデータ駆動型の技術を用いて高次元のメタ安定系を研究するための、包括的かつマルチスケールなフレームワークを提供すると主張している。その主な意義は以下の点にある：

• 方程式および幾何学に依存しない解析： ドリフトベクトル場の明示的な形式や準ポテンシャルを必要とせずに、盆地境界を特定し、遷移統計を推定できる能力。
• 平衡および非平衡の統一的扱い： 本手法は可逆系および非可逆系の両方の確率系に適用可能であり、標準的なLDTの結果とデータ駆動型のスペクトル解析との間の溝を埋めるものである。
• エッジ状態への新たな知見： 脱出を媒介する「エッジ状態」としての緩和を記述する特定のコープマンモードの特定は、従来のラージ・デビエーション論ではほとんど考慮されてこなかった、遷移の媒介者としての新しい視点を提供する。

著者らは、本研究を高次元の複雑な系を解析するための第一歩と位置づけている。現在の研究は対称的な2次元モデルに焦点を当てているが、カーネル化されたEDMDや残差DMDなどの手法を通じて、原理（例えば、劣支配モードのゼロ等位線を用いて盆地境界を近似することなど）は非対称な系や高次元の系にも一般化できることが期待されており、「次元の呪い」を克服できる可能性があるとしている。本研究は、複雑なフラクタル状の盆地境界を解決すると主張しているのではなく、コープマンモードがそのような境界の「ノイズにより正則化された」推定値を提供することを示唆している。