Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

기술 요약: 중이온 충돌에서의 베이지안 바리온 수 서셉터빌리티(Susceptibilities) 추출

문제 정의
보존된 전하, 특히 바리온 수의 요동은 양자 색역학(QCD) 상전이 구조, 특히 추정된 임계점(Critical Point, CP) 근처의 물리적 특성을 탐구하는 데 민감한 프로브 역할을 한다. 제일원리 격자 QCD 계산은 바리온 화학 포텐셜이 0($\mu_B = 0$)인 경우 바리온 수 서셉터빌리티($\chi_B^n$)를 정확하게 결정하지만, 유한 밀도에서의 거동은 외삽 또는 대안적인 이론적 접근 방식에 의존한다. 실험적으로는 RHIC의 BES(Beam Energy Scan) 프로그램이 Au-Au 충돌에서 순 양성자(net-proton) 수의 이벤트별 요동을 측정하여 이 영역을 조사한다. 그러나 이러한 데이터의 정량적 해석은 여러 비평형 및 실험적 효과로 인해 복잡해진다. 즉, 측정된 양성자는 전체 바리온 수의 일부만을 반영하며, 실험적 수용 영역(acceptance)은 운동량 공간에서 제한적이고, 정확한 전역 바리온 수 보존 법칙은 그랜드 캐노니컬(grand-canonical) 기대치와 비교하여 요동을 비자명하게 변형시킨다. 기존의 전략들은 밑바닥의 서셉터빌리티를 직접 추론하기보다는 특정 모델 계산(예: Hadron Resonance Gas, 배제 부피를 포함한 유체역학, 또는 수송 모델)을 데이터와 대조하는 방식을 취해왔다.

방법론
본 연구는 저자들의 지식 범위 내에서, 중이온 충돌 데이터로부터 바리온 수 서셉터빌리티($\chi_B^2, \chi_B^3, \chi_B^4$)를 직접 추출하는 최초의 베이지안 추출법을 제시한다. 이 프레임워크는 이론적 서셉터빌리티로부터 측정된 양성자 컴퓰런트(cumulants)로 가는 순방향 사상(forward map)을 네 단계로 구축한다:

1. 유체역학적 배경(Hydrodynamic Background): 3D 초기 조건과 일정한 에너지 밀도($\epsilon_{sw} = 0.26$ GeV/fm$^3$)에서의 입자화(particlization)를 사용하는 MUSIC 코드를 사용하여 이벤트 평균 유체역학적 초표면(hypersurface)을 생성한다.
2. 최대 엔트로피 동결(Maximum-Entropy Freeze-out, MaxEnt): 초표면 상의 국소 바리온 서셉터빌리티를 결합 바리온-반바리온 컴퓰런트로 변환한다. 본 프레임워크는 국소 서셉터빌리티 입력을 이상적인 Hadron Resonance Gas(iHRG) 베이스라인에 대한 비율로 매개변수화한다: $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$. MaxEnt 처리는 유체역학적 장(field)과 그 요동을 해면(hadron) 다중도 기약 상대 컴퓰런트로 매핑한다.
3. 수용 및 필터링(Acceptance and Filtering): Cooper-Frye 수용 및 바리온-양성자 필터링은 독립적인 이항 샘플링을 통해 결합 컴퓰런트를 수용된 양성자와 반양성자로 매핑하며, 이는 STAR 검출기의 운동학적 수용 영역을 고려한다.
4. 정확한 보존(Exact Conservation): SAM-3.0(Subensemble Acceptance Method)은 정확한 전역 바리온 수 보존을 부과하여, 그랜드 캐노니컬 컴퓰런트를 실험 데이터와 비교 가능한 캐노니컬 컴퓰런트로 변환한다.

베이지안 추론은 BES-II 범위의 특정 충돌 에너지($\sqrt{s_{NN}} = 7.7, 9.2, 11.5, 14.6, 17.3, 19.6, 27$ GeV) 내의 0–5% 중심 Au-Au 충돌에 대한 STAR 데이터를 대상으로 수행된다. 두 가지 상보적인 추출이 (i) 양성자 팩토리얼 컴퓰런트 비율($\hat{C}_n/\hat{C}_1$)과 (ii) 순 양성자 컴퓰런트 비율($\kappa_n/\kappa_1$)을 사용하여 수행된다. 서셉터빌리티 비율에는 평평한 사전 분포(flat priors)가 적용되며, 통계적 및 계통적 불확실성을 사용하여 가우시안 우도(likelihood)가 구성된다.

주요 결과

• 2차 서셉터빌리티 ($\chi_B^2$): 분석 결과, iHRG 값으로 정규화된 2차 서셉터빌리티($\chi_B^2/\bar{\chi}_B^2$)에 대해 엄격한 제약을 얻었으며, 상대적 불확실성은 5–10%이다. 추출된 값은 측정된 에너지 전반에 걸쳐 비단조적인 충돌 에너지 의존성을 보인다: 가장 낮은 에너지($\sqrt{s_{NN}} = 7.7$ GeV)에서 증가($\sim 1.06–1.07$), 중간 에너지($11.5–19.6$ GeV)에서 최소값($\sim 0.73–0.82$), 그리고 높은 에너지에서 평탄한 억제($\sim 0.8–0.9$)를 나타낸다.
• 격자 QCD와의 비교: 화학적 동결(chemical freeze-out) 라인으로 매핑했을 때, 추출된 $\chi_B^2$ 값은 $\mu_B \lesssim 300$ MeV 영역에서 격자 QCD 기반 추정치(특히 4D-TExS 상태 방정식)와 정량적으로 일치한다. 그러나 더 큰 $\mu_B$( $\sqrt{s_{NN}} = 7.7$ GeV에 해당)에서는 추출된 값이 격자 밴드에서 상향 이탈하여 최대 $\sim 50\%$의 증가를 보인다.
• 고차 서셉터빌리티: 3차 및 4차 서셉터빌리티($\chi_B^3, \chi_B^4$)는 바리온에서 양성자로의 매핑 효율 손실과 바리온 보존에 의한 왜곡으로 인해 약하게만 제약된다. 4차 비율은 $\sqrt{s_{NN}} = 19.6$ GeV에서 음의 값을 선호하는 경향($\chi_B^4/\chi_B^2 = -33 \pm 14$)을 보이지만, 큰 불확실성으로 인해 이 고차 모멘트들만으로는 QCD 상전이 구조에 대한 확정적인 결론을 내릴 수 없다.
• 상관관계 지배성(Correlation Dominance): 오직 기약 2-바리온 상관관계만을 가정하는 최소 시나리오($\Delta\chi_B^3 = \Delta\chi_B^4 = 0$)는 $\sqrt{s_{NN}} \approx 11$ GeV 근처에서 관찰되는 양성자 팩토리얼 컴퓰런트 비율 $\hat{C}_3/\hat{C}_1$의 비단조적 피크를 성공적으로 재현한다. 이는 관찰된 구조가 기약 3- 또는 4-바리온 상관관계보다는 $\chi_B^2$의 에너지 의존성과 정확한 바리온 수 보존 사이의 상호작용에서 기인함을 시사한다.

의의 및 주장
본 논문은 특정 상태 방정식(EoS)에 대한 가정 없이, 화학적 동결 라인을 따라 바리온 수 서셉터빌리티를 결정하는 최초의 데이터 기반 베이지안 결정법을 제공한다고 주장한다. 저자들은 자신들의 프레임워크가 2차 서셉터빌리티를 성공적으로 분리해냈으며, 낮은 바리온 밀도에서는 격자 QCD와 일치하지만, 높은 밀도($\mu_B \gtrsim 300$ MeV)에서는 비임계적 외삽에 비해 바리온 수 요동이 증가할 가능성을 보여준다고 주장한다.

연구는 관찰된 현상이 2차 이상의 기약 다중 바리온 상관관계 없이도 상당 부분 설명될 수 있다고 겸허히 결론짓는데, 이는 이들의 존재를 부정하는 것이 아니라 현재의 측정치가 주로 2-바리온 기여도에 민감하다는 것을 의미한다. 저자들은 낮은 충돌 에너지에서 관찰된 증가 현상이 향후 고정 타겟(fixed-target) 데이터에 의해 면밀히 검토되어야 한다고 언급하며, 현재의 추출은 고정된 유체역학 및 동결 설정에 조건부적임을 밝힌다. 본 프레임워크는 유한 밀도 QCD 상태 방정식을 제약하기 위한 일반적인 도구로 제시되며, 모델링 불확실성이 더 잘 통제되면 고정 타겟 에너지, 부피 요동 보정 및 기타 보존된 전하들로 확장될 수 있다.