Microscopic study of baryon stopping in low-energy heavy-ion collisions within UrQMD model

技術要約：UrQMDモデルを用いた低エネルギー重イオン衝突におけるバリオン停止の微視的研究

問題提起
低エネルギー重イオン衝突において、バリオン停止（初期核から中間ラピディティ領域へのバリオン数の輸送）は、生成されたファイアボールの正味バリオン密度およびバリオン化学ポテンシャル（$\mu_B$）を決定するための極めて重要なプロセスである。正味陽子のラピディティ・スペクトルに関する実験的測定は、停止現象に対する制約を与えるものの、それら陽子の微視的な起源までは明らかにしない。具体的には、中間ラピディティにおける正味陽子が、初期の陽子のみに由来するのか、あるいは初期の中性子が電荷交換反応（例：$p + \pi^- \leftrightarrow n + \pi^0$）などを通じてアイソスピン変換を起こしたものからかなりの割合で生じているのかを区別することは、実験的には不可能である。この曖昧さは、バリオン停止の解釈や、ゆらぎ解析で使用される化学平衡の仮定（例：北澤・浅川の定式化）の妥当性の抽出を困難にしている。

手法
本研究では、Ultra-relativistic Quantum Molecular Dynamics (UrQMD) トランスポートモデル（バージョン3.4）を用い、広範なビームエネルギー範囲 $\sqrt{s_{NN}} = 2.4 - 17.3$ GeV（$E_{lab} = 1.23A - 158A$ GeV）における最小バイアス Au+Au 衝突を調査する。解析には以下の2つのモードを使用する：

1. 平均場（MFD）モード： 核ポテンシャルを考慮するために硬いスカーム（Skyrme）状態方程式を用いた、$\sqrt{s_{NN}} < 3.2$ GeV（$E_{lab} < 4A$ GeV）以下のエネルギー用。
2. カスケード（Cascade）モード： すべてのエネルギーで使用され、かつ $\sqrt{s_{NN}} > 3.2$ GeV において排他的に使用される。この領域では、クォークやグルーオンの自由度は解像されず、粒子生成はストリングの励起と断片化によってモデル化される。

中心的な手法の革新は、モデルの出力ファイル（f13およびf15）から完全な相互作用履歴を再構成することにより、最終状態の陽子を微視的に分類することにある。陽子は以下のように分類される：

• 生成された陽子（Produced）： 衝突の進化過程におけるストレンジ粒子やメソンに由来するもの。
• 輸送された陽子（Transported）： 初期の核子に由来するもの。これらはさらに以下に細分される：
  • アイソスピン保存型： 初期の陽子に由来するもの。
  • アイソスピン変換型： アイソスピン状態が変化して陽子となった初期の中性子に由来するもの。

本研究では、3つの衝撃パラメータ間隔（0–4.7 fm、4.7–9.3 fm、9.3–11.4 fm、それぞれ中心衝突、中位中心衝突、周辺衝突に対応）について、各エネルギー点につき20,000イベントを解析する。

主要な貢献と結果

• アイソスピン変換の定量化： 本研究は、中間ラピディティにおける輸送された陽子のかなりの割合が、初期の中性子に由来することを明らかにしている。アイソスピン変換率（輸送された陽子のうち、当初は中性子であったものの割合）は無視できない値であり、調査したエネルギー全域にわたって0.4から0.6の間で推移している。この率はエネルギーとともに増加し、$\sqrt{s_{NN}} \approx 4$ GeV 付近で飽和して Au 核の中性子・核子比（$N/A \approx 0.6$）に近づく。これは、中心衝突において系が完全なアイソスピン無秩序化に近づいていることを示唆している。
• エネルギーおよび中心度依存性：
  • 低エネルギー（$\sqrt{s_{NN}} < 4$ GeV）では、パイオン密度の制限により、アイソスピン変換率は低くなる。
  • $\sqrt{s_{NN}} \approx 4$ GeV 以上では、率は飽和する。この飽和点は核の透過性の発現と一致しており、アイソスピンの無秩序化とバリオン停止は結合した現象であることを示している。
  • 変換率は中心度にも敏感であり、周辺衝突ではハドロン相の有限な寿命と電荷交換機会の制限により、より低い値を示す。
• ラピディティ損失とスペクトル： アイソススピン変換された陽子は、アイソスピン保存型の陽子と比較して大きなラピディティ損失を示し、前者に見られる前方ピークを持たない。正味陽子の平均ラピディティ損失は、アイソスピン変換成分によって支配されている。ラピディティ・スペクトルの減少曲率（$C_{red}$）は、低エネルギーにおいて陽子の起源に基づいた明確な順序を示すが、$\sqrt{s_{NN}} \approx 4$ GeV 以上ではこの傾向は消失する。
• 流動係数：
  • 直接流（$v_1$）： $\sqrt{s_{NN}} \approx 4$ GeV 以下において、明確な起源に基づく順序が観察される。アイソスピン保存型の陽子が最も強い流動を示し、次いで全陽子、アイソスピン変換型、生成された陽子の順となる。MFDによる予測は、一般にCascadeモードよりも高い傾きを与える。
  • 楕円流（$v_2$）： 同様の順序が観察される。MFDによる予測は中心衝突における負の $v_2$ を再現するが、Cascadeによる予測は正のままであり、低エネルギーにおける核ポテンシャルの重要な役割を浮き彫りにしている。
• 観測量および定式化との比較：
  • パイオン比： $\pi^-/\pi^+$ 比とアイソスピン変換率の間の直接的な関連性が確立された。$\sqrt{s_{NN}} < 4$ GeV では、比率が上昇する（中性子過剰な系を反映する）一方で変換率は低下する。このエネルギーを超えると、両者はともに飽和傾向を示す。
  • 北澤・浅川（KA）定式化： 本研究では、アイソスピン変換率をKA定式化の $\alpha_N$ および $\alpha_\pi$ パラメータと比較する。$\sqrt{s_{NN}} \lesssim 10$ GeV 以下において、両パラメータの間に顕著な乖離が見られ、化学平衡（$\alpha_N = \alpha_\pi$）の仮定がこのエネルギー領域において疑わしいことを示している。

意義および主張
本論文は、中間ラピディティにおける従来の正味陽子の測定には、本質的にかなりの割合のアイソスピン変換成分が含まれており、その大きさは衝突エネルギーおよび中心度によって変化することを主張している。この知見は、高 $\mu_B$ 領域で作動するFAIRやNICAのような施設における、バリオン停止の解釈方法の再評価を必要とするものである。

本研究は、バリオン停止とアイソスピン力学のハドロン的なベースラインを確立し、アイソスピンの無秩序化とバリオン停止が結合した現象であり、$\sqrt{s_{NN}} \approx 4-6$ GeV 付近の核透過性の発現と同時に飽和することを実証している。著者らは、これらの結果が、より高いエネルギーにおける非決定論的なクォーク・グルーオンプラズマ（QGP）相の寄与を特定するための参照点を提供すると述べている。論文は、アイソスピン無秩序化の飽和が、非決定論的な相の発現に先行するのか、一致するのか、あるいは後れるのかという問いを投げかけており、これは非決定論的な相を取り入れたさらなる調査を必要とする問題であると結論付けている。

限界
著者は、本研究が純粋なハドロン輸送フレームワーク（UrQMD 3.4）に限定されており、非決定論的な相を含んでいないことを認めている。また、陽子の分類はシミュレーションされた相互作用履歴に依存しており、これらは実験では直接アクセスできないため、知見はパイオン比や流動係数などの観測量を通じた間接的な比較に制約される。さらに、UrQMD 4.0には含まれているが3.4には含まれていない軽核生成の影響については、結果に対して最小限であると述べている。