Solution of Canonical Differential Equations for Integrals on Arbitrary Geometries

技術要約：任意の幾何学における積分のための標準的微分方程式の解法

問題提起
量子場における多ループ散乱振幅の評価は、部分積分（IBP）恒等式を介してファインマン積分をマスター積分（MI）の集合へと簡約することに依存することが多い。これらのマスター積分を決定するための強力な手法は、運動学的不変量および質量に関する一次線形微分方程式の系を解くことである。$\epsilon$因子化された標準形（canonical form）の導入は、接続行列が有理係数の$d\log$形式で構成される多重対数関数のケースにおいて、このプロセスを効率化してきた。

しかし、質量の大きな粒子や複数のスケールを含む問題では、基礎となる幾何学が単純な多重対数関数を超えて、楕円曲線、カラビ・ヤウ多様体、および高次種数の構造へと拡張されることが多い。このようなケースでは、標準的な形式を実現するために、超越関数（幾何学の周期）を用いた変換が必要となる。このような形式は存在するものの、実用上のボトルネックとして、要求される特殊関数の数値的評価が困難であることが挙げられる。これは、特殊関数の閉じた形式による表現の欠如や、解析接続の複雑さに起因する。既存の数値的手法は、級数展開（マルチスケールの領域では限定的）に依存するか、あるいは複雑な特殊関数の明示的な実装を必要とするが、これらは普遍的に実現可能ではない。

手法
著者らは、超越関数の明示的な評価を回避するための戦略を提案している。核心となる観察は、標準的な基底への変換に現れる特殊関数は恣意的なものではなく、有理係数の線形微分方程式を満たす元のマスター積分から派生しているという点である。

その結果、著者らは、これらの超越的な構成要素（例：楕円周期、平方根）自体が、有理係数の微分方程式を満たすことを示している。提案された手法は以下の手順で行われる：

1. 拡張（Augmentation）: マスターインテグラルの標準的微分方程式と、特殊関数を制御する補助的な微分方程式を結合させた、拡張された微分方程式の系を構築する。
2. 有理化（Rationalization）: この拡張された系は、有理係数のみを含む。特殊関数の微分は、閉じた形式の式からではなく、それらを定義する有理微分方程式から直接計算される。
3. 数値積分（Numerical Integration）: 系全体を、特異点を避けるために運動学的空間内の複素化された経路に沿って数値的に解く。特殊関数の解析接続は、それらの微分方程式の積分によって自動的に処理されるため、アドホックな解析接続アルゴリズムの必要性が排除される。
4. 境界条件（Boundary Conditions）: この手法は、非特異点における特殊関数の境界値を必要とする。これらの値は、接続行列を簡略化するために使用される代数的関係と整合している限り、任意に選択できる（例：1に設定する）。標準的なマスター積分は、これらの境界値を用いて標準的な基底へと回転される。

主な貢献

• 理論的洞察: 本論文は、超越的な項を持つ標準的な系を数値的に解くことへの障害は、見かけ上のものに過ぎないことを確立している。特殊関数を、独自の有理微分方程式を満たす動的な変数として扱うことで、問題はより大きな、純粋に有理的な系を解くことへと還元される。
• ソフトウェア実装: 著者らは、高性能な数値的評価のために設計されたスタンドアロンのC++パッケージ（diffeqs）を提供する。このソフトウェアは以下の特徴を持つ：
  • C++で記述されており、Boostライブラリのみに依存し、実行時にMathematicaなどの商用ソフトウェアへの依存を避けている。
  • QDライブラリおよびテンプレート型により、任意の浮動小数点精度をサポートしている。
  • ユーザーが接続行列、特殊関数のベクトル場、および積分経路を定義できるモジュール式インターフェースを実装している。
  • 数値的安定性を確保するための複素経路変形機能を含んでいる。
• 現実的な適用: 著者らは、この手法を非自明な物理的事例である、重いクォークのループを介したダイフォトン、ダイジェット、および$\gamma$+ジェットのハドロン生成における2ループ・マスター積分に適用している。この系には、220個のマスター積分、16個の平方根、および楕円周期に関連する超越関数が含まれる。

結果
著者らは、3つの特定の位相空間点に対して実装のテストを行った。

• 精度: 倍精度において、手法はテストされた点に対して$O(10^{-6})$および$O(10^{-9})$の絶対誤差を達成した。四倍精度では、誤差はそれぞれ$O(10^{-8})$および$O(10^{-19})$に減少しており、解析接続の正当性とこのアプローチの数値的安定性を検証している。
• パフォーマンス: 接続行列の評価時間は、精度とハードウェアに依存するものの、サブミリ秒から数ミリ秒の範囲内で効率的であることが確認された。著者らは、四倍精度バージョンは接続行列の評価に関して倍精度バージョンよりも約10倍遅いと述べている。
• 比較: 結果はAMFlowパッケージを使用して得られた値と照合され、一貫性が確認された。パフォーマンスは、AMDプロセッサと比較してApple Siliconチップ上で優れていることが観察された。

意義および主張
本論文は、この戦略により、楕円またはカラビ・ヤウ幾何を含む任意のファインマン積分に対して、標準的な微分方程式が普遍的に適用可能なツールになることを主張している。その意義は、複雑な特殊関数の閉じた形式の実装を必要とせずに、数値的な現象論を実行できる点にある。幾何学的な構成要素をそれら自身の有理微分方程式を通じて伝播させることで、この手法は特殊関数の「ブラックボックス」的な実装とそれに伴う解析接続の困難を回避している。

著者らは、本研究を、複雑な散乱過程に対する高次摂動計算を自動化するための実用的なステップとして位置づけている。彼らは、現在の実装は高度な最適化なしでも機能しているが、式の最適化や部分分数分解の戦略によってさらなる向上が可能であることを強調している。この手法は、従来の級数展開法が苦戦する、複数のスケールや質量の大きな粒子を含むプロセスにおける、大規模なコリジョン現象論への適用準備ができているものとして提示されている。