📚 课程目录第3章灵敏度计算的核心物理量

3. 灵敏度计算的核心物理量

原子磁力计的灵敏度不是单一器件指标，而是若干物理量在同一个噪声模型下共同作用的结果。要真正读懂灵敏度公式、判断哪一项参数是当前系统的瓶颈，就必须先弄清每个核心物理量的定义、量纲、典型取值，以及它在公式中究竟以“信号增益”还是“噪声源”的身份出现。本节围绕旋磁比 γ、弛豫时间 τ（或弛豫率 Γ）、原子数密度 n、极化率 P 和检测带宽 Δf 五个量展开，并给出把它们串起来的灵敏度表达式与一个可复算的示例。

一、五个物理量的定义与物理意义

旋磁比 γ 描述原子磁矩与角动量之间的比例关系，单位是 rad·s⁻¹·T⁻¹（工程上也常用 Hz/T）。对铷-87 的基态，γ ≈ 2π × 7 Hz/μT，即约 7 Hz/nT。它的物理意义是：外磁场 B 引起的拉莫尔进动角频率为 ω_L = γB。因此 γ 决定了“同样的磁场变化能换来多大的频率/相位响应”，是信号链的转换系数。γ 由原子种类决定，碱金属之间差异明显（K、Rb、Cs 各不相同），选工作原子时它就是一个先天参数。

弛豫时间 τ（等价地用弛豫率 Γ = 1/τ 或线宽 Δν = Γ/2π）描述宏观磁化矢量偏离平衡后恢复的速率，来源包括自旋交换碰撞、自旋-破坏碰撞（与缓冲气体、壁面）、磁场梯度引起的退相干等。τ 越长，磁化矢量“记得”越久，对缓变磁场的响应越灵敏，但同时系统带宽变窄、对快速变化磁场的跟随能力下降。τ 是灵敏度与带宽之间最核心的折中参数。

原子数密度 n 是单位体积内参与工作的碱金属原子数，单位 m⁻³ 或 cm⁻³。它决定可被极化的“自旋资源”总量：n 越大，同样极化率下能提供的总磁矩越大，信噪比越高。但 n 增大会加剧自旋交换碰撞与吸收，可能反过来缩短 τ，因此并非越大越好。

极化率 P 是原子系综沿量化轴的平均自旋极化程度，无量纲，取值 0～1。P 直接决定宏观磁化强度 M = n·μ·P（μ 为单个原子磁矩），是信号幅度的线性因子。P 由光抽运速率与弛豫速率之比决定，稳态下近似 P = R_pump/(R_pump + Γ)。

检测带宽 Δf 是测量系统有效响应的频率范围，单位 Hz。它既由弛豫决定的物理带宽（约 Γ/2π 量级）限制，也由探测光功率、电子学滤波和采样率限制。带宽越宽，纳入的白噪声越多，单点噪声谱密度按 √Δf 增长；因此灵敏度通常以 T/√Hz 的噪声谱密度表述，以便与带宽解耦。

二、它们在灵敏度公式中的角色

在光泵浦磁力计的典型噪声模型中，极限灵敏度可写成磁化矢量在单位带宽内的噪声等效磁场。对以弛豫为主导的系综，一个常用的近似表达为：

δB ≈ 1 / (γ · P · √(n · V · τ · Δf))   （量级关系，系数依探测方式而定）

其中：γ  —— 旋磁比，信号转换系数（分母，越大越好）
      P  —— 极化率，信号线性因子（分母，越大越好）
      n  —— 原子数密度，自旋资源（分母根号内）
      V  —— 作用体积
      τ  —— 弛豫时间，相干积累时间（分母根号内）
      Δf —— 检测带宽，噪声积分窗口（分母根号内，越大越差）

从公式结构可以读出三条关键结论：第一，γ 和 P 是“线性增益”，提升它们对灵敏度是线性收益；第二，n、V、τ 出现在根号内，属于“统计增益”，翻倍只带来约 1.41 倍改善；第三，Δf 在分母根号内，意味着带宽越宽灵敏度数值越差，这正是必须用 T/√Hz 归一化的原因。理解这一点，就能在调参时判断该优先动哪个量。

三、示例：用一组假设参数估算灵敏度并定位瓶颈

下面给出一组教学假设参数（非实测），用于演示如何把五个量代入并判断瓶颈。设工作原子为 Rb-87，γ = 2π × 7 Hz/μT；作用体积 V = 0.1 cm³；原子数密度 n = 1×10¹³ cm⁻³；极化率 P = 0.5；弛豫时间 τ = 10 ms；检测带宽 Δf = 100 Hz。先算弛豫决定的物理线宽 Δν = 1/(2πτ) ≈ 15.9 Hz，说明 100 Hz 的检测带宽已超出物理带宽，实际有效带宽受弛豫限制，此时应把 Δf 取为约 16 Hz 而非 100 Hz，否则会高估噪声。

按量级关系 δB ≈ 1/(γ·P·√(n·V·τ·Δf)) 代入（统一 SI 单位：n = 1×10¹⁹ m⁻³，V = 1×10⁻⁷ m³，τ = 0.01 s，Δf = 16 Hz），可得 n·V·τ·Δf ≈ 1×10¹⁹ × 1×10⁻⁷ × 0.01 × 16 = 1.6×10¹⁰，开方约 1.26×10⁵；γ·P ≈ 4.4×10⁷ × 0.5 ≈ 2.2×10⁷；两者相乘约 2.8×10¹²，取倒数约 3.6×10⁻¹³ T，即约 0.36 pT 量级。这一结果只用于说明数量级与参数敏感度，真实系统还需计入探测光散粒噪声、电子学噪声等。

敏感性分析更有价值：若把 P 从 0.5 提到 0.8，灵敏度改善约 1.6 倍；若把 τ 从 10 ms 提到 40 ms，仅改善约 2 倍（根号效应），但带宽同时收窄 4 倍，对快变场测量不利；若把 n 提高一个数量级，理论上改善约 3.16 倍，但自旋交换碰撞会压缩 τ，实际收益往往被抵消。由此可以形成调参判断：先确认 γ 选型是否合适、P 是否已接近抽运-弛豫平衡上限，再在 τ 与 Δf 之间按被测磁场的时间尺度做折中，最后才考虑提高 n 或 V。

判断瓶颈的实用方法：分别单独改变一个参数（如把抽运光功率提高使 P 上升），观察噪声谱密度是否随之下降。若某参数变化后灵敏度几乎不动，说明它已不是限制项，应转向下一个量。

注意单位一致性：γ 用 rad·s⁻¹·T⁻¹、n 用 m⁻³、τ 用 s、Δf 用 Hz 时结果才是 T。工程中常混用 cm⁻³、μT、Hz/T，换算错误是灵敏度算错的首要原因。

把五个量放在同一张关系图上（见上方 SVG），可以清楚看到：γ 与 P 位于信号增益链的前端，是线性杠杆；n、V、τ 是统计杠杆；Δf 是噪声窗口，且与 τ 存在此消彼长的折中。掌握这一结构，后续无论是读灵敏度公式、做参数标定，还是排查“灵敏度达不到指标”的问题，都能快速定位到具体物理量，而不是笼统地归因于“系统噪声大”。