Rates of convergence of finite element approximations of second-order mean field games with nondifferentiable Hamiltonians

この論文「Rates of convergence of finite element approximations of second-order mean field games with nondifferentiable Hamiltonians（非微分可能なハミルトニアンを持つ 2 次メーンフィールドゲームの有限要素近似の収束率）」は、Yohance A. P. Osborne と Iain Smears によって執筆されたものです。以下に、論文の技術的な要約を問題設定、手法、主要な貢献、結果、そして意義の観点から日本語で詳細に記述します。

1. 問題設定 (Problem)

• 対象: 一般の有界多面体リプシッツ領域 $\Omega \subset \mathbb{R}^d$ 上で定義された、定常状態の 2 次メーンフィールドゲーム（MFG）システム。
• 特徴: ハミルトニアン $H(x, p)$ が非微分可能である場合を扱う。
  • 従来の MFG の数値解析では、ハミルトニアンの微分可能性が仮定されることが多く、これは Kolmogorov-Fokker-Planck (KFP) 方程式の移流項に現れるためである。
  • しかし、最小退出時間問題やバン・バン（bang-bang）制御などの多くの応用において、ハミルトニアンは非微分可能であり、最適フィードバック制御が一意でない場合がある。
• 定式化: 非微分可能性を扱うため、KFP 方程式は偏微分包含（PDI: Partial Differential Inclusion）として定式化される。
  $$\begin{aligned} -\nu\Delta u + H(x, \nabla u) &= F[m](x) \quad \text{in } \Omega, \\ -\nu\Delta m - G(x) &\in \text{div}(m \partial_p H(x, \nabla u)) \quad \text{in } \Omega, \end{aligned}$$
  ここで、$u$ は価値関数、$m$ はプレイヤー密度、$\partial_p H$ はハミルトニアンの偏微分（部分微分）である。
• 課題: 非微分可能なハミルトニアンを持つ MFG PDI に対する有限要素法（FEM）の近似解の**収束率（エラーの定量的評価）**を確立すること。従来の研究では収束性（定性的な結果）は示されていたが、具体的な収束率は不明であった。

2. 手法 (Methodology)

著者らは、以下の 3 つのステップを組み合わせた解析手法を採用している。

1. 正則化アプローチ (Regularization):

   • 非微分可能なハミルトニアン $H$ を、Moreau-Yosida 正則化 $H_\lambda$ で近似する。$H_\lambda$ は $C^{1,1}$ 級（微分可能で微分がリプシッツ連続）であり、$\lambda \to 0$ で $H$ に収束する。
   • これにより、PDI ではなく通常の偏微分方程式（PDE）システム（正則化された MFG）が得られる。
   • 正則化パラメータ $\lambda$ とメッシュサイズ $h$ を独立に制御する。

2. 安定化有限要素法 (Stabilized FEM):

   • 正則化された PDE システムおよび元の PDI システムに対して、安定化項を含む有限要素法を適用する。
   • 離散最大値原理（DMP）を満たすような安定化項 $D_k$ を選択し、数値解の非負性や安定性を保証する。
   • 価値関数 $u$ に対して $H^1$ ノルム、密度 $m$ に対して $L^2$ ノルムでの誤差評価を行う。

3. 誤差の三角不等式による分解:

   • 真の解 $(u, m)$ と離散解 $(u_k, m_k)$ の誤差を、以下の 3 つの項に分解して評価する。
     1. 正則化誤差: 真の解 $(u, m)$ と正則化された連続解 $(u_\lambda, m_\lambda)$ の差。
     2. 離散化誤差: 正則化された連続解 $(u_\lambda, m_\lambda)$ と正則化された離散解 $(u_{k,\lambda}, m_{k,\lambda})$ の差。
     3. 正則化と離散化の整合性誤差: 正則化された離散解 $(u_{k,\lambda}, m_{k,\lambda})$ と元の PDI の離散解 $(u_k, m_k)$ の差。

3. 主要な貢献と技術的工夫 (Key Contributions & Technical Innovations)

• 非微分ハミルトニアンの定量的解析:
  • ハミルトニアンが非微分可能であるため、部分微分集合 $\partial_p H$ からの可測選択間の定量的制御が困難である。著者らは、正則化されたハミルトニアンの微分 $\frac{\partial H_\lambda}{\partial p}$ のリプシッツ定数が $\lambda^{-1}$ として発散する性質を考慮しつつ、誤差評価を導出した。
• 半滑らか性（Semismoothness）の活用:
  • ハミルトニアンの非微分性を克服するため、半滑らか性（Lemma 2.1）を利用し、価値関数の誤差を密度関数の誤差と関連付ける補題（Lemma 3.1）を証明した。
• 双対性アプローチと密度の正則性:
  • 密度 $m$ の正則性が低い（$H^1$ 以下の場合がある）という困難に対処するため、双対性議論（Aubin-Nitsche 型のアプローチ）を工夫し、$L^2$ ノルムでの誤差評価を確立した。
• パラメータの最適化:
  • 正則化パラメータ $\lambda_k$ とメッシュサイズ $h_k$ の関係を最適化することで、全体の収束率を最大化した。具体的には $\lambda_k \sim h_k^{2\gamma/3}$ と選ぶことで、最適なバランスを得ている。

4. 主要な結果 (Results)

• 収束率の定理 (Theorem 5.1):
  • 適切な仮定（安定化条件 H1, H2、強単調な結合項など）の下で、有限要素近似解 $(u_k, m_k)$ と真の解 $(u, m)$ の誤差は以下のように評価される。
    $$\|u - u_k\|_{H^1(\Omega)} + \|m - m_k\|_{L^2(\Omega)} \lesssim h_k^{\gamma/3}$$
  • ここで、$\gamma \in (0, 1]$ は領域 $\Omega$ におけるポアソン方程式の楕円正則性（$L^2$ 源項に対する解の $H^{1+\gamma}$ 正則性）に依存する定数である。
    • 凸領域の場合 $\gamma=1$ となり、収束率は $h^{1/3}$ となる。
    • 非凸多角形領域では $\gamma < 1$ となり、収束率はそれより遅くなる。
• 正則化された離散解の誤差 (Corollary 5.2):
  • 実際の計算では、非微分な包含問題よりも正則化された PDE 系（(4.13)）を解く方が容易である。
  • 正則化された離散解 $(u_{k,\lambda_k}, m_{k,\lambda_k})$ に対しても、誤差評価 $\lesssim \lambda_k^{1/2} + \lambda_k^{-1} h_k^\gamma$ が成立し、$\lambda_k$ を適切に選べば同様の収束率が得られることを示した。

5. 意義 (Significance)

• 理論的進展:
  • 非微分可能なハミルトニアンを持つ MFG 問題に対して、初めて定量的な収束率（エラーバウンド）が証明された。これにより、数値解の精度を理論的に保証できるようになった。
• 実用的価値:
  • 非微分問題の直接解法は数値的に困難（不連続性のため）であることが多い。本研究は、正則化を用いた実用的な計算手法が、理論的な誤差制御の下で有効であることを示し、実装における指針を提供している。
• 応用範囲の拡大:
  • バン・バン制御や最小時間問題など、ハミルトニアンが非微分となる重要な応用分野における数値解析の基盤を築いた。
• 手法の一般性:
  • 安定化有限要素法と正則化手法を組み合わせたアプローチは、他の非滑らかな偏微分包含問題の解析にも応用可能な可能性を示唆している。

要約すると、この論文は非微分可能な MFG 問題の数値解析において、正則化手法と安定化有限要素法を巧みに組み合わせることで、理論的な収束率を初めて導出した画期的な研究である。