Non-perturbative determination of the QCD Equation of State up to the electroweak scale

この論文は、シフト境界条件と非摂動的に定義された結合定数のランニングに基づく定物理線決定法を組み合わせることで、3 GeV から電弱スケールまでの広範な温度領域において、3 質量ゼロクォークを持つ QCD の状態方程式を非摂動的に決定し、摂動論だけでは記述できない非摂動効果の重要性を明らかにしたものである。

要約

1. 何をしたのか？「宇宙のレシピ」を完成させた

この研究の目的は、「クォークとグルーオン（物質の最小単位）」が、超高温でどう振る舞うかを調べる「状態方程式（EoS）」というレシピを完成させることです。

• 対象： 温度が**30 億度（3 GeV）から 1 兆度（165 GeV）**という、とてつもなく高温の領域。
  • これは、太陽の中心の温度（約 1500 万度）の何万倍もの熱さです。
  • 宇宙が誕生してから、電弱相互作用（素粒子の力が分かれる瞬間）が起きるまでの「超高温の宇宙」の状態を再現しています。
• 成果： これまでの計算では「高温になれば、理論（摂動論）で簡単に計算できるはず」と思われていましたが、実は**「そう単純ではない」ことが分かりました。この研究は、その複雑な振る舞いを、理論だけでなく、「第一原理（基本法則）から直接計算する」**ことで、0.5%〜1% という驚異的な精度で描き出しました。

2. 使った魔法の道具：2 つの工夫

この超高温のシミュレーションを成功させるために、研究者たちは 2 つの「魔法の道具」を組み合わせて使いました。

① 「動く窓」から景色を見る（シフトされた境界条件）

通常、箱の中で熱いガス（クォーク・グルーオン）をシミュレーションする場合、箱の壁（境界）で粒子が跳ね返るのを計算する必要があります。しかし、高温になるとこの計算が非常に重く、誤差が出やすくなります。

• アナロジー： 窓の外を走る列車の景色を撮影するとします。
  • 普通の方法： 窓を固定して、外の景色が流れていくのを待つ（静止した箱）。
  • この研究の方法： カメラ自体を動かす（動く窓）。
  • 研究者は、箱の「時間方向」の壁を少しずらす（シフトする）ことで、**「箱ごと移動している」**ような状態を作りました。
  • これにより、**「エントロピー（無秩序さの量）」という重要な値を、余計な計算（引き算）を一切せず、「直接」**測ることができるようになりました。まるで、動く車の中で直接風速を測るようなものです。

② 「段々階段」で高所へ登る（一定物理の線）

高温になるほど、計算に必要な「格子（コンピューターのメモリ）」を細かくし、かつ箱を大きくする必要があります。これは、**「小さな砂粒で巨大な山を正確に描く」**ようなもので、計算量が爆発的に増えます。

• アナロジー： 高い山（電弱スケール）に登るには、麓（ハドロンスケール）から一歩ずつ登る必要があります。
  • 従来の方法だと、麓の「石の大きさ」を基準にすると、山頂では石が小さすぎて描けなくなります。
  • この研究では、「山の標高（エネルギー）」に合わせて、登るペース（格子の細かさ）を自動的に調整する「一定物理の線」という地図を使いました。
  • これにより、低い温度から高い温度まで、**「つなぎ目のない滑らかな道」**で、計算の精度を保ちながら登りきることができました。

3. 発見した驚きの事実：理論は「不完全」だった

これまで物理学者たちは、「温度が非常に高くなれば、物質は理想気体のように振る舞い、既存の理論式で正確に予測できる」と信じていました。

• 予想： 高温になるほど、理論と実験（シミュレーション）の結果はピタリと合うはず。
• 現実： 合いませんでした。
  • 電弱スケール（1 兆度）に近づいても、既存の理論式だけでは説明できない「見えない要素」がまだ残っていました。
  • アナロジー： 天気予報で「明日は晴れ」と言われていたのに、実際には「予報にはない突然の雷雨」が降ったようなものです。
  • この研究は、**「既存の理論式（摂動論）には、高温でも無視できない『非摂動的な（計算しにくい）要素』が潜んでいる」**ことを証明しました。

4. なぜこれが重要なのか？

この研究は、単に数値を正確にしただけではありません。

1. 宇宙の歴史の解明： 宇宙の初期、重力波（時空のさざ波）がどう生成されたかを理解するために、この「高温の物質の性質」が不可欠です。この研究は、そのピースを正確に埋めました。
2. 新しい計算手法の確立： 「動く窓」と「段々階段」の組み合わせは、今後、4 つや 5 つのクォークを扱う場合（より現実的な宇宙のモデル）にもそのまま使えます。
3. 理論の限界の明示： 「どんなに高温でも、単純な理論式では説明できない」という事実を突き止め、物理学の新しい探求の道を開きました。

まとめ

この論文は、**「超高温の宇宙という『見えない世界』を、新しい計算テクニックを使って、初めて鮮明な写真として撮り上げた」**と言えます。

それまで「理論で推測できる範囲」と思われていた領域に、**「理論では説明できない複雑なドラマ」**が隠れていることを発見し、その詳細を 99% 以上の精度で描き出した、画期的な成果なのです。

技術概要

論文要約：電弱スケールまでの QCD 状態方程式の非摂動的決定

論文タイトル: Non-perturbative determination of the QCD Equation of State up to the electroweak scale
著者: Michele Pepe (INFN, Milano-Bicocca)
発表: LATTICE2025 (2025 年 11 月)

1. 背景と課題 (Problem)

量子色力学（QCD）の状態方程式（EoS）は、クォークとグルーオンの微視的ダイナミクスと、宇宙論的な巨視的進化を結びつける重要な物理量です。特に、初期宇宙の進化や重イオン衝突実験におけるクォーク・グルーオンプラズマ（QGP）の理解に不可欠です。

従来の研究における課題は以下の通りです：

• 温度範囲の限界: これまでの格子 QCD による第一原理からの EoS 決定は、ハドロンスケール（$T \sim 150$ MeV）から約 1 GeV 程度までに限られていました。
• 摂動論の限界: 1 GeV 以上の高温領域では、強結合定数 $g$ に関する摂動展開に依存してきました。しかし、摂動級数の収束性は非常に悪く、圧力については $O(g^6)$ のオーダーで、本質的に非摂動的な「ウルトラソフトモード」の影響を受け始めます。
• 電弱スケールまでの未解明: 電弱スケール（$T \sim 100$ GeV）に近づいても、摂動論の予測は完全な結果から大きく乖離している可能性が示唆されており、ステファン・ボルツマン限界への滑らかなアプローチの背後に、非摂動的な寄与が隠れている可能性があります。

本研究は、3 GeV から 165 GeV（電弱スケール）までの広範な温度領域において、$N_f=3$ の質量ゼロクォークを含む QCD の状態方程式を、完全に非摂動的に決定することを目的としています。

2. 手法と理論的枠組み (Methodology)

本研究は、従来の手法の限界を克服するために、以下の 2 つの革新的な戦略を組み合わせることで実現されました。

2.1 物理定数線（LCPs）の決定とスケール設定

高温領域では、ハドロン観測量を用いた従来のスケール設定は計算コストが膨大になり実用的ではありません。そこで、以下のアプローチを採用しました：

• シュレーディンガー汎関数（SF）結合: 有限体積で定義された非摂動的な結合定数 $\bar{g}^2_{SF}$ を用います。
• ステップスケーリング: 格子間隔 $a$ を変化させながら、結合定数のランニングを追跡することで、ハドロンスケールから電弱スケールまで連続的に物理定数線（LCPs）を定義します。これにより、高温シミュレーションに必要なバースパラメータを高精度で決定できます。

2.2 移動座標系とシフト境界条件

従来の積分法（Trace Anomaly の積分）では、紫外発散の除去に異なる温度スケール間での引き算が必要となり、数値的に困難を伴います。これを回避するため、**シフト境界条件（Shifted Boundary Conditions）**を導入しました：

• 移動座標系: ユークリッド時間方向に空間的なシフト $\xi$ を加えることで、理論を移動する座標系として定式化します。
• エントロピー密度の直接測定: この幾何学的設定により、エントロピー密度 $s(T)$ を温度の引き算なしに直接導出できます。
  • 関係式: $s/T^3 = \frac{1+\xi^2}{\xi_k T^4} \frac{\partial f_\xi}{\partial \xi_k}$
  • ここに $f_\xi$ は移動座標系での自由エネルギー密度です。
• この手法は、真空の寄与を大幅に削減し、数値計算の安定性を向上させます。

2.3 数値計算の具体的実施

エントロピー密度の計算は、以下の 2 つの寄与に分解して行われました：

1. 純ゲージ部分（無限重クォーク極限）: 結合定数 $g_0^2$ に関する積分を行い、SU(3) 純ゲージ理論の寄与を評価。
2. ダイナミッククォーク部分: 質量ゼロの 3 種クォークの寄与を、クォーク質量 $m_q$ に関する積分（スカラー密度 $\langle \bar{\psi}\psi \rangle$ の期待値の積分）を通じて評価。

計算は、$L_0/a = 4, 6, 8, 10$ の 4 つの時間方向格子解像度と、$L/a = 144$ の空間サイズで行われ、有限体積効果を抑制しました。

3. 主要な結果 (Key Results)

• 高精度な EoS の決定: 3 GeV から 165 GeV の温度範囲で、エントロピー密度 $s/T^3$ を 0.5%〜1.0% の相対誤差で決定しました。
• 摂動論との比較:
  • 決定された非摂動的なデータは、既知の摂動論展開（$O(g^6)$ まで）や、硬熱ループ（HTL）再総和法による予測とは一致しません。
  • 特に、$O(g^6)$ の項以降に現れる非摂動的な寄与を考慮しないと、電弱スケール付近のデータは説明できません。
  • 摂動論の収束性は、$T \sim 100$ GeV であっても十分ではなく、高次項や非摂動的効果が依然として重要です。
• パラメトリゼーション:
  • 非摂動的なデータを記述するために、既知の摂動論係数に $O(g^6)$ および $O(g^7)$ の項（非摂動的寄与を含む）を追加した多項式展開を用いてフィッティングを行いました。
  • このフィッティングは、500 MeV 付近の既存の格子データ（$N_f=2+1$）とも整合性があり、$T \ge 500$ MeV 全体を記述する有効なパラメトリゼーションとなりました。
• 熱力学ポテンシャル: エントロピー密度から一貫して圧力 $p(T)$ とエネルギー密度 $e(T)$ を導出し、熱力学的一貫性を検証しました。

4. 意義と結論 (Significance and Conclusions)

• 画期的な到達点: 本研究は、QCD の熱力学をハドロンスケールから電弱スケールまで、第一原理（格子 QCD）に基づいて非摂動的に決定した世界初の試みです。
• 摂動論の限界の明確化: 電弱スケールに近づいても QCD は単純な摂動的な理想気体に近づいておらず、非摂動的な構造が重要であることを実証しました。これは、初期宇宙の進化や重力波スペクトルへの影響を評価する上で重要な知見です。
• 汎用性の高い手法: 提案された「シフト境界条件」と「シュレーディンガー汎関数結合を用いたスケール設定」の組み合わせは、$N_f=4$ や $N_f=5$ の質量を持つクォークを含む標準模型の完全な熱力学進化の記述にも適用可能な汎用的な枠組みを提供します。
• 将来展望: 現在、エネルギー・運動量テンソルの非摂動的な重整化定数の決定が進められており、これによりクォーク・グルーオンプラズマの輸送係数（粘性など）の高精度決定が可能になり、重イオン衝突および宇宙論的文脈での理解がさらに深まることが期待されます。

総じて、この研究は高温 QCD の理解に新たな基準を設け、摂動論の限界を超えた非摂動的な QCD の性質を明らかにした重要な成果です。