[論文レビュー] Efficient Evaluation of Cosmological Angular Statistics

Valentin Assassi, Marko Simonović+1 more|arXiv (Cornell University)|May 14, 2017

Cosmology and Gravitation Theories参考文献 28被引用数 12

ひとこと要約

この論文は、運動量空間の統計をべき関数基底に分解することで、ベッセル関数の解析的積分を可能にする、FFTlogに基づく高速な方法を提示している。この手法により、高振動的な球ベッセル関数の高価な数値積分を回避し、とくに角二乗スペクトルにおいて計算時間を大幅に削減できる。この方法は多重極 $\sigma$ に依存しない高速化を達成する。

ABSTRACT

Angular statistics of cosmological observables are hard to compute. The main difficulty is due to the presence of highly-oscillatory Bessel functions which need to be integrated over. In this paper, we provide a simple and fast method to compute the angular power spectrum and bispectrum of any observable. The method is based on using an FFTlog algorithm to decompose the momentum-space statistics onto a basis of power-law functions. For each power law, the integrals over Bessel functions have a simple analytical solution. This allows us to efficiently evaluate these integrals, independently of the value of the multipole $\ell$. In particular, this method significantly speeds up the evaluation of the angular bispectrum compared to existing methods. To illustrate our algorithm, we compute the galaxy, lensing and CMB temperature angular power spectrum and bispectrum.

研究の動機と目的

• 高振動的な球ベッセル関数に起因する、宇宙論的観測量の角統計量の評価における計算ボトルネックを解消すること。
• 銀河、弱引力レンズ、CMB非等方性などの観測量について、角スペクトルおよび角三乗スペクトルを高速かつ正確に評価できる手法を開発すること。
• 運動量空間の統計量をべき関数基底に変換することで、ベッセル積分が解析的に解けるようにし、数値的コストを低減すること。
• さまざまな宇宙論的観測量および相関関数に適用可能な汎用的かつ再利用可能なフレームワークを提供すること。
• 標準的な数値積分と比較して、とくに角三乗スペクトルにおいて顕著な高速化を実証すること。

提案手法

• FFTlogアルゴリズムを用いて、運動量空間の統計量（例えば、スペクトル関数や三乗スペクトル関数）を $k$ のべき関数基底に射影する。
• 各べき関数成分について、球ベッセル関数 $j_\ell(k\chi)$ に対する積分を解析的に評価し、数値台形則による数値積分を回避する。
• 全角相関関数は、これらの解析的に解ける成分の重み付き和として計算され、問題が滑らかな時間積分に還元される。
• 運動量空間における $n$ 点相関関数の可分形を活用し、$f_i(k_i, z_i)$ の因数分解された依存性を仮定する。
• 高 $\ell$ における数値的安定性を確保するため、特殊関数 $\mathbf{I}_\ell(\nu, t)$ の再帰的評価を用い、精度を保つために再初期化を実施する。
• 変換された一般化超幾何関数を用いることで、$t^p \mathbf{I}_\ell(\nu, t)$ を含む積分を効率的に評価し、収束が速い級数展開を実現する。

実験結果

リサーチクエスチョン

• RQ1 高振動的なベッセル関数の直接的数値積分を回避することで、宇宙論的観測量の角スペクトルおよび角三乗スペクトルをより効率的に計算できるか？
• RQ2 観測量の運動量空間統計量をどのように分解すれば、球ベッセル関数積分を解析的に評価できるか？
• RQ3 この手法は、標準的な数値積分と比較して、とくに角三乗スペクトルにおいてどの程度の性能向上を達成するか？
• RQ4 この手法は、銀河数密度、弱引力レンズポテンシャル、CMB温度非等方性など、多様な観測量に一貫して適用可能か？
• RQ5 高多重極 $\ell$ においても、速度を犠牲にせずに安定性と正確性を保証するための数値的戦略は何か？

主な発見

• この手法により、多重極 $\ell$ に依存しない高速な角スペクトルおよび角三乗スペクトルの評価が可能となり、ベッセル関数の $\ell$ 依存性に伴う計算時間の増加が生じない。
• 銀河トモグラフィーにおいて、$N_\ell = 200$ の多重極を用いて、ラップトップで30秒で角スペクトルが計算可能である。
• CMBレンズスペクトルは $N_\ell = 200$ で5秒で計算可能であり、低 $\ell$ 観測量に対して高い効率性を示している。
• CMB非等方性の角三乗スペクトルは、$10^5$ 個の三角形、$N_\nu = 70$、$N_\chi = 200$、$N_r = 80$ を用いて10分で計算可能である。
• 変換された超幾何関数の形を用いることで、特に高 $\ell$ においても収束が速く、再帰的評価による誤差制御により高い精度が達成されている。
• この手法は、運動量空間における $n$ 点関数が可分形であるあらゆる観測量に一般化可能であり、クロス相関を含む応用も可能である。