Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}_2$ from lattice QCD at physical quark masses

本研究は、物理的なクォーク質量を持つ格子 QCD 計算を用いて、大運動量有効理論（LaMET）の枠組みにおいて陽子のアイソベクトル・ヘリシティ部分子分布関数（PDF）を NNLO 精度で導出するとともに、初めて$\overline{\rm MS}$Scheme においてねじれ 3 項のモーメント$\tilde{d}_2$を決定し、これにより陽子内のクォークが経験する平均的なカラーローレンツ力を評価したものである。

要約

🧊 1. 研究の目的：「陽子の回転」の謎を解く

まず、陽子は原子核の材料ですが、実は「止まっている石」ではなく、激しく回転しているものです。
昔の科学者たちは、「陽子の回転は、中に入っている『クォーク（陽子の部品）』の回転で説明できるはずだ」と思っていました。しかし、実験してみると、クォークの回転だけでは陽子の回転の 3 割程度しか説明できず、残りの 7 割は「どこから来たのか？」という大きな謎（「陽子のスピン問題」）が残りました。

この論文は、その**「残りの 7 割」を埋めるために、クォークがどの方向に回転しているか**を、より詳しく、より正確に計算しようとしたものです。

🔍 2. 使った道具：「ラティス QCD」という「超解像カメラ」

通常、クォークは光（レーザーなど）で直接見ることはできません。そこで、この研究チームは**「格子 QCD（ラティス QCD）」**という手法を使いました。

• アナロジー：
  想像してください。陽子の内部を、**「3 次元のドット絵（格子）」で表現しています。
  このドット絵の上で、クォークがどう動き、どう回転しているかを、スーパーコンピューターを使って「シミュレーション」**します。
  今回使ったドット絵は、非常に細かく（解像度が高い）、かつ、現実のクォークの重さ（質量）を正確に再現したものでした。これにより、以前よりもはるかに鮮明な「内部の動き」が見えるようになりました。

🚀 3. 重要な工夫：「高速で走らせる」ことで中身を見る

ここで大きな問題があります。格子の上で静止している陽子を見ると、中身（クォーク）がどう動いているかは見えないのです。
そこで、**「陽子を光の速さ近くまで加速（ブースト）させる」**というアイデアを使いました。

• アナロジー：
  止まっている電車の中を見ると、乗客の動きは分かりません。でも、電車がものすごい速さで走っている時は、窓から乗客の顔や動きがはっきり見えますよね？
  この研究では、陽子を**「時速 150 万 km（1.53 GeV）」という超高速で走らせ、その状態で中身をスキャンしました。
  これにより、静止している状態では見えない「クォークの回転（ヘリシティ）」の分布を、「準 PDF（クオシ-PDF）」**という形で捉えることに成功しました。

🧩 4. 2 つの大きな発見

この研究では、2 つの重要なことを明らかにしました。

① 回転の分布マップ（ヘリシティ PDF）の作成

「陽子の回転」を、クォークが持っているエネルギーの割合（x）ごとに分けて描いた**「分布マップ」**を作りました。

• 結果： 中程度のエネルギーを持つクォーク（x=0.25〜0.75）の回転が、これまでの予想よりも強く寄与していることが分かりました。
• 意味： 「残りの 7 割」の謎を解く鍵が、この「中程度のエネルギーを持つクォーク」にある可能性が高まりました。

② 「ねじれた」力の測定（ねじれ 3 成分のモーメント）

これは少し高度な話ですが、クォークが回転する時、単に回るだけでなく、「グルーオン（クォークを結びつける接着剤のようなもの）」と絡み合って、ねじれた力が働いています。
これを**「ねじれ 3 成分（Twist-3）」**と呼びます。

• 発見： この研究で初めて、この「ねじれた力」の平均値を、理論的に非常に正確な方法（MS 法）で計算しました。
• 結果： 驚いたことに、この「ねじれた力」の値はほぼゼロでした。
• 意味： クォークとグルーオンの絡み合いによる「ねじれ」は、予想以上に小さかったのです。これは、陽子の回転が「単純なクォークの回転」と「グルーオンの回転」で説明しやすいことを示唆しています。

🧱 5. 限界と未来：「端っこの部分」はまだ見えない

今回のシミュレーションは、陽子の回転の「中心部分（中程度のエネルギー）」は非常に鮮明に見えましたが、**「端っこの部分（エネルギーが極端に低い、または高いクォーク）」**はまだ完全には見えませんでした。

• アナロジー：
  高速道路の真ん中はよく見えますが、**「入口（ゼロ）」や「出口（1）」**の部分は、カメラの焦点が合いにくく、ぼやけています。
  この「端っこの部分」を補うために、研究者たちは「理論的なモデル」を使って推測し、全体像を完成させました。

🌟 まとめ：なぜこれが重要なのか？

この論文は、**「陽子という複雑な箱の回転の正体」**を、実験室（シミュレーション）から初めて、非常に高い精度で描き出した画期的な研究です。

• 何が分かった？
  • クォークの回転分布が、予想より中盤で強い。
  • 「ねじれた力」はほとんどない。
• 未来への影響：
  この結果は、将来建設される**「電子・イオン衝突型加速器（EIC）」**などの実験計画にとって、非常に重要な「基準（ゴールポスト）」となります。実験結果がこの計算と合致すれば、私たちは陽子の構造を完全に理解したことになります。

つまり、**「宇宙の物質の 99% を構成する陽子の、回転という『魂』の正体に、一歩近づいた」**と言える研究です。

技術概要

この論文は、格子量子色力学（Lattice QCD）を用いて、物理的なクォーク質量条件下で陽子のアイソベクトル・ヘリシティ・パートン分布関数（PDF）を高精度に計算し、さらにツイスト 3 のモーメント $\tilde{d}_2$ を初めて MS Scheme で抽出した研究です。以下に、問題提起、手法、主要な貢献、結果、および意義について詳細な技術的サマリーを記します。

1. 研究の背景と問題提起

• 陽子のスピン問題: 陽子のスピン構造を理解することは、QCD の基本的な課題です。特に、クォークのスピンが陽子のスピンにどの程度寄与するか（ヘリシティ PDF $g_1(x)$）は、Jaffe-Manohar 分解や Ji のゲージ不変な分解の検証に不可欠です。
• 計算の困難さ: ヘリシティ PDF は光円錐上の非局所演算子で定義されるため、ユークリッド時空の格子 QCD では直接計算できません。
• 既存の課題:
  • 従来の LaMET（Large Momentum Effective Theory）アプローチでは、高運動量状態が必要ですが、計算コストと系統誤差（幂次補正）のバランスが課題でした。
  • ツイスト 3 の分布関数 $g_T(x)$ やそのモーメント $\tilde{d}_2$（クォークが受ける平均的な色ローレンツ力に関連）の計算は、演算子の混合や発散の扱いが複雑で、MS Scheme での直接抽出は行われていませんでした。
  • 多くの計算は物理的なクォーク質量と連続極限の両方を同時に満たすことが困難でした。

2. 手法と計算設定

本研究は、以下の高度な手法と設定を組み合わせています。

• 格子設定:
  • HotQCD コラボレーションによって生成された $N_f=2+1$ の HISQ 作用を用いたアンサンブルを使用。
  • 格子間隔 $a = 0.076$ fm、物理的な軽クォーク・ストレンジクォーク質量（パイオン質量 $\approx 140$ MeV）。
  • 格子サイズ $64^3 \times 64$、全 350 構成。
  • 価クォークにはツリーレベルのタッドポール改善された Clover-Wilson 作用を使用。
• 運動量と演算子:
  • 陽子を $P_z = \{0, 0.25, 1.02, 1.53\}$ GeV にブーストした状態を計算。
  • 非局所演算子（Wilson 線を含む）を測定し、ヘリシティ（$\Gamma = i\gamma_5\gamma_z$）とツイスト 3 横方向（$\Gamma = i\gamma_5\gamma_x$）の行列要素を抽出。
  • 励起状態の汚染を除去するため、2 状態フィッティングと源 - 吸収器間の時間分離（$t_{sep}/a = 6, 8, 10, 12$）を適切に使用。
• 正規化とマッチング:
  • ハイブリッド正規化: Wilson 線の線形発散を除去するため、距離 $z$ が小さい領域では比スキーム、大きい領域では線形発散を減算するハイブリッドスキームを採用。
  • レノルマロン再総和（LRR）: 線形発散の減算に含まれるレノルマロン（IR 曖昧性）を、主値（PV） prescription を用いた再総和により正則化し、MS スキームへの整合的な変換を可能にしました。
  • OPE と RG 不変比: 行列要素の紫外発散を相殺するため、RG 不変比を構成し、OPE（演算子積展開）を用いてモーメントを抽出。
  • LaMET マッチング: 準 PDF（Quasi-PDF）を光円錐 PDF に変換するため、NNLO（次々次次リードオーダー）のマッチング核を使用し、RG 再総和と閾値再総和を適用して理論誤差を低減。
• SDF（Short-Distance Factorization）による端点領域の制約:
  • LaMET が $x \to 0, 1$ の端点領域で破綻するため、座標空間の行列要素を直接フィッティングする SDF 手法を用いて、端点領域のパラメータ（$x^a(1-x)^b$ の係数）を制約し、全 $x$ 領域の PDF を構築しました。

3. 主要な貢献と結果

A. ヘリシティ PDF のモーメントと $\tilde{d}_2$ の抽出

• モーメントの抽出: RG 不変比の OPE 解析により、ヘリシティ PDF のモーメント $\langle x^n \rangle$ を $\mu = 2$ GeV で抽出しました。
  • $\langle x \rangle^{u-d} = 0.291(7)$
  • $\langle x^2 \rangle^{u-d} = 0.101(10)$
  • これらの値は、JAM22 グローバルフィッティングの結果よりもやや大きい傾向を示しました。
• $\tilde{d}_2$ の初回抽出: ツイスト 3 のモーメント $\tilde{d}_2$ を、短距離行列要素の $\lambda$ 依存性から MS スキームで初めて直接抽出することに成功しました。
  • 結果: $\tilde{d}_2^{u-d}(2 \text{ GeV}) = 0.0024(46)$
  • この値は統計的にゼロと一致しており、アイソベクトルチャネルにおける Genuine twist-3 のクォーク - グルオン相関が強く抑制されていることを示唆しています。これは、局所演算子に基づく従来の格子計算や実験的制約と整合的です。

B. $x$ 依存性を持つヘリシティ PDF の決定

• 中 $x$ 領域 ($0.25 \le x \le 0.75$):
  • ハイブリッド正規化と再総和を施した準 PDF を Fourier 変換し、NNLO マッチングを適用して光円錐 PDF を得ました。
  • 結果、中 $x$ 領域（$x \in [0.4, 0.7]$）での分布は、既存のグローバルフィッティング（JAM22, BDSSV24）と比較してより高い値を示す傾向が見られました。この不一致は、より高い運動量や格子系統誤差のさらなる制御が必要であることを示唆しています。
• 端点領域のモデル化:
  • SDF 手法を用いて端点領域のパラメータをフィッティングし、LaMET 結果と結合して全 $x$ 領域の PDF を構成しました。
  • フィッティングにより得られたパラメータは、端点での振る舞いを物理的に制約しつつ、LaMET の予測を補完するものとなりました。

4. 意義と結論

• 物理的意義:
  • 物理的なクォーク質量と微細な格子間隔を用いた計算により、陽子のスピン構造に関する第一原理からの高精度な制約を提供しました。
  • ツイスト 3 のモーメント $\tilde{d}_2$ を MS スキームで初めて抽出し、クォーク - グルオン相関の大きさを定量的に評価しました。その結果、このモーメントは非常に小さいことが確認されました。
• 方法的な進展:
  • レノルマロン再総和を組み込んだハイブリッド正規化スキームの実践により、線形発散と IR 曖昧性を統一的に処理する枠組みを確立しました。
  • LaMET と SDF を組み合わせることで、LaMET の適用範囲（中 $x$ 領域）と SDF の強み（端点領域の制約）を補完し、より信頼性の高い全領域 PDF を構築する手法を提示しました。
• 将来展望:
  • 中 $x$ 領域でのグローバルフィッティングとの不一致を解消するため、より高いハドロン運動量（$P_z$）へのアクセスや、複数の格子間隔を用いた連続極限への外挿、励起状態のさらなる制御が必要であるとしています。

総じて、この論文は格子 QCD におけるハドロン内部構造の計算において、理論的精度（NNLO マッチング、レノルマロン処理）と物理的現実性（物理質量）を両立させた重要なマイルストーンであり、特にツイスト 3 物理の MS スキームでの直接抽出という点で画期的な成果です。