Interpretable and Equation-Free Response Theory for Complex Systems

技術要約：複雑系における解釈可能かつ方程式を用いない応答理論

問題提起
統計力学における応答理論は、外部からの摂動に対してシステムの統計的性質がどのように変化するかを予測するための枠組みを提供する。しかし、この理論を複雑なマルチスケール系に適用する場合、大きな課題が生じる。従来の定式化は、多くの場合、明示的な発展方程式に依存しているが、高次元のシステムにおいては、それらは未知であるか、あるいは計算が困難である。さらに、標準的な応答公式は解釈性に欠けることが多く、変動の支配的なモードや、システムを支配する明確なタイムスケールを明確に分離できていない。揺動散逸定理（FDT）は平衡状態や特定の確率論的システムには適用されるが、散逸的なカオス系や、非線形解析を必要とするほど摂動が大きい場合には破綻する。これらのギャップを埋めるための既存の手法（機械学習サロゲートなど）は、多くの場合「ブラックボックス」として機能し、物理プロセスの階層的な影響を分離する能力に欠けている。

手法
本論文は、有限状態の離散時間マルコフ連鎖に対する線形および非線形応答公式を導出するためのフレームワークを提案する。このアプローチは、基礎となる微分方程式ではなく、システムのデータ駆動表現（マルコフ状態モデルなど）のみに依存する「方程式を用いない（equation-free）」設計となっている。

核となる手法は、以下の3つの要素を統合している：

1. マルコフ状態モデル (MSM): 複雑なシステムの連続的な位相空間を、有限個のメタステーブル（準安定）状態へと粗視化する。システムのダイナミクスは、観測またはシミュレーションされたデータから構築された確率遷移行列 $M$ によって近似される。
2. コープマン形式 (Koopman Formalism): 著者らは、観測量の無限次元空間におけるダイナミクスを線形化するために、コープマン作用素（遷移行列の転置 $K = M^T$ として定義される）を利用する。$K$ のスペクトル分解を行うことで、応答は特定の固有モード（コープマンモード）とその固有値に関連する寄与の和として表現される。
3. 摂動展開: 著者らは、時間依存の摂動に対する観測量および時間遅れ相関の応答に関する明示的な公式を導出する。摂動された確率測度 $\nu(n)$ を、小さなパラメータ $\epsilon$ のべき級数として展開する。これにより、$l$ 次の応答項に関する再帰的な関係式が得られ、これは未摂動の行列 $M$ のスペクトル特性を用いて解かれる。

本フレームワークは、あらゆる次数の摂動に対するグリーン関数の明示的な表現を提供する。線形応答については、グリーン関数はモード固有の係数で重み付けされた指数の和へと分解される。非線形応答については、公式はこれらのモードの畳み込みを含み、異なるタイムスケールがどのように相互作用するかを明示的に示す。

主な貢献

• 方程式を用いない応答公式: 本論文は、一意のエルゴード的不変測度を持つ任意のマルコフ連鎖に適用可能な、行列ベースの一般的な線形および非線形応答公式を導出している。これらの公式は、遷移行列、摂動行列、および観測量のみを必要とし、基礎となる発展方程式が未知であっても適用可能である。
• スペクトル分解による解釈性: 応答理論をコープマン論と結びつけることで、システムの応答が特定のコープマンモードによる明確な寄与に分解できることを示している。これにより、異なるタイムスケールの役割が明確になり、特定の強制力に対してどのモードが応答を支配するかを特定することが可能になる。
• 相関と非線形性への拡張: 本研究は、単純な観測量を超えて、時間遅れ相関および高次の（非線形）応答へと応答理論を拡張している。相関の感度に関する明示的な公式が提供されており、応答を「進化則の変化」と「不変測度の変化」に関連する項へと分解している。
• 時間依存の強制力への一般化: 静的な摂動に焦点を当てることが多かった従来の多くの研究とは異なり、本フレームワークは一般的な時間依存の変調 $f(n)$ を扱い、動的な強制シナリオに対する応答の予測を可能にする。

結果
本手法は、3つのメタステーブル状態を持つポテンシャルを備えた2次元ランジュバン方程式を用いた概念実証的な応用によって検証されている。

• モデル削減: システムはUlam法（625状態）を用いて離散化され、さらに支配的なコープマンモードに基づく3状態のマルコフモデルへと削減された。
• 線形応答: 観測量 $x$ および $y$ に対するグリーン関数を計算した。結果として、応答は各ケースにおいて単一のコープマンモード（強制力と観測量の特定の対称性に対応）によって支配されていることが示され、分解の解釈可能性が確認された。低次モデル（3状態）は、井戸内のダイナミクスを無視しているにもかかわらず、高解像度モデルの応答を高精度に再現した。
• 相関応答: 摂動に対する時間遅れ相関の感度が計算された。応答を4つの明確な項（観測量の期待値に関連するもの2つ、測度の変化に関連するもの1つ、進化則の変化に関連するもの1つ）に分解することで、ダイナミクスを捉えることに成功した。ここでも、低次モデルはフルシステムと強い一致を示した。
• 非線形応答: 2次のグリーン関数を計算した。対称性により、多くの交差項は消失し、固有値に対する単純な単項式の依存性が残った。結果は、非線形応答においても、低次モデルが本質的なダイナミクスを捉えていることを裏付けており、高次元のUlam離散化は、これとほぼ同一の結果を与えた。

意義と主張
本論文は、離刻時間信号を指数の和に適合させる統計的手法であるプロニー法（Prony method）に対し、「動的な基礎」を提供することを主張している。著者らは、グリーン関数のスペクトル分解がプロニー法の構造を説明すると論じている。すなわち、指数の数はマルコフ表現における離散状態の数に対応し、減衰率はコープマン作用素の固有値に対応している。

本研究の意義は、応答理論を以下のものにすることにある：

1. 解釈可能: 応答を変動の特定のモードに明示的に結びつけること。
2. 方程式に依存しない: 微分方程式ではなく、データ駆動型のマルコックスモデルに依拠すること。
3. スケーラブル: 低次モデルが複雑な応答を正確に予測できることを示すことで、高次元システムに伴う「次元の呪い」を回避できる可能性があること。

著者らは、本フレームワークが、微視的な発展方程式の完全な知識を必要とせずに、特に臨界転移付近における複雑系の感度を分析するための経路を提供するものであると結論づけている。また、現在の研究はマルコフ的表現を前提としているが、メモリ（記憶）や隠れた状態を持つシステムへの将来的な拡張が可能であることも述べている。