[论文解读] Testing Higher--Order Lagrangian Perturbation Theory Against Numerical Simulations -- 2. Hierarchical Models

Adrian L. Melott, Thomas Buchert+1 more|ArXiv.org|Apr 8, 1994

Cosmology and Gravitation Theories参考文献 1被引用 14

一句话总结

本文针对负功率谱指数的分层宇宙学模型，将二阶拉格朗日扰动理论（LPT）与N体模拟进行对比。结果表明，二阶LPT在动力学精度方面显著优于泽尔多维奇近似（ZA），尤其在探测引力动力学方面表现优异，且在后期和小尺度上仍具鲁棒性——使其成为建模冷暗物质、热暗物质及混合暗物质宇宙学模型的高计算效率理想选择。

ABSTRACT

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann--Lema\^ıtre cosmogonies investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993), Buchert \& Ehlers (1993), Buchert (1994), Ehlers \& Buchert (1994), is compared with numerical simulations. In this paper we study the dynamics of hierarchical models as a second step. In the first step (Buchert, Melott and Weiß 1994) we analyzed the performance of the Lagrangian schemes for pancake models, i.e., models which initially have a truncated power spectrum. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non--linear regime............We find that for spectra with negative power--index the second--order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low--order statistics like the power--spectrum. In cases with much small--scale power the gain from the higher--order schemes is small, but still measurable. However, in contrast to the results found for pancake models, where the higher--order schemes get worse than TZA at late non--linear stages and on small scales, we here find that the second--order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these

研究动机与目标

• 评估高阶拉格朗日扰动理论在泽尔多维奇近似之外模拟宇宙结构形成的表现。
• 确定在截断功率谱的叠瓦模型中观察到的改进是否可推广至小尺度功率变化的分层模型。
• 优化初始条件滤波（滤波类型与尺度），以最大化与N体模拟的互相关性。
• 评估二阶LPT在深度非线性区域的鲁棒性，特别是针对标准暗物质宇宙学模型。
• 确立二阶LPT在大尺度结构模拟与星系环境建模中的计算与物理优势。

提出的方法

• 采用Buchert（1989–1994）推导的三阶拉格朗日扰动理论解，聚焦一阶与二阶方案。
• 对初始密度场与速度场应用多种滤波类型（高斯型、截断k型），以截断功率谱并优化解析模型性能。
• 针对功率律指数介于-3至+1之间的分层模型执行N体模拟，使用不同滤波尺度的初始条件。
• 通过互相关统计量（包括功率谱与动力学诊断）比较解析LPT预测与N体模拟结果。
• 通过最大化LPT与N体结果之间的互相关系数，优化滤波尺度。
• 分析二阶解在后期与小尺度的行为，与ZA及更高阶方案进行比较。

实验结果

研究问题

• RQ1 在负功率谱指数的分层模型中，二阶拉格朗日扰动理论是否在捕捉动力学方面优于泽尔多维奇近似？
• RQ2 在深度非线性区域，二阶LPT的性能如何演化，特别是在后期与小尺度上？
• RQ3 二阶LPT对ZA的改进在初始条件滤波的不同滤波类型与尺度下是否具有鲁棒性？
• RQ4 在小尺度功率较高的模型中，二阶LPT能否保持精度，即使非线性效应强烈？
• RQ5 二阶方案在标准宇宙学模型（如CDM、HDM与MDM）中在多大程度上仍具有效性？

主要发现

• 二阶拉格朗日扰动理论在动力学精度方面显著优于泽尔多维奇近似，尤其在探测引力动力学的统计量中表现更优。
• 二阶修正带来的改进在深度非线性阶段与小尺度上依然稳健——与叠瓦模型中高阶方案性能退化的情况相反。
• 对于负幂律指数的谱，二阶LPT在动力学统计量中显著优于ZA，而在低阶统计量（如功率谱）中仅略有提升。
• 即使在小尺度功率较高的模型中，二阶修正带来的增益依然可观，尽管相对较小。
• 二阶截断拉格朗日模型计算效率高，CPU开销与泽尔多维奇近似相当。
• 该方法特别适用于建模标准暗物质宇宙学模型（冷暗物质、热暗物质、混合暗物质）中的大尺度结构，可下探至星系群质量尺度（~10^13 M☉），且在截断尺度接近大尺度模拟中的奈奎斯特频率时仍具有效性。