Fully discrete approximation of the semilinear stochastic wave equation on the sphere

技術要約：球面上における半線形確率波動方程式の完全離散近似

問題設定
本論文は、乗法的ガウスノイズによって駆動される単位球面 $S^2$ 上の半線形確率波動方程式に関する数値解析に取り組んでいる。方程式は以下のように与えられる：
$$\partial_{tt}u(t) = \Delta_{S^2}u(t) + f(u(t)) + g(u(t))\dot{W}(t),$$
ここで $\Delta_{S^2}$ は球面ラプラシアン（ラプラス・ベルトラミ演算子）であり、$f$ と $g$ は大域リプシッツ連続な非線形関数、$W$ は球面上の等方的な $Q$-ウィーナー過程である。著者らは、ユークリッド領域におけるSPDE（確率偏微分方程式）の数値解析は確立されているものの、曲面、特に球面における時間依存のSPDEに関する研究は限られていると指摘している。既存の文献は主に線形方程式または加法的ノイズに焦点を当てているが、本研究はより複雑な、乗法的ノイズを伴う半線形方程式を対象としている。

手法
著者らは、主に2つの近似技術を組み合わせた完全離散数値スキームを提案している：

1. 空間離散化（スペクトル・ガラーキン法）：
   解は、球面調和関数 $Y_{\ell,m}$ に基づくスペクトル・ガラーキン法を用いて近似される。解の無限級数展開は、有限モード $\kappa$ で打ち切られる。これにより、問題は最初の $\kappa$ 個の球面調和関数によって張られる有限次元部分空間へと投影される。射影演算子 $P_\kappa$ はラプラス・ベルトラミ演算子と可換であり、元の確率発展方程式の構造を保持する半離散定式化を可能にする。

2. 時間離散化（確率三角関数積分器）：
   時間ステップにおいては、著者らは明示的な確率三角関数積分器（確率指数積分器とも呼ばれる）を採用している。このスキームは、方程式の線形部分によって生成される厳密な半群 $E(h)$ を用いて、状態ベクトル $X_n = (u_n, v_n)^T$ を更新する：
   $$X^\kappa_n = E(h)X^\kappa_{n-1} + hE(h)P_\kappa F(X^\kappa_{n-1}) + E(h)P_\kappa G(X^\kappa_{n-1})\Delta W_{n-1}.$$
   このアプローチは、球面上の波動方程式に関連する厳密な余弦および正弦演算子を利用しており、標準的なオイラー・マルヤマ・スキームで見られるような分裂誤差を回避している。

主要な貢献と理論的結果
論文では、厳密な解のウェルポーズドネス（適切性）を確立し、提案された完全離散スキームの厳密な収束率を導出している。

• ウェルポーズドネスと正則性： 著者らは、球面上の抽象的な確率発展方程式の厳密解の存在と一意性を証明している。彼らは、初期データの正則性とノイズの滑らかさ（共分散演算子 $Q$ によって特徴付けられる）に依存する、ソボレフ空間 $H^\gamma(S^2)$ における解のモーメント有界性と正則性の評価を確立している。
• 空間収束（定理9）： スペクトル・ガラーキン近似の強収束率を証明している。誤差界は $\kappa^{-\min(1, \beta, \delta)}$ のオーダーであることが示されており、ここで $\beta$ は初期データの正則性を、$\delta$ はノイズの滑らかさを表す。ほとんど確実に（almost sure）の収束率も導出されている。
• 時間収束（定理11）： 明示的な確率三角関数積分器の強収束率を証明している。半離散解と完全離散解の間の誤差は $O(h^{\min(1, \beta, \delta)})$ で抑えられる。
• 完全離散収束（系12）： 三角不等式を介して空間および時間の解析結果を組み合わせることで、離散化パラメータが適切にバランスされている場合、完全離散スキームが強収束率 $O(N^{-\min(1, \beta, \delta)})$ （ここで $N$ は時間ステップ数およびモード数）を達成することを実証している。

数値実験
理論的知見は、加法的および乗法的ノイズの両方のケースをカバーする数値実験を通じて検証されている。

• 加法的ノイズ： 実験は、異なるノイズ滑らかさパラメータに対応する予測された空間収束率（$1, 1/2, 1/4$）を裏付けている。時間収束テストでは、提案された三角関数積分器を半陰的なオイラー・マルヤマ・スキームと比較している。三角関数積分器は、滑らかなノイズ（$\delta=1$）に対しては一貫して最適なレート1を、ノースムースなノイズ（$\delta=1/2$）に対しては$1/2$を達成し、後者のケースにおいてオイラー・マルヤマ・スキームを凌駕している。
• 乗法的ノイズ： 非線形として $g(u) = \sin(u)$ および $g(u) = (1+u)/(1+u^2)$ を用いてスキームのテストが行われた。結果は、理論的な強収束率およびパスワイズ収束率を裏付けており、乗法的ノイズのシナリオに対する本手法の堅牢性を実証している。

意義と主張
著者らは、本研究を、乗法的ノイズによって駆動される球面上の半線形確率波動関数に対する、完全離散法の最初の数値解析として位置づけている。先行研究では、曲面上の線形方程式や加法的ノイズについては扱われてきたが、本論文はより一般的な半線形乗法的ケースへと理論的枠組みを拡張している。主な意義は以下の通りである：

1. この特定のクラスの問題に対する、最初の強収束およびほとんど確実な収束の解析を提供したこと。
2. 明示的な確率三角関数積分器をスペクトル法と組み合わせることが、球面上における波動方程式の幾何学的構造を保持しつつ、最適な収束率を達成するために有効であることを示したこと。
3. ユークリッド設定を超え、線形および加法的制約を越えた、多様体上のSPDEの数値的取り扱いに関する文献の空白を埋めたこと。

論文は、提案された手法が理論的に健全であり、数値的に効率的であり、実験によって導出された境界と一致する収束率が確認されたと結論付けている。