星系团磁场中的轴子-光子转换

看不见的轻粒子，怎样让来自远方的光谱出现变化？ 磁场能让光子变成很轻的轴子，再变回来，所以远方星系的一部分光会在路上消失。消失的份额随光子能量变化，会在X射线谱上留下细小起伏。观测到的光谱足够平滑，就能限制轴子与光耦合的强度；恒星冷却给出另一条独立限制。

转动磁畴，改变气体密度和轴子质量，看光子在星系团视线上和轴子来回转换。

到达的光子 $88.3\,\%$

变成轴子的光子 $11.7\,\%$

到达时的线偏振度 $11.6\,\%$

氦燃烧寿命剩余 $100\,\%$

让X射线穿过磁化的星系团

光子能量

3\,\mathrm{keV}

E

轴子-光子耦合

3.16\times10^{-12}\,\mathrm{GeV^{-1}}

g_{a\gamma}

垂直于光束的磁场

5\,\mu\mathrm{G}

B_\perp

气体电子密度

0.00316\,\mathrm{cm^{-3}}

n_e

轴子质量

1\times10^{-13}\,\mathrm{eV}

m_a

在场景里拖动磁畴箭头，可改变磁场方向和强度。在光谱图上横向拖动，可选择上方显示的能量。

这个耦合低于该质量下适用的全部已发表上限。

磁畴、模型各项与检查

磁畴数目

24

N

每个磁畴的长度

8\,\mathrm{kpc}

\ell_B

所选能量下第一个磁畴的各项，单位为千秒差距的倒数：

$\Delta_{a\gamma}$

$\omega_{\rm pl}$

高于临界能量时，混合不再依赖能量；远低于它时，转换被压低。

$E_{\rm crit}=\frac{|m_a^2-\omega_{\rm pl}^2|}{2g_{a\gamma}B_\perp}$

$\ell_{\rm osc}=2\pi/\Delta_{\rm osc}$

每个磁畴使用精确的幺正传递矩阵，总概率除舍入误差外保持为一：

$|P_{\gamma\gamma}+P_{\gamma\to a}-1|$

导出文件保存输入、单位、磁场方向、两条曲线及其标签，以及沿程曲线。

模型假设

这是单条视线的简化模型，没有拟合任何观测。星系团磁场被写成一串等长磁畴；每个磁畴内，垂直于光束的磁场均匀，强度相同，方向随机或一致。气体电子密度取常数。每个磁畴内，光子与一种类轴子粒子的混合方程被精确求解。模型略去了光子吸收、真实磁场的湍流功率谱、密度与磁场的径向分布、宇宙背景光和探测器响应。虚线是理想参照：均匀磁场的两态公式，或对随机磁场方向求平均的结果，不是从实线反推出来的。恒星冷却标尺用的是粒子数据组综述里水平分支星的另一条方程。编号上限照录所引论文的数值，只在各自的质量范围和假设下成立。移动的示踪点只用来示意计算出的比例。

穿过星系团后留下的X射线光子，按能量排列

实线：这条模拟视线虚线：随机磁场方向的理想平均

纵轴是到达光子的比例。两条曲线都来自模型，不是望远镜数据。

所选能量下沿视线的变化

P_{\gamma\gamma}

P_{\gamma\to a}

\Pi_L

点线是磁畴边界。每个磁畴只有一个光子偏振方向参与混合，所以线偏振度会增长。

恒星给出另一条限制

氦燃烧寿命因子

水平分支星计数 $g_{a\gamma}<6.6\times10^{-11}\,\mathrm{GeV^{-1}}$ $m_a\lesssim1000\,\mathrm{eV}$ 置信度95%；PDG 2025 review 89, eq. (89.16); Ayala et al. 2014, arXiv:1406.6053v2, eq. (3)
CAST太阳轴子望远镜 $g_{a\gamma}<5.8\times10^{-11}\,\mathrm{GeV^{-1}}$ $m_a\lesssim0.02\,\mathrm{eV}$ 置信度95%；PDG 2025 review 89, eq. (89.20)
NGC 1275的Fermi-LAT伽马谱 $g_{a\gamma}<5\times10^{-12}\,\mathrm{GeV^{-1}}$ $5\times10^{-10}\lesssim m_a/\mathrm{eV}\lesssim5\times10^{-9}$ 置信度95%；Ajello et al. 2016, arXiv:1603.06978v1, abstract
NGC 1275的Chandra X射线谱 $g_{a\gamma}<7.9\times10^{-13}\,\mathrm{GeV^{-1}}$ $m_a\lesssim1\times10^{-12}\,\mathrm{eV}$ 置信度99.7%；Reynolds et al. 2020, arXiv:1907.05475v4, eq. (5), field model B

灰色虚线标记在当前质量下不适用。上限照录论文数值，不在此计算。

核心问题: 同样的耦合，为什么有时给出平滑光谱，有时给出明显起伏？
一句话直觉: 只有当耦合项压过光子与轴子色散项之差时，混合才会强。磁畴方向随机时，光子比例的平均值趋向 $2/3$ ，而任何一条具体视线都会在这个平均值附近保留随能量变化的起伏。

核心数学模型

一个磁畴里的混合

i\frac{d}{dz}\begin{pmatrix}A_\perp\\A_\parallel\\a\end{pmatrix}=-\begin{pmatrix}\Delta_{\rm pl}&0&0\\0&\Delta_{\rm pl}&\Delta_{a\gamma}\\0&\Delta_{a\gamma}&\Delta_a\end{pmatrix}\begin{pmatrix}A_\perp\\A_\parallel\\a\end{pmatrix}

只有平行于横向磁场的光子偏振分量参与混合。其中 $\Delta_{a\gamma}=g_{a\gamma}B_\perp/2$ ， $\Delta_{\rm pl}=-\omega_{\rm pl}^2/(2E)$ ， $\Delta_a=-m_a^2/(2E)$ 。

均匀磁场中的两态转换

P_{\gamma\to a}=\left(\frac{2\Delta_{a\gamma}}{\Delta_{\rm osc}}\right)^2\sin^2\!\left(\frac{\Delta_{\rm osc}L}{2}\right),\qquad \Delta_{\rm osc}^2=(\Delta_{\rm pl}-\Delta_a)^2+4\Delta_{a\gamma}^2

这是沿磁场偏振的光子。非偏振光的转换概率减半。Lab逐个磁畴精确求解，并在均匀磁场下与此公式核对。

对随机磁场方向求平均

\langle P_{\gamma\to a}\rangle_N=\frac{1}{3}\left[1-\prod_{j=1}^{N}\left(1-\frac{3}{2}P_j\right)\right]

每个 $P_j$ 是第 j 个磁畴对平行偏振的转换概率。磁畴足够多时，平均值趋向三分之一，即两种光子偏振态与一种轴子态均分。

另一条恒星方程

\frac{t_{\rm HB}(g)}{t_{\rm HB}(0)}\approx\frac{80}{80+30\,g_{10}^2},\qquad g_{10}=\frac{g_{a\gamma}}{10^{-10}\,\mathrm{GeV^{-1}}}

氦燃烧核心里的普里马科夫辐射与核燃烧竞争。如果球状星团中的水平分支星比预期少，就说明有轴子冷却。这条方程不使用视线模型。

常见难点

磁场更强不一定转换更多

容易误解磁场加倍，转换份额就加倍。

正确理解概率随距离振荡。一个磁畴里超过半个振荡周期后，在给定能量处加强磁场反而可能减少光子损失。

虚线不是拟合曲线

容易误解虚线是从模拟曲线反推出来的。

正确理解它是独立的理想公式。在缠绕磁场里，它是对所有磁场方向的平均，具体视线会在它附近散开。

动手实验

01
先看默认视线
让缠绕磁场示例运行，看示踪点在两条道之间漂移。
你会看到： 光子道一格一格地变稀又回升，下方的光谱随能量起伏不齐。
02
转动一个磁畴
把链中间一个磁畴的箭头拖转半圈。
你会看到： 沿程曲线从这个磁畴开始改变，望远镜处的偏振方向也跟着转。
03
让磁场变均匀
选“均匀磁场里的振荡”，再拖动能量竖线。
你会看到： 实线与虚线重合，因为精确传递矩阵复现了两态公式。
04
再和恒星比较
把耦合加大到恒星标尺上的第一条编号上限。
你会看到： 这时氦燃烧缩短约百分之十四，而星系团视线已接近完全混合。在各自的磁场假设下，X射线谱能探到比恒星计数弱得多的耦合。

星系团磁场中的轴子-光子转换

跟着X射线光子穿过磁化的星系团

核心数学模型

一个磁畴里的混合

均匀磁场中的两态转换

对随机磁场方向求平均

另一条恒星方程

常见难点

磁场更强不一定转换更多

虚线不是拟合曲线

动手实验

先看默认视线

转动一个磁畴

让磁场变均匀

再和恒星比较

相关实验与文章

文章

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