0617-宇宙学扰动理论

暴涨与量子泡沫

在暴涨宇宙学中， $\delta N$ 是一个非常核心的概念，尤其在处理 大尺度结构的起源 和 非高斯性计算 时常常被用到。我们来从多个角度解释它的含义与作用。

---

🔹 一句话定义：

$$\boxed{\delta N = \mathcal{R}}$$

** $\delta N$ 是不同宇宙区域之间的暴涨 e-folding 数（膨胀倍数）的差异，直接等价于曲率扰动 $\mathcal{R}$ 。**

---

📚 全称： $\delta N$ 公式（ $\delta N$ formalism）

这是由 Sasaki & Stewart（1995）提出的一种用于计算暴涨产生扰动的方法，它基于这样的想法：

“宇宙中不同区域因标量场涨落导致暴涨时间不同 → 最终空间曲率不同。”

---

🧠 物理解读：

✅ 暴涨持续时间的空间变化

• 总的膨胀量（e-folding number）定义为：

  $$N = \int_{t_*}^{t_e} H dt$$
• 如果场 $\varphi$ 有涨落，那么某些区域比其他区域多膨胀一点：

  $$N(\mathbf{x}) = \bar{N} + \delta N(\mathbf{x})$$
• 这个 $\delta N$ 就代表这个位置与平均膨胀量的偏差。

---

📐 几何解释：

在“共动切片（初始时）”与“均匀密度切片（末态）”之间，某一点的膨胀量是：

$$\delta N = N(\varphi + \delta\varphi) - N(\varphi)$$

不同点由于初始 $\delta\varphi$ 不同 → 导致 $\delta N$ 不同 → 最终造成空间曲率扰动。

---

📈 与曲率扰动 $\mathcal{R}$ 的关系：

在超视界尺度上，有非常重要的等式：

$$\boxed{\mathcal{R} = \delta N}$$

这表示大尺度曲率扰动直接由暴涨期间每个区域的膨胀时间差造成。

---

✍️ 实用表达式

如果是单场慢滚暴涨：

$$\delta N \approx \frac{dN}{d\varphi} \delta\varphi \approx \frac{H}{\dot\varphi} \delta\varphi$$

代入 $\delta\varphi \sim \frac{H}{2\pi}$ ，得到经典结果：

$$\mathcal{R} \sim \delta N \sim \frac{H^2}{2\pi \dot\varphi}$$

---

🚀 应用：

1. 计算标量扰动功率谱：

   $$P_\mathcal{R}(k) = \left( \frac{H^2}{2\pi \dot\varphi} \right)^2$$
2. 多场暴涨的扰动传播与非高斯性估计：

   $$\zeta = N_{,\phi^i} \delta\phi^i + \frac{1}{2} N_{,\phi^i\phi^j} \delta\phi^i \delta\phi^j + \cdots$$
3. 非高斯性计算中的：

   $$f_{\text{NL}} \sim \frac{N_{,\phi^i\phi^j}}{(N_{,\phi^i} N_{,\phi^i})}$$

---

📌 小结表格：

项目	含义
$\delta N$	某点相对于背景膨胀数量的变化
本质	膨胀时间不同导致的空间拉伸差异
意义	与曲率扰动 $\mathcal{R}$ 等价
用途	计算功率谱、非高斯性、多场扰动演化等

---

如果你希望我展示如何从 $\delta N$ 推导出 $P_\mathcal{R}(k)$ 、或者在多场暴涨下如何展开 Taylor 展式，也可以继续。需要吗？

Response:

非常好，我们就从最根本的起点，系统性地讲清楚整个 $\delta N$ 理论，从最初动机到最终物理应用，尤其是暴涨宇宙学中如何通过它来描述扰动的来源和演化。

---

🧭 0. 背景问题：结构从哪来？

我们知道，宇宙中的星系、星系团等大尺度结构是由早期宇宙中的微小密度扰动生长而来。那么这些初始扰动是怎么来的？

暴涨理论告诉我们：

这些扰动起源于暴涨时期量子涨落，最初是微观的 $\delta\varphi$ ，然后被膨胀拉伸到宇宙尺度，并最终变成我们观测到的宇宙结构种子。

---

🚀 1. 什么是 e-folding 数 $N$ ？

膨胀时期宇宙的膨胀程度通常通过 e-folding 数 $N$ 表示：

$$N(t_1, t_2) = \ln \frac{a(t_2)}{a(t_1)} = \int_{t_1}^{t_2} H(t) dt$$

• 表示从时间 $t_1$ 到 $t_2$ ，宇宙膨胀了多少个“e倍”。

• 在慢滚暴涨中， $H$ 几乎常数，所以 $N$ 和暴涨持续时间成正比。

---

🎯 2. 暴涨的关键思想：不同位置膨胀时间不一样

暴涨场（inflaton） $\varphi$ 有涨落 $\delta\varphi(\mathbf{x})$ ，意味着不同位置的初始条件略有不同。

→ 这导致：

• 不同空间点 $\mathbf{x}$ 的暴涨时间不同

• 每个 patch 的 e-folding 数变成空间依赖量： $N(\mathbf{x}) = \bar{N} + \delta N(\mathbf{x})$

这就是：

$$\boxed{\delta N(\mathbf{x}) \equiv N(\mathbf{x}) - \bar{N}}$$

---

🌀 3. 从 $\delta N$ 到曲率扰动 $\mathcal{R}$

空间膨胀得多的区域在暴涨后密度就更稀疏，对应负曲率扰动；反之密度更高，对应正扰动。

于是我们有：

$$\boxed{\mathcal{R}(\mathbf{x}) = \delta N(\mathbf{x})}$$

这就是所谓的 ** $\delta N$ formalism** 的核心等式。它告诉你：曲率扰动其实就是每块区域膨胀时间的差异。

---

🧪 4. 如何实际计算 $\delta N$ ？

我们从暴涨动力学出发：

$$N(\varphi) = \int_{\varphi_{\text{end}}}^{\varphi} \frac{H}{\dot{\varphi}} d\varphi$$

那么：

$$\delta N = \frac{dN}{d\varphi} \delta\varphi = \frac{H}{\dot{\varphi}} \delta\varphi$$

而 $\delta\varphi \sim \frac{H}{2\pi}$ 是由量子涨落导出的结果。

代入后得到：

$$\mathcal{R} = \delta N = \frac{H}{\dot{\varphi}} \cdot \frac{H}{2\pi} = \boxed{\frac{H^2}{2\pi\dot{\varphi}}}$$

这就是暴涨生成的标量扰动大小。

---

📈 5. 应用 1：功率谱

标量扰动的功率谱就是：

$$\boxed{P_\mathcal{R}(k) = \left( \frac{H^2}{2\pi\dot{\varphi}} \right)^2 \Bigg|_{k=aH}}$$

这是计算 CMB 初始扰动的基石。

---

🌌 6. 应用 2：非高斯性计算

对于多场暴涨，我们将 $\delta N$ 在多个场上展开：

$$\mathcal{R} = \delta N = N_{,i} \delta\phi^i + \frac{1}{2} N_{,ij} \delta\phi^i \delta\phi^j + \cdots$$

于是可以计算非高斯性参数 $f_{\text{NL}}$ ，例如：

$$f_{\text{NL}} \sim \frac{N_{,ij} N_{,i} N_{,j}}{(N_{,k}N_{,k})^2}$$

这正是现代精密宇宙学分析 Planck 数据时的标准计算方式。

---

🧩 7. 几何解释补充

你可以从两个切片（slices）之间理解：

• 初始是“共动切片”： $\delta\varphi \neq 0$

• 结束时是“均匀密度切片”： $\delta\rho = 0$

$\delta N$ 就是这两种切片之间的“时间间隔”的扰动。

这也是为什么有时写成：

$$\mathcal{R} = N(\rho = \text{const}) - N(\varphi = \text{const})$$

---

✅ 总结表格：

概念	内容
$N$	每个点的暴涨膨胀量
$\delta N$	每个点相对平均膨胀的差异
$\delta N = \mathcal{R}$	与空间曲率扰动等价（超视界上）
$\delta N = \frac{H}{\dot{\varphi}}\delta\varphi$	实际计算公式（单场）
应用	功率谱计算，非高斯性分析，多场演化追踪

---

如果你接下来想从 $\delta N$ 扩展到：

• 非共动切片如何定义扰动？

• 张量扰动是否也能用 $\delta N$ 方法？

• 如何从 $\delta N$ 连接到 Mukhanov-Sasaki 方程？

我可以继续讲下去。你想深入哪个部分？

Response:

Mukhanov–Sasaki 方程是描述暴涨时期标量扰动演化的基本方程，它是早期宇宙微扰理论的核心内容之一。这个方程的解给出了我们今天在宇宙微波背景（CMB）中看到的温度扰动的起源。

我们从头解释：

---

🧭 背景动机：暴涨扰动怎么演化？

我们关心的是**暴涨场 $\varphi$ 的扰动 $\delta\varphi$ ** 如何在弯曲时空中传播，并最终变成我们看到的密度扰动或曲率扰动 $\mathcal{R}$ 。

这需要把：

• **场论扰动 $\delta\varphi$ **、

• 时空扰动（引力） 统一考虑 → 这就是Mukhanov–Sasaki变量登场的地方。

---

🧮 Step 1：定义 Mukhanov–Sasaki 变量

暴涨场扰动 $\delta\varphi$ 和时空扰动（如 $\Phi, \mathcal{R}$ ）纠缠在一起，不好分开分析。

于是引入一个新的变量：

$$\boxed{v \equiv a \left( \delta\varphi + \frac{\varphi'}{\mathcal{H}} \Phi \right)}$$

或者在共动切片下（ $\delta\varphi = 0$ ），也可以写成：

$$\boxed{v = z \mathcal{R}} \quad \text{其中 } z \equiv \frac{a\varphi'}{\mathcal{H}}$$

这个变量 $v$ 就是 Mukhanov–Sasaki 变量，它是规范不变的，并且可以像自由场一样写拉氏量和方程。

---

🧪 Step 2：写出方程本身

Mukhanov–Sasaki 方程的形式如下：

$$\boxed{v_k'' + \left(k^2 - \frac{z''}{z} \right) v_k = 0}$$

• $v_k$ ：是 Fourier 模式的扰动变量

• $k$ ：是波数

• $z''/z$ ：是膨胀背景引入的“有效质量项”

它形式上就像一个变质量的简谐振子，这是其美妙之处：时空膨胀 + 标量扰动 → 凝聚成一个单一变量的方程。

---

🔬 Step 3：物理含义

这个方程控制着标量扰动从量子涨落→冻结→经典扰动的全过程：

• 在 $k \gg aH$ （亚视界）：像自由波一样震荡

• 在 $k \ll aH$ （超视界）：解变为常数 + 渐弱模式，扰动“冻结”

---

📈 Step 4：求解与功率谱

在暴涨时期，我们可以近似 $z''/z \approx 2/\eta^2$ （ $\eta$ 是共形时间），解为：

$$v_k(\eta) \approx \frac{1}{\sqrt{2k}} e^{-ik\eta} \left(1 - \frac{i}{k\eta} \right)$$

然后代入 $\mathcal{R} = v/z$ ，计算得到功率谱：

$$\boxed{P_\mathcal{R}(k) = \left( \frac{H^2}{2\pi \dot\varphi} \right)^2}$$

这正是我们在 $\delta N$ 方法中看到的结果，说明两种方法是一致的。

---

🧩 Mukhanov–Sasaki 方程总结表格

项目	内容
变量定义	$v = z\mathcal{R} = a(\delta\varphi + \frac{\varphi'}{\mathcal{H}}\Phi)$
方程形式	$v_k'' + (k^2 - z''/z)v_k = 0$
适用阶段	暴涨期间的标量扰动演化
特点	可量子化、规范不变、适用于功率谱推导
结果	解释初始扰动幅度、谱指数、非高斯性等

---

✅ Mukhanov–Sasaki 方程 vs $\delta N$ 方法

方法	优点	缺点
Mukhanov–Sasaki 方程	精确，适用于任意模式，推导完整动力学	数学更复杂，不直观
$\delta N$ 方法	简单直观，易计算非高斯性	仅适用于超视界、慢滚近似

---

是否希望我为你推导 Mukhanov–Sasaki 方程的来源（从 action 到 Euler–Lagrange）？还是更想看张量扰动下的类似结构？

Response:

非常好的问题：“为什么扰动变量 $v$ 满足 Mukhanov–Sasaki 方程？”——

我们现在从变分原理出发，从头推导这个方程，解释它为何成立。我们只做线性（一级）扰动分析，处理的是暴涨时的标量扰动。我们将分三步来讲清楚：

---

🧱 Step 1：从 action 出发 —— 二阶扰动拉氏

我们从 Einstein–Hilbert + 标量场的总 action 出发，对背景膨胀宇宙中的小扰动展开到二阶。

Mukhanov 在 1985 年计算了这个 action，并发现：

$$\boxed{ S^{(2)} = \frac{1}{2} \int d\eta\, d^3x\, \left[ (v')^2 - (\nabla v)^2 + \frac{z''}{z} v^2 \right] }$$

这是扰动变量 $v(\eta, \mathbf{x})$ 的有效作用量。这个变量是通过选取规范不变组合得到的，称为 Mukhanov 变量：

$$v = a \left(\delta\varphi + \frac{\varphi'}{\mathcal{H}} \Phi\right)$$

• $a$ 是膨胀因子

• $\Phi$ 是 Newtonian gauge 下的引力势扰动

• $\varphi'$ 是背景场的共形时间导数

• $\mathcal{H} = a'/a$

此时我们已经把引力+标量场的系统简化为一个自由场在时变背景下的行为，惊人的简洁。

---

🧮 Step 2：变分导出运动方程

从该 action 对 $v$ 做变分，得到 Euler–Lagrange 方程：

$$\frac{\partial}{\partial \eta} \left( \frac{\partial \mathcal{L}}{\partial v'} \right) + \nabla^2 v - \frac{z''}{z} v = 0$$

得到：

$$\boxed{ v'' - \nabla^2 v - \frac{z''}{z} v = 0 } \quad \Rightarrow \quad \boxed{ v_k'' + \left( k^2 - \frac{z''}{z} \right) v_k = 0 }$$

这是我们要的 Mukhanov–Sasaki 方程！

---

🧠 Step 3：物理理解

• $v$ ：代表一个有效自由场，既包含了暴涨场的扰动，也包含了引力扰动（因为我们使用的是 gauge-invariant 的组合）

• $z''/z$ ：来自时空膨胀背景的“有效势能项”，控制扰动增长/冻结

• $k^2$ ：空间传播项，控制扰动震荡

这就好比一个“变质量”简谐振子。扰动如何演化完全取决于 $k^2$ 与 $z''/z$ 的对比：

情况	描述	物理行为
$k^2 \gg z''/z$	模式在视界内	震荡，如平面波（亚视界）
$k^2 \ll z''/z$	模式超视界	解冻结，不再变化（扰动定格）

---

🧩 所以 Mukhanov–Sasaki 方程为什么成立？

✅ 因为它是从暴涨系统的总 action 出发，对一阶扰动做精确变分得到的 ✅ 它描述的是 gauge-invariant 变量的演化，因此物理、规范无关、可以量子化 ✅ 它的解决定了暴涨扰动的起源、功率谱和 CMB 各向异性

---

🧪 扩展阅读（可选）

如果你想进一步深入推导，我可以展示：

1. 从完整拉氏：

   $$S = \int d^4x \sqrt{-g} \left[ \frac{1}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \varphi \partial_\nu \varphi - V(\varphi) \right]$$

   展开到二阶并做变量替换。

2. 显式写出背景 Friedmann 方程对 $z = a\varphi'/\mathcal{H}$ 的控制关系。

---

是否希望我继续从完整 action 展开、推导这个二阶 action？这涉及 ADM 分解和大量微分计算，但非常有价值。

Response:

非常好，我们现在从完整 action 出发一步步推导 Mukhanov–Sasaki 方程，讲清楚为什么时空+标量扰动可以被写成一个单变量 $v$ 的二阶方程。这是现代宇宙扰动理论的基石。

---

✅ 总路线图

我们要从下面的引力+标量场系统开始：

$$S = \int d^4x\, \sqrt{-g} \left[ \frac{1}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \varphi \partial_\nu \varphi - V(\varphi) \right]$$

目标是推导出对某个变量 $v$ 的有效二阶 action：

$$S^{(2)} = \frac{1}{2} \int d\eta\, d^3x\, \left[ (v')^2 - (\nabla v)^2 + \frac{z''}{z} v^2 \right] \quad \text{with } v = z\mathcal{R}$$

---

1️⃣ 背景设置：FLRW时空 + 扰动

我们考虑一个平坦FLRW背景加上标量扰动和时空扰动：

$$ds^2 = -(1+2\Phi)\,dt^2 + a^2(t)(1-2\Psi)\,\delta_{ij}dx^i dx^j$$

标量场：

$$\varphi(\vec{x}, t) = \bar{\varphi}(t) + \delta\varphi(\vec{x}, t)$$

其中：

• $\Phi, \Psi$ ：是度规扰动（标量型）

• $\delta\varphi$ ：是标量场的扰动

• 全部只保留到一阶微扰

我们将对整个系统在二阶展开，构造一个 action。

---

2️⃣ 使用 ADM 分解法

这是做引力系统扰动最标准的方式。我们把时空度规写成：

$$ds^2 = -N^2 dt^2 + h_{ij}(dx^i + N^i dt)(dx^j + N^j dt)$$

其中：

• $N$ ：lapse function（测量时间）

• $N^i$ ：shift vector（空间–时间混合扰动）

• $h_{ij}$ ：三维空间 metric

然后我们将标量场和引力的 action 插入这个分解形式，对变量 $N, N^i, h_{ij}, \delta\varphi$ 展开。

关键观察：

• $N$ , $N^i$ 是无动力变量（non-dynamical），可以通过 constraint equation 解掉

• 只留下1个物理自由度

---

3️⃣ 动力学变量提取

经过完整展开后（见 Mukhanov 1988, Maldacena 2002），你会发现：

• 引力和场的耦合可以归结为一个自由变量（自由度数 = 1）

• 最佳选择是一个规范不变变量

定义：

$$\boxed{ v \equiv a \left( \delta\varphi + \frac{\varphi'}{\mathcal{H}} \Phi \right) } \quad \text{或} \quad \boxed{v = z\mathcal{R}}, \quad z \equiv \frac{a\varphi'}{\mathcal{H}}$$

这个变量 $v$ 是一个规范不变的组合，代表暴涨期间真正的自由扰动。

---

4️⃣ 得到二阶 action

使用上述变量，完整的二阶 action 精确地写成：

$$\boxed{ S^{(2)} = \frac{1}{2} \int d\eta\, d^3x\, \left[ (v')^2 - (\nabla v)^2 + \frac{z''}{z} v^2 \right] }$$

这是非常了不起的结果：

• 整个暴涨扰动系统（引力 + 场）被**统一描述为一个“准自由”场 $v$ ** 在“背景势” $z''/z$ 中传播。

• 这个 action 是规范不变的、可以直接进行量子化。

---

5️⃣ 从 action 导出 Mukhanov–Sasaki 方程

从变分原理得到 Euler–Lagrange 方程：

$$\frac{\delta S}{\delta v} = 0 \Rightarrow v'' - \nabla^2 v - \frac{z''}{z}v = 0$$

在动量空间：

$$\boxed{ v_k'' + \left(k^2 - \frac{z''}{z} \right)v_k = 0 }$$

这就是 Mukhanov–Sasaki 方程。

---

🔬 为什么这是最终方程？

• 因为在处理引力和暴涨场扰动时，你只需要考虑 一个自由度（单场暴涨）

• Mukhanov–Sasaki 方程提供了该自由度完整的演化方程

• 解这个方程就能得到功率谱、谱指数、甚至非高斯性等所有重要结果

---

📌 总结

步骤	关键点
Start	从引力+标量场 action 出发
ADM	拆解度规 → 识别非动力变量
Perturb	做二阶扰动展开
Canonical variable	定义 $v = z\mathcal{R}$ ，简化结构
Final action	得到自由场 action + 背景势 $z''/z$
Euler–Lagrange	得到 Mukhanov–Sasaki 方程

---

是否希望我配图解释这个推导结构？或展示在慢滚条件下 $z''/z$ 的近似形式和解的行为？

Response:

非常好，我们现在来解释：Mukhanov–Sasaki 方程究竟讲了什么物理内容？

这不是一个普通的数学方程，而是揭示宇宙结构从量子涨落生成的全过程的核心方程。在现代宇宙学中，它与 Friedmann 方程、Einstein 方程并列为最根本的几条之一。

---

📌 方程本身：

$$\boxed{ v_k'' + \left(k^2 - \frac{z''}{z} \right) v_k = 0 } \quad \text{with} \quad v = z\mathcal{R}, \quad z = \frac{a\varphi'}{\mathcal{H}}$$

这是 Fourier 模式下的 Mukhanov–Sasaki 方程。

---

🌌 它在讲什么？

一句话总结：

这个方程描述了宇宙暴涨期间，量子涨落如何演化成经典的大尺度宇宙结构种子。

我们把这句话拆开细讲：

---

🧠 解构这个方程

1. 主角： $v_k$ 是什么？

• $v_k$ 是 Mukhanov–Sasaki 变量在动量空间的表示

• 它本质上是：

  $$v = z\mathcal{R}$$

  是 **共动曲率扰动 $\mathcal{R}$ ** 的“放大版本”，更方便用于计算和量子化

---

2. 它是一个什么类型的方程？

结构上它像一个“变质量的谐振子”：

$$v_k'' + [k^2 - m^2(\eta)]v_k = 0 \quad \text{with } m^2(\eta) = \frac{z''}{z}$$

• $k^2$ ：对应扰动的空间传播（震荡）

• $z''/z$ ：由背景宇宙决定的“时变势能”或“质量项”

---

🔍 演化行为解析：两种尺度

🟦 1. **亚视界阶段 $k^2 \gg z''/z$ **

• 方程简化为：

  $$v_k'' + k^2 v_k \approx 0 \Rightarrow v_k \sim e^{\pm ik\eta}$$
• 解是震荡的平面波 → 类似真空波动

• 表示扰动还没“冻结”，像自由粒子一样在振动

---

🟥 2. **超视界阶段 $k^2 \ll z''/z$ **

• 方程变成：

  $$v_k'' - \frac{z''}{z}v_k \approx 0$$
• 解变为一个增长模式和一个衰减模式：

  $$v_k \sim A_k z + B_k z \int \frac{1}{z^2} d\eta \Rightarrow \mathcal{R}_k = \frac{v_k}{z} \to \text{const}$$
• 表示扰动冻结在一个常数值上

这解释了为什么大尺度上的曲率扰动 $\mathcal{R}$ 是时间不变的：它在模式出视界时“冻结”下来。

---

📈 功能和作用

功能	说明
✨ 预测初始扰动功率谱	$\mathcal{R} = v/z \Rightarrow P_\mathcal{R}(k) = \left( \frac{H^2}{2\pi\dot\varphi} \right)^2$
🧮 模拟扰动演化	可数值解算，处理非慢滚等情况
🧪 统一引力和场论扰动	以一个变量描述整个系统演化
🎯 描述量子→经典过渡	模式出视界后震荡终止，扰动冻结，变成经典曲率种子

---

🧩 总结一句话

Mukhanov–Sasaki 方程揭示了暴涨扰动从量子起源、震荡传播、穿越视界、最终冻结成大尺度结构种子的完整物理过程。

---

如果你希望，我可以：

• 给你画出它在不同 $k$ 下的解的行为图

• 讲解在慢滚暴涨时 $z''/z$ 的近似

• 对比张量扰动（引力波）的对应方程结构

你想继续哪部分？

---

Powered by ChatGPT Exporter

非常好，我们现在重新深入展开第一部分：暴涨背景与量子扰动的起源，和第二部分：Mukhanov–Sasaki 方程的推导与含义，讲清楚每个变量的物理意义、推导路径和逻辑联系。

---

🧮 第一部分：暴涨与量子扰动的起源

---

✅ 1. 暴涨是什么？

暴涨（Inflation）是早期宇宙的一段极快速指数膨胀过程：

$$a(t) \propto e^{H t}, \quad H = \text{ nearly constant}$$

其目的是解释：

• 宇宙的平坦性问题（Flatness Problem）

• 宇宙各向同性问题（Horizon Problem）

• 起源于量子涨落的宇宙结构种子（Structure Formation）

---

✅ 2. 暴涨由谁驱动？

通过引入一个**标量场 $\phi(t)$ **（称为 inflaton），能量密度以势能为主：

$$\rho \approx V(\phi), \quad \text{if } \frac{1}{2}\dot\phi^2 \ll V(\phi)$$

---

✅ 3. 背景演化方程（慢滚条件）

Friedmann 方程：

$$H^2 = \frac{1}{3 M_{\rm Pl}^2} \left( \frac{1}{2} \dot\phi^2 + V(\phi) \right)$$

Klein–Gordon 方程：

慢滚条件：

$$\epsilon = -\frac{\dot H}{H^2} \ll 1, \quad \eta = \frac{\ddot\phi}{H \dot\phi} \ll 1$$

在这些条件下，宇宙接近de Sitter状态，但有微小演化，才允许谱不完全平坦（非严格scale-invariant）。

---

✅ 4. 量子扰动从哪来？

暴涨时，空间中充满真空涨落（量子起伏）：

• 标量场涨落： $\delta \phi(\vec{x}, t)$

• 度规涨落（扰动引力场）： $\delta g_{\mu\nu}$

由于暴涨将微观尺度扩展到宇宙尺度，这些量子扰动被“冻结”并拉大，成为后来的大尺度结构和CMB起源。

---

✳️ 第二部分：Mukhanov–Sasaki 方程

---

✅ 1. 为何需要 Mukhanov–Sasaki 方程？

我们不能直接用 $\delta \phi$ 来描述扰动，因为它依赖坐标选择。 所以我们使用**规范不变变量（gauge-invariant variable）**来描述物理扰动。

---

✅ 2. 标量扰动的规范不变变量

引入：

• 曲率扰动 $\mathcal{R}$ ：在共动坐标下的空间曲率扰动；

• 它是与 $\delta \phi$ 、时空扰动 $\psi$ 组合后得出的物理变量。

然后定义 Mukhanov–Sasaki 变量：

$$v = z \mathcal{R}, \quad z = a \frac{\dot\phi}{H}$$

---

✅ 3. Mukhanov–Sasaki 方程推导（标量）

从爱因斯坦方程和 Klein–Gordon 方程线性化并组合可得：

$$v_k'' + \left( k^2 - \frac{z''}{z} \right)v_k = 0$$

• $k$ ：波数；

• $\tau$ ：共形时间（ $d\tau = dt/a$ ）；

• $v_k$ ：傅里叶空间中的模式函数。

这本质上是一个变系数谐振子方程。

---

✅ 4. 张量扰动（引力波）

张量扰动只涉及度规的横波部分 $h_{ij}$ （TT gauge）：

定义 $u = a h$ ，得到类似方程：

$$u_k'' + \left( k^2 - \frac{a''}{a} \right)u_k = 0$$

• 解 $u_k$ ，再除以 $a$ 得到引力波模式函数 $h_k$ 。

---

✅ 5. 方程的物理行为分析

子视界尺度 $k \gg aH$ ：

$$v_k'' + k^2 v_k \approx 0 \Rightarrow v_k \sim e^{\pm i k \tau}$$

• 表明扰动在深处视界内是自由平面波。

视界穿越 $k \sim aH$ ：

• 有效势 $z''/z$ 开始起作用，扰动幅度开始冻结。

超视界尺度 $k \ll aH$ ：

$$v_k \sim A_k z + B_k z \int \frac{d\tau}{z^2} \Rightarrow \mathcal{R}_k = \frac{v_k}{z} \to \text{const}$$

• 曲率扰动冻结，不再演化，为后续宇宙留下“印记”。

---

✅ 小结：为什么它重要？

• Mukhanov–Sasaki 方程连接了：

  • 暴涨时期的量子涨落；

  • 后期宇宙中结构形成的初始条件；

• 解这个方程后，能得到功率谱、谱指数、非高斯性等可观测量；

• 也适用于数值暴涨模型（非解析解情形）。

---

🧩 第一、二部分总结逻辑链：

【暴涨场 φ】
    ↓ (真空涨落 δφ + 时空扰动 δg)
【构造规范不变变量：𝓡】
    ↓
【Mukhanov–Sasaki 方程】
    ↓（k 模演化）
【冻结的超视界扰动】
    ↓（计算功率谱）
【P_𝓡(k), n_s, r】
    ↓
【CMB 温度/极化功率谱】
    ↓
【Planck/BICEP 实验约束】

---

接下来如果你需要，我可以继续第三部分（功率谱定义和慢滚推导）展开讲讲，或者直接给你 Python 解 MS 方程的模板。需要我怎么继续？

Response:

很好，我们现在进入整体理论链条的第三部分：从 Mukhanov–Sasaki 方程的解出发，如何精确计算功率谱、光谱指数 $n_s$ 、张量比 $r$ ，以及它们如何从理论连接到实际观测量。

---

📊 第三部分：功率谱的定义、计算与慢滚推导

---

✅ 1. 标量扰动功率谱的定义

定义功率谱 $\mathcal{P}_\mathcal{R}(k)$ 为：

$$\langle \mathcal{R}_{\vec{k}} \mathcal{R}_{\vec{k}'} \rangle = (2\pi)^3 \delta^{(3)}(\vec{k} + \vec{k}') \cdot \frac{2\pi^2}{k^3} \mathcal{P}_\mathcal{R}(k)$$

从 Mukhanov–Sasaki 变量 $v_k$ 得到：

$$\mathcal{P}_\mathcal{R}(k) = \frac{k^3}{2\pi^2} \left| \frac{v_k}{z} \right|^2$$

在视界穿越点 $k = aH$ 时评估。

---

✅ 2. 慢滚近似下的解析表达

在 quasi-de Sitter 情况下， $z''/z \approx \nu^2 - 1/4$ 型势可精确解得 $v_k$ ，解的长波极限（ $k \ll aH$ ）为：

$$|v_k| \sim \frac{1}{\sqrt{2k}} \cdot \frac{1}{k^{\nu - 1/2}} \propto \frac{1}{k^{3/2}} \Rightarrow \mathcal{P}_\mathcal{R}(k) \propto k^{n_s - 1}$$

具体慢滚表达为：

$$\mathcal{P}_\mathcal{R}(k) \approx \frac{1}{8\pi^2 M_{\rm Pl}^2} \cdot \frac{H^2}{\epsilon} \bigg|_{k=aH}$$

---

✅ 3. 光谱指数 $n_s$

谱指数定义为：

$$n_s - 1 = \frac{d \ln \mathcal{P}_\mathcal{R}(k)}{d \ln k}$$

用链式法则：

$$\frac{d}{d \ln k} = \frac{1}{H} \frac{d}{dt}$$

所以：

$$n_s - 1 = \frac{d \ln H^2}{H dt} - \frac{d \ln \epsilon}{H dt} = -2\epsilon - \eta$$

• 第一项来自哈勃变化；

• 第二项来自 inflaton 速度变化。

---

✅ 4. 张量扰动功率谱

张量扰动 $h_k$ 的定义与标量类似，使用重定义变量 $u_k = a h_k$ ，得到：

$$\mathcal{P}_h(k) = \frac{2}{\pi^2} \cdot \frac{H^2}{M_{\rm Pl}^2} \bigg|_{k=aH}$$

（因两个偏振态引入前因子 2）

---

✅ 5. 张量-标量比 $r$

两者比值为：

$$r = \frac{\mathcal{P}_h}{\mathcal{P}_\mathcal{R}} = 16 \epsilon$$

可直接用于连接理论与观测。

---

✅ 6. 总结表达式（慢滚模型下）：

物理量	表达式	物理意义
$\mathcal{P}_\mathcal{R}$	$\frac{1}{8\pi^2 M_{\rm Pl}^2} \cdot \frac{H^2}{\epsilon}$	标量扰动功率谱
$\mathcal{P}_h$	$\frac{2 H^2}{\pi^2 M_{\rm Pl}^2}$	张量扰动功率谱
$n_s - 1$	$-2\epsilon - \eta$	光谱倾斜度
$r$	$16 \epsilon$	张量-标量比

---

📡 第四部分预告：与观测的直接联系

现在我们有：

• 理论预测： $n_s, r$ ，取决于暴涨势 $V(\phi)$

• CMB 观测：Planck, BICEP, ACT 提供 $n_s^{\text{obs}}, r^{\text{obs}}$

我们可以画出模型预测的轨迹（如 Starobinsky, chaotic inflation）在 $(n_s, r)$ 平面上，然后用观测的椭圆区域做拟合比较。

---

---