[論文レビュー] Effective Mori-Zwanzig equation for the reduced-order modeling of stochastic systems

Yuanran Zhu, Huan Lei|arXiv (Cornell University)|Feb 2, 2021

Model Reduction and Neural Networks参考文献 40被引用数 4

ひとこと要約

本稿は、ホイポエルリプティック半群推定を用いて、確率的系における有効モーリ＝ツワンツィグ（EMZ）方程式の厳密な解析的枠組みを構築し、観測可能統計量の事前推定を可能にするとともに、データ駆動型記憶核近似の収束を証明する。モンテカルロシミュレーションを用いたフェルミ＝パスタ＝ユリの鎖（FPU）モデルの検証により、記憶核が指数関数的に減衰することと、低次元モデルの精度が高く保たれることを示した。

ABSTRACT

Built upon the hypoelliptic analysis of the effective Mori-Zwanzig (EMZ) equation for observables of stochastic dynamical systems, we show that the obtained semigroup estimates for the EMZ equation can be used to drive prior estimates of the observable statistics for system in the equilibrium and non-equilibrium state. In addition, we introduce both first-principle and data-driven methods to approximate the EMZ memory kernel, and prove the convergence of the data-driven parametrization schemes using the regularity estimate of the memory kernel. The analysis results are validated numerically via the Monte-Carlo simulation of the Langevin dynamics for a Fermi-Pasta-Ulam chain model. With the same example, we also show the effectiveness of the proposed memory kernel approximation methods.

研究の動機と目的

• ホイポエルリプティック半群解析を用いて、確率的力学系における観測可能統計量（時間自己相関関数、記憶核、ゆらぎ力など）の厳密な事前推定を確立すること。
• 直交ダイナミクスの正則性推定を用いて、EMZ記憶核のデータ駆動型パrametrizationスキームの収束を証明すること。
• 1次元FPU鎖モデルにおけるランジュヴィンダイナミクスのモンテカルロシミュレーションを通じて、理論枠組みの数値的検証を実施すること。
• 半群推定を用いて、一般的に仮定されている指数関数的減衰記憶核を理論的裏付けから正当化すること。

提案手法

• モーリ型射影作用素とコルモゴロフ生成子を用いて、確率的微分方程式に対する有効モーリ＝ツワンツィグ（EMZ）方程式を導出する。
• 直交ダイナミクス半群 $ e^{-t\mathcal{Q}\mathcal{K}\mathcal{Q}} $ に対してホイポエルリプティック解析を適用し、元のSDEの性質から正則性と混合性を確立する。
• 半群推定を用いて、時間自己相関関数および記憶核の事前バウンドを導出し、平衡状態および非平衡状態の両方で指数的減衰を示す。
• 第一原理的およびデータ駆動型手法を導入してEMZ記憶核を近似し、記憶核の正則性推定を用いて収束を証明する。
• コルモゴロフ作用素 $ \mathcal{K}^n $ に対する観測可能量への作用を計算するための再帰的多項式展開技術を用い、展開係数の正確な計算を可能にする。
• 1次元FPU鎖におけるランジュヴィンダイナミクスのモンテカルロシミュレーションを用いて、予測された統計量とシミュレーション結果を比較することで、枠組みの妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ホイポエルリプティック解析からの半群推定を用いて、確率的系における観測可能統計量の事前バウンドを導出可能か？
• RQ2 平衡状態および非平衡状態におけるEMZ記憶核は指数関数的に減衰するか？その理論的正当化は可能か？
• RQ3 直交ダイナミクスの正則性推定を用いて、EMZ記憶核のデータ駆動型パrametrizationスキームの収束を証明可能か？
• RQ4 第一原理的およびデータ駆動型記憶核近似は、低次元モデルにおける全系の統計的挙動をどの程度正確に再現できるか？
• RQ5 提案された枠組みは、確率的FPU鎖における時間自己相関関数およびゆらぎ力の予測にどの程度有効か？

主な発見

• ホイポエルリプティック解析からの半群推定により、時間自己相関関数およびEMZ記憶核が、平衡状態および非平衡状態の両方で指数関数的に減衰することを示す事前バウンドが得られた。
• 記憶核の第一原理的パラメータ化手法は、EMZ方程式の理論的構造と整合的で、高い正確性を示した。
• データ駆動型記憶核近似スキームは収束が証明され、収束速度は半群推定から導かれた記憶核の正則性に依存する。
• FPU鎖の数値的シミュレーションにより、第一原理的およびデータ駆動型手法の両方が、全系統計と密接に一致する高精度な低次元モデルを生成した。
• 本枠組みにより、低次元モデルにおける広く用いられている指数関数的減衰記憶核の仮定が、厳密な解析的裏付けから正当化された。
• 明示的な多項式展開 $ \mathcal{K}^n r_j $ を用いることで、記憶核係数 $ k_n $ の体系的計算が可能となり、実装の容易さが向上した。