付妍妍12张月竹12 魏益焕1 张承民2禺少华2 潘元月2郭元旗2 王德华23

(1渤海大学数理学院物理系锦州121000)

(2中国科学院国家天文台北京100012)

(3贵州师范大学物理与电子科学学院贵阳550001)

中子星吸积加速与引力辐射影响∗

付妍妍1,2†张月竹1,2 魏益焕1 张承民2‡禺少华2 潘元月2郭元旗2 王德华2,3

(1渤海大学数理学院物理系锦州121000)

(2中国科学院国家天文台北京100012)

(3贵州师范大学物理与电子科学学院贵阳550001)

研究了中子星双星系统吸积过程中的中子星自旋加速,得到其周期随着吸积质量的变化关系.并且通过对中子星磁场与周期演化分析,对比脉冲星的观测,获得了理论与观测类似的结果.在此基础上,研究了引力辐射作用对中子星自旋加速过程的影响,并且推导了吸积过程中中子星自旋周期变化率,探讨了吸积与引力辐射影响相当的临界角速度Ωcr,评估了中子星的自旋在引力辐射下的影响.

脉冲星:普通,恒星:中子,双星:普通,X射线:双星,恒星:演化

1 引言

中子星是大质量恒星演化的终端产物之一.一般对于质量大约为8-15M⊙的恒星演化到生命末期的时候,恒星内部的核能就被耗尽,星体内部的辐射压和热压将无法与自身引力平衡,此时恒星就开始塌缩,引起超新星爆发.核心形成的致密星体依靠中子的简并压或核力与自身的引力保持平衡,此致密星体就是中子星[1].

观测发现,大部分中子星是孤立的,还有一些存在于双星系统中[2],截止到2014年, ATNF(Australia Telescope National Facility)脉冲星数据网已经记录了2 300多颗脉冲星,其中有230多颗处于中子星双星系统[3].在中子星双星系统的吸积初期,中子星的磁场极强(∼1012Gs)[4],吸积物质与中子星磁极形成X-ray源[5].随着吸积质量的增加,积累的物质通过磁极不断堆积到中子星表面,导致磁极区域面积扩大,其磁力线总数保持不变,磁力线密度不断变小,从而引起中子星磁场衰减[6].除此之外,还有许多因素会导致磁场衰减,如壳层屏蔽[7]和欧姆耗散[8−9]等.在中子星双星系统中,伴星可以为中子星提供足够的吸积质量,这导致中子星自旋不断加快,当吸积质量达到∼0.1-0.2M⊙时,其自旋周期可加速到毫秒量级(∼20 ms),即形成毫秒脉冲星[4,10].由于吸积的质量不同,中子星可能演化为周期分布为毫秒到几十毫秒的再生脉冲星(recycled pulsar)[11−12].

假设中子星球体形变率为ε,由于轻微形变导致的质量四极矩而产生引力辐射,这将消耗自转动能,对中子星的自旋产生影响[1].引力辐射会消耗吸积得到的角动量,从而限制了毫秒脉冲星自旋的最大速率.

另外,毫秒脉冲星研究是目前在建的500 m口径球面射电望远镜(FAST)的重要科学目标之一[13−14],而脉冲星与引力辐射研究也是FAST的重要关注点[15−17].因此,本文研究内容将为FAST巡天的脉冲星样本特性提供演化信息.

本文主要内容如下:第2部分介绍吸积中子星磁场,以及自旋周期演化;第3部分详细讨论了引力辐射对吸积中子星自旋周期影响,并且推导出中子星临界自旋周期Ωcr随中子星形变率ε的变化公式;文章最后一部分给出了吸积中子星演化和形成的结论及其相应的讨论.

2 吸积中子星自旋周期演化

2.1 中子星磁场衰减模型

中子星通过吸积伴星质量导致磁场衰减的理论模型由Zhang和Kojima提出(以下简称为ZK06)[18].假设中子星壳磁场是冻结的，当磁场为1012Gs时,吸积的物质会通过开放的磁力线落到磁偶极区域;随着吸积物质不断堆积增多,极冠范围扩大,磁力线的密度会被稀释,这导致磁场减弱[19−20].通过这一系列过程,当中子星所吸积的质量超过∼0.1-0.2 M⊙时,中子星的磁球半径将逐渐接近星体表面[21−23],磁场将最终达到底磁场∼108Gs,即磁场可达到最低值约为108Gs[24−25].磁场衰减公式表达为:

式中y为质量参数y=2∆M/(7Mcr),Mcr为中子星的壳质量∼0.2 M⊙,吸积质量与吸积率及吸积时间t满足关系∆M=,参数B0是初始时刻(t=0)时的中子星磁场,Bf被定义为中子星的底磁场,其物理意义是当磁球半径RM到达中子星表面时的磁场值[26−27],磁球半径的表达式为:

RA为阿尔文半径,也叫磁球特征半径,其中参数ϕ代表中子星的磁球半径和阿尔文半径的关系,一般情况下取ϕ=0.5[1,28].

这样我们得到底磁场的表达式为:

当中子星所吸积物质的质量达到约0.2 M⊙时,磁场将达到底磁场Bf[29],

2.2 吸积中子星自旋周期的变化

在中子星进行再循环演化的早期阶段中,其初始自旋周期很长,中子星就会被加速[30],但是当初始周期较短时,中子星又会被减速.中子星自旋周期P的变化与自身磁场B的关系式[28]是:

式中I45是中子星以1045g·cm2为单位的转动惯量,L37为以1037erg/s为单位的中子星光度,n(ωs)为无量纲力矩,其中ωs是快因子(fastness parameter)[29].

如果加入引力辐射之后,中子星自旋周期的变化率就会有所不同.假设中子星有着轻微偏离球形的形变,如均匀分布的椭球星体,引力辐射引起的周期变化率为:

2.3 吸积中子星的自旋周期图(B-P图)

2.3.1 平衡周期线方程(加速线方程)与死亡线方程

假设中子星磁场的对称轴与其自旋轴重合,当磁层半径处物质的频率加速到开普勒频率,我们就会得到自旋周期和磁场的对应关系,成为中子星的加速线方程,也叫做平衡周期线方程[2,31].表达式为:

式中Peq是平衡自旋周期,B9是以109Gs为单位的偶极磁场,M为中子星质量,˙M为中子星吸积率,˙MEdd为爱丁顿吸积率(1018g/s),R6是以106cm为单位的中子星半径.当式中中子星吸积率 ˙M等于爱丁顿吸积率,即 ˙M=1018g/s时,得到的平衡周期线即为爱丁顿加速线.

我们知道在中子星的磁球层中,充满了高能粒子.当极冠的电压高于一定边界值时,极冠中不断产生的正负电子对就会被加速,当极冠电压低于这个值时,这种脉冲行为就会停止.这种脉冲行为停止的界限我们可以用边界磁场和自旋周期的关系来表达:Bs/P2,根据理论预言这个边界电压值B12/P2约为0.2,其中中子星磁场B12是以1012Gs为单位,这一结果也与脉冲星的观测分布大体相符,这个方程就是脉冲星的死亡线方程[32].

2.3.2 自旋周期图(B-P图)

首先,利用ATNF pulsar catalogue的脉冲星数据[3],我们可以画出所有脉冲星的自旋周期及磁场B-P分布图,如图1,并且可以通过观察得出普通脉冲星及毫秒脉冲星相对于加速线和死亡线的分布关系.通过观察我们得出,普通脉冲星的周期一般在0.5 s左右,磁场一般在1012Gs左右;而毫秒脉冲星的周期一般在5 ms左右,其磁场一般集中在108-109Gs左右.这种分布与理论上吸积过程会使脉冲星加速、周期减小、磁场减小的模型相一致.

图1 脉冲星自旋周期及磁场分布图,数据取自ATNF脉冲星数据网.Fig.1The diagram of magnetic fi eld versus spin period,and the data are taken from the ATNF pulsar catalogue.

根据上面的中子星磁场和周期演化的基本方程,假设这一观察过程中吸积率是不变的,中子星的初始周期P0和吸积加速之后的周期P,进行数值计算,我们就可以得到相应的B-P分布图.当初始周期P0=10 s,初始磁场B0=5×1012Gs时,将吸积率=1017g/s,=1018g/s,=1019g/s 3种情况的B-P分布图放到一起比较,如图2,并且观察它们与平衡周期加速线之间的关系.当假定吸积率为 ˙M=1019g/s时,我们通常称为超爱丁顿情况下,所得出的B-P分布图位置靠近爱丁顿加速线.从图中我们可以看出,此种情况的B-P分布图在周期较大的时候,会出现自旋周期位于爱丁顿吸积加速线上方,也说明了吸积率对中子星吸积过程的影响.

在上述计算过程中假设了吸积率是不变的,但是实际情况下许多中子星双星中的吸积过程是变化的(如表现为暂现源),即使对于稳定的吸积过程,在末期也会经历吸积率逐渐减小的过程[25],这些都会影响吸积过程最终周期的大小.

图2 吸积中子星B-P演化图,其中初始周期P0=10 s,初始磁场B0=5×1012Gs,吸积率分别取=1017g/s,=1018g/s,=1019g/s,蓝实线为爱丁顿加速线(见(9)式).Fig.2 The B-P diagram for the accreting NSs with the initial periodP0=10 s,the initial magnetic fi eldB0=5×1012Gs,and the accretion rate is 1017g/s,1018g/s,and 1019g/s,respectively.The blue solid line is the Eddington spin-up line(see Eq.(9)).

在吸积过程中,中子星自旋周期和磁场的变化还可以通过数据画出.当吸积率给定时,如=1018g/s,中子星自旋周期P与磁场B随吸积质量∆M的变化关系,如图3.通过图3我们可以观察出,吸积中子星的周期随着吸积质量的增大,从初始的10 s减小到10 ms左右,说明星体的自旋在加速;同时中子星自身的磁场也从1012Gs减小到108Gs左右.B-P演化轨迹与观测的分布是相符合的.

柳红头上蒙着一条红丝巾。苏石从她的头上剥了下来，红丝巾打着一个结，是个死结。他捏着红丝巾的那个结，硬梆梆的。苏石从大棚里爬了出来，瘫坐在地上，点上一支烟，他边抽烟边借着大棚里透出朦胧的灯光，审视着手上的红丝巾，好像不认识似的。

当中子星吸积加速过程中加入引力辐射的影响,在吸积演化的晚期,中子星最后能达到的最小自旋周期将会受到引力辐射的影响和制约.由于引力辐射的减速效应,会使毫秒脉冲星的最小自旋周期变长.当中子星吸积率分别取=1018g/s和=1019g/s,中子星形变率分别取ε=10−4和ε=10−5时,画出正常吸积过程与加入引力辐射的吸积过程中中子星自旋周期P与磁场B的演化对比图,如图4所示.

从图4中可以看出,在正常吸积过程中,中子星可以达到的最小自旋周期为1.74 ms.而加入引力辐射的吸积过程中,当中子星形变率ε=10−4时,中子星达到的最小自旋平衡周期为44 ms,此时周期变化率==3.2×10−15s/s;当中子星形变率ε=10−5时,中子星达到的最小自旋平衡周期为17.5 ms,此时周期变化率==5.1×10−16s/s.可见引力辐射在中子星吸积加速过程的后期不仅使中子星自旋周期的减小速度减慢,同时也会影响延长中子星能达到的最小自旋周期.

根据目前所观测到的脉冲星周期分布的拟合来看,孤立的毫秒射电脉冲星与双星系统中的毫秒脉冲星的周期分布大体一致,稍微有点区别.孤立的毫秒射电脉冲星的周期分布峰值为4.7 ms左右,而双星中的毫秒脉冲星周期分布峰值为3.5 ms左右.一般认为,孤立毫秒脉冲星的形成也是来自双星吸积,但是其伴星在演化晚期的轨道较小,由于太靠近毫秒脉冲星,导致辐射蒸发伴星[33−34].

图3 中子星磁场(周期)相对吸积质量图,图中吸积率=1018g/s.Fig.3 The diagram of magnetic fi eld(period)versus accretion mass,in which the accretion rate is1018g/s

图4 吸积中子星B-P演化图,演化曲线考虑引力辐射效应,其中初始周期P0=10 s,初始磁场B0=5×1012Gs.图中虚线代表=1019g/s,点虚线代表=1018g/s,黑色实线代表=1018g/s,ε=10−5,点点虚线代表=1018g/s,ε=10−4.Fig.4 The B-P diagram for the accreting NSs with the initial periodP0=10 s and the initial magnetic if eldB0=5×1012Gs,in which the e ff ect of gravitational radiation is taken into account.In the B-P diagram,the dashed line stands for=1019g/s,the dot-dashed line stands for=1018g/s,the black solid line stands for=1018g/s,ε=10−5,and the dot-dot-dashed line stands forM˙=1018g/s,ε=10−4.

3 引力辐射对吸积中子星周期的影响

中子星自旋周期的变化不仅受吸积过程影响,还会受到引力辐射的影响.假设中子星是一种有着轻微形变,并且均匀分布的椭球星体,引力辐射引起的能量变化率为[1]:

可以得到动能损失率星体减速关系:

在中子星双星系统中,吸积加速率由动量守恒给出:

式中Rs为史瓦西半径,且Rs=2GM/c2,一般取∼4.2 km.

当引力波引起的减速率与吸积加速率相当时,中子星的自转速度达到平衡值.所以我们可知由(10)式和(12)式,可得

即为中子星临界角速度Ωcr随中子星形变率ε的变化公式,其中ε−5是以10−5为单位.

在中子星吸积加速的过程中,当自旋周期减小到10 ms左右的时候,根据(14)式可知,中子星形变率ε∼10−5,引力辐射对中子星自旋周期变化的影响开始明显.这意味着,引力波减速效应只有当快速转动才起作用.

中子星形变率与强磁场存在着一定的正比关系,当中子星磁场强度达到1015Gs时,形变率ε会达到10−4[35].但是据目前中子星的实际观测情况来看,中子星磁场一般不会达到1015Gs,所以形变率的数值不会高于理论值.此外,中子星形变率与其他物理量之间的关系还有待于探索与研究.

当ε=10−5时,引力辐射引起的能量变化率=1.76×1021(erg/s)Ω6;吸积引起的能量变化率=1.19×1035(erg/s)Ω.为了研究两者的加速和减速效果,我们将两种情况画到同一个图中进行观察.

当Ωcr定义域取10-10000时,所得的图像如图5所示.图中的交点为引力辐射引起的能量变化率与正常吸积引起的能量变化率相等的情况,即==6.4×1037erg/s,此时Ωcr=583 rad/s,所对应的周期P=1.08×10−4s,即10.8 ms.

从图5中可以看出,当中子星的自旋加速到周期为P=10.8 ms之前,引力辐射对中子星周期改变的影响不太大;当中子星的自旋周期减小到P=10.8 ms以后,引力辐射开始对中子星周期改变起到主导作用,而中子星吸积过程对其周期变化的影响开始慢慢减弱.

图5 中子星能量变化率相对自旋角速度图,其中中子星形变率ε=10−5.图中虚线为吸积过程中的能量变化率(见(12)式),实线为考虑引力辐射效应时的能量变化率(见(10)式).Ωcr是以103rad/s为单位的角速度.Fig.5 The diagram of energy gradient versus angular velocity,in which the neutron star deformation rate is ε=10−5.The dashed line stands for the accretion process only(see Eq.(12)),and the solid line stands for the case when the e ff ect of gravitational radiation is considered(see Eq.(10)).Ωcris the palstance in units of 103rad/s.

4 结论及其讨论

在中子星吸积过程中,我们考虑吸积加速与引力波引起辐射减速两种情况,研究对中子星自旋周期变化的影响.我们推导了两种情况引起的转动动能变化率相同时,中子星临界角速度Ωcr随中子星自身形变率ε的变化公式(见(14)式).当中子星初始周期P0=10 s,初始磁场B0=5×1012Gs时,加速到达周期约10 ms,可得到中子星的形变率ε在10−5量级左右.只有在快速转动时(周期低于10 ms),引力辐射对中子星自旋周期变化的影响开始明显,这可以解释为什么毫秒脉冲星的自旋周期没有加速到亚毫秒量级.

当ε=10−5时,中子星自旋周期P=10.8 ms,自旋角速度Ωcr=583 rad/s,此时能量变化率==6.4×1037erg/s.说明在自旋周期减小到10.8 ms之前,中子星吸积过程对其自身周期的改变起主导作用;当吸积过程中的能量变化率与引力辐射引起的能量变化率相等时,引力辐射作用开始明显增加,使中子星自旋周期加速缓慢.

在本文的计算中,中子星质量假定为1.4 M⊙,这个中子星质量值是针对46颗具有可靠测量的中子星系统的统计而得到的平均值,其中包括射电脉冲双星系统和X射线中子星系统,这一统计结果也与许多理论值相吻合[36−37].

此外,一般而言双星系统中的脉冲星通过吸积加速过程,最终都会演化到爱丁顿平衡周期线的下方,但是根据ATNF提供的脉冲星双星系统的数据观察会发现,也有一些双星脉冲星位于爱丁顿平衡周期线的上方.我们认为,这可能是由于存在高吸积率的脉冲星双星系统所导致的.计算表明,在超爱丁顿情况下(吸积率=1019g/s)自旋周期及磁场演化位置可以出现在B-P图上方(见图2).

[1]Shapiro S L,Teukolsky S A.Black Holes,White Dwarfs,and Neutron Stars.New York:John Wiley, 1983

[2]Bhattacharya D,van den Heuvel E P J.PhR,1991,203:1

[3]Manchester R N,Hobbs G,Teoh A,et al.AJ,2005,29:1993

[4]Alpar M A,Cheng A F,Ruderman M A,et al.Nature,1982,300:728

[5]Liu Q Z,van Paradijs J,van den Heuvel E P J.A&A,2007,469:807

[6]Cheng K S,Zhang C M.A&A,1998,337:441

[7]Zhang C M.A&A,1998,330:195

[8]Cheng K S,Zhang C M.A&A,2000,361:1001

[9]Geppert U,Rheinhaedt M,Zannias T.A&A,2003,412:L33

[10]Burderi L,D’Antona N,Burgay M.ApJ,2002,574:325

[11]Camilo F,Thorsett S E,Kulkarni S R.ApJ,1994,421:L15

[12]Kaspi V M.PNAS,2010,107:7147

[13]Nan R D.ScChG,2006,129:148

[14]Smits R,Lorimer D R,Kramer M,et al.A&A,2009,505:919

[15]Yue Y L,Li D,Nan R D.IAUS,2013,291:577

[16]Li D,Nan R D,Pan Z C.IAUS,2013,291:325

[17]Nan R D,Li D,Jin C J,et al.IJMPD,2011,20:989

[18]Zhang C M,Kojima Y.MNRAS,2006,336:137

[19]Melatos A,Phinney E S.PASA,2001,18:421

[20]Melatos A,Payne D.ApJ，2005,623:1044

[21]Li X D,Wang Z R.A&A,1996,307:5

[22]Li X D,Wang Z R.ApJ,1999,513:845

[23]Zhang C M,Wang J,Zhao Y H,et al.A&A,2011,527:83

[24]Phinney E S,Kulkarni S R.ARA&A,1994,32:591

[25]Tauris M,Langer N,Kramer M.MNRAS,2012,425:1601

[26]Taam R E,van den Heuvel E P J.ApJ,1986,305:235

[27]Tauris M,Langer N,Kramer M.MNRAS,2011,416:2130

[28]Ghosh P,Lamb F K.ApJ,1979,234:296

[29]Pan Y Y,Wang N,Zhang C M.Ap&SS,2013,346:119

[30]Wang J,Zhang C M,Zhao Y H,et al.A&A,2011,526:88

[31]van den Heuvel E P J.NYASA,1977,291:39

[32]Ruderman M A,Sutherland P G.ApJ,1975,196:51

[33]蔡炎,阿里塔尼,赵永恒,等.天文学报,2011,52:449

[34]Cai Y,Ali Taani,Zhao Y H,et al.ChA&A,2012,36:137

[35]Bonazzola S,Marck A.A&A,1994,44:655

[36]程争,张承民,赵永恒,等.天文学报,2013,54:514

[37]Cheng Z,Zhang C M,Zhao Y H,et al.ChA&A,2014,38:294

Accreting Neutron Star Acceleration and E ff ects of Gravitational Radiation

FU Yan-yan1,2ZHANG Yue-zhu1,2WEI Yi-huan1ZHANG Cheng-min2YU Shao-hua2PAN Yuan-yue2GUO Yuan-qi2WANG De-hua2,3
(1 Department of Physics,School of Mathematics and Physics,Bohai University,Jinzhou 121000)
(2 National Astronomical Observatories,Chinese Academy of Sciences,Beijing 100012)
(3 School of Physics and Electronic Science,Guizhou Normal University,Guiyang 550001)

In this paper we study the neutron star’s spin acceleration in binary systems,and the relation how the spin period changes with the accreted mass.We analyze the evolution of magnetic fi elds and the spin periods of neutron stars,and compare the modeled results with the observational data of pulsars,to show that they are consistent with each other.Considering the e ff ect of gravitational radiation on the spin-up process of neutron stars,we derive the changing rate of the neutron star’s spin periods in the accretion process.We estimate the critical angular velocity Ωcr,at which the accretion torque is balanced by that of gravitational radiation,and discuss the in fl uence of gravitational radiation on neutron star’s spin evolution.

pulsars:general,stars:neutron,binaries:general,X-rays:binaries,stars: evolution

P145;

A

10.15940/j.cnki.0001-5245.2015.03.003

2014-11-24收到原稿,2015-01-14收到修改稿

∗国家自然科学基金项目(11173034)及国家重点基础研究发展计划(2012CB821800)资助†fychao@bao.ac.cn

‡zhangcm@bao.ac.cn