4. 灵敏度理论极限公式推导

在上一节中，我们把原子磁力计的灵敏度归结为“能分辨的最小磁场变化”。要回答这个“最小”究竟由什么决定，就必须回到量子层面。原子磁力计的灵敏度理论极限并非来自电路噪声或激光强度噪声，而是来自两个不可回避的量子涨落源：自旋噪声（原子系综自身磁矩的统计涨落）与量子投影噪声（测量过程对自旋态投影引入的随机性）。本节从这两个源头出发，推导灵敏度基本表达式，并解释标准量子极限的物理内涵。

一、自旋噪声：系综磁矩的统计涨落

考虑一个包含 N 个自旋-1/2 原子的系综，沿测量轴（设为 z 轴）的集体自旋算符为 J_z = (ħ/2)Σσ_z^(i)。当原子处于相干自旋态时，横向自旋分量 J_x 与 J_y 的期望值最大，而 J_z 的期望值为零，但其方差不为零。对于 N 个独立原子，J_z 的方差为：

⟨ΔJ_z²⟩ = N·(ħ/2)² = Nħ²/4

这就是自旋噪声的根源：即使没有任何外场扰动，集体磁矩的 z 分量也在 ±ħ/2·√N 量级上随机涨落。若把磁矩换算为等效磁场，这一涨落对应一个“自旋噪声磁场” δB_sn。由于磁矩 μ = γJ（γ 为旋磁比），自旋噪声等效磁场满足：

δB_sn = √(⟨ΔJ_z²⟩)·(2/(γ·N·ħ)) = 1/(γ√(N·T₂))

其中 T₂ 为横向弛豫时间。注意这里已把“单次测量”的时间尺度纳入：在 T₂ 时间内，原子自旋保持相干，测量才能积累信号。因此自旋噪声随 1/√N 下降，这是系综平均的经典结果。

二、量子投影噪声：测量本身引入的随机性

仅有自旋噪声还不够。任何对 J_z 的测量都不是无扰动的：测量把自旋态投影到 J_z 的本征态上，单次测量结果是一个随机变量，其方差至少为 ⟨ΔJ_z²⟩。若用一束探测光连续测量，单位时间内获得的光子数为 ṅ，则测量引入的投影噪声等效磁场为：

δB_pn = 1/(γ·√(N·ṅ·T₂))

当探测光子数足够大（ṅ ≫ 1）时，投影噪声被充分平均，此时灵敏度由自旋噪声主导。反之，若探测光太弱，投影噪声会成为限制因素。两者共同构成总噪声：

δB_total = √(δB_sn² + δB_pn²)

三、灵敏度基本表达式与标准量子极限

把自旋噪声与投影噪声合并，并考虑测量带宽 Δf = 1/(2T₂)，得到原子磁力计灵敏度的基本表达式（单位 T/√Hz）：

δB = (1/(γ√(N·T₂)))·√(1 + 1/(ṅ·T₂))

在理想探测（ṅ·T₂ ≫ 1）下，第二项可忽略，灵敏度退化为：

δB_SQL = 1/(γ√(N·T₂))

这就是标准量子极限（Standard Quantum Limit, SQL）。它的物理内涵十分清晰：灵敏度受限于原子系综的投影噪声，而非技术噪声；要提高灵敏度，只能增加原子数 N、延长相干时间 T₂，或选择旋磁比 γ 更大的原子。值得注意的是，SQL 并非绝对下限——通过自旋压缩（spin squeezing）可让 ⟨ΔJ_z²⟩ 低于 Nħ²/4，从而突破 SQL，进入“海森堡极限”区域，但这需要额外的量子资源。

四、从公式到工程判断：一个示例估算

为把公式落到可操作的层面，设一个示例场景（以下参数为教学假设，非实测）：铷-87 原子，γ ≈ 7×10⁹ rad/(s·T)，原子数 N = 10¹²，相干时间 T₂ = 1 ms。代入 SQL 表达式：

δB_SQL = 1/(7×10⁹ × √(10¹² × 10⁻³))
       = 1/(7×10⁹ × √(10⁹))
       = 1/(7×10⁹ × 3.16×10⁴)
       ≈ 4.5×10⁻¹⁵ T/√Hz
       ≈ 4.5 fT/√Hz

这个量级与文献中高灵敏度 SERF 磁力计的报道范围一致，说明公式抓住了主要矛盾。若把 N 提高到 10¹³，SQL 改善约 √10 ≈ 3.2 倍，降至约 1.4 fT/√Hz；若 T₂ 从 1 ms 延长到 10 ms，同样改善约 3.2 倍。工程上，增加 N 受气室体积与加热均匀性限制，延长 T₂ 受自旋交换碰撞与磁场梯度限制，二者往往相互制约，这正是灵敏度优化的核心权衡。

五、标准量子极限的物理内涵

SQL 的本质是投影噪声与信号积累之间的竞争：测量越强（光子越多），投影噪声越小，但测量反作用会扰动自旋；测量越弱，扰动小但投影噪声大。在理想情况下，最优测量恰好使两者平衡，灵敏度达到 1/(γ√(N·T₂))。这也解释了为什么 SQL 与 √N 成反比——它本质上是 N 个独立量子投影的统计平均结果，而非任何技术缺陷。理解这一点，才能在评估实测灵敏度时正确区分“量子极限”与“工程极限”，避免把技术噪声误判为物理下限。