[论文解读] 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)$ 依赖结构。
• 通过递归计算特殊函数 $\mathbf{I}_\ell(\nu, t)$ 来处理高 $\ell$ 下的数值稳定性，必要时重新初始化以保证精度。
• 通过变换后的广义超几何函数形式，高效计算涉及 $t^p \mathbf{I}_\ell(\nu, t)$ 的积分，其级数具有快速收敛性。

实验结果

研究问题

• RQ1 能否通过避免对振荡贝塞尔函数进行直接数值积分，更高效地计算宇宙学可观测量的角功率谱和三粒子谱？
• RQ2 如何对可观测量的动量空间统计量进行分解，以实现球贝塞尔函数积分的解析求解？
• RQ3 与标准数值积分相比，该方法在角三粒子谱方面性能提升如何？
• RQ4 该方法能否统一应用于多种可观测量，如星系数密度、弱引力透镜势和CMB温度各向异性？
• RQ5 采用何种数值策略可在不牺牲速度的前提下，确保高多极 $\ell$ 下的稳定性和精度？

主要发现

• 该方法实现了与多极 $\ell$ 无关的角功率谱和三粒子谱的快速计算，计算时间不随贝塞尔函数的 $\ell$-依赖性而增长。
• 对于星系分层观测，使用 $N_\ell = 200$ 个多极时，角功率谱在笔记本电脑上仅需30秒即可完成计算。
• CMB弱引力透镜功率谱的计算耗时仅5秒，$N_\ell = 200$，表明该方法在低-$\ell$可观测量中具有极高效率。
• CMB各向异性角三粒子谱的计算耗时10分钟，处理 $10^5$ 个三角形，参数为 $N_\nu = 70$，$N_\chi = 200$，$N_r = 80$。
• 通过采用变换后的超几何函数形式，该方法实现了高精度计算，尤其在高 $\ell$ 时收敛迅速，且通过递归计算实现误差控制。
• 该方法具有通用性，适用于任何在动量空间中具有可分离 $n$-点函数的可观测量，包括交叉相关函数。