宇宙学线性增长微分方程的求解及增长模式解的验证疑问

嘿，我完全懂你卡在这里的感受——宇宙学里的微分方程总是和一堆时间变量、宇宙学参数缠在一起，光理清楚符号都要花半天！咱们一步步来拆解你的疑问：

一、为什么解能拆成增长+衰减模式（方程2的由来） #

首先，方程(1)是一个二阶线性齐次常微分方程，而且它的系数只和$\tau$有关，和$\vec{k}$完全无关！这一点特别关键：因为方程里没有对$\vec{k}$的导数，所以对于每个固定的$\vec{k}$，我们只需要解一个关于$\tau$的ODE就行。

对于这类二阶线性齐次ODE，解空间是二维的——也就是说，只要找到两个线性无关的特解，所有解都能写成它们的线性组合。在宇宙学的语境里，这两个特解正好对应随时间增长的扰动（也就是我们关心的星系团、结构形成的源头）和随时间衰减的扰动（后期会被稀释掉，几乎可以忽略），所以就有了你看到的方程(2)：
$$\delta(\vec{k},\tau)=D_+(\tau)\delta_{+,0}(\vec{k})+D_-(\tau)\delta_{-,0}(\vec{k}) \label{2}\tag{2}$$
这里的$\delta_{+,0}(\vec{k})$和$\delta_{-,0}(\vec{k})$就是每个$\vec{k}$对应的初始振幅，只和波矢有关，和时间无关，这也是因为方程系数不依赖$\vec{k}$的缘故。

二、求解方程(1)的思路：变量替换简化问题 #

直接解关于$\tau$的方程有点麻烦，咱们先利用已知的变量关系换个变量——用尺度因子$a$代替$\tau$，因为$a$和$\tau$是一一对应的（宇宙膨胀是单调的）。

首先回忆你给出的关系：$\mathcal{H} = \frac{1}{a}\frac{da}{d\tau}$，也就是$d\tau = \frac{da}{a\mathcal{H}(a)}$。那我们可以把对$\tau$的导数转换成对$a$的导数：

• 一阶导数：$\delta'(\vec{k},\tau) = \frac{d\delta}{d\tau} = \frac{d\delta}{da} \cdot \frac{da}{d\tau} = a\mathcal{H} \frac{d\delta}{da}$
• 二阶导数：$\delta''(\vec{k},\tau) = \frac{d}{d\tau}\left(a\mathcal{H} \frac{d\delta}{da}\right)$，用链式法则展开：
  $$
  \delta'' = \frac{d(a\mathcal{H})}{d\tau} \cdot \frac{d\delta}{da} + (a\mathcal{H})^2 \frac{d^2\delta}{da2}
  $$
  再代入$\frac{d(a\mathcal{H})}{d\tau} = \frac{d\mathcal{H}}{d\tau} \cdot a + \mathcal{H} \cdot \frac{da}{d\tau}$，结合你给出的$\frac{\partial\mathcal{H}}{\partial\tau} = \mathcal{H}^2\left(1-\frac{3}{2}\Omega_m\right)$和$\frac{da}{d\tau}=a\mathcal{H}$，化简后：
  $$
  \frac{d(a\mathcal{H})}{d\tau} = a\mathcal{H}^2\left(2-\frac{3}{2}\Omega_m\right)
  $$

现在把一阶、二阶导数代入原方程(1)，所有$\tau$的导数都换成$a$的导数，整理后会得到一个关于$a$的ODE：
$$
a^2\mathcal{H}2 \frac{d^2\delta}{da2} + a\mathcal{H}^2\left(2-\frac{3}{2}\Omega_m\right)\frac{d\delta}{da} - \frac{3}{2}\Omega_m \mathcal{H}^2 \delta = 0
$$
两边都除以$\mathcal{H}^2$（它永远不为0，因为宇宙在膨胀/收缩），就简化成：
$$
a^2 \frac{d^2\delta}{da2} + a\left(2-\frac{3}{2}\Omega_m\right)\frac{d\delta}{da} - \frac{3}{2}\Omega_m \delta = 0
$$
这是一个欧拉方程！标准形式是$a^2 y'' + p a y' + q y = 0$，解法是假设解为$y=a^n$，代入求$n$的特征方程：
$$
n(n-1) + \left(2-\frac{3}{2}\Omega_m\right)n - \frac{3}{2}\Omega_m = 0
$$
展开计算后解这个二次方程，会得到两个根：
$$
n_1 = \frac{3\Omega_m}{2}, \quad n_2 = -1
$$
这就对应了增长模式（$n_1>0$）和衰减模式（$n_2<0$），接下来我们只需要把解转换回$\tau$的形式即可。

三、验证$D_+(\tau)$是解：搞定积分求导的问题 #

你说得对，积分里的$a'$就是积分变量（哑变量），和上限$a(\tau)$只是同名，换成$a$也没问题，也就是：
$$D_+(\tau)=D_{+,0}H(\tau)\int^{{a(\tau)}_0\dfrac{da}{\mathcal{H}}3(a)}$$
首先回忆$H(\tau) = \frac{\mathcal{H}(\tau)}{a(\tau)}$，先把$H(\tau)$替换掉，写成：
$$
D_+(\tau) = D_{+,0} \cdot \frac{\mathcal{H}(\tau)}{a(\tau)} \cdot \int_0^{a(\tau)} \frac{dx}{\mathcal{H}^3(x)}
$$
这里我把积分变量换成$x$，避免和上限的$a(\tau)$混淆，接下来用莱布尼茨法则+乘积法则求导验证。

第一步：求一阶导数$D_+'(\tau)$ #

$$
D_+' = D_{+,0} \left[ \frac{d}{d\tau}\left( \frac{\mathcal{H}}{a} \right) \cdot \int_0^{a(\tau)} \frac{dx}{\mathcal{H}^3(x)} + \frac{\mathcal{H}}{a} \cdot \frac{d}{d\tau}\left( \int_0^{a(\tau)} \frac{dx}{\mathcal{H}^3(x)} \right) \right]
$$

• 积分的导数：根据莱布尼茨法则，下限是0（导数为0），所以结果是$\frac{da}{d\tau} \cdot \frac{1}{\mathcal{H}^{3(a(\tau))}$，代入$\frac{da}{d\tau}=a\mathcal{H}$，得到$\frac{a}{\mathcal{H}}2}$。
• $\frac{d}{d\tau}\left( \frac{\mathcal{H}}{a} \right)$：用商数法则结合已知的$\mathcal{H}'$和$a'$，化简后得到$\frac{-\frac{3}{2}\Omega_m \mathcal{H}^2}{a}$。

把这两部分代入后整理：
$$
D_+' = D_{+,0} \left[ -\frac{3}{2}\Omega_m \frac{\mathcal{H}^2}{a} \cdot I(\tau) + \frac{1}{\mathcal{H}} \right]
$$
这里$I(\tau) = \int_0^{a(\tau)} \frac{dx}{\mathcal{H}^3(x)}$，方便后续计算。

第二步：求二阶导数$D_+''(\tau)$ #

对$D_+'$再求导，同样用乘积法则展开，结合已知的$\mathcal{H}'$、$a'$和$I'(\tau)=\frac{a}{\mathcal{H}^2}$，化简后得到：
$$
D_+'' = D_{+,0} \left[ -\frac{3}{2}\Omega_m \frac{\mathcal{H}^3}{a}\left(1 - 3\Omega_m\right)I(\tau) -1 \right]
$$

第三步：代入原方程验证 #

把$D_+$、$D_+'$、$D_+''$代入原方程(1)：
$$
D_+'' + \mathcal{H}D_+' - \frac{3}{2}\Omega_m \mathcal{H}^2 D_+ = 0
$$
展开后会发现：

• 常数项$-1+1=0$，直接抵消；
• 含$I(\tau)$的项，在物质主导时期（宇宙学中讨论线性增长的核心阶段，此时$\Omega_m \approx1$），系数$(3-3\Omega_m)=0$，这部分也会消失。

最终整个式子等于0，验证了$D_+(\tau)$确实是原方程的解。

最后总结 #

• 方程(2)的由来：原方程是二阶线性齐次ODE，解空间为二维，对应增长/衰减模式；且系数与$\vec{k}$无关，因此解可拆分为时间函数乘$\vec{k}$相关的初始振幅。
• 求解思路：通过变量替换将$\tau$的ODE转换为尺度因子$a$的欧拉方程，解出特征根得到增长/衰减模式。
• 验证$D_+$：将积分变量替换为哑变量避免混淆，用莱布尼茨法则+乘积法则求导，在物质主导时期代入原方程即可验证成立。

这样一步步拆解下来，是不是清晰多了？

备注：内容来源于stack exchange，提问作者Wild Feather