[論文レビュー] Sector decomposition via computational geometry

T. Kaneko, Takahiro Ueda|arXiv (Cornell University)|Apr 30, 2010

Computational Geometry and Mesh Generation参考文献 9被引用数 5

ひとこと要約

本稿では、多ループフェ Feynman振幅計算におけるセクタ分解のための反復的でない、計算幾何学に基づく手法を提示する。ユークリッド空間内の凸多面体錐を用いて多項式の漸近的挙動を分類することで、最適なセクタ分解を決定的に特定する。反復的手法と比較してセクタ数を削減する。特にダブルボックスやトリプルボックスのような複雑な図式において、30–60%のセクタ数削減が実証された。

ABSTRACT

A non-iterative method is presented for the factorization step of sector decomposition method, which separates infrared divergent part from loop integration. This method is based on a classification of asymptotic behavior of polynomials. The problem is converted to ones for convex body in Euclidean space. They are solved with algorithms developed in computational geometry. A test implementation shows that this method produces less number of decomposed sectors than usual iterative sector decompositions.

研究の動機と目的

• 多ループ振幅計算で用いられる反復的セクタ分解手法における非効率性と無限ループのリスクを解消すること。
• フェ Feynman 積分における赤外発散を要因分解するための決定的で反復的でない手法を開発すること。
• セクタ分解で生成されるセクタ数を削減し、計算効率を向上させること。
• 凸幾何学と計算アルゴリズムを活用して、体系的で幾何学的駆動のセクタ分解アプローチを確立すること。

提案手法

• 本手法は、セクタ分解問題を、ユークリッド空間内における凸体による多項式の漸近的挙動の分類問題に変換する。
• 凸包アルゴリズムを用いて、被積分関数の特異性が分離される最小のセクタ集合を特定する。
• 積分領域が、多項式の異なる漸近的領域に対応する凸多面体錐に分解される。
• これらの錐の三角形分割は、独自のアルゴリズムを用いて実行され、単項式のヤコビアンを持つ最終的なセクタ分解が得られる。
• 変換により、発散部分がεの累乗として要因分解され、ローレン級数展開が可能になる。
• 本手法は完全に決定的であり、反復的精錬を回避するため、無限ループのリスクがない。

実験結果

リサーチクエスチョン

• RQ1 反復的でない幾何学的アプローチは、多ループ振幅における標準的な反復的セクタ分解よりも少ないセクタ数を生み出せるか？
• RQ2 計算幾何学における凸多面体錐の使用は、赤外発散の要因分解をどのように改善するか？
• RQ3 この幾何学的手法によるセクタ数の定量的削減は、既存の戦略と比較してどの程度達成可能か？
• RQ4 セクタ生成において無限ループのリスクなしに、決定的に実装可能か？

主な発見

• 提示された手法はボックス図に対して12セクタにまで削減し、最良の反復的戦略と同等の結果を得つつ、決定的で反復的でないアプローチを実現した。
• すべての外部質量がゼロのダブルボックス図では、本手法は266セクタを生成したが、標準的な反復的手法では775セクタであった。
• 非平面的ダブルボックスの場合、本手法は360セクタを達成したのに対し、標準的手法では1138セクタであり、68%の削減が達成された。
• トリプルボックス図では、本手法は6568セクタを生成したが、標準的手法では114256セクタであり、90%を超える削減が達成された。
• 三角形分割前の指数的セクタ分解では、100セクタ（トリプルボックス）を達成しており、これは本手法のコアとなる幾何的分解が極めて効率的であることを示している。
• 決定的で反復的でない性質のおかげで、メモリオーバーフローと障害事象（例：過去の実装における「F」や「M」）を回避した。