[论文解读] A refined machinery to calculate large moments from coupled systems of linear differential equations

J. Blümlein, Peter Marquard+1 more|arXiv (Cornell University)|Jan 1, 2019

Numerical methods for differential equations参考文献 67被引用 2

一句话总结

本文提出优化算法，大幅减少计算耦合线性微分方程组大阶矩时所需的初始值数量——此类问题在高能物理计算中常见。通过利用系统结构中的多项式因子及迭代求解策略，该方法可高效计算高达 8,000 阶的矩，进而猜测并求解最小递推关系，最终以嵌套和（如调和和）的形式表达物理量。

ABSTRACT

The large moment method can be used to compute a large number of moments of physical quantities that are described by coupled systems of linear differential equations. Besides these systems the algorithm requires a certain number of initial values as input, that are often hard to derive in a preprocessing step.Thus a major challenge is to keep the number of initial values as small as possible. We present the basic ideas of the underlying large moment method and present refined versions that reduce significantly the number of required initial values.

研究动机与目标

• 解决高能物理中基于矩的方法计算初始值时计算瓶颈过高的问题。
• 减少计算耦合线性微分方程组大阶矩所需初始值的数量。
• 实现对描述物理量的此类系统，高达 µ = 8000 阶矩的计算。
• 促进最小递推关系的发现及其以嵌套和形式的求解，这对以闭式表达物理结果至关重要。
• 通过剔除冗余或非嵌套和贡献，支持复杂物理振幅（如极化异常维数和质量项形式因子）的计算。

提出的方法

• 应用优化的大矩方法，计算耦合线性 ODE 解的 ε 展开中系数函数的矩 Fj,k(n)。
• 对计算出的矩使用递推猜测工具（例如 [42] 中的方法），为每个 Fj,k(n) 推导出最小阶线性递推关系。
• 利用系数矩阵中的多项式公共因子 p(x) 降低递推阶数，从而减少所需初始值的数量。
• 实施迭代求解策略（精化 3），逐步计算 fλ(x,ε), fλ−1(x,ε), ..., f1(x,ε)，使用齐次项依赖性降低的低阶微分方程。
• 应用递推求解工具，将猜测的递推关系转换为不定嵌套和（如调和和）的表示形式。
• 使用 Zürcher 算法或 Gauss 消去法（通过 OreSys 实现）将系统解耦为 f1(x,ε) 的单个标量 ODE，随后通过回代求解其他分量。

实验结果

研究问题

• RQ1 在耦合线性微分方程组中，如何最小化大矩计算所需的初始值数量？
• RQ2 系统的哪些结构特性（如系数中的多项式因子）可被利用以降低递推阶数并减少初始值数量？
• RQ3 与全系统方法相比，迭代求解策略（如先求解 fλ，再 fλ−1，依此类推）在多大程度上可降低矩计算的复杂度？
• RQ4 优化的大矩方法是否能成功计算复杂物理振幅（如极化三圈异常维数）的高阶矩（例如 µ = 8000）？
• RQ5 该方法能否可靠地从主积分解中产生的大量复杂特殊函数（如椭圆或超几何函数）中筛选并仅表达相关的嵌套和分量？

主要发现

• 通过引入一个 13 次多项式因子 p(x)，将 F1,−3(n) 的递推阶数从 17 降低至 4，显著减少了所需初始值的数量。
• 该方法使质量项三圈形式因子的计算达到 8,000 阶矩，实现了完整的递推猜测与嵌套和形式的求解。
• 对于极化三圈异常维数，该方法使用高达 µ = 6,000 阶矩成功复现了结果，验证了方法的正确性与高效性。
• 质量项三圈算符矩阵元 A(3)Qg 的所有 ε⁻¹ 项及常数项的大部分均已成功表达为不定嵌套和的形式。
• 通过使用精化 2 和精化 3，特别是利用公共多项式因子与迭代求解，每个分量 Fj,k(n) 所需的初始值数量显著减少。
• 该方法成功剔除了非嵌套和分量，仅在最终结果中保留了物理上相关的调和和表示形式。