[論文レビュー] Generalized Boltzmann hierarchy for massive neutrinos in cosmology

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

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

ひとこと要約

本稿では、運動量依存性を高次速度加重積分を用いて統合することで、宇宙論的質量ニュートリノのモデル化のための新しい手法、一般化されたボルツマン階層（GBH）を導入する。この手法により、問題は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）を用いた2パラメータ展開を導入する。
• 運動量の正確な統合を伴う、フーリエ空間における無限個の常微分方程式系として一般化されたボルツマン階層（GBH）を導出する。
• 高次速度加重モーメントPnおよびωnを用いて、相対論的領域でも流体的性質を表現する。
• nおよびlの切り捨てスキーム（nmaxおよびlmax）を用い、第一原理から導出された精度制御を行う。
• トイラコードを用いてGBHを検証し、さまざまなニュートリノ質量および波数に対して、高精度なCLASSソルバーと比較する。
• 自由拡散スケール以下の小スケールでは、GBHと同一の切り捨てスキームを用いた流体近似（GBH+FA）を適用する。

実験結果

リサーチクエスチョン

• RQ1 宇宙論的摂動理論において、ニュートリノ分布関数の運動量依存性を正確に統合し得るか？
• RQ2 得られた2パラメータ階層（l, n）は、標準的なボルツマンソルバーと同等の精度のニュートリノ伝達関数を生成するか？
• RQ3 GBHは流体近似および確立されたコード（CLASSやCAMB）と比較して、どのように性能を発揮するか？
• RQ4 nおよびlの切り捨てが、さまざまな赤方偏移およびスケールでGBHの精度に与える影響は何か？
• RQ5 GBHは最新のボルツマンソルバーに効率的に実装可能であり、高精度を維持できるか？

主な発見

• GBHは、すべての赤方偏移および自由拡散スケールまでのすべてのスケールで、CLASSと比較してニュートリノ密度対比伝達関数において0.5%未塔の精度を達成する。
• ODE系を解いた後、数値的運動量積分の必要性が完全に排除され、計算効率が顕著に向上する。
• 小スケールで流体近似（GBH+FA）を組み合わせた場合、GBHはすべてのスケールでCLASSのデフォルト精度（CLASS-DPS）と一致する。
• CAMBのデフォルト精度（CAMB-DPS）は、小スケールでGBH+FAおよびCLASS-DPSを上回り、GBHの流体近似に改善の余地があることを示唆している。
• GBHの切り捨てスキームは、CLASSの流体近似と同等の精度を示すが、CLASSの最適化されたスキームは全体的に優れている。
• フレームワークは、ニュートリノ質量（m = 0.1 eV）および波数（k = 0.008 Mpc⁻¹）の変動に対しても安定しており、相対論的および非相対論的領域で一貫した性能を示す。