[論文レビュー] The PHOBOS Glauber Monte Carlo

B. Alver, M. Baker+2 more|ArXiv.org|May 28, 2008

High-Energy Particle Collisions Research被引用数 136

ひとこと要約

本論文は、重イオン衝突における初期状態の幾何構造をシミュレートするためのPHOBOSモンテカルロゲーラーモデルを提示する。確率的核子配置と衝突パrameterのサンプリングを用いて、$N_{\text{part}}$、$N_{\text{coll}}$、および非対称性を計算する。主な貢献は、詳細なイベントごとの観測量を備えた、公開可能で使いやすい実装であり、RHICおよびLHCにおける相対論的重イオン物理学におけるフラクチュエーション駆動効果の正確な研究を可能にする。

ABSTRACT

``Glauber'' models are used to calculate geometric quantities in the initial state of heavy ion collisions, such as impact parameter, number of participating nucleons and initial eccentricity. The four RHIC experiments have different methods for Glauber Model calculations, leading to similar results for various geometric observables. In this document, we describe an implementation of the Monte Carlo based Glauber Model calculation used by the PHOBOS experiment. The assumptions that go in the calculation are described. A user's guide is provided for running various calculations.

研究の動機と目的

• 重イオン衝突における初期状態幾何構造の詳細で再現可能なモンテカルロゲーラーモデルの実装を提供すること。
• PHOBOS解析に適した、$N_{\text{part}}$、$N_{\text{coll}}$、および非対称性（$\epsilon_{\text{part}}$、$\epsilon_{\text{RP}}$）を含むイベントごとの幾何観測量の正確な計算を可能にすること。
• 実験的および理論的応用の両方を想定した、使いやすく拡張可能なコードベースに加え、チュートリアルおよびドキュメントを提供すること。
• 参加核子配置のイベントごとのサンプリングを可能にすることで、初期状態のフラクチュエーションに関する研究を支援すること。
• PHOBOS特有のゲーラーモデル実装を標準化し文書化することで、解析の一貫性を確保するとともに、他の実験との共同作業を支援すること。

提案手法

• 電子散乱データから導かれる半径確率分布を用いて、核内の核子の確率的配置を実施。電荷密度にはフェルミ分布またはガウス分布を採用。
• 物理的に不自然なクラスタリングを避けるために、最小核子間隔（$d_{\text{min}}$）を導入。デフォルト値は0.4 fm。
• 最大衝突パラメータ $b_{\text{max}} \approx 20$ fm まで、$b$-重み付き分布から衝突パラメータ $b$ をランダムにサンプリング。核はそれぞれ $(-b/2, 0, 0)$ および $(b/2, 0, 0)$ に配置。
• 二体衝突は横方向距離に基づき決定：二つの核子が互いに「ボール径」$D = \sqrt{\sigma_{\text{NN}} / \pi}$ よりも近い場合に衝突とみなす。
• アイケオナル近似の適用：核子はビーム方向に直線的に進行し、縦方向の力学的過程は無視。
• 参加核子数 $N_{\text{part}}$、衝突回数 $N_{\text{coll}}$、および非対称性計算に必要な参加核子位置のモーメントを、イベントごとに計算。

実験結果

リサーチクエスチョン

• RQ1 PHOBOSモンテカルロゲーラーモデルは、重イオン衝突におけるイベントごとの初期状態幾何構造をどのように生成するか？
• RQ2 核子位置の確率的配置および最小間隔の影響は、$N_{\text{part}}$ および $N_{\text{coll}}$ の分布にどのような影響を及えるか？
• RQ3 参加核子非対称性 $\epsilon_{\text{part}}$ と反応面非対称性 $\epsilon_{\text{RP}}$ は、異なる衝突系において $N_{\text{part}}$ に応じてどのように変化するか？
• RQ4 イベントごとの幾何フラクチュエーションは、RHICおよびLHCエネルギー下でのAu+Au、Cu+Cu、Pb+Pb衝突の初期状態にどの程度影響を及えるか？
• RQ5 標準化された公開可能なゲーラーモデル実装は、重イオン実験間での一貫性ある解析をどのように支援するか？

主な発見

• 本モデルは、RHICおよびLHCエネルギー下でのCu+Cu、Au+Au、Pb+Pb衝突において、$N_{\text{part}}$ および $N_{\text{coll}}$ の期待される分布を図2（10,000イベント）に示すように正確に再現している。
• 低 $N_{\text{part}}$ において、参加核子非対称性 $\epsilon_{\text{part}}$ は $\epsilon_{\text{RP}}$ よりも顕著に大きいことが判明し、初期状態に強い非ガウス的フラクチュエーションが存在することを示唆している。
• すべての衝突系において、$N_{\text{part}}$ が増加するにつれて $\epsilon_{\text{part}}$ は減少し、中央衝突における幾何構造の平滑化が顕著に現れている。
• コード実装は、$\epsilon_{\text{part}}$ の正確な計算を可能にし、式 $\epsilon_{\text{part}} = \sqrt{(\text{VarY}-\text{VarX})^2 + 4\text{VarXY}^2}/(\text{VarY}+\text{VarX})$ を用いる。イベントごとのモーメントは出力ntupleに格納されている。
• 合計核-核断面積は Run() 関数内で内部的に計算され、GetTotXSect() を通じてアクセス可能であり、イベントサンプル間で一貫した値を示している。
• $N_{\text{part}}$ および $N_{\text{coll}}$ の分布を理論的期待値と比較することでモデルの妥当性を検証した。コードはスタンドアロンまたは埋め込み解析ワークフローの両方で使用可能な公開コードである。