[論文レビュー] Efficient use of the Generalized Eigenvalue Problem

B. Blossier, Georg von Hippel+3 more|arXiv (Cornell University)|Aug 7, 2008

Quantum Chromodynamics and Particle Interactions参考文献 12被引用数 8

ひとこと要約

本稿では、格子ゲージ理論における一般化固有値問題（GEVP）の最適化された使用法を提案し、エネルギー準位および行列要素を、体系的に制御された指数的に抑制された誤差で抽出することを目的としている。適切に時間間隔 $ t_0 $ を選ぶことで、補正項が $ \mathcal{O}(e^{-\Delta E_n t}) $ のように減少することが保証され、$ t > 0.6\,\text{fm} $ において統計誤差が 1 MeV 未塔の HQET における B メソン質量および崩壊定数の高精度抽出が可能になる。この手法は、$ a = 0.07\,\text{fm} $ および $ 0.1\,\text{fm} $ の格子上で数値的に検証された。

ABSTRACT

We analyze the systematic errors made when using the generalized eigenvalue problem to extract energies and matrix elements in lattice gauge theory. Effective theories such as HQET are also discussed. Numerical results are shown for the extraction of ground-state and excited B-meson masses and the ground-state decay constant in the static approximation.

研究の動機と目的

• 格子ゲージ理論においてエネルギー準位および行列要素を抽出する際の一般化固有値問題（GEVP）における系統的誤差を分析・制御すること。
• GEVPにおけるエネルギー準位および行列要素への補正が時間とともに指数的に減少することを示し、励起状態の精度向上に寄与すること。
• 時間間隔 $ t_0 \geq t/2 $ をとることで、安定で解析的に取り扱いやすい補正が得られる、洗練されたGEVPアプローチを提案すること。特にHQETのような有効理論において有効である。
• 静的近似におけるBメソンの基底状態および励起状態の質量、および基底状態の崩壊定数について、この手法を数値的に検証すること。

提案手法

• GEVPは $ C(t)v_n(t,t_0) = \lambda_n(t,t_0)C(t_0)v_n(t,t_0) $ の形で定式化され、$ C(t) $ はインターオペーラーの相関行列である。
• 有効エネルギー準位は $ E_n^{\text{eff}}(t,t_0) = \frac{1}{a}\log\left(\frac{\lambda_n(t,t_0)}{\lambda_n(t+a,t_0)}\right) $ により抽出され、$ \mathcal{O}(e^{-\Delta E_n t}) $ の補正項とともに真のエネルギーに収束する。
• 行列要素は $ p_{0n}^{\text{eff}} = \frac{\sum_j \langle P(t_1)O_j(0)\rangle (v_n)_{j} }{\|C(t_2)v_n\|} \cdot \frac{\lambda_n(t_0 + t_2/2, t_0)}{\lambda_n(t_0 + t_1, t_0)} $ により計算され、便宜上 $ t_1 = t_2 = t $ とおく。
• 系統的補正は $ \varepsilon_n(t,t_0) = \alpha_N e^{-\Delta E_{N+1,1} t} $ としてモデル化され、$ \Delta E_{N+1,1} $ は次の状態へのギャップを表し、データにフィットされる。
• 数値的実装では、$ C(t_i) $ の主要固有ベクトルを用いて 8x8 行列から $ N \times N $ 相関行列に切り詰める。$ t_i \approx 0.2\,\text{fm} $ を想定。
• シミュレーションは $ 2L \times L^3 $ の格子上で実施され、$ L \approx 1.5\,\text{fm} $ であり、HYP2 および O(a)-改質ウィルソン作用法を用い、$ \beta = 6.0219 $ および $ 6.2885 $ で実行された。

実験結果

リサーチクエスチョン

• RQ1 一般化固有値問題（GEVP）における系統的誤差は時間間隔 $ t $ とともにどのようにスケーリングされ、解析的に制御可能か？
• RQ2 系統的誤差を最小限に抑えつつ数値的安定性を保つために、GEVPにおける $ t_0 $ の最適な選択肢は何か？
• RQ3 GEVPを用いてBメソンの崩壊定数のような行列要素を、制御された指数的減少補正とともに抽出できるか？
• RQ4 有効理論（例：HQET）において補正項はどのように振る舞い、フルQCDと同様に指数的抑制を示すか？

主な発見

• 有効エネルギー準位 $ E_n^{\text{eff}}(t,t_0) $ は、$ \varepsilon_n(t,t_0) = \mathcal{O}(e^{-\Delta E_n t}) $ の補正項とともに真のエネルギーに収束し、Lüscher および Wolff の理論的予想を裏付ける。
• $ t > 0.6\,\text{fm} $ および $ N \geq 4 $ の条件下で、基底状態エネルギーの系統的誤差は 1 MeV 未塔であり、$ t \leq 1\,\text{fm} $ では統計誤差が 3 MeV 未塔であった。
• 有効崩壊定数は $ N=5 $ の場合、$ t_0 \approx 0.4\,\text{fm} $ でプラトーを示し、励起状態補正は 1% 未塔、統計誤差は約 0.7% であった。
• 関数形 $ \alpha_N e^{-\Delta E_{N+1,1} t} $ は、$ t \approx 0.6\,\text{fm} $ まで極めて良好に補正項をフィットし、指数的減少モデルの妥当性を確認した。
• 補正係数 $ \alpha_N $ は理論で予測されたように $ t - t_0 $ とともに減少し、解析的フレームワークの妥当性が裏付けられた。
• この手法は二つの格子スケール（$ a = 0.07\,\text{fm} $ および $ 0.1\,\text{fm} $）で安定に機能し、一貫性のある結果と小さな統計揺らぎを示した。