原初引力波与阿里探测计划.pdf

原初引力波原初引力波与与阿里阿里探测探测计划计划 张新民 苏萌 李虹 李明哲 蔡一夫 2016 年 3 月公布的我国“十三五”规划纲要中列入了“强化宇宙演化、物质结构、生命起源、脑与认知等基础前沿科学研究” ，从国家科技发展的战略高度肯定了宇宙演化等基础前沿科学的重要性不久前的 2016 年 2 月 11 日，LIGO 实验组和美国自然科学基金委员会联合宣布探测到来自于十三亿年前由两个黑洞并合产生的引力波,这是人类第一次直接探测到引力波，可谓一项里程碑式的发现 ▍一、引力波源及其探测方式一、引力波源及其探测方式 引力波源大体可分为两种，天体物理起源和宇宙学起源对应不同的波源，相应的探测方式也不一样，如图 1 所示 图图 1 1，引力波的源及相应的探测方式，引力波的源及相应的探测方式 （图片来自（图片来自 NASA Goddard Space Flight CenterNASA Goddard Space Flight Center）） ▍（一）天体物理起源的引力波及其探测 天体物理起源包括以下三类： （（1 1））中子星、恒星级黑洞等致密天体（几十个太阳质量左右）中子星、恒星级黑洞等致密天体（几十个太阳质量左右）组成的致密双星系统组成的致密双星系统的合并过程的合并过程。

这类引力波的频率处于 10 赫兹- 1000 赫兹量级的高频段，相应的探测手段是地面激光干涉仪与迈克耳孙干涉仪一样， 在引力波激光干涉仪中激光被分裂成两束后在两个垂直的臂中传播并发生干涉当引力波经过时，局部的时空发生变形，两个臂的相对长度会发生改变，相应地激光干涉条纹就会移动此类实验最具代表性的就是 LIGO，利用长达 4 千米的两个互相垂直的探测臂首次探测到了引力波信号 （（2 2））大质量黑洞大质量黑洞并合并合过程的后期、过程的后期、银河系内的白矮双星系统银河系内的白矮双星系统频率为10−5赫兹 -1 赫兹，这类引力波信号可通过空间卫星阵列构成的干涉仪来探测，如欧洲的 eLISA 计划 （（3 3）） 超大质量黑洞 （数百万到数亿太阳质量）超大质量黑洞 （数百万到数亿太阳质量） 并合并合 频率为10−9赫兹 -10−6赫兹， 探测手段是脉冲星计时， 即利用地面上的大型射电望远镜，监视校准后的若干毫秒脉冲星如果其附近有大质量黑洞并合时发出的引力波，这些毫秒脉冲星的脉冲频率会有变化国际上 20世纪 70 年代就开始这方面的研究，90 年代已获得诺贝尔奖，我国在这方面有计划运行的 FAST 实验。

综上所述，天体物理过程产生的是高频引力波，相应的探测装置覆盖的频率范围在10−9赫兹以上 ▍（二）宇宙学起源的原初引力波及其探测 除了天体物理起源以外， 在宇宙的早期剧烈的量子涨落会产生充满整个宇宙空间的引力波，称之为原初引力波自 1940 年代以来，经典热大爆炸宇宙学取得了巨大的成功，其预言的宇宙轻元素丰度、宇宙微波背景辐射等均被实验证实， 其中关于微波背景辐射的研究成果还两度获得诺贝尔物理学奖但是大爆炸宇宙学本身并不完善，存在平直性疑难、视界疑难以及大尺度结构起源等困惑为了解决这些困惑人们提出了暴胀、反弹宇宙、循环宇宙等多种关于早期宇宙的理论或模型 这些理论模型的共同点是在辐射为主时期以前宇宙经历了一段极度非平凡的时期以暴胀为例，在暴胀期间宇宙急剧膨胀，支持暴胀的物质场的真空量子涨落得到增强并迅速扩张到视界外形成大尺度上的原初密度扰动 暴胀结束后宇宙进入减速膨胀的辐射为主时期，并按照经典大爆炸宇宙学理论那样演化至今存在于视界外的原初密度扰动又重新进入视界，在物质的引力不稳定性作用下，生成今天的大尺度结构和微波背景辐射的各向异性 但是暴胀宇宙学的一个重要预言是除了标量型的原初密度扰动以外， 宇宙时空本身的真空量子涨落也会迅速扩张到视界外形成张量型的原初扰动， 即原初引力波。

暴胀结束后，一部分原初引力波又逐渐进入视界因为暴胀发生在整个可观测宇宙中，原初引力波在宇宙中处处存在，形成引力波背景并遗留至今原初引力波与暴胀期间的物质成分没有直接关系，是纯时空涨落，而且是一种量子效应因此,探测原初引力波一方面有利于研究暴胀这样一个极高能标的早期宇宙动力学过程， 另一方面它是引力（即时空）的量子产物，有助于推进人们对量子引力这些基本物理问题的理解 原则上原初引力波在各个频段都有分布，但是随着宇宙的膨胀，不同频段的原初引力波的演化行为并不一样 它决定于如下的引力波方程： 方程左边的ℎ𝑗𝑘表示引力波，𝑎是衡量宇宙大小的标度因子，τ是共形时间，𝐻是牛顿引力常数，撇号表示对时间的导数，后面将要提到的宇宙膨胀率（即 Hubble 参数）为𝐼 = 𝑎′/𝑎2右边的引力波的源𝛱𝑗𝑘(𝜏, x)是物质部分的各向异性应力张量因为波动方程是线性的，引力波在不同尺度上的演化是独立的，适合对它做傅里叶分析经过傅里叶变换以后，波动方程变为： 这里𝑙是波数， 与频率的对应关系是𝑓 = 𝑙/(2𝜌𝑎) 这里可以看出对于固定的波数，频率以1/𝑎的方式逐渐红移，到今天（标度因子归一）相应的观测频率为𝑓 = 𝑙/(2𝜌)。

在暴胀期间，宇宙中不存在各向异性的应力张量，因此上述波动方程的右边为 0，引力波完全由真空涨落生成由上述方程我们可以发现，居于视界外（𝑙 ≪ 𝑎𝐼）的引力波保持不变，而进入视界内（𝑙 > 𝑎𝐼）的引力波类似于阻尼振子，阻尼来自于宇宙膨胀本身， ℎ𝑙具有振荡行为， 但是振荡的幅度在衰减刻画原初引力波的重要物理量主要有两个，第一个是功率谱： 它衡量的是引力波在不同尺度上的幅度如前面所述，暴胀模型还预言了密度扰动的功率谱，目前已经有很精确的测量，因此人们经常用引力波的功率谱与密度扰动功率谱的比值，即张标比𝑟，来表示引力波的幅度大小另一个描述引力波的重要物理量是能量谱： ， 其中 是宇宙临界能量密度，𝜍𝑔𝑤表示引力波的能量密度能量谱与功率谱之间的关系为 暴胀结束后引力波的功率谱的演化可以用下式来描述： 其中Δℎ2(𝑙,𝜏𝑗)是暴胀刚结束时的原初功率谱，𝑇ℎ(𝑙,𝜏)是转移函数一般的慢滚暴胀模型预言了近标度不变的原初功率谱，也就是说Δℎ2(𝑙,𝜏𝑗)与尺度或波数𝑙几乎没有依赖关系暴胀结束后不同尺度的引力波演化行为包含在转移函数中 它主要依赖于三种物理因素：（1）进入视界后由于宇宙膨胀带来的引力波幅度衰减， 也就是前面所说的类似阻尼振荡行为； （2）重新进入视界时宇宙的状态方程 w； （3）由中微子等相对论粒子的自由流（free streaming）产生的各向异性应力张量，这进一步压低视界内小尺度上的引力波幅度。

这些因素叠加在一起我们会发现进入视界的原初引力波的功率谱以及能量谱会被极大地压低不同波长的引力波重新进入视界的时间不同，波长越短（频率越高）的引力波进入视界越早，压低效果就越明显最后得到的能量谱如图 2 所示（𝑟 = 0.1） 图图 2 2，原初引力波的能量谱，原初引力波的能量谱 根据一般的慢滚暴胀模型的预言，遗留到今天的原初引力波在10−15赫兹以上的频段中都会被压低得很厉害，能量谱Ω𝑔𝑤基本上都在10−17以下在频率小于10−15赫兹的长波段，能量谱才会显著地大起来前面所述的引力波探测方式包括地面干涉仪、空间干涉仪和脉冲星计时只能探测频率在10−9赫兹以上的引力波在实验精度方面，目前或不久的将来升级的 LIGO、 Virgo 等地面干涉仪只能探测到最小能量谱为10−9的引力波 将来的 eLISA 空间干涉仪能探测到最小能量谱为10−10的引力波 采用脉冲星计时方法的 FAST 实验的探测精度能到达Ω𝑔𝑤~10−11，SKA 计划的精度能再提高两个量级但是这些工作在高频段的实验的探测能力比起一般慢滚暴胀模型预言的引力波能量谱还差好几个量级， 因此目前或不久的将来这几类实验是探测不到原初引力波的。

当然有些非常规的暴胀模型预言的原初引力波功率谱为标度依赖的蓝谱，它的能谱在高频段比较大，只有在这种情况下将来的高频引力波实验才有希望探测到原初引力波 因此， 就目前来说探测原初引力波最好的方式是宇宙微波背景辐射（CMB）的偏振（或极化）实验与高频引力波实验不一样，CMB的偏振实验探测的引力波频率范围处于10−15赫兹 -10−17赫兹， 从图2 可以看出，在这个频率范围内原初引力波的能量谱与高频段相比要高出好几个量级衡量 CMB 的探测能力最好用功率谱或者张标比𝑟，下面我们就简单地介绍一下 CMB 实验探测原初引力波的方法 CMB 是大爆炸遗留的辐射，在大爆炸宇宙早期，物质处于等离子体状态， 称为重子-光子流体， 其中光子与带电粒子发生频繁的碰撞由于宇宙膨胀，温度在逐渐降低当温度降到 0.1 电子伏特时，原子核与电子复合成中性原子，这个时期被称为复合期之后宇宙变得透明，遗留的光子在宇宙中自由传播（除了第一代恒星形成后发生的重电离过程以外） ，这就是 CMB复合期非常短暂，从垂直于视线方向来看可以看成一个很薄的球面，称为最后散射面，意味着光子在复合期结束时与电子发生最后的散射后开始自由传播。

CMB 充满了整个宇宙，并携带了大量的关于宇宙早期的信息从各个方向发射而来的CMB 具有差不多的温度以及具有10−5量级的涨落除了温度以外，光子还有偏振在复合期以前，由于光子被不停地散射，表现为没有偏振的自然光复合期结束时，自然光经历最后一次散射就形成了偏振光，并遗留下来 从观测的角度来看，CMB 的偏振图像可分解为两种独立的模式，一种是 E 模式，一种是 B 模式某一点附近的不同偏振模式的图像如图 4 所示，中心周围的线段表示偏振取向明显可以看出 B 模式偏振是“有旋”的，因而它在空间反射作用下具有与 E 模式不同的性质在图 4 中沿着通过中心并垂直于纸面的轴作空间反射， 我们将会发现E 模式图像保持不变，而 B 模式的 B>0 和 B 100效应比较明显而原初引力波产生的 B 模式偏振谱的峰值在𝑚 = 100左右人们正是通过寻找大尺度上的 CMB的 B 模式偏振信号来探测原初引力波 图图 6 6，标量扰动（左）和引力波（右）产生的，标量扰动（左）和引力波（右）产生的 CMBCMB 角功率谱角功率谱 （图片来源于（图片来源于 Challinor & Peiris, arXiv:0903.5158Challinor & Peiris, arXiv:0903.5158）） 这里有必要澄清一个观点，人们常说激光干涉仪比如 LIGO 是用来直接探测引力波的，而把 CMB 探测原初引力波的方式归为间接。

其实只要仔细比较两者的探测原理就可发现 CMB 探测原初引力波的方式与激光干涉仪实验一样直接，都是通过看光信号来探测引力波，没有谁比谁更直接的区别引力波经过时会产生时空变形，都是作为背景存在对于激光干涉仪，这种变形使得两个垂直的臂发生伸长和缩短，从而使得两臂交汇处的激光干涉条纹发生改变人们通过观测干涉条纹的移动来记录引力波，看的其实还是光信号对 CMB 来说，宇宙中的原初引力波经过时使得散射前电子周围的时空变形， 散射出来后就产生了CMB光子的B模式偏振， 同样是通过光信号来记录引力波因此对原初引力波探测来说 CMB 是最直接的方式,不比激光干涉仪更为间接不同的地方在于激光干涉仪测的是高频引力波，目前针对的是天体物理波源；CMB 探测的原初引力波的频率非常低,波长是宇宙学尺度，来源于早期宇宙激光干涉仪利用的是光的干涉，实现干涉的主要实验仪器是人造的；CMB 利用的是光的散射，实现散射的“实验装置”是天然的 相对来说原初引力波有一个优点是作为信号源它非常稳定，因为原初引力波充满整个宇宙空间，在同样充满整个宇宙空间的 CMB 上留下的印记也长久存在，可重复验证而 LIGO 探测的黑洞并合这样的事件偶然性比较大， 而且同一事件无法由后续其它实验来验证。