Scalar-Tracking SAV Schemes with Pullback Corrections for Gradient Flows

技術要約：勾配流のためのプルバック補正を用いたスカラー追跡型SAVスキーム

問題提起
$\partial_t \phi = G \frac{\delta E}{\delta \phi}$ によって支配される勾配流は、フェーズフィールド動力学、材料科学、および最適化における散逸進化の基礎的なモデルである。その数値解析における中心的な課題は、離散的なエネルギー散逸則（無条件エネルギー安定性）を維持しつつ、高い軌道精度と計算効率を両나める時間離散化を構築することである。Scalar Auxiliary Variable (SAV) 法は、非線形勾配流を、修正されたエネルギーを散逸させる線形陰的なスキームへと変換する枠組みとして普及してきた。しかし、標準的なSAV法では、すべての非線形力の情報を単一のスカラー補助変数に圧縮するため、状態方程式に対してランク1の補正が適用される。Multiple-SAV (MSAV) 法は、複数のスカラー変数を導入してエネルギー成分を追跡することで精度を向上させるが、これは補助系の次元を増大させる。著者らは、既存の手法における結合された限界を特定している。すなわち、成分ごとの補正（MSAVに見られるもの）によって得られる改良された計量情報が、スカラー補助変数の数と、状態方程式に適用される補正のランクと通常は結びついている点である。本論文は、単一のスカラー補助変数のみを使用して、成分ごとの補正による改良された計量情報を保持できるか否かを問うている。

手法
著者らは、Pullback-Corrected SAV (PB-SAV) ファミリーを導入する。核心となる革新は、スカラーエネルギー追跡メカニズムと、状態方程式に適用される補正のランクをデカップリング（分離）することにある。

1. 単一スカラー・トラッカー: PB-SAVは、単一のスカラー変数 $q(t)$ が全非線形エネルギー $E_1[\phi] + C$ を追跡するという標準的なSAVの定式化を維持する。これにより、修正エネルギー則とスカラーの簿記（bookkeeping）は同一となる。
2. プルバック補正: 全非線形力から導かれるランク1の補正（標準的なSAV）の代わりに、PB-SAVは、非線形エネルギーの許容可能な分解 $E_1[\phi] + C = \sum_{\alpha=1}^r (E_{1,\alpha}[\phi] + C_\alpha)$ によって誘導される高ランクの補正を採用する。
3. 数学的定式化:
   • $R(\phi) = (Q_1(\phi), \dots, Q_r(\phi))$ とし、$Q_\alpha = \sqrt{E_{1,\alpha} + C_\alpha}$ とする。全エネルギー・トラッカーは $Q(\phi) = \|R(\phi)\|_2$ である。
   • 標準的なSAVは、補正 $B_1^n = 2 DQ(\phi^n)^T DQ(\phi^n)$ を適用し、これはランク1である。
   • PB-SAVは、補正 $B_r^n = 2 DR(\phi^n)^T DR(\phi^n)$ を適用する。これは各成分に対応するランク1項の和として解釈される。これは、成分空間のユークリッド計量を、スカラーノルムへと収束させる前に状態空間へとプルバック（引き戻し）するものと解釈される。
4. ソルバーの実装: 結果として得られる状態方程式は、標準的な半陰的な作用素に対する低ランク摂動を含む。著者らは、Sherman–Morrison–Woodbury 公式を利用してシステムを効率的に解く。これは、時間ステップサイズ $\Delta t$ に依存せず、各タイムステップにつき基本線形ソルバーを $r+1$ 回適用するだけで済む。
5. ステップ依存の分解: PB-SAVの重要な特徴は、補正に使用される成分分解が、安定性を損なうことなく毎タイムステップ変化できることである。これは、スカラー変数が全エネルギーのみを追跡するためである。

主な貢献

• 構成と安定性: 著者らは、第1次PB-SAVスキームを構築し、それが標準的なSAVと同一の修正エネルギー散逸則を満たすことを証明した。これは、成分分解がステップごとに変化する場合でも成立する。
• 理論的比較: 彼らは、PB-SAVの補正行列 $B_r^n$ が、標準的なSAVの補正 $B_1^n$ に対して正定値順序において優位である（$B_r^n \succeq B_1^n$）ことを確立した。その差は、正規化された成分方向微分の重み付き分散として明示的に特徴付けられる。これは、成分分割を精緻化（ランクを増加）させると、補正の質が厳密に向上または維持されるという単調性の結果を与える。
• 最適化の解釈: 有限次元において、PB-SAVの補正は、非線形エネルギーの最小二乗表現のガウス・ニュートン行列として特定される。これにより、複数の補助変数を導入することなく、局所的なヘッセ行列構造をより良く近似することで、なぜ軌道精度を向上させることができるのかが明確になる。
• 緩和戦略: 補助変数 $q_n$ が定義値 $Q(\phi_n)$ から乖離すること（これが最適化の文脈で停滞を引き起こす可能性がある）に対処するため、著者らは、散逸則を維持しながら一貫性を保つエネルギー許容的な緩和策を提案している。

結果
有限次元の勾配流、Allen–Cahn動力学、および非局所Cahn–Hilliardモデルに関する数値実験は、2つの異なる性能領域を示している：

1. 誤差定数の減少: 滑らかな領域（例：Allen–Cahnの滑らかな初期条件）において、PB-SAVはSAVおよびMSAVと同じ第1次の時間収束率を示す。これらのケースでは、主な利点は誤差定数の減少であり、軌道精度自体は標準的なSAVと同程度である。
2. 大幅な軌道改善: 非線形エネルギーが異なるスケールや空間的位置を持つ明確な成分（例：グループ化された曲率を持つ有限次元問題、または特定のスペクトル内容を持つ非局所モデル）を持つ問題において、高ランクの補正は精度を大幅に向上させる。
   • 2つの曲率クラスを持つ凸二次形式のテストにおいて、ランク2のPB-SAVは、SAVと比較して目的関数のギャップを5分の1に減少させた。
   • 正則化最小二乗問題において、ランク2の補正はヘッセ行列を機械精度まで復元し、緩和されたSAVよりも数桁速い収束を実現した。
   • 非局所Cahn–Hilliardモデルでは、補正ランクを1から4に増やすことで、最終的な $L^2$ 誤差が約7.3倍減少し、活動的なスペクトル内容が捉えられるところで飽和した。
   • 適応ランク戦略が実証され、能動的なエネルギー成分に基づいてランクを動的に圧縮することで、固定の高ランクスキームの精度を、平均ランク $\approx 2$ で達成した。

意義
本論文は、PB-SAVが「スカラーエネルギー・トラッカーの数」と「補正のランク」をうまく分離していると主張している。これにより、実務者は、単一のスカラー補助変数、単純なスカラー修正エネルギー則、および時間依存しない線形ソルバーの低ランク摂動というSAVの実用的な利点を保持しながら、MSAVの優れた挙動の根源である成分ごとの計量情報を取り入れることが可能となる。この手法は、補正ランクを問題の幾何学的構造に合わせて調整または適応できる柔軟な枠組みを提供し、無条件のエネルギー安定性を維持しながら、特定の領域において大幅な精度向上を実現する。