네 개 탑 쿼크 생산 3렙톤 채널에서 NLO QCD 보정 완전 분석

초록

본 연구는 LHC에서 네 개의 탑 쿼크가 동시에 생성되는 희귀 과정(pp→tt̄tt̄)을 3렙톤(세 전자·뮤온) 최종 상태로 고려하고, 생산과 붕괴 단계 모두에 대해 좁은 폭 근사(NWA)를 이용해 차세계(QCD) 차수(NLO) 보정을 수행하였다. QCD 보정이 미치는 영향을 조사하기 위해 두 가장 강한 라이트 제트의 불변 질량 차이(Mjj) 제한 파라미터 Qcut을 변동시키며, 전체 및 미분 단면에 대한 K‑factor와 스케일 불확실성을 상세히 분석하였다.

상세 분석

이 논문은 네 개 탑 쿼크 동시 생산 과정(pp→tt̄tt̄)을 3렙톤 채널에서 NLO QCD 정확도로 계산한 최초의 전면적 연구 중 하나이다. 저자들은 좁은 폭 근사(NWA)를 적용해 생산 단계와 붕괴 단계(각 탑→Wb) 를 명확히 분리하고, 세 가지 NWA 구현 방식을 제시한다. 첫 번째인 “NLO full”은 생산·붕괴 모두에 NLO 보정을 적용하고 탑 폭을 고정값(NLO 정확도)으로 사용한다. 두 번째인 “NLO exp”(또는 단순히 NLO)는 동일하게 보정을 적용하되, 탑 폭을 αs 전개식으로 확장해 일관성을 유지한다. 세 번째인 “NLO LO dec”은 생산 단계에만 NLO 보정을 적용하고, 붕괴는 LO 폭을 사용한다. 이러한 구분은 각각의 단계에서 발생하는 하드 방사와 스핀 상관관계가 어떻게 변하는지를 정량화하는 데 필수적이다.

계산은 HELAC‑NLO 프레임워크를 이용해 √s=13.6 TeV에서 수행되었으며, PDF는 MSHT20(N)LO, 스케일은 동적 선택 µ0=ET/4 (ET는 최종 상태 탑과 반탑의 transverse mass 합) 로 설정하였다. 사건 선택 기준은 최소 4개의 b‑제트, 최소 2개의 라이트 제트, 정확히 3개의 전자·뮤온, 그리고 누락된 에너지(ν)이다. NLO에서 추가 라이트 제트가 발생하면 기존 W→jj 제트가 재조합될 수 있어, Mjj가 W 질량과 크게 벗어나는 경우가 발생한다. 이러한 비물리적 위상을 억제하기 위해 |Mjj−mW|<Qcut 조건을 도입했으며, 기본값 Qcut=25 GeV 로 설정하고 5 GeV부터 100 GeV까지 변화를 조사하였다.

결과적으로, LO 단면은 Qcut에 무감각하지만 NLO 단면은 Qcut에 강하게 의존한다. Qcut를 100 GeV까지 확대하면 K‑factor가 1.45에 달하고, 스케일 불확실성도 +42%/‑30% 로 급증한다. 반면 Qcut=25 GeV에서는 σNLO=44.9 fb (스케일 +18%/‑23%) 로 비교적 안정적인 결과를 얻는다. 또한, NWA 구현 방식 간 차이는 Qcut이 작을 때(≤15 GeV) 최대 40%까지 벌어지며, 이는 스케일 불확실성을 초과한다. Qcut=25 GeV에서는 모든 구현이 10‑12% 이내 차이로 수렴한다.

미분 분포에서는 pT, j2와 ΔRj1j2가 특히 민감하게 변한다. LO에서는 pT, j2가 200 GeV에서 급격히 꺼지지만 NLO에서는 추가 방사로 인해 고pT 영역이 채워진다. ΔRj1j2는 LO에서 ≈1에 집중되지만 NLO에서는 뒤쪽(π) 피크가 나타나, 두 제트가 뒤로 튀는 이벤트가 증가함을 보여준다. 붕괴 단계에 NLO 보정을 포함하지 않은 “NLO LO dec”와 “NLO exp”를 비교하면, ΔRj1j2에서는 최대 22% 차이가 발생해 스케일 불확실성을 넘어선다. 반면 Δφℓ1ℓ2(두 가장 강한 전자·뮤온 사이의 방위각)에서는 차이가 10% 이하로, 스케일 밴드 안에 머문다. 전반적으로, 붕괴 단계까지 NLO 보정을 적용하면 전체 스케일 불확실성이 감소하고, 스핀 상관관계가 정확히 전달된다.

마지막으로 저자들은 향후 파트론 샤워와 매트릭스 엘리먼트 보정(MEC)을 포함한 NLO+PS 매칭을 수행해, 현재 고정된 고정계산과의 차이를 정량화하고, 실험 분석에 바로 적용 가능한 시뮬레이션 툴을 제공할 계획임을 밝힌다.