Stochastic numerical approximation for nonlinear Fokker-Planck equations with singular kernels

技術要約：特異なカーネルを持つ非線形フォッカー・プランク方程式に対する確率的数値近似

問題設定
本論文は、以下の形式の非線形フォッカー・プランク偏微分方程式（PDE）の解の数値近似を扱う：
$$\partial_t u = \Delta u - \nabla \cdot (u K * u)$$
ここで、相互作用カーネル $K$ は特異である。このような方程式は、ケラー・セーゲル・ケモタキシス・モデル、静電気学や重力におけるクーロン・ポテンシャル、および2次元ナビエ・ストークス方程式のビオ・サバール則など、様々な物理的または生物学的文脈で発生する。

主な課題は、ドリフト項 $K * u$ が特異なカーネルを含むことにある。標準的な数値的手法、例えば、関連するマッケアン・ヴラソフ確率微分方程式（SDE）にオイラー・丸山（EM）スキームを直接適用する方法は、二つの困難に直面する：

1. 密度の計算： EMスキームは、各タイムステップにおいて解の密度を評価する必要があるが、これは直接的には利用できない。
2. 特異性とウェル・ポーズドネス（適切設定性）： タイムステップ $h \to 0$ となるにつれ、相互作用カーネルの特異性が粒子系のウェル・ポーズドネスを損なう可能性がある。

これらを克服するために、著者は、経験測量（empirical measure）を介して解を近似する $N$ 個の相互作用粒子系を検討する。しかし、特異なドリフトのため、数値的安定性と理論的な有界性を確保するために、モリファイア（平滑化関数）列 $(V^N)$ と滑らかなカットオフ関数 $F_A$ を用いて、系を正則化している。

手法
本論文では、正則化および離散化された粒子系（式 1.4/2.2）から、非線形フォッカー・プランク方程式の解への収束を分析する。その手法は以下の通りである：

1. 正則化と離散化：

   • 特異なカーネル $K$ は、モリファイア $V^N$（スケーリングパラメータ $\alpha$）によって畳み込み演算される。
   • ドリフトを一様に有界にするために、真のドリフト $K * u_t$ の $L^\infty$ 有界性と整合するように、カットオフ関数 $F_A$ が適用される。
   • 系は、タイムステップ $h$ を用いたオイラー・丸山スキームによって離散化される。

2. 収束解析戦略：
   解析では、モリファイ化された経験測量 $\tilde{\mu}^{N,h}_t$（および単一粒子の密度）と真の解 $u_t$ の間の誤差を、以下のいくつかの成分に分解する：

   • 初期誤差： 初期モリファイ化経験測量と $u_0$ の差。
   • カオス伝播誤差： 正則化された粒子系とマッケアン・ヴラソフSDE（非離散化）の差。
   • 離散化誤差： EMスキームによって導入される誤差（時間離散化）。
   • 確率的畳み込み： 粒子系のノイズから生じる項。

3. 技術的ツール：

   • セミグループ・アプローチ： 著者は、熱核 $g_2(t, \cdot)$ を含むフォッカー・プランク方程式の弱形式（mild formulation）を利用する。
   • 伊藤の公式： 粒子系の密度を熱核の周りで展開するために適用される。
   • モーメント有界性： 著者は、$L^p$ 空間におけるモリファイ化された経験測量の $N$ に関する一様な有界性を確立する。
   • グロンウォールの不等式： 誤差分解から生じる積分不等式を扱うために使用される。
   • 正則性仮定： カーネル $K$ は、特定の可積分条件（仮定 AK）および、$L^1 \cap L^r$ ノルムに関する畳み込み演算 $K * f$ のヘルダー型の正則性条件を満たすものと仮定される。

主要な貢献と結果
本論文は、粒子数極限（$N \to \infty$）およびタイムステップ極限（$h \to 0$）における二つの対象の明示的な収束レートを導出している：

1. モリファイ化された経験測量の収束（定理 2.5）：
   本論文は、モリファイ化された経験測量 $\tilde{\mu}^{N,h}_t$ が、$L^1 \cap L^r$ ノルムにおいてフォッカー・プランク方程式の解 $u_t$ に収束することを確立する。収束レートは以下のように与えられる：
   $$\mathbb{E}\left[ \sup_{s \in [0,t]} \| \tilde{\mu}^{N,h}_s - u_s \|^m_{L^1 \cap L^r} \right]^{1/m} \leq C \left( \mathbb{E}[\| \tilde{\mu}^{N,h}_0 - u_0 \|^m_{L^1 \cap L^r}]^{1/m} + N^{-\rho + \epsilon} + N^{v_2} h^{v_3} \right)$$
   ここで、指数はカーネルの特性およびモリファイ化パラメータ $\alpha$ に依存する。

2. 粒子密度の収束（定理 2.6）：
   本論文は、$i$ 番目の粒子の密度 $u^{i,N,h}_t(x)$ がマッケアン・ヴラフス密度 $u_t(x)$ へ向かう点ごとの収束レートを提供する。この境界にはガウス重み項が含まれる：
   $$|u^{i,N,h}_t(x) - u_t(x)| \leq C (\dots) \left( \int_{\mathbb{R}^d} g_{c/p}(t, x-z) u_0(dz) \right)^{1/p}$$
   この結果は、テスト関数がディラック測量であっても適用可能な、EMスキームの弱誤差を定量化する。

3. 特定のカーネルに対する明示的なレート（セクション 4）：
   著者は、一般的な結果を以下の特定の場合に適用している：

   • 有界リプシッツ・カーネル： $O(N^{-\frac{1}{d+2}} + h^{1/2})$ 次のレートを回収する。
   • リースのポテンシャルおよびクーロン・カーネル： $d \geq 2$ の場合、レートは $O(N^{-\frac{1}{2(d+1)}} + N^{\frac{d}{2(d+1)}} h^{1/2})$ である。
   • クリロフ・ロックナー型カーネル： $K \in L^p$ かつ $p > d$ を満たすカーネルに対して、同様のレートが導出される。

意義と主張
著者は、著者の知る限り、特異な相互作用カーネルを持つ相互作用粒子系に適用されるオイラー・丸山スキームに対して、定量的な収束レートを提供する最初の研究であると主張している。

• 特異性の処理： 以前の研究の多くは、ワッサースタイン距離におけるドリフトのリプシッツ連続性を必要としていたが（これは特異なカーネルでは成立しない）、本研究では、$L^1 \cap L^r$ 空間における畳み込み演算 $K * f$ の特定の正則構造を利用している。
• 離散化の解析： 非離散化（連続時間）粒子系に関する定量的な収束は最近 Olivera ら [35] によって確立されたが、本論文はこれを離散化されたEMスキームへと拡張し、粒子数 $N$ とタイムステップ $h$ の間のトレードオフを明示的に定量化している。
• 堅牢性： 本手法は、標準的なクリロフ・ロックナー条件を満たさないカーネル（例：$\alpha \geq 2$ である非常に特異なケラー・セーゲル・カーネル）に対しても、本論文で導入された特定の構造的仮定を満たす限り、堅牢であることが示されている。

著者は、導出されたレートが $N$ および $h$ に対して多項式的に依存すること、そして解析が次元 $d$ およびカーネルの特異性の強さへの依存関係を注意深く追跡していることを強調している。本研究は、理論的なカオス伝播の結果と、特異な相互作用モデルの実際的な数値実装との間の溝を埋めるものである。