[논문 리뷰] An efficient algorithm for the computation of the trace of the symmetrized product of an arbitrary number of Dirac matrices with two indices

Fernando Izaurieta, R. Ramirez-Gomez외 1명|arXiv (Cornell University)|2011. 06. 08.

Black Holes and Theoretical Physics인용 수 1

한 줄 요약

이 논문은 2n차원 시공간에서 n개의 딜라 클라크 행렬의 대칭화된 곱의 트레이스를 계산하는 매우 효율적인 알고리즘을 제시한다. 이는 카른-시몬스 중력 이론을 구성하는 데 핵심적인 단계이다. 트레이스를 n의 정수 분할에 대한 합으로 표현하고, 계산 가능한 계수를 사용하며, 재귀 관계를 활용함으로써, n ≤ 25일 경우 1분 이내에 결과를 산출할 수 있다. 이는 기존 도구들에 비해 극적으로 뛰어난 성능을 발휘한다.

ABSTRACT

A genuine gauge theory for the Poincare, de Sitter or anti-de Sitter algebras can be constructed in (2n-1)-dimensional spacetime by means of the Chern-Simons form, yielding a gravitational theory that differs from General Relativity but shares many of its properties, such as second order field equations for the metric. The particular form of the Lagrangian is determined by a rank n, symmetric tensor invariant under the relevant algebra. In practice, the calculation of this invariant tensor can be reduced to the computation of the trace of the symmetrized product of n Dirac Gamma matrices \Gamma_{ab} in 2n-dimensional spacetime. While straightforward in principle, this calculation can become extremely cumbersome in practice. For large enough n, existing computer algebra packages take an inordinate long time to produce the answer or plainly fail having used up all available memory. In this talk we show that the general formula for the trace of the symmetrized product of 2n Gamma matrices \Gamma_{ab} can be written as a certain sum over the integer partitions s of n, with every term being multiplied by a numerical coefficient \alpha_{s}. We then give a general algorithm that computes the \alpha-coefficients as the solution of a linear system of equations generated by evaluating the general formula for different sets of tensors B^{ab} with random numerical entries. A recurrence relation between different coefficients is shown to hold and is used in a second, algorithm to greatly speed up the computations. Runtime of the minimal algorithm stays below 1 min on a typical desktop computer for up to n=25, which easily covers all foreseeable applications of the trace formula.

연구 동기 및 목표

• 큰 n에 대해 딜라 클라크 행렬의 대칭화된 곱의 트레이스를 계산하는 데 발생하는 계산 블로킹 문제를 해결하기 위해.
• n > 10일 경우 실패하거나 실용적으로 너무 느려지는 기존 컴퓨터 대수 시스템의 한계를 극복하기 위해.
• 게이지 이론 기반 중력 모델의 응용에 적합한 확장 가능한 알고리즘을 개발하기 위해.
• 트레이스 공식의 계산 복잡도를 감소시키기 위해 계수 간의 재귀 관계를 유도하기 위해.

제안 방법

• 트레이스는 n의 정수 분할 s에 대한 합으로 표현되며, 각 항은 수치적 계수 αs에 의해 가중된다.
• 계수 αs는 랜덤 텐서 입력 B^{ab}를 사용해 트레이스 공식을 평가함으로써 유도된 선형 방정식계를 풀어 계산된다.
• 계수 αs 간의 재귀 관계가 도출되어 계산 복잡도를 감소시키는 데 활용된다.
• 알고리즘은 모든 n ≤ 25에서 높은 성능을 유지하는 최소 형태로 구현된다.
• 기존의 기호 간소화가 아닌 수치 평가와 선형 대수를 기반으로 하여 기호 폭발을 피한다.
• 알고리즘은 랜덤 텐서 구성에 대한 테스트를 통해 검증되었으며, 기존 알려진 결과와의 일관성을 확인했다.

실험 결과

연구 질문

• RQ1 큰 n에 대해 n개의 딜라 클라크 행렬의 대칭화된 곱의 트레이스를 어떻게 효율적으로 계산할 수 있는가?
• RQ2 트레이스 공식의 분할 기반 전개에서 계수 αs의 구조는 어떠한가?
• RQ3 계수 αs 간의 재귀 관계를 도출하여 계산을 가속화할 수 있는가?
• RQ4 n > 10일 경우 제안된 알고리즘이 기존 컴퓨터 대수 시스템에 비해 성능 면에서 어떻게 비교되는가?
• RQ5 이 방법을 사용해 신뢰성 있고 신속하게 트레이스를 계산할 수 있는 최대 n의 값은 얼마인가?

주요 결과

• 표준 데스크톱 컴퓨터에서 모든 n ≤ 25에 대해 딜라 클라크 행렬의 대칭화된 곱의 트레이스를 1분 이내에 계산할 수 있다.
• 계수 αs는 랜덤 텐서 입력 B^{ab}로부터 유도된 선형 시스템을 풀어 결정되며, 이는 수치적 안정성과 정확성을 보장한다.
• αs 계수 간의 재귀 관계가 도출되어 계산 시간을 크게 감소시켰다.
• 기존 컴퓨터 대수 시스템에서 큰 n에 대해 발생하는 메모리 오버플로우 및 성능 저하 문제를 이 방법은 피한다.
• 이 방법은 강건하고 확장 가능하여 카른-시몬스 중력 이론의 모든 예상 응용을 커버한다.
• 런타임 벤치마크를 통한 성능 검증을 통해 이론 물리학 응용 분야에서 실용적인 사용 가능성을 입증했다.