在超强磁场中修正的相对论电子压强*

董爱军 高志福 杨晓峰 王娜 刘畅 彭秋和

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

2)(中国科学院新疆天文台,乌鲁木齐 830011)

3)(南京大学天文与空间科学学院,南京 210000)

4)(贵州省射电天文数据处理重点实验室,贵阳 550001)

当前脉冲星领域一个重要的研究热点是磁星.本文在朱翠等(Zhu C,Gao Z F,Li X D,Wang N,Yuan J P,Peng Q H 2016Mod.Phys.Lett.A 31 1650070)工作的基础上,重新研究了磁星超强磁场下(B≫Bcr,Bcr 是电子的量子临界磁场)电子朗道能级的稳定性及其对电子压强的影响.首先,对弱磁场极限下(B≪Bcr)中子星内部电子压强进行必要的回顾;然后,通过引入电子朗道能级稳定性系数gν 和Dirac-δ 函数,推导出在超强磁场下修正的相对论电子压强Pe 的表达式,给出表达式适用条件:物质密度ρ≥107 g·cm—3和Bcr ≤B <1017 G(1 G=10—4 T).超强磁场通过修正相对论电子的相空间,提高了电子数密度ne,而ne 的增加意味着Pe 的增加.利用修正的电子压强表达式,讨论了超强磁场下费米子自旋极化现象、电子磁化现象以及超强磁场对物态方程的修正.最后,本文的结果与其他类似工作进行对比,并对未来的工作进行展望.本文的研究将为磁星以及强磁化白矮星的物态方程和热演化的探索提供极有价值的参考,将为普通射电脉冲星等离子磁层数值模拟、高磁场脉冲星辐射机制等相关研究提供有用的信息.

1 引言

磁星是主要由磁场提供辐射能量的一类脉冲星,属于近些年来被关注的特殊类的致密天体[1,2],其表面磁场超过电子的临界磁场Bcr(Bcr=4.41×1013G(1 G=10—4T)),在这样强的磁场下,电子朗道能级被强烈地量子化[3,4].磁星的典型特征之一是X 射线波段和软γ射线波段有着超过爱丁顿光度的短期爆发;部分磁星会发射出中等或巨型闪耀,后者在不到0.5 s 内能释放出高达1046erg(1 erg=10—7J)的能量.软γ射线重复暴和反常X射线脉冲星被认为是磁星候选者.磁星和太阳耀发常被解释为由剧烈磁重联引起,但是磁星在宁静状态下X 射线辐射常被归因于超强磁场的衰减[3-7].

中子星内部处于β平衡中的完全简并的相对论电子气体,其分布函数f(Ee)服从费米-狄拉克统计.当温度T→ 0 时,电子化学势µe称为电子费米能.电子费米能是致密星物态方程中一个非常重要的物理量.近年来,国内外不少作者对致密星环境下的物态方程、中微子辐射和相对论电子进行关注和研究[8-14].其中,刘晶晶等[15-19]对包括修正的Urca 过程、β衰变、电子俘获以及中微子吸收等进行长期的研究,并取得了丰硕的成果.高志福等[20]通过引入超强磁场下的Dirac-δ函数,给出了超强磁场下简并的相对论电子压强Pe(若无特殊强调,本文Pe单位均为dyns·cm—2)的一个特解:

其中Ye是电子丰度,B是磁场强度,ρ是物质密度,饱和核密度ρ0=2.8×1014g·cm—3.方程(1)的适用范围为ρ≥107g·cm—3和Bcr≤B<1017G.文献[20]的主要结论有:磁场强度越大,电子压强越大,磁星的总压强总是各向异性的,与普通射电脉冲星相比,磁星可能是密度更大的中子星;若考虑磁能对物态方程的贡献,磁星的质量可能更大.

朱翠等[21]通过引入电子朗道能级稳定性系数,研究了磁星内部电子朗道能级的稳定性,推导出超强磁场下相对论的特解,并讨论了特解的适用范围.1965 年,Kubo[22]给出关于强磁化电子气体的微观状态数Npha的计算,得出在超强磁场中随B的增加而减低的结论,彭秋和等[23]对这种计算方法有所质疑.我们发现文献[22]人为地引入了一个错误的假设,即将非相对论电子回旋运动方程的解应用到相对论电子运动中,给出超强磁场中电子气体在z轴方向pz→pz+dpz的动量区间能态密度,得出∝1/B的结论,其中 ℏ 表示被约化的普朗克常数.在过去的50 年中,其他作者反复地引用文献[22],这可能导致人们错误地理解强磁化中子星的物理.当然,由于地球上缺乏中子星量级的强磁场,人们无法在实验室中直接验证与B的关系.我们认为在超强磁场下处理量子化的电子朗道能级应该引入Dirac-δ函数,如果B≥1017G,窄长的朗道柱面将会简化为一维或二维的线性链,文献[20,21]中的结论将不再适用.

由于篇幅的局限性,朱翠等[21]没有对超强磁场中电子压强Pe和磁场强度B、物质密度ρ及电子丰度Ye之间关系进行研究,因此没有给出关于Pe与B的表达式,而这样的表达式无疑对于修正中子星的物态方程至关重要.第2 节对中子星内部弱磁场近似下的电子压强进行回顾;第3 节推导出超强磁场中修正的电子压强表达式,并与其他工作进行对比;第4 节讨论超强磁场对电子相空间的修正、超强磁场下费米子极化和磁化问题;第5 节在相对论平均场模型下给出超强磁场对物态方程影响的相关估算,讨论中子星内部压强的各向异性;第6 节进行总结与展望.

2 中子星弱磁场极限下的电子压强

如果中子星内部磁场太弱,磁场对的影响可以忽略.在弱磁场极限下B≪Bcr,和Pe均由ρ和Ye共同决定,Pe可写成多方形式:Pe=KρΓ,其中K和Γ为常数,在以下两种情况下取值不同.

1)对于非相对论性电子,ρ≪107g·cm—3,

2)对于相对论性电子,ρ≥107g·cm—3,

这里,A和Z分别表示核子数和质子数.对于Z和A都给定的原子核,电子与质子的丰度相等:Ye=Yp=Z/A;对于理想的中子-质子-电子(n-p-e)气体,Ye=Yp=ne/(np+nn)≈ne/nN,其中np,nn和nN分别表示质子数密度、中子数密度和核子数密度.将电子数密度ne=NAρYe(NA为阿伏伽德罗常数)代入方程(2)和方程(3)中,得到中子星壳层中Pe随ne的关系Pe2.45×,该式由李新虎等[24]文章中的方程(14)可以很容易算出.

图1 给出了弱磁场极限下中子星内部相对论电子的Pe随ne的变化.可以看出,Pe随ne的增加而增大,点-虚线(dM/dρ＞ 0)表示星体的质量随着物质密度的增加而增大,具有这样质量的中子星是稳定的;实线(dM/dρ<0)表示星体质量随密度的增加而减小,dM/dρ<0 这样的中子星不稳定.

图1 中子星内部弱磁场极限下相对论电子压强Pe 随电子数密度ne 的变化Fig.1.Relativistic electron pressurePe with electron number densityne in the limit of weak magnetic field inside a neutron star.

3 超强磁场中相对论电子的压强

3.1 超强磁场中电子朗道能级稳定性

为了得到超强磁场中修正的电子压强表达式,需要回顾朱翠等[21]的工作.他们首次引入超强磁场中电子的朗道能级稳定系数gν,并假定gν具有幂率指数的形式:

式中ν=n+1/2+σ是电子占据的朗道能级数,n是主量子数,σ=±1/2 是电子自旋值,自旋向下(↓)的电子占据的能级数ν=0,1,2,3,···,自旋向上(↑)的电子占据的能级数ν=1,2,3,···;g0和α分别表示电子的基态能级稳定性系数和能级稳定性指数.当ν=1 时,g0=g1,即基态和第一激发态能级具有相同的稳定性.为简化起见,假定g0=g1=1.由量子力学可知,占据能级数n较大的电子具有较高的能量,容易发生激发跃迁,从而占据ν较低的能级,并且ν越大能级稳定性越低,即gν＜gν-1＜gν-2,因此α为负数.在超强磁场中,电子的费米面变成了狭长的朗道柱面.定义无量纲磁场B*B/Bcr,电子最大的朗道能级数为[21]

其中Int[·]表示对括号内的值取最大整数,σ=±1/2 是电子自旋,pz为电子在z轴方向的动量,me是电子的质量,c是光速.不难看出,νmax随磁场增强而减小.

3.2 超强磁场中电子能态密度

将电子磁矩假定在一个沿着z轴方向的强磁场B*中,电子的费米能为

其中,δ[p⊥/mec-(2B*)1/2] 为Dirac-δ函数,物理意义是:由于朗道能级强量子化,在第n个和第n+1 朗道环面之间不存在任何微观量子态.Dν表示自旋简并度,当ν=0 时,Dν=1,ν≥1,Dν=2,则方程(8)写为

由方程(5)和方程(6)可知,当pz=0 时,则有

于是方程(9)被化简为

将方程(10)代入方程(11)后,则有

引入比率qI2/I1,其中

当νmax≥6 时,比率q≈ 1,方程(11)中的求和用积分来代替,具体为

引入变量t=pzc/,于是方程(12)被简化为

式中λe=h/(mec)表示电子的康普顿波长.

3.3 超强磁场中相对论电子的压强

根据泡利不相容原理的要求,电子能态密度等于电子数密度,于是得到

其中,

求解方程(16)给出

即始终是由ne确定,而后者由B,ρ和Ye共同决定.高志福等[20]已经证明:在平行和垂直于磁场方向电子的最大动量都等于电子的费米动量,

将电子动量进行体积分,得到电子压强:

式中φ(xe)是多项式,

方程(23)适用于ρ≥107g·cm—3和Bcr

图2 不同磁场下中子星内部电子压强Pe随物质密度ρ 的变化Fig.2.Relation between electron pressurePe and matter densityρ in neutron stars with different magnetic fields.

中子星壳层以电子简并压为主导,特别是在年轻的强磁化中子星(如磁星),星体壳层中磁化率、热导率、比热等物理量振荡行为可能与费米面附近的电子在强磁场中的活动性有关,这些行为本质上体现了电子朗道能级的不稳定性,未来我们将深入地讨论上述问题.近30 年来,强磁场对电子朗道能级的影响早已引起人们广泛的广注,研究结果已成功地应用于依靠电子简并压来抵抗引力的白矮星研究之中[25-31].下面将以三篇典型性的工作[25-27]为例,对强磁场中电子气体的压强的计算方法进行简介,将前人的工作及与本文的工作进行对比.

比较发现方程(24)与方程(7)一致.超强磁场中电子数密度ne为

方程(27)和方程(23)是Pe在强磁场中两种不同的表达式.两者存在相同点和差异性,相同点包括:1)电子看成理想的气体,并且忽略了电子反常磁矩;2)电子压强都与磁场强度B*、电子数密度ne和电子费米能有关,并且ne和都是包含B*的复杂函数;3)给定ne和,Pe随B*的增加而增大;4)随着磁场的增加,由方程(27)拟合的压强随密度变化曲线会出现不规则的突起或波动,这是因为在第ν个能级电子能态由部分填充到完全填充转变,或电子由第ν个到第ν+1 个能级跃迁引起,这种现象被认为与费米面附近的电子在强磁场中行为有关[25,26],其本质上体现了朗道能级不稳定性.

两者的差异性在于:1)方程(27)是一个包含能级数ν、复杂的求和公式,并且νmax取值非常小,通常情况下0 ≤νmax≤ 3[26-28];费米面附近的电子主要集中在基态ν=0 和第一激发态ν=1,有意思的是在给定电子费米能和物质密度时,B*越强,νmax越小,并且ν值小的能级所对应的电子压强较大;在强磁场下,随着电子能量值的增加,电子态密度下降到非常小的值,使得压强梯度变得非常陡峭,增加了电子的向外简并压,从而解释了超钱德拉塞卡质量极限的强磁化白矮星;2)本文给出的方程(23)是一个没有包含能级数ν、简洁的积分表达式,但在推导过程中要求νmax取值较大,通常情况下νmax≥6,使得求和变为积分,在超强磁场下,电子态密度和电子能量值同时增加,电子向外简并压也随之增加,不需要考虑电子压强的梯度,同样支持强磁化白矮星模型.产生上述差异性的原因是:在以往的工作中没有引入电子朗道能级稳定性系数和Dirac-δ函数.

图3 是本文工作与强磁场中电子数密度和电子压强其他研究的对比.图3(a)显示强磁化白矮星中电子压强Pe随ρ的变化[28].白矮星由碳元素12C 组成,压强P=Pe+PL,这里Pe是按照方程(27)计算的电子压强,PL是晶格压强(受磁场影响忽略不计);此外实线、长虚线、短虚线和点线分别表示星体中心磁场B*=0,10,100 和200 时的拟合曲线.比较发现:当0 ≤B*≤ 10 时,磁场对白矮星的状态方程和星体结构的影响非常小;当B*≫1 时(如B*≥100),强磁场的效应则不能忽略.图3(a)中,由于PL≪Pe,因此不需要考虑PL.图3(a)还表明在相同的物质密度下,Pe随B*的增加而增大,且曲线出现多处突起,说明电子的朗道能级不稳定.

图3 本文与其他强磁场中电子数密度和电子压强研究的对比(a)强磁化白矮星中电子压强Pe随ρ 变化关系;(b)中子星壳层电子数密度ne随ρ 变化关系;(c)磁化白矮星中(最大电子费米能量EFmax=20mec2)电子压强Pe随ρ 变化关系;(d)两种不同的理论模型下白矮星中电子压强Pe随ρ 变化关系Fig.3.Study of electron number density and electron pressure in strong magnetic fields by other authors and their comparison with this work:(a)Relationship between electron pressurePe andρ in a strongly magnetized white dwarf(WD);(b)relationship between the electron number densityne andρ in the crust of a neutron star;(c)electron pressurePe as a function ofρ in a magnetized WD with maximum electron Fermi energyEFmax=20mec2;(d)electron pressurePe as a function ofρ in a magnetized WD under two different theoretical models.

与强磁化白矮星对比,图3(b)给出了中子星内壳层ne在不同磁场下随ρ的变化[29],点线、虚线和实线分别表示壳层磁场强度B*=0,103和104时的拟合曲线,后两条曲线由方程(23)拟合得到.对于B*=103或更小的场强,大量的朗道能级被电子填充,电子数密度与零磁场的结果几乎没有区别.对于较强磁场,例如B*=104,电子分布在第0 个(ν=0)朗道能级或某个质量密度范围内的几个能级,在这种情况下与零磁场情况相比,电子数密度显著增强.而本文的计算给出关系式:(B*)1/6ne,其中和ne分别为超强磁场下和零磁场下的电子数密度.这个关系式应用到中子星壳层,给出的电子数密度变化趋势与图3(b)给出的变化趋势基本一致,因为强磁场通过对晶格能的修正会改变原子核的序列,提高质子丰度和核子平衡密度,电子的化学势也得到提升,相关的工作可以参考作者在2013 年的工作[20].然而在密度较高时,在图3(b)中没有看到电子丰度的增加,反而是电子数密度变化随磁场的增加而减低,这与本文给出的结论相反.文献[29]给出有意思的解释:在强磁场下,在密度较高区域从原子核中流出的中子减少.在中子星高密度区域,如外核与内核区,由于中子简并压主导中子星的物态方程,文献[29]对电子压强不再讨论.需要强调的是,文献[29]给出B*=104这样的强磁场超出本文电子压强表达式的适用范围,给出的低磁场B*≤ 103与本文给出的低磁场B*≤ 1 相差较大.

为了和文献[27]的结果进行比较,图3(d)给出两种不同的理论模型下白矮星中电子压强Pe随物质密度ρ的变化,密度范围为ρ=1.2×109—1.17×1010g·cm—3.黑色实线表示B*=0时Pe随ρ的变化,根据Pe=2.45×10—17ne4/3=2.45×10—17(NAρYe)4/3得到Pe的变化范围为6.31×1026—1.35×1028dyns·cm—2,红色点-虚线表示在EFmax=20mec2和B*=199.5 情况下Pe随ρ的变化.为了简化计算,取νmax=1,由方程(27)得到Pe的变化范围为8.46×1027—3.27×1028dyns·cm—2;蓝色点-虚线表示在B*=199.5时Pe随ρ的变化,由方程(23)得到Pe的变化范围为1.24×1028—2.56×1029dyns·cm—2.对比DM2012 模型[27],发现在相同的磁场下电子压强随密度增加较快,所对应的最大电子费米能远大于20mec2,这是因为两种不同的理论模型给出电子数密度、电子费米能和电子压强的计算方法存在着一定的差异.

4 超强磁场下量子电动力学现象

4.1 超强磁场对电子相空间的修正

超强磁场通过修正相对论电子的相空间,会导致简并的相对论电子重新分布.根据泡利不相容原理,电子将从朗道能级的最低量子能态填充到最高能态.随着磁场强度B的增加,越来越多的电子占据更高的朗道能级,但随着电子朗道能级数的增加,电子朗道能级稳定性系数会减小.在超强磁场中电子数密度的增加意味着电子简并压强的增加.

在中子星内部,简单的中子衰变和连续的电子俘获同发生,然而超强磁场可能会使前者的速率大于后者.由于更多的中子转化为质子,质子的丰度增加,根据电中性要求,电子丰度相应增加,因此ne也会增加;另外,由核物理知识可知,质子的丰度反映核物质的非对称性,质子丰度的值与核物质的对称能、体积束缚能等参量密切相关[30-35];超强磁场可能会增加核物质的非对称性,提高质子的丰度,因此核物质的平均电子数密度也相应地增加.当然,这些理论上的可能性有待于实验验证.

Table 1.Partial c alculations ofnN,,Pe,PandMina relativisticmeanfieldmodel withtheTMA parameter set.表1 在相对论平均场TMA 参数模型下nN,,Pe,P和M 的部分计算值

Table 1.Partial c alculations ofnN,,Pe,PandMina relativisticmeanfieldmodel withtheTMA parameter set.表1 在相对论平均场TMA 参数模型下nN,,Pe,P和M 的部分计算值

4.2 超强磁场下费米子自旋极化现象

最近,不少工作[35-39]对中子星内部零温(T=0)和强磁场下费米子系统自旋极化现象进行研究.本节考虑自然单位制,费米子数密度为[36]

式中,无量纲变量(ε(B)是单粒子的能量)γf的表达式为

从方程(32)可以看出,当γf→ ∞时,对应零磁场下粒子自旋完全非极化场景,有n↑=n↓=n/2;对于γf→ 1,对应超强磁场下费米子自旋完全极化场景,有n↑=n和n↓=0.由于当γf<1 时,n↓的值变为负值,因此γf的最小值为1.对应粒子自旋完全极化场景,中子星物质所能维持的最大磁场强度即磁场强度饱和值Bs为

图4 中子星内部费米子完全极化场景下饱和磁场强度Bs 随粒子数密度n 的变化关系(a)质子/电子完全极化下Bs vs.ne/np;(b)中子完全极化下Bs vs.nB(nB 为重子数密度)Fig.4.Relationship between the saturated magnetic field strengthBs and the particle number densityn in a fully polarized neutron star fermion matter:(a)Bs vs.ne/np in a fully polarized scenario for proton/electron matter system;(b)Bs vs.nB in a fully polarized scenario for the neutron matter system(nB is the baryon number density).

4.3 超强磁场中电子磁化现象

由4.2 节可知,当中子星内部出现电子自旋极化现象时,自旋平行于和反平行于磁场方向的电子数密度不再相等,电子系统被磁化,中子星内部出现一个宏观上的感应磁矩,感应磁矩就会导致一个感应磁场,后者的大小取决于电子系统的磁化率χ或磁化强度M.根据最新研究[39],磁星内部可能包括化石起源的原始磁场和顺磁磁化产生的感应磁场[38].后者为磁星提供制动力矩,影响磁星辐射特性和内部热演化.电子磁化由两种成分构成:一是电子内禀磁矩在外磁场中的顺向分布引起的泡利顺磁部分;二是在外场中电子轨道运动量子化引起的朗道抗磁部分.高志福等[40]采用霍尔斯坦-普里马可夫变换等方法,对外磁场中的铁磁体、反铁磁及亚铁磁体的自旋波谱讨论并计算了临界磁场.由电动力学可知,磁化强度M=χB及相对磁导率µr=1/(1—µ0χ),µ0为真空磁导率,当µ0χ→ 1时,磁场变得越强,出现临界磁化现象.强磁场中相对论电子的磁化系数为[41]

对于生鲜类冷冻冷藏产品，采用全程冷链关键。从供应商的出库到电商的冷库，中间不能离开冷藏，仓库在验收入库时要检查送货车辆以及温度是否达标。水果的新鲜度和成熟度也是验收时的重点关注内容，如果拟入库商品在短期内要出库或者销售状况好、市场需求大，就适合选择成熟度比较高的商品，但是如果拟入库商品有较长的销售期，就适合选择成熟度比较低的商品，这样才能避免配送到消费者手上的时候出现腐烂变质的情况。不过成熟度偏低的产品在口感和品质上会有一定程度的欠缺，所以在采购入库时要尽量避免这一类。

利用方程(35),图5 给出不同强磁场下相对论电子的磁化率χ随ne的变化,水平的点-虚线表示临界磁化率(在高斯单位制下µ0=1).当ne一定时,χ随B的增加而增大,当B一定时,χ随ne的增加而增大.在中子星高密度区域会出现临界磁化现象,甚至磁化系数大于1.本文不考虑质子和中子磁化现象,因为核子磁化所需要的磁场更高(B≥5×1018G),在这样高的磁场下本文得到的电子压强修正表达式不再成立.在中子星内部简并电子系统的磁化将扮演重要角色.对于原生中子星,在以中微子辐射主导的冷却过程中,随着温度的降低,一方面物质分布的层次性使得中子星内部磁化不均匀,导致新的磁场诱导项产生;另一方面顺磁化的电子系统相当于减弱了磁场的扩散,从而出现等效扩散系数等于零的临界情况,在中子星磁场演化过程中增加了相变可能性.随着扩散的减弱,磁化强度会随着磁场的增加发生振荡,出现De Haasvan Alphen 不稳定性,中子星表现出磁星的活动性,如耀斑和外暴[1-5].

图5 不同磁场下中子星内部相对论电子的磁化率χ 与电子数密度ne 的变化关系Fig.5.Relation between the magnetic susceptibilityχ and number density of relativistic electronsne in neutron stars with different magnetic field strengths.

5 超强磁场对物态方程的修正

5.1 超强磁场下相对论平均场理论模型

相对论平均场理论(relativistic mean field theory,RMFT)已成为研究有限核物质性质的标准方法.TMA 参数组与相对论Brueckner-Hartree-Fock 理论预测相一致,能很好地描述中子星状态方程[42-48].将结合RMFT 框架下的TMA 参数组及本文所得到的超强磁场中电子压强的表达式,求解广义相对论流体静力学平衡方程(TOV 方程),研究超强磁场如何影响中子星的物态方程.主要讨论由核子(n,p)、轻子(e,µ)组成的传统中子星,记作npeµ物质.强子之间的相互作用是由介子(σ,ω,ρ)来传递.核子数密度nN为中子数密度nn和质子数密度np之和.强磁场环境下,总的相互作用拉格朗日密度为[46,47]

其中ψN和ψl分别为核子和轻子的Dirac 旋量场;σµν=i/2[γµ,γν];κN和κl分别为核子和轻子的反常磁矩系数,κp=1.7928µN,κn=—1.9130µN,κe=1.15965×10—3µB,κµ=1.16952×10—3µB[48,49];τ是核子同位旋,τ0是τ的第三分量;σ,ω和ρ分别表示标量介子、矢量介子、等矢量-矢量介子场,ωµν=∂µων—∂νωµ和Rµν=∂µRν—∂νRµ分别表示ω和ρ介子的场张量;U(σ)=(1/3)g2σ3+(1/4)g3σ4是σ介子非线性自相互作用项;电磁场张量Fµν=∂µAν—∂νAµ,Aµ=(0,0,Bx,0);m和ml核子分别表示核子和轻子的质量;mσ,mω和mρ分别表示σ,ω 和ρ 介子的质量;gσ,gω和gρ分别是σ,ω 和ρ 介子与核子之间的耦合常数;有关TMA 参数组核物质饱和特性和耦合常数,请参考文献[45]中表4.本文采用一种广泛应用的、从表层到中心随密度变化的中子星磁场模型[44],具体形式如下:

其中,β=0.001,γ=3,Bsurf和Bcent分别表示中子星表面和中心磁场.根据电子极化限制和磁星表面磁场观测,选取Bcent=1.0×1016G及Bsurf=1.0×1014G.在如此强度的变化磁场下反常磁矩对中子星物态方程影响较小,故而暂不考虑核子和轻子的反常磁矩.图6 表示中子星内部磁场B随密度ρ的变化关系.

图6 中子星内部磁场B 随物质密度ρ 的变化关系Fig.6.Relation of the magnetic fieldB and matter densityρ in a neutron star.

把(36)式代入Euler-Lagrange 方程,得到核子、轻子和介子在强磁场中运动方程:

每个组分的费米能与它们的费米动量之间关系为

其中m*=mN—gσσ是核子有效质量,量子数s=+1 和—1 分别表示粒子自旋向上和自旋向下.核子标量密度ρS分别为

物质系统中的中子、质子、缪子与电子的粒子数密度与费米动量的关系分别为

其中εp,εn,εe和εµ分别为质子、中子和电子和缪子的能量密度,具体形式为

考虑电磁场对能量密度的贡献,总的能量密度ε和热力学压强P分别为

利用中子星结构的TOV 方程来研究星体的宏观性质,具体形式为[50,51]

其中M(r),ε(r)和P(r)分别为以半径为r的球面内中子星的质量、能量密度和压强.以物态方程(51)作为输入量,在边界条件P(R)=0和M(R)=0情况下可以获得TOV 方程的解.由于中子星中心的重子数密度由nN=(M/mN)/(4πr3/3)来决定,得到中子星质量M与nN的变化关系.表1 中的第2—5 列表示在零磁场中的计算值;第6—9 列表示在密度依赖的中子星强磁场模型下的计算值.由表1 可以看出,在TMA 参数模型下,电子费米能、电子压强Pe和压强P都随星体中心的nN的增加而增大,在相同nN的情况下,、Pe和压强P都随着磁场的增加而增大.当重子数密度nN=0.9156 fm—3时,中子星质量达到最大值,当B≪B*时,M的最大值为1.9916M⊙(M⊙为太阳质量),当B＞B*时,M的最大值为2.034M⊙,增加的质量来源于磁场能的贡献.本文模型能够解释两倍太阳质量的中子星[52],这表明本文选择的TMA 参数组比较可靠.

5.2 超强磁场下压强的各向异性

5.1 节在球对称空间下求解了TOV 方程,忽略了压强的各向异性.在实际的中子星环境下,超强磁场会改变空间的球对称性,在平行于磁场的方向的压强P//与垂直于磁场方向的压强P⊥存在着一定的差异.在本文所选择的密度依赖的强磁场模型下(B≤ 1016G),核子的反常磁矩忽略不计.系统总的能量-动量张量Tµν 由物质部分和电磁场部分所组成,分别表示为

其中磁化张量由M=—∂εm/∂B来决定,磁场张量BμBμ-B2.方程(52)改写为

垂直于磁场方向的压强和平行于磁场方向的压强由下列方程决定:

将方程(51)代入方程(54)中,得到

为了讨论压强的各向异性,定义两个方向上的压强之差:

2013 年我们对磁星内部压强的各向异性进行详细地讨论[20],当B≪1020G 时,星体内部磁化方向与外磁场方向相反,即M<0,ΔP=P⊥—P//＞ 0,或P⊥＞P//,沿磁场方向磁张压使中子星发生形变,变成类似地球的旋转椭球星;若B≥1020G 时,会出现中子磁矩顺磁磁化,则M＞ 0,故P⊥

表2 相对论平均场模型下nN,ρ,B,ne,|MB|,ΔP和P//的部分计算值,这里选择TMA 参数组和密度依赖的中子星强磁场模型Table 2.Partial calculations ofnN,ρ,B,ne,|MB|,ΔP和P//in a relativistic mean field model.TMA parameter set and a density-dependent magnetic field model for a neutron star are selected.

6 总结与展望

在我们以往的工作基础上,通过引入电子朗道能级稳定性系数和Dirac-δ函数,本文推导出在超强磁场下修正的相对论电子压强的表达式;讨论了超强磁场对电子相空间的修正、超强磁场下费米子自旋极化现象、超强磁场中电子系统磁化现象以及超强磁场对物态方程的修正.磁星内部可能包括化石起源的原始磁场和顺磁磁化产生的超强感应磁场.后者提供制动力矩,影响磁星辐射特性和内部热演化.彭秋和等[53]认为:磁星的超强磁场可能起源于各向异性的3P2中子超流在低温下顺磁所产生的感应磁场,最大的偶极磁场强度为Bd为3×1015—4×1015G.本文所得的结论适用于文献[53]的磁星模型.方仁洪等[54]从均匀磁场中单个大质量狄拉克费米子的朗道能级出发,构造了系统的配分函数,通过配分函数可以导出包括压强在内的所有的热力学量(如矢量流、能量-动量-张量等);与本文不同的是,方仁洪等[54]是在有限温度下的狄拉克费米子系统的热力学(也可得到温度趋近零时情形);文献[54]与本文的结果大致相同,但处理问题的方法和侧重点不同,前者为强磁场,后者为一般磁场,得到热力学压强后,同样可求出狄拉克费米子的物态方程.

由于篇幅有限,只选择传统的中子星内部的npeµ物质系统为代表,讨论了超强磁场对物态方程的影响,没有讨论量子电动力学效应对弱相互作用过程和辐射机制的影响,也无法对超强磁场在中子星其他物质区域的量子电动力学效应逐个进行讨论,这些将是我们未来研究工作的方向之一.本文研究为磁星及依靠电子简并压来抵抗引力的强磁化白矮星的物态方程和热演化的研究提供了参考,可为射电脉冲星等离子磁层数值模拟[55-57]、高磁场脉冲星及磁星的辐射机制、限制中子星引力波的上限[58,59]等相关研究提供有用的信息.