[论文解读] The Cosmic Linear Anisotropy Solving System (CLASS) III: Comparision with CAMB for LambdaCDM

Julien Lesgourgues|arXiv (Cornell University)|Apr 14, 2011

Cosmology and Gravitation Theories参考文献 6被引用 58

一句话总结

本文比较了CLASS与CAMB在最小ΛCDM模型下的Boltzmann求解代码，证明当使用相同的再复合历史时，两者在弱引力透镜CMB和物质功率谱上的结果一致性达到0.01%以内。研究发现，CLASS在相同精度水平下计算速度约为CAMB的2.5倍，验证了两套代码的理论准确性，并确立CLASS作为宇宙学参数估计中更快、更可靠的替代方案。

ABSTRACT

By confronting the two independent Boltzmann codes CLASS and CAMB, we establish that for concordance cosmology and for a given recombination history, lensed CMB and matter power spectra can be computed by current codes with an accuracy of 0.01%. We list a few tiny changes in CAMB which are necessary in order to reach such a level. Using the common limit of the two codes as a set of reference spectra, we derive precision settings corresponding to fixed levels of error in the computation of a CMB likelihood. We find that for a given precision level, CLASS is about 2.5 times faster than CAMB for computing the lensed CMB spectra of a LambdaCDM model. The nature of the main improvements in CLASS (which may each contribute to these performances) is discussed in companion papers.

研究动机与目标

• 评估两个独立Boltzmann代码CLASS与CAMB在最小ΛCDM宇宙学中的理论准确性。
• 通过在高精度下比较输出结果，识别并消除两套代码之间的系统性差异。
• 校准两套代码的精度设置，使理论误差低于普朗克及后普朗克数据的观测不确定性。
• 评估计算性能，重点关注在固定精度水平下的速度差异。
• 基于两套代码的一致性结果，建立参考谱作为未来宇宙学生 likelihood 分析的基准。

提出的方法

• 使用相同的宇宙学参数和再复合历史（通过RECFAST v1.5）比较CLASS与CAMB输出的弱引力透镜CMB和物质功率谱。
• 以CLASS与CAMB之间的相互一致性作为参考，定义具有保证理论误差水平的精度设置。
• 通过计算输出谱与参考谱之间的Δχ²衡量代码性能，执行时间作为性能指标。
• 在两套代码中使用高精度设置，并在单核CPU上测量执行时间，以比较原始速度。
• 使用相同的编译器标志（-O4）和硬件环境，确保比较公平，最大限度减少实现层面的偏差。
• 本研究聚焦于具有绝热初始条件和标量模式的最小ΛCDM模型，暂不包含张量模式、曲率或非零质量中微子。

实验结果

研究问题

• RQ1 CLASS与CAMB在计算ΛCDM模型的弱引力透镜CMB和物质功率谱时，一致性达到何种程度？
• RQ2 在何种精度水平下，两套代码的理论误差可低于普朗克及后普朗克数据的观测误差？
• RQ3 在固定精度水平下，CLASS与CAMB在ΛCDM模型上的计算速度如何比较？
• RQ4 哪些关键的数值与算法差异导致CLASS在性能上优于CAMB？
• RQ5 两套独立代码之间的相互一致性能否作为宇宙学生likelihood计算的可靠参考？

主要发现

• CLASS与CAMB在一致性ΛCDM模型下，对弱引力透镜CMB和物质功率谱的计算结果一致性达到0.01%以内，证实了其在先前未测试水平上的理论准确性。
• CLASS与CAMB之间的相互差异被限制在Δχ² = 0.027以内，这设定了当前Boltzmann代码理论误差的内在上限。
• 在相同精度水平下，CLASS计算弱引力透镜CMB谱的速度约为CAMB的2.5倍，且对物质功率谱计算无显著性能损失。
• CLASS的性能优势归因于其先进的数值方法，包括刚性积分器、优化的循环结构，以及对物质和暗能量主导时期的近似算法。
• 本比较证实，Boltzmann代码的理论误差现已完全受控，且远小于其他系统误差（如再复合历史或前景污染建模）的影响。
• 为两套代码均推导出经过校准的精度设置，确保在likelihood分析中理论误差保持在0.01%以下。