📚 课程目录第2章单光子探测物理基础

2. 单光子探测物理基础

在时间相关单光子计数（TCSPC）系统中，探测器面对的从来不是“稳定连续的光流”，而是随机到达的离散光子。理解这种随机性，是理解计数率上限、计数损失、寿命测量偏差以及系统死时间的起点。本节从光子统计出发，逐步建立泊松过程、探测概率、死时间与堆积效应之间的定量联系，并给出可用于工程估算的示例参数。

一、光子为什么服从泊松统计

对于理想的稳定光源（如连续波激光经强衰减后），在给定时间窗口内到达探测器的光子数并不是固定值，而是围绕平均值波动的随机变量。若满足以下条件：光子到达彼此独立、单位时间平均到达率恒定、任意两个光子在同一瞬间到达的概率可忽略，则时间窗口 T 内到达的光子数 N 服从泊松分布：

P(N = k) = (λT)^k · e^(−λT) / k!

λ：平均光子到达率（photons/s）
T：观测时间窗口（s）
λT：窗口内平均光子数

泊松分布有两个关键特征：均值等于方差，即 Var(N) = λT；相邻光子到达时间间隔服从指数分布，其均值为 1/λ。这意味着“光子间隔”本身也是随机的，短间隔出现的概率并不为零——这正是后面死时间与堆积效应的物理根源。

二、探测概率与量子效率

光子到达探测器并不等于被探测到。单光子探测器（SPAD、PMT、SNSPD 等）存在量子效率 η，即入射光子产生可探测雪崩/脉冲的概率。若入射平均光子率为 λ，则有效探测事件率近似为：

μ = η · λ

η：量子效率（0～1）
μ：探测事件平均率（counts/s）

在时间窗口 T 内，探测到 k 个事件的概率同样服从泊松分布，只需把 λT 替换为 μT。工程上常关心“至少探测到一个光子”的概率：

P(N ≥ 1) = 1 − e^(−μT)

例如，若 μ = 1×10⁶ counts/s，T = 1 ns，则 μT = 0.001，P(N≥1) ≈ 0.001，即约千分之一的窗口能记录到事件。这个数量级说明：在高速 TCSPC 中，绝大多数激发周期是“空”的，系统必须依靠大量重复周期累积统计。

三、死时间：探测器“失明”的窗口

SPAD 在每次雪崩后需要一定时间恢复，这段时间内无法响应新光子，称为死时间 τ_d。典型值：

硅基 SPAD：20～100 ns（被动淬灭），部分主动淬灭可到 10 ns 量级；
PMT：1～20 ns；
SNSPD：几 ns 到几十 ns。

死时间分为两类：非扩展型（non-paralyzable）——死时间内到达的光子被忽略，死时间不延长；扩展型（paralyzable）——每个到达光子都重新触发死时间，可能使探测器永久饱和。SPAD 通常接近非扩展型。

对非扩展型探测器，若真实入射事件率为 μ，测得计数率为 m，则：

m = μ / (1 + μ τ_d)

反解：μ = m / (1 − m τ_d)

该式说明：当 μ 接近 1/τ_d 时，m 趋于饱和值 1/τ_d。例如 τ_d = 50 ns，则理论最大计数率约 2×10⁷ counts/s；若测得 m = 1×10⁷ counts/s，则真实 μ ≈ 1×10⁷ / (1 − 0.5) = 2×10⁷ counts/s，计数损失高达 50%。

四、堆积效应：脉冲重叠与计数偏差

堆积（pile-up）与死时间密切相关但侧重点不同。死时间描述“探测器不能响应”，堆积描述“多个光子事件在时间上重叠，导致系统把它们当成一个事件或产生错误时间标记”。在 TCSPC 中，堆积最典型的后果是：高计数率下，晚到达的光子被早到达光子的死时间“屏蔽”，导致荧光衰减曲线整体向短寿命方向偏移。

判断堆积是否可忽略，常用“计数率准则”：

若 μ · τ_d ≪ 1，则堆积可忽略。
工程经验：μ · τ_d < 0.01～0.05 时，寿命偏差通常小于 1%。

例如 τ_d = 50 ns，要求 μ·τ_d < 0.01，则 μ < 2×10⁵ counts/s。若实际需要更高计数率，就必须做死时间/堆积校正，或降低单周期激发概率（降低激光功率、增大衰减）。

五、参数关系总览

下面用一张关系图把本节核心量串起来：入射光子率 λ、量子效率 η、探测事件率 μ、死时间 τ_d、测得计数率 m、堆积程度 μτ_d 之间的因果与反馈关系。图中同时标出“可忽略区”和“需校正区”的边界，便于工程判断。

六、工程判断与操作建议

在实际 TCSPC 系统中，建议按以下顺序检查光子统计与堆积：

估算 μ：由激光重复频率 f_rep、每周期平均光子数 n̄ 得 μ = f_rep · n̄ · η。例如 f_rep = 40 MHz，n̄ = 0.01，η = 0.5，则 μ = 2×10⁵ counts/s。
计算 μτ_d：若 τ_d = 50 ns，则 μτ_d = 0.01，处于可忽略边界，寿命偏差通常小于 1%。
观察计数率线性：逐步增大激光功率，记录 m 随 μ 的变化。若 m 开始偏离线性、趋于饱和，说明进入堆积区。
校正或降功率：若无法降功率，用 μ = m/(1−mτ_d) 反解真实率，或采用多周期/多探测器并行降低单通道计数率。

教学示例参数说明：以上 τ_d = 50 ns、η = 0.5、f_rep = 40 MHz 等均为典型量级示例，用于说明计算过程，不代表某一具体型号的实测指标。实际系统应以探测器数据手册和实测计数率为准。

总结来看，光子统计决定了 TCSPC 的“随机底噪”，泊松过程给出了探测概率的定量描述，死时间设定了计数率上限，堆积效应则把“高计数率”转化为“寿命偏差”。四者环环相扣：只有把 μτ_d 控制在合理范围，后续的寿命拟合与系统校准才有意义。