4、暗计数基础：暗计数的来源、温度依赖性与DCR统计特性
SPAD（单光子雪崩二极管）在完全没有入射光子的条件下，仍然会产生雪崩脉冲，这些脉冲被称为暗计数（Dark Count）。暗计数是SPAD噪声的主要来源，直接决定了探测器的暗计数率（DCR，Dark Count Rate），是衡量SPAD性能的核心指标之一。理解暗计数的物理来源、温度依赖规律以及统计特性，是进行死时间优化和系统信噪比提升的前提。本章将从这三个维度展开，帮助工程师建立对暗计数噪声的系统认知。
4.1 暗计数的物理来源与产生机制
SPAD的暗计数并非单一机制所致，而是多种物理过程共同作用的结果。在工程实践中，不同来源的暗计数对温度、偏压的响应规律不同，因此需要分别理解其物理本质，才能有针对性地进行优化。
4.1.1 热激发（SRH复合中心辅助）
热激发是室温下SPAD暗计数最主要的来源。在PN结耗尽区内，由于晶格缺陷或杂质的存在，会形成Shockley-Read-Hall（SRH）复合中心。这些缺陷能级位于禁带中间附近，充当了电子-空穴对产生的"跳板"：价带电子先被热激发到缺陷能级，再从缺陷能级跃迁到导带，形成自由载流子。该过程所需的激活能约为禁带宽度的一半（硅约0.55eV），因此对温度极其敏感。
在SPAD结构中，耗尽区宽度越宽、缺陷密度越高，热激发产生的载流子数量就越多。实际器件中，工艺引入的金属杂质（如Fe、Cu）和晶格位错是主要缺陷来源。从工程角度看，降低热激发暗计数的手段包括：优化外延层质量、采用"优区"结构（在倍增区周围设置保护环）、以及通过降温运行来指数级压低热激发速率。
4.1.2 带间隧穿（Band-to-Band Tunneling）
当SPAD工作在较高的过偏压（VEX）时，PN结内部电场强度可达10⁶ V/cm量级。在如此高的电场下，价带电子可以直接隧穿禁带进入导带，形成电子-空穴对，这就是带间隧穿效应。与热激发不同，带间隧穿几乎不依赖温度，而是对电场强度（即偏压）极其敏感。
带间隧穿电流密度可用以下简化公式描述：
J_BTBT ≈ A · E · exp(-B / E)其中E为结电场强度，A和B为材料相关常数。从公式可以看出，电场强度E出现在指数项的分母中，意味着电场微小的增加会导致隧穿电流的急剧上升。因此，在追求高光子探测效率（PDP）而提高过偏压时，带间隧穿暗计数会呈超线性增长，成为限制SPAD最大工作偏压的关键因素。
4.1.3 陷阱辅助隧穿（Trap-Assisted Tunneling, TAT）
陷阱辅助隧穿是介于热激发和带间隧穿之间的机制。载流子先被热激发到缺陷能级，然后在电场作用下隧穿出陷阱进入导带（或价带）。TAT的激活能低于SRH热激发，通常为禁带能量的1/3到1/2，因此其温度依赖性介于纯热激发和纯隧穿之间。在中等偏压和中等温度区间，TAT往往占据主导地位。
从器件设计角度看，TAT与工艺缺陷密度直接相关。采用"浅结"结构或降低倍增区掺杂浓度，可以有效降低峰值电场，从而抑制TAT。此外，在SPAD外围增加"虚拟环"（Virtual Guard Ring）结构，可以避免边缘电场集中导致的TAT增强。
4.2 暗计数率的温度依赖性
DCR随温度的变化规律是SPAD系统设计中最关键的工程参数之一。理解这一依赖关系，不仅有助于选择合适的工作温度，还能指导制冷方案的功耗预算和温控精度设计。
4.2.1 理论模型与经验公式
热激发主导的DCR温度依赖性可以用Arrhenius公式描述：
DCR(T) = C · T² · exp(-Ea / kT)其中C为常数，T为温度，Ea为激活能（硅SPAD典型值0.5~0.6eV），k为玻尔兹曼常数（8.617×10⁻⁵ eV/K）。T²项来源于载流子热速度和有效态密度的温度依赖。在实际工程中，更常用的是"每10°C倍增因子"（M10）来快速估算：
M_10 = exp( Ea · ln(10) / (k · T · (T+10)) )在室温（300K）附近，Ea=0.55eV时，M10≈2.0，即温度每升高10°C，DCR翻倍。这一经验法则在-20°C到+60°C范围内具有较好的工程精度。
4.2.2 低温平台的产生
当温度降低到一定程度（通常低于-40°C），热激发速率指数下降，此时带间隧穿和TAT成为暗计数的主要来源。由于这两种机制几乎不依赖温度，DCR随温度降低而趋于饱和，形成"低温平台"。这意味着单纯依靠制冷无法将DCR无限降低，存在一个由隧穿机制决定的"地板值"。
对于需要极低暗计数的应用（如量子密钥分发QKD、单分子检测），设计者需要在器件选型时关注其低温平台DCR水平。例如，深制冷到-60°C时，优秀硅SPAD的DCR可低至1 cps以下，而InGaAs/InP SPAD在同样温度下平台DCR通常在100~1000 cps量级。
4.2.3 温度依赖性的工程测量方法
在实验室中准确测量DCR的温度依赖性，需要遵循以下步骤：
- 暗室环境准备：确保SPAD完全避光，使用遮光罩和暗箱，排除环境光子干扰。
- 温度控制：使用TEC（半导体制冷器）或恒温槽，温度稳定度优于±0.1°C，避免温度波动导致DCR测量误差。
- 偏压设置：固定过偏压（如VEX=3V），记录击穿电压VBD随温度的变化（典型硅SPAD约0.3~0.5V/°C），确保过偏压恒定。
- 数据采集：每个温度点至少采集60秒以上，统计总暗计数脉冲数，除以有效探测时间得到DCR。
- 数据拟合：将DCR-T数据按Arrhenius公式拟合，提取激活能Ea，判断暗计数主导机制。
4.3 暗计数率的统计特性
暗计数在时间轴上的分布并非完全随机，而是包含泊松随机分量和相关性分量（后脉冲）。理解这些统计特性，对于设计死时间策略、信号处理算法和系统误码率评估至关重要。
4.3.1 泊松分布与时间均匀性
在排除后脉冲影响的理想情况下，暗计数事件在时间上服从泊松分布。这意味着：在任意时间间隔Δt内，出现k个暗计数的概率为：
P(k) = (λ·Δt)^k · exp(-λ·Δt) / k!其中λ为平均暗计数率（即DCR）。泊松分布的关键特性是：方差等于均值（σ²=μ），且事件之间相互独立。在工程上，这意味着暗计数的时间间隔服从指数分布，其概率密度函数为：
f(t) = λ · exp(-λ·t)这一特性为死时间设计提供了理论依据：如果死时间设置为τD，则因暗计数导致的"有效探测时间损失"比例约为λ·τD（当λ·τD << 1时）。
4.3.2 后脉冲效应与相关性
后脉冲（Afterpulsing）是SPAD特有的相关性噪声。在一次雪崩过程中，部分载流子被陷阱能级捕获，在随后（微秒到毫秒量级）被释放，触发新的雪崩。后脉冲的概率PAP随死时间τD的增大而指数下降：
P_AP(τ_D) = P_AP0 · exp(-τ_D / τ_trap)其中τtrap为陷阱寿命（硅SPAD典型值50~500ns），PAP0为初始后脉冲概率（与过偏压和工艺相关）。后脉冲的存在使得暗计数在时间上呈现"聚集"效应，破坏了泊松分布的独立性假设。
4.3.3 DCR统计特性的测量与评估
为了全面评估SPAD的暗计数统计特性，建议采用以下测量方案：
| 测量参数 | 测量方法 | 评估指标 | 工程意义 |
|---|---|---|---|
| 平均DCR | 长时间计数（≥60s） | DCR (cps) | 系统噪声基底 |
| 时间间隔分布 | 时间相关单光子计数（TCSPC） | 指数拟合优度R² | 判断是否纯泊松 |
| 后脉冲概率 | 双脉冲相关测量 | PAP (%) | 死时间优化依据 |
| 计数方差/均值比 | 分段计数统计 | Fano因子F=σ²/μ | 偏离泊松的程度 |
4.4 案例拆解：一款硅SPAD的暗计数优化实战
下面通过一个完整的工程案例，展示如何利用暗计数知识进行系统优化。
4.4.1 项目背景与初始现象
某激光雷达（LiDAR）项目采用单光子探测方案，使用商用硅SPAD（有效直径25μm，击穿电压VBD=24.5V@25°C）。系统要求在10m距离内实现90%的探测概率，同时保持误警率低于1kHz。初始样机测试发现：在25°C环境温度下，DCR实测为8.5kcps，远高于数据手册标称的2kcps（@25°C），导致误警率超标。
4.4.2 问题分析与诊断过程
工程师首先进行了变温实验，结果如下：
- 15°C时DCR=4.2kcps，25°C时DCR=8.5kcps，35°C时DCR=17.8kcps
- 计算M10 = 17.8/8.5 ≈ 2.09，符合热激发主导特征
- 但25°C时DCR仍比手册值高4倍以上
进一步检查发现：样机PCB布局中，SPAD附近有一个功率电阻（用于LED驱动），其表面温度实测为62°C，通过热传导使SPAD局部温度比环境温度高约8°C。修正后SPAD实际结温约33°C，对应DCR应为8.5kcps，与实测吻合。
4.4.3 优化措施与结果
针对热耦合问题，采取了以下措施：
- 热隔离：在SPAD与功率电阻之间增加隔热槽（PCB开槽），并在SPAD底部增加铜散热块连接至外壳。
- 温度补偿：增加NTC热敏电阻紧贴SPAD封装，实时监测结温，动态调整过偏压以保持恒定VEX。
- 死时间优化：根据后脉冲测量结果（PAP=1.2% @τD=100ns），将死时间从100ns调整为80ns，在保证后脉冲概率不显著增加的前提下，提高了最大计数率。
优化后，在25°C环境温度下，SPAD实际结温降至27°C，DCR降至3.1kcps，误警率降至0.6kHz，满足系统指标。同时，由于死时间缩短，最大可测光强提升了25%。
4.4.4 复盘与经验总结
本案例揭示了暗计数优化中两个容易被忽视的要点：
- 热管理优先于电气优化：在许多系统中，SPAD的"异常"暗计数往往源于局部热源，而非器件本身缺陷。在投入复杂的算法或电路优化前，应先进行热成像检查。
- 死时间并非越大越好：死时间增加虽然降低后脉冲，但会牺牲最大计数率和增加光子丢失概率。最优死时间应基于后脉冲概率曲线和系统最大光通量需求综合确定。
4.5 暗计数优化的常见误区与设计建议
在SPAD系统设计中，暗计数相关的误区往往导致设计反复或性能不达标。以下总结了几类高频问题：
| 误区类型 | 典型表现 | 正确做法 |
|---|---|---|
| 忽视击穿电压温度漂移 | 固定偏压导致过偏压随温度变化 | 实时监测温度，动态补偿VBD |
| 混淆DCR与暗电流 | 用直流暗电流推算DCR | DCR必须通过脉冲计数测量 |
| 忽略后脉冲的统计影响 | 假设暗计数纯泊松分布 | 测量后脉冲概率，纳入噪声模型 |
| 过度制冷 | 追求极低温度忽略功耗和冷凝 | 根据低温平台DCR选择合适工作点 |
| 死时间设置一刀切 | 所有像素使用相同死时间 | 根据各像素后脉冲特性差异化设置 |
在系统设计层面，建议将暗计数优化视为一个"温度-偏压-死时间"三维联合优化问题。首先通过变温实验确定器件的最优工作温度区间，然后在该温度下扫描过偏压，找到DCR与光子探测效率（PDP）的最佳折中点，最后根据后脉冲测量结果确定死时间。这一流程虽然耗时，但能确保系统在真实工作条件下达到最优性能。
此外，对于阵列型SPAD系统，建议在固件中实现"暗计数实时校准"功能：在每次测量前插入一段短时间（如1ms）的暗计数采样，动态估计当前DCR，并据此调整触发阈值或进行背景扣除。这种方法能有效应对温度漂移和器件老化带来的DCR变化，提升系统的长期稳定性。