SIS100 고밀도 영역에서 화학동결선의 체계적 분석

초록

본 연구는 SIS100 에서 진행되는 Au+Au 충돌(√sₙₙ = 3–5 GeV)의 화학동결 파라미터를 UrQMD 시뮬레이션과 Thermal‑FIST HRG 모델을 이용해 체계적으로 조사한다. 캐스케이드와 CMF 포텐셜을 포함한 두 가지 방정식 상태(EoS)를 비교하고, 포함된 하드론 종류(경량 핵·반양성자 포함 여부)에 따라 추출된 온도(T)와 화학퍼텐셜(μ_B), 그리고 평균 바리온 밀도(n_B)가 어떻게 변하는지를 분석한다. 결과는 하드론 선택과 EoS에 따라 T와 μ_B가 수십 MeV까지 차이날 수 있음을 보여주며, 화학동결선과 QCD 위상도(특히 임계점) 연결 시 주의가 필요함을 강조한다.

상세 분석

이 논문은 고밀도·저온 영역(μ_B/T ≈ 5–7)에서 화학동결선이 실제 물리적 위상 전이와 얼마나 연관될 수 있는지를 검증하기 위해 두 단계의 계산 체계를 구축한다. 첫 번째 단계는 Ultra‑relativistic Quantum Molecular Dynamics(UrQMD v4.0)를 이용해 Au+Au 충돌을 0–10 % 중심성으로 시뮬레이션하고, 두 가지 상이한 방정식 상태(EoS)를 적용한다. ‘Cascade’ 모드에서는 포텐셜이 없으며 매우 연한 EoS를 제공하고, ‘CMF’ 모드에서는 parity‑doublet chiral mean‑field 포텐셜을 포함해 압축성을 크게 높인 강직한 EoS를 구현한다. 이 두 시나리오는 특히 스트레인징(kaon)과 반양성자(anti‑p) 생산에 뚜렷한 차이를 만든다(그림 1). CMF EoS에서는 압축이 억제돼 K⁺와 anti‑p가 현저히 감소한다.

두 번째 단계는 Thermal‑FIST 프레임워크를 이용해 UrQMD에서 얻은 전산 하드론 수율을 HRG(Hadron Resonance Gas) 모델로 피팅한다. 피팅은 Grand‑Canonical Ensemble을 사용하고, γ_s 파라미터로 스트레인징 억제를, Q/A = 0.4의 전하‑바리온 비율을 고정한다. 중요한 점은 ‘누락된 하드론’ 효과를 검증하기 위해 세 가지 하드론 집합을 정의했다. Set 1은 π±, K±, p, Λ, anti‑p, 그리고 중성자(데우터론)까지 포함하고, Set 2는 데우터론을 제외, Set 3는 데우터론과 anti‑p까지 모두 제외한다. 모든 경우에 대해 10 %의 통계적 오차를 가정해 χ²/도프(dof)를 계산했으며, 전반적으로 χ²/도프 < 1을 기록해 매우 좋은 피팅 품질을 확인했다(그림 2). 특히 Set 1에서 낮은 에너지 구간에서 χ²가 약간 상승하는데, 이는 단순한 코얼레선스 모델이 저에너지에서 데우터론 생산을 충분히 설명하지 못하기 때문이다.

피팅 결과를 T–μ_B 평면에 나타내면(그림 3) 세 집합 모두 기존 실험 데이터와 파라미터화된 화학동결선(⟨E⟩/N ≈ 1 GeV) 근처에 위치한다. 그러나 하드론 선택에 따라 온도와 화학퍼텐셜이 체계적으로 변한다. anti‑p를 포함하면 T가 약 10–15 MeV 상승하고, 데우터론을 포함하면 추가로 5 MeV 정도 상승한다. 이는 반양성자가 초기 고온 단계에서 주로 생성되며, 경량 핵이 포함될 경우 전체 수율을 맞추기 위해 모델이 높은 온도를 요구하기 때문이다. 한편, CMF EoS를 사용할 경우 μ_B가 약 30–40 MeV 상승한다. 이는 강직한 EoS가 압축을 억제해 동일한 바리온 수를 달성하려면 더 높은 화학퍼텐셜이 필요하기 때문이다.

또한, HRG 모델을 이용해 피팅된 (T, μ_B)에서 바리온 밀도 n_B를 계산한 결과를 T–n_B 평면에 표시하였다(그림 4). CMF EoS에서는 동일 온도에서 n_B가 약 10 % 낮게 나오며, 이는 엔트로피/바리온(e/N) 값이 감소했기 때문이다. 반면 Cascade EoS는 더 높은 n_B를 예측한다. 이러한 차이는 실제 실험에서 추출된 n_B와 비교했을 때, EoS 선택이 물성 추정에 미치는 영향을 명확히 보여준다.

전체적으로 논문은 (1) UrQMD 자체가 비평형 동역학을 구현함에도 불구하고 HRG 피팅이 좋은 χ²를 제공한다는 점, (2) 피팅에 사용되는 하드론 종류와 시뮬레이션 EoS가 화학동결 파라미터에 수십 MeV 수준의 편향을 일으킨다는 점, (3) 따라서 화학동결선을 QCD 위상도(특히 임계점)와 직접 연결할 때는 이러한 시스템적 불확실성을 반드시 고려해야 함을 강조한다. CBM@FAIR 실험에서 기대되는 정밀도와 비교했을 때, 현재 수준의 모델링 불확실성은 임계점 탐색에 제한 요인으로 작용할 수 있다. 향후 연구에서는 (i) 더 정교한 코얼레선스 모델링, (ii) 다양한 중앙성 및 충돌 시스템에 대한 다중 시뮬레이션, (iii) 실험 데이터와의 직접적인 비교를 통한 γ_s 및 부피 파라미터의 제약이 필요하다.