S40954 65喷流的高分辨率观测研究.pdf
《S40954 65喷流的高分辨率观测研究.pdf》由会员分享,可在线阅读,更多相关《S40954 65喷流的高分辨率观测研究.pdf(12页珍藏版)》请在咨信网上搜索。
1、第 42 卷第 2 期天文学进展Vol.42,No.22023 年 6 月PROGRESS IN ASTRONOMYJun.,2023doi:10.3969/j.issn.1000-8349.2023.02.06S4 0954+65 喷喷喷流流流的的的高高高分分分辨辨辨率率率观观观测测测研研研究究究王 玕1,2,3,路如森1,4,沈志强1,4,江 悟1,4,黄 磊1,5,赵杉杉1,闫 玺1,3,程 远1,2,3(1.中国科学院 上海天文台,上海 200030;2.上海科技大学,上海 201210;3.中国科学院大学,北京100049;4.中国科学院 射电天文重点实验室,南京 210008;5.
2、中国科学院 星系宇宙学重点实验室,上海 200030)摘要:S4 0954+65是一个典型的BL Lac天体,观测表明,其射电、光学和 射线波段的辐射具有准同步耀变现象。利用美国甚长基线干涉阵列VLBA 对该源分别在2018 年6 月和2019 年11月进行多波段观测,发现其射电核的流量密度谱有明显变化。通过同步自吸收机制,得到射电核心磁场强度为(3.30.9)104T。通过分析2017年5月2021年5月期间VLBA 43GHz观测的归档数据,发现2 个明显的喷流运动成分和靠近核心的一个静止成分。对两个运动成分分析得到喷流的视向夹角分别为5.5和6.9,张角为1.8和2.4。利用VLBA 4
3、3GHz归档数据以及Fermi 卫星4 期源表(4FGL)中该源的 射线流量,发现射电核心和 射线在2018 年10 月2021年5月间变得更加活跃,且射电核心流量变化更强。关键词:BL Lac 天体;喷流结构;核心磁场;喷流成分中中中图图图分分分类类类号号号:P157.7文文文献献献标标标识识识码码码:A1引言S4 0954+65 是一个典型BL Lac 天体,具有相对论性单边喷流1。其与视线方向夹角为1.59.324。Becerra等人5通过测量MgII2800A 发射线,得到它的红移z=0.36940.0011。Fan和Cao6通过测定H的宽度,得到其中心超大质量黑洞的质量MBH3.3
4、108M(M为太阳质量)。S4 0954+65的光变十分剧烈。在射电和光学波段,时标为几小时的光变被探测到很多次,特别是处于耀变状态时713。Morozova等人12发现,在2011年3月至4月间,该源在43GHz 射电、光学和 射线波段表现出准同步耀变现象;他们还发现三个射电成分的收稿日期:2022-04-27;修回日期:2022-05-16资助项目:国家自然科学基金(11933007);中科院基础前沿科学研究计划从0 到1 原始创新项目(ZDBS-LY-SLH011);上海市基础研究特区计划(JCYJ-SHFY-2021-013)通讯作者:王玕,258天 文 学 进 展42 卷258天 文
5、 学 进 展42 卷258天 文 学 进 展42 卷出现与 射线波段耀变相关。在这期间,光学R 波段的流量也有显著变化,其中最明显的变化发生在2011 年3 月9 日,此时,流量增加了0.7mag。在2011 年2 月27 日至2011年3 月19 日期间,光学R 波段流量一直上升,偏振位置角平滑变化300。2015年2 月,该源在VHE(100 GeV)射线首次被探测到,被认为与新生成的射电成分有关联14。Marscher等人15认为耀变是由于一个新生成的射电成分在喷流的底部被加速,并沿着螺旋磁场运动,当成分穿过明亮的射电核心时,就会产生光学、射线以及准同步射电耀变。该源不同波段耀变会表现出
6、延迟。Ahnen等人14发现,2015年2 月,X射线耀变相对于 射线耀变存在约1周的延迟。Raiteri等人16发现,在2019年8月,射电耀变相对于光学耀变有约3周延迟,并认为射电辐射区比光学辐射区更靠近喷流底部。Wehrele等人17认为,高频观测是来自非热光子的逆康普顿散射,而射电辐射来自电子同步辐射,所以两者存在时间延迟。本文第2章介绍我们的观测和数据处理。在3.1节,我们根据多频率VLBA观测结果计算喷流的物理性质;3.2节给出S4 0954+65在2017年5月2021年5月喷流的运动学结果;在3.3 节,我们讨论了射电核心流量变化与 射线光变的联系。第4 章进行总结。在本文中,
7、我们使用的宇宙学常数为:H0=71 km s1 Mpc1,m=0.27,=0.7318。2观测数据及处理分析2.1多频率VLBA观测我们对S4 0954+65 观测是在对近邻低光度活动星系核M8119的VLBA 观测期间作为校准源而进行,观测共分为两个历元。在2018 年6 月10 日,VLBA观测分别在8GHz、22GHz和43GHz交替进行,在每个频率共观测21次扫描(scan),每次扫描为0.5min,总的时间跨度为8h。在2019 年11月2日和11月4日,观测分别在22GHz 和5GHz 及8GHz进行。在每个频率,S4 0954+65 分别观测8次扫描,每次扫描时长为3min,每个
8、频率总的时间跨度为8h。我们采用美国国立射电天文台开发的AIPS软件对S4 0954+65VLBA观测数据进行校准。校准过程采用标准化流程,包括手动相位校准(manual phase cal),全局条纹拟合(globalfringe fitting)以及幅度校准。其中对于幅度校准,我们分别在5GHz和8GHz考虑了电离层改正,在22GHz和43GHz考虑了大气不透明度改正。2.2VLBA 43GHz 档案数据为进一步研究S4 0954+65 的喷流运动学,我们使用波士顿大学研究团队20在43GHz对该源的VLBA监测数据。该项目利用美国VLBA阵在43GHz对数十个明亮的 射线耀变体进行长期监
9、测3。我们对2017年5月2021年5月间的归档数据重新进行了成图,成图以及成分认证采用Jorstad等人3,21的方法。感谢NASA 通过费米客座研究员计划(Fermi Guest Investigator Program)资助的VLBA-BU Blazar 监测计划38(BEAM-ME和VLBA-BU-BLAZAR),为本文研究提供了便利。2 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2592 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2592 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2592.3 射线观测数据Fermi 卫星上的大面积望远镜(L
10、AT)每3h 扫描一遍全天空。本文使用Fermi 卫星4 期源表(4FGL)数据22。我们选用S4 0954+65采样周期为7d的 射线光变曲线数据,每个数据点的不确定度是2,观测时间范围为2017年05月13日2021年05月28日(MJD5788759363)。3结果与讨论3.1射电喷流的物理性质我们按时间顺序把2018年6月命名为第一个历元,把2019年11月命名为第二个历元。第二个历元43GHz数据取自VLBA归档数据。图1中展示了第二个历元S4 0954+65在4个频率的图像。在秒差距的尺度上,S4 0954+65 表现出单边的喷流,喷流的初始方向由北偏向西北并逐渐向西,与千秒差距尺
11、度上的喷流形成近180的夹角23。图2a)展示了射电核在两个历元的谱。由于C3成分(参见3.2节)只有在43GHz能够分解开,所以我们把C3和核成分合并来考虑。由图2a)可知,在第一个历元,射电核成分在我们的频率覆盖范围内表现出平谱,通过相应的射电核成分流量密度变化可知,此时射电核心并不活跃。在第二个历元,核心的谱在低频端表现出同步自吸收特征,为了确定其峰值频率,我们采用Casadio等人242016年5月在86GHz的观测数据,并通过假设其与43GHz之间谱指数不变的条件外推到我们观测的时间。我们利用同步辐射自吸收(SSA)来解释核心在第二个历元的倒转谱(见图2a)。利用SSA谱公式进行拟合
12、25,流量密度可表示为:S=Sm(m)t1 em(/m)t1 em,(1)其中,m是以GHz为单位的倒转频率;Sm是倒转频率处的流量密度,以Jy为单位;t是光学厚发射区的谱指数,为2.5;是光学薄处的谱指数;m=1.5(1 8/3t 1)是倒转频率的光深。通过拟合得到倒转频率m=(452)GHz,Sm=1.3Jy,光学薄处的谱指数=0.1,拟合结果如图2a)所示。我们注意到同步自吸收未能很好地拟合核心成分的谱。一方面这可能是由于低频观测不能很好地分解开核心成分和C3 成分;另一方面,通过86GHz的总流量密度外推得出的核心成分在该频率的流量密度具有较大的不确定性。利用这些物理参数,带入SSA磁
13、场强度公式,公式由Marscher26给出:BSSA=105b()b45mS2m(1+z)(104T),(2)其中,b()是与谱指数有关的函数,由Marscher26列出表1可知,b(0.1)=1.1。核心成分大小为b=0.11mas。多普勒因子可由3.2节表2取8.3,红移z=0.367。将上述参数代入式(2),核心磁场BSSA=(3.3 0.9)104T。260天 文 学 进 展42 卷260天 文 学 进 展42 卷260天 文 学 进 展42 卷注:等强度轮廓线最小值分别为0.24,0.11,0.41和0.76mJy,并以2倍大小递增,总强度的峰值流量密度分别为0.73,0.77,0.
14、80和1.24Jy。每张图左下角的灰色椭圆表示洁束大小。图 1S4 0954+658在5GHz,8GHz,22GHz和43GHz的总强度图像图 2a)中的红色和蓝色圆点分别代表第一个历元和第二个历元的射电核心流量密度,蓝线是第二个历元核心的SSA谱拟合;b)是两个历元射电核心的亮温度2 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2612 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2612 期王玕,等:S4 0954+65 射电喷流的高分辨率研究261另一种计算磁场的方法是假设磁场能量和粒子能量均分,Feng等人27给出磁场与其他物理量的关系式:Beq=4.5(1+k
15、)f(,107,1011)LR32/7,(3)其中,k 是重粒子能量与电子能量比值,通常取k=100,这是由于电子和正电子是在质子流与射电源内的气体和尘埃发生碰撞后产生的次级粒子28。R 是以cm 为单位的成分大小。L是源的同步辐射光度,近似可以写成4d2LSmm,dL是以Mpc 为单位的光度距离,dL=1966Mpc1。f(,107,1011)是谱指数函数,f(0.1,107,1011)=0.6107。式(3)可以表示为如下形式:Beq=5.37 1012Smmd2LR32/7(104T).(4)将前面的参数代入可得核心磁场Beq 0.3G。SSA磁场与均分磁场结果相近。VLBI成分在源静止
16、坐标系中的表观亮度温度由下式得出29:Tb=1.22 1012S2b2(1+z)K.(5)由图2b)可知,第一个历元射电核心亮温度与第二个历元射电核心亮温度相近,亮温度不随频率显著变化,并保持在1011 1012K之间。3.2射电喷流的运动学为了研究喷流的运动学,我们分析了在2017年5月2021年5月观测期间的共37个历元43GHzVLBA归档数据。我们将流量密度小于20mJy极弱成分剔除,得到2个明显的运动成分C1 和C2,和一个靠近核心的静止成分C3。成分C1,C2和C3如图3a)所示。C1存在的时间为2017年8月2018年7月;C2存在的时间为2019年2月2021年5 月;C3 存
17、在时间为2018年4月2021年5月(见图3b)。成分C1,C2和C3在天空平面上运动轨迹如图3a)所示,原点是射电核心,我们假设射电核心静止30;图3b)和3c)横坐标是修正儒略日,纵坐标分别是赤经(图3b)和赤纬(图3c)。通过图3a)可知,C1和C2这两个成分在运动过程中都在正北方向距离核心约0.4mas 处产生转折,并通过图3b)和3c)在转折之前,两个成分都是近似沿正北方向运动;转折后,在赤经方向近似做匀速运动(转折后近似是一条直线);通过在赤纬方向运动轨迹与时间的关系可知,这两个成分在转折前后赤纬方向速度不改变。由此可推断,在转折处,两个成分在赤经方向短暂加速,并继续匀速运动。转折
18、区可能是激波生成区,成分加速可能由激波导致31。由图4可知,C1通过转折区时,流量和成分大小都没有显著变化;而C2通过转折区之前流量突然升高,成分在不断变大,这可能是由于C2穿过C3时,与C3发生碰撞。C3成分类似于Morozova 等人12观测得到的A0,在靠近核心处静止。C1和C2成分的其他结果见表1。我们把37个历元所有成分的分布展示在图5a)中;图5b)展示了Morozova等人12观测的11个运动成分(K1K11)和一个离核心很近的静止成分A0的运动轨迹。在这两个时间范围内,喷流的运动成分表现出两条明显路径,两条路径夹角几乎没有变化,约为30。将262天 文 学 进 展42 卷262
19、天 文 学 进 展42 卷262天 文 学 进 展42 卷图 3a)3个成分在天空平面上的运动轨迹;b)和c)成分的相对赤经和相对赤纬随时间的变化图 4C1,C2和C3的成分的流量密度和尺寸大小随时间的变化2 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2632 期王玕,等:S4 0954+65 射电喷流的高分辨率研究2632 期王玕,等:S4 0954+65 射电喷流的高分辨率研究263表 1S4 0954+65的喷流中观测到的新的射电成分的特征成分历元数平均位置角(PA)/()平均大小(FWMH)/mas最大流量/mJy视速度app/c生成时间C11013.7 4.10.14
20、0.0381 710.4 0.32017.05.23 9C21914.2 6.40.23 0.08471 288.2 0.22018.12.02 8图3a)与图5b)对比可知,除K2和A0成分,Morozova等人12观测结果的其余成分与本文的C1和C2成分很类似,都是先沿着正北方向运动,再向西转折。根据Morozova等人12观测结果可知,K2和K3的存在时间分别约为500d和700d,生成时间相差约3个月,这表明两个成分在两条路径上同时运动。图 5a)本次工作所有历元的成分分布,不同历元的成分用不同颜色表示;b)Morozova等人观测得到的12个成分的运动轨迹12S4 0954+65 中
21、心超大质量黑洞的质量为MBH=3.3 108M6,质量与史瓦西半径的关系rs=2GM/c2。图3 是运动成分在天空平面的二维投影,1mas=5.08pc,利用视向夹角范围为1.59.3(本节稍后给出),可算出转折点到射电核心的距离为7.8 10440.3 104rs,转折区在邦迪半径附近。在Xu等人32的工作中,BL Lac天体S5 0716+714也存在类似于S4 0954+65的双路径结构。Kim等人33发现S5 0716+714的喷流的运动成分的分布在距离核心0.53mas处发生转折,转折点到射电核心的距离约为邦迪半径。双路径喷流结构和运动路径的转折区发生在邦迪半径附近,这可能与黑洞影响
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- S40954 65喷流的高分辨率观测研究 65 喷流 高分辨率 观测 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。