roAp 星观测和理论研究综述

2021-10-15 11:55史芳菲张华伟付建宁
天文学进展 2021年3期
关键词:对流脉动振幅

史芳菲,张华伟,付建宁

(1.北京大学 物理学院 天文系,北京 100871;2.北京大学 科维理天文与天体物理研究所,北京 100871;3.北京师范大学 天文系,北京 100875)

1 引言

一般认为,恒星的元素丰度主要取决于其所诞生分子云的化学组成,并在恒星演化的过程中发生变化。在恒星形成与演化的过程中各种元素的丰度与太阳相应的元素丰度相比都保持着相近的比例[1]。但是在某些恒星中,由于一些物理过程使某些元素的丰度增高或降低,从而造成了化学丰度的异常。在早B 到早F 型星中,有大约10% 的恒星自转缓慢,在其对流层中,一些元素(主要是稀土元素)在很强的辐射压作用下被带到恒星表面,造成这些元素丰度增丰的现象。这类星被称作化学特殊星——CP (chemically peculiar)星。

CP 星主要有4 个亚群:(1)Am 星,Ca 和Sc 含量低,Fe 族元素和重元素含量高;(2)Ap 星,具有增丰的Si,Cr,Sr,Eu 或稀土元素,因普遍具有强磁场,也称为磁Ap 星;(3)HgMn 星,具有高的Hg,Mn 等其他重金属元素丰度;(4)具有弱He 线的He-weak 星[2]。其中,Am 星与Ap 星有相似的光谱特征,两者的主要区别在于Ap 星普遍存在强磁场。然而,在低分辨率光谱中无法观测到磁效应引起的光谱特征,很难识别这种差异。因此,在一些统计研究中,很容易因混淆而引入误差。

强磁场是Ap 星最重要的属性,而且普遍存在[3]。典型的有效磁场强度范围为(1~2)×10−4T,最大可达2×10−3T。与其他没有磁场特征的CP 星相比,Ap 星的自转会由于强磁场导致的制动而变得缓慢[4],有些Ap 星的自转周期甚至长达几个世纪[5]。强磁场还可以抑制对流,从而为许多具有吸收特征的离子提供一个稳定的环境。在此环境中,辐射压可以使离子抵抗重力作用而上升到表层。然而,其他元素,尤其是He,会沉入到底层。这样一来,由于磁场存在下的元素分层效应,一些化学元素在表面上分布不均匀[6],造成某些化学元素,如La,Pr,Si,Cr,Sr,Eu 和Nd 等元素含量高于“正常”的A 型星。在某些情况下,这些元素丰度可能会达到太阳丰度的100 万倍。

恒星在靠近磁极处磁场更强,元素分层会表现得更加明显,使得一些元素被辐射压带到表层并在恒星大气中形成斑点。这些斑点会造成该处辐射的吸收再辐射,而辐射的波段和流量不同造成了与周围区域的亮度差异。当这颗星自转时,亮度会随着自转发生周期性的变化。这种变化很稳定,并且可以用倾斜旋转模型来解释:恒星具有偶极磁场,其轴倾斜于自转轴。当恒星旋转时,这些斑点会周期性地产生光谱和光度变化(见图1)[7]。通过光度变化也可以推测表面“斑点”的信息,例如:TYC 2488-1241-1 光变曲线只有一个类正弦函数,推测表面有可能只有一个“斑点”;而KIC 10685175 的光变曲线有多个类正弦函数,因此表面可能有多个“斑点”。

图1 TYC 2488-1241-1 和KIC 10685175 光变曲线相位图

人们对在赫罗图上位置与经典脉动不稳定带重叠的Ap 星进行观测,发现一些Ap 星表现出了脉动,它们被称为roAp 星(快速振荡的Ap 星)。尽管Ap 星的名字中有“A”,但它们实际上覆盖了从晚B 型到早F 型的光谱类型(15 000>Teff>7 500 K),恒星质量约为2M⊙,表面具有几百到几千高斯大小的磁场[9]。其自转周期很长,从几天到几十年不等;脉动周期很短,约为4~15 min;脉动模式为非径向高阶p 模式。它们在赫罗图上位于经典脉动不稳定带的下部与主序正上方区域的交叠处。在这个区域,可以发现另外两类多模式脉动变星:g 模式脉动的γDor 型变星,脉动周期大约是0.25~3 d[10];p 模式脉动的δScuti 型变星,脉动周期大约是几小时[11]。相较于同样具有p 模式脉动的δScuti型变星,roAp 星的脉动周期更短,大约4~15 min,振幅非常低(ΔB<0.034 mag)。结合倾斜转动模型,可以很好地解释roAp 星在观测上表现出的脉动频率的自转分裂以及自转造成的脉动振幅和相位的调制现象。有的roAp 星还有不止一种模式的脉动,甚至还同时具有δScuti和γDor 模式的脉动。更详细的观测性质将在第4 章介绍。

通过对脉动的研究,利用星震学可以进一步研究roAp 星的自转倾角、磁偏角、磁场大小、半径、质量和年龄等性质,能够更精确地确定恒星年龄和演化阶段;结合有效温度、光度和金属丰度等信息构建星震学的理论模型并计算得到本征脉动模式;通过与观测得到的多重模式脉动信息进行对比,获得恒星内部结构和演化的信息,帮助我们理解恒星内部的物理过程和产生脉动的物理机制。

本文第2,3 章将介绍roAp 星的观测手段和观测性质;第4 章介绍roAp 星的脉动激发理论;第5 章介绍目前研究存在的问题;最后是总结与展望。

2 roAp 星的观测手段

自从Kurtz[12]发现第一颗roAp 星以来,截至2020 年12 月,通过地面和空间观测共发现了77 颗roAp 星(见表2)。脉动快、振幅小,导致了对roAp 星观测的困难。很长一段时间以来,人们通过地面望远镜的时序测光来研究roAp 星的脉动、几何结构和星震学。例如,在早期阶段,Kurtz 等人[12]均是通过这种方式发现了大多数roAp 星。随着Kepler 和TESS 等太空望远镜的使用,人们可以探测到振幅远低于地面观测探测极限的脉动。

Kepler 望远镜观测有两种模式:长段(long cadence,LC)观测模式,积分曝光时间为29.43 min;短段(short cadence,SC)观测模式,积分曝光时间为58.85 s。在其4 年的正常任务中,Kepler 望远镜对绝大多数恒星都采用LC 模式观测,而在SC 模式下观测的恒星数量很少。利用Kepler SC 模式观测的数据,人们已经发现了5 颗roAp 星。Murphy 等人[13]指出,虽然LC 模式会有超Nyquist 频率①Nyquist 频率:fNyq=1/2ΔTobs,ΔTobs为总观测时间。的问题,但是通过对Kepler 数据的时间项进行日心校正可以解决这一问题。Hey 等人[14]利用该方法发现了6 颗roAp 候选星。然而,史芳菲等人[15]在对同时被Kepler LC 模式和TESS观测的roAp 星KIC 10685175 进行研究后发现,虽然Hey 等人[14]通过Kepler LC 模式的数据已证认了这颗星为roAp 星,但并没有找到真正的脉动频率,而是把脉动频率与Nyquist 频率的叠加频率当成了脉动频率,这说明LC 模式中仍然存在超Nyquist 频率的问题。

TESS 的观测策略[16]是将天球的每个半球分为13 个区域,每个区域观测27 d。除了30 min 的积分曝光模式外,在每个区域选择20 000 颗恒星进行2 min积分曝光模式观测。相较于Kepler LC 模式,虽然积分时间减少会造成信噪比降低,但是对于研究roAp 星十几分钟周期的脉动显然更加合适。TESS 打开了一个发现roAp 星的新窗口。到目前为止TESS发现了很多多模式脉动roAp 星。通过对大多数已知的roAp 星的重新观测,TESS 也将发现roAp 星更多的脉动模式。同样地,TESS 在已知的或新发现的roAp 星上发现的多重自转分裂的频率,将对这些恒星的自转倾角和磁倾角提供更多约束。

在研究roAp 星时,地面和空间测光观测主要有两个不同之处:(1)地面可探测的脉动振幅一般在5×10−4~10×10−3mag 之间,而Kepler 望远镜可以观测到几微星等。(2)因为对于roAp 星来说,B 波段的脉动振幅最大[67],地面观测通常在B 波段进行。而Kepler 望远镜和TESS 的观测是在更宽的波段(Kepler 望远镜,约400~900 nm;TESS,约600~1 100 nm)中进行的(在这个波段,脉动振幅可能是B 波段的[22])。

除了时序测光,对roAp 星也可以通过时序光谱来观测研究。实际上,在探测roAp 星脉动方面,光谱探测比测光更为灵敏,因为某些谱线(大多数来自稀土元素)的脉动视向速度振幅可以达到每秒几千米。然而,光谱观测很难在保证光谱的高时间分辨率和高信噪比的同时兼顾效率。因此,除了光谱之外,时序测光观测仍是必要且有效的。

3 roAp 星的观测特征

通过研究roAp 星的观测特征发现,大部分roAp 星都表现出共同性质,仅个别roAp星表现出的性质较为特殊。总结roAp 星的共同性质,可以帮助我们了解roAp 星相较于其他变星的特性;而一些roAp 星之间的差异也为我们提供了研究特殊物理过程的机会。

3.1 一般观测性质

在大多数roAp 星的频谱图中(如KIC 10685175,见图2),不止一个,而是一组脉动频率。实际上,这一组频率在误差范围内是等间隔的,且其间隔等于自转频率。根据脉动模式的不同,脉动频率有时会分裂成三重频率(三分裂)、五重频率(五分裂)甚至更高重的分裂。如图2 就是一个被自转频率“五分裂”之后的脉动频率的例子。

图2 KIC 10685175 频谱图

Kurtz[9]在倾斜自转模型的基础上提出倾斜脉动模型:在这个模型中,磁轴与自转轴有一个夹角,脉动轴与磁轴之间的夹角可以忽略,一般看作重合,且模型同时考虑了磁场和自转,并且假定磁场对于频率分裂的影响要远大于自转。同时,这一模型也可以解释左右分裂的频率振幅不相等,并给出振幅比值A与磁偏角β、倾角i之间的关系[9]:

式中,A 的上标和下标分别代表星震学中的球谐度ℓ和方位角阶数m,具体在频谱图中的示例如图2 所示。

通过倾斜转动模型也可以测量内部磁场。以最简单的偶极磁场(ℓ=1)为例,在微扰近似下,磁场造成微扰的脉动频率为:

如果对一颗roAp 星观测到了不止一组频率,则可以通过星震学分析进行脉动模式的证认,得到恒星相关参数。例如通过自转分裂的频率间隔大小可以估算自转速度,自转分裂模式的振幅相对变化可以估算旋转轴倾角,大频率间隔(相同球谐度ℓ、相邻径向节点数n、m=0 的本征模的频率差)可以用来确定恒星平均密度等。

以HR3831 为例,如果把观测数据每50 个脉动周期分成一小段,对每一段计算脉动的振幅和相位,会发现脉动振幅和相位随着时间发生变化。把这种变化按照自转相位叠加,变化的规律性会更加明显,且周期刚好为自转周期。这种变化叫做自转对脉动的调制现象(见图3),最早被Kurtz 等人[30]发现。这也可以用倾斜转动模型解释。图4 给出了不同自转相位处的脉动示意图。在0 相位处,磁轴北极刚好处在视线方向上,此时由于脉动轴与磁轴几乎重合,越靠近脉动两极脉动振幅越大,因此此时脉动振幅最大;在0.25 相位处,磁轴翻转,脉动反相,振幅沿赤道对称,此时脉动振幅为0,相位π跳变;在0.5 相位处,虽然磁轴南极处在视线方向上,但是与视线方向夹角很大,观测到的脉动主要来自赤道,因此此时脉动振幅虽然达到极大值,但比0 相位处的振幅要小。

图3 HR 3831 脉动振幅、相位变化图[30]

图4 不同自转相位处脉动示意图[68]

3.2 特殊表现的roAp

有一些roAp 星的某些参数偏离平均值比较多,如有效温度、自转周期、脉动周期等。这些参数与物理性质相关,因此值得深入研究。roAp 星在经典脉动不稳定带上的分布使我们关注它的有效温度。目前发现最热的roAp 星是HD 43226[22],其有效温度为8 293 K,这颗星处在经典脉动不稳定带的红边缘与蓝边缘之间。最冷的roAp 星是HD 258048[61],其有效温度为6 600 K,处在红边缘之外。与此相近的还有HD 216641 (6 640 K)[22]和HD 154708 (6 800 K)[48]。自转周期与磁场大小相关,roAp 星中自转周期长于1 000 d 的有HD 166473 (10.5 a)[69]和HD 201601 (4 190 d,大于11 a)[56]。脉动周期与演化阶段有关,脉动周期最长的是HD 116114 (21 min)[39]和HD 177765 (23.6 min)。脉动振幅最大的是2MASS J19400781-4420093 (0.034 mag)[61]。

HD 42659[25]是已知的唯一一颗处于双星系统中的roAp 星。双星系统对于roAp 星的脉动、演化、磁场等的影响需要进一步研究。KIC 8677585[64]是唯一一颗同时具有δScuti 脉动、γDor 脉动和roAp 脉动的变星(见图5)。什么机制能够同时激发两种不同类型的脉动,这种变星是否常见也值得探讨。KIC 10195926[65]是一颗已经演化过主序的roAp 星,而通常roAp 星都在主序阶段。通过对HR 3831[30,70]的长时间观测,还发现一些脉动模式的寿命很短,HD 12932[21]也出现了这种模式变化,曾经观测到的脉动频率在多年之后的观测中消失,反而观测到新的脉动模式。然而这些只是目前发现的一些个例,还没有办法进行更深入和系统的研究。

图5 KIC 8677585 频谱图[71]

一些roAp 星只有一个脉动频率,如HD 161459[49]和HD 196470[55],但是其他一些roAp 星则有多个脉动模式,如HD 63087[18,22](见图6)有6 个模式,其中包含3 个单频脉动模式、2 个三分裂和1 个五分裂。也有人指出,只有一个脉动模式的roAp 星仍可能存在未被探测到的脉动模式[72]。随着探测能力的增强,会发现更多的多模式脉动roAp 星。如HD 134214,在Kreidl 与Kurtz[73]及Kurtz 等人[74]两次测光观测的分析中只发现了一个脉动频率;但在Gruberbauer 等人[75]的文章中发现,至少存在两个m值不同的脉动模式。在赫罗图上,Ap 星与经典脉动不稳定带重叠的部分,除了roAp 星外还有一部分Ap 星至今没有探测到脉动。同理,有些低振幅的脉动无法被探测到,也有可能是受观测的限制。时序光谱观测可以探测到脉动振幅更低的roAp 星。因此,人们利用这种方法发现了一些测光观测没有发现的roAp 星,如HD 115226[76],HD 116114[39]等在内共16 颗roAp 星。

图6 HD 63087 频谱图

一部分roAp 星的磁场会造成脉动模式的畸变。原本脉动振幅沿恒星表面的分布是按赤道轴对称的,但是由于磁场的影响,脉动振幅不再这样分布,而是产生一些偏差。这就造成观测上的一些特殊表现,如在自转0.25 相位左右,脉动振幅不再为0,脉动相位的跳变不再是π。很多roAp 星都有这样脉动畸变的现象,如KIC 10685175[15](见图7),HD 99563[35]等。Saio[77]在理论模型中引入了磁场的影响,可以用来拟合得到最佳的自转倾角和磁偏角。假定一颗roAp 星没有受到磁场畸变的影响,也可以通过自转分裂的振幅比计算自转倾角和磁偏角。如果这两种方法得到的结果偏差很大,则表明该roAp 星的脉动模式的确受到了磁场的影响。例如KIC 10685175 (见图8),可以看到两种情况计算的自转倾角和磁偏角不同,且不考虑磁场的情况计算的结果对观测数据的拟合效果不如考虑磁场的模型拟合得好。

图7 a) KIC 10685175 脉动振幅随自转相位的变化;b) 脉动相位随自转相位的变化[15]

图8 KIC 10685175 有无磁场计算的自转倾角和磁偏结果对比[15]

4 roAp 星的脉动激发机制

关于roAp 星的脉动激发模式存在很激烈的争论,是一个至今没有完全解决的问题。由于roAp 星在赫罗图上的位置与δScuti 的分布有所重合,人们首先想到的是δScuti 的脉动激发模式——在He II 电离区激发的κ机制。然而,由于在roAp 星内部大部分的He 没有得到足够的辐射压力支撑,而在重力的作用下沉入更深的内部,残存的He II 电离区能否激发κ机制令人质疑。但Cox 等人[79]在低He 丰度、没有对流区域的δScuti 的线性径向脉动模型中验证金属特殊性和脉动可以同时出现。这种情况可以发生在靠近经典脉动不稳定带的红边缘附近的200~500 K 范围内,范围的大小与He 丰度有关。Dolez 等人[80]在自己的模型中引入了星风的影响。Dolez 和Gough[81]通过计算发现,在磁极附近的元素分层效果最明显,这一区域的He 应该最先耗尽,然而观测显示在靠近磁极处的脉动振幅最大。于是Dolez 等人[80]在原本的原子扩散的模型中加入了星风的影响,使得He 得以在磁极附近积累,从而激发脉动。除了He II 电离区,Matthews[82]认为在极区Si 丰度足够的Si IV 电离区也可以激发κ机制,并且提出Si 丰度可以被用来区分roAp 星和没有脉动的Ap 星。

目前更为广泛接受的模型是H I 区激发的脉动,因为这些模型可以更好地拟合观测到的高阶脉动模式。Gautschy 等人[83]认为roAp 星可能有色球层,并且其特殊的温度-光深关系导致在靠近表面的低光深处温度反而升高,这使得一些高阶p 模式的振动可以稳定存在。Balmforth 等人[84]认为roAp 星中存在的强磁场抑制了磁极周围部分区域的对流,而赤道附近对流依然存在,因此建立了仅磁赤道区域保留对流的“斑点”模型。模型的结果显示,κ机制可以在极区的H I 区被激发;而在赤道附近,由于对流和湍流的影响,使得κ机制没有办法发挥作用。这一情况与太阳振动类似。Saio[77]的非绝热近似的恒星脉动模型中,恒星包层内的对流被完全抑制,磁场对振荡的直接影响被考虑在内。并且他发现,所有磁场造成畸变的高阶脉动模式都可以在磁场的作用下稳定存在,因此畸变的偶极脉动模式和四极模式在roAp 星中更容易被激发。Theado 等人[72]考虑了包层对流有没有抑制,以及不同的金属丰度和金属丰度剖面的不同模型。不同模型给出的理论脉动不稳定带的红边缘与蓝边缘相似,但比标准模型给出的红边缘更热,这实际上无法满足目前的观测结果。

除了κ机制,也有人从别的角度解释脉动。1983 年,Shibahashi[85]提出了另一种可能性:磁过稳定对流。在这种情况下,脉动过稳定是由磁力线(“粘”在等离子体上,像橡皮筋一样)的回复力所引起,它在大气的超绝热区阻止对流运动。每个质量单元与其周围环境之间的热交换使振动逐渐增长。这种振动的周期要比径向基模的周期短得多,因此这种振动是相对于磁轴对称的非常高的泛音。

要检验和限制roAp 星的理论模型,需要通过观测搜寻更多的roAp 星,在赫罗图上更为准确地划分roAp 星的分布区域和边界,从而更全面深入地认识roAp 星的脉动性质。

5 尚待解决的问题

关于roAp 星至今还有很多问题没有解决。例如:Cunha 等人[86]计算了roAp 星脉动频率的上限值——截止频率,高于截止频率的震波会直接穿透恒星表面而消散掉,不能稳定存在(见图9)。但是有观测发现,一些roAp 星存在着高于截止频率的脉动模式,比如HD 24355[23]和HD 42659[25],其脉动频率远远高于截止频率。因此,高于截止频率的脉动在什么样的条件下能够稳定存在仍然是一个问题。

图9 roAp 星在Teff-νL/M 平面上的分布[23]

此外,roAp 星在赫罗图上的位置,相较于Cunha[87]计算的理论脉动不稳定带,比较偏红边缘,而在蓝边缘没有roAp 星,却有一些roAp 星超出了理论不稳定带(见图10)。这意味着在红边缘外的其他区域,在特定的物理条件下,某些脉动模式也可能被激发。

图10 roAp 星在赫罗图上的分布[22]

人们关于roAp 脉动模式的选择也进行了一些思考。roAp 星脉动径向节点数n的分布很广泛:15~75。如果只有一种激发机制,很难产生这么大范围的径向节点数n。此前,人们认为磁场会抑制低阶模式的脉动,但是近年来逐渐探测到了一些具有低阶脉动的roAp星。Murphy 等人[88]的理论模型认为低阶模式的脉动可以在氦不足的情况下被激发,且与磁场大小有关。当极区磁场Bp大于0.15 T 时任何模式的脉动都可以被激发,Bp大于0.4 T 都不能被激发,而Bp在0.15~0.4 T之间时只有径向脉动和低阶的g 模式脉动可以被激发。Cunha 等人[86]通过构建仅在极区抑制对流的脉动模型,提出roAp 星脉动的激发机制有两种:低阶模式的脉动由κ机制激发,高阶模式的脉动由对流激发。理论的脉动不稳定带是基于对流被磁场完全抑制来计算的。如果引入湍流压的激发机制,脉动不稳定带也会有所扩展。

要想解决这些问题,首先需要扩大roAp 星的样本,最重要的就是对Ap 星无偏采样。然而,目前的地面观测发现的都是冷Ap 星。因此,下一阶段要重点搜寻热的Ap 星中是否存在roAp 星;其次,要保证在赫罗图上位置的准确性就要有准确的恒星参数——半径、有效温度、光度等。TESS 作为正在运行的空间望远镜,其观测数据的数量和质量都很适合做roAp 星的搜寻和星震学研究。

6 总结与展望

roAp 星是处于主序阶段、光谱型为晚B 到早F 型、具有4~15 min 周期的快速脉动、有强磁场并表现出某些金属元素丰度增丰的脉动变星。由于它们的自转轴与磁轴一般不重合,因而在磁轴与视线方向存在一定夹角的情况下,其自转会造成亮度变化。这种变化来源于靠近磁极处出现的元素增丰斑点,且这种亮度变化很稳定,可以用来准确地得出自转周期。自转也会对脉动造成影响,比如脉动频率的分裂、脉动振幅和相位的调制,等等。

对roAp 星的脉动既有时序测光也有时序光谱观测。光谱观测可以探测到更弱的脉动,但测光观测可以获得更多的星震学信息。曾经很长一段时间内,对roAp 星的测光只能通过地面观测进行。随着Kepler 和TESS 的升空,空间观测在roAp 星的探测上发挥了很大作用,大大扩充了成员星的样本。截至2020 年12 月,已经有77 颗roAp 星被发现,而随着TESS 的继续运行,有望发现更多的roAp 星。

roAp 星的脉动激发机制普遍被认为是在H I电离区激发的κ机制。但一些roAp 星的特殊性使得人们对已有的理论进行进一步的思考和补充研究:为什么roAp 星只集中在脉动不稳定带的红边缘,为什么一部分roAp 星超出了理论的红边缘?为什么有的roAp 星有多个脉动频率,而有的只有一个?为什么有的roAp 星可以存在比截止频率还高的脉动频率?个别roAp 星中低阶模式的脉动和δScuti 型脉动又是如何被激发的?这些问题促使人们更加深入地研究roAp 星脉动的激发机制,并在此基础上推测roAp 星脉动的激发机制有两种:低阶模式的脉动是由κ机制激发,高阶模式的脉动由对流激发。

为了进一步研究roAp 星的脉动,需要更多的无偏且具有准确恒星参数的样本。随着探测精度的不断提高,有望发现更多的roAp 星,从而更为准确地确定roAp 星在赫罗图上的分布区域;也可以对已知的roAp 星进行重新观测和深入研究。如果能探测到更多的脉动模式,就可以进行星震学研究,为理论研究提供更丰富的信息。

猜你喜欢
对流脉动振幅
齐口裂腹鱼集群行为对流态的响应
RBI在超期服役脉动真空灭菌器定检中的应用
十大涨跌幅、换手、振幅、资金流向
十大涨跌幅、换手、振幅、资金流向
超临界压力RP-3壁面结焦对流阻的影响
十大涨跌幅、换手、振幅、资金流向
沪市十大振幅
基于ANSYS的自然对流换热系数计算方法研究
浅谈我国当前挤奶机脉动器的发展趋势
地脉动在大震前的异常变化研究