Kinetic Freeze-Out Conditions and Net Baryon Density in Au+Au Collisions at $\sqrt{s_{NN}} = 7.7$--$39$ GeV within a Collective Flow Fireball Model

RHIC ビームエネルギー・スキャンにおける Au+Au 衝突の STAR データを用いた統計的ファイアボールモデルの解析により、縦方向の集団的流れを考慮することで凍結温度が上昇し、特に大きな縦流速では格子 QCD が予測するクォーク・グルーオンプラズマからハドロン相への転移温度を超えることが示され、このエネルギー領域では大きな縦流速が物理的に不利である可能性が指摘されるとともに、ネットバリオン密度の最大値が$\sqrt{s_{NN}} \lesssim 11.5$ GeV 付近に存在し、動的流れがバリオン圧縮を最大約 20% 増大させることが明らかになった。

要約

1. 実験の舞台：「超高温の火の玉」

まず、実験の状況をイメージしてください。
金（Au）の原子核同士を、光速に近い速さで正面衝突させます。すると、一瞬にして**「クォーク・グルーオンプラズマ（QGP）」**という、原子核の部品がバラバラになった超高温・高密度の「火の玉」が生まれます。

この火の玉は、すぐに膨張して冷えていき、やがて「ハドロン（陽子や中性子など）」という粒子の塊に戻ります。この**「火の玉が冷えて、粒子がもうそれ以上相互作用しなくなる瞬間」を、物理学者は「運動学的凍結（Kinetic Freeze-out）」**と呼びます。

この研究は、**「その凍結した瞬間、火の玉はどれくらい熱かったのか（温度 T）」と「どれくらい粒子が詰まっていたのか（バリオンの密度）」**を、より正確に測ろうとしたものです。

2. 従来の考え方の限界：「静止している風船」

これまでの研究では、この火の玉を**「静止しているお風呂の湯」**のように扱ってきました。
「湯の温度が一定で、粒子がその中をランダムに動いている」と仮定して、実験データ（飛び散る粒子の速さの分布）に当てはめて、温度を計算していました。

しかし、現実はもっとダイナミックです。
火の玉は**「風船が破裂して勢いよく膨張している」**ような状態です。

• 横方向への膨張（横流）： 風船が横に広がる動き。
• 縦方向への膨張（縦流）： 風船が前後に伸びる動き。

この論文の著者たちは、**「この『膨張する動き（流れ）』を無視してはいけない」と考えました。特に、「前後方向（縦方向）への動き」**が、温度の計算にどう影響するかを初めて詳しく調べました。

3. 発見：「動き」と「温度」のトリック

ここで、この研究の最も面白い発見（トリック）があります。

**「風船が前後に速く動いていると、観測者には『もっと熱い』ように見える」**という現象です。

• 例え話：
  あなたが静止しているお風呂（温度 40 度）に入っているとします。
  しかし、もしそのお風呂自体が、前後に猛烈な勢いで揺れ動いていたらどうでしょう？
  粒子（湯気）があなたにぶつかる瞬間、その「揺れの勢い」が粒子の速さに加算されます。
  観測者は「粒子が速い＝温度が高い」と勘違いしてしまいます。

この論文では、**「縦方向の動き（縦流）」と「温度」が、計算上は「ごまかし合い（縮退）」**を起こしていることを発見しました。

• 縦方向の動きを無視すると、温度は低めに出る。
• 縦方向の動きを考慮して計算し直すと、同じデータから導き出される温度は、実際よりも高く（143〜209 メガ電子ボルトの範囲で）出るようになります。

4. 重要な結論：「どれくらいの動きが現実的か？」

この「温度のトリック」から、重要な結論が導き出されました。

• QCD の転移温度（Tc）という壁：
  理論物理学（格子 QCD）によると、物質が「粒子の塊（ハドロン）」から「バラバラの部品（QGP）」に変わる境界線の温度は、約155〜160 メガ電子ボルトです。
  これを超えると、もう「ハドロンガス」としてのモデルは使えません。

• 研究の結果：

  • 縦流が小さい場合（動きが穏やか）： 計算された温度は 150〜180 メガ電子ボルト程度。これは「境界線」の近くか、少し上ですが、まだハドロンガスとして説明がつけられます。
  • 縦流が大きい場合（激しく前後に動く）： 計算された温度は200 メガ電子ボルトを超えてしまいます。これは「境界線」を大きく超えており、**「もうハドロンガスではない（QGP になっているはず）」**ことを意味します。

**つまり、「縦方向に激しく動く（vz=0.4）」という仮定は、このエネルギー領域では物理的にありえない（不自然だ）」**という結論に至りました。
「風船が前後に飛び跳ねるほど激しく動くなら、それはもう『ハドロンガス』のモデルでは説明できない」というわけです。

5. 物質の密度：「圧縮されたスポンジ」

もう一つの発見は、**「バリオンの密度（物質の詰まり具合）」**についてです。

• 衝突エネルギーが低いと、原子核同士がぶつかり合って止まりやすいため、中央部分に物質がギュッと詰まります。
• 衝突エネルギーが高いと、すり抜けてしまうため、物質は薄くなります。

この研究では、「集団的な流れ（膨張）を考慮すると、物質の圧縮は、静止している場合よりも最大で 20% 強く見積もられる」ことがわかりました。
まるで、「静止しているスポンジ」よりも、「勢いよく膨らんでいるスポンジ」の方が、一時的に内部の圧力が強く感じられるようなものです。

まとめ：この研究がなぜ重要なのか？

1. 「温度」の測り直し： 過去の計算では見落としていた「縦方向の動き」を考慮することで、凍結時の温度が実際よりも高い可能性があると示しました。
2. 物理的な制約： 「動きが激しすぎると、温度が物理的にありえない値（QGP 領域）になってしまう」ため、**「このエネルギー領域では、縦方向の動きはそれほど激しくない」**という制約を見つけました。
3. 将来への道しるべ： この結果は、将来のより高度なシミュレーション（流体の動きを詳しく計算するモデル）や、FAIR（ドイツ）や RHIC（アメリカ）の新しい実験にとって、**「正しい基準（ベンチマーク）」**となります。

一言で言うと：
「爆発する火の玉の『前後の揺れ』を無視すると、温度の計算がズレてしまう。そして、激しく揺れすぎると、そのモデル自体が破綻してしまうことがわかった。だから、火の玉はほどほどに揺れているはずだ」という、**「火の玉の動きと温度の正しい関係」**を突き止めた研究です。

技術概要

以下は、提供された論文「Kinetic Freeze-Out Conditions and Net Baryon Density in Au+Au Collisions at √sNN = 7.7–39 GeV within a Collective Flow Fireball Model」の技術的な要約です。

論文概要

タイトル: 集団流ファイアボールモデルにおける Au+Au 衝突（√sNN = 7.7–39 GeV）の運動学的凍結条件と正味のバリオン密度
著者: Sk Noor Alam, Victor Roy
対象: RHIC ビームエネルギー・スキャン（BES）プログラムにおける 0–5% 中心 Au+Au 衝突

---

1. 研究の背景と課題 (Problem)

重イオン衝突実験の主要な目的の一つは、量子色力学（QCD）の相図を温度（$T$）と化学ポテンシャル（$\mu_B$）の空間でマッピングし、臨界点の位置を特定することです。

• 既存の課題: 従来の熱的モデル（ファイアボールモデル）では、ファイアボールを静的な熱源として扱うことが多く、集団流（Collective Flow）の影響を十分に考慮していない場合がありました。
• 具体的な問題点: 粒子の横運動量（$p_T$）スペクトルは、熱的運動と集団流（横方向および縦方向）の両方の影響を受けます。特に、縦方向の集団流（$v_z$）を無視したり、適切に扱わなかったりすると、凍結温度（$T$）や化学ポテンシャル（$\mu_B$）などの凍結条件を抽出する際にバイアスが生じる可能性があります。
• 本研究の目的: 横方向および縦方向の集団流を明示的に組み込んだ共変統計ファイアボールモデルを用いて、RHIC の BES エネルギー領域（7.7–39 GeV）における運動学的凍結条件と正味のバリオン密度（$\rho_B$）を再評価し、流体力学的な膨張が凍結軌道に与える影響を明らかにすること。

2. 手法 (Methodology)

• モデル: 共変統計ファイアボールモデル（Covariant Statistical Fireball Model）を採用。
  • 静止ファイアボールモデル（Randrup & Cleymans の先行研究）を拡張し、局所静止系での熱的分布に、横方向流速 $v_T$ と縦方向流速 $v_z$ を持つ流体要素の集団流を Lorentz 共変的に組み込みました。
  • 粒子の分布関数は、流体の 4 元速度 $u^\mu$ と粒子の 4 元運動量 $p^\mu$ の内積 $p_\mu u^\mu$ に依存する Lorentz 不変な形式で記述されます。
• データ: STAR コラボレーションによる 0–5% 中心 Au+Au 衝突の陽子と正パイオン（$\pi^+$）の横運動量スペクトル（$|y| < 0.1$）を使用。
• フィッティング手順:
  • 縦方向流速 $v_z$ を 3 つの固定値（0.0, 0.2, 0.4）で設定。
  • 各 $v_z$ に対して、横運動量スペクトルに式 (17) を最小二乗法でフィットし、以下のパラメータを同時に抽出：
    • 凍結温度 $T$
    • バリオン化学ポテンシャル $\mu_B$
    • 横方向流速 $v_T$
    • 規格化定数 $\lambda$
• 熱力学的量の計算: 抽出された $T$ と $\mu_B$ を用いて、ハドロン共鳴ガス（HRG）モデルに基づき、正味のバリオン密度 $\rho_B$ とエネルギー密度 $\varepsilon$ を計算し、$(\rho_B, \varepsilon)$ 平面での凍結軌道を再構築しました。

3. 主要な結果 (Key Results)

A. 凍結温度 ($T$) と縦方向流速 ($v_z$) の関係

• 系統的な温度上昇: 縦方向流速 $v_z$ を増加させると、抽出される凍結温度 $T$ が系統的に上昇しました。
  • $v_z = 0.0$: 143–171 MeV
  • $v_z = 0.2$: 150–180 MeV
  • $v_z = 0.4$: 175–209 MeV
• 物理的メカニズム: この温度上昇は、$p_T$ スペクトルの硬化（Hardening）によるものではなく、**運動学的な縮退（Kinematic Degeneracy）**に起因します。
  • 分布関数において、流体の縦方向運動が Lorentz 因子を通じて粒子の実効エネルギーに寄与します。
  • 観測されたスペクトルを再現するために、フィッティングアルゴリズムは縦方向流速の効果を補償するために、より高い温度 $T$ を選択せざるを得なくなります。
• QCD クロスオーバー温度との比較:
  • $v_z = 0.0$ および $0.2$の場合、抽出された温度は格子 QCD による QCD クロスオーバー温度（$T_c \approx 155–160$ MeV）と概ね一致します。
  • しかし、$v_z = 0.4$ の場合、特に高エネルギー側で $T$ が $T_c$ を大幅に超えます（最大 209 MeV）。これは、ハドロン共鳴ガスモデルが適用不可能な領域（クォーク・グルーオンプラズマ領域）を示唆しており、このエネルギー領域において $v_z = 0.4$ という大きな縦方向流速は物理的に好ましくない可能性を示しています。

B. バリオン化学ポテンシャル ($\mu_B$) と正味のバリオン密度 ($\rho_B$)

• $\mu_B$ の振る舞い: $\mu_B$ は $v_z$ に依存せず、衝突エネルギーの増加とともに単調に減少しました（7.7 GeV で約 420 MeV から 39 GeV で約 200 MeV へ）。これは、エネルギーが高いほど核子の停止（Baryon Stopping）が小さくなり、中間ラピディティでの正味のバリオン密度が低下することを反映しています。
• $\rho_B$ の最大値: $(\rho_B, \varepsilon)$ 平面における凍結軌道を再構築した結果、中間エネルギー（$\sqrt{s_{NN}} \lesssim 11.5$ GeV）で正味のバリオン密度に明確な最大値が存在することが確認されました。
• 集団流の影響: 動的な流（Collective Flow）を考慮することで、静的な限界（$v_z=0$）と比較して、推定されるバリオン圧縮（$\rho_B$）が最大で約 20% 増大することが示されました。

C. 粒子収量

• 陽子収量: 中間ラピディティ（$|y| < 0.1$）における陽子収量は、$v_z$ の変化に対してほぼ不変でした。これは、正味のバリオン量が衝突の初期段階（平衡前）のバリオン停止によって決定され、その後の集団的縦方向流によるラピディティのシフトが狭い窓内での積分収量にほとんど影響を与えないためです。
• パイオン収量: 衝突エネルギーの増加とともに単調に増加し、高エネルギー側でのメソン支配への移行を示しています。

4. 貢献と意義 (Significance)

• 凍結条件の精緻化: 集団流（特に縦方向流）を明示的に扱うことで、凍結温度の抽出における系統的なバイアスを明らかにしました。$v_z$ を無視して $T$ を決定すると、物理的に不適切な高温領域に陥るリスクがあることを示唆しています。
• 物理的に妥当なパラメータ範囲の特定: 抽出された温度が QCD クロスオーバー温度 $T_c$ を超えないという制約を課すことで、$v_z = 0.0$ および $0.2$のようなより物理的に妥当な流速範囲が特定でき、$v_z = 0.4$ のような過大評価が排除されました。
• 将来の研究への指針:
  • FAIR や低エネルギー RHIC 実験における衝突エネルギーの選択や、臨界点探索のための実験計画に対して、より信頼性の高い凍結ベンチマークを提供しました。
  • 今後の研究では、ストレンジハドロン（K, $\Lambda$, $\Xi$など）を同時にフィットして化学ポテンシャル（$\mu_S, \mu_Q$）を決定することや、より広いラピディティ範囲でのデータを用いて $v_z$ の縮退を解くこと、そして粘性流体力学モデルとの結合が提案されています。

結論

本研究は、RHIC BES エネルギー領域における Au+Au 衝突の運動学的凍結条件を、横・縦両方の集団流を考慮した共変統計モデルを用いて再評価しました。その結果、縦方向流速と温度の間の運動学的縮退が凍結温度の抽出に大きな影響を与えること、そして物理的に妥当なモデル（HRG）の適用範囲（$T \lesssim T_c$）を考慮すると、大きな縦方向流速はこれらのエネルギーでは支持されないことが示されました。また、集団流を考慮することで、中間エネルギー領域における正味のバリオン密度の最大値がより顕著に増大することが確認され、QCD 相図の探索における重要な知見を提供しました。