[논문 리뷰] Algorithms to solve coupled systems of differential equations in terms of power series

Jakob Ablinger, Arnd Behring외 3명|arXiv (Cornell University)|2016. 01. 01.

Numerical methods for differential equations참고 문헌 36인용 수 9

한 줄 요약

이 논문은 양자장론에서 발생하는 결합된 미분방정식계를 해결하기 위한 두 가지 고도화된 알고리즘 전략을 제시한다. 특히 질량이 있는 3계열 파인만 적분에 대해 적용된다. 거듭제곱 급수 전개와 기호 합산 기법을 활용하여 계수를 초월하는 하이퍼지오메트릭 곱의 중첩 합으로 효율적으로 계산하며, 새로운 전략이 재귀 순서를 감소시키고 초기값 요구 조건이 높은 계산적 블로킹을 최소화한다.

ABSTRACT

Using integration by parts relations, Feynman integrals can be represented in terms of coupled systems of differential equations. In the following we suppose that the unknown Feynman integrals can be given in power series representations, and that sufficiently many initial values of the integrals are given. Then there exist algorithms that decide constructively if the coefficients of their power series representations can be given within the class of nested sums over hypergeometric products. In this article we will work out the calculation steps that solve this problem. First, we will present a successful tactic that has been applied recently to challenging problems coming from massive 3-loop Feynman integrals. Here our main tool is to solve scalar linear recurrences within the class of nested sums over hypergeometric products. Second, we will present a new variation of this tactic which relies on more involved summation technologies but succeeds in reducing the problem to solve scalar recurrences with lower recurrence orders. The article will work out the different challenges of this new tactic and demonstrates how they can be treated efficiently with our existing summation technologies.

연구 동기 및 목표

• 3계열 파인만 적분에서 유도된 결합된 미분방정식계를 해결할 때 발생하는 계산적 블로킹 문제를 해결한다. 특히 초기값 계산이 불가능한 경우를 대비한다.
• 결합되지 않은 시스템에서 유도된 스칼라 재귀식의 재귀 순서를 감소시키는 새로운 알고리즘 전략을 개발하여 필요한 초기값의 수를 줄인다.
• Sigma 및 SolveCoupledSystem과 같은 기존 도구가 고순서 재귀나 메모리 집약적인 분리 단계를 처리할 때 겪는 한계를 극복한다.
• 해결책이 초월적 또는 대수적 매개변수를 포함하는 하이퍼지오메트릭 곱의 중첩 합으로 표현될 수 있는지 알고리즘적으로 결정할 수 있도록 한다.
• 1차 재귀 가능성을 체계적으로 판단하는 프레임워크를 제공하고, 반복적 방법 외의 다른 해법이 필요한 시스템을 식별한다.

제안 방법

• Zürcher의 방법을 구현한 OreSys를 사용하여 결합된 미분방정식계를 분리 알고리즘을 통해 스칼라 재귀식으로 변환한다.
• SolveCoupledSystem 패키지 내부에서 기호 합산 도구인 Sigma 및 SumProduction을 적용하여 하이퍼지오메트릭 곱의 중첩 합으로 스칼라 재귀식을 해결한다.
• 차별화된 전략을 도입하여 미분 연산자를 그 최대공약수로 나누어 재귀 순서를 감소시키며, 필요한 초기값의 수를 크게 줄인다.
• 상한 A를 제어 가능한 범위로 잘라낸 무한합을 사용하여 중첩 합 간의 대수적 의존성을 제거하고, 이질적 항이 사라지도록 보장한다.
• 특정한 악성 분모(예: a≠0,±1인 1−ax 형태)를 가진 항들만 걸러내고 처리하며, 이러한 항들은 대수적 관계에 기인할 가능성이 높다.
• 문제가 되는 합을 제거한 후 A→∞의 극한을 취하여, 대수적으로 독립적인 중첩 합과 처리 가능한 복잡도를 가진 간소화된 표현식을 도출한다.

실험 결과

연구 질문

• RQ1 3계열 파인만 적분에서 유도된 결합된 미분방정식계는 하이퍼지오메트릭 곱의 중첩 합으로 알고리즘적으로 해결될 수 있는가?
• RQ2 결합되지 않은 시스템에서 유도된 스칼라 재귀식의 재귀 순서는 어떻게 감소시킬 수 있으며, 이를 통해 필요한 초기값의 수를 최소화할 수 있는가?
• RQ3 초월적 또는 대수적 매개변수 존재 하에서 중첩 합 간의 대수적 의존성을 효과적으로 제거할 수 있는 기호 합산 기법은 무엇인가?
• RQ4 비정수 또는 복소수 극점에 기인한 이질적 항(불필요한 기여)은 어떤 조건에서 거듭제곱 급수 해에서 사라지는가?
• RQ5 고순서 재귀에 대해 기존 방법에 비해 메모리 사용량과 계산 시간 측면에서 새로운 알고리즘 접근법이 우수한가?

주요 결과

• 새로운 전략은 미분 연산자의 최대공약수로 나누어 재귀 순서를 감소시키며, 기존 표준 방법에 비해 훨씬 낮은 재귀 순서를 달성한다.
• 모든 테스트 사례에서 악성 분모(예: a≠0,±1인 1−ax 형태)를 가진 모든 중첩 합이 대수적 간소화 후 사라지며, 이질적 항의 기여가 제거된다.
• 최종 표현식에서 중첩 합의 수가 다룰 수 있는 수준으로 줄어들었고, 남은 모든 합은 대수적으로 독립적이며 효율적인 기호 계산이 가능하다.
• 이전에는 고순서 재귀(예: 16 이상)로 인해 초기값 계산이 불가능하여 실패했던 시스템도 이 방법을 통해 해결할 수 있다.
• SolveCoupledSystem에 통합된 알고리즘 파ip라인은 합성항에 최대 1000개의 중첩 합을 처리할 수 있으며, 확장성을 입증한다.
• 이 방법은 1차 재귀 가능성을 알고리즘적으로 판단할 수 있게 하여, 반복적 해법 외의 접근이 필요한 시스템을 식별할 수 있다.