Sign of Wilson-Fermion Determinant using Contour-Integral Spectral Projection

技術要約：輪郭積分スペクトル射影を用いたウィルソン・フェルミオン行列式の符号の決定

問題設定
ウィルソン・フェルミオンを用いた格子量子色力学（QCD）において、フェルミオン行列式の符号 $\text{sign}(\det D)$ は極めて重要な量である。$\gamma_5$-エルミート性により行列式は実数であることが保証されているが、単一フレーバーまたは質量非退化のウィルソン・フェルミオンの場合、それが正であることは保証されない。このシグネチャー（符号）は、非エルミート・ウィルソン・ディラック演算子（$D$）の負の実軸上に存在する固有値の数 $N_-$ のパリティによってのみ決定される。このシグネチャーの決定は計算量的に困難であり、特に低エネルギー領域にモードが存在する場合や、「例外的な構成（exceptional configurations）」のように演算子の条件数が悪化しているパラメータ領域において顕著である。スペクトルフロー法（質量関数としてのエルミート演算子 $Q=\gamma_5 D$ の固有値を追跡する手法）や、非エルミート行列に対する直接的な反復固有値ソルバーなどは、しばしば広範な低エネルギー・スペクトル領域を解像する必要があったり、蓄積された誤差を補正するための重み付け因子（$W_{\text{sign}}$）を必要とする多項式近似や有理近似に依存したりすることがある。

手法
著者らは、負の実固有値を直接カウントするために、低エネルギー・スペクトル全体を解する必要を回避し、負の実固有値を直接分離する、輪郭積分に基づく標的型のスペクトル射影法を提案している。この手法の核となる構成要素は以下の通りである：

1. 輪郭積分スペクトル射影：
   本手法は、関心の対象となるスペクトル領域を孤立させるために、実軸に対して対称な楕円形の複素アルガン図形内の輪径 $\Gamma$ を利用する。スペクトル射影体は次のように定義される：
   $$P_\Gamma = \frac{1}{2\pi i} \oint_\Gamma (zI - A)^{-1} dz$$
   ここで $A$ は非エルミート・ウィルソン・ディラック演算子である。この積分は、ノード数 $N_q$、ノード $z_j$、および重み $\omega_j$ を持つ有限の求積則を用いて数値的に近似される。得られる近似射影体 $\hat{P}_\Gamma$ は、演算子を縮小された部分空間へと写像する。

2. $\gamma_5$-エルミート性の活用：
   ウィルソン・ディラック演算子の $\gamma_5$-エルミート性により、左および右の固有系は互いに関連付けられている（$r_k^\dagger \gamma_5$ が $\lambda_k^*$ に対する左固有ベクトルとして機能する）。その結果、本手法は右射影部分空間の求解のみを必要とし、左の固有系を独立して決定する必要がない。

3. シフトされた線形系とソルバーの適応：
   射影プロセスには、各求積ノードにおいてシフトされた線形系 $(z_j I - \hat{D})q_j = Y$ （ここで $\hat{D}$ は $D$ の偶のシュア補行列）を解くことが含まれる。著者らは、もともとフル演算子 $D$ 用に設計された DFL-SAP-GCR ソルバーを、シフトされた偶のシュア補行列 $\hat{D} - zI$ に対してのみ作用するように適応させている。この適応は、シフトされたフル演算子のシュア補行列が、シフトされたシュア補行列と一致しないため、極めて重要である。

4. スペクトル境界のための特異値除外：
   負の実固有値をすべて包含しつつ、計算コストを増大させないように輪径の左端（$X_L$）を決定することが重要な課題である。著者らは、特異値推定を用いる実用的な戦略を導入している。具体的には、一連の点 $s_i$ におけるシフトされた系 $\hat{D} - s_i I$ の最小特異値 $\sigma_{\min}$ を評価することで、固有値を含まない除外区間を特定する。これにより、エルミート部分のレイリー商のみから導かれる過度に保守的な境界を避け、保守的かつ効率的な $X_L$ の決定が可能となる。

5. 部分空間反復：
   試行部分空間ベクトルの生成には FEAST パッケージが用いられる。試行部分空間の次元 $m$ は、包含されると予想される固有値の数を超えるように選択される。本手法は急速な収束を示し、通常、2回の FEAST 反復以内にすべての負の実モードを解明する。

主な結果
本手法は、CLS コンソーシアムによって生成された格子QCDアンサンブルの多様な構成に対してテストされた。主な知見は以下の通りである：

• 堅牢性と精度： 本手順は、負の実解を持たない構成、浅い負の実解を持つ構成、深い負の実解を持つ構成、および解のオーダーが大きく異なる複数の解を持つ構成を含む様々な構成において、負の実解をマシン精度まで正確に再現した。
• パラメータ効率： ほとんどの構成において、最小の求積ノード数 $N_q = 8$ およびシフトソルバーの許容誤差 $10^{-3}$ で十分であった。すべての負の実固有値は、2回目の FEAST 反復後に解明された。
• 極端なケースへの対応： 固有値が2桁の範囲にわたる構成（例：configuration C102r005n56）では、$N_q$ を 12 に増やし、ソルバーの許容誤差を $10^{-5}$ に厳格化することで、大きな計算オーバーヘッドなしにすべてのモードを忠実に捉えることができた。
• 計算コスト： 主なコストは、シフトされた線形系の求解から生じる。輪径および求積ノードの選択に伴うオーバーヘッドは、全評価の数パーセント未満である。特異値除外戦略により、積分領域が効果的に最小化され、エルミート部分のみから導かれる境界を使用する場合と比較して、必要な求積ノード数が削減された。

意義と主張
本論文は、この輪郭積分に基づくスペクトル射影が、望ましくない物理モードの解明や広範なスペクトルフロー解析に頼ることなく、フェルミオン行列式の符号を決定するための系統的かつ効率的な手段を提供すると主張している。固有値ソルバーを行列式のシグネチャーを支配する正確な領域に限定することで、本手法は、多項式変換を併用した従来の反復固有値ソルバーに代わる堅牢な選択肢を提供する。

格子QCDにおける $W_{\text{sign}}$ の評価という特定の適用範囲を超えて、著者らは、このアプローチが、高密度な複素スペクトルの中に埋め込まれた物理的に関連のある固有値を分離するための一般的な戦略を提供することを示唆している。これには、非エルミート量子系、安定性解析、および大規模な固有値問題への応用が考えられ、広範な探索ではなく、スペクトルの特定の部分集合に対する標的型のフィルタリングを提供できる。本研究は、行列式の符号を信頼性高く決定することが、シミュレーションにおける有理近似誤差の補正（重み付けによる）や、単一フレーバーにおける測度の正当性の検証において不可欠であることを強調している。