Diagonal Kenney-Laub Rational Approximation to the Overlap Operator using Wilson and Brillouin Kernel

技術要約：オーバーラップ演算子に対する対角型ケニー・ローブ有理近似

問題設定
オーバーラップ・ディラック演算子は、連続理論のカイラル対称性の性質を保持し、ギンスパルグ・ウィルソン（GW）関係に密接に従う格子ディラック演算子の離散化である。しかし、その実用的な実装は、エルミート核演算子 $X$ の行列符号関数 $\text{sgn}[X]$ を近似する必要があるため、主に大きな計算上の課題に直面している。チェビシェフ多項式近似や単一分数有理近似などの従来の手法は、しばしば高い計算コストを要し、チューニングのために事前のスペクトル情報（固有値推定値）を必要としたり、近似次数の選択を困難にする振動的な収束挙動を示したりすることがある。さらに、核演算子の選択（例：ウィルソン型か改良型か）は、$X$ の条件数、ひいては符号関数の近似効率に大きな影響を与える。

手法
著者らは、以下の3つの特定のコンポーネントを統合したオーバーラップ演算子の定式化を提案している：

1. 対角型ケニー・ローブ（KL）イテレート： 標準的な多項式近似の代わりに、著者らは行列符号関数の近似のための、大域的に収束する有理反復スキームである対角型 $(n, n)$ KLイテレートを採用している。この手法の主要な特徴は、核演算子の極端な固有値の事前推定を必要としないことである。
2. 部分分数（PF）展開： 対角型KLイテレートは、部分分数展開を介して表現される。これにより、近似は $(X^2 + \sigma_i)^{-1}$ の形式の項の和へと変換される。この形式により、すべての項の同時反転にマルチシフト共役勾配（MSCG）アルゴリズムを用いることが可能となり、単一の高次多項式分数を反転する場合と比較して、計算オーバーヘッドが大幅に削減される。
3. 核の選択： 本研究では、標準的なウィルソン・ディラック演算子とブリリュアン演算子の2種類の核演算子を評価している。ブリリュアン演算子は、パラメータフリーで改良されたウィルソン型の演算子であり、理想的なギンスパルグ・ウィルソン・フェルミオンに近いスペクトルを示すことで知られており、潜在的により優れたカイラル特性を提供できる。

著者らは、この定式化をクエンチド格子QCD設定（$\beta=6.2$, $48 \times 24^3$ 格子）において実装し、確立されたチェビシェフ多項式近似と比較・検証している。彼らは、GW関係の破れ、符号関数の収束、および物理的観測量、具体的にはPCACクォーク質量と臨界裸クォーク質量（$m_{bare}^{critical}$）を分析している。

主な貢献

• アルゴリズムの実装： 著者らは、部分分数展開を用いた対角型KLイテレートの実用的な実装を実証しており、これによりスペクトル推定を必要としない堅牢なMSCGソルバーの使用が可能になる。
• 系統的な改善： 著者らは、対角型KL次数 $n$ を増やすことで、符号関数の近似およびGW関係の充足が、系統的かつ単調に改善されることを示している。
• 核の比較： 本研究は、KLフレームワーク内でのウィルソン核とブリリュアン核の詳細な比較を提供し、ブリリュアン核の優れた条件数とその収束率への影響を明らかにしている。
• ベンチマーク： KL-PFアプローチとチェビシェフ多項式法との包括的な比較を行い、符号関数の収束および特定の物理的精度目標に達するために必要な計算コストを評価している。

結果

• 単一分数に対する部分分数の効率性： PF展開は、単一分数（SF）形式よりも圧倒的に優れていることが証明された。高次（例：$n=5$）の場合、SFアプローチは極めて高価になる（ウィルソンの場合、384コア時間を超えるコスト）が、PFアプローチは効率的なままである（約0.22コア時間）。
• 符号関数の収束： KL次数 $n$ の増加は、「符号二乗の破れ」（$\delta_{sgn^2}$）の単調な減少をもたらす。ブリリュアン核は、特に条件数が小さい構成においては、ウィルソン核よりも低い計算コストで与えられた破れの閾値を達成する。
• 物理的観測量： PCAC質量とパイ中間子有効質量は、$n$ が増加するにつれて急速に安定する。収束は実質的に $n \geq 3$ で到達する。臨界裸質量（$m_{bare}^{critical}$）は、$n$ の増加とともに単調にゼロに近づく。
• KL vs. チェビシェフ：
  • 収束挙動： KL法は目標精度に向かって単調な収束を示すが、チェビシェフ法は振動的な挙動を示すため、最適な次数の選択がより困難になる。
  • 計算コスト： $|m_{bare}^{critical}| \approx 4 \times 10^{-4}$ という参照精度に達するために、KL-ブリリュアンの組み合わせは、チェビシェフ-ウィルソンの組み合わせよりも約20%少ない計算コストを必要とする。KL-ブリリュアンの組み合わせは、テストされたすべての手法と核のペアの中で最も効率的である。

意義および主張
本論文は、対角型KL有理近似とブリリュアン核の組み合わせが、オーバーラップ演算子の実用的かつ効果的な実装を提供すると主張している。主な意義は、この手法の計算効率と予測可能性にある。部分分数展開を利用することで、著者らは高次有理近似のスケーリング問題を克服し、高次のKLイテレートを可能にしている。KL法の単調な収束は、振動的なチェビシェフ法と比較して、目標精度を達成するために必要なチューニングの労力を軽減する。さらに、結果は、ブリリュアン核が（1回あたりの適用コストは高いものの）、その優れたスペクトル特性により、KL近似と組み合わされた際に純粋な計算上の利点を提供することを示唆している。著者らは、この組み合わせがベクトル型ゲージ理論の格子QCD計算における実装に値すると結論付けている。