Bjorken sum rule down to the elastic region from the lattice

技術要約：ファインマン・ヘルマン法によるビョルケン総和則の格子決定

問題提起
ビョルケン総和則は、陽子と中性子の偏極構造関数（$g_1$）の第1モーメントの差として定義され、核子のスピン構造の理解および摂動論的量子色力学（pQCD）の検証における基本的な量である。この量は、$Q^2 \approx 1 \text{ GeV}^2$ まで実験的に検証されており、最近の実験によって弾性領域（$Q^2 \to 0$）へとマッピングされているが、摂動論的領域と非摂動論的な弾性／ゲラシュモフ・ドレヘーン（GDH）領域との間の理論的な遷移については、依然として十分に理解されていない。

高次ツイスト寄与（$1/Q^2$ 次の冪展開項）の直接的な格子QCD計算は、低次元演算子との演算子混合に起因する冪発散的な混合のため、極めて困難であることが知られている。したがって、演算子の繰り込みに伴う複雑さを回避しつつ、これらのモーメントや関連する高次ツイスト効果にアクセスするための代替的な第一原理的手法が必要とされている。さらに、ビョルケン総和則は、あらゆる運動量スケールにおける強結合定数 $\alpha_s$ の非摂動論的な定義への潜在的な経路を提供している。

手法
著者らは、四点相関関数の直接計算に伴う計算コストの増大とノイズを回避するため、ファインマン・ヘルマン（FH）定理を用いて偏極前方コンプトン振幅を計算している。

1. 理論的枠組み:

   • 偏福コンプト振幅は、分散関係を介して構造関数 $g_1$ と関連付けられる。格子上では、アクセスは物理的に不自然な領域 $|\omega| < 1$ に制限されており、振幅は $g_1$ のコーンウォル・ノートン・モーメントを用いて展開される。最低次のモーメント（$n=1$）は、直接的にビョルケン総和則に対応する。
   • ベクトルカレントに結合する、外部からの摂動場がQCD作用に導入される。これらの結合に対する核子のエネルギーの二階微分が、コンプト振幅を与える。
   • スピン依存部分を孤立させるため、著者らは背景場が存在する場合の、スピンアップ状態とスピンダウン状態の核子のエネルギー分裂を解析している。

2. 格子設定:

   • 計算は、QCDSFコラボレーションによって生成された、$2+1$ フレーバー（stout-smearingを用いた非摂動的な $O(a)$ 改良ウィルソン・フェルミオン）の単一ゲージアンサンブルを用いて行われた。
   • アンサンブルは、格子間隔 $a \approx 0.068 \text{ fm}$、体積 $48^3 \times 96$ であり、$m_\pi \approx 412 \text{ MeV}$ のピオン質量において $SU(3)$ フレーバー対称点に調整されている。
   • 核子は静止状態（$p=0$）に保たれ、運動量転移 $q$ は $0.5 \lesssim Q^2 \lesssim 8 \text{ GeV}^2$ の範囲をカバーするように変化させている。

3. 技術的革新:

   • 対角的 vs 非対角的摂動: 著者らは、摂動がスピン空間において対角的である場合（単純な比による抽出が可能）と、非対角的である場合（スピン混合を導入する）の両方に対処している。後者の場合、相関関数行列を対角化してエネルギー固有値を抽出するために、**一般化固有値問題（GEVP）**を実装している。
   • LOPE改良: 回転対称性の破れに起因する離散化誤差を軽減するため、著者らは格子特有の三角関数形式を用いて運動学的因子を置き換える、格子演算子積展開（LOPE）改良を適用している。
   • 弾性項の減算: コンプト振幅には弾性と非弾性の両方の寄与が含まれるため、電磁形因子（三点関数から導出）を用いて弾性ポール寄与を計算し、それらを減算することで、総和則に必要な非弾性モーメントを孤立させている。

主要な貢献と結果

• 第一原理計算: 本論文は、広範な $Q^2$ 範囲（$0.5 \text{ to } 8 \text{ GeV}^2$）にわたるビョルケン総和則のモーメントの格子決定を提示している。これは、摂動論的領域と非摂動論的領域の両方をカバーしている。
• 手法の妥当性確認: 本研究は、非対角的なスピン摂動を扱うためにGEVPを用いたFH手法を用いることが、直接的な四点関数計算に対する補完的な枠組みとなることを成功裏に実証した。
• 理論との比較:
  • 高 $Q^2$ ($> 4 \text{ GeV}^2$): 計算されたモーメントは、摂動論的量子色力学の予測（$\text{N}^3\text{LO}$ 精度の精度）およびクォーク・パートン模型と定性的な一致を示しているが、離散化アーティファクトによる系統的なオフセットが観察されている。
  • 中 $Q^2$ ($0.5 - 4 \text{ GeV}^2$): 格子の結果は、$Q^2$ が減少するにつれて純粋な摂動論的予測から逸脱する。計算された弾性寄与を加えると、格子データは理論曲線に近づくが、残留する不一致が残っている。
• 高次ツイスト効果: 摂動論的予測と（弾性減算後の）格子モーメントとの間の残留不一致は、高次ツイスト寄与の信号として解釈されている。しかし、著者らは、大きな $Q^2$ における残留する離散化アーティファクトのため、これらの高次ツイスト項を厳密に定量的に抽出することは現時点では困難であると述べている。

意義と主張
本論文は、ファインマン・ヘルマン法を用いて、非シングレット深非弾性散乱（DIS）のモーメントを格子上で決定することが可能であることを示している。主な意義は以下の通りである：

1. 冪発散的な演算子混合の問題なしに、高次ツイスト行列要素にアクセスするための第一原理的なアプローチを提供すること。
2. ビョルケン総和則における摂動論的領域と非摂動論的領域の間の遷移を研究するための経路を提供すること。
3. 連続極限およびカイラル極限の計算が行われ、かつ高次ツイスト効果が制御されることを前提として、ハドロン観測量から強結合定数 $\alpha_s$ を決定するための基礎を築くこと。

著者らは、現在の精度については謙虚な姿勢を維持しており、単一の格子間隔と重いピオン質量が、高次ツイスト効果と離散化誤差を明確に分離することを制限していることを認めている。彼らは、$1/Q^2$ の冪展開項を信頼性高く推定し、$\alpha_s$ を直接決定するためには、複数の格子間隔とより軽いクォーク質量を用いたさらなる研究が必要であると結論付けている。