Automated NRQCD and NRQED simulations of quarkonium and leptonium production with P-wave states and physical-mass effects

기술 요약: P-파동 상태 및 물리적 질량 효과를 포함한 쿼크로니움 및 레프토늄 생성의 자동화된 NRQCD 및 NRQED 시뮬레이션

문제 정의
무거운 쿼크 결합 상태(heavy quarkonia)와 레프토니아(leptonia, 레프톤의 결합 상태)의 생성 및 붕괴는 비상대론적 영역에서 양자 색역학(QCD)과 양자 전기 역학(QED)을 테스트하기 위한 중요한 실험실 역할을 한다. 비상대론적 QCD(NRQCD) 및 비상대론적 QED(NRQED) 프레임워크는 단거리 계수(short-distance coefficients)와 장거리 행렬 요소(long-distance matrix elements, LDMEs)의 엄격한 인자 분해(factorization)를 제공하지만, 이벤트 생성의 자동화는 역사적으로 제한적이었다. MadOnia 및 HELAC-Onia와 같은 기존 도구들은 주로 S-파동 상태에 집중하거나 임의의 수의 P-파동 상태에 대한 완전한 자동화가 부족했다. 또한, 기존의 섭동 계산은 종종 결합 상태의 질량을 구성 입자 질량의 합($M_B \approx m_Q + m_{\bar{Q}}$)으로 근사한다. 이러한 근사는 특히 운동학적 임계 영역과 동일한 구성 성분을 가지나 물리적 질량이 서로 다른 여러 결합 상태가 포함된 과정(예: $J/\psi + \psi(2S)$)에서 체계적인 불확실성을 초래한다.

방법론
저자들은 MadGraph5_aMC@NLO(MG5_aMC) 프레임워크의 새로운 모듈인 MadSONS(MadGraph Simulations Of NRQFT States)를 제시한다. 이 모듈은 기존에 확립된 S-파동 구현을 확장하여, P-파동 상태를 포함한 임의의 트리 레벨 생성 과정과 물리적 질량 효과를 포함한다.

1. 이중 수(Dual Numbers)를 통한 P-파동 구현:
   S-파동 상태와 달리, P-파동 상태는 원점에서 파동함수가 0이 되므로, 그 생성율은 파동함수의 1차 미분에 의존한다. 운동량 공간에서 이는 상대 운동량에 대한 미분을 포함하는 궤도 각운동량 투영기(projectors)에 해당한다.

• 저자들은 다변수 이중 수(multivariate dual numbers), 즉 자동 미분의 한 형태를 사용하여 이러한 미분을 구현한다. 이 접근 방식은 기호 미분(symbolic differentiation)이나 유한 차분(finite-difference) 근사에 의존하지 않고 정확하고 수치적으로 안정적인 미분을 제공한다.
• 이중 수 형식은 진폭 생성 체인을 통해 미분을 전파할 수 있게 하여, 행렬 요소 계산 내에서 궤도($P_L$) 및 총 각운동량($P_J$) 투영기를 직접 평가할 수 있게 한다.
• 이 방법은 향후 고차 궤도 각운동량(예: D-파동) 및 고차 미분이 필요한 차세대 격차(NLO) 계산을 위한 확장에 용이하도록 설계되었다.

2. 물리적 질량을 위한 운동량 재배열(Momentum Reshuffling):
   물리적 결합 상태 질량과 구성 입자 질량의 합 사이의 불일치를 해결하기 위해, 저자들은 운동량 재배열 절차를 도입한다.

• 위상 공간 생성: 물리적 질량을 사용하여 이벤트를 생성함으로써, 올바른 가용 위상 공간 부피와 운동학적 임계치를 보장한다.
• 행렬 요소 평가: NRQFT 형식에서 요구되는 단거리 계수를 위해, 구성 입자의 온-쉘(on-shell) 조건(구성 입자 질량의 합)을 사용하여 행렬 요소를 평가한다.
• 재배열 알고리즘: 물리적 운동량과 구성 입자 수준의 운동량 사이의 매핑은 두 가지 구별된 반동(recoil) 체계를 통해 이루어진다:
  • 초기 상태 재배열(Initial-state reshuffling): 최종 상태의 운동량을 고정시킨 채로 입사 파톤의 운동량을 수정하여 4-운동량을 보존한다.
  • 최종 상태 재배열(Final-state reshuffling): 물리적 질량에 대한 온-쉘 조건을 만족시키기 위해 출력 운동량(및 잠재적으로 파톤 중심-질량 에너지)을 수정한다.
• 이 전략은 위상 공간 측도(phase-space measure)에 고차 상대론적 보정을 포함시키는 것을 근사하며, 임계 영역 근처에서의 관측량 묘사를 개선한다.

3. UFO 모델 및 자동화:
   이 구현은 업데이트된 UFO (Universal Feynman Output) 모델인 sm_onia에 캡슐화되어 있다. 이 모델은 Fock 상태(P-파동 구성인 $^3P_J$ 및 $^1P_1$ 포함)를 정의하며, 사용자가 입력 카드(onia_card.dat)를 통해 물리적 질량과 LDME를 지정할 수 있도록 한다. 이 프레임워크는 컬러 싱글렛(color-singlet) 및 컬러 옥텟(color-octet) 구성을 모두 포함하는 쿼크로니아(NRQCD)와 레프토니아(NRQED)를 지원한다.

주요 기여 및 결과

• 검증: P-파동 구현은 HELAC-Onia 생성기를 통해 검증되었다.

  • 제곱 진폭(Squared Amplitudes): 개별 위상 공간 지점에서의 비교 결과, 이중 정밀도 산술(double-precision arithmetic) 수준($10^{-14}$에서 $10^{-16}$ 사이의 상대적 편차)에서 일치함을 보였다.
  • 적분 단면적: 단일 및 이중 쿼크로니아 생성에 대한 벤치마크는 $pp$및$e^+e^-$ 충돌에서 몬테카를로 통계적 불확실성 범위 내에서 일치함을 보여준다.
  • 레프토니아: NRQED의 경우, 코드는 광범위한 충돌 에너지 영역에서 포지트로늄(positronium) 및 디타우오니움(ditauonium)의 P-파동 생성에 대한 해석적 단면적 계산을 서브 퍼밀(sub-permille) 수준으로 재현한다.

• 현상론적 응용:

  • 단일 및 이중 쿼크로니아 생성: 저자들은 LHC($\sqrt{s}=13$ TeV) 및 B-팩토리($\sqrt{s}=10.58$ GeV)에서의 다양한 과정에 대한 단면적을 제시한다. 이들은 운동량 재배열의 영향을 정량화하였으며, 단일 생성의 경우 효과가 미미하지만, 큰 질량 차이를 가진 상태들이 포함된 이중 생성 채널(예: $J/\psi + \psi(2S)$)에서는 효과가 상당할 수 있음(최대 ~20-50%)을 언급하였다.
  • $J/\psi + \psi(2S)$ 사례 연구: LHC에서의 연관 $J/\psi + \psi(2S)$ 생성에 대한 상세 연구는 본 프레임워크의 유용성을 입증한다. 물리적 질량과 P-파동 피드다운(feed-down) 기여($\psi(2S)$ 및 $\chi_{cJ}$로부터)를 포함함으로써, 예측된 단면적은 전통적인 처리 방식에 비해 감소한다. 이러한 감소는 LHCb 데이터에 대한 개선된 설명을 제공하며, 잠재적인 이중 파톤 산란(DPS) 기여를 위한 여지를 더 많이 남겨준다. 이 연구는 P-파동 기여가 직접 생성의 약 17%, 전체 파이딜(fiducial) 단면적의 25%를 차지한다는 점을 강조한다.
  • 레프토니아 생성: 이 프레임워크는 포지트로늄 및 디타우오니움과 같은 이색 원자(exotic atoms)의 연구를 가능하게 한다. 저자들은 닫힌 형태(closed-form)의 NRQED 예측을 도출하고, 자동화된 시뮬레이션이 해석적 결과와 일치함을 확인하여 QED 영역에서의 P-파동 투영기를 검증하였다.

의의 및 주장
본 논문은 MadSONS가 결합 상태 생성의 완전한 자동화를 향한 중요한 단계라고 주장한다. LO S- 및 P-파동 과정을 자동화함으로써, 이 도구는 NRQCD에서 $v^7$ 차수까지 현상론적으로 유의미한 모든 Fock 상태를 체계적으로 포함할 수 있게 한다.

저자들은 이중 수가 향-차(NLO) 확장에 필수적인 높은 수치 정밀도를 요구하는 적외선 발산(infrared divergence)의 상쇄를 위해 견고한 토대를 제공한다는 점을 강조한다. 또한, 운동량 재배열 절차는 행렬 요소에 대한 완전한 상대론적 보정을 요구하지 않으면서도 물리적 질량 효과를 포함할 수 있는 실용적인 방법을 제공하여, 운동학적 임계 영역 근처에서의 시뮬레이션 신뢰도를 높인다.

이 작업은 MadSONS를 쿼크로니아 생성 메커니즘을 조사하고, NRQCD 인자 분해를 테스트하며, 이색 레프토니아 시스템을 탐구하기 위한 다재다능한 도구로 위치시킨다. 동시에 Pythia8와 같은 표준 파톤 샤워 및 해드로니제이션 도구와 호환성을 유지한다. 저자들은 LO 예측이 데이터를 개선하지만, 고차 QCD 보정, 상대론적 효과, 그리고 잠재적인 DPS 기여로 인한 불일치가 남아 있을 수 있으며, 이는 향후 조사의 대상임을 겸허히 언급하였다.