2. 原子磁力计的物理基础

原子磁力计之所以能测到远低于地磁百万分之一的弱磁场,根源在于它把“磁场”翻译成了“原子自旋的进动频率”。要理解这台仪器,必须先把四件事串起来:原子能级如何被磁场改写(塞曼效应)、自旋在磁场中如何转动(拉莫尔进动)、什么条件下会发生共振(磁共振条件),以及这些量之间如何换算成可测量的信号。本节以碱金属原子(如铷 87Rb、铯 133Cs)为对象,给出可计算的物理图像。

一、原子能级与磁矩:从零场到有场

碱金属原子最外层只有一个价电子,其基态由电子自旋 S 与核自旋 I 耦合形成总角动量 F。以 87Rb 为例,核自旋 I = 3/2,基态 2S1/2 分裂为 F = 1 与 F = 2 两个超精细能级。在没有外磁场时,同一 F 下的各磁子能级 mF 是简并的;一旦加上外磁场 B,简并解除,能级按 mF 分裂,这就是塞曼效应。

在弱场条件下(B 远小于超精细分裂对应的场,对 87Rb 约为 0.1 mT 量级),能级移动可写成线性形式:

ΔE = g_F · μ_B · m_F · B
g_F = g_J · [F(F+1) + J(J+1) − I(I+1)] / [2F(F+1)]

其中 μB ≈ 9.274×10−24 J/T 为玻尔磁子,gF 是朗德因子。对 87Rb 的 F = 2 态,gF ≈ 1/2;对 F = 1 态,gF ≈ −1/2。相邻磁子能级间隔为 ΔE = gFμBB,对应一个频率:

f_L = g_F · μ_B · B / h
h ≈ 6.626×10^(−34) J·s

代入数值可得一个非常实用的换算:对 gF = 1/2 的态,fL ≈ 7 Hz/nT(即 1 nT 对应约 7 Hz)。这个数字是原子磁力计灵敏度的“汇率”——测频精度越高,测磁分辨率越高。

二、拉莫尔进动:自旋在磁场中的旋转

把原子看成一个磁矩 μ,置于磁场 B 中,它受到的力矩为 τ = μ × B。由于角动量守恒,磁矩不会直接“倒向”磁场方向,而是绕 B 做进动,这就是拉莫尔进动。其角频率为:

ω_L = γ · B
γ = g_F · μ_B / ħ

对 87Rb 的 F = 2 态,γ ≈ 2π × 7 Hz/nT,即约 4.4×107 rad·s−1·T−1。进动方向由右手定则决定:若 B 沿 +z,则磁矩绕 +z 逆时针旋转(从 +z 向 −z 看)。这个方向约定在后续光路与线圈设计中必须保持一致,否则会把符号搞反。

进动本身不产生净磁矩变化,但一旦自旋被“倾斜”到与 B 成一定角度,它就会在横向平面(x-y 平面)产生一个以 fL 旋转的横向磁矩分量。这个旋转的横向磁矩正是原子磁力计真正“读出”的信号源。

三、磁共振条件:让自旋持续翻转

如果只让自旋自由进动,横向磁矩会因弛豫而衰减。要维持或操控它,需要施加一个与进动同步的射频(RF)磁场 B1,方向垂直于 B0。当 RF 频率等于拉莫尔频率时,发生磁共振:

f_RF = f_L = γ · B_0 / (2π)

在旋转坐标系中看,B1 变成一个静止的横向场,自旋绕它进动,从而被连续翻转。翻转角为:

θ = γ · B_1 · τ_p

其中 τp 是 RF 脉冲持续时间。若要让自旋从 +z 翻到 x-y 平面(90° 脉冲),需要 γB1τp = π/2。例如取 B1 = 10 nT,对 87Rb 的 γ ≈ 4.4×107 rad·s−1·T−1,则 τp ≈ 3.6 ms。这个量级说明:原子磁力计的 RF 脉冲通常在毫秒级,而不是微秒级。

共振条件不是“频率相等”这么简单:还必须满足 B1 ⊥ B0,且 RF 的相位与进动相位保持相干。相位失配会直接降低翻转效率,表现为信号幅度下降。

四、从物理量到可测信号:一张关系图

把上述链条整理成一张信号链关系图,可以清楚看到“磁场 → 能级分裂 → 进动频率 → 共振 → 光信号”的传递路径,以及每个环节的关键参数。

原子磁力计物理链:从外磁场到可测光信号 外磁场 B₀ 方向 +z 量级:nT ~ μT 塞曼能级分裂 ΔE = g_F·μ_B·m_F·B₀ m_F 简并解除 拉莫尔进动 ω_L = γ·B₀ f_L ≈ 7 Hz/nT 磁共振条件 f_RF = f_L B₁ ⊥ B₀ 自旋进动与横向磁矩 z (B₀) x y 进动方向 横向磁矩 M_xy • 自旋绕 B₀ 进动,横向分量以 f_L 旋转 • 弛豫时间 T₂ 决定横向信号衰减速度 • RF 脉冲使自旋翻转,产生可观测的宏观磁矩 光与原子相互作用 泵浦光:制备自旋极化 探测光:读出偏振旋转 光电探测器 光强/偏振变化 → 电信号 锁相放大提取 f_L 分量 频率 → 磁场 B = 2π·f_L / γ 分辨率 ≈ 1/γ · Δf 共振驱动

图中最值得注意的两条“换算线”:一是从 B0 到 fL 的 γ 换算,二是从 fL 到 B 的反换算。实际仪器中,我们并不直接测 B,而是测 fL,再用 γ 反推。因此 γ 的标定精度直接决定磁场测量精度。

五、一个可复算的教学示例

假设我们要设计一台 87Rb 原子磁力计,目标测量范围 20 nT ~ 100 nT,要求频率分辨率 0.1 Hz。按 fL ≈ 7 Hz/nT 计算:

参数数值说明
目标磁场 B50 nT示例工作点
拉莫尔频率 f_L350 Hz7 Hz/nT × 50 nT
RF 频率 f_RF350 Hz共振条件 f_RF = f_L
B₁ 幅度10 nT典型 RF 驱动场
90° 脉冲宽度 τ_p≈ 3.6 msγB₁τ_p = π/2
频率分辨率 Δf0.1 Hz锁相放大器带宽决定
磁场分辨率 ΔB≈ 0.014 nTΔf / 7 Hz/nT

这个示例说明:在 50 nT 的工作点,拉莫尔频率只有 350 Hz,属于音频范围,对电子学带宽要求不高;但要达到 0.014 nT 的磁场分辨率,需要把频率测量稳定到 0.1 Hz 以内,这对光源稳定度、气室温度稳定度和磁屏蔽都提出了要求。

常见误区:把 fL 与 RF 驱动频率混为一谈。自由进动时没有 RF,fL 由 B0 决定;只有施加 RF 且频率匹配时,才发生共振翻转。两者在数值上相等,但物理角色不同。

六、弛豫与线宽:为什么信号会“变钝”

自旋进动不是永久的。横向磁矩会因自旋交换碰撞、与气室壁碰撞、磁场梯度等原因衰减,时间常数为 T2。共振线宽与 T2 的关系为:

Δf_FWHM = 1 / (π · T₂)

若 T2 = 10 ms,则线宽约 32 Hz。线宽越窄,频率分辨能力越强,但信号幅度也越低。实际设计中需要在“高灵敏度”与“高信噪比”之间折中:通常通过增大气室尺寸、优化抗弛豫涂层、降低磁场梯度来延长 T2,同时用锁相放大和调制解调技术把微弱信号从噪声中提取出来。