Thermal axion production in QCD from the lattice

技術要約：格子QCDによる熱的アクシオン生成

問題提起
初期宇宙における熱的アクシオンの生成率は、光錐上で評価されるQCDトポロジカル電荷密度スペクトル関数 $\sigma(k, k)$ によって決定される。摂動論的手法は非常に高温では有効であり、カイラル摂動論は転移領域以下で適用可能であるが、クロスオーバー領域（$T \sim 150\text{--}300$ MeV）は強結合であるため、これらの手法は信頼性に欠ける。標準的な格子QCDのアプローチは、ミンコフスキー・スペクトル関数を必要とする一方で、それらはユークリッド相関関数を計算するという根本的な障害に直面している。このギャップを埋めるための過去の試みは解析接続に依存していたが、これは逆問題の性質上、大きな系統誤差を伴いやすいプロセスであった。本論文は、解析接続を行うことなく、ユークリッド格子のデータから直接、アクシオン生成率の重み付き積分を計算する課題に取り組んでいる。

手法
著者らは、ユークリッド相関関数を用いてスペクトル関数の重み付き積分を計算することを可能にする、最近のフォトンの生成率に関する研究（Harvey Meyerらによるもの）に着想を得た手法を提案している。

1. 理論的枠組み:

   • アクシオンの生成率は、スペクトル関数 $\sigma(\omega, \omega)$ に関連している。
   • 著者らは、虚空間運動量を含むユークリッド相関関数 $H_E(\omega_n)$ を定義し、これはカーネル積分を介してスペクトル関数に関連付けられる。
   • $H_E(\omega_n)$ の直接的な評価は、カーネルが距離に対して指数関数的に増大し、積分が紫外（UV）領域および高周波数において発散する可能性があるため、困難を伴う。
   • この発散問題を解決するために、著者らは松原モードの線形結合、具体的には $H_E^{(0,1)}$ および $H_E^{(0,1,2)}$ を構成する。これらの組み合わせは、主要なUV発散をキャンセルし、高周波数での収束を保証する減算カーネルを導入するものである。

2. 格子実装:

   • セットアップ: 本研究は、純粋 $SU(3)$ ヤン・ミルズ理論（クエンチ近似）において、格子サイズ $N_\tau = 16$（アスペクト比 4）、温度 $T = 1.5 T_c$ で行われる。
   • グラディエントフロー: トポロジカル電荷密度演算子 $q(x)$ の乗法的繰り込みと短距離ノイズを処理するために、著者らは改良された（Zeuthen）グラディエントフローを採用している。これにより、演算子は半径 $r \sim \sqrt{8t_f}$ でスミアリングされ、トポロジカル電荷が安定化される。
   • テイル・フィッティング: ユークリッド相関関数 $G(\omega_n, r)$ は大きな距離で指数関数的に減衰する。しかし、重み関数 $\sinh(\omega_n r)/(\omega_n r)$ は指数関数的に増大する。積分を計算するために、著者らは相関関数の「テイル」（大きな $r$）を、コンパクト化された幾何学における理論のスペクトル表現から導かれる指数関数の和（$e^{-mr}/r$）にフィットさせる。そして、精密な格子データ（小さな $r$）とフィットされた関数（大きな $r$）を組み合わせることで、積分を計算する。
   • 外挿: 結果は有限のフロー時間 $t_f$ において得られる。著者らは、$t_f \gtrsim 0.4a^2$ においてスミアリング関数は実質的に $O(4)$ 不変であり、統合された量 $H_E$ のフロー時間依存性は無視できる（$t_f^2$ 以上でスケールする）と主張している。彼らは、異なるフロー時間にわたるデータに対して定数フィットを行い、$t_f \to 0$ への外挿を行っている。

主な貢献

• 解析接続を用いない手法: 本論文は、不安定な解析接続を回避し、ユークリッド格子シミュレーションから直接、アクシオン生成率の重み付き積分を抽出するための実行可能な経路を示している。
• 指数関数的なノイズの処理: 著者らは、積分における指数関数的に増大する重み関数によって引き起こされる信号対雑音比（S/N比）の低下を克服するために、グラディエントフローとテイル・フィッティング技術を組み合わせた堅牢な数値戦略を開発した。
• 減算スキーム: 積分の有限性と格子上での定義可能性を確保する、特定の松原モードの線形結合（$H_E^{(0,1)}$ および $H_E^{(0,1,2)}$）の定式化。

結果

• 数値値: 純粋なグルーオンQCD（$T = 1.5 T_c$）において、著者らは以下の値を正常に抽出することに成功した：
  • $H_E^{(0,1)} = -0.0357(58) T^4$
  • $H_E^{(0,1,2)} = -0.0152(91) T^4$
    負の符号は、トポロジカル電荷密度演算子の反エルミート性に一致している。
• 摂動論との比較: 結果は、次次項（NLO）摂動計算と比較されている（$\alpha_s \approx 0.3$ を使用）。格子の結果はNLOの予測と一致しているが、最高次（LO）の結果は著しく小さすぎる。
• 実現可能性: 本研究は、単一の格子間隔において、適度な計算コスト（約60,000 CPUコア時間）で、少なくとも1つ、あるいはおそらく2つの有限な重み付き積分を抽出することが可能であることを確認した。

意義と主張
本論文は、他の手法が失敗する強結合のQCD領域において、熱的アクシオン生成率を決定するための新しい道を切り開くと主張している。純粋なグルーオン理論においてこれらの積分を計算することに成功したことにより、この手法の実現可能性を確立した。著者らは、この手法を、クロスオーバー付近およびその少し上の温度（$T \in [200, 1000]$ MeV）における完全なQCDへと拡張できると考えている。この温度域は、アクシオン暗黒物質に関する宇宙論的制約にとって極めて重要である。

著者らは、自らの手法が微分的な生成率 $d\Gamma/d^3k$（運動量 $k$ の関数）ではなく、積分のモーメントを与えるものであることを指摘し、即時的な現象論的影響については控えめな姿勢を保っている。また、これらの積分制約を、アクシオンの質量や結合定数の精密な限界へと翻訳するには、さらなる発展が必要であることも認めている。さらに、彼らは、この手法の成功は、最も軽いフレーバーレスの $0^{-+}$ 状態が温度スケールに対して十分に重く、相関関数が成長する重み関数よりも速く減衰するという条件に依存していることを強調している。この条件は、結合が弱くなる非常に高い温度においては、困難になる可能性がある。