读懂Chandra对NGC 1275的轴子上限

我们通读的是arXiv第4版全文，不是期刊版，两者没有逐字比对。我们没有复现这条上限。Lab里“百分之一量级的起伏”预设用发表的耦合值和示意性磁场参数，只显示效应有多大。

图3：光谱有多平滑？

先看左图。上半部分是Chandra两路光栅的一阶谱合在一起，HEG为蓝色，MEG为红色，叠着最佳幂律模型；下半部分是数据除以模型。右图把这些比值汇成直方图，点线标出正负百分之三。

拟合模型是一条幂律，光子指数 $\Gamma=1.890\pm0.009$ ，被银河系气体吸收，柱密度 $1.32\times10^{21}\,\mathrm{cm^{-2}}$ 。拟合统计量 $C=4956$ ，自由度4874。源的观测时长约490千秒，使用1到9 keV波段。

常见误读是把图上的点当成拟合用的数据。图注写明，数据为了展示做了大幅合并，拟合用的是未合并的数据。百分之三那两条线只帮人看清散布有多大。

式(2)和两种磁场模型

电子密度用早先X射线工作给出的双成分剖面，它决定每个半径上的等离子体项：

n_e(r)=\frac{3.9\times10^{-2}}{[1+(r/80\,\mathrm{kpc})^2]^{1.8}}+\frac{4.05\times10^{-3}}{[1+(r/280\,\mathrm{kpc})^2]^{0.87}}\,\mathrm{cm^{-3}}.

磁场由一块块小区搭成。模型A里每块的磁场按 $B_0[n_e(r)/n_e(0)]^{0.7}$ 缩放， $B_0\approx25\,\mu\mathrm G$ ；每块尺寸在3.5到10 kpc之间抽取，概率正比于 $L^{-1.2}$ 。中心10 kpc被排除，所以小块从10 kpc一直排到维里半径1.8 Mpc；只保留星系团法拉第旋转不超过 $2000\,\mathrm{rad\,m^{-2}}$ 的组态。模型B里磁压固定为气体压强的百分之一，25 kpc处磁场约7.5微高斯，按等温星系团处理时随密度的平方根下降。

论文明说，无论小块模型还是高斯随机场，都不是自然界真实存在的磁场；对转换要紧的是横向强度和沿视线的相干长度。本系列的星系团磁场与等离子体中的混合沿模型B剖面算过临界能量：在这个上限耦合下它远高于9 keV，英仙座处在弱混合区。

图4：一个实现，许多参数

两幅图都用模型B的同一个代表性实现。左图质量固定在 $\log_{10}(m_a/\mathrm{eV})=-12.7$ ，耦合从 $\log_{10}(g_{a\gamma}/\mathrm{GeV^{-1}})=-11.9$ 一路降到 $-12.7$ 。右图耦合固定在 $-12.1$ ，质量从 $-12.3$ 变到 $-12.9$ 。

左图里耦合越大凹陷越深，这符合弱混合的预期，转换大致随耦合的平方增长。右图里质量越大，低能端越被压住，这就是临界能量的表现。不要把这些曲线读成平均值：每条都是单个磁场实现，起伏的位置只属于这个实现。

图5：磁场模型与磁场实现

这里取一个参数点， $m_a=10^{-12}$ eV、 $g_{a\gamma}=10^{-12}\,\mathrm{GeV^{-1}}$ ，每种模型画两个实现，模型A为黑色，模型B为红色。4 keV以下两种模型给出的畸变大小相近，模型A的更窄；4 keV以上两者都变得更有周期性，模型A的畸变明显变大。

要点是：同一模型的两个实现之间的差别，可以和两种模型之间的差别一样大。任何上限都必须对实现做边缘化，这是下一步的内容。

式(3)和式(4)：对磁场做边缘化

作者为每种模型生成500个合格的磁场组态，在一张网格上计算存活曲线： $\log_{10}(m_a/\mathrm{eV})$ 从 $-13.6$ 到 $-11.1$ ， $\log_{10}(g_{a\gamma}/\mathrm{GeV^{-1}})$ 从 $-13$ 到 $-10.7$ ，共约26万条曲线。每条曲线乘上被吸收的幂律，再去拟合未合并的光谱，对第 $i$ 个实现得到统计量 $C(m_a,g_{a\gamma},i)$ 。

后验正比于 $\exp(-C/2)$ ，按式(3)在质量、耦合和实现上归一化；按式(4)对实现求和，就得到只关于质量和耦合的后验。先验在 $\ln m_a$ 和 $\ln g_{a\gamma}$ 上是平的。未知磁场正是通过这次求和进入结果的：只有几乎所有可信的磁场组态都没法让某个耦合藏在数据底下，这个耦合才被排除。

图6和式(5)：上限本身

上图是模型A，下图是模型B。粗黑线围出99.7%置信度的排除区，涂浅蓝；黄色是95%的允许区，红点是最佳拟合点。图里还画了更早的上限：M87（Marsh等2017）、SN1987A（Payez等2015）和Fermi，以及IAXO的预期灵敏度。

对 $m_a<10^{-12}$ eV，论文在99.7%置信度下给出：模型A为 $g_{a\gamma}<6.3\times10^{-13}\,\mathrm{GeV^{-1}}$ ，模型B为 $7.9\times10^{-13}\,\mathrm{GeV^{-1}}$ 。式(5)把模型B的值列为主要结果。

最容易误读的是那块黄色的“岛”。在95%水平上，模型B偏好 $g_{a\gamma}\approx4$ 到 $8\times10^{-13}\,\mathrm{GeV^{-1}}$ ，起因是1.2 keV处一个不大的凹陷和3到4 keV的一个宽隆起。作者认为它们很可能来自光栅有效面积的残余定标误差，因此只把结果当作上限。

四类陈述

已发表结果。 Reynolds等2020年对模型A和B给出的99.7%上限是经过审稿的结果。同样属于此类的有Ayala等2014年的水平分支上限 $g_{a\gamma}<6.6\times10^{-11}\,\mathrm{GeV^{-1}}$ ，以及Ajello等2016年用同一星系六年数据得到的Fermi-LAT上限：0.5到5 neV之间 $5\times10^{-12}\,\mathrm{GeV^{-1}}$ 。Fermi那篇我们只读了摘要。PDG轴子综述2025汇总了这些结果，并把X射线上限引作 $8\times10^{-13}\,\mathrm{GeV^{-1}}$ 。

预印本。 这里读的文本都是arXiv版本。Reynolds等的第4版页眉仍写着“Submitted to ApJ”；Matthews等的重新分析读的是arXiv第1版，标注已被同一期刊接收。最终期刊版的措辞有没有改动，我们没有比对。

模型预测。 图6里的IAXO线是预期灵敏度，不是测量。Day与Krippendorf 2018年的X射线偏振上限也是：摘要里给的是方法和预期上限，我们只读了摘要。论文里的每条存活曲线，以及Lab里的每条曲线，都是在某个假定磁场下的模型预测。

开放问题。 上限到底多依赖磁场？Matthews等，2022重新分析同一批光谱，发现在约 $10^{-13}$ eV以下，换用别的磁场模型时上限弱0.1到0.3 dex。他们的结论是：只要气体压强与磁压之比约为100，模型B的上限就稳健；但磁场的强度和结构仍不确定。1.2 keV那个凹陷是否纯属定标问题，也还没有定论。

其他途径有什么不同

吸收。 在X射线里，吸收是基线模型的一部分：银河系气体平滑地压低低能端；Matthews等发现，加一个部分覆盖的吸收体并不改变上限。到了伽马射线，角色反过来：甚高能光子会和河外背景光产生电子对而被吸收，一段路以轴子身份走的光子可以躲过去。PDG综述认为这一可能目前尚无定论。

湍流磁场。 每条星系团上限都继承了磁场模型。PDG综述指出，伽马射线上限依赖的磁场，在那些方向上还没有被直接观测证实。

偏振。 只有沿磁场偏振的光会转换，所以偏振带有独立的信息。PDG综述引了来自磁白矮星光学偏振的上限： $1.7\times10^{-12}\,\mathrm{GeV^{-1}}$ ，置信度95%，这些白矮星的磁场可以通过塞曼分裂测出。

恒星冷却。 恒星根本不需要外部磁场。它们通过能量收支限制同一个耦合，对几keV以下的质量都成立，但比Chandra上限宽松约80倍。

下次读一条上限时的检查清单

用的是哪条光谱、哪个波段、哪个基线模型，吸收怎么处理？
用的是哪种磁场模型，有哪些观测约束它？
用了多少个磁场实现，上限是否对它们做了边缘化？
质量范围、置信度和先验分别是什么？
有没有讨论偏好区域，什么仪器效应可能造成它？
你读的是哪个版本的文本，哪些陈述是预期而不是测量？

到Lab里试试：选“百分之一量级的起伏”，比较实线和虚线平均。