Elliptic flow of charged hadrons in d+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV using a multi-phase transport model

本論文は、AMPT モデルを用いて 200 GeV の d+Au 衝突における荷電ハドロンの楕円流を解析し、初期のパートン散乱がその形成に決定的な役割を果たす一方、後期のハドロン再散乱の影響は限定的であることを示し、実験データとの比較を通じて非対称系における集団的挙動の理解を深めたものである。

要約

1. 実験の舞台：「デコボコな衝突」

通常、原子核同士（金と金など）を正面からぶつけると、円筒形のようなきれいな衝突になります。しかし、今回の実験は**「重水素（d）」**という、プロトン 1 個と中性子 1 個の「小さな双子」のような粒子と、大きな「金（Au）」の原子核をぶつけるものです。

• アナロジー：
  • 大きな「金」の原子核は、**「巨大な丸いクッション」**です。
  • 小さな「重水素」は、**「指先」や「小さな石」**のようなものです。
  • この「指先」が「クッション」に当たると、衝突の形はいつも同じではなく、**「デコボコ」**になります。

この「デコボコ」な衝突から、極限まで熱く高密度になった物質（クォーク・グルーオンプラズマ）が生まれます。研究者たちは、この衝突で粒子がどのように飛び散るかを調べることで、その中身が「液体のように流れているか」を確認しようとしています。

2. 何を測ったのか？「楕円形の流れ（楕円流）」

衝突が起きると、粒子は四方八方に飛び散りますが、実は**「楕円形（ひし形に近い形）」に偏って飛び散ることがあります。これを「楕円流（v2）」**と呼びます。

• アナロジー：
  • 風船を指で押して潰すと、押した方向には縮み、横方向には膨らみますよね。
  • 衝突の瞬間、この「デコボコ」な圧力によって、粒子が**「横に押し広げられるように」**流れる現象です。
  • この「流れの強さ」を測ることで、衝突した物質が「硬い砂」なのか、「滑らかな液体」なのかを判断します。

3. 使った道具：「AMPT モデル」というシミュレーター

実験自体はアメリカの加速器（RHIC）で行われますが、この論文では**「AMPT モデル」**という、コンピューター上の「仮想実験室」を使ってシミュレーションを行いました。

• 2 つのモード：
  1. デフォルトモード： 衝突直後は「粒子（ハドロン）」のまま動き、後から相互作用するシミュレーション。
  2. ストリング・メルトモード： 衝突直後に「粒子」が溶けて「クォーク（素粒子）」のスープ（液体）になり、その後また固まるシミュレーション。

研究者たちは、この 2 つのモードを使い分け、**「どの段階で『液体』のような流れが生まれるのか」**を探りました。

4. 発見した重要なこと

① 「クォークの衝突」が流れを作る鍵

シミュレーションの結果、**「クォーク同士がぶつかる頻度（断面積）」**を変えると、流れの強さが大きく変わることがわかりました。

• 発見： クォーク同士が頻繁にぶつかる（衝突確率が高い）ほど、「液体のような流れ」が強くなることが確認できました。
• 意味： 衝突の直後、物質は「クォークとグルーオンの液体」の状態になっており、その中で激しく相互作用しているからこそ、あのきれいな流れが生まれていることが示唆されました。

② 「最後の衝突」はあまり影響しない

一方、粒子が固まった後の「ハドロン（粒子）」同士の衝突（再散乱）の時間を長くしても、流れの強さはほとんど変わりませんでした。

• 意味： 流れの形を決めるのは、**「衝突の直後（クォークの段階）」**であり、その後の粒子同士のこすれ合いはあまり重要ではないことがわかりました。

③ 実験データとの比較：「見方」によって答えが変わる

実際の実験データ（STAR や PHENIX という実験グループのもの）と比較したところ、面白い結果が出ました。

• 現象： シミュレーションで「流れ」を計算する際、**「どの基準（角度）で測るか」**によって、実験データとの合致度が大きく変わりました。
• アナロジー：
  • 風船を横から見たら丸く見えますが、斜めから見ると楕円に見えます。
  • これと同じで、**「衝突の中心（参加核子）」を基準に測るか、「観測された粒子の分布」**を基準に測るかで、計算結果が異なります。
• 結論： 非対称な衝突（指先とクッション）では、実験データを解釈する際に非常に注意が必要で、単に「合っている・合っていない」ではなく、**「どの視点で見たか」**を考慮しないと誤解を招く可能性があります。

5. まとめ：この研究が伝えたかったこと

この研究は、**「小さな粒子（重水素）と大きな原子核（金）をぶつけただけでも、巨大な原子核同士をぶつけたときと同じように、『クォークの液体』が生まれている可能性がある」**ことを示唆しています。

• 重要なメッセージ：
  • 衝突の直後に、粒子が溶けて「クォークのスープ」になり、そこで激しくぶつかり合うことが、あの美しい「流れ」を生み出している。
  • しかし、非対称な衝突では、データの読み方が非常に繊細で、**「視点（基準）」**を間違えると、液体の性質を過大評価したり過小評価したりしてしまう危険性がある。

つまり、「宇宙の始まりのような極限状態の物質」を、小さな衝突でも再現できる可能性を、コンピューターシミュレーションを通じて詳しく検証した、非常に興味深い研究です。

技術概要

以下は、提供された論文「Elliptic flow of charged hadrons in d+Au collisions at √sNN = 200 GeV using a multi-phase transport model」の技術的な要約です。

論文概要

タイトル: 多相輸送モデルを用いた√sNN = 200 GeV における d+Au 衝突での荷電ハドロンの楕円流
著者: J. Tanwar, I. Aggarwal, L. Kumar (Panjab University), V. Bairathi, S. Kabana (Universidad de Tarapacá)
日付: 2026 年 3 月 6 日（論文記載日付）

---

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

相対論的重イオン衝突実験（RHIC および LHC）において、クォーク・グルーオンプラズマ（QGP）の形成が確認されています。非対称な衝突系（p+Au, d+Au, 3He+Au など）においても、QGP 様の集団的挙動（コレクティビティ）が観測されていますが、そのメカニズムは完全には解明されていません。
特に、d+Au 衝突のような非対称系では、初期幾何学的な揺らぎ（サブ核子レベルの構造）が楕円流（$v_2$）の形成に重要な役割を果たすと考えられています。しかし、実験データ（STAR, PHENIX）と理論モデルの比較において、非フロー効果（non-flow effects）やイベントごとの揺らぎの影響をどう解釈するか、またハドロン化のメカニズム（ストリング融解など）が $v_2$ にどのように寄与するかについて、さらなる検討が必要です。

2. 手法とモデル (Methodology)

本研究では、RHIC における√sNN = 200 GeV の d+Au 衝突をシミュレートするために、**多相輸送モデル（AMPT: A Multi-Phase Transport model）**を使用しました。

• モデル設定:
  • AMPT-DEF（デフォルトモード）: パーティオンがストリングと再結合し、その後ストリングがハドロンに断片化するモード。
  • AMPT-SM（ストリング融解モード）: 初期状態のストリングがまずパーティオンに融解し、クォークの合体（coalescence）モデルを経てハドロン化するモード。QGP 様の相を明示的に含みます。
  • 1 億個の最小バイアス（minimum bias）d+Au 衝突イベントを生成しました。
• 解析手法:
  • 楕円流 ($v_2$) の測定: 中間ラピディティ（$|\eta_{c.m.}| < 0.9$）の荷電ハドロンについて、$p_T$ 依存性を解析。
  • イベント平面の定義:
    1. $\eta$-サブイベント平面法 ($\psi_2^{ep}$): 異なる擬ラピディティ領域（$-1.8 < \eta < -1.0$ と $1.0 < \eta < 1.8$）を用いて非フロー効果を抑制。
    2. 参加者平面法 ($\psi_2^{pp}$): 衝突に関与する核子の位置から定義される真の対称面。
  • パラメータ変異: パーティオン間散乱断面積（$\sigma_{qq}$: 1.5, 3.0, 6.0 mb）とハドロン再散乱時間（$t_{had}$: 0.6, 30 fm/c）を変化させて、$v_2$ への依存性を系統的に調査しました。
• 比較対象: STAR および PHENIX 実験の 0-10% および 0-5% 中心度での実験データと比較。

3. 主要な結果 (Key Results)

A. $p_T$ 依存性とイベント平面の違い

• $\psi_2^{ep}$（イベント平面）: $p_T$ の増加とともに $v_2$ が単調増加し、集団運動の寄与を示唆。
• $\psi_2^{pp}$（参加者平面）: $p_T \sim 2.0$ GeV/c まで増加した後、飽和する傾向を示す。
• 差異の原因: $v_2^{ep}$ と $v_2^{pp}$ の大きな差異は、d+Au 衝突におけるイベントごとの大きな揺らぎ（最終状態の荷電粒子数と参加核子数の変動）に起因すると考えられます。

B. 中心度依存性

• $v_2^{ep}$: $p_T < 2.0$ GeV/c では中心度依存性は弱いですが、$p_T > 2.0$ GeV/c では対称系（Au+Au）と同様に、中心衝突から外側衝突へ向かって $v_2$ が増加する傾向が見られました。
• $v_2^{pp}$: 全 $p_T$ 範囲で中心度依存性を示しますが、その傾向は対称系とは逆転しています。これもイベントごとの揺らぎの影響が強いことを示しています。

C. 粒子種依存性と NCQ スケーリング

• 低 $p_T$ 領域での質量順序付けは弱く、中間 $p_T$ 領域でのバリオン・メソン分離（baryon-meson splitting）は明確には観測されませんでした。
• 構成クォーク数（$n_q$）によるスケーリング（NCQ scaling）も確認されませんでした。これは、d+Au 衝突で形成される系において、対称系のような強いパートオン相互作用やクォーク合体によるハドロン化が支配的ではない可能性を示唆しています。

D. パーティオン散乱断面積とハドロン再散乱の影響

• パーティオン散乱断面積 ($\sigma_{qq}$): $\sigma_{qq}$ を増大させると $v_2$ の大きさは増加します。特に、AMPT-SM モデルで $\sigma_{qq} = 6.0$ mb とすると、STAR 実験の参加者平面に基づくデータとよく一致します。
• ハドロン再散乱時間 ($t_{had}$): $t_{had}$ を 0.6 fm/c から 30 fm/c に変化させても、$v_2^{ep}$ への影響はほとんど見られませんでした。一方、$v_2^{pp}$ は $t_{had}$ の増加に伴って増加しました。
• 実験データとの比較:
  • STAR データ (0-10%): AMPT-SM (6.0 mb) が参加者平面 ($\psi_2^{pp}$) に対して実験値を良く再現。
  • PHENIX データ (0-5%): AMPT-DEF モデルおよび AMPT-SM (1.5 mb) がイベント平面 ($\psi_2^{ep}$) に対して実験値を良く再現。
  • この不一致は、非フロー効果の除去手法やイベント平面の定義の違いに起因しており、非対称系における実験結果の解釈には注意が必要であることを示しています。

4. 結論と意義 (Significance)

本研究は、AMPT モデルを用いて d+Au 衝突における楕円流の特性を詳細に解明しました。

• 初期状態の重要性: 初期段階のパートオン相（特にパーティオン散乱断面積）が $v_2$ の形成に決定的な役割を果たしていることが示されました。
• 非対称系の特性: 対称系とは異なり、イベントごとの大きな幾何学的揺らぎが $v_2$ の測定値（特にイベント平面と参加者平面の差）に大きな影響を与えることが確認されました。
• モデルの妥当性: AMPT モデル（特にストリング融解モード）は、適切なパラメータ設定により、非対称系における $p_T$ 依存性をよく記述できることを示しました。
• 実験解釈への示唆: 実験データと理論モデルの比較において、非フロー効果の除去やイベント平面の定義が結果に大きく影響するため、非対称系におけるコレクティビティの解釈には慎重なアプローチが必要であるという重要な示唆を与えています。

この研究は、RHIC における非対称衝突系における QGP 様の物質の性質理解と、ハドロン化メカニズムの解明に貢献するものです。