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

Thomas Buchert, Adrian L. Melott+1 more|arXiv (Cornell University)|Sep 30, 1993

Geological Modeling and Analysis参考文献 1被引用 10

一句话总结

本文在幂律指数 n = -1 的叠瓦模型中，将三阶拉格朗日扰动理论（LPT）与数值模拟进行对比，表明在弱非线性区（σ ≲ 2）中，高阶修正显著提升了对 Zel'dovich 近似（ZA）的精度。研究证实，三阶 LPT 能够降低相位场与密度场的误差，尤其在结合初始功率谱截断时效果更佳，为宇宙大尺度结构形成提供了更精确的解析框架。

ABSTRACT

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasilinear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993a), Buchert \& Ehlers (1993), Buchert (1993b), Ehlers \& Buchert (1993), is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work (Coles \etal 1993, Melott \etal 1993, Melott 1993) the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973, hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of $σ\approx 2$). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-index $n=-1$) using cross-correlation statistics employed in

研究动机与目标

• 评估高阶拉格朗日扰动理论（LPT）在超越 Zel'dovich 近似（ZA）的宇宙大尺度结构形成建模中的表现。
• 评估位移场中的三阶修正是否能提升在幂律指数 n = -1 的叠瓦模型中，弱非线性区（σ ≲ 2）的预测精度。
• 通过互相关统计量，将三阶 LPT 的预测能力与数值模拟进行比较。
• 研究初始功率谱截断对分层聚类模型中解析近似精度的影响。
• 为高阶 LPT 建立基准，作为宇宙模拟中第一阶方法的更优替代方案。

提出的方法

• 为爱因斯坦-德西特尘埃宇宙学构建三阶拉格朗日扰动理论，扩展了先前的一阶与二阶解。
• 将该理论应用于具有截断功率谱（n = -1）的单一叠瓦模型，以抑制小尺度噪声并改善收敛性。
• 执行 N 体数值模拟，生成用于比较的参考密度场与速度场。
• 利用解析预测与模拟输出之间的互相关统计量，量化不同尺度下的精度。
• 将三阶 LPT 结果与 Zel'dovich 近似（ZA）及其他先前研究中的解析近似进行比较。
• 评估相位场与密度场误差作为线性均方根密度对比度 σ 的函数，重点关注 σ ≲ 2 的范围。

实验结果

研究问题

• RQ1 三阶拉格朗日扰动理论在预测叠瓦模型中密度场与速度场时，相较于 Zel'dovich 近似有何改进？
• RQ2 初始功率谱截断在弱非线性区中，对高阶 LPT 精度的提升程度如何？
• RQ3 与一阶方法相比，使用三阶 LPT 在相位场与密度场重构中可实现多大程度的误差降低？
• RQ4 解析预测与数值模拟之间的互相关统计量如何验证三阶 LPT 的性能？
• RQ5 在 ZA 开始失效的大尺度区域（弱非线性区）中，高阶 LPT 是否仍能保持精度？

主要发现

• 在弱非线性区（σ ≲ 2）中，三阶拉格朗日扰动理论相较于 Zel'dovich 近似，显著降低了相位场与密度场的误差。
• 引入三阶修正后，解析预测与数值模拟之间的互相关性得到提升，尤其在高密度对比度区域更为明显。
• 初始功率谱截断提升了 ZA 与高阶 LPT 的精度，但其对三阶 LPT 的提升更为显著，原因在于对非线性模耦合的更好控制。
• 本研究证实，三阶 LPT 在将初始条件映射到最终密度结构方面，优于一阶方法。
• 该方法在密度场的相位与振幅上均与模拟结果更加一致，尤其在叠瓦与丝状结构的形成中表现更优。
• 结果验证了高阶 LPT 作为建模宇宙大尺度结构的稳健解析工具，尤其在弱非线性区中具有显著优势。