---
title: "读懂Chandra对NGC 1275的轴子上限"
description: "逐图读Reynolds等人2020年的论文：光谱显示了什么，磁场模型怎样决定上限，伽马射线、偏振和恒星又有何不同。"
canonical_url: "https://sandboxphysics.com/zh/stories/axion-photon-mixing/reading-the-chandra-ngc-1275-limit/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/axion-photon-mixing/reading-the-chandra-ngc-1275-limit/"
last_modified: "2026-11-12"
---

[用宇宙磁场找轴子](https://sandboxphysics.com/zh/stories/axion-photon-mixing/)

# 读懂Chandra对NGC 1275的轴子上限

Sandbox Physics · 2026-11-12 · L4 · 论文阅读

把Christopher Reynolds等人的arXiv:1907.05475v4打开放在旁边，这篇论文后来于2020年发表在《天体物理学报》（ApJ）上。整个论证落在图3到图6和式(2)到式(5)上：一条平滑的X射线谱，一个英仙座磁场模型，一个存活曲线库，以及质量和耦合上的后验分布。读完之后，你应该能说清$g_{a\gamma}<7.9\times10^{-13}\,\mathrm{GeV^{-1}}$这条上限靠什么撑着，以及伽马射线、偏振和恒星为什么以不同方式限制同一个耦合。[Reynolds等，2020](https://arxiv.org/abs/1907.05475v4)

我们通读的是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微高斯，按等温星系团处理时随密度的平方根下降。

论文明说，无论小块模型还是高斯随机场，都不是自然界真实存在的磁场；对转换要紧的是横向强度和沿视线的相干长度。本系列的[星系团磁场与等离子体中的混合](https://sandboxphysics.com/zh/stories/axion-photon-mixing/mixing-in-cluster-fields-and-plasma/)沿模型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年](https://arxiv.org/abs/1406.6053v2)的水平分支上限$g_{a\gamma}<6.6\times10^{-11}\,\mathrm{GeV^{-1}}$，以及[Ajello等2016年](https://arxiv.org/abs/1603.06978v1)用同一星系六年数据得到的Fermi-LAT上限：0.5到5 neV之间$5\times10^{-12}\,\mathrm{GeV^{-1}}$。Fermi那篇我们只读了摘要。[PDG轴子综述2025](https://pdg.lbl.gov/2025/reviews/rpp2025-rev-axions.pdf)汇总了这些结果，并把X射线上限引作$8\times10^{-13}\,\mathrm{GeV^{-1}}$。

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

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

**开放问题。** 上限到底多依赖磁场？[Matthews等，2022](https://arxiv.org/abs/2202.08875v1)重新分析同一批光谱，发现在约$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倍。

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

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

[到Lab里试试](https://sandboxphysics.com/zh/labs/axion-photon-conversion/)：选“百分之一量级的起伏”，比较实线和虚线平均。

## 来源

1. [Reynolds等，2020](https://arxiv.org/abs/1907.05475v4)
2. [Ayala等2014年](https://arxiv.org/abs/1406.6053v2)
3. [Ajello等2016年](https://arxiv.org/abs/1603.06978v1)
4. [PDG轴子综述2025](https://pdg.lbl.gov/2025/reviews/rpp2025-rev-axions.pdf)
5. [Day与Krippendorf 2018年](https://arxiv.org/abs/1801.10557v2)
6. [Matthews等，2022](https://arxiv.org/abs/2202.08875v1)

## 继续阅读

用宇宙磁场找轴子 · 阅读路径

### 从两个混合态读到一条已发表上限

1. [文章 **看不见的轻粒子，怎样让来自远方的光谱出现变化？**](https://sandboxphysics.com/zh/stories/axion-photon-mixing/invisible-particles-and-distant-spectra/)
2. [文章 **把光子和轴子看成两个会混合的态**](https://sandboxphysics.com/zh/stories/axion-photon-mixing/photon-and-axion-as-two-states/)
3. [文章 **星系团磁场与等离子体中的混合**](https://sandboxphysics.com/zh/stories/axion-photon-mixing/mixing-in-cluster-fields-and-plasma/)
4. [文章 **随机磁畴、偏振与恒星冷却**](https://sandboxphysics.com/zh/stories/axion-photon-mixing/random-domains-polarization-and-stars/)
5. [文章 **用保存的输入复算一条星系团视线**](https://sandboxphysics.com/zh/stories/axion-photon-mixing/recompute-a-cluster-sightline/)
6. [当前页面 **读懂Chandra对NGC 1275的轴子上限**](https://sandboxphysics.com/zh/stories/axion-photon-mixing/reading-the-chandra-ngc-1275-limit/)