[论文解读] Dynamics of gravitational clustering II. Steepest-descent method for the quasi-linear regime

Patrick Valageas|ArXiv.org|Jul 6, 2001

Cosmology and Gravitation Theories参考文献 15被引用 60

一句话总结

本文提出一种非微扰的最速下降法，用于计算引力聚集准线性 regime 中密度对比 $\mathcal{P}(\delta_R)$ 的概率密度函数（PDF）。通过在 $\sigma \to 0$ 极限下处理作用量的路径积分，该方法得到渐近精确的结果，修正了以往对 $n < 0$ 功率谱在高密度尾部的错误，并提供了一种无需依赖流体近似、比微扰流体动力学方法更严格且更直观的非微扰框架。

ABSTRACT

We develop a non-perturbative method to derive the probability distribution $P(δ_R)$ of the density contrast within spherical cells in the quasi-linear regime. Indeed, since this corresponds to a rare-event limit a steepest-descent approximation can yield asymptotically exact results. We check that this is the case for Gaussian initial density fluctuations, where we recover most of the results obtained by perturbative methods from a hydrodynamical description. Moreover, we correct an error which was introduced in previous works for the high-density tail of the pdf. This feature, which appears for power-spectra with a slope $n&lt;0$, points out the limitations of perturbative approaches which cannot describe the pdf $P(δ_R)$ for $δ_R \ga 3$ even in the limit $σ o 0$. This break-up does not involve shell-crossing and it is naturally explained within our framework. Thus, our approach provides a rigorous treatment of the quasi-linear regime, which does not rely on the hydrodynamical approximation for the equations of motion. Besides, it is actually simpler and more intuitive than previous methods. Our approach can also be applied to non-Gaussian initial conditions.

研究动机与目标

• 在不依赖微扰流体动力学近似的情况下，推导球形区域内密度对比 $\mathcal{P}(\delta_R)$ 的概率密度函数，适用于准线性 regime。
• 在 $\sigma \to 0$ 极限下，为计算 PDF 提供一种非微扰、严格的框架，此时微扰理论发散。
• 修正先前未被考虑的、针对 $n < 0$ 功率谱在 PDF 高密度尾部的错误，该问题在微扰方法中即使在 $\sigma \to 0$ 极限下也难以准确描述。
• 将该方法扩展至非高斯初始条件，如合作者工作所展示。

提出的方法

• 对密度对比的生成泛函使用最速下降近似，该近似在 $\sigma \to 0$ 时渐近成立。
• 将该方法应用于由无碰撞 Boltzmann 方程和泊松方程在共动坐标系中导出的作用量泛函。
• 通过鞍点近似构建相关生成函数 $\overline{\psi}(y)$，并引入非线性映射 $\delta_R = \mathcal{G}(\tau)$。
• 引入稀释因子 $1/(1 + \delta_R)$，以考虑共动坐标系中高密度区域的收缩，从而校正背景膨胀的影响。
• 通过最速下降法计算 $\overline{\psi}(y)$ 的逆拉普拉斯变换，得到 $\mathcal{P}(\delta_R)$ 的闭式表达式。
• 通过与球形坍缩模型对比验证结果，发现其与形式 $\mathcal{P}(\delta_R) \propto \frac{1}{1 + \delta_R} \frac{d\nu}{d\delta_R} e^{-\nu^2/2}$ 完全一致。

实验结果

研究问题

• RQ1 如何在不依赖微扰流体动力学近似的情况下，严格计算准线性 regime 中密度对比 $\mathcal{P}(\delta_R)$ 的概率密度函数？
• RQ2 为何即使在 $\sigma \to 0$ 极限下，微扰方法仍无法准确描述 $\delta_R \gtrsim 3$ 时的 $\mathcal{P}(\delta_R)$ 高密度尾部，特别是对 $n < 0$ 功率谱而言？
• RQ3 基于最速下降近似的非微扰方法是否能在 $\sigma \to 0$ 极限下对高斯初始条件给出精确结果？
• RQ4 背景宇宙的膨胀如何影响 PDF 的高密度尾部？这一效应如何在路径积分表述中一致地纳入？
• RQ5 该方法能否推广至非高斯初始条件，从而克服标准微扰方法的局限性？

主要发现

• 最速下降法在 $\sigma \to 0$ 极限下对 $\mathcal{P}(\delta_R)$ 给出渐近精确结果，为先前的微扰结果提供了严格的理论基础。
• 该方法修正了以往未被考虑的、针对 $n < 0$ 功率谱在 PDF 高密度尾部的错误，该问题在微扰方法中即使在 $\sigma \to 0$ 极限下也难以准确描述。
• 所导出的 PDF 形式为 $\mathcal{P}(\delta_R) \simeq \frac{1}{\sqrt{2\pi}\sigma} \frac{1}{1 + \delta_R} \frac{1}{|\mathcal{G}'(\tau)|} e^{-\tau^2/(2\sigma^2)}$，与球形坍缩模型完全一致。
• 引入稀释因子 $1/(1 + \delta_R)$ 可有效描述由于背景膨胀导致的共动坐标系中高密度区域收缩，这一物理效应在早期微扰处理中被忽略。
• 该方法比以往方法更简洁、更直观，且不依赖流体近似，因此适用于非高斯初始条件。
• PDF 中的指数截断是精确的，而前因子则通过最速下降法一致导出，验证了 Valageas (1998) 提出的球形模型方法。