3. 二阶关联函数定义与物理内涵
在光子统计测量中，我们真正关心的不是“某一时刻探测到几个光子”这一孤立事件，而是“在某一时刻探测到一个光子之后，另一时刻再探测到一个光子的概率有多大”。二阶关联函数正是把这种“成对探测”的统计规律定量化的工具。本节从归一化二阶关联函数 g²(τ) 的定义出发，推导它与聚束、反聚束、相干态三类光子统计的对应关系，并说明时间差 τ 的计算方式与归一化因子的物理意义。
一、从联合探测概率到 g²(τ) 的定义
设光场在 t 时刻的强度为 I(t)。理想单光子探测器在 t 时刻的探测概率正比于 I(t)。那么“在 t 时刻探测到一个光子，并在 t+τ 时刻再探测到一个光子”的联合概率，正比于强度关联函数：
G²(τ) = ⟨ I(t) · I(t+τ) ⟩其中尖括号表示对时间（或对系综）求平均。G²(τ) 的量纲是强度平方，其绝对数值依赖于光强、探测效率、光路损耗等实验条件，无法在不同实验之间直接比较。为了得到与装置无关、只反映光场本身统计性质的量，需要做归一化。最常用的归一化方式是除以两个单点强度平均值的乘积：
g²(τ) = ⟨ I(t) · I(t+τ) ⟩ / ⟨ I(t) ⟩²这就是归一化二阶关联函数的定义。若用光子数算符表示，在稳态光场下可写成 g²(τ) = ⟨a†(t)a†(t+τ)a(t+τ)a(t)⟩ / ⟨a†(t)a(t)⟩²。两种写法本质一致：分子是“先探测一个、再探测一个”的联合概率，分母是“两次独立探测”的参考概率。归一化因子 ⟨I⟩² 的物理意义正在于此——它把联合探测概率折算成“相对于随机独立探测的倍数”。因此 g²(τ) 是一个无量纲的相对量：g²=1 表示两次探测彼此独立，g²>1 表示成对出现的倾向增强，g²<1 表示成对出现的倾向被抑制。
二、三类光子统计的判据
根据 g²(τ) 在 τ=0 附近的取值，可以把光场的光子统计分为三类，这是二阶关联测量最核心的判据。
- 热光（聚束）：g²(0)=2，g²(τ) 随 |τ| 增大按 1+exp(−|τ|/τc) 的形式衰减到 1。光子倾向于“扎堆”到达，称为光子聚束。典型对象是热光源、混沌光、以及多模热态。
- 相干态：g²(τ)≡1，与 τ 无关。相干态的光子到达服从泊松分布，前后探测事件完全独立，是“随机基准”。理想单模激光在阈值以上近似相干态。
- 单光子源（反聚束）：g²(0)=0，g²(τ) 在 τ=0 处出现凹陷，随 |τ| 增大回升到 1。一个单光子源发射一个光子后需要重新激发，无法同时发射两个光子，因此“同时探测到两个光子”的概率被压制，称为光子反聚束。g²(0)<0.5 是判定单光子源的重要实验门槛。
需要强调，判据看的是 τ=0 附近的行为，而不是任意 τ 处的值。对于脉冲式单光子源，g²(τ) 会在每个脉冲重复周期处出现一系列峰，但 τ=0 处的凹陷才是反聚束的本质特征。
三、时间差 τ 的计算方式
在 HBT（Hanbury Brown–Twiss）型测量中，光场被分束器分成两路，分别送入两个单光子探测器，两路探测器的输出电脉冲送入时间相关单光子计数器（TCSPC）或符合计数器。τ 的定义必须与硬件记录方式严格对应：
τ = t₂ − t₁其中 t₁ 是“起始通道”（start）探测到光子的时刻，t₂ 是“停止通道”（stop）探测到光子的时刻。若起始通道接探测器 A、停止通道接探测器 B，则 τ>0 表示 B 比 A 晚探测到光子，τ<0 表示 B 比 A 早。由于两个探测器在物理上是对称的，理想情况下 g²(τ) 应满足 g²(τ)=g²(−τ)，因此实际分析中常把正负延迟合并，或只取 |τ| 作图以提升信噪比。但要注意：若两路存在固定的电子延迟差或光纤长度差，曲线整体会沿 τ 轴平移，这个平移量必须在标定中扣除，否则 τ=0 的位置会判断错误，进而误判 g²(0)。
时间轴的分辨率由探测器抖动（jitter）和 TCSPC 的时间通道宽度共同决定。若时间通道宽度为 Δt，则每个 τ 的取值落在 [kΔt, (k+1)Δt) 区间内，直方图的 bin 宽度就是 Δt。对于反聚束测量，若单光子源的寿命远小于 Δt，τ=0 处的凹陷会被时间分辨“抹平”，测得的 g²(0) 会偏高，因此 Δt 应显著小于光源的特征时间尺度。
四、归一化因子的物理意义与实验估计
归一化因子 ⟨I⟩² 在实验上并不是直接测得的强度平方，而是通过“远延迟参考”来估计。当 |τ| 远大于光场的相干时间（或脉冲源的重复周期）时，两次探测事件彼此独立，此时 G²(τ→∞)=⟨I⟩²。因此实际数据处理中常用长延迟区的符合计数平均值作为分母：
g²(τ) = N_符合(τ) / [ N₁ · N₂ · Δt · T ]其中 N_符合(τ) 是延迟 τ 处的符合计数，N₁、N₂ 是两路的单路计数率，Δt 是时间通道宽度，T 是总测量时间。分母 N₁·N₂·Δt·T 正是“若两路完全独立，在 τ 通道内期望的偶然符合数”。用这个分母归一化，g²(τ) 就直接表示“实测符合数相对偶然符合数的倍数”：等于 1 说明符合纯属偶然，大于 1 说明存在超出偶然的成对倾向，小于 1 说明成对事件被压制。这也解释了为什么 g²(0) 可以小于 1 甚至接近 0——它衡量的是“同时到达”相对于随机基准被压低的程度。
五、一个教学示例
假设某脉冲单光子源的重复频率为 80 MHz（周期 12.5 ns），探测器时间抖动约 300 ps，TCSPC 通道宽度设为 16 ps。测量得到：单路计数率 N₁=N₂=5×10⁴ s⁻¹，总测量时间 T=600 s，τ=0 通道内符合计数为 12，而长延迟区（|τ|>50 ns）平均每通道符合计数为 480。则按长延迟归一化，g²(0)=12/480=0.025，远小于 0.5，可判定该源具有良好的单光子性。若把通道宽度改为 512 ps，τ=0 通道会吞入更多真实符合，g²(0) 可能升至 0.2 左右——数值变化并非光源变差，而是时间分辨改变所致。这个例子说明，报告 g²(0) 时必须同时给出时间通道宽度与归一化方式，否则数值不可复现。
综上，g²(τ) 的定义把“成对探测”折算为相对独立探测的倍数，其 τ=0 附近的取值直接对应聚束、相干、反聚束三类统计；τ 的计算依赖 start/stop 通道约定与系统延迟标定，归一化因子则决定了 g² 数值的参照基准。理解这三者的对应关系，是后续符合计数处理与 g²(0) 判读的基础。