[논문 리뷰] The PHOBOS Glauber Monte Carlo

B. Alver, M. Baker외 2명|ArXiv.org|2008. 05. 28.

High-Energy Particle Collisions Research인용 수 136

한 줄 요약

이 논문은 레이저 중이온 충돌에서 초기 상태 기하학을 시뮬레이션하기 위한 PHOBOS 몬테카를로 글라우버 모델을 제시한다. 무작위로 뉴클론 위치를 정하고 영향 매개변수를 샘플링하여 $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}}$보다 유의미하게 크며, 이는 초기 상태에서 강한 비정규 분포의 변동성이 있음을 시사한다.
• 모든 시스템에서 $\epsilon_{\text{part}}$는 $N_{\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}}$ 분포를 기대치와 비교하여 검증되었으며, 코드는 단독 또는 통합 분석 워크플로우에서 사용 가능한 공개된 자료이다.