[論文レビュー] Convergence of Random Batch Method for interacting particles with disparate species and weights

Shi Jin, Lei Li+1 more|arXiv (Cornell University)|Mar 25, 2020

Markov Chains and Monte Carlo Methods参考文献 42被引用数 6

ひとこと要約

本稿は、異なる種類および重みを持つ相互作用粒子系に対するランダムバッチ法（RBM）の収束を確立し、$O(\sqrt{\tau})$ の強い誤差と $O(\tau)$ の弱い誤差を、$N$ に対して一様に示した。従来の同一粒子に対する結果を拡張するもので、非交換性を扱うために重み付き誤差平均化と洗練された補助補題を導入した。

ABSTRACT

We consider in this work the convergence of Random Batch Method proposed in our previous work [Jin et al., J. Comput. Phys., 400(1), 2020] for interacting particles to the case of disparate species and weights. We show that the strong error is of $O(\sqrtτ)$ while the weak error is of $O(τ)$ where $τ$ is the time step between two random divisions of batches. Both types of convergence are uniform in $N$, the number of particles. The proof of strong convergence follows closely the proof in [Jin et al., J. Comput. Phys., 400(1), 2020] for indistinguishable particles, but there are still some differences: since there is no exchangeability now, we have to use a certain weighted average of the errors; some refined auxiliary lemmas have to be proved compared with our previous work. To show that the weak convergence of empirical measure is uniform in $N$, certain sharp estimates for the derivatives of the backward equations have been used. The weak convergence analysis is also illustrating for the convergence of Random Batch Method for $N$-body Liouville equations.

研究の動機と目的

• 同一粒子から異なる種類および重みを持つ粒子系へ、ランダムバッチ法（RBM）の収束解析を拡張すること。
• 粒子数 $N$ に対して一様な強い収束と弱い収束の速度を確立すること。
• 重み付き／多様種粒子系における非交換性の欠如を克服するため、重み付き誤差平均化技術を導入すること。
• 弱収束が $N$ に対して一様となるように保証するため、後向き方程式の導関数に対する鋭い推定値を導出すること。
• 弱収束解析を通じて、RBM が $N$-体リウヴィル方程式へ応用可能であることを示すこと。

提案手法

• 各タイムステップ $\tau$ で粒子をランダムなバッチにグループ化し、相互作用をバッチ内に限定することで、RBM フレームワークを適応する。
• 異なる種類や重みによる非交換性に対応するため、誤差の重み付き平均を用いる。
• 弱収束を解析するため、後向きコルモゴロフ方程式のアプローチを適用し、解の導関数に対する鋭いバウンドを必要とする。
• 後向き方程式の導関数推定値に対してグローナルの不等式を適用し、$N$ に対して一様に誤差の増大を制御する。
• 高階導関数（$\nabla_{x_i}u$, $\nabla_{x_i}\nabla_{x_j}u$ など）の再帰的推定値を導出し、$k$ 階導関数に対して $O(1/N^k)$ の減衰を示す。
• テンソル表記と縮約規則を用いて、後向き方程式の生成子における複雑な導関数表現を管理する。

実験結果

リサーチクエスチョン

• RQ1 異なる重みと種類を持つ相互作用粒子系に対するランダムバッチ法の強い収束速度は何か？
• RQ2 異種粒子を有する $N$-粒子系において、経験的測度の弱収束はどのように振る舞うか？
• RQ3 非交換性を持つ粒子系において、強い誤差と弱い誤差の両方が $N$ に対して一様に収束する解析は可能か？
• RQ4 弱収束の一様性を保証するため、後向き方程式の導関数に必要な洗練された推定値は何か？
• RQ5 この手法は $N$-体リウヴィル方程式へどのように拡張可能か？

主な発見

• ランダムバッチ法の強い誤差は、非交換性を持つ異なる重み・種類の粒子に対しても、$O(\sqrt{\tau})$ のオーダーであり、$N$ に対して一様である。
• 弱い誤差は $O(\tau)$ のオーダーであり、$N$ に対して一様である。これは、平均場極限の一貫性にとって不可欠である。
• 強い収束の証明には、非交換性の欠如に起因し、同一粒子の場合とは異なり、誤差の重み付き平均化が必要である。
• 後向き方程式の導関数に対する鋭い推定値が得られ、$\|\nabla_{x_i}u\|_{\infty} = O(1/N)$ および $i \neq j$ のとき $\|\nabla_{x_i}\nabla_{x_j}u\|_{\infty} = O(1/N^2)$ であることが示された。これは弱収束に不可欠である。
• 解の高階導関数は、$k$ 階導関数に対して $O(1/N^k)$ のオーダーで減衰し、収束の一様性を支持する。
• 弱収束解析は、RBM が $N$-体リウヴィル方程式に適用される際の収束に基礎を提供し、理論的妥当性を拡張する。