Event-by-Event Multiplicity Fluctuations in Heavy-Ion Collisions Using Modified HIJING Monte Carlo Generator

この論文は、修正版 HIJING モンテカルロ生成器を用いて相対論的重イオン衝突における事象ごとの多重度揺らぎを解析し、生成された物質の状態（高温または低温）の識別、エネルギー損失モデルの検証、および QCD 相図上の臨界点探索における第一相転移の兆候検出への有効性を示しています。

要約

🍳 1. 研究の目的：宇宙の「最初のスープ」を作ってみる

宇宙が生まれた直後（ビッグバン直後）は、ものすごい熱と圧力で、原子の核（陽子や中性子）さえもバラバラになり、**「クォーク・グルーオンプラズマ（QGP）」**という、超高温・高密度の「スープ」のような状態だったと考えられています。

現代の科学者たちは、巨大な加速器（重イオン衝突実験）を使って、原子核同士を光の速さ近くまで加速してぶつけ、一瞬だけこの「スープ」を再現しようとしています。

でも、問題はここからです。
「本当にスープ（QGP）ができたのか？それともただの温かいお粥（通常の物質）なのか？」を、一瞬で消えてしまう現象から見分けるのは至難の業です。

そこでこの論文の著者たちは、**「出来上がった料理の『揺らぎ（ムラ）』」**に注目しました。

🌊 2. 核心のアイデア：料理の「ムラ」が語る真実

料理を作るとき、鍋の中が均一に温まっている場合と、一部が沸騰して泡立っている場合では、鍋の揺れ方や泡の立ち方が違いますよね？

この論文では、衝突して生まれた粒子の数（マルチピシティ）が、「イベントごと（一回の衝突ごと）」にどれくらい変動するかを詳しく分析しました。

• 通常の物質（冷たいお粥）： 粒子の数は比較的安定している。
• 新しい物質（熱いスープ/QGP）： 粒子の数が激しく揺れ動く。
• 相転移（状態変化）の瞬間： 氷が水に変わるときのように、物質の状態が混ざり合う瞬間には、「揺らぎ」が異常に大きくなるはずだ、と予想されます。

著者たちは、この「揺らぎの大きさ」を測ることで、**「今、鍋の中で何が起こっているか（QGP ができているか、相転移が起きているか）」**を見極めようとしています。

🎮 3. 使われたツール：HIJING という「料理シミュレーター」

実際に実験をするには、何百万回も衝突させる必要がありますが、それは現実的には不可能です。そこで、コンピューター上で衝突を再現する**「モンテカルロ・シミュレーター」**を使いました。

彼らが使ったのは**「HIJING」**という有名なシミュレーターですが、これには少し改造を加えました。

• 改造のポイント： 衝突した瞬間に、粒子が「熱いスープ（QGP）」の中を通過するときに、エネルギーを失う（減速する）仕組みを、よりリアルに組み込みました。
  • 熱い中を走る場合： 粒子は「放射」と「衝突」の両方でエネルギーを失います（まるで、熱い油の中で走ると、熱と摩擦でエネルギーを奪われるイメージ）。
  • 冷たい中を走る場合： 主に「放射」だけでエネルギーを失います。

この改造により、コンピューターの中で「QGP ができた場合」と「通常の物質の場合」を、まるでパラメータを切り替えるようにシミュレーションできるようになりました。

📊 4. 発見されたこと：揺らぎは「状態変化」のサイン

シミュレーションの結果、以下のようなことがわかりました。

1. 揺らぎは温度に敏感：
   衝突のエネルギーを上げて温度が上がり、QGP ができると、粒子数の「ムラ（揺らぎ）」が急激に大きくなりました。これは、「揺らぎの大きさ」が、鍋の中の温度や状態を測る「温度計」の役割を果たすことを示しています。

2. 相転移のサイン（1 次相転移）：
   物質が「冷たい状態」から「熱いスープ状態」へ劇的に変わる瞬間（1 次相転移）には、揺らぎが特に顕著に現れました。

   • 特に、「3 番目の揺らぎ」を「2 番目の揺らぎ」で割った値を見ると、この変化がより鮮明に現れることがわかりました。
   • これは、**「鍋の中で氷と水が混ざり合っている瞬間」**を、揺らぎのパターンから見分けることに成功したようなものです。

3. 見る範囲を広げるとより明確に：
   観測する粒子の範囲（角度や運動量）を広げると、相転移のサインがよりはっきりと現れました。

🏁 5. まとめ：なぜこれが重要なのか？

この研究は、単に「粒子が揺れた」という事実を報告しただけではありません。

• QCD 相図の地図作り： 物質の状態がどう変わるかを示す「地図（QCD 相図）」において、**「臨界点（Critical Point）」**と呼ばれる、相転移の性質が変わる重要な場所を探し出すための、非常に鋭い道具（揺らぎの分析）を提供しました。
• 実験の指針： 今後の実験（MPD や FAIR など）で、実際に「相転移のサイン」を見つける際に、どのデータに注目すればいいかの手がかりを与えています。

一言で言うと：
「原子核をぶつけてできた『宇宙の最初のスープ』が、本当にできているかどうかは、**『粒子の数のムラ（揺らぎ）』**を詳しく見ることで、鍋の中で何が起きているか（相転移かどうかも）見分けることができる！」という、新しい「料理の味見方法」を提案した論文です。

技術概要

以下は、提示された論文「Event-by-Event Multiplicity Fluctuations in Heavy-Ion Collisions Using Modified HIJING Monte Carlo Generator（修正 HIJING モンテカルロ生成器を用いた重イオン衝突におけるイベントごとの多重度揺らぎ）」の技術的サマリーです。

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

• クォーク・グルーオンプラズマ（QGP）の探索: 近年の高エネルギー実験の主要な目標の一つは、QCD（量子色力学）の漸近的自由性に由来する QGP の存在確認と、その相転移の秩序（一次相転移か、クロスオーバー転移か）の特定です。
• 実験的観測の難しさ: 有限の正味バリオン密度における QCD 物質の性質を第一原理から計算することは困難であり、相転移の順序や臨界点（CEP）の特定は実験に依存しています。
• 課題: 相転移（特に一次相転移）に敏感であり、かつ実験的に測定可能な観測量の確立が必要です。特に、臨界点近傍での混合相（ハドロンと QGP の気泡）の形成に伴う大きな揺らぎを検出できる手法が求められています。

2. 手法 (Methodology)

本研究では、重イオン衝突で生成される物質の状態を診断するツールとして、「イベントごとの荷電ハドロン多重度揺らぎ」を分析しました。

• シミュレーションツールの改変:
  • 標準的な輸送モデルであるHIJING (v.1.41) モンテカルロ生成器を改変しました。
  • 改変の核心は、ホットな媒体（QGP）とコールドな媒体における部分子（パトン）のエネルギー損失モデルを統合して実装した点にあります。
• エネルギー損失モデル:
  • ホットな媒体: 放射によるエネルギー損失（BDMPS 定式化）と、媒体粒子との衝突によるエネルギー損失の両方を考慮しました。
    • 放射損失：$dE/dz|_{radiative}$（有限・無限媒体のサイズ依存性を考慮）。
    • 衝突損失：$dE/dz|_{collisional}$（温度 $T$、デバイ遮蔽長 $\mu$ などに依存）。
  • コールドな媒体: 放射によるエネルギー損失のみを考慮。
• 一次相転移のシミュレーション:
  • 一次相転移を仮定し、ビナodal（二相共存）エネルギー点 $(\sqrt{s_{NN}})_b$ 付近での相転移確率 $\omega$ を導入しました。
  • 転移確率は、密度や温度などの揺らぎの強度 $\sigma^2$ をパラメータとする誤差関数（erf）分布で記述し、イベントごとに相転移の有無を確率的に決定しました。
• 分析対象:
  • 中心衝突（0-5%）の Au+Au 衝突を $\sqrt{s_{NN}} = 20 \sim 200$ GeV の範囲でシミュレーション（各エネルギー点で 350 万イベント）。
  • 観測量として、荷電ハドロン（$\pi, K, p$ など）の多重度 $N$ に関する累積量（Cumulant）を計算しました。
    • 平均 $\langle N \rangle$
    • 第 2 累積量 $C_2^N = \langle (\delta N)^2 \rangle$
    • 第 3 累積量 $C_3^N = \langle (\delta N)^3 \rangle$
    • 比：$C_2^N / C_1^N$（分散/平均）、$C_3^N / C_2^N$（歪度関連）
  • 運動学的カット：$|y| < 0.5$ または $1.0$、$p_T$ 範囲を変えて分析。

3. 主要な貢献 (Key Contributions)

• 揺らぎと媒体状態の相関の解明: 部分子のエネルギー損失モデルを HIJING に組み込むことで、生成される物質の種類（ホットかコールドか）と揺らぎの振る舞いの関係を定量的に示しました。
• 一次相転移の検出手法の提案: 累積量の比（特に $C_3/C_2$）が、ビナodal 点や一次相転移の存在に対して非常に敏感なシグナルとなることを示しました。
• 運動学的領域の影響評価: 観測する運動学的領域（$y$ や $p_T$ の範囲）を広げることで、相転移近傍での揺らぎの差異がより顕著になることを明らかにしました。

4. 結果 (Results)

• エネルギー依存性:
  • 荷電ハドロン多重度の揺らぎ（特に累積量の比）は、温度の上昇（すなわち、ホットで高密度な物質における部分子のエネルギー損失の増大）とともに増加することが確認されました。
  • この効果は、衝突エネルギーが高くなるほど顕著になります。
• モデル間の比較:
  • 「完全なエネルギー損失モデル（inf.+finite）」や「BDMPS モデル」は、「コールド物質」や「クエンチングなし（no quenching）」のケースと比較して、揺らぎが大幅に増大しました。
  • 一次相転移（1PT）を仮定した場合、特に $C_3^N / C_2^N$ の比において、ビナodal 点付近で明確な応答（ピークや急激な変化）が観測されました。
• 累積量の比較:
  • 一次相転移の検出には、第 2 累積量と第 1 累積量の比だけでなく、第 3 累積量と第 2 累積量の比（$C_3/C_2$）の方がより鋭敏な応答を示すことが分かりました。
• 運動学的領域の影響:
  • 観測範囲を拡大する（例：$|y|$ や $p_T$ の下限を広く取る）と、相転移近傍における揺らぎの挙動の違いがさらに明確になりました。

5. 意義と結論 (Significance)

• QCD 相図の描画への寄与: 本研究で示された「イベントごとの多重度揺らぎ」は、QCD 相図をマッピングし、臨界点（Critical Point）や一次相転移の存在を探索する実験データ（PHENIX, STAR, MPD, FAIR など）を解釈する上で極めて重要な観測量となります。
• モデルの妥当性検証: 揺らぎの大きさや性質を分析することで、エネルギー損失モデルの妥当性を検証し、衝突で生成された媒体がホットな QGP なのかコールドなハドロンガスなのかを識別する手段を提供します。
• 将来の実験指針: 第 3 累積量と第 2 累積量の比、および広範な運動学的領域での測定が、相転移のシグナルを捉えるために有効であることを示唆しており、今後の実験設計やデータ解析の指針となります。

要約すると、本論文は改変された HIJING 生成器を用いて、部分子エネルギー損失と相転移モデルを統合し、荷電ハドロン多重度の揺らぎ（特に高次累積量の比）が、QGP 生成や一次相転移の検出に対する極めて感度の高いプローブであることを実証したという点に最大の意義があります。