[论文解读] Generalized Boltzmann hierarchy for massive neutrinos in cosmology

Caio Nascimento|arXiv (Cornell University)|Apr 1, 2021

Astrophysics and Cosmic Phenomena参考文献 33被引用 2

一句话总结

本文提出广义玻尔兹曼层级（GBH），通过高阶速度加权积分精确积分动量依赖性，将宇宙学中质量中微子的建模方法简化为一个双参数无穷阶常微分方程组，实现了与CLASS相比在大尺度和中等尺度上中微子转移函数的<0.5%精度，提供了一种计算效率更高的标准玻尔兹曼求解器替代方案，同时保持高精度。

ABSTRACT

Boltzmann solvers are an important tool for the computation of cosmological observables in the linear regime. In the presence of massive neutrinos, they involve solving the Boltzmann equation followed by an integration in momentum space to arrive at the desired fluid properties, a procedure which is known to be computationally slow. In this work we introduce the so-called generalized Boltzmann hierarchy (GBH) for massive neutrinos in cosmology, an alternative to the usual Boltzmann hierarchy, where the momentum dependence is integrated out leaving us with a two-parameter infinite set of ordinary differential equations. Along with the usual expansion in multipoles, there is now also an expansion in higher velocity weight integrals of the distribution function. Using a toy code, we show that the GBH produces the density contrast neutrino transfer function to a $\lesssim 0.5\%$ accuracy at both large and intermediate scales compared to the neutrino free-streaming scale, thus providing a proof-of-principle for the GBH. We comment on the implementation of the GBH in a state of the art Boltzmann solver.

研究动机与目标

• 开发一种计算高效的替代方案，用于宇宙学中质量中微子的标准玻尔兹曼求解器。
• 在中微子分布函数中精确积分动量依赖性，避免昂贵的数值动量积分。
• 在大尺度和中等尺度上，实现中微子密度扰动转移函数的亚百分之一精度。
• 使该框架能够高效集成到最先进的玻尔兹曼求解器中。
• 将GBH框架与现有流体近似方法及标准求解器（如CLASS和CAMB）进行比较。

提出的方法

• 引入双参数展开：中微子分布函数的多极展开（l）和速度加权矩展开（n）。
• 在傅里叶空间中推导广义玻尔兹曼层级（GBH），作为无穷阶常微分方程组，实现动量积分的精确处理。
• 使用更高阶速度加权矩Pn和ωn来表示质量中微子的流体特性，即使在相对论性区域也适用。
• 采用n和l的截断方案，其中nmax和lmax控制精度，该方案基于第一性原理推导。
• 通过一个简易代码验证GBH，将结果与高精度的CLASS求解器在不同中微子质量和波矢下进行比较。
• 在自由流标度以下的小尺度上应用流体近似（GBH+FA），并采用与GBH相同的截断方案。

实验结果

研究问题

• RQ1 在宇宙学微扰理论中，中微子分布函数中的动量依赖性能否被精确积分？
• RQ2 由此产生的双参数层级（l, n）是否能产生与标准玻尔兹曼求解器相当的准确中微子转移函数？
• RQ3 GBH与流体近似方法及成熟代码（如CLASS和CAMB）相比表现如何？
• RQ4 n和l的截断对GBH在不同红移和尺度上的精度有何影响？
• RQ5 GBH能否在现代玻尔兹曼求解器中高效实现，同时保持高精度？

主要发现

• GBH在所有红移和尺度（直至自由流标度）上，与CLASS相比，中微子密度扰动转移函数的精度优于<0.5%。
• 该方法在求解常微分方程组后，完全消除了对数值动量积分的需求，显著提升了计算效率。
• 当在小尺度上结合流体近似（GBH+FA）时，该方法在所有尺度上均达到CLASS默认精度（CLASS-DPS）。
• CAMB在默认精度（CAMB-DPS）下在小尺度上优于GBH+FA和CLASS-DPS，表明GBH流体近似仍有优化空间。
• GBH截断方案的精度与CLASS的流体近似相当，尽管CLASS经过优化的方案整体表现更优。
• 该框架在不同中微子质量（m = 0.1 eV）和波矢（k = 0.008 Mpc⁻¹）下均表现稳健，且在相对论和非相对论区域均保持一致性能。