Numerical stationary states for nonlocal Fokker-Planck equations via fixed points of consistency maps

技術要約：一貫性写像の不動点を用いた非局所フォッカー・プランク方程式の数値定常状態

問題設定
本論文は、以下の形式を持つ非局所フォッカー・プランク型方程式の定常状態の数値計算について扱っている：
$$\partial_t u = \nabla \cdot (u\nabla(H'[u] + \Phi[u](x) + V(x))),$$
ここで、$u(x,t)$ は密度、$H(u)$ は内部エネルギー、$V(x)$ は外場ポテンシャル、そして $\Phi[u]$ は非局所的（通常は非線形的）な相互作用を表現する。これらの式は、流体力学、相分離、細胞接着、生物学的凝集など、多様な分野で発生する。

主な課題は、これらの系における長時間ダイナミクスにある。標準的な時間依存型ソルバー（例：有限体積法）は、動的に安定な状態を検出することには長けているが、不安定な定常状態や複雑な分岐枝を特定することはしばしば困難である。これは、メタ安定性や複数のタイムスケールを持つ系において顕著である。さらに、既存の数値継続ソフトウェアは、このクラスの方程式に見られる一般的な非線形拡散や非局所項の特定の構造的要件への対応に苦慮することが多い。

手法
著者らは、微分作用子の直接的な離散化を回避する、不動点に基づいた数値フレームワークを提案している。代わりに、定常問題を、元の偏微分方程式とその非局所的相互作用項から導かれる非線形不動点方程式 $T u - u = 0$ として再定式化する。

1. 不動点への再定式化:
   $\partial_t u = 0$ とし、質量の保存則と正値性を利用することで、定常状態 $u^*$ はその支持集合上で一致条件 $H'[u] + \Phi[u] + V = c_0$ を満たすことが示される。$H'$ が可逆であると仮定すると、定常解は写像の不動点として以下のように表される：
   $$T u := (H')^{-1}[c_0(u) - \Phi[u](x) - V(x)].$$
   線形拡散（$H[u] = \sigma u(\log u - 1)$）の場合、この逆関数は指数関数となる。

2. 数値解法戦略:
   著者らは、離散化された系 $F_n(u) = T_n(u) - u = 0$ を Jacobian-free Newton-Krylov (JFNK) 法を用いて解く。

   • 外側ループ: ニュートン反復によって解を更新する。
   • 内側ループ: クライロフ部分空間法（具体的にはGMRES）を用いて、線形システム $J \delta u = -F(u)$ を解く。
   • 行列フリー・アプローチ: この手法は、ヤコビ行列を明示的に構築しない。代わりに、ベクトル $\phi$ に対するフレシェ微分 $DT[u](\phi)$ の作用を計算する。
   • 微分近似: 著者らは、フレシェ微分の2つのアプローチを比較検討している：
     • 解析的（Analytic）: 写像 $T$ の微分の正確な閉形式の式を導出する。
     • 有限差分（Finite-Difference）: 中心差分を用いて微分を近似する。

3. 離散化:
   非局所項（畳み込み）は、周期領域における高速フーリエ変換（FFT）または適切な求積則を用いて処理される。これにより、精度の負担は微分ステンスルから、畳み込みと積分規則の精度へと移行する。

主要な貢献

• フレームワークの開発: 著者らは、マッケーン・ヴラスコフ方程式、クッカー・スマレス・フロッキング・モデル、およびニューラル・フォッカー・プランク方程式という3つの代表的なモデルに対して、堅牢な不動点フレームワークを開発し、検証した。
• 解析的微分と離散的微分の比較: 本論文は、解析的なフレシェ微分を使用する場合と、有限差分近似を使用する場合の直接的な比較を提供している。解析的なアプローチは、保証された精度と計算効率を提供する一方で、有限差分アプローチも解析的な形式が困難な場合には実行可能であることを示している。
• 分岐解析: 本手法は、既知の分岐図を正確に再現し、安定および不安定な定常状態の両方（超臨界および亜臨界の分岐を含む）を検出することに成功している。
• 新たな挙動の発見: 本フレームワークは、特にマルチモード・カーネルや2次元領域において、これまでこれらの問題では観察されていなかった複雑な分岐構造を明らかにする。

結果
著者らは、既存の解析結果および既知の分岐理論に対して自身のアプローチを検証している：

• 精度: クルカモ・型のカーネルを持つマッケーン・ヴラスコフ方程式において、数値解は解析的な参照解に対して誤差減衰率 $O(N_n^{-4})$（ここで $N_n$ は格子点数）で収束する。臨界分岐閾値は高い精度（$\sim 10^{-5}$）で回収される。
• マルチモード・カーネル: バイモーダルな相互作用カーネルを用いることで、本手法は、主要な分岐点から派生する二次的な分岐を含む、複数の解の枝を特定する。線形作用素が2次元の核（カーネル）を持つ場合（退化した分岐）、本手法は標準的な局所理論が特徴付けに苦慮する超越的（transcritical）な挙動を検出する。
• 安定性: 不動点ソルバーと摂動解析のための時間依存型ソルバーを組み合わせることで、各枝の安定性をマッピングしている。超臨界の枝は一般に安定であるが、亜臨界の枝は不安定であり、異なるプロファイルを持つ枝同士の交差などの複雑な相互作用が発生することを観察している。
• 汎用性: 本アプローチは以下に成功裏に拡張されている：
  • 2次元領域: トーラス $\mathbb{T}^2$ 上での適用可能性を示す。
  • 異なる領域: 無限領域（クッカー・スマレス）や、無流束境界を持つ半無限領域（ニューラル・フォッカー・プランク）の扱い。
  • 様々なカーネル: トップハット、三角形、および吸引・反発カーネルを用いたテストにより、二次分岐を含む豊かな解構造を明らかにしている。

意義
本論文は、この不動点アプローチが非局所フォッカー・プランク方程式の定常解を探索するための優れたツールであることを主張している。その主な意義は、**ダイナミカルな安定性に依存しない（agnosticism to dynamical stability）**点にあり、時間発展型ソルバーが見逃してしまう不安定な状態や複雑な分岐枝の検出を可能にしている。偏微分方程式の構造的特性を利用して問題を不動点写像へと再構成することで、微分作用子を直接離散化することに伴う困難を回避している。著者らは、現在の研究は線形拡散に焦点を当てているものの、必要な構造的条件が満たされる限り、本フレームワークは非線形拡散（例：ポロウス・ミディアム）や相互作用する個体群の系へと拡張可能であると述べている。