跳到主要内容
Sandbox Physics

中微子 · D03 · SNO 与消失的味道

太阳中微子味道观测台

走进重水探测器，回放三种反应，采集模拟曝光，再让各通道的通量约束相交。一个只产生电子味的暗太阳，能解释同一份记录吗？

交互模型太阳中微子味道观测台
推断电子味通量—\text{—}
推断其他活性味通量—\text{—}
推断总通量—\text{—}

三种响应 · 同一轮冻结曝光

消失的味道，留下了指纹。

CC

带电流

只响应电子味

νe+d→p+p+e−\nu_e+d\to p+p+e^-
记录计数—\text{—}
当前检查点预测，含本底—\text{—}
NC

中性流

响应全部活性味

να+d→p+n+να\nu_\alpha+d\to p+n+\nu_\alpha
记录计数—\text{—}
当前检查点预测，含本底—\text{—}
ES

弹性散射

响应加权的味道混合

να+e−→να+e−\nu_\alpha+e^-\to\nu_\alpha+e^-
记录计数—\text{—}
当前检查点预测，含本底—\text{—}

实心棒为记录计数，细棒为当前通量检查点预测；三通道使用同一计数标度。这里理想区分通道，真实 SNO 通过联合统计分析提取相互重叠的事件群。

让通道约束相交

点选通量，或聚焦图表后按方向键。色带显示每个已选通道的一倍计数涨落尺度。

两轴单位均为 106 cm−2 s−110^6\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}

CCNCES最佳拟合
当前通量检查点—\text{—}
响应矩阵、精确计数与光学信号

每单位曝光、通量的计数响应。三行依次为 CC、NC、ES；各行本底另列。这些是教学系数，不是实测 SNO 效率。

A=(3000100100609.24),B=bT(20358)A=\begin{pmatrix}300&0\\100&100\\60&9.24\end{pmatrix},\qquad B=bT\begin{pmatrix}20\\35\\8\end{pmatrix}

对照历史实验结果

SNO 2002 年论文在指定硼-8 能谱假设下报告电子味与总活性通量。以下发表值是独立参考；本模拟器不拟合它们，也不复现合作组分析。

Φe=1.76±0.05 (stat.)±0.09 (syst.)ΦNC=5.09−0.43+0.44 (stat.)−0.43+0.46 (syst.)\begin{aligned}\Phi_e&=1.76\pm0.05\,\mathrm{(stat.)}\pm0.09\,\mathrm{(syst.)}\\\Phi_{\mathrm{NC}}&=5.09^{+0.44}_{-0.43}\,\mathrm{(stat.)}^{+0.46}_{-0.43}\,\mathrm{(syst.)}\end{aligned}

106 cm−2 s−110^6\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}

SNO · 2002 年原始论文 ↗

继续探索

从味道到证据

传播预测味概率，实验还需要哪些环节，才能把相互作用变成可用记录？

物理教程

同一个来源，三种不同的问题。

背景放射化学实验发现太阳电子中微子偏少。SNO 用重水多问了一句：不管味道，总共有多少活性中微子到达？比较带电流、中性流和弹性散射响应，就能分开较暗的来源与味道改变。

为什么值得理解本实验把推断直接画出来。每个通道在二维通量平面上约束一条带；单一电子味敏感的速率无法分开源变弱与变味，不同敏感性才能补足信息。

先抓重点

核心问题
相同的电子味通道计数，能否来自总活性通量不同的两个太阳？
一句话直觉
中性流把活性味一起计数；电子味速率低，不一定意味着总活性通量低。

核心数学模型

两个通量分量

Φx=Φμ+Φτ,Φ=Φe+Φx\Phi_x=\Phi_\mu+\Phi_\tau,\qquad\Phi=\Phi_e+\Phi_x

第二个分量合并其他两种活性味。本模型不能分开 μ 子味与 τ 子味，也不包含惰性态。

三种理想通道响应

λCC=300TΦe+20bTλNC=100T(Φe+Φx)+35bTλES=60T(Φe+0.154Φx)+8bT\begin{aligned}\lambda_{\mathrm{CC}}&=300T\Phi_e+20bT\\\lambda_{\mathrm{NC}}&=100T(\Phi_e+\Phi_x)+35bT\\\lambda_{\mathrm{ES}}&=60T(\Phi_e+0.154\Phi_x)+8bT\end{aligned}

通量数值以每平方厘米每秒百万个为单位，曝光为相对量。系数与已知本底都是教学选择；固定散射权重用于说明依赖能谱的响应，不是普适比值。

冻结计数记录

Ni∼Poisson⁡(λi)N_i\sim\operatorname{Poisson}(\lambda_i)

三个理想区分的通道独立采样。真实 SNO 拟合重叠事件群并处理相关误差；光学示例不用来给这里的总计数标记通道。

只拟合选中的响应

Q=∑i∈S[Ni−Bi−T(AΦ)i]2max⁡(Ni,1)Q=\sum_{i\in\mathcal S}\frac{[N_i-B_i-T(A\boldsymbol\Phi)_i]^2}{\max(N_i,1)}

无约束加权最小二乘用观测计数估计方差，下限为一。局部协方差与残差分数只描述该近似，都不是标定的置信区域或显著性。

光子到达显示的传感器球

tj=n ∣rj−rv∣ct_j=\frac{n\,|\mathbf r_j-\mathbf r_v|}{c}

场景射线都终止于已存的传感器坐标，距球心 9 米。时间从有效电子发光开始，水折射率取 1.33；省略亚克力折射、色散、散射与中子俘获延迟。

常见难点

中性流闪光不是中微子的直接图像

容易误解NC 场景中的光是中子发出的。

正确理解释放的中子被俘获，产生的伽马射线散射电子，再由电子发出切伦科夫光。本方向平均示意从有效发光点开始，省略俘获延迟与扩散。

负拟合值是近似的警告

容易误解其他味估计为负，意味着物理通量为负。

正确理解这里特意采用无约束拟合。接近零时，计数涨落会使估计超出物理区域；实际推断需要物理边界处理，本实验不提供该区间分析。

动手实验

  1. 01

    比较三种反应

    在探测器上方选择 CC、NC 和 ES，拖动光飞行时间轴，检查侧视与相互作用特写。

    应该观察到： CC 和 ES 显示展宽的电子光锥；NC 显示方向平均的俘获闪光。中子不带电，不是发出切伦科夫光的粒子。
  2. 02

    让两个太阳给出相近的电子味计数

    先采集暗太阳，再采集已经变味，用公共标度比较 CC 与 NC。

    应该观察到： 两个预设的电子味分量相同，但总活性通量较大时 NC 增多；电子味减少不等于把来源整体调暗。
  3. 03

    移除信息

    先关掉 NC，再关掉 ES，只留下 CC 拟合，随后恢复第二通道。

    应该观察到： 单一速率只能约束一条线，不能确定唯一交点。在指定响应下，CC 加 ES 可分辨两分量，但其他味误差比 CC 加 NC 大。
  4. 04

    检查纯电子味解释

    检查最好的纯电子味假设，比较其预测棒与记录计数。

    应该观察到： 该假设在选中通道上单独优化。记录需要其他活性味分量时，满足单一通道无法同时满足另外的通道。
  5. 05

    分开真值和推断

    揭示模拟来源，改变随机种子，再采集长曝光并导出 CSV。

    应该观察到： 真值标记独立于拟合。计数会涨落，最佳拟合未必正好命中来源；显示的计数、响应、示意信号与传感器坐标都可导出。