时雪草 潘 翔 窦立明 王建国 姜 鹏 杨臣威 周宏岩

(1 中国科学技术大学天文学系合肥230026)

(2 中国极地研究中心上海200136)

(3 广州大学天体物理中心广州510006)

(4 中国科学院云南天文台昆明650216)

1 引言

星际介质是星系的重要组成部分, 很大程度上决定了恒星和行星的形成, 因而对于天体起源研究尤为关键[1].尘埃是星际介质的主要组分之一, 与星际环境频繁相互作用:(1)对入射光进行吸收、散射、再发射, 从而改变辐射场(如星光、吸积盘、超新星等产生的辐射); (2)通过耗散效应, 沉积周围气态物质元素, 催化促进分子生成; (3)碰撞、聚集使尘埃尺寸逐渐增长, 为天体形成创造条件[2].

消光曲线, 即尘埃对光子的吸收强度随波长的变化特征, 是观测上探测和研究星际尘埃的重要手段.尘埃消光在紫外-光学波段最为显著, 而在这一波长范围内最重要的观测特征为2175波长处的宽吸收峰.在本星系群中, 该特征在质量较大、金属丰度较高的银河系中尤为显著[3], 因而也常称遥远宇宙中的2175吸收体为类银河吸收体.除银河系外, 大麦哲伦云[4]及某些遥远漩涡星系[5]的部分视线上同样观测到2175吸收峰,可见该特征在星际尘埃中较为普遍.但是2175吸收峰的载体尘埃到底是什么化学成分、具有怎样的物理性质目前尚无定论.

近期, 我们在7个类星体中发现疑似与蓝移宽吸收线(外流气体)相关联的类银河尘埃吸收特征, 暗示2175吸收体可能存在于类星体外流气体中[8].在另一个类星体SDSS J170542.91+354340.2 (简称J1705+3543)的系统红移处, 我们同样发现了显著的类银河吸收峰, 结合残余莱曼辐射、吸收线光变等特征分析, 该吸收体很可能内禀于类星体[9], 而遗憾的是并没有获得J1705+3543在X射线波段的观测数据, 因而难以对该类星体在高能波段的辐射特性进行限制.本文给出了类星体SDSS J091613.60+292106.1 (简称J0916+2921, 系统红移zem= 1.1418 ± 0.0018)系统红移处的类银河2175消光特征及其光学和X射线波段观测, 弥补了之前研究中的部分不足.

如无明确说明, 本文中误差为1σ置信水平; 同时在红移距离计算中, 本文假定标准ΛCDM宇宙学模型, 参数为哈勃常数H0= 70 km·s−1·Mpc−1、暗能量密度参数ΩΛ= 0.7、物质密度参数Ωm= 0.3.

2 光学观测与数据分析

斯隆数字巡天(Sloan Digital Sky Survey, SDSS)于2003年12月21日对类星体J0916+2921进行了测光观测, 并于2005年1月18日对其进行了光谱观测[10].类星体J0916+2921在SDSS的u、g、r、i、z波段星等分别为19.99 ± 0.04、19.55 ± 0.02、18.74 ± 0.02、18.53 ± 0.02、18.45 ± 0.04; SDSS光谱波长覆盖为3800–9200, 分辨率R ∼2000.在近红外波段, 英国红外望远镜深场巡天(UKIRT Infrared Deep Sky Survey, UKIDSS)测量的Y、J、H、K波段星等分别为17.77 ± 0.02、17.43 ± 0.02、16.86 ± 0.04、16.10 ± 0.03, 观测时间在2007年2月5日至2009年4月3日之间[11], 无显著光变.

随后, 我们使用6.5 m多镜面望远镜(Multiple Mirror Telescope, MMT)的蓝通道光谱仪(Blue Channel Spectrograph, BCS)于2008年2月11日和3月30日对J0916+2921进行了两次光学光谱观测.观测中分别使用1200 grid·mm−1和300 grid·mm−1光栅, 光谱分辨率分别约为2000和1000, 波长覆盖范围分别为4160–5480和3180–8384.两条MMT光谱中, 前者(1200 grid·mm−1)具有较高分辨率, 同时覆盖了绝大多数吸收线, 使得吸收线柱密度的测量更为可靠; 后者(300 grid·mm−1)相对于SDSS光谱向蓝端延伸, 使得静止系紫外波段消光曲线测量更为可靠.我们参照银河系尘埃分布[12], 对所有光谱、测光数据进行了银河系消光改正.所有光谱流量与SDSS测光流量大体吻合.为提高信噪比, 我们将所有光谱合并, 并测量消光曲线, 而吸收线测量则使用高分辨率的MMT 1200 grid·mm−1光谱和SDSS光谱.

2.1 消光曲线

如之前所介绍, 我们合并J0916+2921的所有光学光谱, 以研究其消光曲线.假定高光度类星体合成谱[13]为J0916+2921的无吸收模板.将类星体模板在K波段流量调整到J0916+2921的观测值后(图1中天蓝色实线), 直观对比可以看出, J0916+2921的测光(绿色圆点)和光谱(黑色实线)流量在短波长处受到强烈抑制, 同时在静止系2200附近存在明显的类银河吸收峰.随后, 我们将在K波段与J0916+2921归一的类星体合成谱模板插值到观测波段, 得到无吸收模型单位波长(λ)的流量结合观测流量计算得到实测消光曲线(图1子图中黑色实线和绿色圆点).这条观测消光曲线整体较为平滑, 表明我们选用的合成谱模板很好地对目标类星体J0916+2921的无吸收谱进行了近似, 即使在发射线附近也没有很强的局部特征, 因而对于消光曲线的分析影响不大, 尤其是附近没有强发射线的2175吸收峰特征.另一方面, 吸收线附近区域呈现很窄的尖峰, 由于特征很窄, 同样不会对平滑的连续谱消光曲线产生干扰.

图1 J0916+2921的静止系紫外–光学–近红外测光(绿色圆点)、光谱数据(黑色实线)与类星体合成谱(天蓝色实线)的对比(在K波段归一化).以合成谱为模板计算得到消光曲线(子图中黑色实线与绿色圆点), 并得到最优拟合(红色实线).Fig.1 The UV–optical–NIR photometric (green dots) and spectroscopic (black solid lines) data of J0916+2921 are shown, together with the X-shooter quasar composite (sky blue solid line) scaled to the flux of J0916+2921 at the K band.The apparent extinction curve is obtained assuming the quasar composite as an extinction-free model (black solid lines and green dots in the insert panel), and the best–fit model is shown with red solid lines.

鉴于J0916+2921的消光曲线与银河系类似, 我们选用银河系中常用的参数化模型:

其中λ为波长, 以µm 为单位, x为波数, 单位为µm−1.该分段模型由3种成分组成: (1)紫外波段的线性消光成分, 斜率为c2, 标度为c1; (2)以Drude轮廓描述的2175吸收峰: 中心波数x0, 吸收峰宽γ, 表达式为

系数为c3; (3)远紫外的2阶消光项, 以c4(强度)和c5(2阶项起始位置)为参量.即便吸收线、发射线对消光曲线测量的干扰较小, 我们仍然对强吸收和强发射线中心区域进行了屏蔽, 随后用参数化模型拟合J0916+2921的观测消光曲线, 并用自举法重复1000次测量以计算模型参数的误差.最后的拟合参数为c1= −0.13±0.15、c2= 0.18±0.06、c3=1.78±0.56、x0=4.72±0.03、γ =1.53±0.20、c4=0.40±0.15和c5=5.03±0.23.吸收峰中心波数略高于银河系的均值, 且宽度较宽, 类似于大部分类星体视线上实测的中间插入2175吸收体[14].吸收峰强度为Abump= πc3/2γ =1.89±0.09.

2.2 吸收线与柱密度

鉴于吸收线较窄, 我们可采用局部平滑、迭代屏蔽的方式对观测光谱进行归一化处理[14].在光谱中, 我们探测到两组MgII吸收线系统, 其中一组为中间插入吸收, 红移为0.8376, 由于与本文研究目标没有直接相关性, 因而不作讨论; 另一组吸收线系统红移为zabs=1.1413±0.0002, 与系统红移zem=1.1418±0.0018一致, 为类星体关联吸收体.上节中显著的2175尘埃吸收峰表明类星体附近应该存在显著的低电离吸收线系统, 而该近邻吸收体是系统红移处唯一的强吸收线系统, 因而无疑两者有物理联系, 所以我们在后文中称该吸收体为2175尘埃吸收体.该尘埃吸收体的主要强吸收线轮廓见图2,其中灰色方块为SDSS光谱数据, 黑色散点为MMT光谱数据, 两次观测吸收轮廓接近,且MMT数据质量有明显提升.可以看到在尘埃吸收体静止系(zabs=1.1413±0.0002)中,吸收线主要覆盖的速度区间为[−200, +200]km·s−1, 因而我们考察这些证认的吸收线在这一范围的光深以及等值宽度, 通过视光深法(Apparent Optical Depth, AOD[15])和生长曲线法(Curve Of Growth, COG[16])分别估计离子的柱密度.在使用生长曲线法之前, 我们需要得到尘埃吸收体的等效速度展宽b.考虑到我们探测到多个FeII跃迁的吸收,其中有FeII λ1608这样非常微弱的吸收线, 也有FeII λ2382这样吸收接近见底的强线, 通过对FeII等值宽度的模拟, 可以准确限定生长曲线.最终拟合结果见图3, 卡方分布图表明等效展宽为b=34.5±3.0 km·s−1.基于该生长曲线,我们同样对其他离子柱密度进行测量.视光深法和生长曲线法测得的柱密度见表1中lg NAOD及lg NCOG列, 对于弱线, 两者测量结果较为接近, 而对于强线, 由于存在潜在的饱和区域, 视光深法测量并不可靠,因而我们采用生长曲线法的结果作为离子柱密度.通过累加同种元素的不同离子柱密度, 我们对元素柱密度进行了估计(表1中lg NX列).鉴于锌的升华温度低难以附着尘埃,因而常用于表征吸收体气体元素原始丰度.若以太阳丰度[17]为基准, 则可以计算得到吸收体中Al、Cr、Fe相对锌元素的丰度[X/Zn]([X/Zn]≡lg(NX/NZn)−其中lg NX为实测J0916+2921的元素丰度, lg为太阳元素丰度[17], 此定义中太阳的元素相对值丰度为0), 结果为: [Al/Zn]=−1.68±0.10、[Cr/Zn]=−0.49±0.10与[Fe/Zn]= −0.81±0.18.这种气体中元素大量缺失的现象是由于金属离子附着尘埃形成, 因而称之为尘埃耗散.该吸收体中的尘埃耗散严重表明系统中尘埃丰富, 进一步验证了气体吸收线与观测的显著尘埃消光特征之间的关联.

图2 J0916+2921近邻吸收体中主要吸收线的归一化轮廓, 已经投影到速度空间, 横轴∆V 表示相对于吸收线系统速度(红移zabs = 1.1413±0.0002)的偏离, 负数为蓝移.其中黑色散点为MMT数据, 灰色方块为SDSS数据, 绿色点线为归一化续谱, 红色虚线标记[−200, +200]km·s−1范围.Fig.2 The absorption line profile for major absorption lines in the proximate absorption line system of J0916+2921.The relative velocity (∆V) in the panels is calculated in the rest-frame of the absorber zabs = 1.1413±0.0002, and the negative value means blue shift.Black dots are derived using the MMT data, and grey squares correspond to the SDSS data.The green dot lines are the normalized continuum,and the red dashed lines mask the [−200, +200]km·s−1 velocity range.

图3 右图: 利用多普勒参数b和多种FeII离子的柱密度拟合的生长曲线的约化卡方轮廓, 红点为最佳拟合, 从内到外的实线分别对应1σ、2σ、3σ置信区间.其中NFeII为FeII离子柱密度, b为谱线多普勒速度展宽; 左图: 各吸收线在最佳拟合生长曲线(黑色实线)上的位置, 其中灰色实线表示生长曲线1σ的误差范围.其中N为离子柱密度, f为谱线振子强度, W 为谱线等值宽度.Fig.3 Right panel: The contour of the reduced χ2 in the fitting of COG using the doppler parameter b and FeII multiplets.The best-fit result is indicated with a red dot, and the rings from inside out correspond to 1σ, 2σ, and 3σ confidence levels.Here, NFeII is the column density of FeII.Left panel: The positions of different absorption lines in the best-fit COG (black solid line) are shown, where N is the column density, f is the oscillator strength of the absorption lines, and W is the absorption line equivalent widths.The grey solid lines show the 1σ uncertainty of the COG.

3 X射线观测与数据分析

表1 J0916+2921近邻吸收体中主要吸收线及测量结果Table 1 Major proximate absorption lines in J0916+2921 and measurement results

对于抽取得到的光谱, 我们以每个单元最低20个计数进行分割后合并.鉴于MOS1和MOS2仪器性质接近, 且实测光谱相吻合, 我们将MOS1、MOS2光谱作合并处理.对于所抽取得到的PN及MOS光谱, 我们用XSPEC (12.6.0版本)软件同时进行拟合,拟合误差为90%置信度水平.拟合过程中采用同样的模型, 并对PN与MOS光谱的参数(除归一化参数外)进行绑定.首先,我们尝试了最简单的幂律谱模型(XSPEC的zpowerlw,并取zem= 1.1418 ± 0.0018), 同时考虑银河系中性吸收(柱密度为= 1.76 ×1020cm−2[18]).结果见图4上, 以卡方/自由度表示的最佳拟合优度为χ2/DOF =14.09/28, 拟合得到幂指数为Γ = 1.74 ± 0.13, 同时银河吸收改正后的0.3–10.0 keV X射线亮度为1.37 ×10−13erg·s−1cm−2(PN)及1.32× 10−13erg·s−1·cm−2(MOS).鉴于J0916+ 2921中存在显著的内禀气体吸收和尘埃消光, 我们尝试加入本征吸收体对光谱进行拟合.在幂律谱的银河吸收基础上, 我们添加zabs= 1.1413±0.0002处的中性吸收体.由于之前所测量的吸收线均来自金属离子而没有氢、氦元素的柱密度估计, 我们对吸收体的金属丰度并不了解.考虑到吸收体很可能与类星体或其寄主星系相关联以及对类银河消光曲线的探测, 本征吸收体很可能类似银河系星际介质, 即大约为太阳丰度, 因此我们同样假定X射线的本征吸收体为太阳丰度.最终拟合优度与之前相仿,为χ2/DOF = 14.09/27, 因而光谱质量无法精确计算吸收体性质.本征吸收模型拟合得到的幂指数为与之前一致, 同时本征吸收体柱密度限制为1021cm−2, 置信度为90%.若假定太阳丰度, 这一柱密度值对应的Zn元素柱密度限制为实测Zn元素柱密度在X射线模拟的估计范围内, 这一结果表明之前的假定, 即吸收体可能对X射线产生吸收以及吸收体为太阳丰度, 在现有数据基础上是合理的.拟合得到的卡方对幂指数、本征吸收体柱密度的依赖关系见图4中的子图.

图4 上图: J0916+2921的XMM–Newton EPIC光谱(散点)以及最佳拟合结果(阶梯实线).其中PN和MOS的数据及模型分别由黑色和红色表示.子图显示了卡方对拟合参数X射线谱指数Γ及本征吸收体柱密度的依赖, 从内到外的环线分别表示1σ、2σ、3σ置信区间.下图: 上图中最佳拟合的卡方.Fig.4 Upper panel: The XMM–Newton EPIC spectrum of J0916+2921 and the best–fit model of a simple power law.The upper and lower spectra correspond to the PN and MOS data.The inset shows the confidence contours (1σ, 2σ, and 3σ for the two interesting parameters) for the photon index Γ and the intrinsic neutral absorption column density.Lower panel: residuals of the best–fit model as in the upper panel.

4 讨论: 类星体高能辐射环境下的类银河系尘埃?

结合测光流量, 我们估算得到J0916+2921在消光改正后的光学波段5100单色光度为L5100= 4.46×1045erg·s−1.另一方面, 在X射线波段, J0916+2921的观测系下0.3–10 keV光度为L0.3−10keV=2.10×1045erg·s−1.X射线–光学光度比为L0.3−10keV/L5100= 0.47, 对于常用类星体SED[19], 比值为0.50, 两者相一致.注意到我们对光学流量进行了消光改正, 而受限于数据质量, 对X射线流量并未修正潜在吸收体的吸收, 可见J0916+2921的实际X射线辐射强度高于直接观测结果, J0916+2916的高能辐射很可能相对一般类星体更强.按照一般AGN中热光度Lbol与L5100经验关系[20], 可以得到Lbol= 4.60 × 1046erg·s−1= 1.20 × 1013L⊙.J0916+2921光度较高, 远高于普通星系(如银河系、仙女座星系光度∼2×1010L⊙).若假定尘埃升华温度为T1500K= T/1500K,则J0916+2921的尘埃升华半径为= 1.6 pc[21], 其中LUV,46为类星体紫外光度, 单位是1046erg·s−1.

类星体J0916+2921具有较高的光度, 同时具有比一般类星体相对更强的高能X射线辐射.在这样的高能辐射环境下类银河2175吸收体是否存在, 是类星体及星际介质的研究中一个十分有趣的问题.J0916+2921中近邻尘埃吸收体有3种可能的起源[15], 与类星体相关性由远及近分别为: (1)处于类星体伴星系中; (2)来源于类星体寄主星系的恒星形成区; (3)来源于类星体核区.首先, 在近邻宇宙中, 2175吸收体普遍存在于漩涡星系(如银河系、仙女座星系)的星系盘中, 其所处的星际环境中通常金属丰度较高且尘埃丰富, 而伴星系(如大、小麦哲伦云)通常质量较低, 且金属元素、尘埃成分贫乏, 因而伴星系不太可能是2175吸收体的载体.因而J0916+2921中2175吸收体更可能位于寄主星系或内禀于类星体核区.对于(2)、(3)两种情况,基于已有的观测数据并不能加以区分.若吸收体位于寄主星系恒星形成区, 考虑到一般漩涡星系的尺寸, J0916+2921中吸收体距离星系核的距离至多为几十个kpc (千秒差距).与宁静星系相比, J0916+2921的恒星形成区会受到较强的来自核区的高能电离辐射, 同时类星体对寄主星系的外流反馈作用也会对星际介质产生扰动.这类环境下星际介质的性质会受到怎样的影响以及怎样进一步影响恒星形成过程都值得进一步研究.若吸收体内禀于类星体, 吸收体无法在尘埃升华半径以内存在, 因而距离中央引擎须超过1.6 pc.与之相关联的类星体核区结构可能有尘埃环、窄发射线区或者外流.尘埃环和窄发射线区由于受类星体引力场束缚, 其速度场有较大可能性为维里化状态, 因而吸收体红移很可能与系统红移非常接近,同时吸收线轮廓也偏向于对称结构; 而外流对应的吸收体通常相对类星体系统整体蓝移(吸收体红移低于系统红移), 同时吸收轮廓通常更为复杂.由于吸收线较窄, 我们的光谱观测分辨率不足以对吸收线轮廓进行分辨, 同时波长覆盖范围内也没有很好的系统红移处指针(如窄发射线或星光吸收线), 因而目前无法对这3种起源加以区分.但无疑, 尘埃吸收体的内禀起源局限于类星体外围结构, 因而了解吸收体的性质可能对研究类星体外区的结构、演化等过程提供帮助.

J0916+2921中探测到近邻的类银河吸收体对2175吸收载体能否在高能辐射场中存在提出了挑战.即使在赛弗特星系中, 理论上即使距中央引擎千秒差距的距离上, 核区的高能X射线和极紫外光子也可以轻易瓦解PAH分子[22], 而在光度更高的类星体J0916+2921中更是如此.而实际观测中, 90%的星爆星系中探测到的PAH发射特征,只出现在约10%的活动星系中[23], 在某种程度上印证了活动星系中高能辐射摧毁PAH分子的理论预期.但另一方面, 我们也不难看出, PAH发射特征在活动星系中并非没有探测, 一定数量的例外仍然存在.同时我们也发现了一批类星体近邻的2175吸收体, 如引言中提到的J1705+3543以及部分宽吸收线类星体.实际上, 尘埃是类星体的重要结构—尘埃环的主要成分, 而碳基小颗粒尘埃则是大量尘埃颗粒的根源之一, 因而不可或缺.那么这些小颗粒尘埃如何在活动星系的辐射场中生存, 很大程度上可能源于遮蔽效应[22], 即在某些视线上, 类星体电离辐射到达尘埃之前便已经被足够厚的物质所吸收,使这些受遮蔽的小颗粒尘埃免受高能射线的解离.这些存活的小颗粒尘埃如何进入类星体周围, 它们的进一步演化以及与周围物质的相互作用都是很有意义的课题, 值得后续的观测研究和理论分析来进行解答.

5 总结与展望

本文中, 我们对一个特别类星体J0916+2921的近邻类银河2175尘埃吸收体进行了详细的光学、X射线波段的联合诊断.通过分析光学测光和光谱数据, 我们对其消光曲线进行测量, 探测到强度为Abump= 1.89±0.09的类银河2175吸收峰.其吸收峰相比银河平均值较宽, 与之前发现的河外2175吸收特征类似.在光学光谱中,我们发现系统红移附近唯一的一个近邻MgII吸收系统.我们测量了吸收体的气体元素柱密度, 以Zn元素为基准, 发现Al、Cr、Fe元素均存在严重的尘埃耗散现象.种种迹象均表明该吸收线系统与2175吸收峰相关联.通过对X射线和光学测光流量进行分析, 计算得到J0916+2921的X射线光度为L0.3−10keV= 2.10× 1045erg ·s−1, 热光度为Lbol= 1.20×1013L⊙, 具有典型的类星体高能辐射特征.若考虑X射线波段存在潜在的近邻吸收体吸收, J0916+2921的X射线辐射相较一般类星体更强, 远高于正常星系.与之前的近邻2175吸收体类似, J0916+2921的类银河吸收体可能来自寄主星系恒星形成区, 也可能内禀于类星体.无论是哪种可能性, J0916+2921中的强X射线辐射都对近邻2175吸收体的存在提出了更严重的挑战.对J0916+2921这样本身具有很强高能辐射而兼有近邻2175吸收的目标, 后续多波段观测及大样本分析对于研究类星体环境下尘埃的性质、尘埃对类星体的演化、类星体与寄主星系的相互作用和共同演化都具有独特的意义.