[论文解读] Perturbative Lagrangian Approach to Gravitational Instability

F. R. Bouchet, S. Colombi+2 more|arXiv (Cornell University)|Jun 6, 1994

Cosmology and Gravitation Theories被引用 9

一句话总结

本文提出了一种微扰拉格朗日方法，用于建模宇宙结构中的引力不稳定性，通过展开粒子轨迹而非密度对比度。该方法在二阶时显著优于泽尔多维奇近似，且在高阶下保持更高精度，为从幂律初始条件出发分析红移畸变和大尺度结构形成提供了稳健框架，并通过N体模拟得到验证。

ABSTRACT

This paper deals with the time evolution in the matter era of perturbations in Friedman-Lemaitre models with arbitrary density parameter $Ω$, with either a zero cosmological constant, $Λ= 0$, or with a non-zero cosmological constant in a spatially flat Universe. Unlike the classical Eulerian approach where the density contrast is expanded in a perturbative series, this analysis relies instead on a perturbative expansion of particles trajectories in Lagrangian coordinates. This brings a number of advantages over the classical analysis. In particular, it enables the description of stronger density contrasts. Indeed the linear term is the famous Zel'dovich approximate solution (1970). We present here a systematic and detailed account of this approach. We give analytical results (or fits to numerical results) up to the third order. We then proceed to explore the link between the lagrangian description and statistical measures. We show in particular that Lagrangian perturbation theory provides a natural framework to compute the effect of redshift distortions, using the skewness of the density distribution function as an example. Finally, we show how well the second order theory does as compared to other approximat- ions in the case of spherically symmetric perturbations. We also compare this second order approximation and Zel'dovich solution to N-body simulations in the description of large-scale structure formation starting from a power law (n=-2) power spectrum of Gaussian perturbation. We find that second order theory is both simple and powerful.

研究动机与目标

• 开发一种在物质时代基于拉格朗日轨迹而非欧拉密度对比度的引力不稳定性系统性微扰框架。
• 将泽尔多维奇近似扩展至更高阶（最高至三阶），以提高对非线性结构形成的描述精度。
• 建立一种自然的形式化方法，用于计算统计量（如红移畸变效应）并利用密度分布的偏度进行分析。
• 通过N体模拟验证二阶拉格朗日理论在幂律初始谱（n = -2）下的准确性。
• 评估拉格朗日方法在捕捉线性区域之外强密度对比度区域时的表现。

提出的方法

• 该方法采用拉格朗日坐标中粒子轨迹的微扰展开，一阶项对应泽尔多维奇近似。
• 通过拉格朗日微扰理论形式化系统推导高阶项（最高至三阶），实现对相空间中粒子运动的精确追踪。
• 通过直接演化轨迹避免欧拉微扰理论在高密度对比度下失效的问题。
• 利用拉格朗日框架计算密度场的偏度等统计量，以建模红移畸变。
• 将理论预测与初始幂律谱（n = -2）及不同Ω值下的N体模拟结果进行比较。
• 采用解析拟合与数值解法，验证不同宇宙学模型（Ω为零或非零、空间平坦）中二阶理论的适用性。

实验结果

研究问题

• RQ1 拉格朗日微扰方法在建模非线性结构形成方面，相较于泽尔多维奇近似有何改进？
• RQ2 与其它近似相比，二阶拉格朗日理论在描述球对称扰动时的准确性如何？
• RQ3 拉格朗日框架能否通过偏度等统计量自然地描述红移畸变？
• RQ4 对于幂律初始功率谱（n = -2），二阶拉格朗日理论与N体模拟结果的匹配程度如何？
• RQ5 拉格朗日方法在描述物质时代强密度对比度区域时的有效范围是什么？

主要发现

• 二阶拉格朗日理论在描述球对称扰动演化方面显著优于泽尔多维奇近似，尤其在高密度区域表现更优。
• 拉格朗日方法通过计算密度分布的偏度成功建模了红移畸变，显示出其在统计分析中的实用性。
• 对于幂律初始谱指数n = -2的情况，二阶拉格朗日理论与N体模拟结果高度一致，即使在非线性区域亦然。
• 由于基于轨迹的公式化，该方法在密度对比度高于传统欧拉微扰理论时仍保持有效。
• 该框架提供了最高至三阶的解析表达式，可在无需完整数值模拟的情况下实现精确预测。
• 该理论在不同宇宙学模型中表现良好，包括具有非零宇宙学常数的平坦空间模型。