핵방정식 연구를 위한 클러스터·하이퍼핵 프로빙: PHQMD 모델과 STAR 데이터 비교

초록

PHQMD 전이동역학 모델을 이용해 √sₙₙ=3 GeV Au+Au 충돌에서 측정된 STAR의 양성자·Λ·경량 클러스터·하이퍼핵 데이터와 세 가지 핵 방정식(EoS)을 비교하였다. 압축계수 K가 200 MeV인 연성(momentum‑dependent) EoS가 대부분의 관측량(수율, p_T·rapidity 분포, v₁·v₂ 흐름)과 가장 잘 일치했으며, 경성(EoS‑hard)도 유사한 경향을 보였다.

상세 분석

본 논문은 고밀도 핵물질의 상태 방정식(EoS)을 실험적으로 제약하기 위해, 최신 STAR 협력팀이 제공한 √sₙₙ=3 GeV Au+Au 충돌 데이터를 PHQMD(Parton‑Hadron‑Quantum‑Molecular Dynamics) 전이동역학 모델과 정밀 비교하였다. PHQMD는 QMD 기반의 입자 위치·운동량을 Gaussian 파동함수로 기술하고, 파트론·하드론 전이와 클러스터·하이퍼핵 형성을 동시에 다룰 수 있는 N‑body 접근법이다. 모델 내부에서 핵 평균장(mean‑field)은 세 가지 형태로 구현되었다. 첫 번째는 전통적인 Skyrme‑type 정적 포텐셜로, 압축계수 K=200 MeV(‘soft’)와 K=380 MeV(‘hard’) 두 버전을 사용하였다. 두 번째는 위 두 정적 포텐셜에 momentum‑dependence를 추가한 ‘soft‑MD’ 형태이며, 이는 p‑A 산란 실험에서 추출된 광학포텐셜 U_opt(p)를 기반으로 파라미터화하였다. 저자들은 U_opt(p)의 고에너지 영역(>1 GeV/c)에서 데이터가 부족함을 인식하고, 세 가지 서로 다른 외삽(parametrization I, II, III)을 도입해 민감도 검증을 수행하였다.

시뮬레이션 결과는 다음과 같이 정리된다. (1) 양성자와 Λ의 전산량, rapidity 및 p_T 스펙트럼은 ‘soft‑MD’ EoS가 가장 정확히 재현했으며, ‘hard’ EoS도 전반적인 형태는 맞추지만 절대값에서 약 10‑15% 차이를 보였다. (2) A≤4 클러스터(디튬, 트리톤, ³He, ⁴He)의 수율과 스펙트럼 역시 ‘soft‑MD’가 실험값과 거의 일치했으며, ‘hard’는 특히 고p_T 구간에서 과다 예측하였다. (3) 하이퍼핵(³ΛH, ⁴ΛH) 역시 ‘soft‑MD’가 적절히 설명했지만, 통계가 제한적이어서 확정적인 결론을 내리기엔 부족했다. (4) 흐름 계수 v₁과 v₂는 압축계수 K에 매우 민감한데, ‘soft‑MD’가 양성자·Λ 모두에서 STAR 측정값을 재현했고, ‘hard’는 v₁의 기울기를 과도하게 크게 만들었다. 특히 v₂(elliptic flow)의 경우, momentum‑dependent 포텐셜이 없을 때는 과소평가되지만, ‘soft‑MD’는 실험 데이터와 거의 일치하였다.

또한, 저자들은 momentum‑dependent 포텐셜의 고에너지 외삽이 관측량에 미치는 영향을 검증하였다. 파라미터화 I(지수형 감소)와 II(고p에서 상수) 사이의 차이는 주로 초기 고p_T 입자들의 분포에 영향을 주었으며, 최종 수율·흐름에는 큰 차이를 만들지 않았다. 이는 충돌 후 7 fm/c 이후 시스템이 거의 열평형에 가까워지면서, 초기 고p_T 부분이 평균장에 미치는 영향이 감소함을 의미한다.

결론적으로, 압축계수 K≈200 MeV에 momentum‑dependence를 포함한 ‘soft‑MD’ EoS가 현재 √sₙₙ=3 GeV 영역에서 가장 신뢰할 수 있는 핵 방정식임을 제시한다. 이는 기존의 ‘soft’와 ‘hard’ 구분만으로는 충분히 설명되지 않는 흐름·클러스터 데이터들을 동시에 만족시키는 최초의 전이동역학 연구라 할 수 있다. 향후 더 높은 통계와 다양한 입자(특히 다중 하이퍼핵) 측정이 이루어지면, momentum‑dependent 포텐셜의 정밀한 형태와 고밀도 물질의 비선형 응답을 더욱 구체화할 수 있을 것으로 기대된다.