怎样读一篇薄盘光谱论文

把第 3.2 节、式 14–17、附录 E 和第 4.2 节放在手边。先读模型曲线的坐标与图注，再追到发射和传播公式，最后看拟合还需要哪些外部测量。仅凭光谱“更硬”，自旋还定不下来。

先分清纵轴数的是光子还是能量

图 5 的横轴是远方观察者测得的光子能量，单位千电子伏特；纵轴是光子数谱通量密度，单位为每千电子伏特、每平方厘米、每秒的光子数。两条轴都取常用对数。横轴坐标为零，对应一千电子伏特；向右走一个单位，能量就增大十倍。

把光子能量记为 $E$，光子数谱记为 $N_E$。若统一用焦耳、平方米和秒，$N_E$ 就是每平方米、每秒、每焦耳内的光子数。乘上每个光子的能量，得到能量谱通量密度：

同样多的光子，高能光子送来的能量更多。换成频率谱时，用 $E=h\nu$，其中 $h$ 是普朗克常量，$\nu$ 是频率。同一小段光谱携带相同的能量，满足 $F_\nu\,d\nu=F_E\,dE$，于是

这个能量因子会改变曲线形状。数光子的峰和按对数频率区间数能量的峰，可以落在不同位置。拿两张光谱图作比较，先把各自的纵轴写下来。

试着换一次单位：在一千电子伏特处，光子数谱通量密度若为每千电子伏特、每平方厘米、每秒一个光子，相应的 $\nu F_\nu$ 是多少？一千电子伏特等于 $1.602176634\times10^{-16}\,\mathrm{J}$，一平方厘米等于 $10^{-4}\,\mathrm{m^2}$。把光子能量和能量区间一起换算，结果是

这里不能只乘能量，漏掉“每千电子伏特”的区间单位。电子伏特换算和普朗克常量可查 NIST 的 SI 常量表。

上下两幅各只改了什么

先读上幅图注。论文用几何单位写无量纲自旋 $a/M$，依次取 $0,0.5,0.9,0.999$。观察方向与盘面法线夹角固定为 $30^\circ$。其余条件是质量十个太阳质量、距离十千秒差距、向内吸积率 $10^{19}\,\mathrm{g\,s^{-1}}$、零内缘力矩，以及颜色硬化因子 $f_{\rm col}=1$。第 3.2 节还说明，图 5–8 都假定气体局部向各方向等强度发光。

自旋越大，曲线的高能端越向右伸。在论文的模型里，盘内缘位于最内稳定圆轨道；顺行自旋越快，这条轨道越靠近中心，盘便能延伸到更热的内区。第 3.2 节同时讨论了返回盘面的光：引力把一部分内区辐射弯回盘面，连较远环带的供光也会受影响。解释这组曲线，需要同时跟着盘的发射和光的传播往下看。

下幅固定 $a/M=0.9$，把视角依次改为 $0^\circ,40^\circ,70^\circ,85^\circ$，其余图注参数不变。低能处，盘越倾斜，曲线越低；到了高能端，排序变了。假如视角只把可见面积乘一个余弦，整条曲线会按同一比例缩放，曲线不该换序。作者把高能增强归因于多普勒束射与引力聚焦；这正是简单投影图像需要补上相对论计算的地方。

“更硬”只说明高能光子的相对贡献增大。温度、视角和其他模型假设都可能改变这种相对贡献。先说清这次只改了哪个条件，才有具体的物理解释。

沿附录 E 从天空找到发光的气体

式 E1 把观察者天空上各方向的强度相加。这里 $F_E$ 是收到的能量谱通量密度，$I_E$ 是每单位立体角的比强度，$d\Omega$ 是天空上的一小块立体角，单位球面度。式 E2–E3 再引入频率比 $g$，也就是收到的光子能量除以发射时的光子能量：

读这个式子，要把三次方因子和括号里的能量一起看。观察者在某个能量收到的光，出发时可能处在另一个能量区间；只改强度系数，会漏掉这次频移。

接着读式 E4–E6。有效温度由盘面发出的总通量定义，颜色温度决定热谱的形状；颜色硬化之后，再用四次方稀释恢复原来规定的总通量。式 E7 把这些关系放进同一个天空积分。每个像方向可能遇到不同温度，也可能有不同频移，最后按它在天空所占的立体角相加。这个立体角已经包含盘面呈现在眼前的几何，不能再随手乘一次平盘的投影因子。

可对照 Zimmerman 等，astro-ph/0408209v1 的第 2 节，看牛顿薄盘怎样从零内缘力矩得到局部温度。本课程的例子沿用这条供能律，再接史瓦西传播；Li 等研究的是 Kerr 薄盘，并包括返回盘面的辐射。课程里的时空固定为不自转，改变规定的内缘，可以比较供光范围和可见区域怎样变化。若要把内缘解释成自旋，还需要前后一致的 Kerr 发射与传播模型。

读到第 4.2 节，再谈自旋推断

连续谱拟合利用实测光谱中平滑的热辐射部分，估计发光内区的尺度。要把这个尺度连到自旋，还得假定盘内缘与最内稳定圆轨道有怎样的关系。质量 $M$ 决定引力半径 $GM/c^2$ 有多长，其中 $G$ 是引力常量、$c$ 是真空光速；距离 $D$ 决定同一份接收通量对应多大的本征供光；倾角 $i$ 则影响投影面积和相对论传播。每个量都需要独立测量及其不确定度。

第 4.2 节可以逐句练习这种读法。前一节比较归一化时，作者固定自旋，让距离参与拟合；进入自旋比较后，作者采用独立的距离估计，再让自旋和吸积率变化。先看“这次放开了哪些参数”，再看拟合出的自旋。XSPEC 的 KERRBB 官方说明也把质量、距离、倾角和颜色硬化因子分别列为输入，并单列内缘力矩、返回辐射及临边变暗的选择。

盘大气也是解释的一部分。吸收与散射会改造气体发出的谱形；用一个颜色硬化因子概括这些过程，是一种近似。论文第 4.2 节专门讨论了它怎样影响自旋估计。图 5 取 $f_{\rm col}=1$，本课程例子取 $1.7$，这两套条件不能直接当成同一组自旋比较。真正拟合观测时，还会用到仪器响应、沿途吸收和其他光谱成分；图 5 本身给出的是控制参数后的模型曲线。

铁线反射先问光从哪里来

连续谱拟合追踪盘自身的热辐射。反射谱还要跟着另一束光：高温日冕的X射线照到盘上，经过吸收、荧光和散射后离开。铁荧光在约$6\text{–}7\,\mathrm{keV}$附近留下谱线，轨道运动和引力又改变线形。拟合时，气体电离状态与照明方式也要写清。RELXILL作者说明把局部反射与相对论传播接在一起。

打开García等人的arXiv:1312.3231v2，2013年12月12日版，读第1至4.1节，再看PDF第8页的图7。左图比较四种观察倾角下的归一化能量谱；中图用保留发射角的预测，除以通常的角度平均卷积结果。这是两种模型的比值，不是观测残差。右图列出各个局部发射方向对收到的光贡献多少。光路会弯，一个观察者便能收到朝多个方向发出的光。

试着判断：中图的比值随能量变化，只调整整条谱的高低能消除它吗？不能。读自转拟合时，把角度处理、倾角、电离、铁丰度和照明方式一并记下。课程中的热谱例子帮助比较发射与传播，反射计算由所引模型另行提供。

换一张图，自己试着读

光子数谱和能量谱的峰不在同一处，哪张错了？ 单凭这一点判断不了。先换成同一种谱量和单位，再比较峰的位置。

两条倾角曲线交叉，一个投影面积因子够吗？ 常数因子会把所有能量同时按比例缩放，解释不了交叉。接着找依赖光子能量的多普勒频移，以及引力怎样改变盘的像。

拟合把较热的光谱解释成较快自旋，要查什么？ 先列出固定的参数，再查它们来自哪些外部测量、误差多大。内缘位置、内缘应力和大气假设也要放在同一页看。

把中央图像裁得更紧，完整光谱会变吗？ 只裁显示画面，原来的积分不变；实际收光孔径缩小，收集到的光才会改变。比较图像和光谱前，先问清积分覆盖哪里。

下一次读论文，留一张小纸条。记下版本与页码、坐标量及其单位、变化参数、固定物理假设和积分范围，再写清拟合用了哪些外部输入。把这些写全，再比较曲线，就能说清图中哪一步支持了结论。