李清心,李兆升
(湘潭大学 物理与光电工程学院,湘潭411105)
小质量X射线双星(low-mass X-ray binary,LMXB)由致密星(中子星或黑洞)和一颗小质量主序星(M 通过RXTE,Swift和NICER卫星的观测数据,发现LMXB源中约20%的I型X射线暴中观测到了X射线暴振荡现象,振荡的频率范围约为11∼620 Hz,振荡的均方根振幅在5%∼20%之间,在有些源中,均方根振幅可接近50%。在吸积毫秒脉冲星中,观测到的XBO频率接近星体的自转频率,差异在几赫兹之内。在I型X射线暴的上升和衰减尾部,XBO的频率有几赫兹的漂移,通常在尾部渐近中子星的旋转频率,且XBO更容易在I型X射线暴的尾部探测到。 Strohmayer等人[1]认为X射线暴振荡的信号来自于中子星的自转,即在I型X射线暴爆发时,可能在中子星表面形成热斑,导致不均匀的温度分布,由星体自转在X射线能段产生周期性信号,但I型X射线暴的爆发过程中热斑如何形成和演化等问题尚未解决。表面模式模型(如Rossby模式)认为I型X射线暴可激发中子星海洋的振荡,在星体表面不同高度形成温度梯度,通过星体自转调制,形成XBO。这些模型只能解释部分观测的XBO现象,有待进一步的完善。 本文分别从I型X射线暴,X射线暴振荡的观测以及振荡理论分别进行介绍。 I型X射线暴是发生在中子星LMXB中的不稳定核燃烧,主要有H燃烧、He燃烧、C燃烧,其中C燃烧在观测上为超暴,由Grindlay等人[4]在4U 1820—30中首次发现。典型的I型X射线暴的光变曲线呈快速上升和指数下降形状,上升时间约1∼10 s(见图1),持续时间10∼100 s,爆发间隔时间为几小时到几天,能谱一般由黑体谱拟合,特征温度为1∼3 keV,辐射的总能量1032∼1033J。目前在115个中子星LMXB源中观测到了I型X射线暴①https://p ersonal.sron.nl/~jeanz/bursterlist.html。 图1 I型X射线暴的光变曲线及其功率密度谱[1] I型X射线暴的核燃烧过程敏感地依赖于中子星的局部吸积率,定义局部爱丁顿吸积率为: 其中,mp为质子质量,c为光速,X为大气中的H丰度,σTH为汤姆孙散射截面,R为中子星半径。对于太阳丰度的大气(X=0.73)且半径为10 km的中子星,其局部爱丁顿吸积率为8.8×104g·cm−2·s−1。 中子星表面温度较低时,H燃烧由冷CNO循环主导,少数情况会发生热CNO循环。而中子星表面温度足够较高时,H燃烧主要是rp过程。冷CNO循环过程如下: 循环过程中,C,N,O元素捕获质子,并在反应的最后一步释放α粒子,C,N,O的生成速率与消耗速率相同,是循环过程的催化剂。冷CNO的反应速率依赖于捕获质子的速率,尤其是14N(p,γ)15O反应。当吸积率低于0.01m˙Edd,H经冷CNO循环过程燃烧,反应速率非常快,这个过程释放的能量与T16(T为环境温度)成正比,反应过程不稳定,同时会触发He燃烧,观测上表现为H/He混合燃烧I型X射线暴。 He丰度较高的壳层可以在稳定的H燃烧(吸积率0.01m˙Edd∼0.1m˙Edd)的壳层之下形成,或者通过吸积He白矮星、大气成分为He的主序星形成。在He壳层中,当吸积率为0.1m˙Edd∼1m˙Edd,环境温度较高(T>108K),简并的电子密度达到106g/cm3时,会触发不稳定的He燃烧(反应速率∝T30,此过程极其不稳定),此时He燃烧过程是α+α+α→12C+γ。这个过程分两个阶段,第一阶段是两个α粒子结合成8Be,8Be的寿命只有2.6×10−16s;此时会源源不断地形成8Be,从而有足够的8Be粒子与α结合成C的激发态12C⋆。12C⋆可以衰变成8Be并释放α粒子,或者辐射γ光子变成C基态,这是3α过程的第二个阶段。观测上表现为He燃烧I型X射线暴,与H/He混合燃烧I型X射线暴相比,其光变曲线的上升下降时间更短。 在很少的几个源中(例如4U 1820−30),观测到一类特殊的X射线暴,也称为超暴(superburst),与上述两类I型X射线暴相比,超暴释放的总能量高约103倍,持续时间长约103倍,发生率约为这类暴可能是中子星大气中更底层的C壳层的燃烧。目前,超暴的发生条件尚不清楚,可能需要额外的加热机制(例如,中子星壳层的超密态核反应),以克服C聚变的库仑势垒。 对XBO的探测,需要高能X射线望远镜具有高时间分辨率、大有效面积、光子堆积和死时间效应不显著等特性。自1996年RXTE首次发现中子星LMXB的XBO起,目前共在30个源中发现了这类现象,XBO频率见表1,其中确定暴振荡频率的源按脉冲星和非脉冲星分类,XBO频率有待进一步观测验证的源单独分出。除RXTE外,人们还从Swift/BAT和NICER卫星的观测中搜索到了XBO信号。下面,将分别介绍XBO信号的搜索方法、信号的概率估计和X射线望远镜对XBO的观测。 表1 探测到XBO的源[7,9–11]Hz (续表) 搜索周期性的XBO信号,可通过快速傅里叶变换(fast fourier transformation,FFT)或统计方法[5],它们分别针对并道数据及事例数据①事例数据记录了每个光子的到达时间,且可转换为并道数据;并道数据记录每一个时间间隔中的光子计数率。。 对于FFT,将总持续时间T的光变曲线均匀划分为N个并道,且N=2m(m是整数),时间序列记为x k(k=0,1,2,3,···,N−1),表示第k个并道中的光子数,Leahy归一化功率密度谱P j为[6]: 其中,νj是傅里叶频率(νj=j/T,j=0,1,2,3,...,N/2),Ntot为总光子数。功率密度谱的最小频率为νmin=1/T,最大频率(即奈斯奎特频率)为νmax=N/(2T),频率间隔为∆ν=1/T。 其中,n为谐波数,ϕj为光子相位: 其中,t0是参考时间,t j是光子相对于t0的到达时间,ν(t)是频率模型,当ν(t)不随时间变化时,ϕj=2πνt j。 除了振荡频率,功率谱还包含有关脉冲幅度的信息。对于Leahy归一化的功率谱,均方根分数(root mean square,rms)振幅定义为: 括号项用来校正背景因素影响,Nγ是总光子数,B是背景光子数①一般取暴前的光子平均计数率乘以计算功率密度谱所对应的时间段宽度,作为背景光子数,P s是信号功率,如果使用方法,则将公式(5)中的P s替换为 3.2.1 估计快速傅里叶变换的真实信号概率 由FFT方法计算I型X射线暴的功率密度谱,低频包含红噪声和泊松噪声,而高频部分主要由泊松噪声组成。泊松噪声产生的功率密度具有两个自由度的χ2分布[7]: 其中,R n和I n分别为噪声傅里叶变换系数的实部和虚部,根据Groth[8]解析表达式,测量功率P m和信号功率P s的概率分布: 其中,I是第一类的修正贝塞尔函数,而k是谐波数(一般取1)。对于仅有噪声的情况,P s=0,k=1,公式(7)为具有两个自由度的χ2分布的概率密度函数(probability densityfunction,pdf),等于 对于k>1,p k(P m,0)满足自由度为2k的χ2分布。 图2 时频窗口可视化[43] 假设噪声为泊松分布,由噪声产生噪声分布δ的测量信号概率Pnoise为: 根据Pnoise≤1%,Muno等人[44]所定义的三个标准来判断重要信号: (1)将N=N t×Nν=16×10=160次检验和Pnoise≤1%代入式(9),可以得到噪声分布δ≤7×10−5。 (2)计算暴上升阶段第1 s,即N t=1,Pnoise≤1%时由公式(9)得到噪声分布δ≤10−3。 (3)一个时频窗口周围有两个相邻的频率窗口以及一个时间窗口,计算所选时频窗口和其中一个相邻窗口的噪声概率:两个窗口单独包含噪声的概率由两个窗口独立的噪声概率的乘积给出,Prob1,2=P1(N1,δ1)×P2(N2,δ2),因可取3个相邻窗口计算噪声概率,所以 其中一个Nν=3,可近似表示为: 当Prob1,2≤10−2时,将其代入式(10),则有δ1δ2≤1.3×10−6。 其中,k是谐波数(一般k=1),而I k−1是第一类修正的Bessel函数。给定Z s,通过相关的累积分布函数获得介于0和Z m之间的测量功率的概率是: 在给定Z m的情况下,真实信号功率在0到Z s间的概率是: 如果在整个爆发中均未找到振荡信号,则选择最大的功率作为信号的上限。 3.2.3 蒙特卡洛模拟 蒙特卡洛模拟方法可用来评估无信号时噪声功率的分布,且能精确模拟实际观测的数据,如设定的时间窗口等。首先选择合适的模型拟合I型X射线暴的光变曲线,然后以此为基础生成大量具有泊松计数统计信息但没有任何周期性信号的光变曲线样本。从这些模拟出来的光曲线中获取功率谱,即可得出仅由噪声构成的功率密度谱。在测得的振荡信号附近统计功率分布,此分布即为噪声分布,可判断出测量到的XBO信号来自于噪声的概率。 Strohmayer等人[1]将整个I型X射线暴的光变曲线应用FFT方法,发现功率密度谱中出现了多峰结构(见图1),这是因为在整个爆发期间XBO的振荡频率随时间演化。为了有效地搜索和追踪XBO信号,一般用动态功率密度谱方法分析I型X射线暴。 图3给出了产生动态功率谱的示意图,图3a)为I型X射线暴的光变曲线,可将光变曲线划分为多段,这里以三段为例,∆t0,∆t1,∆t2是移动窗口的时间长度,相邻的窗口之间有重叠,其中T0,T1和T2为窗口的中间时刻。图3b)分别计算每段光变曲线的功率密度谱,蓝色代表低功率,紫色代表中功率,红色代表高功率。图3c)将图3b)的功率密度谱和相对应的中间时刻结合,给出时间-功率-频率三维图。最后在图3d)中将图3c)功率的等高线投影在光变曲线上,得到动态功率密度谱,可简洁地展示XBO出现的时间、XBO信号的频率和强度随时间的演化等信息。 图3 动态功率密度谱[9] RXTE卫星于1995年12月30日发射升空,在长达16年的运行后,于2012年1月5日正式退役,该卫星是到目前为止研究XBO最重要和最成功的高能时变探测器。RXTE卫星采用事例格式记录高时间分辨率数据,2019年,Bilous和Watts[7]用FFT方法系统搜索和分析了所有RXTE观测的I型X射线暴的XBO信号,搜索了所有57个中子星LMXB源的2118次爆发,列出了爆发峰值的信噪比,提取了2∼60 keV能量范围内的光变曲线,搜索时间窗口为2 s,移动步长为0.5 s,提供了暴振荡频率以及在爆发中出现的位置等信息,并计算了振荡的分数振幅,对已经确定振荡频率的源进行搜索,确定了已知的17个XBO频率与之前测的频率或者脉冲星自传频率相等;对之前一些论文中所搜索到的爆发振荡但还有待检验的源进行检验,其中一些源并没有搜索到暴振荡信号,另一些源搜索到的振荡信号过于微弱,不排除是噪声信号(关于XBO频率有待证实的源的具体信息,可见表1);还发现了SAX J1810.8−2609在531 Hz的爆发振荡信号以及IGR J17473−2721频率在600 Hz左右的振荡信号,该源不同暴的振荡频率之间相差在3 Hz以内,很可能是一个新的XBO源,其准确振荡频率还需要进一步检验。 图4是RXTE观测到的典型的XBO动态功率密度谱,其中图4a)SAX J1808.4−3658和图4b)XTE J1814−338是两个持续吸积脉冲星,图4c)HETE J1900.1−2455以及图4d)Aql X−1为间歇吸积脉冲星,而图4e)4U 1636−536、图4f)4U 1728−34都是非脉冲星。 图4 不同类型源的动态功率密度谱[9] 研究发现SAX J1808.4−3658在暴开始时就出现振荡信号,并且随着暴进行向上或者向下有几Hz的漂移,在暴峰值振荡突然消失,之后振荡会在稍高或稍低频率处重新出现并保持稳定;XTE J1814−338的振荡频率变化在一个傅里叶频率分辨率之内,振荡信号也不会在暴峰值期间消失,而且经常出现在暴的尾部;持续吸积脉冲星振荡频率与自转频率非常相近,在几Hz之内。HETE J1900.1−2455只有一个时间窗口检测到爆发振荡信号的峰值,其余频率比较低的振荡信号会持续一段时间;Aql X−1和HETE J1900.1−2455等间歇吸积脉冲星的爆发振荡频率,在整个爆发过程中有缓慢的漂移,渐近最大频率且比自转频率低约1 Hz。4U 1636−536爆发振荡信号存在于暴上升阶段;4U 1728−34的振荡信号存在于整个爆发过程中;4U 1636−536和4U 1728−34等非脉冲星所观察到的爆发振荡,随着爆发的进行,与间歇吸积脉冲星相似,频率约1 Hz向上漂移,在暴的尾部渐近最大频率。 RXTE还从中子星LMXB的超暴中发现了XBO现象。Strohmayer和Markwardt[23]在2001年观测到的4U 1636-53的超暴光变曲线和动态功率谱见图5,在动态功率谱中还可以看到明显的频率漂移,来自于中子星在其双星轨道上的多普勒运动。超暴振荡频率比一般I型X射线暴的振荡频率高出约0.4 Hz,比正常爆发振荡更接近中子星的真实自转频率。另外,在LMXB 4U 0614+09[38]中也发现了超暴的爆发振荡信号。 图5 4U 1636-53超暴的光变曲线和动态功率密度谱[23] 2017年6月3 日,NASA的中子星内部成分探索器(NICER)发射成功,它的目标是通过精确测量中子星的质量和半径来研究中子星的物质状态和内部结构。NICER提供了56个准直X射线聚光器阵列,每个都配有一个硅漂移探测器,这些探测器工作在0.2∼12 keV能段,在1 keV处的能量分辨率约为85 eV,在1.5 keV处的有效面积超过2000 cm2,时间分辨率小于300 ns,显著地优于RXTE。NICER成为继RXTE之后研究XBO最主要和最优秀的设备。 3.5.1 SAX J1808.4−3658 Bult等人在2019年[12]用NICER对SAX J1808.4−3658源进行观测,发现了预期的401 Hz频率处的XBO,且没有明显的漂移。他们还给出了XBO的脉冲轮廓和持续辐射的脉冲轮廓,两种波形如图6所示,可以看出XBO的轮廓和幅度与持续吸积脉冲相似,这表明两者都源于中子星表面的一片受限制区域的燃烧(热斑,详见第5.1节),但是XBO领先脉冲34◦±7◦,XBO的波形与能量的相关性与持续吸积脉冲星不同,一般在低能谱测得暴发振荡信号。 图6 暴振荡波形与脉冲波形相比较图[12] 3.5.2 4U 1728−34 从开始工作到2019年2月间,NICER探测到了4U 1728−34的7次I型X射线暴,Mahmoodifar等人[27]发现其中的3次爆发有XBO现象,且都出现在其暴衰减尾部(如图7所示),频率非常接近,与之前的观测结果一致。 图7 NICER观测到4U 1728−34的光变曲线上的动态功率密度谱[27] 比较这三个在暴尾部的振荡,表明暴在高能带中可以更强,并且可以具有不同的幅度和特性。其中的两个振荡只在高能段(E>6 keV)中被检测到,并且它们具有非常大的分数均方根振幅(48±9)%和(46±9)%。热斑模型(第5.1节)和表面模式振荡模型(第5.2节)都难以解释如何在爆发衰减尾部出现的大振幅的振荡。如果中子星表面的所有位置都以相同的方式加热和冷却,则无法在中子星表面产生足够高的温度梯度;而冷却尾流(cooling wake)模型[45](详见5.3节)表明,如需要在中子星表面形成很高的温度梯度,则要中子星表面先燃烧部分比后燃烧部分更早冷却,在燃烧到赤道点时火焰会迅速蔓延至整个赤道,进而使整个中子星表面的温度梯度增加,所以一般温度梯度存在于赤道以及附近的区域[27]。 吸积毫秒脉冲星(accreting millisecond X-ray pulsars,AMXPs)是吸积供能(通过吸积伴星的物质,将引力势能转换为X射线辐射)的快速转动中子星,正是在AXMPs中观测到了XBO,才证实了XBO信号来自于中子星自转。研究AXMPs的XBO可以更好地了解爆发振荡的物理机制。这类源的自转频率可以被独立测量,而且,它们是可以用来研究星体磁场影响以及燃料分布不均匀作用的唯一一类源。 AMXPs的XBO频率与其自转频率非常接近,通常都是在几Hz以内,但是不同AMXP的XBO频率与自转频率的间隔存在差异。例如XTE J1814−338的XBO频率非常接近自转频率;而SAX J1808.4−3658和IGR J17511−3057的暴振荡在暴上升阶段快速增加,之后略微超过自转频率,最后在暴尾部会稳定在自转频率为10 mHz左右[2,15]。AXMPs与非脉冲星的XBO现象也有差异。对于非脉冲星的I型X射线暴,当源处于软态(高吸积率)时,才会出现XBO现象[19]。相比之下,即使在硬态(低吸积率),AXMPs中也能检测到XBO现象。AXMPs和非脉冲星的暴振荡持续时间也存在着不同,对于SAX J1808.4−3658和XTE 1814−338,暴振荡现象会在整个暴中持续;对于非脉冲星,虽然有时会在整个爆发中检测到爆发振荡,但是大多暴振荡都出现在暴的尾部[19],这点和间歇吸积供能脉冲星一样。IGR J17511−3057是上面两种现象的过渡:在较低吸积率下,整个暴中都可以检测到振荡,但是随着吸积率的增加,爆发的振荡信号会从暴的上升阶段中消失[15]。AXMPs的XBO信号比间歇脉冲星和非脉冲星的谐波含量更高[2,15,46],并且它们振幅与能量的关系也不一样,AXMP的XBO振幅随能量的上升而下降,而间歇脉冲星和非脉冲星的XBO振幅随能量的上升而上升[47,48]。XBO振幅在不同的源之间也有很大的区别,大多数情况下都在爆发时吸积供能脉冲幅度的1%到10%[2,15];但是对于过渡性AMXP射电脉冲星IGR J18245−2452是一个例外:对于检测到的其中一个暴,爆发的振荡幅度远高于爆发由吸积动力产生的脉冲幅度[49]。目前,对造成这些脉冲星和非脉冲星暴振荡之间特性差异的原因尚不清楚,动态磁场可能起主要作用。 在XBO的观测中发现,振荡信号在大多数暴中频率会向上漂移(频率增加),一般出现在暴的尾部,在一些暴中振荡信号有着很高的振幅,吸积脉冲星的振荡频率与脉冲星自转频率的偏差在几Hz之内,研究爆发过程中的点火和火焰蔓延方式对了解暴振荡的物理机制起着重要的作用。目前提出了热斑模型、表面模式振荡模型(如图8)和冷却尾流(cooling wake)模型来解释这些XBO现象,Cavecchi和Spitkovsky[50]在2020年模拟了爆发过程中不稳定火焰前段的3D模型,不稳定的火焰前端导致整体燃烧速度加快并出现明显火焰涡流,这些涡流也与振荡频率的漂移有关,它们很可能是激发表面振荡模式的原因。 图8 爆发振荡简化示意图[9] 现在一些暴的振荡信号还很难观测到,未来由中国主导的eXTP卫星(the enhanced X-ray timing and polarimetry mission),有更强的X射线时变观测能力,将为研究暴振荡机制提供更多的条件与机会。 为解释暴振荡频率向上漂移和振荡主要出现在暴尾部的现象,人们提出了两种热斑模型:热斑传播即热斑在暴上升阶段中临时产生,火焰从点火点扩散到整个中子星;另一种是将燃烧限制在中子星表面的一小片区域内,热斑持续存在。 Strohmayer等人[51]首先提出,火焰在中子星上扩散时的水平运动可能会导致振荡频率的漂移,为进一步解释频率总是逐渐上升的现象,提出另一种假设:在加热燃烧层的垂直膨胀过程中保持角动量守恒,快速膨胀会使热斑旋转得更慢,从而使观测到的频率低于自转速率,并且随着层的冷却和收缩频率会升高,逐渐达到自转频率。在大多数暴尾部,拟合暴能谱所获得的黑体半径与中子星的半径相似,这说明火焰已经在整个中子星表面燃烧[19]。但是积累爆发物质需要数小时甚至数天,而热核失控不到1 s[52],为了同时满足整个中子星表面的点火条件,热状态必须保持一致[53];而吸积中存在一些不对称性现象(例如由于磁通道或赤道边界层的存在),使得这种情况(尤其是在AXMP中)不太可能发生。所以,如果火焰在某一点开始燃烧,那么火焰需要在整个中子星表面扩散,但是目前还难以解释燃烧物质在这时可以达到所需要的热均衡状态。 散布热斑模型:随着整个暴能谱黑体半径增加,热斑尺度增大,振荡幅度下降[22,46,51,54];可以解释暴的上升阶段中存在爆发振荡,但是无法解释大多数暴振荡信号都出现在暴的尾部。其中一个解释是因为观测效应:暴上升阶段持续时间很短,很难检测[9]。XBO探测难易程度的另一个因素是点火处的纬度,在快速旋转的中子星上,赤道处的有效重力较低,吸积率最大,这样更容易积累点火所需要的物质[55]。研究发现,赤道处的I型X射线暴,会更容易探测到爆发振荡[56]。但是,在边界层附近,吸积物质的核反应过程介于稳定燃烧与不稳定燃烧之间,可能在高纬度处点火[57]。对磁场较高的源,由于局部密度过大或者过度加热,也可能在磁极处优先点火[58]。火焰如何扩散对暴振荡的可探测性同样重要,因为火焰必须能够在不对称方向上持续传播。多年以来,控制火焰扩散的过程在I型X射线暴中一直是一个悬而未决的问题,它受热传导、湍流和对流等各种传热机制的影响[9]。Spitkovsky等人[55]在2002年指出流体动力学效应的重要性,特别是科里奥利力与热燃烧物质膨胀之间的相互作用,科里奥利力的作用是减缓火焰传播,从而保留热斑,这种约束对于控制火焰的扩散非常重要,因为热物质的快速膨胀和横向扩散可能会导致火焰失速[59]。但是,科里奥利力约束对于IGR J17480−2446(自转频率仅为11 Hz)无效[33]。 同样暴尾部的高振幅振荡现象也没有准确的解释,第一种机制是热斑得以保留,燃烧火焰并未扩散到整个中子星上。让热斑保留下来的一种方法是通过中子星磁场限制燃料。吸积物质沿磁力线进入中子星磁极,压力足够大时场线变形,吸积物质才能越过磁力线。Brown和Bildsten[60]的结果表明,要使物质约束在正常I型X射线暴的点火深度内,其场强至少应为106T,这要远远高于LMXB的磁场强度[61,62],而磁流体动力学的不稳定性也可能使磁约束无效[63],因此,这种机制的可能性并不大。另一种机制是限制火焰前端:燃烧一旦开始,蔓延一段距离后停止传播,这需要以某种方式阻止其热传递。Cavecchi等人[33]认为燃烧物质与磁场的径向分量相互作用,充分放大磁场的水平分量后,也能阻止热传导。对于低磁场的LMXB,这些过程难以出现,除非在爆发过程中有某些机制(例如对流发电)暂时增加磁场[64]。 表面模式振荡模型认为,I型X射线暴点火以及火焰在中子星周围的扩散可能会激发中子星海洋(由离子和电子的等离子体组成,分布在中子星外层到等离子体可以结晶成固体的深度)中的大尺度波浪,海洋中不同高度处温度有差异,以致X射线的亮度不同,通过星体转动,表现为XBO。人们借鉴了地球海洋表面模式振荡的研究方法,写出了旋转球体振荡的流体动力学连续性方程(质量守恒),动量和能量守恒方程,并简化为拉普拉斯潮汐方程式[65–67]。这个系统的恢复力是浮力和科里奥利力。拉普拉斯潮汐方程的非轴对称模式解依赖于exp(imφ)的方位角φ,其中方位角特征值m是一个整数,给定一个中子星的旋转模式具有方位角数量m和频率νr的模式,观测到的频率为: 其中,νs是中子星自转频率,而νr的符号是正号还是负号,取决于模式是沿自转方向(向东)还是逆自转方向(向西)传播。如果νr取决于海洋的热状态,就可以解释XBO频率漂移现象,即振荡频率会随着爆发的进行和尾部的冷却而改变。 拉普拉斯潮汐方程式有三种模式。 (1)庞加莱模式。向逆方向(向西)传播,并在不旋转的临界处变为纯重力波动(仅靠浮力恢复),因此通常被称为g模式[68,69]。 (2)开尔文模式。沿与自转相同方向(向东)传播,在不旋转的临界处变为纯重力波动。但与庞加莱模式不同的是,开尔文模式处于地转平衡状态(水平气压梯度力与科里奥利力相平衡),因此只涉及到方位角运动,并只限于在赤道上进行火焰传播。 (3)Rossby模式。沿逆方向(向西)传播,也称为r模式[70],在零浮力的极限下,Rossby模式是由科里奥利力的纬度变化驱动的。在不旋转的临界处变为零频率的一般解。如有不同燃烧物质之间的浮力存在,则也将起到恢复模式的作用,因此也通常称为浮力r模式。 Heyl[69]认为浮力r模式是最有可能的,例如:逆向传播的模式中|νr|仅有几赫兹,给出的观测频率略低于自转频率;随着燃烧层冷却,|νr|变小,观测到的频率将趋向于自转频率。而庞加莱模式(或g模式)预计的|νr|要大很多,所以被排除[9,69]。开尔文模式被排除原因是频率将随着海洋冷却而下降,这与大多数观测到的频率漂移相反,开尔文模式和庞加莱模式也都将燃烧严格限制在快速旋转的中子星上的赤道区域,这将降低整体暴振荡幅度。浮力r模式预言燃烧层冷却时的频率漂移约10 Hz,比观测值高一个量级。Piro和Bildsten[71]提出振荡出现在燃烧层底部与中子星壳层之间的区域,振荡频率漂移约4 Hz,依然高于AMXP源的测量值。Heyl upciteHey04指出浮力r模式如果限制在光球层,则漂移频率可能更小,细节有待研究。一些暴的光度会暂时超过爱丁顿光度,热能转换为动能和引力势能,使中子星最外层发生膨胀,称为光球半径膨胀(photospheric radius expansion,PRE),通常会从PRE暴中检测到爆发振荡信号,但是在PRE阶段本身很少检测到爆发振荡。为了解释这种在峰值阶段没有检测出振荡信号的现象,Cooper[72]考虑燃烧层高效对流能量传输机制,分析了单区域模型,发现对流确实抑制了振荡的发生,可以解释一些暴振荡在暴尾部的现象。 Chambers和Watts[73]考虑浮力r模式中的参考系拖曳效应和引力红移效应,计算发现频率漂移约1∼3 Hz,与观测结果相符,但是给出的TBO频率与自转频率的偏差仍然比观测值大。 目前在一些X射线暴的尾部检测到振幅很大的爆发振荡信号,因为缺少明显的高梯度温度差,普通的热斑模型无法解释这种现象。Cumming和Bildsten[53]认为可能是由于中子星表面冷却而形成温度梯度导致了暴尾部的振荡现象,这即为冷却尾流模型。 Mahmoodifar和Strohmayer[74]提出了三种冷却尾流模型,其中“规范”冷却尾流模型分为两种:一种假设中子星上每一片区域以相同的方式燃烧和冷却,中子星表面一旦从某点开始点燃燃烧,火焰扩散到每个区域之后该区域温度会由指数上升和下降,使得中子星表面出现一定的温度梯度,温度梯度和中子星自旋相影响,使得暴尾部出现爆发振荡信号。第二种“规范”模型称为纬度依赖模型,因为该模型假设中子星表面的冷却时间随纬度变化而变化,在快速旋转的中子星表面,有效重力加速度沿赤道方向减小,在极点处最大,所以质量吸积率和冷却时间在赤道处最大、极点处最小。虽然这个模型考虑了冷却时间变化因素,但保持了自旋轴对称性。第三种为非对称冷却唯象模型,冷却时间变化不再满足中子星自旋轴对称性,也不考虑每个单独区域的温度变化,而是只考虑已燃烧区域(热区)和剩下区域(冷区)的温度。火焰燃烧到整个中子星表面时,时间很短(用时大概为1 s),中子星表面温度梯度可忽略,当开始冷却时,先燃烧的区域最先冷却,假设冷却与燃烧以相同方式扩散,但是时间更长(用时大概为10 s),因此先燃烧区域不仅先冷却,而且还要比其他区域冷却时间短很多,此过程可以增加暴尾部振荡信号的振幅。 精确测量中子星的质量和半径对于我们认识致密物质状态方程和中子星的内部结构非常关键。LMXB的中子星质量和半径可由多种独立的方法测定,例如,由光球层膨胀暴的时间分辨能谱[75],测量质量和半径;已知距离,拟合宁静态的热谱[76],测量质量和半径;由中子星表面特征谱线[77]测量质量和半径之比等。拟合脉冲轮廓也能测量中子星的质量和半径。X射线脉冲星因表面热斑的周期性自转,表现出脉冲形式的辐射。在强引力场的影响下,热斑的辐射束弯曲,因此观测到的脉冲轮廓,与脉冲星的质量与半径之比有关,同时也与热斑形状、与视线方向的夹角等因素相关。Psaltis等人[78]研究了快速转动的吸积X射线脉冲星(自转频率180∼620 Hz),此时星体因高速旋转而形变,明显改变了脉冲轮廓,这个效应独立于光线弯曲,给出了另外的质量半径的限制关系,从而可以不用测量源的距离而测量出星体的质量和半径。NICER卫星的主要科学目标就是研究转动供能X射线脉冲星(脉冲星的辐射能量来源于星体的自转能)的物态,目前,NICER累积了大量的观测数据,测量了PSR J0030+0451的质量和半径[79],约束了致密物质的状态,此方法依赖于辐射光子来自于中子星表面而不是中子星附近的磁层。在X射线脉冲星中应用的脉冲轮廓拟合方法也适用于XBO,高精度质量和半径的测量,需要探测到106个光子。未来eXTP卫星将会搭载大有效面积的探测器,通过测量中子星LMXB的XBO,累积多次观测得到的脉冲轮廓,拟合其脉冲轮廓,将对中子星物态给出更加严格的限制。现在已经观测到了30颗中子星LMXB(吸积供能)的XBO现象,观测样本数远远多于转动供能X射线脉冲星,更多的中子星质量和半径的测量将在更宽的质量谱上约束致密物质状态①例如,可能发现大质量(>2M⊙)或小质量(<1M⊙)的中子星,而小质量中子星的半径如能准确测量,将有效地区分奇异星和中子星模型。。结合其他方法测量到的星体质量和半径,将有效地减小质量和半径测量的误差。 多年来,人们通过多波段观测搜索各类中子星,研究中子星的自转频率分布和演化,因为它们是理解恒星和双星系统演化的重要因素。观测表明,中子星LMXB的XBO频率接近于星体的自转频率,其误差在几赫兹之内,由XBO的观测,独立发现了较多吸积供能中子星的自转频率。通过对中子星LMXB 4U 1608-522的XBO观测,发现其自转频率约为620 Hz,这是由X射线能段测量到转动最快的中子星。发现更多快速自转的中子星,对理解射电毫秒脉冲星的再循环机制至关重要,尤其是吸积演化过程中的物质转移、磁场衰减和吸积加速等过程。另外一个重要的问题是中子星可以达到的最大自转频率②是否存在亚毫秒中子星,亚毫秒周期的中子星将严格而干净地约束致密物质的状态方程[80–82]。而且,在快速旋转中子星上,是否会观测到暴振荡,这也是目前尚未解决的问题。 本文介绍了中子星LMXB中的XBO现象。人们提出了多个理论模型解释XBO的观测现象,但这些模型难以解释所有的观测现象,尚待完善。目前,NICER探测到了XBO现象,并将观测到更多的I型X射线暴。未来,中国主导的eXTP卫星也将观测中子星作为主要的研究目标,我们期待通过观测XBO,发现亚毫秒脉冲星,精确测量中子星的质量和半径并严格限制致密物质状态,建立更加完善的XBO理论模型。 附录A 搜索Swift/BAT的X射线暴数据中的X射线暴振荡信号,以Swift在2006年观测到的4U 0614+09的爆发数据为例,步骤如下。 (1)下载对应观测数据。 (2)在数据所在目录,终端输入命令 对下载数据进行处理,输出文件按类型排放在输出目录中。 (3)在目录00234849000-results/events,得到能量范围为13∼20 keV的事例数据转变为时间分辨率为0.25 s的并道数据,终端输入: 其中detmask是探测器的质量图(不良像素图)文件名称,用来反映BAT望远镜中科学数据质量的任何操作问题。 (4)画出上面得到并道数据的光变曲线,查看X射线暴所在时间范围,爆发时间内光变曲线取每4 s进行快速傅里叶变换,每次移动步长为0.25 s,得到功率密度谱,本例用Python对数据进行相关操作,具体代码见代码块(1)。 (5)得到功率密度谱在414 Hz附近有明显的振荡信号,取频率范围413∼416 Hz,13∼20 keV的事例数据,做检验,将=22,26,30,34,38的等高线投影在光变曲线上,得到的动态功率谱与Strohmayer等人[38]对4U 0614+09爆发数据的处理结果一致,振荡信号频率为414.75 Hz,具体操作见代码块(2)。2 I型X射线暴
2.1 H燃烧
2.2 He燃烧
2.3 C燃烧
3 X射线暴振荡的观测
3.1 搜索XBO的方法
3.2 如何估计真实信号概率
3.3 XBO的动态功率密度谱
3.4 RXTE对XBO的观测
3.5 中子星内部成分探索器(NICER)探测XBO
4 吸积毫秒X射线脉冲星与XBO的关系
5 XBO理论解释
5.1 热斑模型
5.2 表面模式振荡模型
5.3 冷却尾流(cooling wake)模型
6 暴振荡的应用
6.1 测量中子星的质量和半径
6.2 测量中子星自转
7 结论与展望