[論文レビュー] Loaded Dice in Monte Carlo: importance sampling in phase space integration and probability distributions for discrepancies

A. van Hameren|ArXiv.org|Jan 10, 2001

Mathematical Approximation and Integration参考文献 14被引用数 4

ひとこと要約

本学位論文は、多粒子QCD過程における位相空間積分のための新規重要度サンプリングアルゴリズムを開発し、運動量のアンテナ極構造を活用することで、次元mに対してO(m)のスケーリングを達成し、拒否サンプリングと比べて著しく優れた性能を発揮する。また、量子場理論的手法を用いて2次的不均一性の生成関数を導出し、準モンテカルロ積分の均一性に関する厳密な解析を可能にする。

ABSTRACT

Discrepancies play an important role in the study of uniformity properties of point sets. Their probability distributions are a help in the analysis of the efficiency of the Quasi Monte Carlo method of numerical integration, which uses point sets that are distributed more uniformly than sets of independently uniformly distributed random points. In this thesis, generating functions of probability distributions of quadratic discrepancies are calculated using techniques borrowed from quantum field theory. The second part of this manuscript deals with the application of the Monte Carlo method to phase space integration, and in particular with an explicit example of importance sampling. It concerns the integration of differential cross sections of multi-parton QCD-processes, which contain the so-called kinematical antenna pole structures. The algorithm is presented and compared with RAMBO, showing a substantial reduction in computing time. In behalf of completeness of the thesis, short introductions to probability theory, Feynman diagrams and the Monte Carlo method of numerical integration are included.

研究の動機と目的

• 運動量のアンテナ極構造を持つ多粒子QCD過程における位相空間点の生成のため、効率的で均一にスケーラブルなアルゴリズムの開発。
• 2次的不均一性の確率分布を用いて、準モンテカルロ積分に用いられる点集合の均一性特性を分析すること。
• 特に経路積分とインスタントンを含む量子場理論の技術を用いて、不均一性分布の生成関数を計算すること。
• 高次元位相空間積分における標準的拒否サンプリングの計算効率の良い代替手法を提供すること。
• 高度な数学的物理道具を用いて、不均一性に基づく誤差推定の理論的基盤を確立すること。

提案手法

• 運動的制約によって定義されるm次元ポリトープ内に一様に分布する運動量を生成するユニタリなアルゴリズム形式を構築する。
• アンテナ関数の構造に基づく統合領域の再帰的分解を用い、ベータ分布および正弦分布による効率的サンプリングを実現する。
• 三角関数積分およびベータ関数の恒等式を用いてm次元ポリトープの体積V_m(k)を導出し、閉形式の表現を得る。
• 量子場理論の生成関数技法を応用し、LegoおよびL₂*不均一性を含む2次的不均一性の確率分布を計算する。
• インスタントン解析およびリーマン面構造の計算を実施し、不均一性分布の漸近的挙動を検討する。
• 有限標本補正および強収束極限を導出し、実用的な数値応用における精度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 多粒子QCD過程における位相空間積分を、次元の増加に伴って均一に効率的に行うにはどうすればよいか？
• RQ2 高次元点集合における2次的不均一性の確率分布は何か？ そして、それが解析的にどのように計算できるか？
• RQ3 経路積分やインスタントンを含む量子場理論の技術は、数値積分における不均一性分布の分析に応用可能か？
• RQ4 提案されたアルゴリズムは、高次元位相空間において、標準的拒否サンプリングと比べてどのように性能を発揮するか？
• RQ5 有限標本補正とは何か？ そして、それらは準モンテカルロ積分の精度にどのように影響するか？

主な発見

• 提案されたアルゴリズムは、m次元位相空間点を生成するにあたりO(m)回の呼び出しで済み、次元にかかわらず均一な効率性を達成する。これに対して拒否サンプリングは指数関数的O(2^m)のスケーリングを示す。
• m ≥ 4の場合、アルゴリズムのCPU時間はmに対して線形に増加するが、拒否サンプリングは指数関数的に増加し、m=10からm=20に至る間に20倍の時間増加を示す。
• アルゴリズムは相関のない運動量を生成し、任意の二つのx_i間の相関係数が正確に1/2となるため、強固な妥当性の検証が可能である。
• L₂*不均一性の生成関数は経路積分法を用いて導出され、その確率分布を正確に計算可能となる。
• Legoおよびダイアホニー不均一性に対して、有限標本補正および強収束極限が導出され、有限Nの状態でも精度が向上する。
• cos(πx/2)密度を持つ非一様分布への拡張においても、均一な効率性が維持され、重み付き体積はO(m(8/π²)^{m/2})のスケーリングを示す。