3. 原子荧光与散射基础

中性原子阵列的成像，本质上是在“数光子”：原子被近共振激光照射后不断散射光子，其中一小部分被物镜收集并落到相机上，形成荧光斑。要判断一个原子处于亮态还是暗态，或者判断某个格点是否有原子，最终都归结为相机像素上收集到的光电子数是否超过噪声。理解共振散射、饱和散射率、光子收集效率、信噪比与量子效率之间的关系，是设计成像方案和估计读出保真度的基础。

一、共振散射：从散射截面到散射率

考虑一个二能级近似下的原子，基态 |g⟩、激发态 |e⟩，共振频率为 ω₀，自然线宽为 Γ（对应激发态寿命 τ = 1/Γ）。当一束频率为 ω_L、强度为 I 的激光照射原子时，原子在激光驱动下做拉比振荡，同时被自发辐射衰减。稳态下，原子处于激发态的布居为

ρ_ee = (s₀ / 2) / (1 + s₀ + (2Δ/Γ)²)

其中 s₀ = I / I_sat 为共振饱和参数，
I_sat = π h c Γ / (3 λ³) 为饱和强度，
Δ = ω_L − ω₀ 为激光失谐。

原子单位时间散射的光子数，即散射率，为

R_sc = Γ · ρ_ee = (Γ / 2) · s₀ / (1 + s₀ + (2Δ/Γ)²)

这个式子包含三个关键信息。第一，当 Δ = 0 且 s₀ ≫ 1 时，R_sc → Γ/2，这就是饱和散射率的极限：原子最多每 2τ 散射一个光子，因为激发态布居最多到 1/2。第二，失谐会显著压低散射率，Δ = Γ/2 时散射率约降为共振时的一半。第三，弱光下 s₀ ≪ 1，R_sc ≈ (Γ/2)·s₀，散射率与光强成正比，此时原子像一个线性散射体。

以 ⁸⁷Rb 的 D2 线为例，λ ≈ 780 nm，Γ ≈ 2π × 6.07 MHz，τ ≈ 26.2 ns，I_sat ≈ 1.67 mW/cm²（对圆偏振、均匀光场）。若用 s₀ = 10 的共振光照射，R_sc ≈ Γ/2 × 10/11 ≈ 0.45Γ ≈ 2π × 2.76 MHz，即每秒约 2.76 × 10⁶ 个光子。若把光强降到 s₀ = 0.1，散射率降到约 0.045Γ，每秒约 2.76 × 10⁵ 个光子。可见光强直接决定“光子预算”。

二、光子收集效率：几何与光学系统的乘积

原子向 4π 立体角各向同性辐射，物镜只能收集其中一部分。若物镜半孔径角为 θ，数值孔径 NA = n sinθ（空气中 n ≈ 1），则几何收集效率为

η_geo = (1 − cosθ) / 2 = (1 − √(1 − NA²)) / 2

NA = 0.5 时 η_geo ≈ 6.7%；NA = 0.75 时约 11.7%；NA = 0.95 时约 19.5%。这是单透镜的物理上限。实际收集效率还要乘上光学系统的透过率 T（窗口、滤光片、分束器、物镜镀膜等），以及原子与物镜焦点的对准效率。因此

η_coll = η_geo · T · η_align

典型高 NA 物镜加滤光片系统，T 约 0.7～0.9，对准效率接近 1，总收集效率约 5%～15%。这意味着即使原子每秒散射 2.76 × 10⁶ 个光子，落到相机上的可能只有每秒 1.4 × 10⁵～4 × 10⁵ 个光子。曝光时间 t_exp 内收集到的光子数为 N_ph = R_sc · η_coll · t_exp。

三、量子效率与信噪比：从光子到电子

相机（sCMOS 或 EMCCD）把入射光子转换为光电子，量子效率 QE 表示单个光子产生一个可探测电子的概率。于是曝光时间内产生的光电子数为

N_e = N_ph · QE = R_sc · η_coll · t_exp · QE

信号就是 N_e。噪声来源主要有三类：光子散粒噪声 √N_e、暗电流噪声 √(D·t_exp)（D 为暗电流率）、读出噪声 σ_r。总噪声为

σ_total = √(N_e + D·t_exp + σ_r²)

信噪比 SNR = N_e / σ_total。在散粒噪声主导时，SNR ≈ √N_e，即信号越强，相对噪声越小，但改善速度是平方根关系。若读出噪声主导（弱信号、短曝光），SNR ≈ N_e / σ_r，此时提高 QE 或延长曝光比单纯增加光强更有效。

以一个具体示例说明：假设 R_sc = 2.76 × 10⁶ s⁻¹，η_coll = 10%，QE = 0.8，t_exp = 10 μs。则 N_ph = 2.76 × 10⁶ × 0.1 × 10 × 10⁻⁶ = 2.76，N_e ≈ 2.2。散粒噪声 √2.2 ≈ 1.5，若读出噪声 σ_r = 1.5 e⁻，则 σ_total ≈ √(2.2 + 2.25) ≈ 2.1，SNR ≈ 1.05。这个信噪比不足以可靠区分“有原子”和“无原子”。若把曝光延长到 100 μs，N_e ≈ 22，散粒噪声 4.7，σ_total ≈ √(22 + 2.25) ≈ 4.9，SNR ≈ 4.5，单次判断的保真度就明显提高。若再提高收集效率到 20%，N_e ≈ 44，SNR ≈ 6.4。

四、读出判据与保真度

实际判断一个格点是否有原子，通常设定一个阈值 N_th。若测得光电子数 N_e > N_th，判为“有原子”（亮）；否则判为“无原子”（暗）。由于散粒噪声和读出噪声的存在，亮态和暗态的光电子数分布会有重叠，重叠区域导致误判。假设亮态均值 N₁、暗态均值 N₀（通常接近 0），两者标准差分别为 σ₁、σ₀，则阈值取在 (N₁ + N₀)/2 附近时，误判率由高斯尾部积分给出。若 N₁ = 22、σ₁ = 4.9、N₀ = 0、σ₀ = 1.5，阈值取 11，则亮态误判为暗的概率约为 Φ(−(22−11)/4.9) ≈ Φ(−2.24) ≈ 1.3%，暗态误判为亮的概率约为 Φ(−11/1.5) ≈ Φ(−7.3) ≈ 10⁻¹³，单次读出保真度约 98.7%。若把曝光延长到 200 μs，N₁ ≈ 44，σ₁ ≈ 6.7，阈值取 22，亮态误判率降到约 0.05%，保真度超过 99.9%。

这里有一个重要的工程权衡：延长曝光可以提高 SNR，但会降低成像帧率，也可能引入原子加热、丢失或光位移。提高光强可以增加 R_sc，但接近饱和后收益递减，而且过强光会导致原子加热和碰撞损失。提高收集效率需要更大 NA 的物镜，但工作距离和视场会受限。提高 QE 依赖相机选型，sCMOS 在 600～800 nm 的 QE 通常为 70%～90%，EMCCD 更高但读出噪声和成本也更高。实际方案往往是在这些约束下取折中。

总结来看，共振散射决定了原子的“光子发射上限”，饱和散射率给出了这个上限的物理边界；光子收集效率把各向同性辐射压缩到物镜孔径内；量子效率完成光子到光电子的转换；散粒噪声、暗电流和读出噪声共同决定信噪比；最终由信噪比和阈值判据决定读出保真度。这条链路中每一环都可以独立优化，但彼此耦合，理解它们的定量关系，才能在原子阵列成像实验中做出合理的参数选择。