[논문 리뷰] Automation of 2-loop Amplitude Calculations

Stephen Jones|arXiv (Cornell University)|2016. 08. 12.

Particle physics theoretical and experimental studies인용 수 3

한 줄 요약

이 논문은 양자장이론에서 두 루프 진폭을 자동으로 계산하는 프레임워크를 제시한다. 특히 QCD에서 NLO 수준의 힉스 보손 쌍 생성을 중심으로 다룬다. 섹터 분해와 준 몬테카를로(QMC) 적분을 조합하여 마스터 적분을 수치적으로 평가함으로써, GPU 가속화 및 위상공간 샘플링 최적화된 워크플로우를 사용해 $gg \to HH$의 미분 단면적에 대한 첫 번째 완전 수치적 계산을 가능하게 하였다.

ABSTRACT

Some of the tools and techniques that have recently been used to compute Higgs boson pair production at NLO in QCD are discussed. The calculation relies on the use of integral reduction, to reduce the number of integrals which must be computed, and expressing the amplitude in terms of a quasi-finite basis, which simplifies their numeric evaluation. Emphasis is placed on sector decomposition and Quasi-Monte Carlo (QMC) integration which are used to numerically compute the master integrals.

연구 동기 및 목표

• QCD에서 NLO 초과 정밀도 계산을 위한 자동화되고 수치적으로 안정적인 프레임워크 개발.
• 해석적 결과가 제공되지 않을 경우 다중 스케일을 가진 2루프 진폭을 계산하는 데 도전.
• 정확한 단면적 평가를 위해 필요한 위상공간 점 수를 최소화하기 위한 최적화된 샘플링 기법 적용.
• 섹터 분해와 QMC 적분을 통한 마스터 적분의 효율적 수치 평가 보장.
• GPU 가속화와 다중 스레딩을 활용한 가상 진폭의 고성능 계산 달성.

제안 방법

• 우주론적 및 적외선 발산을 다루기 위해 $D = 4 - 2\epsilon$를 사용하는 전통적인 차원 정규화(CDR)를 적용.
• 루프 적분을 분리하기 위해 형식 인자 $F_1, F_2$를 프로젝터 $P_1^{\mu\nu}, P_2^{\mu\nu}$와 수축.
• REDUZE를 통한 적분 감소를 적용하여 독립적 적분 수를 감소시키나, 비플레인어 4점 적분은 별도의 처리 필요.
• 해석적 감소가 불가능한 적분에 대해 섹터 분해를 적용하여 특이성을 분리하고 수치 적분을 위한 준비.
• R1SL 및 VEGAS 알고리즘을 사용한 준 몬테카를로(QMC) 적분을 구현하여 분해된 적분을 고정밀도로 평가.
• OpenCL을 활용해 CPU 및 GPU를 통한 수치 적분을 병렬화하여 고처리량 및 코어 수에 비례하는 선형 스케일링 달성.

실험 결과

연구 질문

• RQ1 해석적 해가 없는 다중 스케일을 가진 2루프 진폭을 수치적으로 어떻게 계산할 수 있는가?
• RQ2 고차원 위상공간에서 섹터 분해된 준유한 적분에 QMC 적분을 적용했을 때 성능과 정밀도는 어떠한가?
• RQ3 정밀도를 훼손하지 않고도 위상공간 샘플링을 얼마나 최적화할 수 있는가?
• RQ4 수치적으로 집약적인 2루프 진폭 계산에 대해 GPU 가속화 및 다중 스레딩 구현이 얼마나 효과적으로 확장되는가?
• RQ5 일반적인 2루프 과정에 대해 적분 감소, 섹터 분해, 수치 적분을 통합한 완전 자동화 프레임워크를 구축할 수 있는가?

주요 결과

• 이 프레임워크를 통해 QCD에서 NLO 수준의 힉스 보손 쌍 생성($gg \to HH$)에 대한 미분 단면적의 첫 번째 수치 계산을 달성하였다.
• 사건 샘플링으로 인한 통계적 불확도 0.3%와 마스터 적분의 수치 적분으로 인한 추가 불확도 0.1%로 총 단면적을 계산하였다.
• 위상공간 점당 중앙값 GPU 소요 시간은 2시간이며, 벽시계 시간 6일 동안 총 4680 GPU 시간을 소비하였다.
• 물리적 코어 수 한계까지 CPU 코어 수에 비례하여 선형 확장성을 확보하였으며, 단일 Nvidia Tesla K20Xm GPU를 사용해 상당한 속도 향상을 달성하였다.
• VEGAS 및 R1SL QMC 알고리즘 모두 양호한 성능을 보였으며, R1SL은 더 큰 함수 차원에서 더 우수한 확장성을 보였고, VEGAS는 거의 최적의 몬테카를로 스케일링을 제공하였다.
• 위상공간 격자 없이 666개의 위상공간 점을 성공적으로 계산하여, 다양한 컷, PDF 또는 질량 중심계 에너지에 대해 재평가가 가능하게 하였다.