Appearance
引力波天文学
概述
引力波天文学是利用引力波观测研究宇宙的新兴学科。通过探测时空涟漪，我们开启了观测宇宙的全新窗口，不再仅仅依赖电磁波，而是能够直接"聆听"宇宙中最剧烈的事件。
一、引力波探测原理
1.1 激光干涉仪
迈克尔逊干涉仪原理
现代引力波探测器的核心是激光干涉仪，基于迈克尔逊干涉仪原理。一束激光被分束器分成两束，分别沿着两条垂直的臂传播，在末端被反射镜反射回来重新汇合，产生干涉图样。
当引力波经过时，会拉伸和压缩空间，导致两条臂的长度发生微小变化。设臂长为 $L$，引力波引起的应变 $h$ 定义为：
$$ h = \frac{\Delta L}{L} $$
对于典型的引力波信号，$h \sim 10^{-21}$，这意味着对于 $L = 4\thinspace\text{km}$ 的臂长，长度变化仅为：
$$ \Delta L = h \times L = 10^{-21} \times 4 \times 10^3\thinspace\text{m} = 4 \times 10^{-18}\thinspace\text{m} $$
这比质子半径（约 $10^{-15}\thinspace\text{m}$）还要小三个数量级！
两束光重新汇合时的相位差为：
$$ \Delta \phi = \frac{4\pi}{\lambda} \Delta L $$
其中 $\lambda$ 是激光波长。通过精确测量干涉图样的变化，可以探测到引力波信号。
1.2 噪声来源
探测挑战
引力波探测面临极其严峻的噪声挑战。探测器必须在 $10^{-18}\thinspace\text{m}$ 的尺度上测量长度变化，这要求对各种噪声源进行严格控制。
地震噪声
地面振动是最主要的低频噪声源。地震、交通、海浪甚至微风都会引起地面振动，其幅度远大于引力波信号。地震噪声在低频段（$< 10\thinspace\text{Hz}$）占主导地位。
热噪声
探测器的反射镜和悬挂系统由原子组成，这些原子的热运动会引起镜面位置的随机涨落。热噪声在中频段（$10-100\thinspace\text{Hz}$）最为显著。根据涨落-耗散定理，热噪声功率谱密度与温度 $T$和机械损耗$\phi$ 成正比：
$$ S_x(f) \propto \frac{k_B T \phi}{f} $$
散粒噪声
激光由离散的光子组成，光子到达探测器的统计涨落称为散粒噪声（shot noise）。这是一种量子噪声，在高频段（$> 100\thinspace\text{Hz}$）占主导。散粒噪声与激光功率 $P$ 的关系为：
$$ S_h \propto \frac{1}{\sqrt{P}} $$
因此增加激光功率可以降低散粒噪声。
辐射压噪声
然而，增加激光功率会带来另一个问题：辐射压噪声。光子撞击反射镜时会传递动量，大量光子的统计涨落会引起镜面位置的随机运动。辐射压噪声在低频段显著，与散粒噪声形成权衡。
标准量子极限
散粒噪声和辐射压噪声之间存在基本权衡，两者共同决定了探测器的标准量子极限（Standard Quantum Limit, SQL）。这是量子力学对连续位置测量精度的基本限制。
1.3 噪声抑制
隔震系统
为了抑制地震噪声，现代引力波探测器采用多级隔震系统。LIGO 使用主动隔震和被动隔震相结合的方案：
- 主动隔震：使用传感器和执行器实时抵消地面振动
- 被动隔震：使用多级摆锤悬挂系统，利用摆的高频滤波特性
LIGO 的反射镜悬挂系统由四级摆组成，从粗到细逐级隔离振动，最终将高频段的振动抑制到 $10^{-18}\thinspace\text{m}$ 量级。
低温技术
为了降低热噪声，可以将探测器冷却到极低温度。KAGRA 探测器将反射镜冷却到约 $20\thinspace\text{K}$，显著降低了热噪声。低温技术是第三代探测器（如 Einstein Telescope）的关键技术之一。
压缩态光
量子非破坏测量
压缩态光（squeezed light）是一种量子态，其某个正交分量的噪声低于真空涨落，而另一个分量的噪声相应增加。通过注入适当相位的压缩真空态，可以在不增加辐射压噪声的情况下降低散粒噪声，或反之。
LIGO 和 Virgo 已经成功应用了频率相关的压缩态技术，在宽频带内降低了量子噪声。这是量子光学在引力波探测中的重要应用。
信号循环
信号循环（signal recycling）技术在干涉仪输出端放置一个额外的反射镜，形成信号循环腔。通过调节腔的参数，可以优化探测器在特定频段的灵敏度。这种方法可以"调谐"探测器，使其对特定类型的引力波源更敏感。
二、地面引力波探测器
2.1 LIGO
激光干涉引力波天文台
LIGO（Laser Interferometer Gravitational-Wave Observatory）是目前最成功的引力波探测器，由美国国家科学基金会资助，加州理工学院和麻省理工学院联合运行。
LIGO 有两个探测器：
- Hanford 探测器：位于华盛顿州 Hanford
- Livingston 探测器：位于路易斯安那州 Livingston
两个探测器相距约 $3000\thinspace\text{km}$，臂长均为 $4\thinspace\text{km}$。两个探测器的联合观测可以：
- 确认信号：排除局部噪声引起的假信号
- 定位源：通过到达时间差确定引力波源的天空位置
- 测量偏振：验证广义相对论对引力波偏振的预言
历史性发现：GW150914
2015 年 9 月 14 日，LIGO 首次直接探测到引力波信号 引力波与黑洞，命名为 GW150914。这个信号来自两个黑洞的并合：
- 黑洞质量：$36\thinspace{}M_\odot$和$29\thinspace{}M_\odot$
- 并合后黑洞质量：$62\thinspace{}M_\odot$
- 释放能量：$3\thinspace{}M_\odot c^2 \approx 5.4 \times 10^{47}\thinspace\text{J}$
这是人类历史上首次直接探测到引力波，验证了爱因斯坦 1916 年的预言，开启了引力波天文学的新纪元。2017 年，LIGO 的三位创始人（Rainer Weiss、Kip Thorne、Barry Barish）因此获得诺贝尔物理学奖。
2.2 Virgo
Virgo 探测器位于意大利比萨附近，臂长 $3\thinspace\text{km}$，由欧洲引力波组织（EGO）运行。Virgo 与 LIGO 联合观测可以：
- 改善定位：三个探测器的基线可以更好地确定源的位置
- 测量偏振：Virgo 的取向与 LIGO 不同，有助于测量引力波的偏振态
- 减少假信号：多个独立探测器提高信号可信度
在 O3 观测运行中，Virgo 的加入显著提高了引力波事件的定位精度，从数百平方度缩小到数十平方度，为多信使天文学观测提供了更好的目标。
2.3 KAGRA
KAGRA（Kamioka Gravitational-Wave Detector）位于日本飞驒市神冈矿山地下 $200\thinspace\text{m}$ 处，由日本国立天文台运行。KAGRA 的特色包括：
- 地下探测器：地下环境显著降低了地震噪声和人类活动噪声
- 低温技术：反射镜冷却到约 $20\thinspace\text{K}$，降低热噪声
- 量子技术：计划采用压缩态光和量子非破坏测量
KAGRA 于 2020 年加入观测网络，虽然目前灵敏度低于 LIGO 和 Virgo，但其技术创新为未来探测器提供了宝贵经验。
2.4 观测运行
LIGO 和 Virgo 进行了多次观测运行 相对论观测：
O1 运行（2015 年 9 月 - 2016 年 1 月）
- 首次探测到引力波 GW150914
- 探测到第二个黑洞并合事件 GW151226
- 确认了引力波探测的可行性
O2 运行（2016 年 11 月 - 2017 年 8 月）
- 探测到多个黑洞并合事件
- 首次探测到双中子星并合 GW170817
- 开启多信使天文学时代
O3 运行（2019 年 4 月 - 2020 年 3 月）
- Virgo 加入，三探测器联合观测
- 探测到约 90 个引力波事件
- 包括黑洞并合、中子星并合等多种类型
截至 O3 运行结束，LIGO-Virgo 网络共探测到约 100 个引力波事件，绝大多数是双黑洞并合，少数是双中子星并合和中子星-黑洞并合。
三、空间引力波探测器
3.1 LISA
激光干涉空间天线
LISA（Laser Interferometer Space Antenna）是欧空局（ESA）计划的空间引力波探测器，预计 2030 年代发射。LISA 由三个航天器组成，形成边长为 250 万公里的等边三角形，绕太阳运行。
LISA 的优势：
- 臂长巨大：$2.5 \times 10^6\thinspace\text{km}$，远大于地面探测器
- 低频探测：敏感频段为 $0.1\thinspace\text{mHz} - 1\thinspace\text{Hz}$，覆盖地面探测器无法触及的低频引力波
- 无地震噪声：空间环境没有地震噪声
LISA 的科学目标包括：
- 探测超大质量黑洞并合（$10^4 - 10^7\thinspace{}M_\odot$）
- 观测极端质量比旋进（EMRI）
- 探测银河系内的双白矮星
- 寻找原初引力波
LISA 面临的技术挑战包括：
- 精确测量 $250$ 万公里的臂长变化
- 消除太阳辐射压和宇宙射线的影响
- 维持航天器的无拖曳运动
3.2 太极计划
太极计划是中国科学院主导的空间引力波探测计划，设计与 LISA 类似，但臂长略长（约 300 万公里）。太极计划分为三个阶段：
- 太极一号（2019 年发射）：技术验证卫星，测试高精度惯性传感器和激光干涉测量技术
- 太极二号：双星系统，进一步验证关键技术
- 太极三号：完整的三臂引力波探测器
太极一号已于 2019 年成功发射，验证了多项关键技术，为中国空间引力波探测奠定了基础。
3.3 天琴计划
天琴计划由中山大学主导，采用与 LISA 不同的轨道方案。天琴的三个航天器位于地球轨道上，以地球为中心，臂长约为 17 万公里。
天琴计划的优势：
- 臂长较短：技术难度相对较低
- 地球轨道：通信和供电更方便
- 分阶段实施：可以先发射部分航天器进行科学观测
天琴计划的目标频段为 $0.1\thinspace\text{mHz} - 1\thinspace\text{Hz}$，与 LISA 互补，可以观测中等质量黑洞并合、双白矮星等引力波源 天体物理。
四、脉冲星计时阵列
4.1 基本原理
宇宙级引力波探测器
脉冲星计时阵列（Pulsar Timing Array, PTA）利用毫秒脉冲星作为精确时钟，探测纳赫兹频段的引力波。这是一种天然的、银河系尺度的引力波探测器。
毫秒脉冲星是快速旋转的中子星，自转周期极其稳定（精度可达 $10^{-15}$），可以与原子钟媲美。当引力波经过地球和脉冲星之间时，会拉伸和压缩时空，导致脉冲到达时间（Time of Arrival, TOA）发生微小变化。
对于一颗脉冲星，引力波引起的到达时间残差为：
$$ R(t) = \int_0^t \frac{\Delta \nu}{\nu}(t')\thinspace{}dt' $$
其中 $\Delta \nu / \nu$ 是脉冲频率的相对变化。
通过监测多颗脉冲星的到达时间残差，可以寻找引力波的特征信号——Hellings-Downs 曲线，这是一种特定的角相关函数，是引力波背景的独特指纹。
4.2 主要项目
NANOGrav
NANOGrav（North American Nanohertz Gravitational-Wave Observatory）利用阿雷西博望远镜和绿岸望远镜监测北美天空的毫秒脉冲星。
EPTA
EPTA（European Pulsar Timing Array）利用欧洲多个大型射电望远镜（包括 Effelsberg、Nançay、Lovell 等）监测欧洲天空的脉冲星。
PPTA
PPTA（Parkes Pulsar Timing Array）利用澳大利亚 Parkes 64 米射电望远镜监测南天脉冲星。
IPTA
IPTA（International Pulsar Timing Array）是 NANOGrav、EPTA 和 PPTA 的国际联合项目，整合全球数据，提高灵敏度。此外，中国也开展了CPTA（Chinese Pulsar Timing Array）项目，利用 FAST（500 米口径球面射电望远镜）监测脉冲星。
4.3 纳赫兹引力波
重大突破：2023 年
2023 年 6 月，NANOGrav、EPTA、PPTA 和 CPTA 联合宣布发现了纳赫兹引力波背景的强证据。这是脉冲星计时阵列多年观测的里程碑式成果。
观测到的信号特征：
- 频段：$1 - 100\thinspace\text{nHz}$（周期约 1 - 30 年）
- 应变幅度：$h_c \sim 10^{-15}$
- 角相关：符合 Hellings-Downs 曲线的预期
最可能的解释是超大质量双黑洞背景：
- 宇宙中大量超大质量双黑洞（$10^8 - 10^{10}\thinspace{}M_\odot$）在旋进过程中发射纳赫兹引力波
- 这些引力波形成随机背景，被脉冲星计时阵列探测到
其他可能的解释包括：
- 原初引力波（宇宙学起源）
- 宇宙弦
- 其他新物理
这一发现开启了纳赫兹引力波天文学的新窗口，为研究超大质量黑洞的形成和演化提供了全新手段。
五、引力波源
5.1 双黑洞并合
最常见的引力波源
双黑洞并合是目前探测到的最常见的引力波源。恒星级双黑洞由大质量恒星的演化形成，质量范围通常为 $5 - 100\thinspace{}M_\odot$。
双黑洞并合过程分为三个阶段：
旋进阶段（Inspiral）
- 双黑洞在引力辐射反作用下逐渐靠近，轨道频率逐渐增加
- 引力波频率 $f_{\text{GW}} = 2 f_{\text{orb}}$，其中 $f_{\text{orb}}$ 是轨道频率
- 波形可以用后牛顿近似精确描述
并合阶段（Merger）
- 双黑洞最终并合形成一个黑洞
- 强场动力学，必须用数值相对论模拟
- 引力波振幅达到最大
铃宕阶段（Ringdown）
- 新形成的黑洞通过发射引力波"铃宕"到稳态
- 波形由黑洞的准正模（quasi-normal modes）决定
- 频率和衰减时间由黑洞的质量和自旋决定
双黑洞并合的波形特征被称为啁啾信号（chirp signal），频率和振幅随时间增加。啁啾质量（chirp mass）定义为：
$$ \mathcal{M} = \frac{(m_1 m_2)^{3/5}}{(m_1 + m_2)^{1/5}} $$
啁啾质量决定了引力波波形的演化速率，是从观测数据中提取的最精确参数。
5.2 双中子星并合
多信使天文学的里程碑：GW170817
GW170817 是 2017 年 8 月 17 日探测到的首个双中子星并合事件，也是首次同时观测到引力波和电磁对应体的事件，开启了多信使天文学的新时代。
GW170817 的关键信息：
- 距离：约 $40\thinspace\text{Mpc}$（1.3 亿光年）
- 中子星质量：$1.36 - 1.60\thinspace{}M_\odot$
- 并合产物：可能形成大质量中子星或黑洞
千新星（Kilonova）
双中子星并合后，抛射出大量中子富集物质，通过 r 过程核合成形成重元素。这些重元素的放射性衰变加热抛射物，产生短暂的光学/红外辐射，称为千新星（kilonova）。
GW170817 的电磁对应体：
- 短伽马暴（GRB 170817A）：并合后约 1.7 秒探测到，持续时间约 2 秒
- 千新星（AT 2017gfo）：在光学和红外波段观测到，持续数天
- 余辉：在 X 射线和射电波段观测到，持续数月
r 过程核合成
双中子星并合是宇宙中重元素（如金、铂、铀等）的主要产地。r 过程（快速中子捕获过程）在中子富集环境中进行，可以在短时间内合成大量重元素。GW170817 的观测证实了双中子星并合确实是 r 过程核合成的场所。
5.3 中子星-黑洞并合
中子星-黑洞并合是第三种被探测到的致密双星并合类型。这类事件的特点：
- 潮汐效应：黑洞的潮汐力可以在并合前撕裂中子星
- 物质抛射：潮汐撕裂导致物质抛射，可能产生千新星
- 依赖参数：是否发生潮汐撕裂取决于黑洞质量、自旋和中子星物态方程
LIGO-Virgo 在 O3 运行中探测到几个中子星-黑洞并合候选事件，如 GW200105 和 GW200115。
5.4 超新星爆发
非对称爆发
大质量恒星（$> 8\thinspace{}M_\odot$）在生命末期会发生核心坍缩，形成超新星爆发。如果爆发是非对称的，会发射引力波。
核心坍缩超新星的引力波来源：
- 反弹激波：核心坍缩到核密度时反弹，产生强激波
- 对流不稳定性：激波传播过程中的对流
- 中微子加热：中微子加热引起的不对称性
- 旋转不稳定性：快速旋转的核心可能发生非轴对称不稳定性
超新星引力波信号的特征：
- 频段：$100 - 1000\thinspace\text{Hz}$
- 持续时间：约 $0.1 - 1\thinspace\text{s}$
- 波形：复杂的爆发信号，依赖于爆发机制
目前尚未明确探测到超新星引力波，但未来的第三代探测器有望实现这一目标。
5.5 连续引力波
持续不断的信号
连续引力波是由快速旋转的中子星发射的持续引力波信号。如果中子星存在不对称性（如山脉、椭率），就会持续发射引力波。
连续引力波的特征：
- 频率：$f_{\text{GW}} = 2 f_{\text{spin}}$（对于椭率驱动）或 $f_{\text{GW}} = f_{\text{spin}}$（对于单极子振荡）
- 持续时间：理论上可以持续数年
- 振幅：极其微弱，$h_0 \sim 10^{-26} - 10^{-25}$
中子星的椭率 $\epsilon$ 定义为：
$$ \epsilon = \frac{I_{xx} - I_{yy}}{I_{zz}} $$
其中 $I_{xx}, I_{yy}, I_{zz}$是主转动惯量。中子星地壳的"山脉"高度上限约为几厘米，对应的椭率$\epsilon \sim 10^{-6}$。
连续引力波探测的挑战：
- 信号极其微弱，需要长时间积分
- 需要精确知道脉冲星自转参数
- 计算量巨大，需要使用分布式计算（如 Einstein@Home）
5.6 随机引力波背景
宇宙的引力波噪声
随机引力波背景是由大量独立的引力波源叠加形成的随机信号，类似于宇宙微波背景辐射。
随机引力波背景可以分为两类：
天体物理起源
- 双黑洞背景：宇宙中所有双黑洞并合的叠加
- 双中子星背景：双中子星并合的叠加
- 超新星背景：超新星爆发的叠加
宇宙学起源
- 原初引力波：暴胀期间产生的引力波
- 宇宙弦：宇宙弦振荡和断裂发射的引力波
- 相变引力波：早期宇宙相变产生的引力波
- 预加热：暴胀结束后的预加热过程
随机引力波背景的探测对于理解宇宙早期物理具有重要意义。脉冲星计时阵列已经探测到纳赫兹频段的随机背景，LISA 和地面探测器有望探测到更高频段的随机背景。
六、波形分析
6.1 后牛顿近似
弱场慢速展开
后牛顿近似（Post-Newtonian Approximation）是处理双星系统旋进阶段的解析方法。它将爱因斯坦场方程按 $v/c$和$GM/(rc^2)$ 展开，适用于弱场和慢速运动的情况。
后牛顿展开的参数为：
$$ x = \left(\frac{GM\omega}{c^3}\right)^{2/3} \sim \left(\frac{v}{c}\right)^2 $$
其中 $\omega$是轨道角频率。后牛顿展开到$n$PN 阶意味着保留到 $O(x^n)$ 的项。
目前已知的高阶后牛顿结果：
- 3.5PN：完整的旋进波形（包括辐射反作用）
- 4PN：部分结果
后牛顿近似的适用范围：
- 适用于旋进阶段的早期，$v/c \ll 1$
- 在并合附近失效，需要数值相对论
引力波频率的演化方程为：
$$ \frac{df}{dt} = \frac{96}{5}\pi^{8/3}\left(\frac{G\mathcal{M}}{c^3}\right)^{5/3}f^{11/3} $$
这表明频率随时间增加，形成啁啾信号。
6.2 数值相对论
强场动力学的唯一工具
数值相对论（Numerical Relativity）是直接在计算机上求解爱因斯坦场方程的方法，用于模拟双黑洞并合的强场动力学阶段。
数值相对论的挑战：
- 坐标选择：需要选择合适的坐标条件避免奇点
- 稳定性：长时间模拟需要数值稳定性
- 计算资源：需要超级计算机和大规模并行计算
2005 年，三个研究组（Pretorius; Campanelli et al.; Baker et al.）独立实现了稳定的双黑洞并合模拟，这是数值相对论的重大突破，为引力波波形模板的构建奠定了基础。
数值相对论模拟的输出：
- 波形：引力波应变 $h(t)$ 作为时间的函数
- 最终黑洞参数：质量、自旋、反冲速度
- 辐射能量：以引力波形式辐射的总能量
6.3 有效单体方法
连接旋进和并合
有效单体方法（Effective One-Body, EOB）是一种半解析方法，将双体问题映射为一个粒子在有效度规中运动的问题。EOB 方法可以连接旋进阶段和并合阶段，提供完整的波形。
EOB 方法的核心思想：
- 将双体动力学映射为一个质量为 $\mu = m_1 m_2 / (m_1 + m_2)$的粒子在有效度规$g_{\mu\nu}^{\text{eff}}$ 中运动
- 有效度规包含后牛顿信息和数值相对论的校准参数
- 通过调节参数，使 EOB 波形与数值相对论波形匹配
EOB 波形的优势：
- 计算速度快，适合参数估计
- 物理图像清晰
- 可以外推到任意质量比
6.4 参数估计
匹配滤波
最优探测方法
匹配滤波（Matched Filtering）是从噪声中提取已知波形信号的最优方法。其基本思想是将观测数据与理论模板进行互相关：
$$ \rho = \frac{\langle d | h \rangle}{\sqrt{\langle h | h \rangle}} $$
其中 $d$ 是观测数据，$h$ 是理论模板，内积定义为：
$$ \langle a | b \rangle = 4 \text{Re} \int_0^\infty \frac{\tilde{a}^*(f) \tilde{b}(f)}{S_n(f)}\thinspace{}df $$
$S_n(f)$是噪声功率谱密度。信噪比$\rho$ 越大，信号越显著。
贝叶斯分析
参数估计使用贝叶斯推断方法。给定观测数据 $d$，参数 $\theta$ 的后验概率为：
$$ p(\theta | d) = \frac{p(d | \theta) p(\theta)}{p(d)} $$
其中：
- $p(d | \theta)$ 是似然函数
- $p(\theta)$ 是先验概率
- $p(d)$ 是证据（evidence）
对于高斯噪声，似然函数为：
$$ p(d | \theta) \propto \exp\left(-\frac{1}{2}\langle d - h(\theta) | d - h(\theta) \rangle\right) $$
参数后验分布
通过马尔可夫链蒙特卡洛（MCMC）或嵌套采样（Nested Sampling）方法，可以从后验分布中抽样，得到参数的估计和不确定性。
典型的引力波事件参数包括：
- 质量参数：$m_1, m_2, \chi_1, \chi_2$（质量和自旋）
- 轨道参数：距离、天空位置、轨道倾角
- 时间参数：并合时间、参考相位
参数估计的结果通常以角分布（corner plot）的形式展示，显示各参数的边际分布和相关性。
七、多信使天文学
7.1 引力波 + 电磁波
多信使天文学的黄金时代
GW170817 是多信使天文学的里程碑事件，首次同时观测到引力波和电磁波信号，开启了多信使天文学的新时代。
短伽马暴
双中子星并合被认为是短伽马暴（持续时间 $< 2\thinspace\text{s}$）的中心引擎。并合形成黑洞或磁星，驱动相对论性喷流，产生伽马射线暴。GW170817/GRB 170817A 证实了这一图像。
千新星
如前所述，双中子星并合抛射的物质通过 r 过程核合成形成重元素，产生千新星辐射。GW170817 的千新星 AT 2017gfo 在光学和红外波段被广泛观测。
余辉
并合产生的喷流与周围介质相互作用，产生从 X 射线到射电波段的余辉辐射。GW170817 的余辉持续了数月，提供了关于喷流结构和周围环境的重要信息。
7.2 引力波 + 中微子
限制引力波速度
引力波和中微子的联合观测可以限制引力波的传播速度。根据广义相对论，引力波以光速传播。
GW170817 的观测：
- 引力波到达时间：$t_{\text{GW}}$
- 伽马射线到达时间：$t_{\gamma} = t_{\text{GW}} + 1.7\thinspace\text{s}$
考虑到源的距离 $D \approx 40\thinspace\text{Mpc}$，可以限制引力波速度与光速的差异：
$$ \frac{|v_{\text{GW}} - c|}{c} < 10^{-15} $$
这是对广义相对论的强有力验证，也排除了许多修改引力理论。
IceCube 中微子天文台也在寻找与引力波事件相关的中微子信号，但目前尚未发现明确的关联。
7.3 科学成果
哈勃常数测量
双中子星并合提供了独立测量哈勃常数 $H_0$的方法。引力波信号直接给出光度距离$D_L$，电磁对应体给出红移 $z$，从而可以确定 $H_0$：
$$ v = H_0 D_L $$
GW170817 给出的 $H_0 = 70^{+12}_{-8}\thinspace\text{km/s/Mpc}$，虽然精度不高，但提供了一种独立的方法，有助于解决当前 $H_0$ 的"哈勃张力"问题。
中子星物态方程
双中子星并合的引力波信号包含关于中子星内部结构的信息。潮汐形变参数 $\Lambda$ 与中子星物态方程密切相关：
$$ \Lambda = \frac{2}{3} k_2 \left(\frac{c^2 R}{GM}\right)^5 $$
其中 $k_2$ 是 Love 数，$R$是中子星半径。GW170817 的限制$\Lambda \lesssim 800$ 排除了某些"硬"物态方程。
广义相对论检验
引力波观测提供了在强场和动态时空检验广义相对论的独特机会：
- 引力波速度：与光速一致
- 引力波偏振：符合广义相对论预言的张量偏振
- 后牛顿系数：与广义相对论一致
- 黑洞无毛定理：铃宕频率符合克尔黑洞的预言
八、未来展望
8.1 第三代探测器
灵敏度提升 10 倍
第三代地面引力波探测器计划将灵敏度提高约 10 倍，探测体积增加 1000 倍，可以观测到宇宙边缘的引力波事件。
Einstein Telescope
Einstein Telescope（ET）是欧洲计划的第三代探测器，位于地下，采用三角形几何（三条 10 km 臂）。ET 的特点：
- 地下位置：显著降低地震噪声
- 低温技术：反射镜冷却到约 $10 - 20\thinspace\text{K}$
- XYLO 配置：不同温度探测器覆盖不同频段
- 灵敏度：$h \sim 10^{-25}$（在 $1 - 1000\thinspace\text{Hz}$）
Cosmic Explorer
Cosmic Explorer（CE）是美国计划的第三代探测器，臂长 $40\thinspace\text{km}$，是 LIGO 的 10 倍。CE 的特点：
- 超长臂：直接提高灵敏度 10 倍
- 先进技术：压缩态光、更高功率激光
- 灵敏度：与 ET 相当
ET 和 CE 计划协同运行，形成第三代引力波观测网络。
8.2 新物理
探索未知物理
引力波天文学为探索基础物理提供了全新途径，有望发现超越标准模型的新物理。
原初引力波
暴胀期间量子涨落产生的引力波称为原初引力波。原初引力波的背景谱依赖于暴胀模型：
$$ \Omega_{\text{GW}}(f) \propto f^{n_t - 1} $$
其中 $n_t$ 是张量谱指数。探测原初引力波可以直接验证暴胀理论，确定暴胀能标。
LISA 和空间探测器有望探测到 $f \sim 10^{-4} - 10^{-1}\thinspace\text{Hz}$ 的原初引力波背景。
宇宙弦
宇宙弦是早期宇宙相变产生的拓扑缺陷。宇宙弦振荡和断裂会发射引力波，形成特征性的随机背景。探测宇宙弦引力波可以直接验证大统一理论和超弦理论。
相变引力波
早期宇宙的强一阶相变会产生引力波。相变引力波的频率和幅度取决于相变参数：
- 相变温度 $T_*$
- 相变强度 $\alpha$
- 相变速率 $\beta/H$
LISA 敏感频段对应电弱能标的相变，探测到相变引力波将直接证明电弱相变是强一阶的，这对解释宇宙重子不对称性具有重要意义。