Implications of \textit{SARAS3} data for Coulomb-like interacting dark matter

本文通过分析SARAS3在55.5–84.4 MHz频段未探测到21厘米信号的结果，通过自洽地建模气体冷却与结构形成抑制，对类库仑相互作用的暗物质施加约束，最终发现相互作用暗物质相较于标准冷暗物质并无统计显著的偏好性，同时确立了全球21厘米信号振幅的有意义上限。

要点

以下是该论文的通俗化解释，借助类比帮助理解其中的概念。

宏观图景：聆听宇宙的“婴儿啼哭”

想象早期的宇宙是一个巨大而黑暗的育婴室。在大爆炸后约 1 亿到 2 亿年，第一代恒星刚刚开始诞生。这些恒星发出的光与充满宇宙的氢气发生相互作用，产生了一种特定的无线电信号，即21 厘米信号。

可以将这种信号比作宇宙黎明时的“婴儿啼哭”。如果我们能清晰地听到它，就能得知气体的温度是冷是热，以及第一代恒星形成的速度有多快。

长期以来，科学家们希望能听到这种啼哭。然而，该信号极其微弱，就像试图在飓风中听到一声低语。这场“飓风”由我们银河系、地球大气层以及射电望远镜本身产生的无线电噪声构成。

谜团：暗物质的秘密对话

我们知道，宇宙的大部分由暗物质组成，这是一种不发光的神秘物质。标准理论认为暗物质是“冷”且“懒惰”的——它只是静止不动，仅通过引力与普通物质（如气体）发生相互作用。

但如果暗物质更像一只“社交蝴蝶”呢？如果它会与普通气体粒子碰撞并交换热量，就像两个人握手并传递体温一样？这就是**相互作用暗物质（IDM）**的概念。

本文作者希望测试一种特定的“社交型”暗物质，其相互作用方式类似于库仑力（类似于电荷之间的吸引或排斥）。他们问道：如果暗物质真的如此，它将如何改变“婴儿啼哭”（即 21 厘米信号）？

两步效应：冷却与延迟

本文指出，如果暗物质与气体相互作用，会导致两个主要变化，作者对此进行了仔细建模：

1. “冰袋”效应（冷却）：
   通常情况下，随着宇宙膨胀，气体会缓慢冷却。但如果暗物质比气体更冷，它就会像冰袋一样，从气体中吸走热量。这使得气体变得比预期更冷。

   • 结果： 更冷的气体会产生更深、更响亮的“啼哭”（即更强的吸收信号）。

2. “交通堵塞”效应（恒星延迟）：
   当暗物质与气体碰撞时，会产生摩擦（阻力）。这会减缓气体的运动，使其更难坍缩并形成恒星。

   • 结果： 恒星形成被推迟。由于恒星提供热量和光线，最终会使气体变暖，因此“啼哭”发生的时间会延后，且强度会比恒星按时形成时更弱。

作者意识到，以往的研究往往只关注“冰袋”（冷却）效应，而忽略了“交通堵塞”（恒星延迟）效应。本文首次同时建模了这两种效应，以呈现完整的图景。

侦探工作：SARAS3 实验

为了验证这一理论，研究团队查看了SARAS3 实验的数据。

• 设置： 与地面上的其他望远镜不同，SARAS3 是漂浮在湖面上的天线。水面充当了完美且均匀的背景，有助于过滤掉来自地面的部分“噪声”。
• 结果： SARAS3 在特定频率范围内寻找“婴儿啼哭”，但并未发现它。他们只看到了静电噪声。

调查：“一无所获”揭示了什么？

通常，当科学家说“我们没找到它”时，这感觉像是一个死胡同。但作者将这种“零结果”（一无所获）视为一条线索。

他们构建了一个复杂的计算机模型，模拟了：

1. 基于其暗物质理论的“婴儿啼哭”（21 厘米信号）。
2. “噪声”（如银河系无线电波等前景干扰）。

随后，他们使用统计方法（贝叶斯推断）来观察他们的“暗物质 + 噪声”模型是否能解释 SARAS3 的数据。

发现：

• 信号被掩盖： 数据过于嘈杂，无法确定暗物质的确切质量或其相互作用的强度。这就像站在狂风中试图猜测一根羽毛的确切重量；风（噪声）太大，无法分辨。
• “太响亮”的规则： 然而，他们可以说明信号不是什么。数据证明，在他们观测的频率范围内，“婴儿啼哭”不可能极其深或响亮。具体来说，在某个特定时间点（红移 23.6），信号深度不能超过**-277.6 毫开尔文**。如果暗物质的相互作用强到足以使信号达到这个深度，SARAS3 本应观测到它。既然他们没有观测到，那些特定的强相互作用就被排除了。
• 暗物质与标准模型： 作者将他们的“社交型暗物质”模型与标准的“懒惰型暗物质”模型进行了比较。他们问道：数据是否更倾向于社交版本？
  • 裁决： 不。数据尚无定论。这是一场平局。虽然“社交版本”的胜算略高（1.7 比 1），但不足以断定其绝对正确。本质上，这是一场平局。

结论

本文是一堂关于如何聆听寂静的课程。尽管 SARAS3 没有发现信号，但作者了解到：

1. 我们尚不能排除暗物质与气体相互作用的想法，但我们知道它不能过于强烈（否则信号会太响亮而无法被忽略）。
2. 要解开这个谜团，我们需要来自REACH等未来实验的更好数据（更小的风噪、更清晰的信号）。
3. “社交型暗物质”理论依然成立，但尚未被证实。

简而言之：宇宙仍在低语，而我们仍在试图弄清楚这低语是来自一个标准的幽灵，还是一个健谈的幽灵。SARAS3 告诉我们这个幽灵没有大喊大叫，但它尚未告诉我们它究竟在低语什么。

技术摘要

以下是论文《SARAS3 数据对类库仑相互作用暗物质的启示》的详细技术总结。

1. 问题陈述

暗物质（DM）的本质仍是宇宙学中最重要的未解之谜之一。虽然标准冷暗物质（CDM）模型假设暗物质是无碰撞的，但各种异常现象（例如矮星系动力学、宇宙射线正电子过剩）暗示了暗物质与重子之间可能存在非引力相互作用。具体而言，**类库仑相互作用暗物质（IDM）**模型提出，暗物质粒子携带微小的分数电荷，或通过一种与相对速度四次方成反比的力进行相互作用（$\sigma \propto v^{-4}$）。

宇宙黎明时期的 21 厘米信号（即第一代恒星形成的时期）是探测星系际介质（IGM）热历史的高灵敏度探针。先前的研究通常依赖具有争议的 EDGES 探测结果，认为 IDM 可以通过过量的重子冷却来解释异常深的吸收谷。然而，SARAS3 实验（背景射电频谱整形天线测量）在 55.5–84.4 MHz 频段未探测到如此深的特征，以 95% 的置信度排除了 EDGES 的轮廓。

核心挑战： 先前利用 21 厘米数据对 IDM 的约束往往不完整，主要关注重子冷却效应，而忽略了动量交换效应。动量交换会抑制结构形成（晕的形成），从而延迟恒星形成，并减少 Ly$\alpha$、X 射线和电离背景。本文旨在利用 SARAS3 数据提供首个自洽的 IDM 分析，在严格的贝叶斯框架内同时模拟冷却效应和延迟恒星形成效应。

2. 方法论

A. 理论框架

作者利用 ECHO21 Python 包对全球 21 厘米信号进行建模，纳入了 IDM 的两种不同物理机制：

1. 过量重子冷却： 暗物质与重子之间的热交换使气体温度（$T_k$）低于绝热极限，从而加深 21 厘米吸收谷。
2. 延迟恒星形成： 动量交换抑制物质功率谱（$P(k)$），降低晕质量函数（HMF）。这会延迟恒星形成的开始，进而延迟 Ly$\alpha$ 耦合和 X 射线加热。这种效应倾向于将吸收谷移至更低的红移并减小其深度。

该模型求解以下耦合微分方程：

• 气体动力学温度（$T_k$）
• 暗物质温度（$T_\chi$）
• 暗物质与重子之间的相对速度（$v_{b\chi}$）
• 电离分数（$x_e$）

B. 参数空间

该模型由 7 个自由参数定义：

• IDM 物理： 暗物质粒子质量（$m_\chi$）和相互作用截面参数（$\sigma_{45}$）。
• 天体物理： Ly$\alpha$ 背景缩放（$f_{Ly}$）、X 射线背景缩放（$f_X$）、X 射线谱指数（$w$）、电离光子逃逸分数（$f_{esc}$）以及恒星形成的最小维里温度（$T_{vir}^{min}$）。
• 注： 由于与其他参数存在简并性，恒星形成效率（$f_*$）被固定为 0.1。

C. 数据与推断

• 数据集： SARAS3 天线温度数据（55.5–84.4 MHz）。
• 前景建模： 使用 6 阶对数 - 对数多项式（7 个系数）来模拟银河系/河外前景和系统误差。
• 统计方法： 使用 Polychord 嵌套采样算法进行联合贝叶斯推断。
  • 总模型为 14 维（7 个信号参数 + 7 个前景系数）。
  • 先验： 所有参数均使用均匀（无信息）先验以避免偏差。
  • 似然函数： 基于观测天线温度与建模天线温度之差的 Gaussian 似然函数，并按 SARAS3 噪声水平加权。
• 分析指标：
  • 贝叶斯证据： 用于比较 IDM 与 CDM 模型（通过贝叶斯因子）。
  • Kullback-Leibler (KL) 散度： 用于衡量从先验到后验的信息增益（约束能力）。
  • 高斯模型维度（GMD）： 用于量化有效约束参数的数量。

3. 主要贡献

1. 自洽的 IDM 建模： 这是第一项通过同时包含冷却效应（加深信号）和动量转移效应（抑制恒星形成并减弱信号）来联合约束 IDM 参数的研究。
2. 严格的贝叶斯框架： 与先前近似分析不同，本研究使用嵌套采样对前景和天体物理不确定性进行了完全边缘化处理。
3. 函数约束： 作者不仅约束了单个参数，还将参数后验分布转化为红移范围内 21 厘米信号幅度（$T_{21}$）的函数约束。
4. 模型比较： 利用 SARAS3 数据对 IDM 情景与标准 CDM 情景进行了直接的贝叶斯证据比较。

4. 结果

A. 参数约束

• 参数约束微弱： 鉴于 SARAS3 数据的噪声水平以及前景与宇宙学信号之间的简并性，数据并未显著约束单个 IDM 或天体物理参数。
  • 7 维信号参数空间的 KL 散度极低（$D_{KL} \approx 0.007$）。
  • 高斯模型维度可忽略不计（$d_G \approx 0.01$）。
  • $m_\chi$、$\sigma_{45}$ 和天体物理参数的后验分布与先验分布基本一致。

B. 信号幅度界限

虽然单个参数未受约束，但数据确实排除了观测频段内的深吸收特征：

• 最强的约束出现在 $z = 23.6$（$\nu \approx 57.7$ MHz）处，此处 KL 散度达到峰值（$0.52$）。
• 在该红移处，数据对信号幅度施加了 $3\sigma$ 下限：
  $$T_{21}(z=23.6) \gtrsim -277.6 \text{ mK}$$
• 这有效地排除了在 SARAS3 频段内预测吸收特征深于该限制的 IDM 模型。

C. IDM 与 CDM 的偏好

• 作者计算了 IDM 和标准 CDM 模型的贝叶斯证据。
• 对数贝叶斯因子： $\ln B = \ln Z_{IDM} - \ln Z_{CDM} \approx 0.53$。
• 贝叶斯因子： $B \approx 1.7$。
• 结论： 根据 Jeffreys 标度，这代表非决定性证据。基于当前的 SARAS3 数据，在统计上没有显著偏好 IDM 胜过 CDM。

5. 意义与启示

• 零结果具有信息量： 尽管 SARAS3 未探测到信号，但这一非探测结果为类库仑 IDM 模型的全局 21 厘米信号幅度提供了有意义的上限。它排除了会产生深吸收谷的最极端冷却情景。
• 动量交换的重要性： 该研究强调，忽略动量交换（延迟恒星形成）会导致高估冷却效应。自洽模型表明，动量抑制可以抵消冷却作用，使信号变浅，更难与 CDM 区分。
• 未来前景： 目前的约束受限于 SARAS3 的噪声水平和高红移天体物理（如 X 射线和 Ly$\alpha$ 背景）的巨大不确定性。作者建议，未来的实验需要更宽的带宽（例如 REACH），并结合数据集（例如 21 厘米功率谱、紫外光度函数），以打破简并性并对 IDM 参数施加更严格的约束。
• 工具发布： 作者已更新并发布 ECHO21 Python 包，纳入了类库仑 IDM 建模，以促进该领域的未来研究。

总之，本文建立了一个稳健的框架，用于利用全局 21 厘米数据检验相互作用暗物质，证明了虽然当前数据无法区分 CDM 和 IDM，但它成功排除了预测极深吸收特征的模型。