The QCD crossover temperature and curvature coefficient from unbiased exponential resummation at physical quark masses in (2+1)-flavor lattice QCD

技術要約：(2+1)フレーバー格子QCDにおける無偏向指数関数的再総和（Unbiased Exponential Resummation）

問題提起
非零の温度およびバリオン化学ポテンシャル（$\mu_B$）における強相互作用物質の熱力学を決定することは、量子色力学（QCD）の中心的目標である。主要な課題は、$\mu_B = 0$ 付近におけるQCDクロスオーバー線とその曲率を特定することであり、これはQCD相図の制約および重イオン衝突実験の解釈において不可欠である。標準的な格子QCDのアプローチは、有限の$\mu_B$における符号問題によって阻害されている。その結果、ほとんどの第一原理研究は、$\mu_B = 0$ の周りでの熱力学的観測量のテイラー展開、あるいは虚数$\mu_B$からの解析接続に依存している。この「無偏向指数関数的再総和（unbiased exponential resummation）」の枠組みは、有限密度表現のQCD分配関数を再構成するために開発されているが、擬臨界観測量、具体的には擬臨界温度（$T_{pc}$）および主要な曲率係数（$\kappa_2^B$）の決定への適用については未探索のままである。

手法
著者は、物理的なクォーク質量を持つ(2+1)フレーバー格子QCDシミュレーションに対し、$32^3 \times 8$ の格子上で、Highly Improved Staggered Quark (HISQ) アクションを用いて、無偏向指数関数的再総和法を適用している。温度範囲は $T \in [135, 176]$ MeV である。

手法の核心は、$N$次まで打ち切られた、有限の$\hat{\mu}_B$における分配関数の比 $R_N$ の再構成にある：
$$R_N = \frac{Z_N(\hat{\mu}_B)}{Z(0)} = \left\langle \exp \left[ \sum_{n=1}^N \frac{D_n}{n!} \hat{\mu}_B^n \right] \right\rangle$$
ここで、$D_n$ は $\hat{\mu}_B=0$ におけるフェルミオン行列式の対数の微分である。著者は、$N=4$ を用いて、$\mathcal{O}(\hat{\mu}_B^2)$ まで正確にテイラー展開を再現しつつ、指数構造を通じて高次の項を生成する、$\mathcal{O}(\hat{\mu}_B^4)$ までのバリオン数感受率 $\hat{\chi}_2^B$ を再構成する。

$T_{pc}$ と $\kappa_2^B$ を決定するために、2つの異なるアプローチが採用されている：

1. 熱的微分解析（Thermal Derivative Analysis）: $\hat{\chi}_2^B$ の $\hat{\mu}_B \in [0, 0.5]$ における温度依存性を再構成する。2つの熱的微分観測量、$\hat{\chi}_T = T \partial \hat{\chi}_2^B / \partial T$ および $\hat{\chi}_C = T_c \partial \hat{\chi}_2^B / \partial T$ （ここで $T_c = 145$ MeV）を構成する。$\hat{\mu}_B = 0$ におけるこれらの微分のピーク位置が、擬臨界温度 $T_{pc}^T$ および $T_{pc}^C$ を定義する。
2. スケーリング関係と $\hat{\chi}_4^B$: 再構成された $\hat{\chi}_2^B$ の $\hat{\mu}_B$ 依存性から、4次のバリオン数感受率 $\hat{\chi}_4^B$ を抽出する。$\hat{\chi}_4^B$ のピークは、3番目の推定値 $T_{pc}^\chi$ を提供する。
3. 曲率の決定:
   • スケーリングに基づく手法: $\hat{\chi}_4^B$ と $\hat{\chi}_2^B$ の熱的微分の間のスケーリング関係を用いて、曲率を $\kappa_2^B = \frac{1}{6} \frac{\hat{\chi}_4^B}{\hat{\chi}_C}$ として推定する。この手法は、スケーリング領域における非普遍的なスケーリングパラメータの相殺に依拠している。
   • 直接フィット（Direct Fit）: 擬臨界線の $\hat{\mu}_B$ 依存性の明示的な解析として、ピーク位置（$T_{pc}^T$ および $T_{pc}^C$）を $\hat{\mu}_B \in [0, 0.3]$ の範囲で $T_{pc}^x(\hat{\mu}_B) = T_{pc}^x(0) [1 - \kappa_{2,x}^B \hat{\mu}_B^2 + \mathcal{O}(\hat{\mu}_B^4)]$ の形式にフィットさせる。

主な結果

• 擬臨界温度: $\hat{\mu}_B = 0$ において、熱的微分解析は以下を与える：
  • $T_{pc}^T = 160.2(4)(3)$ MeV
  • $T_{pc}^C = 158.1(5)(3)$ MeV
    $\hat{\chi}_4^B$ のピークは、$T_{pc}^\chi = 159.4(10)(12)$ MeV（フィット法）および $159.4(9)(8)$ MeV（固定法）を与える。これらの値は互いに一致しており、従来のテイラー展開の結果（例：HotQCDコラボレーション）とも一致している。
• 曲率係数 ($\kappa_2^B$):
  • スケーリング関係: スケーリングに基づく推定値は、異なる $T_{pc}$ の定義間で安定している：
    • $\kappa_2^B(T_{pc}^T) = 0.02082(14)$
    • $\kappa_2^B(T_{pc}^C) = 0.02064(15)$
    • $\kappa_2^B(T_{pc}^\chi) = 0.02075(14)$
  • 直接フィット: $\hat{\mu}_B$ 依存性の直接解析は、より大きな中心値を示すが、有意に大きな不確かさを伴う：$\kappa_{2,T}^B = 0.051(35)$ および $\kappa_{2,C}^B = 0.051(34)$。
• 整合性: 直接フィットによる不確かさが大きいにもかかわらず、結果はスケーリングに基づく推定値と統計的に互換性がある。スケーリングに基づく曲率の値は、同等の温度におけるテイラー展開法から得られた値とも一致している。

意義および主張
本論文は、物理的なクォーク質量におけるQCD擬臨界温度および曲率の決定に対する、無偏向指数関数的再総和の初適用を提示していると主張している。主な意義は、この再総和の枠組みが、期待されるスケーリング挙動および従来のテイラー展開に基づく格子QCDの結果と整合していることを実証することにある。

具体的には、著者は以下のことを確立している：

1. 無偏向指数関数的再総和は、すべての高次テイラー係数を明示的に計算することなく、有限密度の感受率およびその熱的微分を正常に再構成できる。
2. この手法から導出された観測量は、曲率係数に必要なスケーリング関係を満たしており、再総和された記述に対する非自明なクロスチェックを提供している。
3. 本手法は、小さな有限バリオン密度におけるQCDクロスオーバーを調査するための、テイラー展開に対する補完的なアプローチとして機能する。

著者は、直接的な擬臨界線による曲率決定の精度については、現在の結果がスケーリングに基づく推定値よりも「大幅に制約が弱い」と述べ、謙虚な姿勢を保っている。彼らは、この枠組みが実行可能なツールとして確立された一方で、有限の $\hat{\mu}_B$ における擬臨界温度の定義の改善、統計量の増加、および連続体外挿を行うためのさらなる作業が必要であると結論づけている。