[论文解读] A direct method to compute the galaxy count angular correlation function including redshift-space distortions

J.E. Campagne, S. Plaszczynski+1 more|arXiv (Cornell University)|Mar 8, 2017

Galaxies: Formation, Evolution, Phenomena参考文献 28被引用 4

一句话总结

本文提出了一种直接、数值高效的计算方法，用于计算红移空间中的星系计数角相关函数 ξ(θ,z₁,z₂)，通过一种新颖的球谐函数加法定理，绕过了中间的功率谱计算过程。该方法实现了亚百分之一的精度，并在速度上比现有代码（如 CLASSgal）快10至20倍，使得无需假设基准模型即可实现大规模结构数据的快速、宇宙学无关分析。

ABSTRACT

In the near future, cosmology will enter the wide and deep galaxy survey area allowing high-precision studies of the large scale structure of the universe in three dimensions. To test cosmological models and determine their parameters accurately, it is natural to confront data with exact theoretical expectations expressed in the observational parameter space (angles and redshift). The data-driven galaxy number count fluctuations on redshift shells, can be used to build correlation functions $C(θ; z_1, z_2)$ on and between shells which can probe the baryonic acoustic oscillations, the distance-redshift distortions as well as gravitational lensing and other relativistic effects. Transforming the model to the data space usually requires the computation of the angular power spectrum $C_\ell(z_1, z_2)$ but this appears as an artificial and inefficient step plagued by apodization issues. In this article we show that it is not necessary and present a compact expression for $C(θ; z_1, z_2)$ that includes directly the leading density and redshift space distortions terms from the full linear theory. It can be evaluated using a fast integration method based on Clenshaw-Curtis quadrature and Chebyshev polynomial series. This new method to compute the correlation functions without any Limber approximation, allows us to produce and discuss maps of the correlation function directly in the observable space and is a significant step towards disentangling the data from the tested models.

研究动机与目标

• 开发一种数值高效、直接计算红移空间中角相关函数 ξ(θ,z₁,z₂) 的方法，避免依赖中间的功率谱计算。
• 通过球谐函数展开和加法定理，提供 ξ(θ,z₁,z₂) 的闭式解析表达式，简化推导与实现。
• 在不使用 Limber 近似的情况下，准确计算包含红移空间扭曲（RSD）、密度、多普勒和引力透镜效应的相关函数。
• 开发一个公开的软件工具（Angpow），在宇宙参数推断中比 CLASSgal 等现有代码具有更高的速度和内存效率。
• 通过直接在可观测量空间（角度与红移）中工作，实现数据驱动的、宇宙学无关的 ΛCDM-GR 宇宙学范式检验。

提出的方法

• 利用球谐函数展开和加法定理，推导出 ξ(θ,z₁,z₂) 的闭式表达式，实现对多重极的直接求和。
• 将相关函数表示为勒让德多项式上的无穷级数，避免因有限 ℓ_max 截断引起的振荡。
• 通过 3C-算法高效评估，无需显式计算 Cℓ(z₁,z₂)，从而从角功率谱重建相关函数。
• 在谐函数空间公式中引入额外的导数和算子，系统性地纳入相对论效应，如 RSD、多普勒频移和引力透镜。
• 使用 ℓ_max = 3500 且无加窗函数的基准宇宙学模型，与 CLASSgal 代码进行方法验证。
• 在公开的 Angpow 软件中实现该方法，针对多核系统进行了低内存和高速计算的优化。

实验结果

研究问题

• RQ1 能否在不计算中间 Cℓ(z₁,z₂) 功率谱的情况下，推导出 ξ(θ,z₁,z₂) 的直接闭式表达式？
• RQ2 如何系统性地将红移空间扭曲、多普勒频移和引力透镜等相对论效应纳入相关函数形式化框架？
• RQ3 新方法是否能在数值效率和精度上超越 CLASSgal 等现有方法？
• RQ4 该方法在多大程度上通过直接在可观测量空间工作，消除了数据处理中对基准宇宙学模型的依赖？
• RQ5 红移选择函数的选择（狄拉克型 vs 高斯型）在多大程度上影响观测到的相关函数结构，特别是 BAO 振荡特征？

主要发现

• 新方法与 CLASSgal 代码在不同红移和角度尺度上均实现了亚百分之一的精度一致，验证了其准确性。
• 该实现比 CLASSgal 快 10 至 20 倍，在 16 个线程下仅用 18 秒完成，内存占用仅 1 GB，而 CLASSgal 在 8 个线程下需 3 分钟且占用 100 GB 内存。
• 由于高阶 Cℓ 随多极数迅速衰减，该方法天然抑制了高频效应，无需使用加窗或加apodization 处理。
• 正确再现了 RSD 和透镜效应导致的中心相关峰收缩以及在 θ ≈ 0 和 Δz ≈ ±0.05 处负相关性的出现。
• 高斯型红移选择函数在 θ ≈ 0 附近模糊了 BAO 振荡，但在更高角度和小 Δz 区域仍能保留这些特征，与理论预期一致。
• Angpow 软件已公开发布于 https://gitlab.in2p3.fr/campagne/AngPow，支持在大尺度结构研究中实现快速、宇宙学无关的参数推断。