On the role of positivity preservation for high order approximations of the Dean--Kawasaki equation

技術要約：ディーン–カワサキ方程式の高次近似における正値性保持について

問題提起
ディーン–カワサキ（Dean–Kawasaki）方程式は、$N$ 個の独立なブラウン粒子による経験的測度の進化を記述する。これは微視的なレベルでは原子的なランダム測量としてのマルチンゲール解として正確であるが、非リプシッツな平方根ノイズ係数とスケーリングの超臨界性により、密度値を持つ確率偏微分方程式（SPDE）として解釈するにはあまりに特異である。この粒子系を連続体モデルで近似するためには、正則化が必要となる。

この近似における決定的な課題は、**正値性の保持（positivity preservation）**である。経験的密度は非負である必要があるため、理想的な数値的または解析的な近似は、この性質を保持すべきである。しかし、既存の高次手法は以下のトレードオフに直面している：

1. 有限差分アプローチ（例：[CF23, CFIR26]）は、高密度領域において高次の弱収束を達成するが、発散演算子の非局所的な性質により、解が負の値を取ることを許容するため、正値性を保持できない。
2. 正値性保持アプローチ（例：[DKP24, DJP26]）は、非負性を保証するためにノイズを切断するか、あるいは特異性を正則化するが、収束率が弱くなるか、近似の質を制限する強いノイズ切断を必要とする。

本論文では、正値性保持の欠如が弱近似誤差に与える具体的な影響、特に初期密度が局所的に消失し得る低密度領域における影響を調査する。具体的には、正値性を保持しない近似が超多項式収束率を達成できるのか、また初期密度がゼロから離れていない場合に誤差がどのように振る舞うのかを問う。

手法
著者らは、以下のスペクトル正則化を用いたディーン–カワサキ方程式を分析する：
$$du_t = \frac{1}{2}\Delta u_t dt + \frac{1}{\sqrt{N}} \nabla \cdot \left( K_\varepsilon * (f(u_t^+) dW_t) \right)$$
ここで、$K_\varepsilon$ はフーリエ変換がコンパクトな台を持つ平滑化関数であり、$u_t^+ = \max(u_t, 0)$ である。$f$ は平方根ノイズ（例：$f(x)=\sqrt{x}$）を模したものとする。これまでの正値性保持スキームとは異与なり、ウィーナー過程 $W$ は切断されず、スペクトル正則化はノイズ係数全体に対して作用する。この選択により、解 $u_t$ が負の値を取ることを許容するが、正則化をテスト関数への畳み込みへと移動させ、テスト関数の正則性を活用できるという利点がある。

分析は以下に基づいている：

1. ハミルトン–ヤコビ–ベルマン（HJB）方程式との双対性： [KLvR19] の双対性の結果を複素数値のテスト関数へと拡張し、SPDE 解のモーメント母関数と無限次元コルモゴロフ後退方程式の解との関係を関連付ける。
2. ラプラス変換解析： ラプラス変換の差 $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$ を通じて弱誤差を測定する。これはコーシーの積分公式により、あらゆる次数のモーメントを制御する。
3. 誤差の分解： 全誤差を、初期近似誤差、フーリエカットオフ誤差、および解の負の部分（$u_t^-$）から特異的に生じる誤差の3つに分解する。
4. 下限の構成： 初期密度が部分区間で消失する特定の一次元例を構成することで、近似誤差の下限を導出し、非正値性保持スキームの限界を示す。

主要な貢献と結果

1. 時間一様弱誤差界：
   論文は、ラプラス変換を通じて測定される時間一様弱誤差界を確立している。誤差界は以下の3つの成分からなる：

   • 初期近似誤差。
   • フーリエカットオフ誤差（スペクトル分解能 $\varepsilon$ に依存）。
   • 負の部分誤差： 解が負の値を取る確率と大きさに比例する項。

2. 高密度領域 ($\rho_{\min} > 0$)：
   初期密度がゼロから離れている場合（$\rho_0 \ge \rho_{\min} > 0$）、解の負の部分は指数関数的に抑制される。その結果、スペクトル近似は滑らかなテスト関数 ($\phi \in C^\infty$) に対して超多項式収束率を達成する。具体的には、最適な $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$ の選択において、誤差は任意の $\alpha > 0$ に対して $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$ で減衰する。この結果は、密度がノイズ強度に対して十分に高い限り、厳密な正値性保持がなくても成立する。

3. 低密度領域 ($\rho_{\min} = 0$)：
   初期密度が消失することを許容する場合、負の部分の指数関数的抑制は消失する。上界は多項式レート $O(N^{-(3\alpha+d)/(2\alpha+d)})$ へと悪化する。
   極めて重要なことに、著者らは $\rho_{\min}=0$ の一次元例に対する下限を提示している。彼らは、この領域において弱誤差が超多項式的に減衰することはないと証明している。負の部分の解が相対的に大きくなり、ゆらぎを支配して近似の質を制限するためである。これは、正値性を保持しない近似は、低密度領域においては本質的に精度が制限されることを示唆している。

4. 既存手法との比較：

   • 次元 $d=1, 2$ において、非正値性保持スキームの低密度領域における上界は、[DKP24] における正値性保持SPDE近似よりも劣っている。
   • これらの結果は、ハイブリッドスキーム（高密度領域での高次スペクトル近似と、低密度領域での粒子ベースまたは正値性保持手法を組み合わせたもの）が最も効率的な戦略である可能性を示唆している。

5. 数値的検証：
   数値実験は理論的知見を裏付けている。シミュレーションによれば：

   • 高密度領域では、スペクトルガラーキン近似は滑らかな観測量に対して超多項式収束を示す。
   • 低密度領域では、解が負になる確率が著しく増加し、理論的な誤差の下限と相関している。
   • $\rho_{\min}$ が大きいときは正値性の喪失は稀であるが、$\rho_{\min} \approx 0$ のときは頻繁に発生する。

意義と主張
本論文は、低初期密度を許容するディーン–カワサキ方程式の非正値性保持近似に関する初の分析を提供したと主張している。その主要な意義は、正値性保持の欠如によるコストを定量化したことにある：

• 初期密度が消失する場合、正値性保持が超多項式収束率を達成するために不可欠であることを実証した。
• スペクトル正則化は高密度領域では正値性を厳密に保持せずに高次の精度を達成できるが、低密度領域では「負への逸脱（negative excursions）」が支配的な誤差源となり、超多項式収束を妨げることを明らかにした。
• 本研究は、無限次元のマルコフ構造とHJB方程式との双対性を利用することで、有限次元のアプローチ（[CF23] など）とは一線を画し、時間一様な誤差界と複素数値テスト関数の取り扱いを可能にしている。

著者らは、これらの結果の普遍性については慎重な姿勢を保っており、下限が特定の一次元例から導出されていること、また、粒子間の相互作用や外部ポテンシャル（これらは動的に低密度領域を作り出す可能性がある）との関係についてはさらなる研究が必要であることを述べているが、現在の分析はそのようなシナリオにおける近似の限界に関する基礎的な理解を提供するものである。