Probing the onset of hydrodynamization in peripheral p-Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV

要約

あなたが、まだ流体として振る舞うことができる最小の水たまりの大きさを突き止めようとしていると想像してください。巨大な海があれば、それは容易に流れます。しかし、一滴だけあれば、そこにとどまったり、バラバラになったりするかもしれません。では、その境界線はどこにあるのでしょうか。水分子の集まりが、流体として振る舞うのをやめ、個々の混沌とした粒子として振る舞い始めるのは、どの大きさでしょうか。

この論文は、クォーク・グルーオンプラズマ（QGP） におけるその正確な「転換点」を見つけることを目的としています。

QGP とは何か

QGP を宇宙の「原始のスープ」と考えてください。これはビッグバン直後のわずかな瞬間に存在した物質の状態です。この状態では、原子の構成要素（クォークとグルーオン）が溶け合い、超高温・超高密度の液体のように自由に流れています。

通常、科学者たちはこのスープを作るために、鉛のような重い原子を光速に近い速度で衝突させます。しかし最近、科学者たちは奇妙なことに気づきました。鉛のような重い原子同士を衝突させるのではなく、単一の陽子が鉛原子核に衝突する（p-Pb 衝突）といった、はるかに小さなものを衝突させた際にも、この「液体スープ」の兆候が現れるのです。

大きな疑問は、それは実際に液体なのか、それとも単に無秩序に跳ね回っている粒子の集まりに過ぎないのか、という点です。

実験：陽子を鉛に衝突させる

この論文の著者たちは、この「スープ」を流体力学（流れる液体を記述するための数学）の法則で記述できる最小の大きさを突き止めようとしていました。

彼らはJETSCAPEと呼ばれる大規模なコンピュータシミュレーションを使用しました。このシミュレーションは、衝突プロセス全体を 4 つの段階で再現するハイテクなビデオゲームエンジンと考えることができます。

1. 準備（TRENTo）：舞台を設定し、陽子と鉛原子核を初期位置に配置します。
2. 試合前（Freestreaming）：「液体」が形成される前に、粒子がごく短時間、自由に飛び回ります。
3. 流れ（MUSIC）：これが流体力学の部分です。シミュレーションは粒子を流れる流体として扱おうとします。
4. 結末（iSS + SMASH）：スープが冷えるにつれて、粒子は検出器で観測可能な実際の陽子、パイオン、その他の粒子へと凍結します。

テスト：スープはどの程度「液体」か

スープが実際に流体のように振る舞っているかどうかをテストするために、科学者たちは楕円流と呼ばれるものを観察しました。

アナロジー：2 台の車が正面衝突すると想像してください。もしそれらが完全に丸く、真ん中で衝突すれば、破片は円形に飛び散ります。しかし、少しずれて（すれ違い様に）衝突すれば、破片はより楕円形（フットボールのような形）に飛び散ります。

• もし内部の物質が完全な流体として振る舞えば、その楕円形に強く絞り出されます。
• もし物質が単に無秩序に跳ね回っている粒子であれば、楕円形は弱いか、あるいは存在しないでしょう。

科学者たちは、「周縁」衝突（陽子と鉛原子核の重なりが小さいすれ違い様の衝突）のシミュレーションを実行しました。彼らは問いかけました：流体としての振る舞いが崩壊する前に、この重なりはどの程度小さくできるか？

意外な展開：「緩和時間」のノブ

実際の流体では、流体を押してから反応するまでの間に遅れがあります。物理学ではこれをせん断緩和時間と呼びます。

著者たちはある手を使いました。この「緩和時間」のノブを極端な設定に回したのです。

• 彼らは問いかけました：「もし流体の反応が非常に鈍ければどうなるか？非常に速ければどうなるか？」
• 彼らはこれらの極端な条件下で楕円流（楕円形）を観察しました。

発見：転換点

彼らは、より一層「すれ違い様」な衝突（つまり、関与する物質の量、すなわちdN/dyが小さくなる衝突）をシミュレーションするにつれて、流体の振る舞いを観察しました。

• 結果：物質の量が単位ラピディティあたり約**7 粒子（dN/dy ≈ 7）**まで減少すると、流体の振る舞いは突然揺らぎ始め、崩壊し始めました。
• 比喩：人々が流体のように移動しようとしている群衆を想像してください。100 人いれば、スムーズに流れます。10 人でも、まだ流れるかもしれません。しかし、7 人まで減ると、彼らは個々にぶつかり始め、滑らかな「流れ」は消えてしまいます。

この論文は、彼らが研究したエネルギーにおける陽子 - 鉛衝突について、系が約 7 粒子より小さくなると流体力学は機能しなくなると結論付けています。それ以下では、「スープ」は液体として振る舞うには小さすぎ、単に個々の粒子の集まりに過ぎません。

なぜこれが重要なのか

これは、科学者が自然の根本的な限界を理解する助けとなります。それは、物質の「液体」状態が魔法ではなく、最小の大きさの要件を持っていることを教えてくれます。系が小さすぎれば、流体力学の法則はもはや適用されず、代わりに個々の粒子を見る必要があります。

著者たちはまた、彼らの結果が、以前の大規模な衝突（鉛 - 鉛衝突など）に関する研究とはわずかに異なっていたと指摘しました。おそらく、今回使用されたコンピュータモデルの方が安定しており、「試合前」の段階を異なった方法で処理していたためでしょう。

要約すると：彼らは、まだ「流体」と呼ぶことができる最小のクォーク・グルーオンプラズマの水たまりを見つけました。そして、その水たまりがまとまりを保つためには、少なくとも約 7 個の粒子を含んでいる必要があることがわかりました。

技術概要

技術的概要：周辺 p-Pb 衝突における流体力学的化の開始の探求

問題提起
相対論的流体力学は、大型の重イオン衝突（Pb-Pb、Au-Au）におけるクォーク・グルーオンプラズマ（QGP）のモデル化に成功しているが、小型衝突系（高多重度プロトン - プロトンおよびプロトン - 原子核）への適用性は激しい議論の対象となっている。中心的な問いは、これらの小型系で観測される集団的流れのシグネチャーが、真の QGP 様の熱い原子核物質効果に起因するのか、それとも冷たい原子核物質効果に起因するのかである。具体的には、著者らは流体力学的記述が破綻する最小の系サイズ（粒子多重度 $dN/dy$ で定量化される）を決定することを目的としている。著者らの先行研究では、周辺 Pb-Pb および Au-Au 衝突において $dN/dy \approx 10$ で流体力学的化の開始が確立されたが、本研究はこれを $\sqrt{s_{NN}} = 5.02$ TeV のより小さな p-Pb 系へと拡張したものである。

手法
本研究は、高エネルギー衝突のソフトセクターとハードセクターの両方を一貫してモデル化するように設計されたモジュール型シミュレーションツールであるJETSCAPEフレームワーク（バージョン 3.6.4）を利用している。ソフトセクターのシミュレーションチェーンは以下の通りである：

1. 初期条件： TRENToモデルが、核厚さ関数の一般化された平均に基づき、核子揺らぎを取り入れて初期エネルギー密度を生成する。
2. 非平衡期： 衝突なしのボルツマン方程式に基づくFreestreamingモジュールが、初期状態と流体力学フェーズの間のギャップを埋める。
3. 流体力学： MUSICソルバーが、2 次粘性流体力学（BRSSS 形式）を実装する。非流力学的モードを調節する 2 次輸送係数であるせん断緩和時間（$\tau_\pi$）が主要な変数である。
4. ハドロンアフターバーナー： iSS（粒子化）およびSMASH（ハドロン輸送）モジュールが、凍結とそれに続くハドロン再散乱を処理する。

流体力学的化の開始を探るため、著者らはせん断緩和時間（$\tau_\pi$）を極端な値まで変化させる。$\tau_\pi$ は $\tau_\pi(T) = b_\pi \frac{\eta}{s(T)} \frac{1}{T}$ によってせん断粘性とエントロピー密度の比（$\eta/s$）に関連しているため、2 つの異なる変化戦略を採用する：

1. 温度依存の $\eta/s$ を固定したまま、正規化定数 $b_\pi$ を 1 から 20 まで変化させる。
2. $b_\pi$ を固定したまま、温度依存の $\eta/s$ の範囲（下限の KSS 限界と上限のベイズ限界を使用）を変化させる。

流体力学的振る舞いの破綻は、楕円流（$v_2$）とその揺らぎを監視することで特定される。理論的な議論によれば、非流力学的モード（$\tau_\pi$ によって調節される）の寄与が流力学的モードの寄与を超えると、流体記述は失敗し、流れの観測量に急激な逸脱が生じる。

主要な貢献と結果

• ベンチマーク： このモデルは、中心性（$N_{track}$ で定義される）の関数としての荷電粒子の rapidity 分布、横運動量（$p_T$）分布、および平均横運動量に関する ALICE、ATLAS、CMS の実験データを成功裡に再現している。
• 揺らぎ解析： 周辺 p-Pb 衝突をシミュレーションし、$\tau_\pi$ を変化させることで、著者らは系が中心衝突から周辺衝突へと移動するにつれて楕円流の揺らぎが増加することを観測した。
• 開始の決定：
  • $b_\pi$ を変化させた場合、中心性範囲 75–85% から 70–80% の間で楕円流の揺らぎが急激に増加する。
  • $\eta/s$ を変化させた場合、揺らぎは 70–80% の中心性範囲以下で増加し始める。
  • これらの揺らぎは、$dN/dy \approx 7$ の荷電粒子多重度に対応する。
• 先行研究との比較： 著者らの以前の大型系に関する研究（開始が $dN/dy \approx 10$ で見つかった）とは異なり、p-Pb における遷移はより低い多重度で現れる。著者らは、推定された開始以下の逸脱が以前の研究よりも急激ではないと指摘しており、これは以前使用された ECHO-QGP に比べてより安定した MUSIC 流体力学ソルバーを使用していることによる可能性があると述べている。

意義と主張
本論文は、5.02 TeV における p-Pb 衝突で低次流体力学によって満足にモデル化できる最小の QGP 体積の見積もりを提供すると主張している。主要な結論は、流体力学的振る舞いがおよそ$dN/dy \approx 7$で破綻するということである。

著者らは、この遷移の性質について控えめなトーンを維持しており、これが鋭い開始を表すのか、それとも滑らかな遷移領域を表すのかは依然として未解決の問題であることを認めている。また、本研究はソフトセクターを分離していることを強調しており、開始の決定的な確認には、ハード（ジェットクエンチングなど）とソフトの観測量が流体力学的化の決定にどのように相互に影響するかを見るために、ハードセクターの同時調査が理想的であるとしている。この研究は、小型衝突系における流体力学モデルの有効性の限界を制約するための、非流力学的モード解析の具体的な応用として機能する。