Comparision of modeling of prompt photon production in proton-proton collisions at energies NICA $\sqrt{s}$ =10 GeV in PYTHIA and calculations using FeynCalc

기술 요약: NICA 에너지 영역의 양성자-양성자 충돌에서의 프롬프트 광자 생성 모델링 비교

문제 제기
본 연구는 NICA(Nuclotron-based Ion Collider fAcility) 및 SPD(Spin Physics Detector) 실험과 관련된 에너지 영역인 중심 질량 에너지 $\sqrt{s} = 10$ GeV에서의 양성자-양성자($pp$) 충돌 내 프롬프트 광자 생성의 모델링 및 이론적 계산을 다룬다. 프롬프트 광자는 파톤 분포 함수(PDF)와 강한 산란 메커니즘을 이해하는 데 중요한 탐침 역할을 하지만, 특정 부과정(subprocess)의 기여와 초기 상태 편극(initial-state polarization)의 효과와 관련하여 미분 단면적의 차이를 이해하기 위해 현대적인 몬테카를로 이벤트 생성기와 엄밀한 해석적 계산을 비교할 필요가 있다.

방법론
저자들은 프롬프트 광자 생성의 세 가지 주요 부과정을 모델링하기 위해 이중 방법론을 채택하였다:

1. 몬테카를로 모델링 (PYTHIA 8.316): 본 연구는 초기 입자 상태, 강한 산란 행렬 요소, 후속 파톤 샤워(초기 및 최종 상태 복사 - ISR/FSR), 다중 파톤 상호작용(MPI), 그리고 룬드 스트링 모델(Lund string model)을 통한 강입자화(hadronization)를 포함한 전체 이벤트 체인을 시뮬레이션하기 위해 PYTHIA 8.316을 활용하였다. 이 접근 방식은 비섭동적 효과와 소프트/콜리니어 복사를 포착한다.
2. 해석적 계산 (FeynCalc): FeynCalc를 기반으로 한 병렬 계산이 페인만 다이어그램법을 사용하여 수행되었다. 이 계산은 다음의 특정 부과정에 대한 최고 차수(Leading Order, LO) 행렬 요소에 집중하였다:
   • 쿼크-글루온 컴프턴 산란 ($qg \to q\gamma$).
   • 쿼크-반쿼크 쌍 소멸 ($q\bar{q} \to g\gamma$).
   • 쿼크의 브렘스트랄룽(Bremsstrahlung, 제동 복사) ($q \to q\gamma$).
     해석적 접근 방식은 LHAPDF 라이브러리(버전 6.5.5)를 통해 파톤 분포 함수(PDF)를 통합하였으며, 편극된 PDF와 헬리시티 의존적 행렬 요소를 사용하여 종방향 편극된 양성자로 확장되었다.

연구진은 통계적 분포(광자 에너지, 라피디티, 불변 질량)와 충돌 에너지($\sqrt{s}$), 횡운동량($p_T$), 산란각 코사인($\cos\theta$), 그리고 라피디티($y$)에 따른 미분 단면적을 분석하였다.

주요 기여 및 결과

• 부과정 지배력 (해석적 계산): 특정 해석적 프레임워크(FeynCalc) 내에서 부과정이 격리되고 파톤 캐스케이드가 제외된 경우, **컴프턴 산란($qg \to q\gamma$)**이 가장 지배적인 메커니즘으로 나타났으며, 전체 단면적의 **51.27%**를 차지하였다. 쿼크-반쿼크 소멸($q\bar{q} \to g\gamma$)은 **47.22%**를 기여하였고, 쿼크 브렘스트랄룽은 **0.003%**로 무시할 만한 수준이었다. 논문은 PYTHIA 8.316이 (캐스케이드와 강입자화를 포함한) 전체 광자 생성 메커니즘을 시뮬레이션하므로, 이 특정 백분율 분해는 전체 PYTHIA 이벤트 시뮬레이션이 아닌 격리된 파톤 수준의 계산에 적용된다는 점을 명시한다.
• 통계적 분포:
  • 생성된 광자의 에너지($E_\gamma$) 분포는 고에너지 쪽으로 비대칭적으로 평탄해지는 특성을 가진 **로그-노멀 분포(log-normal distribution)**를 따르는 것으로 밝혀졌다.
  • 라피디티($y$)에 따른 분포는 중앙 라피디티를 중심으로 대칭인 **정규 분포(Gaussian distribution)**를 따랐으며, 이는 대부분의 강한 광자가 상호작용하는 파톤들이 동일한 운동량 분율을 가질 때 생성됨을 나타낸다.
• PYTHIA와 FeynCalc의 비교:
  • 질적 유사성: 두 방법 모두 $p_T$, $\cos\theta$, $y$에 대한 미분 단면적의 함수적 의존성에 대해 유사한 결과를 생성하였다. 예를 들어, 두 방법 모두 $p_T$가 증가함에 따라 단면적이 급격히 감소하고, 각도 및 라피디티 의존성에서 "U자형" 분포(전방/후방 방향에서 피크)를 보였다.
  • 양적 차이: 상당한 수치적 차이가 관찰되었다. PYTH레이야 결과는 FeynCalc의 "순수(bare)" LO 계산에 비해 일반적으로 더 매끄러운 분포와 다른 정규화(normalization)를 보여주었다.
  • 차이의 원인: 저자들은 이러한 차이를 PYTHIA에 포함된 파톤 샤워, 고유 횡운동량($k_T$), 그리고 강입자화 때문이라고 설명하였다. 구체적으로, FeynCalc 계산은 낮은 $p_T$ 및 위상 공간 경계에서 조절되지 않은 소프트 및 콜리니어 발산(divergence)을 보인 반면, PYTHIA의 수다코프(Sudakov) 폼 팩터와 샤워링 메커니즘은 이러한 영역을 자연스럽게 억제하였다.
• 편극 효과:
  • 비편극 양성자의 경우, PYTHIA와 FeynCalc 사이의 적분 단면적 비율은 약 101-102%로 안정적으로 유지되었으며, 이는 운동학적 재분배에도 불구하고 총 정규화가 보존됨을 시사한다.
  • 종방향 편극 양성자의 경우, 비율이 1에서 크게 벗어났으며 에너지에 대한 강한 비선형 의존성을 보였다. 저자들은 이것이 표준 이벤트 생성기(PYTHIA 등)의 이론적 한계를 강조한다고 주장하는데, 이러한 생성기들은 종종 스핀 평균화되거나 근사된 분할 커널(splitting kernels)에 의존하여, 엄밀한 해석적 처리를 수행하는 FeynCalc에 비해 다중 방출을 통한 정확한 양자 헬리시티 상태를 완전히 보존하지 못한다는 점을 지적한다.

의의 및 주장
본 논문은 이 비교가 중간 에너지 영역(NICA)에서의 프롬프트 광자 생성 메커니즘을 명확히 하는 데 필요한 자료를 제공한다고 주장한다.

• 상호 보완적 접근: 저자들은 해석적 계산(FeynCalc)이 특정 운동학적 변수와 스핀 구조에 대한 투명성과 제어를 제공하여 특정 부과정의 기여를 격리할 수 있다고 상정한다(예: 격리된 계산에서 나타난 51.27%의 컴프턴 지배력). 반대로, PYTHIA는 실험적 예측에 필수적인 비섭동적 효과를 포함한 전체 이벤트 환경에 대한 현실적인 시뮬레이션을 제공한다.
• 도구의 검증: 본 연구는 현대적인 도구인 PYTHIA가 복잡한 QCD 효과를 모델링함으로써 실험 예측의 정확도를 높이지만, 방사 및 강입자화 처리 방식의 차이로 인해 순수 이론 계산과는 수치적으로 다를 수 있음을 입증한다.
• 스핀 물리학: 이 작업은 스핀 의존 프로세스의 모델링의 복잡성을 강조하며, 편극된 충돌에 대한 해석적 결과와 생성기 기반 결과 사이의 불일치가 파톤 샤워를 통해 정확한 스핀 밀도 행렬을 추적하는 것의 어려움을 드러낸다고 언급한다.

저자들은 두 방법의 사용이 필수적이라고 결론짓는다. 즉, 섭동 QCD 및 PDF 민감도를 테스트하기 위해서는 FeynCalc가 필요하며, 비섭동적 효과가 중요한 데이터를 분석하고 탐지 기술을 개발하기 위해서는 PYTHIA가 필요하다.