The nonlocal attraction-repulsion transport equation with power kernels

技術要約：べき乗カーネルを用いた非局所的引力・斥力輸送方程式

問題設定
本論文は、与えられた背景測度 $\omega$ への引力と、粒子間の自己斥力の競合によって駆動される、$\mathbb{R}^d$ 上の確率密度 $\phi_t$ の進化を記述する非局所的な連続方程式を調査している。ダイナミクスは以下の式によって支配される：
$$\partial_t \phi_t = \text{div}\left(\phi_t (\nabla V - K_r * \phi_t)\right), \quad t > 0,$$
ここで、背景ポテンシャルは $V = \psi_a * \omega$ であり、相互作用カーネルはべき関数 $\psi_s(x) = |x|^{1+s}$ および $K_s(x) = \nabla \psi_s(x) = (1+s)|x|^{s-1}x$ で定義される。指数 $a$（引力）および $r$（斥力）は $[0, 1)$ の範囲にある。

本研究は、$a \neq r$ である「不一致（discordant）」領域、および任意の質量 $\omega(\mathbb{R}^d)$ を持つ等指数 $a=r$ のケースに焦点を当てている。この方程式は、形式的には相互作用エネルギー
$$E(\phi) = \int_{\mathbb{R}^d} (\psi_a * \omega) d\phi - \frac{1}{2} \iint_{\mathbb{R}^d \times \mathbb{R}^d} \psi_r(x-y) d\phi(x) d\phi(y)$$
のワッサースタイン・勾配流である。主な数学的課題は、原点におけるカーネルの特異性（$s=0$ の場合）およびリプシッツ連続性の欠如にあり、これにより、標準的な $\lambda$-凸性を必要とする勾配流理論の直接的な適用が妨げられている。

手法
著者らは多角的な解析的手法を用いている：

1. 二乗半径正則化（Squared-Radius Regularization）： 特異性とリプシッツ連続性の欠如に対処するため、著者らは正則化 $\psi_{s,\epsilon}(x) = (|x|^2 + \epsilon)^{(1+s)/2}$ を導入している。畳み込みによる平滑化とは異なり、これは凸性と一方向の符号推定を正確に保持する。彼らは正則化された系における大域的なウェル・ポーズドネス（解の存在と一意性）を確立し、コンパクト性議論を用いて $\epsilon \to 0$ の極限へと移行する。
2. ラグランジュ的枠組み： 解は、特性曲線となる常微分方程式 $\dot{X}_t = -\nabla V(X_t) + K_r * \phi_t(X_t)$ を満たすフロー写像 $X_t$ による初期データの押し出し（push-forward）として構成および解析される。
3. 自由境界による特徴付け： 定常状態については、フーリエ解析と斉次分布の理論を用いる。著者らは、分数ラプラシアンを用いて畳み込み演算子の逆演算を行い、定常条件 $\phi \nabla F_\phi = 0$ を、分数型外部ディリクレ問題に関わる自由境界問題へと還元する。
4. エネルギーおよびモーメント推定： 正則化された系におけるグロンウォールの不等式型の議論を通じて、一様 $L^\infty$ 境界およびモーメント推定が導出される。これらは、解のタイトネス（tightness）および引力領域におけるサポートの閉じ込めを証明するために不可欠である。

主要な貢献と結果

• 大域的ウェル・ポーズドネス： 本論文は、有限のモーメントを持つ任意の $L^1 \cap L^\infty$ 初期データに対する大域的なラグランジュ的ウェル・ポーズドネスを確立している。解は、一様な $L^\infty$ およびモーメント境界を満たすラグランジュ的分布解のクラス内で一意的である。さらに、著者らは滑らかな初期データに対する $W^{n,\infty}$ 正則性の伝播を証明している。
• サポートの一様コンパクト性： 引力が支配的な領域（$a > r$、または $a=r$ かつ $\omega(\mathbb{R}^d) > 1$）において、初期データがコンパクトなサポートを持つ場合、その解のサポートはすべての時間において一様に有界であることを証明した。これは、遠距離において斥力に対して引力が優勢であるため、サポートの直径が無限に増大できないことを示すことで達成される。また、$a=r > 1$ の場合にはこの性質が成立しないことを示す反例も提示されている。
• 定常状態の特徴付け： 非二次的な引力支配領域（$0 \le r \le a < 1$）において、本論文はゼロフラックス定常状態を自由境界問題として特徴付けている。具体的には、定常密度 $\phi$ は $\phi = Q_{a,r} - L_r F$ と表される。ここで $L_r$ は分数ラプラシアンに関連する逆畳み込み演算子であり、$F$ は分数型外部ディリクレ問題を満たす。
• 明示的な定常状態の例： 得られた特徴付けを用いて、特定の幾何学的形状に対する明示的な定常測度を構成している：
  • 1次元： リーツ・ポテンシャルを含む区間解。
  • 2次元： 分数ラプラシアンのBoggio-Kelvinグリーン核を利用した円盤解。
  • 3次元： 一様球およびガウス背景における解。ここでは、特定の指数（$a=r=0$）により、特徴付けが局所微分演算子（ラプラシアン）に帰着する。
• 長時間収束： 境界付きエネルギーと一様モーメント境界を持つすべての大域的解は、ゼロフラックス定常状態に収束する部分列を持つことを証明している。$\omega$-極限集合が唯一の定常状態を含む場合、完全な収束が確立される。

意義と主張
著者らは、本研究を、任意の次元 $d \ge 1$ および不一致な領域 $a \neq r$ への、従来の一次元および等指数結果の拡張として位置づけている。

• 理論的進展： 本論文は、エネルギーが $\lambda$-凸ではない（この条件は標準的な勾配流の存在理論を通常排除する）特異な引力・斥力流を解析するための厳密な枠組みを提供している。二乗半径正則化の使用が、必要な構造的特性を保持するための重要な技術的革新として強調されている。
• モデリングの洞察： これらの結果は、「不平等なアクセス」における資源分配のプロトタイプモデルを提供する。引力支配の場合、平衡集団 $\phi_\infty$ は資源分布 $\omega$ に引き寄せられるが、内部的な斥力によって閉じ込められるため、定常状態は $\omega$ とは一致しない。$\phi_\infty$ の明示的な特徴付けにより、富や資源の分布における不平等の定量化が可能となる。
• MMDとの関連： 本論文は、これらの輸送方程式と機械学習で使用される最大平均偏差（MMD）との関係を明確にしている。先行研究では滑らかなカーネルや特定のケースが扱われてきたが、本論文は高次元における特異な負距離カーネルへの解析を拡張し、$a=r$ かつ $\omega(\mathbb{R}^d) \ge 1$ の場合のウェル・ポーズドネスと漸近挙動を確立している。
• 数値的検証： 自由境界による特徴付けから導かれた理論的な定常プロファイルは、数値的な粒子シミュレーションによって検証されており、一次元および二次元の両方で一致が示されている。

著者らは、サポートの一様閉じ込めに関する結果が、$\omega$ がコンパクトに支持されていない場合の $a=r$ かつ $\omega(\mathbb{R}^d)=1$ のケースには適用されないことを明記している（その場合、サポートは $\omega$ に一致するように増大することが予想される）。また、$a=r$ かつ $\omega(\mathbb{R}^d)=1$ のケースは、Rosenzweig, Slepčev, および Wang による並行研究において、平均場収束率に焦点を当てて扱われていることを認めている。