[論文レビュー] 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)を活用して再帰関係の次数を低下させ、これにより必要な初期値の数を削減する。
• 反復的解法戦略（Refinement 3）を実装し、fλ(x,ε), fλ−1(x,ε), ..., f1(x,ε)を段階的に計算することで、非同次項依存性が小さい低次の微分方程式を用いる。
• 再帰関係解法ツールを用いて、推測された再帰関係を不定ネスト和（例：調和和）の表現に変換する。
• ZürcherのアルゴリズムまたはGaussの消去法（OreSysを介して）を用い、系をf1(x,ε)の単一スカラーODEに分離し、その後他の成分に対して後退代入を実行する。

実験結果

リサーチクエスチョン

• RQ1 結合された線形微分方程式系において、大規模モーメント計算に必要な初期値の数を最小化する方法は何か？
• RQ2 再帰関係の次数を低下させ、結果として初期値の数を減らすために、系の構造的性質（例：係数に含まれる多項式因子）をどのように活用できるか？
• RQ3 反復的解法戦略（例：fλから順にfλ−1, ..., を解く）は、全系アプローチと比較して、モーメント計算の複雑さをどの程度軽減できるか？
• RQ4 洗練された大モーメント法は、偏極3ループ異常次元のような複雑な物理振幅について、高次のモーメント（例：µ = 8000）を正しく計算できるか？
• RQ5 マスタインテグラルの解から生じる複雑な特殊関数（例：楕円関数や超幾何関数）の大量の項の中から、非ネスト和の成分を信頼性高く除外し、物理的に関連のあるネスト和表現のみを保持できるか？

主な発見

• 13次多項式因子p(x)を導入することで、F1,−3(n)の再帰関係次数を17から4にまで低下させ、必要な初期値の数を著しく削減した。
• 本手法により、質量のある3ループフォーム因子について最大8,000モーメントの計算が可能となり、ネストされた和の形での完全な再帰関係の推測と解法が実現した。
• 偏極3ループ異常次元の計算において、最大µ = 6,000モーメントを用いて結果を再現し、本手法の正しさと効率性を確認した。
• 質量のある3ループオペレータ行列要素A(3)Qgのε⁻¹項および定数項の大部分が、不定ネスト和の形で正しく表現された。
• 特に共通多項式因子の使用と反復的解法の組み合わせにより、各成分Fj,k(n)の初期値の数が顕著に削減された（Refinement 2およびRefinement 3の活用）。
• 本手法は非ネスト和成分を効果的に除外し、最終的な結果に物理的に関連のある調和和表現のみを保持した。