[論文レビュー] Particles approximations of Vlasov equations with singular forces : Propagation of chaos

Pierre‐Emmanuel Jabin, Maxime Hauray|arXiv (Cornell University)|Jul 19, 2011

Gas Dynamics and Kinetic Theory被引用数 17

ひとこと要約

本稿は、$d \geq 3$次元における $1/|x|^\alpha$ 形式の特異的相互作用力（$α < 1$）を有する粒子系について、平均場極限およびカオスの伝播を確立する。さらに、$N$-依存的なカットオフを導入することで、$α < d-1$ のより強い特異性を持つ力に対しても結果を拡張し、クーロン力や重力相互作用にほぼ完全に近い状況をカバーする。主な貢献は、このような特異的力則の下で $N$-粒子力学の極限として、ボルツマン方程式が厳密に導出されることの確立である。

ABSTRACT

We obtain the mean field limit and the propagation of chaos for a system of particles interacting with a singular interaction force of the type $1/|x|^α$, with $α&lt;1$ in dimension $d \geq 3$. We also provide results for forces with singularity up to $α&lt; d-1$ but with large enough cut-off. This last result thus almost includes the most interesting case of Coulombian or gravitational interaction, but it is also interesting when the strength of the singularity $α$ is larger but close to one, in which case it allows for very small cut-off.

研究の動機と目的

• $d \geq 3$ 次元における特異的力 $F(x) \sim 1/|x|^\alpha$ を有する $N$-粒子系の平均場極限を確立すること。
• このような系に対してカオスの伝播を証明し、粒子の周辺分布がボルツマン方程式の解に収束することを示すこと。
• 原点近傍に $N$-依存的なカットオフを導入することで、$α < d-1$ のより強い特異性を持つ力に対しても解析を拡張すること。
• ハミルトニアン構造やポテンシャルに基づく力に依存しないように一般の特異的カーネルを許容すること。
• プラズマおよび天体物理学における運動方程式の粒子シミュレーションの厳密な基礎を提供すること。

提案手法

• 粒子数 $N \to \infty$ の際に全質量が有限に保たれるように、力項に $1/N$ 要素を含む平均場スケーリングを用いる。
• 存在性と連続性を仮定した上で、積分方程式を用いて $N$-粒子系のグローバル解を定義する。
• 二項分布と指数モーメントを用いた確率的推定により、粒子密度の揺らぎを制御する。
• 位相空間をサイズ $L\varepsilon$ の立方体に分割し、これらの立方体上での粒子数の上界を評価する。
• チェビシェフの不等式と最適な指数モーメント評価を用いて、経験測度の $L^\infty$ ノルムに関する大偏差推定を得る。
• 特性曲線と電場 $E = F * \rho$ のリプシッツ連続性を用いて、ボルツマン方程式の解の存在性と正則性を確立する。

実験結果

リサーチクエスチョン

• RQ1 特異的力 $F(x) \sim 1/|x|^\alpha$ を有する $N$-粒子系が $N \to \infty$ の際にボルツマン方程式に収束する条件は何か？
• RQ2 ポテンシャルから導かれない、またはハミルトニアン構造を満たさない特異的力に対しても、カオスの伝播は確立可能か？
• RQ3 クーロン力や重力ポテンシャルに近い、より強い特異性を持つ力に対しても、解析をどのように拡張できるか？
• RQ4 $N$-依存的なカットオフが相互作用を正則化し、平均場極限を可能にする役割は何か？
• RQ5 このような特異的系において、カオスの伝播の収束速度はどの程度か？

主な発見

• 力が $|F(x)| \leq C/|x|^\alpha$ を満たし、$α < 1$ の場合、原点近傍でポテンシャルが連続的かつ有界であるため、平均場極限が成立する。
• 力が $α < d-1$ の場合、$N^{-m}$ のカットオフを導入することで平均場極限が確立され、$|x| \leq N^{-m}$ に対して $|F_N(x)| \leq N^{m\alpha}$ が成り立つ。これにより、ほぼクーロン的相互作用が可能になる。
• 経験測度に関する大偏差推定を用いてカオスの伝播が証明され、大偏差の確率は $N$ の任意の逆数より速く減少する。
• 指数モーメントと被覆論法を用いて、経験測度の $L^\infty$ ノルム $\|f_N^Z\|_\infty$ の上界を制御し、$\mathbb{P}(\|f_N^Z\|_\infty \geq \beta c_\phi \|f\|_\infty) \leq c_0 N^\gamma e^{-(\beta \ln \beta - \beta + 1)(2L)^n \|f\|_\infty N^{1-\gamma}}$ を得る。
• ボルツマン方程式の解はグローバルに存在し、特性曲線法により与えられる。$α < d-1$ かつ $\rho \in L^\infty$ であるため、$E$ はリプシッツ連続である。
• 分布関数 $f$ の $L^1$ および $L^\infty$ ノルムは、特性曲線に沿って保存され、解の安定性が保証される。