2、单光子探测基础：光子统计特性、单光子探测器类型（SPAD、SNSPD、PMT）及工作原理

单光子探测是量子光学、量子密钥分发（QKD）、激光雷达（LiDAR）、生物荧光寿命成像（FLIM）以及暗物质搜索等前沿领域的核心使能技术。在宏观尺度下，光表现为连续的电磁波；但在单光子层面，光呈现离散的粒子性，其能量量子为 \(E = h\nu\)。要实现对单个光子的有效探测，必须深刻理解光子的统计特性，并依据不同应用场景选择合适的光电探测器件。本章将系统梳理光子统计的数学物理基础，深入剖析SPAD、SNSPD、PMT三类主流单光子探测器的工作原理、关键性能参数与选型依据，并结合实际工程案例给出噪声抑制与信号增强的底层逻辑。

2.1 光子统计特性：从相干光到非经典光

光子统计特性描述了光场中光子数涨落的概率分布规律，是设计单光子探测系统时选择探测器类型、评估噪声基底和优化信号提取算法的理论基础。根据光场的量子态不同，光子数分布呈现三种典型统计行为：泊松分布、超泊松分布和亚泊松分布。

2.1.1 核心原理：三种统计分布与二阶关联函数

对于理想的相干光源（如单纵模激光器），光子数 \(n\) 在探测时间窗口内的概率分布服从泊松分布：

P(n) = (μ^n · e^(-μ)) / n!

其中 μ 为平均光子数。泊松分布的方差等于均值（σ² = μ），其归一化二阶关联函数 g²(0) = 1，表明光子到达事件彼此独立，无时间关联性。对于热光源（如卤素灯、超辐射发光二极管），光子倾向于聚束到达，呈现超泊松分布，方差大于均值，g²(0) > 1。而对于单光子源（如量子点、色心）产生的Fock态，光子数分布呈现亚泊松特性，方差小于均值，g²(0) < 1，这是非经典光场的标志性特征。

在实际单光子探测系统中，我们关注的不仅是信号光子的统计特性，还包括探测器自身的噪声统计。暗计数（无光入射时的计数）通常服从泊松分布，而后脉冲（由前一次雪崩或倍增过程引起的虚假计数）则具有时间相关性，会破坏计数的独立性假设。理解这些统计特性是后续设计噪声抑制算法（如符合计数、时间相关单光子计数TCSPC）的前提。

2.1.2 判断方法：如何评估光子统计特性

在工程实践中，评估光子统计特性的标准方法是测量二阶关联函数 g²(0)。使用Hanbury Brown-Twiss（HBT）干涉仪，将光束分束后分别送入两个单光子探测器，测量符合计数率与单通道计数率的关系：

g²(0) = N_coincidence / (N_1 · N_2 · τ_window)

其中 N_coincidence 为两通道符合计数率，N₁ 和 N₂ 分别为两通道单计数率，τ_window 为符合时间窗口。当 g²(0) < 0.5 时，可判定为单光子源；g²(0) = 1 对应相干光；g²(0) > 1 则表明存在多光子聚束或背景噪声贡献。

2.1.3 案例拆解：QKD系统中光子统计对误码率的影响

在量子密钥分发（QKD）系统中，发送端通常使用弱相干激光脉冲（平均光子数 μ ≈ 0.1-0.5）近似单光子源。以BB84协议为例，若 μ 设置过高（如 μ = 1.0），则多光子脉冲概率显著增加，攻击者可通过光子数分裂（PNS）攻击窃取信息；若 μ 设置过低（如 μ = 0.01），则信号计数率急剧下降，暗计数占比上升，量子误码率（QBER）恶化。实际系统中，需根据信道损耗和探测器暗计数率优化 μ 值。例如，在15 dB信道损耗下，采用暗计数率 100 Hz 的SPAD探测器，最优 μ 约为 0.3-0.4，此时QBER可控制在 3% 以下。

2.1.4 常见误区与总结

误区一：认为所有单光子探测器都适用于任何波长。实际上，SPAD在近红外（>1000 nm）探测效率急剧下降，而SNSPD在可见光波段虽然也能工作，但性价比不如SPAD。误区二：忽视后脉冲效应。后脉冲会导致计数率偏高，尤其在高速计数场景下，需通过设置合适的死时间或后脉冲抑制电路来消除。误区三：将泊松分布直接套用于所有光源。对于脉冲激光器，其光子数分布可能偏离理想泊松模型，需实测验证。

2.2 单光子探测器类型与工作原理

目前工程上最常用的三类单光子探测器分别为：单光子雪崩二极管（SPAD）、超导纳米线单光子探测器（SNSPD）和光电倍增管（PMT）。三者的物理机制、性能边界和适用场景差异显著，选型时需综合权衡探测效率、暗计数、时间分辨率、死时间、工作温度和成本。

2.2.1 核心原理：三类探测器的物理机制

SPAD（单光子雪崩二极管）：基于PN结的盖革模式雪崩倍增效应。当反向偏压超过击穿电压（V_BR）时，单个光生载流子触发自持雪崩，产生宏观可测电流脉冲。雪崩过程需通过淬灭电路（被动淬灭电阻或主动淬灭晶体管）终止，随后进入死时间恢复期。SPAD的增益可达 10⁶ 以上，时间抖动（Timing Jitter）典型值为 50-300 ps，探测效率在可见光波段可达 50-70%。

SNSPD（超导纳米线单光子探测器）：利用超导材料（如NbN、WSi）在临界温度以下处于零电阻态的特性。当光子入射到处于电流偏置的纳米线上时，局部超导态被破坏形成热点（Hotspot），产生电阻区域，从而在偏置电路中形成电压脉冲。SNSPD具有极低的暗计数（<1 Hz）、皮秒级时间抖动（<30 ps）和宽波段响应（可见光至中红外），但需要工作在 2-4 K 的低温环境，系统复杂度高。

PMT（光电倍增管）：基于外光电效应和二次电子倍增。光子入射到光电阴极产生光电子，经聚焦电极加速后撞击倍增极（Dynode），逐级倍增产生大量二次电子，最终在阳极形成电流脉冲。PMT的增益可达 10⁵-10⁷，响应速度快（上升时间 <1 ns），且具有大面积探测能力（直径可达数英寸），但量子效率较低（典型 20-30%），且对磁场敏感。

2.2.2 操作流程：SPAD探测器的选型与偏压设置

以SPAD为例，工程实施的关键步骤包括：

1. 选型评估：根据目标波长选择合适的光敏材料（Si、InGaAs/InP）。Si-SPAD适用于 400-1000 nm，InGaAs-SPAD适用于 900-1700 nm。
2. 偏压设置：将反向偏压设置为击穿电压以上 1-5 V（过压 ΔV）。过压越大，探测效率越高，但暗计数和后脉冲也增加。需通过实验标定最优过压点。
3. 淬灭电路设计：被动淬灭（电阻 100-500 kΩ）结构简单但恢复时间长；主动淬灭（快速比较器+晶体管）可缩短死时间至 10 ns 以下。
4. 信号读出：雪崩脉冲经高速比较器整形为LVTTL/NIM信号，送入时间数字转换器（TDC）或计数卡。
5. 温控与校准：使用TEC制冷至 -20°C 以下，记录暗计数率随温度变化曲线，确保工作点稳定。

2.2.3 案例拆解：1550 nm 长距离QKD系统中SNSPD与InGaAs-SPAD的对比

在某城域量子骨干网项目中，需要实现 100 km 光纤链路的QKD。初始方案采用InGaAs/InP SPAD（暗计数 1 kHz，探测效率 20%，死时间 50 ns），实测QBER为 8%，密钥率仅 1.2 kbps，无法满足业务需求。分析发现，主要瓶颈在于SPAD的暗计数和后脉冲在长距离弱光信号条件下导致误码率过高。

随后改用SNSPD（暗计数 <1 Hz，探测效率 85%，时间抖动 25 ps），配合门控模式（门宽 2 ns）和偏振编码方案。关键操作包括：

• 将SNSPD安装于闭循环制冷机中，温度稳定在 2.5 K ± 0.01 K；
• 优化偏置电流至 0.9 I_c（临界电流），使探测效率达到峰值；
• 采用符合计数法，将门控信号与探测器输出进行时间关联，进一步抑制暗计数。

最终系统QBER降至 1.5%，密钥率提升至 15 kbps，满足业务需求。该案例表明，在长距离、低光功率场景下，SNSPD的极低暗计数优势远大于其系统复杂度带来的代价。

2.2.4 常见误区与总结复盘

误区一：盲目追求高探测效率而忽略暗计数。在弱光探测中，信噪比（SNR）才是决定性指标，SNR = η·N_signal / √(η·N_signal + N_dark·t_window)。误区二：忽视时间抖动对系统的影响。在TCSPC和符合测量中，时间抖动直接决定时间分辨率，需根据实验需求选择合适抖动的探测器。误区三：认为PMT已过时。实际上，在需要大面积探测或超快时间响应的场景（如飞行时间质谱、正电子湮灭寿命谱），PMT仍是不可替代的选择。