Out-of-equilibrium contributions to charm hadrons in a fluid-dynamic approach

기술 요약: 유체역학적 접근법을 통한 참(charm) 하드론의 비평형 기여

문제 제기
참 쿼크는 거대한 질량 덕분에 주로 초기 강한 산란(hard scattering) 과정에서 생성되며, 화염구(fireball)의 수명 내에 열적 생성(thermal production)을 통해 열적 평형에 도달하지 않으므로, 쿼크-글루온 플라즈마(QGP)를 탐사하는 독특한 프로브 역할을 한다. 기존의 유체역학적 프레임워크는 참 쿼크의 역학을 국소 열적 평형에 근접해 있다는 가정하에 성공적으로 모델링해 왔다. 그러나 이러한 모델들은 두 가지 결정적인 단계, 즉 충돌의 초기 순간(자유 스트리밍 단계 이후)과 동결 해제(freeze-out) 초평면에서의 비평형 보정(out-of-equilibrium corrections)을 간과하였다. 동결 해제 보정의 누락은 참 하드론 다중도(multiplicity)가 공간 확산 계수($D_s$)에 대해 비물리적인 의존성을 갖게 만들었으며, 이는 참 수송 계수의 체계적인 결정을 방해했다. 또한, 초기 분포는 종종 국소 열적 평형 상태에 있다고 가정되어, 초기 생성의 비평형적 특성을 무시하였다.

방법론
저자들은 초기 단계와 동결 해제 표면 모두에서 비평형 보정을 명시적으로 계산함으로써, 기존에 개발된 참 쿼크용 유체역학적 프레임워크를 확장한다.

1. 초기 단계 초기화:

   • 이 모델은 참 생성 시간과 유체역학적 초기화 시간($\tau_0$) 사이의 자유 스트리밍(free-streaming) 단계를 도입한다.
   • 저자들은 특성 곡선법(method of characteristics)을 사용하여, 이 단계 동안 직선 궤적을 따라 이동하는 입자들에 대한 운동론적 방정식을 해석적으로 푼다.
   • 이 접근 방식은 초기 참 밀도를 수정하고, 이진 충돌 밀도($n_{coll}$)의 공간적 의존성과 초기 강한 산란의 운동량 분포로부터 유도된 비제로(non-zero) 확산 전류($\nu^\mu_q$)를 $\tau_0$에서 생성한다.
   • 란다우 매칭 조건(Landau matching conditions)이 사용되어, 미분 생성 단면적에 관해 전체 참 분포 함수의 적분을 고정한다.

2. 동결 해제 보정 (다종 접근법):

   • 저자들은 참의 하드론 공명 가스(HRGC) 내에서 다종 접근법을 채택하며, 시스템을 모든 참 하드론 종들이 경량 QGP 성분들과 상호작용하는 유체로 취급한다.
   • 각 하드론 종 $i$에 대한 분포 함수의 비평형 보정($\delta f_{i,p}$)은 모멘트 방법(method of moments)과 나비에-스토크스(Navier-Stokes) 근사를 사용하여 유도된다.
   • 보정은 기약 모멘트(irreducible moments)로 전개되며, 참 하드론의 벌크 압력 및 전단 응력 기여는 경량 성분에 비해 무시할 수 있다고 가정하여 랭크 $l=1$(확산 전류)에서 절단된다.
   • 핵심적인 해석적 결과(식 23)는 보정을 다음과 같이 표현한다:
     $$\delta f_{i,p} = f^{(0)}_{i,p} \frac{1}{P_i} \bar{\kappa}_i \nu^\mu_q p_{\langle \mu \rangle}$$
     여기서 $\bar{\kappa}_i$는 종의 전하와 밀도에 의존하는 정규화된 확산 계수이고, $\nu^\mu_q$는 전체 참 확산 전류이다. 이 공식화는 개별 종의 전류 대신 전체 전류만을 진화시킬 수 있게 한다.

3. 진화 및 스펙트럼:

   • 시스템은 일정한 $D_s T$를 갖는 1+1D 유체역학 프레임워크(FluiduM)를 사용하여 진화된다.
   • 불변 수율(invariant yields)과 횡방향 운동량($p_T$) 스펙트럼은 유도된 비평형 보정을 포함하는 쿠퍼-프라이(Cooper-Frye) 절차를 통해 계산된다.
   • 모델은 FONLL pQCD 예측과 ALICE 측정 참 단면적을 모두 입력값으로 사용하여, ALICE 실험 데이터(Pb-Pb 충돌, $\sqrt{s_{NN}} = 5.02$ TeV, 0–10% 중심도)와 비교 검증된다.

주요 기여 및 결과

• 비물리적 의존성 해결: 동결 해제 시의 비평형 보정을 포함함으로써, 참 하드론 다중도가 공간 확산 계수 $D_s$에 대해 보이는 비물리적인 의존성을 해결하였다. 이러한 보정이 없으면 다중도는 $D_s$에 따라 크게 변하지만, 보정을 적용하면 통합 수율은 테스트된 전체 $D_s$ 범위에서 비확산 사례 및 초기 참 쿼크의 총수와 일치하게 유지된다.
• 초기 확산 전류: 자유 스트리밍 초기화는 $\tau_0$에서 후기 시간에 생성되는 전류보다 한 자릿수 더 큰 확산 전류를 생성한다. 그러나 본 연구는 이 초기 조건이 제로 초기 전류(zero-initial-current) 시나리오와 비교했을 때 후속 참 장(field)의 진화와 최종 불변 수율에 미치는 영향이 미미하다는 것을 발견하였다.
• 모델-데이터 비교:
  • FONLL 단면적을 사용할 경우, 모델은 불확실성 범위 내에서 참 메존의 통합 수율을 재현한다.
  • $\Lambda_c^+$ 바리온 수율에서 텐션(tension)이 관찰되는데, 이는 과소 예측되었다. 기존 문헌을 따라, 저자들은 들뜬 참 바리온의 통계적 가중치를 3배로 높이는 "참-바리온 강화(charm-baryon enhancement)"를 적용한다. 이 조정을 적용하고 ALICE 측정 단면적을 사용하면, 모델은 $p_T \approx 3$ GeV/c까지까지 실험적 수율 및 $p_T$ 스펙트럼과 호환성을 달성한다.
• 적용 한계: 논문은 분포 함수의 양의 성질(positivity)에 기반한 유체역학적 접근법의 근본적인 한계를 식별한다.
  • 비평형 보정 $\delta f_{i,p}$는 (확산 전류의 음의 반경 방향 프로파일로 인해) 음수가 될 수 있다.
  • $|\delta f_{i,p}| \approx f^{(0)}_{i,p}$가 되면, 전체 분포 함수가 음수가 되어 물리적 해석이 불가능해진다.
  • 이러한 $p_T$ 스펙트럼의 "제로 크로싱(zero-crossing)"은 $D_s T$가 증가함에 따라 더 낮은 $p_T$ 값에서 발생한다. $D_s T \gtrsim 0.2$인 경우, 모델은 $p_T \sim 2$ GeV/c에서 음의 스펙트럼을 생성한다.
  • 따라서, 이 프레임워크는 전체 $p_T$ 범위에 대해 $D_s T \lesssim 0.1$인 경우에만 신뢰할 수 있으며, 높은 $D_s T$의 경우 제한된 저-$p_T$ 영역(예: $D_s T \approx 0.24$일 때 $p_T \lesssim 2.5$ GeV/c)에서만 유효하다.

의의
본 논문은 충돌의 초기 및 최종 단계 모두에서 비평형 효과를 고려하는 일관된 참 하드론 유체역학 프레임워크를 구축한다. 수송 계수에 대한 비물리적인 수율 의존성을 교정함으로써, 이 연구는 참 확산 계수의 체계적인 추출(예: 베이지안 분석을 통한)을 위한 필수적인 토대를 마련한다. 그러나 저자들은 현재의 선형화된 접근법이 높은 확산 계수 및 횡방향 운동량에서 섭동 가정($|\delta f| \ll f^{(0)}$)이 붕괴되는 문제로 인해 제한적이라고 겸허히 결론짓는다. 그들은 향-후 연구가 전체 운동학적 범위에서 양의 분포 함수를 보장하기 위해 최대 엔트로피 접근법(maximum entropy approach)과 같은 대안적인 방법을 탐구해야 한다고 제안한다.