光子数分辨能力如何提升多路复用 SNSPD 探测器阵列？

局部分辨能力放宽探测器数量的平方增长

由通断型 SNSPD 构成的多路复用阵列，通过把 N 个光子分散到 D 个像素上来近似实现光子数分辨。由于每个像素只报告有响应或无响应，同一像素内的多光子碰撞会模糊读出结果，而要让这种碰撞概率保持较低所需的元件数量，会随光子数 N 呈平方增长。论文提出的问题是：如果每个单元本身已经携带一定的光子数信息，这一要求能否被放宽。

作者把已有的二进制公式推广到具有任意可分辨态数的探测器。一个 s 态探测器给出 s 个输出类别：0、1、……、s-2 个可分辨光子，以及最后一个饱和类别 (s-1)+，它把占据数等于或大于 s-1 个光子的情况归为一类。论文具体研究了 2 态、3 态和 4 态多路复用系统，其中 3 态探测器区分 0、1 和 2 个及以上光子事件，4 态探测器区分 0、1、2 和 3 个及以上光子事件。

核心结论是：提高局部光子数分辨能力会系统性地改善探测器保真度和贝叶斯光子数重建，并显著放宽纯二进制多路复用阵列所要求的探测器数量与可分辨光子数之间的平方增长关系。这是由响应张量推导出的模型结果，而不是在已制备阵列上的测量结果。

响应张量与二项损耗卷积

核心对象是概率响应张量 P(K(s)|N,D,eta)：当 N 个光子到达由 D 个 s 态探测器组成、量子效率为 eta 的系统时，观测到特定输出向量 K(s) 的概率。输出向量记录有多少个探测器感知到 0、1、2、……(s-1)+ 个光子，其各分量之和为 D。

在理想无损耗情形下，作者用指数生成函数从组合学角度推导该张量，统计 k1 个探测器各接收一个光子、k2 个探测器各接收两个光子，依此类推直到饱和类别的方案数。当 s=2 时，该表达式恰好退化为早先关于把 N 光子 Fock 态均匀分配到 D 个理想通断探测器上的引理。

有限效率通过二项损耗卷积单独处理：由于空间分配是均匀的，含 N 个入射光子的有损系统，其行为等同于注入数量减少为 R 个有效探测光子的理想系统。暗计数被单独评估，作者指出 SNSPD 中典型的亚赫兹暗计数率使其影响相比碰撞驱动和分辨能力驱动的贡献可以忽略，因此数值结果采用无噪声的 4 探测器和 8 探测器系统。

似然矩阵向理想对角线靠拢

为了直接比较读出结果，完整张量被投影到标量估计量 N_s = K(s) . v(s) 上，其中 v(s) = (0, 1, ..., s-1)。这种边缘化得到似然矩阵 P(s)(n_s|N,D,eta)，它覆盖所有给出相同估计量值的输出构型。由于饱和的存在，N_s 是一个下界估计量：一般有 N_s(K) <= N。

在 D=8、eta=0.92 时，论文给出了 N 变化到 16 的似然矩阵。白色对角线标记读出结果等于到达光子数的理想响应。2 态响应被强烈压缩，因为随着 N 增大，被探测到的光子数增长速率下降。3 态似然分布在更接近 N 的 n3 值附近，而 4 态矩阵明显更集中于对角线，因为碰撞和损耗仍然有限。

实际含义是：在探测器数量和效率固定时，随着态数增加，投影后的读出结果对入射光子数更具信息量。这是关于建模似然的陈述，而不是关于某个具体已制备器件的陈述。

保真度提升与 0.55 阈值

论文从响应张量推导出广义 N 光子探测器保真度，并对其进行数值评估。图 2a 报告了在 D=8 时不同效率下保真度随光子数的变化。3 态系统在 N 达到 12 之前给出的保真度明显高于 2 态系统，而无损耗 2 态系统仅在 2 个光子以内优于有损耗 3 态系统。

图 2b 报告了满足最低保真度阈值 0.55 所需的探测器数量。对于效率有限的 2 态系统，所需数量增长快于 e^N，在 7 个光子时已不实用，此时 D 超过 1000。对于 3 态系统，在 5 个光子以内所需数量保持在十个探测器以下，当 eta >= 0.92 时损耗影响很小。

这一比较具体地放宽了早先针对通断探测器的规则，即高保真度要求可分辨特征数量超过 N 的平方。0.55 这一阈值和 D=8 这一设定是这些数字成立的条件，它们并非普适的性能保证。

探测器保真度不等于重建保真度

论文区分了两个容易混淆的指标。N 光子探测保真度问的是探测器输出是否与真实光子数一致。重建保真度问的是从探测器输出中获得了多少关于真实光子数的信息，当光子数估计用于支持后续决策时，它尤其重要。

作者把这一区别具体化：一个与真实光子数一致的饱和结果不会降低探测器保真度，但它产生的后验仍分布在若干光子数上，重建保真度因此下降。任何作用于光子数阈值的协议，只有在给定观测结果时该阈值之上的残余概率质量足够小的情况下才可靠，而重建保真度恰好量化了这一概率。

贝叶斯重建对响应张量求逆，以在给定观测输出向量的情况下推断 N 的似然，并使用一个先验分布。论文还定义了一个系综重建保真度，它对局部保真度在探测器结果和先验分布两者上取平均，并指出先验既隐式地出现在贝叶斯项中，也显式地作为源统计的权重出现。

泊松平均重建受益于高效率

在重建比较中，作者指定入射光子数 N，并用多项分布把它们随机分配到 D 个探测器上；对于等概率探测器，每个探测器接收一个光子的概率为 1/D。他们还表明，低态数系统的输出向量可以通过层级关系从高态数输出向量中还原，例如对于 mu=8、D=6 的泊松源，K(4)=[2,1,2,1] 给出 K(3)=[2,1,3] 和 K(2)=[2,4]。

以探测器效率为函数表示时，泊松平均重建保真度在 eta < 0.9 时小于 0.6，但在 eta 更大时大幅提升，尤其是 D=8 的情形。作者指出，这对多路复用多态 SNSPD 器件尤其有利，因为 SNSPD 可以轻松超过 eta=0.9。先验均值 mu 越小，出于同样的原因效率越高，但 mu=8 的 4 态探测器仍优于 mu=4 的 2 态探测器。

这些是在给定先验和探测器数量下的模型结果。它们描述的是在该框架内重建保真度如何随效率和态数变化，而不是在物理阵列上测得的重建结果。

这些结果的适用范围与边界

该框架在可分辨态数上具有一般性，但数值结果具体针对 2 态、3 态和 4 态多路复用系统，探测器数量为 4 或 8，无噪声运行，效率包括 eta=1 和 eta=0.92。暗计数处理被放在附录中，不属于主要数值结果，理由是亚赫兹的 SNSPD 暗计数相比碰撞和分辨效应可以忽略。

暗计数附录还表明，符合窗口 tau 会强烈改变每个窗口的暗计数概率 epsilon：在 1 Hz 速率下，tau 约为 300 ps 时 epsilon 接近 3E-10，而 tau 约为 100 微秒时 epsilon 接近 1E-4。作者假设 tau 处于 50 至 800 ps 的电阻畴寿命范围内，这正是其暗计数可忽略这一结论成立的条件。

该论文是预印本，其结论是解析和数值性质的。它没有报告已制备的多态阵列，因此保真度和重建数字应被理解为对建模架构的信息价值评估，而不是已演示的器件性能。

优先关注重建保真度而非仅探测器保真度

论文的收尾论点是：在表征光子数分辨探测器时，传统上大量强调 N 光子探测保真度，但重建保真度回答的是一个不同且往往更贴近决策的问题：从输出中实际获得了多少关于真实光子数的信息。一个饱和结果可以与真实光子数一致并保持探测器保真度不变，同时后验仍然分散、重建保真度下降。

因此，对于权衡多态 SNSPD 阵列的读者来说，关键比较是在固定 D 和 eta 下，对 2 态、3 态和 4 态探测器同时按两个指标进行评判。0.55 的保真度阈值，以及 2 态系统在 7 个光子时 D>1000 的要求，标记了二进制多路复用变得不实用的位置；而 3 态系统在 5 个光子以内所需探测器少于十个，当 eta >= 0.92 时损耗影响很小。

最有用的结论取决于具体情景：当光子数估计用于后续阈值决策时，例如评估预报式单光子源是否携带隐藏的多光子污染，重建保真度就是量化残余风险的指标，也是局部光子数分辨能力所改善的指标。

本文引用的文献