平谱射电类星体B3 0307+380的15 GHz射电光变分析*

杨 星，易庭丰，毛李胜

(云南师范大学物理与电子信息学院，云南 昆明 650500)

耀变体(Blazar)是一类喷流与视线方向近似平行的活动星系核，其观测特征包括宽波段电磁辐射、高光度、高偏振、视超光速运动等[1]。耀变体包括平谱射电类星体( Flat Spectrum Radio Quasars, FSRQs)和蝎虎型天体(BL Lacertae Objects)两个子类， 两者的差别主要体现在光谱中有无宽发射线[1]。普遍认为， 耀变体在几乎整个电磁波段表现出各种时标(从几分钟到几年量级不等)的光变，少数耀变体表现出周期性光变现象[2-5]。光变研究一直是探索耀变体内部辐射过程和物理机制的最有效手段之一。由于耀变体喷流方向与观测视线之间的夹角小(θ< 10°)，耀变体的辐射因多普勒效应而增强[1,6]。多普勒因子(δ)可以通过喷流中物质流速度(β)以及视角(θ)两个本征参数定义。然而这两个量都是不可直接观测的量，因此，有必要用新的其他方法估算多普勒因子。在各种估算多普勒因子的方法中，通过射电光变估算多普勒因子(δV)的方法相对比较简单，因而被广泛应用[7-10]。

B3 0307 + 380(RA=47.707 8°，Dec=38.248 3°)是一个被探测到伽马射线辐射的平谱射电类星体[11]，红移z=0.816。本文利用欧文斯谷射电天文台[12]的15 GHz长时间监测数据(观测时长约11.4年，观测时间从2008年1月8日到2019年5月25日)，研究B3 0307 + 380的射电光变，估算射电亮温度和光变多普勒因子，并考察是否存在周期性光变。采用的宇宙学参数为H0=73 km·s-1·Mpc-1，Ωmatter=0.27，Ωvacuum=0.73。

1 数 据

从欧文斯谷射电天文台官网(http://www.astro.caltech.edu/ovroblazars/)获取了B3 0307 + 380的15 GHz长时间测光数据。 去除少量的坏点后， 采用的测光数据包括664个点， 约化儒略日MJD范围约为54 473～58 628。 图1给出了光变曲线，横坐标为MJD，纵坐标为以Jy为单位的15 GHz流量及误差(大部分误差很小，无法在图中分辨)。

图1 B3 0307 + 380的15 GHz光变曲线

Fig.1 The 15GHz light curve of B3 0307 + 380

2 多普勒因子估算

基于光变时标限定的尺度，文[13]给出了耀变体亮温度的计算公式：

(1)

其中，Tb为亮温度；ΔF为以Jy为单位的流量变化值(光变幅度)；tob为以天为单位的光变时标；λ为以cm为单位的观测波长；D为以Mpc为单位的光度距离。

图1显示B3 0307 + 380的射电光变曲线中表现出多个爆发过程。从中挑选12个包含上升和下降阶段的完整爆发过程， 采用文[14]提供的双指数函数公式拟合这12个爆发过程。具体公式如下：

(2)

其中，Fc为基底流量；t0为峰对应的时间；tr和td分别为指数上升时标和下降时标；F0为爆发的幅度。

(3)

图3给出了多普勒因子δV的分布情况。δV的值分布在(9.12 ± 0.58)～(35.38 ± 3.34)，平均值为18.65 ± 7.11 (sd)，中值约为18.14。

图2 12个爆发过程的曲线拟合

Fig.2 Curve fitting of twelve flares

表1 爆发过程的拟合结果

3 周期分析

历史文献中不断有关于耀变体周期性光变的报道[18-19]。耀变体的周期性光变研究， 对理解其内部辐射的物理机制和辐射区的几何性质有重要意义。目前，比较常用的周期分析方法有Lomb-Scargle周期图方法[20-21]和加权小波分析方法[22]等。 利用这两种方法探索了B3 0307 + 380的15 GHz光变曲线是否存在周期光变，结果列于图4。图4(a)为Lomb-Scargle周期图方法的分析结果，图中按照文[23]提供的基于红噪声的蒙特卡洛方法给出了95%、 99%、99.7%置信度曲线。 结果可看出，存在一个约244 ± 11天的 >99%置信度的周期，同时也存在一个约465 ± 40天的 >99%置信度的周期，后者大约为前者的两倍。图4(b)为加权小波分析结果， 图中大致可以看出，存在一个大约250天的周期，并且这个周期基本持续，特别是在观测时间的中段。

图3 光变多普勒因子δV的分布

Fig.3 The distribution of variability Doppler factorδV

4 总结与讨论

耀变体光变曲线的主要特征是在多波段呈现大幅快速光变且存在高偏振。通常认为，相对论性喷流的聚束效应是导致这些极端观测特征出现的原因。 欧文斯谷射电天文台的长时间监测数据显示，平谱射电类星体B3 0307 + 380的15 GHz射电光变曲线中表现出多个爆发过程。利用双指数函数拟合了12个爆发过程，得出光变多普勒因子分布在(9.12 ± 0.58)～(35.38 ± 3.34)范围内，平均值为18.65 ± 7.11 (sd)。 这表明B3 0307 + 380的射电辐射存在显著的多普勒增亮效应，支持被广泛接受的相对论性喷流模型。

图4 (a) Lomb-Scargle周期分析结果；(b) 加权小波周期分析结果

Fig.4 (a) The results of LSP; (b) WWZ period analysis

耀变体的周期光变可能与辐射区域在喷流中沿着螺旋形轨迹作相对论性运动有关[24]。具体来说，由于螺旋形运动导致观测视角作周期性变化，多普勒因子也发生相应的变化，从而导致多普勒增亮的辐射流量受此效应的调制。图5给出了多普勒因子在12个爆发过程之间的变化情况。

图5 光变多普勒因子δV的变化。 每个点对应的约化儒略日值取拟合参数t0及其误差，δV取上升和下降阶段对应多普勒因子的平均值及其误差

Fig.5 The variability ofδV. The MJD of each point adopts the fitting parametert0and its error, andδVadopts the average value of Doppler factors during the rise and decay processes and its error

从图5可以发现，多普勒因子随时间呈现一定的振荡。 我们猜测B3 0307 + 380的射电周期光变可能与螺旋喷流模型有关。需要指出的是，B3 0307 + 380约244天的射电周期置信度仅为大于2σ水平，还有待进一步确认。另外，关于周期的解释，螺旋喷流模型仅是可能的原因，还有很多其他的机制值得参考[18,25]。

致谢：感谢欧文斯谷射电天文台的数据支持。