소형 핵 충돌을 위한 최신 Glauber 모델 예측: 산소·네온 시스템

초록

본 논문은 TGlauberMC v3.3의 업데이트 내용을 소개하고, LHC에서 5.36 TeV(OO, NeNe) 및 9.62 TeV(pO)에서 진행될 산소‑산소, 네온‑네온, 양성자‑산소 충돌에 대한 기하학적 단면, 참여핵수(Npart), 충돌수(Ncoll), 편심도(εₙ) 등 초기 상태 양자를 예측한다. 최신 핵밀도 프로파일(3pF, HO, NNLOsat)과 개선된 핵자‑하위구조 모델을 적용해 중간 급속도 다중생성량의 중심도 의존성을 제시하고, 코드와 사용자 가이드를 공개한다.

상세 분석

TGlauberMC v3.3은 기존 Monte‑Carlo Glauber 프레임워크에 두 가지 주요 개선을 도입한다. 첫째, 작은 핵(⁶⁶O, ²⁰Ne)의 전하밀도 데이터를 기반으로 3‑parameter Fermi(3pF), 수정된 조화진동자(HO) 및 최신 NNLOsat 계산에서 얻은 정확한 밀도 프로파일을 구현한다. 표 I에 제시된 파라미터들은 핵반지름 R, 확산도 a, 변형 파라미터 β 등을 포함해, 기존 ‘Odat’, ‘Opar’ 등과 비교해 rms 반지름 차이가 5 % 이내임을 확인한다. 둘째, 핵자‑하위구조를 보다 현실적으로 묘사하기 위해 여러 전이 함수(Γ‑parameterization, TRENTo, HIJING, PYTHIA‑Monash)를 도입한다. ω(0–1) 혹은 Gaussian 폭 w(0.4–0.7) 파라미터를 조정함으로써 핵자‑핵자 충돌 확률 p(bₙₙ)를 하드‑스피어 모델에서 벗어나 부드러운 형태로 전환한다(Fig. 1). 이러한 전이 함수 선택은 Ncoll·σₙₙ 관계와 초기 에너지 밀도 분포에 직접적인 영향을 미치며, 특히 작은 시스템에서의 통계적 변동성을 크게 확대한다.

논문은 OO와 NeNe 충돌에 대해 총 단면 σₕₐ𝑑𝑟𝑜𝑛을 5.36 TeV에서 각각 2.9 b와 3.2 b(±2 %) 수준으로 예측한다. 초기 상태 양자 Npart, Ncoll, ε₂, ε₃ 등은 중앙성(impact parameter) 의존성을 상세히 제시했으며, 특히 ε₂(NeNe)/ε₂(OO) 비율이 1.15 ~ 1.20 정도로, 이는 실험적으로 관측되는 v₂ 비율과 직접 연결될 수 있다. α‑클러스터링 효과를 포함한 Trajectum 기반 프로파일(TR_OXYGEN_V15, TR_NEON_V14)도 비교했으며, 클러스터링이 ε₂를 약 5 % 상승시키는 것으로 나타났다.

입자 생성에 대해서는 다중 부분입자 상호작용(MPI) 모델을 적용해, 기존 pp 및 pPb 데이터에 기반한 파라미터를 OO·NeNe에 확장한다. 중앙성에 따른 중간 급속도 다중생성량 ⟨dNch/dη⟩는 Npart에 거의 선형적으로 스케일링하며, 예측된 중앙성 구간(0–5 %)에서 ⟨dNch/dη⟩≈1200(OO), ≈1500(NeNe) 정도가 된다. 이는 기존 PbPb(5.02 TeV)와 비교해 약 30 % 낮은 값이지만, 작은 시스템에서도 QGP‑like 흐름을 관찰하기에 충분한 입자 밀도를 제공한다는 점을 강조한다.

코드 배포는 HepForge(버전 v3.3.2)에서 이루어지며, 사용자 가이드와 추가 라이브러리(Trajectum Nucleus Generator, TrNucGen)도 포함한다. 이는 실험 협업팀이 실시간으로 초기 조건을 교정하거나, 다양한 핵밀도 모델을 테스트할 수 있게 한다. 전체적으로 본 연구는 작은‑중간 규모 핵 충돌에서 초기 기하학과 입자 생성 메커니즘을 정밀히 연결함으로써, QGP 형성 임계조건 탐색에 필수적인 기준점을 제공한다.