超强磁场下中子星壳层的电导率和磁星环向磁场欧姆衰变*

陈建玲 王辉 贾焕玉 马紫微 李永宏 谭俊

1) (运城学院物理与电子工程系,运城 044000)

2) (西南交通大学物理科学与技术学院,成都 610031)

3) (运城学院数学与信息技术学院,运城 044000)

1 引 言

中子星是研究致密天体物理的最佳实验室,它们有很强的磁场.已知的中子星磁场或者是相对稳定的、或者随时间缓慢地变化[1,2].磁场在中子星自旋减慢演化及磁层活动性中发挥重要作用,因此得到了广泛的研究[3−7].

磁星是指主要由磁场提供辐射能量的一类脉冲星.对磁星的观测和理论研究是当前脉冲星领域一个重要的热点[2,3].磁星大致分为软伽玛射线重复爆(soft Gamma-ray repeaters,SGRs)和反常X射线脉冲星(anomalous X-ray pulsars,AXPs),被发现的磁星及其候选体数量已经增长到约29个.磁星表面偶极磁场可以由转动周期和一阶周期导数来估计,其量级约为1014—1015G (1 G=10–4T)[3−7].磁星表现出一系列广泛的X射线活动性,包括短暴、长暴、耀斑和准周期振荡,通常伴随着一系列有趣的到达时间行为,如增大的自转减慢、周期跃变和反周期跃变[8−13].这些行为一般被解释为超强磁场的衰变,或者是由于表面磁张压增加导致磁星壳层破裂,这会加剧磁层的扭曲[14].近年来,国内不少专家对中子星(包括磁星)的磁场演化、冷却及辐射机制做了大量的、深入的研究[15],取得了丰富的成果,但是涉及到磁星的环向磁场衰变和表面热辐射方面的研究很少.

我们知道,具有纯粹的极向磁场构型的中子星是极其不稳定的,同样地,由于泰勒不稳定性[16],具有纯粹的环向磁场构型的中子星也是不稳定的.而一个稳定的中子星磁场结构要求在中子星内部具有极向磁场和环向磁场组成的混合磁场,并且壳层磁场表现为偶极磁场(占主导)和高阶多极磁场(对壳层磁场贡献较小)[17].在霍尔漂移(Hall drift)的作用下,一个纯粹的极向磁场首先会在星体内部产生一个环向磁场,然后与环向磁场紧密相关的极向电流将极向磁场推压到磁极冠区,从而增大磁偶极矩[18,19].一个新诞生的中子星(原中子星)在壳层形成之前由于对流形成包含多级成分的极向磁场[20].在核塌缩后,由于较差自转极向磁场不断地被扭曲,因此,中子星磁场可能存在一个很强的环向成分,极向磁偶极矩在径向较差自转的作用下会产生环向的磁四极场[21].

以往关于中子星磁场的欧姆衰变(又称欧姆耗散)的研究主要集中于对普通中子星壳层磁场衰变的研究,这包括了对壳层磁场欧姆衰变本征模的计算[22]、星体表面磁场欧姆衰变的自相似解[23]、欧姆耗散方程[24,25]以及由于多极磁场的欧姆衰变导致的磁能湮灭率[26].1994年Geppert和Urpin[27]首次研究了在吸积中子星中也会由于欧姆衰变导致磁场演化.有两个因素可以减缓磁场的衰减: 引力红移效应以及空间内禀的弯曲几何[28].由于广义相对论效应,磁场衰变时标会增加,但将与平直时空中磁场衰变时标保持相同的量级[29−32].最近,王辉、高志福、王娜等[33](以下称WGW19)在广义相对论框架下推导出普通中子星磁场欧姆衰变的本征值方程.他们发现: 靠旋转供能的脉冲星PSR J1640-4631内部可能发生环向磁场的欧姆衰变,但是释放的磁场能不足以提供持续的软X射线光度,并且讨论了该源可能的、各向异性的软X射线辐射机制.

WGW19虽然推导出在广义相对论框架下旋转供能脉冲星的偶极环向磁场欧姆衰变本征方程,但是磁场通过欧姆衰变再以焦耳热的方式释放热能这一现象普遍地存在于包括磁星在内的年轻的中子星内部,WGW19推导出的方程可能同样地适用于磁星模型,主要的原因分析如下.1)对于包括磁星在内的强磁化的中子星,磁场能以欧姆耗散焦耳加热在其演化的早期非常重要,WGW19推导出的偶极环向磁场欧姆衰变本征方程是基于极向磁场强度不低于约1013G、温度不高于约108K的中子星的模型,该模型同样地适用于磁星.2)相比于依靠旋转供能的脉冲星(简称普通中子星),磁星理论上可能具更高的内部多极磁场、扭曲的磁层和几种不同的磁场起源机制[20],但是观测上没有证明磁星的磁场结构位型、演化方式与中子星之间存在本质上的差别,即磁星和普通中子星应当同样地具有极向磁场和环向磁场分量、同样地拥有类似的边值条件.3)由于较高的磁能密度贡献,磁星的物态方程可能比普通中子星的稍硬、质量稍大.这个物态方程上的理论假设并没有得到观测上的支持.在描述磁星星体结构时,通常认为: 磁星和普通中子星一样具有等离子体的磁层、固态的壳层和流体的核,核内具有核子超流与质子超导.在壳层中,电阻率(电导率)主要是由于电子-声子和电子-杂质散射过程[23−27],导致比在流体内部更有效的欧姆耗散.4)按照电动力学,中子星内部可能存在的径向超导电流能产生环向磁场,由于核内超导与外部的真空,磁场能的湮灭最为有效方式是: 在具有电阻的中子星的壳层以欧姆耗散的方式产生焦耳热,欧姆耗散率由组成物质的有限电导率决定[33].磁场能以欧姆耗散的方式产生焦耳热的情况只会发生在具有电阻的壳层.5) WGW19以广义的麦克斯韦方程组为出发点,采用一个球对称的中子星模型和时空几何,并假定一个无力场(force-free)磁场位型,在该模型下磁场被限制在壳层,即核内由于超流排斥作用而不存在磁场.本文假定磁星具有无力场磁场位型,并采用球对称的中子星模型和时空几何.我们有理由相信,WGW19在广义相对论下推导出的偶极环向磁场欧姆衰变本征方程也同样地用于计算磁星的欧姆衰变率、磁场能释放率以及软X射线光度.

观测表明: 部分宁静状态下的的磁星存在持续的软X射线辐射,其光度的典型值LX为1034—1036erg/s (1 erg/s=10–7J/s).磁星在爆发过程中[34],软X射线光度急剧增加,可达1038erg/s或更高; 在外爆衰退期,磁星会出现热斑缩小、射电辐射、X射线谱变软及脉冲轮廓变得简单,这些与磁星磁层扭曲程度和磁层电流发生变化有关.软X射线光子可能存在热起源,可以用黑体谱(幂律+黑体 或者双黑体谱)来拟合,对应的温度约为2×106—6×106K,这比依靠转动能损提供能量的脉冲星的典型温度值高很多,这可能是由于磁场衰变提供额外的加热[30].Kaspi和Beloborodov[9]讨论了磁星的磁场欧姆衰变和软X射线光度之间的关系,他们认为: 欧姆衰变可以提供磁星软X射线辐射所需要的热能,首先,相比于磁星的年龄,磁场的欧姆衰变时标相当的短,这就要求磁场在一个小尺度范围内发生变化; 其次,磁场变化幅度δB≈1016G.Beloborodov和Li[31]也认为磁星软X射线光度可能与磁场的欧姆衰变有关,但他们没有在超强磁场中进行电导率的计算,也没有考虑广义相对论效应对环形磁场的欧姆衰变的影响.

本文将在第2节给出广义相对论下霍尔感应方程和核物质状态方程; 第3节计算在超强磁场下壳层的电导率、环向磁场衰变率及磁能衰变率;第4节讨论磁星的旋转能损率Lrot、环向磁场衰变磁能释放率LB及软X射线光度LX之间的关系以及软X-射线可能的各向异性的起源,第5节进行总结与讨论.

2 广义相对论下的霍尔感应方程和核物质状态方程

2.1 霍尔感应方程

为方便起见,我们假设一个球对称的中子星模型和时空几何,考虑到广义相对论效应,中子星的时空几何形状可以由下式描述

其中x0=ct,Vij是空间坐标xi的函数,Φ(r) 是引力势[32].结合磁流体动力学平衡与欧姆定律,我们得到广义相对论下中子星磁场霍尔感应方程

其中s代表导电率,eijk代表完全反对称的Levi-Civata张量的分量,D表示变量微分算子.在非奇异的、静态的、球对称几何背景下,(1)式中的空间部分的坐标可以用下式表示:

其中hi=hix1,x2,x3是标度因子,它们分别为

其中m(r)=GM(r)/c2,M(r) 是半径为r的球体内部物质总质量.消除(2)式中磁场的坐标分量,得到更为简化的感应方程

其中Z=eΦ≡(1−2M(R)/R)1/2是红移因子[29,33].

2.2 磁场的三分量

轴对称条件下,中子星的磁场B可以分解成极向磁场分量Bp=Brer+Bθeθ与环形磁场分量Bt=Bφeφ两个部分,其中er,eθ和eφ分别代表r方向、q方向和φ方向的单位矢量.在广义相对论下磁场三分量由下式给出:

这里B是中子星磁极冠区表面处偶极极向磁场(与l=1,m=0相对应);Fl,Tl和Hl是待定函数,Ylm是l阶勒让德多项式.为了简易起见,我们仅考虑极向磁场的偶极分量,即l=1,m=0,于是得到 dYlm/dθ=−sinθ与 Ym=cosθ,再利用分离分量法[29],很容易得到

以下为了方便表达,我们一律省去了脚标1,于是(6)式简化为

这里r为无量纲的距离.以下星体半径表示为RNS→R.为了求解F(r,t) 需要两个边界条件: 磁场局限于从壳层表面到壳核边界处,或是局限于从核半径Rcore到星体半径R处.磁场的内部边界条件: 当r→Rcore时,F(r,t)→0.其次,在中子星表面处,壳层与真空中磁场是连续的,得到磁场的内部边界条件

其中y=R/RS,RS≡2GM/c2是史瓦西半径.首先,由于中子星的核可能具有超导特性,核内磁场存在的可能性非常小,在无力场磁场位型下磁场满足

其中µ是与磁场曲率有关的一个参数,被解释为斯托克斯函数F(r,t) 的波数,µ的值取决于中子星的核物质物态方程.为了简单起见,对于(10)式中第一个等式,我们仅考虑µ为常数的解,这样第二个等式就会自动地满足.为了得到(10)式第一等式中内部磁场的一个通解,我们选择限制条件:

保证这个解能满足B在r和θ方向的两个分量之间平衡.当B被限制在壳层中调整µ使r分量在壳-核边界处消失,得到

其中Rcore表示壳-核边界处到球心距离.另外还有在两种中子星磁场位型: TC1与TC2[28],在这两种磁场结构位型中,环形磁场都仅仅被限制在中子星壳层内部,而在无力场磁场位型下,环形磁场可以穿越壳层渗透到星体外部.利用勒让德函数及球形贝塞尔函数,可以求解在无力场磁场位型下的高阶环形磁场的欧姆衰变本征方程.

2.3 中子星核物质状态方程

为了得到参数µ的值,需要采用较为实际的中子星结构参量,即采用较为实际的核物质状态方程.预计的中子星最大质量是依赖于不同模型下的核物质状态方程[34,35].中子星的最大质量与最小质量的数值都具有争议性.观测的中子星质量约在(1—2)M⊙的范围内(http://www.stellarcollapse.org/).尽管从状态方程的角度来看,中子星的质量可以小于一倍M⊙,但很难从超新星爆发的机理来解释其成因.在外壳层区域,我们采用Baym-Pethick-Sutherland物态方程[36]; 在内壳层区域,我们采用可压缩的液滴模型,即Baym-Bethe-Pethich物态方程(BBP模型)[37]; 在中子星的核内,我们采用相对论平均场(relative mean-field,RMF)理论模型,RMF是在具有修正效应的有效耦合常数的基础上建立起来的,是研究有限核的一个标准方法.作为对比,我们选用三个具有代表性的RMF模型——NL3[38],GM1[39]和TMA[40,41].

表1中r0表示饱和核密度;E0,K0,m*,K′,J,分别表示对称核物质饱和密度下的结合能、不可压缩系数、无量纲的有效核子质量、偏斜系数、对称能、对称能的坡度、对称能曲率、对称能偏斜系数和与体积相关的同位旋不可压缩系数.TMA模型在高密度下的物态方程与Dirac-Brueckner-Hartree-Fock 理论所预测的行为一致,因此TMA是当前最为成功的RMF参数组之一.从表1可以看出,相比于GM1和NL3模型,TMA模型具有较软的对称能,因此给定中子星最大质量值也较小.图1所示为在NL3,GM1和TMA模型中中子星的质量和半径的关系.NL3 模型给出中子星最大质量Mmax=2.78M⊙,GM1 模型给出中子星最大质量Mmax=2.45M⊙,而TMA模型给出最大质量Mmax=2.032M⊙,这个值非常接近目前观测到的中子星的最大质量[42],因此,TMA模型较为实际和可靠.在以下计算中,我们将一律采用TMA模型的物态方程.采用由文献[7]中更为恰当的M-R-I关系近似表达式,

表1 在NL3,GM1和TMA模型下饱和核物质特性.Table 1.Saturation properties of nuclear matter in the parameterizations for NL3,GM1 and TMA models.

图1 在NL3,GM1和TMA模型下中子星的质量和半径的关系Fig.1.Relationships between mass and radius of neutron stars in NL3,GM1 and TMA model.

图2 在TMA模型中磁星的转动惯量I随质量m和半径R的关系Fig.2.Relationship of moment of inertial I to mass M and radius R for magnetars in TMA models.

将内壳层的厚度近似地看成为中子星总的壳层厚度Rc=0.97 km (BBP模型),结合TMA模型和(9),(10)式,可以计算在TMA模型中对应任意最大质量的磁星的µ和I的值,并将计算结果列于表2.尽管磁星的I随质量的增加而增大,但是µ的值比较稳定(µ=1.67—1.68).为了比较在无力磁场结构位型下弯曲时空中磁场三分量之间的关系,我们选择一个典型质量M=1.45M⊙的磁星,TMA模型中,对应µ=1.676及x=r/R=0.917—1.0.利用(8)—(13)式,得到归一化的磁场三分量:Br/(Bcosθ) ,Bθ/(Bsinθ) 及Bϕ/(Bsinϕ)随x的变化,其中B为在中子星磁极表面处(x=1)极向偶极磁场.如图3所示,Br与Bj都随x的增加而增大,但是后者增加比前者快得多; 在星体表面处,Bϕ/(Bsinϕ) 达到极大值,这时Bj的强度高于Br约一个量级.

表2 在TMA模型中磁星的m,R,Rcore/R,µ和I的部分值Table 2.Partial values of m,R,Rcore/R,µ and I for magnetars in TMA model.

图3 在无力磁场结构位型下壳层归一化磁场分量Br/(Bcosθ)(红线),Bθ/(Bsinθ) (蓝线),及Bϕ/(Bsinϕ)(黄线)与归一化径向坐标x的关系(选取µ=1.676,对应在TMA模型下的M=1.45M⊙,R=11.77 km及I=1.45×1045 g·cm2)Fig.3.Normalized magnetic field components of the crustal confined for the force-free field: Br/(Bcosθ) (red line),Bθ/(Bsinθ)(blue line),and Bϕ/(Bsinϕ) (yellow line) vs.normalized radial coordinate x.Here we assume the parameter µ=1.676,corresponding to M=1.45M⊙,R=11.77 km and I=1.45×1045 g·cm2 in the TMA model.

3 磁星环向磁场欧姆衰变

3.1 环向磁场的本征方程

WGW19为了研究旋转供能脉冲星PSR J1640-4631壳层磁场演化,推导出在广义相对论框架下偶极环向磁场欧姆衰变本征方程,引入一个矢量因此,(5)式中B的j分量可以写为

对于任意一个矢量A,它的旋度服从以下规律

于是得到K的旋度在j,q和r方向的分量,并利用分离变量法,得到.

通过比较(16)式与文献[29]中(3.9)式,WGW19发现H(r,t) 与F(r,t) 有相同的演化方式.结合(8),(10)和(11)式,得到

由(17)式可知,在广义相对论框架下,磁场三分量均可以由标量的斯特克斯函数F(r,t) 来表示.磁场在φ方向和r方向上的强度之比值关系为

这里R是无量纲的中子星半径.通过耦合爱因斯坦方程与流体的能量-动量-张量,得到关于m(r) 和Φ(r)的微分方程以及流体静力学的方程P(r)[29].将m(r) ,Φ(r)和P(r) 代入文献[29]中方程(3.9),得到

(19)式必须需要满足相对论的斯托克斯流函数F(t,r).为了解出方程(19),将F(t,r) 展开成多项式形式

其中n=1,2,···,x=r/R; ∑是对与所有本征模 求和;Xn(x) 满足边值条件

L表示本征算符[29].通过引入一阶球形贝塞尔函数,得到

和弯曲时空下环向磁场本征方程

作为对比,文献[28]忽略了广义相对论效应,通过引入标量函数在平直空间将磁场的三分量表示如下:

由以上分析可知,在无力磁场位型下,由于磁场在 (r,θ) 方向必须满足相同的平衡条件和在(r,φ)方向满足相同的限制条件,在平直空间和弯曲空间中,磁场的三个分量都可以由一个分量函数来表示,并且环向磁场和极向磁场的强度之比、磁能密度之比是相同的.尽管如此,由于广义相对论效应,平直空间中的(x,t) 和弯曲空间中的F(r,t) 表达形式和演化并不相同.为了探索出磁场的演化形式,必须要确定电导率s的限制范围.

3.2 强磁场下中子星壳层的电导率

壳层内导电率分别来源于电子-声子散射与电子-杂质散射对电导率的贡献.电导率大小非常依赖温度T与密度r,后者跨越6个或6个以上数量级.决定导电率s的第三个参量是介质的不纯净度,由于在低密度的外壳层区域电导率非常低,磁场欧姆衰变时标特别短(τOhm≈ 10—102a),因此外壳层将不再考虑.根据BBP模型,内壳层密度范围约为: 从中子滴出密度ρ∼4.66×1011g·cm−3到壳-核边界ρ∼1.30×1014g·cm−3,由于随着核子数Z与质量数A单调而任意地增加,BBP模型的物态方程在较高密度ρ≥1.72×1014g·cm−3遭到Shapiro和Teukolsky[43]的质疑,我们将停止在相应的较高密度区的计算.在超强磁场下,壳层中重子数密度nB可以被看成是一个不变的量,电子丰度Ye=ZnN/nB,其中nN是原子核数密度,则电子数密度ne=nBYe=ZnN.本文将利用由文献[44]开发的、较为实际的电导率的程序(公开下载的网页为http://www.ioffe.ru/astro/conduct),并且结合BBP模型,来计算磁星壳层的电导率.由于文献[44]考虑到强磁场因素,较之以往的计算电导率程序,这个程序提供的结果更能反映中子星壳层电导率的真实情形.

我们计算出在BBP模型中不同的T和Q下磁星壳层电导率,部分计算结果列于表3.表3上半部分对应Bp=5.0×1014G,表3下半部分对应Bp=3.0×1015G.可以看出,在给定T,Q和Bp的情况下,s随着r增加而增大; 在给定T,r和Bp的情况下,s随着Q增加而减小; 在给定Q,r和Bp的情况下,s随着T增加而减小; 在给定Q,r和T的情况下,s随着Bp增加而增大,磁场在低密度区域对s的影响比在高密度区影响大,在低密度区域,s相对增加率 ∆σ/σ⩽5% ; 在高密度区域,由于 ∆σ/σ⩽10−3,磁场对s的影响忽略不计.为了简单起见,在表3中我们假定不纯净度参数Q是一个不变的量,这与中子星壳层的实际情形可能相差很大,必须给出一个合理的Q值范围.

我们知道,在较高的温度和较低的密度下,晶格声子限制了电子的运动,因此热和电荷的输运由电子-声子散射(碰撞)主导,而在高密度环境下,热和电荷的输运由电子-杂质散射主导.文献[16]对于以电子-声子散射主导电荷输运的低密度层进行研究,给出一个小的Q值范围:Q∼ 0.001—0.1;最近,对于以电子-杂质散射主导的更深的高密度层的研究预示不纯净度参数更高:Q> 1.通过对中子星磁热演化的研究,文献[45]给出更高不纯净度参数:Q> 1.通过对中子星磁热演化的研究,文献[46]给出中子星内壳层不纯净度的范围:Q∼1—100; 磁星内壳层温度比星体表面温度高出1—2个量级,但是最高温度不能超过各异性中子超流的临界温度,我们选取一个比表3更为合理的磁星壳层温度范围T∼5.0×107— 2.0×108K ,在较低密度层ρ∼4.66×1011— 1.0×1013g·cm−3对应Q∼ 0.001—0.1; 在高密度层ρ∼1.0×1013—1.30×1014g·cm−3对应Q∼ 1—25; 选取一个典型超强磁场Bp=5.0×1015G,给出在一定磁场下电导率随温度和不纯净度的变化关系图,如图4所示.由于r,T及Q来决定s的值确实存在很大的不确定性,为了方便计算,我们给出合理的参数范围: 对于低密度层电导率由电子-声子散射主导,选取Q∼ 0.1,对于较高密度层ρ∼1.0×1013—5.0×1013g·cm−3电导率由电子-杂质散射主导,选取Q∼ 1,对于更高密度层电导率仍由电子-杂质散射主导,Q值可能增大,选取Q∼ 2; 对应壳层温度稍低的磁星,T=6.0×107K,得到壳层s的变化范围s∼ 8.65×1023—8.75×1024s–1; 对应壳层温度稍高的磁星,T=2.0×108K,得到壳层s的变化范围s∼ 1.09×1023—2.52×1024s–1.

3.3 超强磁场下的磁能释放率

选择典型质量为M=1.45M⊙的磁星,并且采用BBP 和TMA物态方程,对应R=11.77 km及µ=1.676.将(22)和(23)式进行求导,得到环向磁场和极向磁场的变化率:

表3 在不同温度和不同纯净度参数下磁星壳层电导率的部分值(采用BBP模型)Table 3.Partial values of electrical conductivity for different temperatures and impurity parameters in the crust of magnetars.Here we use the equation of station (EOS) of BBP model.

图4 磁星壳层电导率随密度、温度及不纯净度参数的变化 (a)电导率由电子-声子散射主导; (b)电导率由电子-杂质散射主导; 物态方程一律采用BBP 模型Fig.4.Relationship of s to r,Τ and Q in the inner crust for magnetar: (a) The conductivity due to electron-phonon scattering; (b) the conductivity due to electron-impurity scattering.The EOS of BBP model is used.

磁场能释放率也是时间的函数,可以用下式来进行估算:

其中dV=4πr2dr,Rc=0.98 km.磁能释放率主要是由环向磁场所主导:

选择一个典型的强磁场Bp(0)=2.0×1015G,计算了磁场的衰变率和磁能释放率,部分结果列在表4中.

将σ=8.75×1024s−1和σ=2.52×1024s−1分别代入方程组(22),(23),(26),(27)中,用数值模拟方法得到在Bp(0)=3.0×1015G和Bp(0)=5.0×1014G两种情况下Bp,dBp/dt,Lp,Bt,dBt/dt,Lt和LB随时间的变化,如图5所示.由图5可以看出,Bp和Bt都同样地经历缓慢衰变和快速衰变的过程.

4 与磁星的观测对比

4.1 磁星旋转能损率和软X射线光度

从McGill磁星数据网站(http://www.physics.mcgill.ca/~pulsar/magnetar/main.html)上,可以看到有29颗磁星及其候选体,其中22颗源:CXOU J164710.2–455216 (缩写CXOU J164710),CXOU J0100430.1–721134 (缩写CXOU J01004),1RXS J170849.0–400910 (缩写1RXS J170849),1E 2259+586,1E 1048.1–5937,1E 1841–045,4U 0142+61,SGR 0418+5729,SGR 0526–66,SGR 1900+14,SGR 1806–20,1E 1547.0–5408,XTE J1810–197,CXOU J171405.7–381031 (缩写CXOU J171405),SGR 1627–41,Swift J1822–1606,Swift J1834.9–0864,SGR J1745–2900,PSRJ1622–4950,PSR J1846–025,暂变源AX J1845.0–0258及3XMM J185246.6+003377 (缩写3XMM J185246),它们具有软X射线光子.周期约为7 s的X射线的源AX J1845.0–0300,与超新星遗迹G29.6+0.1成协,X射线亮度较低,光谱较软,表明它可能是制动的X射线脉冲星[47].然而,由于缺乏精确稳定的周期导数值,人们无法对其旋转能量损失率进行评估.因此AX J1845.0–0258的磁场衰变、旋转能损率与软X-射线光度的关系将不再被考虑.3XMM J185246的旋转周期P=11.5587 s,周期导数的上限假定一个经典的偶极制动模型,给出偶极磁场的上限1013G及旋转能损率上限Lrot< 4.75×1030erg/s,该源与超新星遗迹Kes 73成协,超新星遗迹(supernova remnant,SNR)年龄4.7 ka,宁静状态下的软X射线光度上限为:初始磁场可以选取与SGR 0418–5729相同,因为它和SGR 0418–5729一样,属于低偶极磁场、低光度、暂变的X射线源(X射线流量仍然在缓慢下降,没有得到相对稳定的值).在表5中,名称依次为: 源名、自转周期(P)、自转周期导数特征年龄真实年龄估计、可能成协的物体(如SNR、大质量分子云、大质量星团、氢II区、大、小麦哲伦云、银河系中心)、实际年龄确定方法、软X射线光度、旋转能损率.动力学年龄是指磁星从最初位置(如星团)移动到现在所处位置的时间,通过测量SNR的年龄、自行年龄或特征年龄来估算.有关与超新星遗迹的距离和年龄信息请参考SNR网站http://www.physics.umanitoba.ca/snr/SNRcat/.被观测到的磁星软X射线光度是指经过引力红移的、各向同性的软X射线光度,由下式进行估算:

表4 当Bp(0)=2.0×1015 G时Bp,dBp/dt,Lp,Bt,dBt/dt,Lt和LB的部分值(假定一个中等质量的磁星M=1.45M⊙,R=11.77 km,Rc=0.98 km,对应着I=1.47I45和 µ=1.676 ; 表格上和下半部分分别对应着σ=8.75×1024s−1和 σ=2.52×1024s−1)Table 4.Partial values of Bp,dBp/dt,Lp,Bt,dBt/dt,Lt and LB when Bp(0)=2.0×1015 G.Here we assume a mediummass magnetar M=1.45M⊙,R=11.77 km,Rc=0.97 km,corresponding to I=1.47I45 and µ=1.676 ,respectively.The top and bottom parts correspond to σ=8.75×1024s−1 and σ=2.52×1024s−1 ,respectively.

其中D是源到地球的距离是指观测到的软X射线谱流,即在10–19J)范围内没有被星际介质吸收软X射线谱流).中子星旋转能量损失率被定义为

其中W是角速度是W的导数,n=W/2π是自转频率.本文采用TMA参数组,取中等质量的中子星M=1.45M⊙,对应I=1.45(2)×1045g·cm2.

图5 磁星磁场欧姆衰变的数值模拟 (a) 在x=1处极向磁场Bp随时间t的变化; (b) 在x=1处极向磁场Bt随时间t的变化;(c) 在x=1处极向磁场衰减率dBp/dt,随时间t的变化; (d) 在x=1处环向磁场衰减率−dBt/dt,随时间t的变化; (e) 极化磁场的能量衰减率Lp随时间t的变化; (e) 环向磁场的能量衰减率Lt随时间t的变化; 在(a)(f)图中红色和蓝颜色的线分别表示σ=2.52×1024s−1和σ=8.75×1024s−1Fig.5.Numerical fitting of Ohmic decay for magnetars: (a) The poloidal magnetic field,Bp,as a function of t at x=1; (b) the toroidal magnetic field,Bt,as a function of t when at x=1; (c) the poloidal magnetic field decay rate,dBp/dt,as a function of t when at x=1; (d) the toroidal field decay rate,dBt/dt,as a function of t when at x=1; (e) the poloidal field energy decay rate,Lp,as a function of t; (f) the toroidal filed energy decay rate,Lt,as a function of t.The red and blue lines in (a)−(f) indicate σ=2.52×1024s−1 and σ=8.75×1024s−1 ,respectively.

在图6中,Radio PSR 表示160颗普通的射电脉冲星[89]; HB表示15颗强磁场脉冲星; XINS表示7颗X射线的孤立中子星[90,91]; CCO表示8颗中心致密天体[92]; 实线是利用文献[89]给出的公式进行拟合得到的图线.点划线是利用文献[88]给出的经验公式

表5 具有软X射线辐射的22颗磁星的到达时间及其辐射特性Table 5.The persistent timing,ages and emission characteristics for 22 magnetars with observed soft X-ray flux.

同样地,对于10颗满足情况2)的磁星,得到拟合公式:

卡方/自由度(χ2/df)的比值反映拟合的准确度,理论上χ2/df的值越接近1越好.对于(32)和(33)式,χ2/df的值分别为8.51/7 和7.98/9,这说明拟合的表达式和数据非常吻合.随着磁星样本的增加,我们会修正拟合的表达式.

图6 在各向同性加热模型下磁星及相关致密天体 -的关系图Fig.6.The -plot for our magnetars and selected objects in isotropic heating models.

图7 在各向同性加热模型下拟合得到的磁星的旋转能损率与软X射线光度的关系Fig.7.Fitting relationship between the soft X-ray luminosity and rotational energy loss rate of magnetars in the isotropic heating model.

4.2 磁星X射线辐射和环向磁场的欧姆衰变

本节只考虑12颗旋转能损率远小于软X射线光度的磁星.假定内部环向磁场的欧姆衰变可以提供磁星各向同性的软X射线辐射.为了估计这12颗磁星的Bp(0)值,我们引用Viganò等[46]的工作:

1)通过引入最先进的动力学系数,并且考虑到霍尔漂移项的重要影响,Viganò等[46]给出中子星磁-热演化二维模拟最新结果,并与包括17颗磁星在内的40个源的观测进行比较.结果发现,仅靠改变初始磁场、质量和包层成分,在该理论模型下,磁星、高磁场射电脉冲星、孤立暗中子星等观测多样性可以得到很好地解释.

2)通过比较具有磁场Bp(0)=1015G的铁壳层中子星的冷却曲线,推断8颗磁星(CXOU J171405,SGR 1900+14,1E 1048.1–5937,SGR 0526–66,CXOU J010043,1RXS J170849,1E 1841–045,SGR 1806–20)能自诞生时就具有量级为几个1015G的磁场,磁场能提供硬X射线辐射,总的光度~1036erg·s–1[46].

3)通过比较具有磁场Bp~1—5×1014G 的中子星群的冷却曲线,得到启示: 9颗磁星(1E 1547.0–5408,SGR 1627–41,SGR 0501+4516,XTE J1810–197,CXOU J164710,1E 2259+586,4U 0142+61,Swift 1822.3–1606,SGR 0418+5729)可能具有初始磁场Bp(0)=3×1014G,这样的磁场能够解释观测的到达时间特性和持续的软X射线光度.

结合文献[25]和当前的磁星观测(软X射线光度、表面热温度和实际年龄等),我们给出12颗磁星(>Lrot)的初始磁场和相关参数,列举在表6中.在中子星内部可能存在着超流涡丝爬行、放射性同位素化学元素衰变[94]、最小冷却等各向同性加热机制,但是这些加热机制热能产生率与磁星高值的X射线光度相比忽略不计[33].我们假定各向同性的软X-射线辐射来自星体表面,给定磁星初始偶极磁场的估计值和实际年龄,计算出在两个典型电导率情况下磁场能释放率,计算结果列于表6中.

表6 12颗旋转能损率远小于软X射线光度的磁星的辐射特性及磁场能衰变率Table 6.The X-ray emission characteristics and magnetic field energy decay rates of 12 magnetars with rotational energy loss rates less than their soft X-ray luminosities.

从表6可以看到,磁星内部磁场由于欧姆衰变,磁场能释放率普遍地高于其软X射线光度1—2个数量级,主导磁能释放率的环向磁场衰变足够提供磁星宁静状态下持续的X射线辐射,表6中LB的计算结果支持了本文的理论模型.客观上讲,我们对LB的理论计算与磁星内部磁场欧姆衰变实际情形还存在一定的差异.这是由于对磁星初始磁场、实际年龄和内壳层厚度的估算还存在一定的不确定性,另外,由于星际介质对X射线的吸收以及磁星距离估计的不确定性,由磁星网站所给出的的值(表5中第8列)也存在一定的不确定性.随着对磁星观测手段、观测设备、观测方法的改进,以及理论方面的深入研究,本文模型也会得到进一步改进,理论结果将更好地符合磁星的高能观测.

由于欧姆衰变过程释放的热能,绝大部分由热中微子带走,还有一部分被星体物质吸收维持热平衡,因此,只有很少部分的湮灭的磁场能转换为热X射线光子.我们定义磁星的软X射线的转换系数为

对于每颗磁星来说,其软X射线的转换系数的大小存在一定的差异,我们利用(34)式算出在两种电导率下软X射线的转换系数的值,分别列于表6中第8和第10列中,发现12颗磁星的η分布范围为10–2—10–1.从理论上讲,一旦确定磁星的η和LB的值,通过下式来估计表面有效温度TS,

这里σS为Stefan-Boltzmann常数.由于引力红移,观测的磁星表面温度小于表面有效温度TS,两者关系:=TS(1−rg/R)1/2,其中rg=GM/c2为史瓦西半径.

5 总结与讨论

磁星的活动性和辐射特征归因于内部可能存在的超强磁场.在WGW19的工作基础上[33],我们计算了超强磁场下的中子星壳层电导率、由于欧姆衰变内部磁场衰变率和磁场能释放率.我们重新研究了22颗具有软X射线辐射的磁星的LX-Lrot关系,得到了新的拟合公式.我们发现,对于LXLrot的磁星,内部环向磁场的欧姆衰变可以提供其高光度的软X射线辐射并维持较高的表面热温度.

关于磁星初始偶极磁场,我们主要参考了文献[43],因为文献[43]首次提出了中子星磁热演化模型,这个模型可以成功地解释年轻的中子星(包括磁星以及高磁场脉冲星在内)X射线的辐射机制和冷却机制.

我们还讨论了磁星软X射线其他可能的各向异性的起源机制,例如磁斑点、热塑性流波加热等.尽管这些加热机制与欧姆衰变不同,但是这些加热机制都要求磁星内部必须存在超强的环向磁场,各向异性的加热机制要求磁星内部存在更高的磁多极场(如磁八极场),并且涉及到复杂的霍尔漂移,这些将成为我们未来关注的焦点.随着对磁星观测手段、观测设备、观测方法的改进,以及理论方面的深入研究,本文模型也会得到进一步的改进,理论结果将更好地符合磁星的软X射线辐射和表面热温度的观测.