Data-driven measurements of the beauty hadron production in p-Pb collisions at \sqrt{s_{\mathrm{NN}}} = 5.02 and 8.16 TeV

기술 요약: $\sqrt{s_{NN}} = 5.02$ 및 $8.16$ TeV p–Pb 충돌에서의 뷰티 헤드론 생성에 대한 데이터 기반 측정

문제 및 동기
중질량 쿼크를 포함하는 헤드론의 생성은 양자 색역학(QCD)을 테스트하는 중요한 시험대 역할을 한다. 짐(charm)과 뷰티(beauty) 생성의 경우 쿼크의 질량이 크기 때문에 섭동적 QCD(pQCD) 계산이 신뢰할 수 있지만, 중이온 충돌에서 쿼크-글루온 플라즈마(QGP)의 효과를 분리해내기 위해서는 냉핵 물질(Cold Nuclear Matter, CNM) 효과를 정밀하게 정량화하는 것이 필요하다. 탈국소화(deconfinement)가 없는 핵 환경에서 발생하는 이러한 효과는 핵자 분포 함수(nPDF)와 운동학적 분포를 수정한다.

뷰티 헤드론은 큰 질량 덕분에 섭동적 계산이 용이하며, 짐에 비해 핵 변조(nuclear modification)에서 뚜렷한 질량 계층 구조를 형성하기 때문에 CNM 효과의 귀중한 프로브가 된다. 그러나 양성자-핵(p–Pb) 충돌에서의 뷰티 생성에 관한 기존 측정들은 붕괴 생성물(비즉각적 $D^0$ 및 $J/\psi$ 메존 또는 붕괴 전자)의 재구성에 의존해 왔다. 이러한 방법들은 모(parent)가 되는 뷰리 헤드론의 운동학에 직접 접근하는 것을 제한하여, 전체 운동학적 위상 공간에 걸쳐 nPDF를 제약하는 능력을 제한한다. 본 연구는 넓은 라피디티(rapidity) 범위에 걸쳐 포괄적인 뷰티 헤드론 생성 단면적 및 핵 변조 인자($R_{p\text{Pb}}$)를 직접적이고 데이터 기반으로 결정해야 할 필요성을 다룬다.

방법론
본 연구는 측정된 비즉각적 딸(daughter) 스펙트럼으로부터 모 뷰티 헤드론의 운동학적 분포를 추론하기 위해 데이터 기반 베이지안 언폴딩(unfolding) 기법을 사용한다. 이 분석은 ALICE 및 LHCb 협력단의 입력을 활용한다:

• 입력값: 중간 라피디티($-0.96 < y < 0.04$)에서의 비즉각적 $D^0$ 메존 스펙트럼과, $\sqrt{s_{NN}} = 5.02$ 및 $8.16$TeV p–Pb 충돌에서 중간 라피디티, 전방($1.5 < y < 4.0$), 그리고 후방($-5.0 < y < -2.5$) 라피디티를 커버하는 비즉각적$J/\psi$ 메존 스펙트럼을 사용한다.
• 언폴딩 프레임워크: 실제 모 뷰티 헤드론이 재구성된 딸 입자를 생성할 확률을 인코딩하는 응답 행렬(response matrix)을 EvtGen 2.0 시뮬레이션으로부터 도출하는 반복적 베이지안 절차(D'Agostini 방식 준수)를 적용한다.
• 분기비(Branching Fractions): 딸 입수의 모 입자 단면적으로의 변환은 포괄적 분기비에 의존한다. $J/\psi$의 경우 PDG 값을 사용한다. $D^0$의 경우, 특정 붕괴 모드에 대한 PDG 값과 LHCb에서 pp 충돌에서 측정한 $b$-헤드론 생성 분율($f_u, f_d$)을 결합하여 포괄적 확률을 계산한다.
• 참조 계산: $\sqrt{s} = 8.16$ TeV에서의 pp 참조 단면적이 존재하지 않으므로, $\sqrt{s} = 5.02$ 및 $13$ TeV에서 언폴딩된 pp 결과를 FONLL(Fixed-Order Next-to-Leading Log) 예측의 비율을 스케일링 인자로 사용하여 외삽함으로써 pp 참조 단면적을 얻는다.
• 불확실성 전파: 계통 불확실성은 $p_T$ 빈(bin) 전반에 걸쳐 완전히 상관된 것으로 취급된다. 소스에는 실험적 입력 불확실성, 분기비 불확실성, 응답 행렬 통계, 반복 의존성, 그리고 저에너지 $p_T$ 외삽이 포함된다.

주요 기여 및 결과
본 논문은 p–Pb 충돌에서 전체 라피디티 범위($-4.0 < y < 4.0$)를 아우르는 최초의 언폴딩 기반 뷰티 헤드론 생성 측정치를 제시한다.

1. 단면적: 본 연구는 개방형 뷰티 생성에 대한 $p_T$ 및 $y$ 미분 단면적을 보고한다. 서로 다른 붕괴 채널($D^0$ 및 $J/\psi$)과 서로 다른 협력단(ALICE 및 LHCb)으로부터 얻은 언폴딩 결과는 공통 운동학 구간으로 매핑했을 때 상호 일관적이다.
2. 핵 변조 인자 ($R_{p\text{Pb}}$): $R_{p\text{Pb}}$는 횡운동량 및 라피디티의 함수로 보고된다. 결과는 기존 실험적 측정치 및 nPDF 세트(EPPS16, nCTEQ15, EPPS16*)에 기반한 이론적 계산과 호환되는 정밀도로 전체 라피디티 범위를 커버한다.
   • 언폴딩된 $R_{p\text{Pb}}$ 값은 기존의 실험적 측정치(비즉각적 $D^0$, $J/\psi$, $B$ 메존) 및 nPDF 세트에 기반한 이론적 계산과 일치한다.
   • 이 분석은 미래의 전역 nPDF 피팅을 위해 보완적인 정보를 제공하는, 데이터로 제약된 $R_{p\text{Pb}}$ 결정을 제공한다.
3. 질량 계층 구조: 중간 라피디티에서 뷰티와 짐 헤드론을 비교하면 뷰티 헤드론에서 체계적으로 더 큰 $R_{p\text{Pb}}$가 나타난다. 이는 짐에 비해 뷰티에 대해 더 약한 CNM 억제가 있음을 나타내며, 이는 냉핵 물질 효과의 예상된 질량 계층 구조와 일치한다.
4. 총 단면적: p–Pb 충돌에서의 총 포괄적 $b\bar{b}$ 생성 단면적을 보고한다:
   • $\sigma_{b\bar{b}}^{\sqrt{s_{NN}}=5.02\text{ TeV}} = 38.33 \pm 1.33 \text{ (stat.)} \pm 2.40 \text{ (sys.) mb}$
   • $\sigma_{b\bar{b}}^{\sqrt{s_{NN}}=8.16\text{ TeV}} = 64.60 \pm 1.81 \text{ (stat.)} \pm 5.87 \text{ (sys.) mb}$
     이 결과들은 이전의 위상 공간 외삽 결과들과 일치한다.
5. 스펙트럼 형태: 평균 제곱 횡운동량 $\langle p_T^2 \rangle$은 두 에너지 모두에서 p–Pb 및 pp 충격 간에 일관되게 나타나며, 이는 현재의 불확실성 범위 내에서 p–Pb 충돌에서의 $p_T$ 스펙트럼 형태에 유의미한 변형이 없음을 나타낸다.

의의
본 논문은 베이지안 언폴딩을 사용하여 반포괄적(semi-inclusive) 붕괴 채널로부터 포괄적 뷰티 헤드론 운동학을 추출하는 견고한 방법론을 확립한다. 높은 정밀도로 전체 라피디티 범위를 커버하는 최초의 $R_{p\text{Pb}}$ 결과를 제공함으로써, 이 연구는 특히 낮은 $x$ 영역에서의 핵자 글루온 분포를 포함한 핵자 분포 함수(nPDF)를 직접적으로 제약한다. 뷰티가 짐보다 약한 억제를 보인다는 관측은 질량 의존적 CNM 효과에 대한 이론적 이해를 강화한다. 이러한 결과는 중이온 충돌에서 QGP의 신호를 분리하기 위한 필수적인 벤치마크 역할을 하며, 향후 전역 nPDF 피팅을 위한 중요한 입력을 제공한다.