3. 单光子探测器件与计数模式
单光子探测是光子计数统计建模的物理起点。探测器把“有没有光子到达”转换成“有没有电脉冲输出”,而统计建模真正关心的,是输出脉冲序列与入射光子序列之间的映射关系。理解 SPAD、PMT、SiPM 三种主流器件的工作机理,以及死时间如何扭曲这一映射,是后续建立计数分布模型、做计数率修正的前提。
三种器件的工作机理
PMT(光电倍增管)是最早成熟的单光子探测器。光子打到光电阴极,通过外光电效应激发出光电子;光电子在打拿极链中被逐级二次电子倍增,典型增益可达 10⁵~10⁷,最终在阳极形成可被甄别的电流脉冲。PMT 增益高、响应快、暗计数低,但体积大、需要数百至上千伏高压、对磁场敏感,且量子效率在长波段下降明显。
SPAD(单光子雪崩二极管)工作在盖革模式,即反向偏压高于击穿电压。单个光生载流子即可触发自持雪崩,输出与光子到达时刻对应的宏观电流脉冲,因此 SPAD 天然具备单光子灵敏度和皮秒级时间分辨能力。但雪崩一旦触发就会持续,必须由淬灭电路(被动或主动)把偏压拉回击穿电压以下才能复位,这段无法响应新光子的时间就是死时间。SPAD 输出是二值化的“有/无”信号,本身不区分光子数。
SiPM(硅光电倍增器)由大量并联的 SPAD 微单元(典型 10³~10⁴ 个,单元尺寸 10~100 μm)集成在同一硅片上,每个微单元串联一个淬灭电阻。光子击中某个微单元时,该单元雪崩并输出一个标准化脉冲;所有单元的输出在公共负载上叠加,因此 SiPM 的脉冲幅度正比于同时被触发的微单元数,具备一定的光子数分辨能力。SiPM 工作电压低(约 30 V)、体积小、对磁场不敏感,但存在光学串扰、后脉冲和较高的暗计数率,且动态范围受微单元总数限制。
死时间效应:非瘫痪型与瘫痪型
死时间是单光子探测中最关键的计数损失机制。按器件在死时间内是否仍能响应新光子,分为两类:
- 非瘫痪型(non-paralyzable):死时间内到达的光子被完全忽略,但不会延长死时间。SPAD 的被动淬灭、多数 PMT 计数模式近似属于此类。
- 瘫痪型(paralyzable):死时间内到达的光子会重新触发并延长死时间,计数率随入射率升高反而下降。部分主动淬灭 SPAD 和某些 PMT 工作模式接近此类。
设真实入射光子率为 n(单位 cps),死时间为 τ,则两种模型的观测计数率 m 分别为:
非瘫痪型:m = n / (1 + n·τ)
瘫痪型: m = n·exp(−n·τ)
两式都表明:当 n·τ 不可忽略时,m 明显小于 n,且随 n 增大逐渐饱和。非瘫痪型在高计数率下趋于 1/τ 的饱和值;瘫痪型则先升后降,存在一个极大值。实际建模时,必须先判断器件属于哪一类,再选择对应修正公式。
计数率修正的实际操作
修正的目标是由观测计数率 m 反推真实光子率 n。对非瘫痪型,可直接解析反解:
n = m / (1 − m·τ)
对瘫痪型,解析反解需要 Lambert W 函数,工程上更常用迭代或查表。以非瘫痪型为例,若某 SPAD 死时间 τ = 50 ns,实测 m = 8×10⁶ cps,则:
n = 8e6 / (1 − 8e6 × 50e-9)
= 8e6 / (1 − 0.4)
= 1.333e7 cps
可见观测值比真实值低了约 40%,若不修正,后续的光子统计分布拟合会出现系统性偏差。死时间 τ 的标定通常用双脉冲法或衰减法:用已知时间间隔的两个光脉冲测量最小可分辨间隔,或改变入射光强拟合 m–n 曲线反推 τ。
从脉冲到计数:三种计数模式
探测器输出脉冲后,计数方式决定了最终统计量的形式:
- 脉冲计数模式:在固定门宽 T 内统计脉冲个数,得到计数值 k,服从泊松或经死时间修正后的分布。这是光子计数统计最常用的模式。
- 时间标记模式(TCSPC):记录每个脉冲的到达时刻,形成时间戳序列,可重建光子到达时间分布,用于荧光寿命、飞行时间测量。
- 电荷积分模式:对 SiPM 输出电荷积分,得到与光子数近似成正比的模拟量,适合高光强下避免计数饱和。
三种模式对死时间的敏感程度不同:脉冲计数模式直接受死时间影响,必须修正;时间标记模式在低计数率下影响较小,但高计数率时仍会丢失时间戳;电荷积分模式则主要受微单元饱和限制,而非时间死时间。
建模时的关键判断
把上述内容落到光子计数统计建模,需要依次确认三件事:第一,器件输出是二值脉冲还是幅度可分辨信号,这决定统计量是计数值还是电荷量;第二,死时间属于非瘫痪型还是瘫痪型,这决定修正公式的形式;第三,计数模式是脉冲计数、时间标记还是电荷积分,这决定死时间影响的作用位置。三者确定后,才能写出正确的观测模型,进而用泊松分布、死时间修正泊松分布或复合分布去拟合实验数据。
实际实验中,建议先用弱光条件标定暗计数率和死时间,再逐步提高光强观察 m–n 曲线的偏离趋势,与理论公式对比确认器件类型。只有完成这一步,后续的光子计数统计分布拟合才具有物理意义。