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 输出电荷积分,得到与光子数近似成正比的模拟量,适合高光强下避免计数饱和。

三种模式对死时间的敏感程度不同:脉冲计数模式直接受死时间影响,必须修正;时间标记模式在低计数率下影响较小,但高计数率时仍会丢失时间戳;电荷积分模式则主要受微单元饱和限制,而非时间死时间。

单光子探测链路:器件 → 死时间 → 计数模式 → 修正 PMT 光电倍增管 光电阴极 → 打拿极链 → 阳极 增益 10⁵~10⁷ 输出:模拟脉冲,可甄别 死时间:约 10~100 ns 非瘫痪型为主 SPAD 单光子雪崩二极管 盖革模式,偏压 > 击穿电压 单载流子触发自持雪崩 输出:二值脉冲,无光子数分辨 死时间:20~100 ns(淬灭决定) 被动淬灭≈非瘫痪型 SiPM 硅光电倍增器 10³~10⁴ 个 SPAD 微单元并联 每单元串淬灭电阻 输出幅度 ∝ 触发单元数 死时间:单元级,约 50~200 ns 存在串扰、后脉冲、暗计数 死时间效应:观测计数率 m 与真实光子率 n 的偏离 非瘫痪型:m = n / (1 + n·τ) 死时间内光子被忽略,不延长死时间 瘫痪型:m = n·exp(−n·τ) 死时间内光子重新触发,延长死时间 脉冲计数模式 门宽 T 内统计脉冲数 k 受死时间影响最大 需计数率修正 时间标记模式 TCSPC 记录每个脉冲到达时刻 重建光子到达时间分布 低计数率下死时间影响小 电荷积分模式 积分 SiPM 输出电荷量 幅度 ∝ 光子数 受微单元饱和限制 修正目标:由 m 反推 n,消除死时间导致的计数损失

建模时的关键判断

把上述内容落到光子计数统计建模,需要依次确认三件事:第一,器件输出是二值脉冲还是幅度可分辨信号,这决定统计量是计数值还是电荷量;第二,死时间属于非瘫痪型还是瘫痪型,这决定修正公式的形式;第三,计数模式是脉冲计数、时间标记还是电荷积分,这决定死时间影响的作用位置。三者确定后,才能写出正确的观测模型,进而用泊松分布、死时间修正泊松分布或复合分布去拟合实验数据。

SiPM 的微单元饱和与死时间是两个不同机制:前者限制单次事件可分辨的最大光子数,后者限制单位时间可响应的光子数。高光强下两者同时起作用,建模时不能只修正其中一个。

实际实验中,建议先用弱光条件标定暗计数率和死时间,再逐步提高光强观察 m–n 曲线的偏离趋势,与理论公式对比确认器件类型。只有完成这一步,后续的光子计数统计分布拟合才具有物理意义。