[論文レビュー] Extrapolation Algorithms for Infrared Divergent Integrals

Elise de Doncker, J. Fujimoto+6 more|arXiv (Cornell University)|Oct 17, 2011

Numerical methods in inverse problems参考文献 13被引用数 5

ひとこと要約

本稿では、量子場の理論における赤外発散を示す1ループ頂点積分の主要係数を計算するために、4倍精度算術を用いた二重外挿法を提示する。正則化パrameter ε および δ を小さくする積分列に対して非線形外挿（ε-アルゴリズムを用いて）を適用することで、積分領域の内部および境界に特異点が存在する場合でも、ほぼ機械精度に近い精度を達成する。

ABSTRACT

This paper describes applications of extrapolation for the computation of coefficients in an expansion of infrared divergent integrals. An extrapolation procedure is performed with respect to a parameter introduced by dimensional regularization. While this treats typical IR singularities at the boundaries of the integration domain, special care needs to be taken in cases where the integrand is singular in the interior of the domain as well as on the boundaries. A double extrapolation is devised for a class of massless vertex integrals. Quadruple precision results are presented, demonstrating high accuracy. The computations are supported by the use of general adaptive integration programs from the QUADPACK package, in iterated integrations with highly singular integrand functions.

研究の動機と目的

• 非可積分な特異点が存在するため、標準的な統合手法が失敗する赤外発散ループ積分の数値的計算における課題に対処すること。
• 特に質量ゼロの頂点図において生じる、領域の境界および内部の両方で特異点を持つ被積分関数の困難さを克服すること。
• 外挿技術を用いて、このような積分のローレンツ展開における係数を計算するための堅牢な数値フレームワークを構築すること。
• 4倍精度算術とQuadpackの適応的統合ルーチンを用いて、高精度の結果を示すこと。
• 特異点構造の事前解析なしに、より自動的かつ信頼性の高い複雑なループ積分の数値計算を可能にすること。

提案手法

• 赤外発散を正則化するため、小さなパrameter ε を用いた次元正則化を用い、積分を ε におけるローレンツ級数に展開する。
• Dqags/Dqag を用いた繰り返し数値統合により、ε の値を小さくした場合の積分近似値を計算する。
• ε → 0 のときの積分近似列の収束を加速するために、ε-アルゴリズムを用いた非線形外挿を実装する。
• 内部特異点を有する積分の場合、分母における小さな虚数項 iδ の導入により、δ に関する第二の外挿を実施して非可積分特異点に対処する。
• 収束性と安定性を向上させるために、ε および δ に対して幾何級数（例：Bulirsch 級数）を用いる。
• Quadpack が複素関数をネイティブにサポートしていないため、複素積分の実部と虚部を別々に計算し、それぞれに外挿を適用する。

実験結果

リサーチクエスチョン

• RQ1 積分領域の境界および内部に特異点が存在する場合、赤外発散頂点積分の数値統合をどのように安定化できるか。
• RQ2 外挿技術は、このような積分のローレンツ展開における係数抽出の精度をどの程度向上できるか。
• RQ3 ε および δ に関する二重外挿は、複雑なループ積分に対してほぼ機械精度の結果を達成するのにどの程度有効か。
• RQ4 幾何級数と調和級数の異なる系列が、外挿プロセスの安定性および収束性に与える影響は何か。
• RQ5 4倍精度算術と適応的統合、および外挿を組み合わせることで、高精度な量子場の理論計算に信頼性のある結果を得られるか。

主な発見

• 二重外挿法により、赤外発散積分の係数が相対誤差 5.84×10⁻¹⁶ まで計算可能であり、4倍精度算術の限界に近い。
• 1つのオンシェルおよび2つのオフシェル粒子の場合、係数 $ C_{-1} $ の外挿された虚部は 1.421026277606×10⁻⁴ に収束し、誤差推定値は 5.84×10⁻¹⁶ である。
• 8回の外挿ステップを経て、$ C_0 $ の実部は 5.24199463865×10⁻⁴ に計算され、誤差推定値は 5.64×10⁻¹⁴ である。
• 虚部の収束が実部よりも速いため、特異点の性質に起因する異なる収束挙動が示された。
• ε-アルゴリズムを用いた固定長（18）の外挿列は安定した結果を提供するが、将来的な改善のためには動的終了基準の導入が提案される。
• 本手法は、複雑な特異点を有する積分に対しても高精度な結果が得られることを示しており、ループ積分の自動化された数値計算の可能性を裏付けている。