[论文解读] Algorithms to solve coupled systems of differential equations in terms of power series

Jakob Ablinger, Arnd Behring+3 more|arXiv (Cornell University)|Jan 1, 2016

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）在处理高阶递推或内存密集型解耦步骤时的局限性。
• 实现算法化判断解是否可表示为超几何乘积的嵌套和（包括调和和分圆和）的形式。
• 提供一种系统性框架，用于判断系统的一阶可约性，并识别出需要采用替代求解方法的系统。

提出的方法

• 通过Zürcher方法（在OreSys中实现）的解耦算法，将耦合微分系统转化为标量递推。
• 在SolveCoupledSystem软件包中应用符号求和工具——Sigma和SumProduction，以超几何乘积的嵌套和形式求解标量递推。
• 实施一种新策略：将微分算子除以其最大公因式（GCD），以降低递推阶数，显著减少所需初始值的数量。
• 使用截断的无穷级数并控制上界A，以消除嵌套和之间的代数依赖性，确保外来的项消失。
• 应用高级简化技术，仅筛选和处理含有“坏”分母（如1−ax，其中a≠0,±1）的子表达式，这些分母易引发代数关系。
• 在消除有问题的和之后，取A→∞的极限，得到一个由代数独立的嵌套和构成的简化表达式，且复杂度可控。

实验结果

研究问题

• RQ1 能否以超几何乘积的嵌套和形式，算法化地求解来自3环费曼积分的耦合微分方程组？
• RQ2 如何在从解耦系统导出的标量递推中降低递推阶数，以最小化所需初始值的数量？
• RQ3 在存在超越或代数参数的情况下，哪些符号求和技术能有效消除嵌套和之间的代数依赖性？
• RQ4 在幂级数解中，外来项（来自非整数或复数极点的不想要的贡献）在何种条件下会消失？
• RQ5 在高阶递推情况下，新算法方法在内存使用和计算时间方面是否优于现有方法？

主要发现

• 新策略通过将微分算子除以其最大公因式，成功降低了递推阶数，所得递推阶数远低于标准方法。
• 在所有测试示例中，所有含‘坏’分母（如1−ax，其中a≠0,±1）的嵌套和在代数简化后均消失，从而消除了外来项的贡献。
• 最终表示中嵌套和的数量被减少到可管理的水平，且所有剩余的和均为代数独立，从而支持高效的符号计算。
• 该方法可求解以往因递推阶数过高（如16阶或以上）而无法求解的系统，此时使用现有工具计算初始值不可行。
• 集成于SolveCoupledSystem中的算法处理管道，成功处理了被求和项中包含多达1000个嵌套和的系统，展现出良好的可扩展性。
• 该方法实现了对一阶可约性的算法化判断，能够识别出需要采用非迭代求解技术的系统。