[論文レビュー] The strong convergence of operator-splitting methods for the Langevin dynamics model

Adam Telatovich, Xiantao Li|arXiv (Cornell University)|Jun 13, 2017

Theoretical and Computational Physics参考文献 34被引用数 6

ひとこと要約

本稿では、加法的ノイズを伴うランジュバンダイナミクスモデルにおける作用素分割法について、きめ細やかな強い収束解析を提示する。著者らは、クニタの解表現とイトウ・テイラー展開を活用し、強い収束次数が最大3に達する新たな高次作用素分割スキームを導出し、標準的な分割法（SVVやナード・スプリット）は唯一の収束次数にとどまることを示している。

ABSTRACT

We study the strong convergence of some operator-splitting methods for the Langevin dynamics model with additive noise. It will be shown that a direct splitting of deterministic and random terms, including the symmetric splitting methods, only offers strong convergence of order 1. To improve the order of strong convergence, a new class of operator-splitting methods based on Kunita's solution representation are proposed. We present stochastic algorithms with strong orders up to 3. Both mathematical analysis and numerical evidence are provided to verify the desired order of accuracy.

研究の動機と目的

• 加法的ノイズを伴うランジュバンダイナミクスモデルにおける広く用いられる作用素分割法の強い収束挙動を分析すること。
• 既存の分割スキーム（例えば、確率的速度バーレット（SVV）やナード・スプリット）の根本的限界を特定し、それらが強い収束次数1にしか達しないことを示すこと。
• 確率微分方程式のクニタの解表現に基づく、新たな高次作用素分割スキームのクラスを構築すること。
• 数学的解析と数値実験を通じて、理論的収束次数を検証すること。
• 提案された高次スキームの過減衰および過不足衰極限における振る舞いを検討すること。

提案手法

• 著者らは、イトウ型SDEのクニタの解表現を用い、ドリフトおよび拡散微分作用素の交換子を含む作用素指数関数として解を表現する。
• 彼らは、クニタ展開を段階的に高次に切り詰めることで、高次数の数値スキームを構築し、ドリフトおよび拡散項に対する微分作用素の反復適用によって係数を導出する。
• この手法では、最大3次の多重イトウ積分を計算する必要があり、例えば $ \bigint_{0}^{\tau} \bigint_{0}^{s_3} \bigint_{0}^{s_2} dW_{s_1} ds_2 ds_3 $ のような項が含まれるが、これらは解析的に評価され、精度を確保する。
• スキームは多重インデックス $ \boldsymbol{\nu} $ に従って項をグループ化され、合計レベル $ l(\boldsymbol{\nu}) + n(\boldsymbol{\nu}) $ によって次数が定まる。ここで $ l $ はウィENER積分の数を数え、$ n $ は微分作用素の適用回数を表す。
• 加法的ノイズの場合、強い収束次数3と次数2.5のスキームは同一となる。これは、すべての追加係数関数 $ g_{\boldsymbol{\nu}} $（$ \boldsymbol{\nu} \notin \tilde{\boldsymbol{\nu}}_{2.5} $）が、$ L^1 $ および $ L^0 $ 操作下で $ g_{(1,0,0)} $ および $ g_{(1,0,0,0)} $ が定数またはゼロであるため、消えるためである。
• 得られたスキームは、明示的かつ積分可能な部分ステップとして実装され、従来の分割法と同様に実装の簡便さを保ちつつ、より高い強い収束次数を達成している。

実験結果

リサーチクエスチョン

• RQ1 なぜ、分子動力学シミュレーションで広く用いられている標準的かつ作用素分割法（SVV やナード・スプリット）は、強収束次数1にしか達しないのか？
• RQ2 加法的ノイズを伴うランジュバンダイナミクスモデルの作用素分割法において、より高い強収束次数を達成できるか？
• RQ3 クニタの解表現とイトウ・テイラー展開を用いた高次スキームの構築の理論的根拠は何か？
• RQ4 強い収束誤差ノルムで評価した場合、新しいスキームは古典的手法と比べてどれほど高い精度を示すか？
• RQ5 新しい高次スキームは、過減衰および過不足衰極限においても正しい振る舞いを保持するか？

主な発見

• 確率的速度バーレット（SVV）およびナード・スプリット法は、部分ステップの正確さにもかかわらず、強い収束次数1にしか達しない。
• クニタの解表現を適用し、イトウ・テイラー展開を高次に切り詰めることで、強い収束次数が最大3に達する新たな作用素分割スキームを著者らが構築した。
• 加法的ノイズの場合、強い収束次数3と次数2.5のスキームは同一である。これは、すべての追加係数関数 $ g_{\boldsymbol{\nu}} $（$ \boldsymbol{\nu} \notin \tilde{\boldsymbol{\nu}}_{2.5} $）が、$ L^1 $ および $ L^0 $ 操作下で $ g_{(1,0,0)} $ および $ g_{(1,0,0,0)} $ が定数またはゼロであるため、消えるためである。
• 新しいスキームは数学的に望ましい強い収束次数に達することが証明されており、数値的証拠によって理論的予測が確認されている。
• これらの手法は過減衰および過不足衰極限の両方の状態に適用可能であり、両領域で一貫した収束挙動が観察された。
• 従来のイトウ・テイラー法とは異なり、力の高次微分の必要がなく、ルンゲ＝クッタ法とは異なり、複数回の力評価も不要であるため、分子動力学への実用的応用が可能である。