计算宇宙学与前沿观测

大学前沿课程。数值宇宙学模拟、引力波天文学、多信使天文学、JWST 最新发现，以及暗物质暗能量的前沿研究。

第一章 数值宇宙学模拟

1.1 为什么需要计算机模拟

宇宙学研究的独特挑战：

• 无法做实验：不能创造或操控宇宙
• 时间尺度巨大：宇宙演化以十亿年为单位
• 唯一性：我们只有一个可观测宇宙
• 非线性结构形成：暗物质晕、星系形成涉及复杂物理

计算机模拟的作用：

• 从初始条件（CMB）出发，模拟宇宙演化
• 预测可观测量的统计分布
• 与真实观测数据对比，检验宇宙学模型

1.2 N 体模拟

核心思想： 将暗物质离散化为大量粒子，计算引力相互作用。

关键算法：

1.3 主流宇宙学模拟项目

1.4 实践：用 Python 做简易 N 体模拟

import numpy as np

class NBodySimulation:
    def __init__(self, n_particles, box_size):
        self.n = n_particles
        self.box = box_size
        # 随机初始位置
        self.pos = np.random.rand(n_particles, 3) * box_size
        self.vel = np.zeros((n_particles, 3))
        self.mass = np.ones(n_particles) * (box_size**3 / n_particles)

    def compute_forces(self, softening=0.01):
        """计算引力（直接求和 + 软化）"""
        forces = np.zeros_like(self.pos)
        for i in range(self.n):
            for j in range(i+1, self.n):
                diff = self.pos[j] - self.pos[i]
                # 周期性边界
                diff = diff - self.box * np.round(diff / self.box)
                dist = np.linalg.norm(diff)
                force_mag = self.mass[i] * self.mass[j] / (dist**2 + softening**2)
                force_dir = diff / dist
                forces[i] += force_mag * force_dir
                forces[j] -= force_mag * force_dir
        return forces

    def step(self, dt):
        """Leapfrog 积分"""
        forces = self.compute_forces()
        acc = forces / self.mass[:, np.newaxis]
        self.vel += acc * dt
        self.pos += self.vel * dt
        # 周期性边界
        self.pos = self.pos % self.box

# 运行模拟
sim = NBodySimulation(n_particles=100, box_size=1.0)
for step in range(1000):
    sim.step(dt=0.001)

第二章 引力波天文学

2.1 引力波的发现

2015 年 9 月 14 日，LIGO 首次直接探测到引力波（GW150914）：

• 两个黑洞合并（36 + 29 → 62 太阳质量）
• 3 个太阳质量转化为引力波能量
• 信号传播了 13 亿年

意义： 开启了引力波天文学新纪元，2017 年诺贝尔物理学奖。

2.2 引力波探测器

2.3 多信使天文学

GW170817（2017 年 8 月 17 日）： 首次引力波 + 电磁波联合观测

• 引力波信号：双中子星合并
• 1.7 秒后：费米卫星探测到伽马射线暴
• 数小时后：光学望远镜观测到千新星（kilonova）
• 证实：重元素（金、铂）在中子星合并中产生

多信使组合：

引力波（LIGO/Virgo）→ 定位天区
    ↓
伽马射线（Fermi/INTEGRAL）→ 确认方向
    ↓
光学/红外（VLT/Hubble/JWST）→ 详细观测
    ↓
射电（VLA/ALMA）→ 余辉观测
    ↓
中微子（IceCube）→ 高能过程

第三章 JWST 与新一代观测

3.1 詹姆斯·韦伯太空望远镜（JWST）

2021 年 12 月 25 日发射，2022 年 7 月发布首批图像。

关键参数：

• 主镜直径：6.5 米（18 块镀金铍镜面）
• 观测波段：0.6 - 28.5 微米（红外）
• 轨道：日地 L2 拉格朗日点（150 万公里）
• 工作温度：~40K（被动 + 主动制冷）

3.2 JWST 重大发现

• 最遥远星系：JADES-GS-z14-0（红移 z=14.32，大爆炸后 2.9 亿年）
• 系外行星大气：首次探测到系外行星大气中的 CO₂、SO₂
• 恒星形成：以前看不到的尘埃云内部结构
• 早期宇宙：发现大爆炸后早期就有大量成熟星系，挑战现有模型

3.3 其他新一代设备

第四章 暗物质与暗能量前沿

4.1 暗物质候选者

4.2 暗能量本质

观测事实： 宇宙加速膨胀（1998 年发现，2011 年诺贝尔奖）

理论解释：

4.3 哈勃常数危机

问题： 不同方法测得的哈勃常数 H₀ 不一致：

5σ 差异 → 可能暗示新物理！

可能解释：

• 早期暗能量
• 中微子新物理
• 系统误差（尚未发现）
• 需要新的宇宙学模型

第五章 研究前沿

5.1 AI + 宇宙学

• 模拟加速：用神经网络替代 N 体模拟（速度提升 1000 倍）
• 参数推断：从观测数据直接推断宇宙学参数
• 异常检测：AI 发现观测数据中的异常模式
• 引力波信号处理：深度学习从噪声中提取引力波信号

5.2 21 厘米宇宙学

• 中性氢的 21 厘米辐射是探测"宇宙黎明"的关键
• SKA 将绘制宇宙大尺度结构的三维地图
• 可能揭示第一代恒星和星系的形成

5.3 量子引力与宇宙起源

• 圈量子宇宙学：大爆炸可能是"大反弹"
• 弦论宇宙学：膜碰撞产生宇宙
• 全息原理：宇宙可能是低维边界的全息投影
• 永恒暴胀：我们的宇宙只是多重宇宙之一

推荐资源

• 书籍：《Cosmological Physics》（Peacock）、《Numerical Cosmology》
• 模拟数据：IllustrisTNG 公开数据、SkySim
• 课程：Coursera "Astrophysics"（ANU）、EdX "Cosmology"
• 工具：Astropy（Python 天文库）、Gadget-4（N 体模拟代码）
• 数据：NASA ADS（论文）、SIMBAD（天体数据库）、NED（NASA/IPAC）