Evaluation of resonances: adaptivity and AAA rational approximation of randomly scalarized boundary integral resolvents

技術要約：適応型AAA有理近似による共鳴の評価

問題設定
本論文は、開いた（open）および閉じた（closed）空洞、あるいはその他の散乱構造における音響および電磁共鳴を評価するという計算上の課題に取り組んでいる。数学的に、これらの共鳴は、領域 $\Omega$ 上の斉次境界条件（ディリクレ、ノイマン、またはザレムバ条件）の下でのヘルムホルツ固有値問題 $\Delta u + k^2 u = 0$ の解 $(u, k)$ に対応する。

• 内部問題の場合、$k$ は実固有値を表す。
• 外部および開いた空洞の問題の場合、$k$ は複素平面の下半平面（$\Im k \le 0$）に位置する複素「散乱極」または「複素固有値」を表す。

この問題は、境界積分演算子の離散化から導出された行列値関数 $F(k)$ に関する非線形固有値問題（NEP）として定式化される。目標は、$F(k)w = 0$ を満たす非ゼロベクトル $w$ が存在するような複素数 $k$ を見つけることである。

手法
提案されるアルゴリズムは、主に3つのコンポーネント、すなわち、ランダムにスカラー化されたレゾルベント、AAA有理近似アルゴリズム、およびセカント法による精緻化段階を備えた適応的探索戦略を組み合わせたものである。

1. スカラー化レゾルベント: フルサイズの行列NEPを直接解く代わりに、本手法はスカラー関数 $S(k) = u^* F(k)^{-1} v$ を評価する。ここで、$u$ と $v$ は固定されたランダムなベクトルである。$S(k)$ の極は、元の問題の固有値に対応する。
2. AAA有理近似: 関数 $S(k)$ は、AAAアルゴリズムを用いて有理関数 $r(k)$ によって近似される。このアルゴリズムは、重心形式の有理近似式を構築するために、補間点と重みを貪欲に選択する。この有理関数 $r(k)$ の極が、共鳴の初期近似として機能する。
3. 適応的探索戦略: 探索領域内のすべての共鳴を確実に捉えるために、アルゴリズムは領域の再帰的な二分分割（複素平面内の矩形、または実直線上の区間）を採用している。
   • アルゴリズム1（基本）: 点の集合で $S(k)$ を評価し、AAA近似を構築し、領域内の極を返す。
   • アルゴリズム2（適応型）: 領域を再帰的に分割する。さらなる分割を行った際に、サブ領域内の極の数が変化した場合、プロセスを繰り返す。これにより、サンプリング密度の不足によって共鳴を見逃すことがないようにする。
4. 精緻化と検証:
   • セカント法: 初期極が特定された後、$1/S(k)$ に対してセカント法の終了段階を適用し、マシン精度に近いレベルまで指数関数的な収束を実現する。このステップは、信頼できる誤差推定も提供し、AAAアルゴリズムが時折生成する可能性のある偽の固有値を排除する役割も果たす。
   • 固有関数: 精緻化された固有値における $F(k)$ の零空間を計算することで、対応する固有関数を復元する。これにはピボット選択を伴うガウス消去法を用いる。

主な貢献

• 直接的な極の抽出: NEPを線形化するために有理近似を使用する従来のAAAベースのアプローチとは異なり、本手法はスカラー化されたレゾルベントの有理近似の極を直接、固有値の推定値として使用する。
• 幾何学的な柔軟性: 輪郭積分に依存するのではなく、有理近似に依存することで、矩形領域や実区間を含む任意の探索領域の使用が可能となる。これは、適応的分割戦略の設計を容易にする（これは、矩形に対する台形則のような従来の輪郭積分法では実装が困難である）。
• 堅牢性と効率性: 適応的探索とセカント精緻化の組み合わせにより、低周波数および高周波数状態の両方において高い精度（マシン精度に近い精度）を確保する。本手法は、様々な波数におけるスカラー化されたレゾルベントの評価のみを必要とし、これは並列化が極めて容易である。
• 開いた弧（Open Arcs）の取り扱い: 本論文は、単層ポテンシャルの法線微分の特定のジャンプ関係を利用して、エッジの特異点を処理し、開いた弧の問題へとレゾルベントの極と固有値を関連付ける理論的枠組みを拡張している。

数値結果
論文では、本手法の性能を示す様々な数値実験が提示されている。

• 一般的なNEP: アルゴリズムは、NLEVPコレクションの「CDプレーヤー」（60個の固有値）および「バタフライ」（256個の固有値）問題からすべての固有値を正常に特定した。
• 低周波散乱: 開口部を持つ凧形および円形の空洞における固有値の正確な計算。本手法は、実固有値（内部問題）と実軸付近の複素極（開いた空洞）の両方を扱った。
• 輪郭積分との比較: ブロックSS法およびBeynの「積分アルゴリズム1」との比較において、Beynの手法は精緻化なしでも高速に収束できるが、提案手法は局所的な精緻化によって同等の精度を達成することを示した。決定的な違いとして、提案手法は、ブロックSS/Beynにおいて $L > 1$ の場合に必要となる、周波数点ごとの複数回の行列計算を回避できるため、LU分解が利用できない大規模な3D問題において潜在的な優位性を持つ。
• ギャップサイズの漸近挙動: 開いた弧が閉じた曲線に近づく際の散乱極の収束の研究に使用された。結果として、減衰率（$\Im k$）は、モードのギャップとの整列に応じて $O(\theta^4)$ または $O(\theta^8)$ の速度で収束し、空間周波数（$\Re k$）は $O(\theta^2)$ または $O(\theta^4)$ で収束することが示された。
• 高周波性能: アルゴリズムは、単位円の区間 $[1, 100]$ において1,244個の異なる内部固有値をマシン精度に近い精度で計算することに成功した。また、キートやロケット型の領域の高周波固有関数を可視化し、開いた空洞の固有関数と平面波照射下での散乱解との間の強い対応関係を示した。

意義と主張
著者らは、提案されたアプローチが単純で実装が容易であり、急速に収束すると主張している。その主要な意義は、輪郭積分法によって課される幾何学的制約を受けることなく、幅広い周波数範囲にわたって、散乱共鳴およびNEPの解に対して非常に高精度な近似を提供できる能力にある。適応的戦略は網羅的な探索能力を保証し、セカント精緻化段階は高い精度と信頼性を保証する。本論文は、本手法を、特に複雑な幾何学形状を含む問題や、繰り返しの行列分解に計算リソースが制限される問題において、既存の輪郭積分ベースの固有値ソルバーに対する柔軟な代替手段として位置づけている。