[論文レビュー] The iEBE-VISHNU code package for relativistic heavy-ion collisions

Chun Shen, Zhi Qiu+4 more|arXiv (Cornell University)|Sep 29, 2014

High-Energy Particle Collisions Research参考文献 77被引用数 27

ひとこと要約

iEBE-VISHNUコードパッケージは、(2+1)次元の粘性流体力学とハドロンカスケードモデルを組み合わせたハイブリッドフレームワークを用いて、相対論的重イオン衝突のイベントごとのシミュレーションを可能にする。モンテカルロシミュレーションのためのモジュラでスケーラブルなパイプラインを提供し、熱光子、ジェットクエンチング、重いクォーク拡散といったレアプローブの数値的に効率的なインターフェースを備え、せん断粘性係数とエントロピー比 η/s などのクォーク- gluonプラズマの輸送特性を定量的に制約可能である。

ABSTRACT

The iEBE-VISHNU code package performs event-by-event simulations for relativistic heavy-ion collisions using a hybrid approach based on (2+1)-dimensional viscous hydrodynamics coupled to a hadronic cascade model. We present the detailed model implementation, accompanied by some numerical code tests for the package. iEBE-VISHNU forms the core of a general theoretical framework for model-data comparisons through large scale Monte-Carlo simulations. A numerical interface between the hydrodynamically evolving medium and thermal photon radiation is also discussed. This interface is more generally designed for calculations of all kinds of rare probes that are coupled to the temperature and flow velocity evolution of the bulk medium, such as jet energy loss and heavy quark diffusion.

研究の動機と目的

• 相対論的重イオン衝突のイベントごとのモデリングのための包括的でモジュラなシミュレーションフレームワークの開発。
• クォーク-グルーオンプラズマの性質を制約するため、η/s などの輸送係数を調整することでモデルとデータの比較を可能にする。
• レアプローブ（熱光子、ジェットエネルギー損失、重いクォーク拡散など）を進化するバルク媒体に結合する、数値的に効率的なインターフェースの実装。
• 初期条件のフラクチュエーションを伴う流体力学の数値的課題と、コペル-フライ采択における負の確率の問題に対処する。
• 統一されたワークフローとSQLiteベースのデータストレージを用いた大規模モンテカルロシミュレーションをサポートする。

提案手法

• パッケージは、KNOスケーリングに基づく多重度フラクチュエーションを用いて初期条件を生成する superMC モジュールを使用し、√sNN = 5.02 TeV における p+Pb の多重度分布を再現する。
• VISHNew は 400×400 グリッド上で (2+1)次元の粘性流体力学を実行し、せん断粘性係数補正を施した相対論的ナビエ-ストークス方程式を解く。
• iSS モジュールはコペル-フライの公式を用いて流体セルから最終状態の粒子を抽出し、負の確率問題を処理するための3つの方法を採用する。
• UrQMD は、freeze-out 後の最終状態相互作用と粒子崩壊をモデル化するハドロンカスケード段階をシミュレートする。
• テンソル分解法により、実験フレームのせん断応力と局所静止フレームの放射テンソルを組み合わせて、熱光子の放射率に粘性補正を加える。
• 繰り返しのローレンツブーストを回避するため、ローレンツ不変スカラー積 (πμνq̂μq̂ν) と射影演算子をインターフェースで使用し、効率的かつポータブルなレアプローブ計算を実現する。

実験結果

リサーチクエスチョン

• RQ1 小スケール系における多重度フラクチュエーションは、KNOスケーリングに整合する形でどのようにモデル化可能か？ これにより、実験的な p+Pb 多重度分布が再現可能か？
• RQ2 初期条件のフラクチュエーションが、重イオン衝突における回転対称性の度合い εn の分布に与える影響は何か？
• RQ3 高空間分解能を伴うフラクチュエーションを伴う初期条件から進化する粘性流体力学を、数値的に安定化するにはどうすればよいか？
• RQ4 イベントごとのシミュレーションにおけるコペル-フライ採択における負の確率問題を効果的に緩和する方法は何か？
• RQ5 熱光子のようなレアプローブは、粘性補正を施した進化する流体力学場に放射率を結合することで、どのように効率的に計算可能か？

主な発見

• superMC モジュールは、MC-グローバー・モデルに KNO スケーリングを組み込むことで、√sNN = 5.02 TeV における p+Pb 衝突の測定済み多重度分布を成功裏に再現した。
• VISHNew は (2+1)次元の粘性流体力学において数値的安定性と収束性を示し、ブースト不変流れに対するガブザーの解析解との整合性を確認した。
• iSS モジュールの3つのサンプリング手法は、負の確率問題の低減に効果的であり、性能ベンチマークではさまざまなシミュレーション設定で安定した粒子生成が確認された。
• 熱光子放射インターフェースは、放射率をローレンツスカラーに分解することで、繰り返しのテンソルブーストを回避し、粘性補正の高速で再利用可能なテーブル化を実現した。
• 1つの 0-5% Pb+Pb イベントについて、ラップトップ（2.4 GHz i5、4GB RAM）で合計約690秒のシミュレーション時間を達成し、流体力学が最も計算コストの高いコンponentであった。
• このフレームワークは、将来の拡張をサポートしており、非平衡ダイナミクス、ボリューム粘性、HBT相関計算が含まれる。新規モジュールは https://u.osu.edu/vishnu/ でリリースされる予定である。