[논문 리뷰] New developments in PJFry

J. Fleischer, T. Riemann외 1명|DESY (CERN, DESY, Fermilab, IHEP, and SLAC)|2012. 10. 15.

Black Holes and Theoretical Physics참고 문헌 29인용 수 6

한 줄 요약

이 논문은 양자장론에서 수치적으로 안정적인 1-loop 텐서 적분 계산을 위한 오픈소스 C++ 패키지인 PJFry를 제시한다. 차원 재귀와 급수 전개를 사용하여 역행렬 그람 행렬식을 피함으로써, 질량이 임의로 설정된 5점 및 6점 함수까지 랭크 6까지의 텐서 적분을 안정적으로 평가할 수 있게 하며, 고에너지 물리학의 진폭 계산 및 자동화된 진폭 계산에 적용 가능하다.

ABSTRACT

We report on recent progress in numerical evaluation of one loop tensor integrals. A public C++ package PJFry implementing algorithms from arXiv:1009.4436 and its extension to hexagons up to rank 6 are presented.

연구 동기 및 목표

• 다양한 레그를 포함한 과정에서 역행렬 그람 행렬식으로 인한 1-loop 텐서 적분의 수치적 불안정성을 해결하기 위해.
• 4점 함수를 초월하는 고랭크 텐서 적분을 위한 강력하고 공개 가능한 도구를 개발하기 위해.
• 임의의 질량과 랭크를 가진 5점 및 6점 함수까지의 텐서 적분 축소 적용 범위를 확장하기 위해.
• 작은 그람 행렬식 영역에서 점근 전개를 사용하여 텐서 형상 인자에 대한 정확한 평가를 가능하게 하기 위해.
• 자동화된 진폭 계산 및 루프 진폭 라이브러리에서 사용 가능한 고성능, 재사용 가능한 프레임워크를 제공하기 위해.

제안 방법

• 패키지는 분모 거듭제곱이 이동된 고차원 스칼라 적분으로 텐서 적분을 표현하기 위해 차원 재귀 관계를 사용한다.
• 식 (6)과 (7)의 재귀 관계를 사용하여 텐서 형상 인자를 고차원에서의 스칼라 적분으로 감소시켜, 그람 행렬식의 직접적인 역행렬을 피한다.
• 작은 그람 행렬식 영역에서는 재귀 관계 (10)에서 유도된 급수 전개를 적용하여 수치적 안정성을 유지한다.
• 조합적 계수와 거듭제곱 상승/하향 연산자를 통해 텐서 형상 인자를 일반화된 스칼라 적분으로 매핑한다.
• 캐시 시스템을 통해 저차원 텐서 적분을 재사용하여 계산을 가속화하고 중복 평가를 줄인다.
• C, FORTRAN, C++, Mathematica, GoSam 등의 다양한 인터페이스를 지원하며, qd 라이브러리를 통한 복소수 질량의 고정밀도 연산도 제공한다.

실험 결과

연구 질문

• RQ1 작거나 0에 가까운 그람 행렬식이 존재하는 상황에서 1-loop 텐서 적분은 어떻게 안정적으로 평가할 수 있는가?
• RQ2 5점 및 6점 함수의 고랭크 텐서 적분을 평가할 때 가장 효율적이고 수치적으로 안정적인 방법은 무엇인가?
• RQ3 차원 재귀와 급수 전개를 조합하여 역행렬 그람 행렬식을 도입하지 않고도 작은 그람 영역에서 정확성을 유지할 수 있는가?
• RQ4 임의의 질량 구성 조건에서 랭크 6까지의 텐서 적분을 지원할 수 있는 일반적인 목적의 고성능 라이브러리를 어떻게 설계할 수 있는가?
• RQ5 텐서 적분 축소 과정에서 중간 결과를 캐싱하고 재사용하는 최적의 전략은 무엇인가?

주요 결과

• PJFry는 작은 그람 행렬식 영역이 존재하는 경우에도 5점 및 6점 함수의 랭크 6까지의 텐서 적분을 완전한 수치적 안정성으로 평가할 수 있다.
• 재귀 기반 감소와 점근 급수 전개를 통해 역행렬 그람 행렬식을 피함으로써, 모든 위상공간 영역에서 강력한 안정성을 확보한다.
• 낮은 차원의 텐서 적분을 재사용하는 캐싱 메커니즘 덕분에 고성능을 달성하여 추가 비용 없이 계산을 가속화한다.
• 6개의 외부 입자와 한 개의 질량 있는 내부 선을 가진 테스트 위상공간 점에서 스칼라 육각형 $F_0$는 $\epsilon=0$일 때 $-30.3579 - 205.8213i$로 평가되며, $\epsilon=0,1,2$에서 일관된 결과를 얻는다.
• 벡터 육각형 성분은 동일한 프레임워크를 사용하여 정확하게 계산되었으며, 복소 평면 상에서 기대되는 행동과 일치한다.
• Mathematica 인터페이스는 LoopTools에서 사용하는 표준 관례를 지원하며, 기호적 진폭 계산에서 직접 사용 가능하다.