2. Rydberg原子物理基础

里德堡原子是量子门控制工程中最核心的物理载体。理解它的第一步,不是记住某个公式,而是建立一幅清晰的“能量—波函数—寿命”图像:一个价电子被激发到主量子数 n 很大的轨道上,远离原子实,感受到的屏蔽越来越强,能级间距越来越小,轨道半径按 n² 增长,寿命也随 n 显著变长。本节围绕四个关键量展开:主量子数 n、量子亏损 δ、里德堡态波函数 ψ、以及寿命与黑体辐射效应。

一、主量子数 n:里德堡态的“尺度旋钮”

主量子数 n 决定里德堡态的基本尺度。对碱金属原子(如 Rb、Cs)的价电子,能量可近似写成

E(n, l, j) = -R_y / (n - δ_{l,j})²

其中 R_y 是里德堡常数,δ 是量子亏损。n 越大,能级越靠近电离阈,相邻能级间距按 1/n³ 缩小。例如 Rb 的 nS 态,n=50 与 n=51 的能量差约几十 GHz,而 n=100 时已降到约几 GHz。这个“能级变密”直接决定了微波操控、态分辨和串扰控制的难度。

轨道半径的标度关系为 r ∝ n²。n=50 时,里德堡电子平均半径约 0.1 μm;n=100 时可达约 0.4 μm。两个这样的原子靠近时,电子云显著重叠,产生强范德华或偶极相互作用,这正是里德堡阻塞与两比特门的物理来源。

二、量子亏损 δ:原子实穿透的修正

量子亏损 δ 描述价电子穿透原子实、感受到未被完全屏蔽的核电荷的程度。它依赖轨道角动量 l:l 越小,电子越容易靠近原子实,δ 越大;l 越大,电子被离心势挡在外面,δ 越小。对 Rb:

态δ(近似)物理含义
nS≈ 3.13穿透最深,能级下移最多
nP≈ 2.64穿透较强
nD≈ 1.35穿透减弱
nF≈ 0.02几乎不穿透,接近氢原子

工程上,δ 不是“小修正”,而是决定有效主量子数 n* = n - δ 的关键参数。计算相互作用强度、阻塞半径、微波耦合矩阵元时,必须用 n* 而不是裸 n。例如 n=60 的 Rb nS 态,n* ≈ 56.87,能量与半径都要按 n* 重新标度。

三、里德堡态波函数:径向节点与角向分布

里德堡态波函数 ψ(r, θ, φ) = R_{n*l}(r) Y_{lm}(θ, φ)。径向部分在原子实外近似为库仑波函数,具有 n* - l - 1 个径向节点,振荡周期随 r 增大而变长。角向部分由球谐函数描述,l 决定电子云形状:s 态球对称,p 态哑铃形,d 态四叶形。

波函数的“尾巴”极其重要:在 r 很大处,R_{n*l}(r) 仍保持可观幅度,衰减长度约 n*² a₀。两个里德堡原子相距 R 时,相互作用矩阵元正比于 ψ(r) 在对方位置的重叠积分。因此,波函数的空间延展直接转化为可调的相互作用强度。

里德堡态能量—波函数—寿命关系图(以 Rb nS 为例) 能级与 n* E n=30, n*≈26.9 n=50, n*≈46.9 n=70, n*≈66.9 n=90, n*≈86.9 n=110, n*≈106.9 间距 ∝ 1/n*³,越往上越密 径向波函数 R(r) 与节点 r / a₀ 节点 节点 节点 节点数 = n* - l - 1;尾部衰减长度 ∝ n*² 寿命与黑体辐射 自发辐射寿命 τ₀ ∝ n*³ n=50: τ₀ ≈ 100 μs n=100: τ₀ ≈ 800 μs 300 K 黑体辐射主导 τ_BBR 可缩短至 ~10 μs 低温 4 K:BBR 抑制 τ 接近 τ₀,适合长门操作 室温:需缩短脉冲或主动补偿 n ↑ → n* ↑ → 半径 ∝ n*² ↑ → 相互作用 ↑ → 阻塞半径 ↑;同时 τ₀ ∝ n*³ ↑,但 BBR 损耗也 ↑

四、寿命与黑体辐射效应

里德堡态寿命由自发辐射和黑体辐射共同决定。自发辐射寿命近似 τ₀ ∝ n*³。以 Rb nS 态为例,n=50 时 τ₀ 约 100 μs,n=100 时约 800 μs。但室温 300 K 下,黑体辐射光子会诱导受激跃迁,把寿命压缩到几十微秒甚至更短。黑体辐射跃迁率随 n 增大而迅速上升,因为相邻里德堡态之间的偶极矩阵元 ∝ n*²,而黑体光子谱在微波/远红外波段有足够密度。

工程上必须区分两种场景:

  • 室温环境:黑体辐射主导,n=70 的 Rb 里德堡态有效寿命可能只有约 10–20 μs。此时门操作时间必须远小于该值,通常控制在 1 μs 以内,或采用低温环境。
  • 低温环境(4 K):黑体辐射被指数抑制,寿命接近自发辐射极限,适合需要长相干时间的量子门或存储实验。
在估算门保真度时,可先用 τ_eff ≈ (1/τ₀ + 1/τ_BBR)⁻¹ 得到有效寿命,再与门时间 T_gate 比较。若 T_gate / τ_eff > 0.1,就需要把寿命损耗纳入保真度预算。

五、参数选择与工程判断

实际选择 n 时,需要同时权衡三件事:相互作用强度、寿命、以及操控可行性。n 太小,相互作用弱,阻塞半径小,两比特门速度慢;n 太大,寿命被黑体辐射严重压缩,且能级过密导致微波串扰和态分辨困难。对 Rb 两比特门,常见工作区间是 n = 50–80,对应阻塞半径约 5–15 μm,有效寿命约 10–100 μs(视温度而定)。

量子亏损 δ 的选取则与态类型有关:nS 态 δ 大,能级下移明显,适合用双光子激发;nP 态常用于单光子激发和微波耦合;nD 态量子亏损小,更接近氢原子,适合研究长程偶极相互作用。实际实验中,先用光谱测量确定目标态的 δ 和 n*,再据此计算相互作用强度与阻塞半径,最后结合寿命预算确定门时间与激光/微波参数。

注意:量子亏损 δ 随 l 和 j 变化,不同精细态的能量差可能达到 GHz 量级。若忽略精细结构,计算出的共振频率会偏离实际,导致激发效率下降。工程上应使用实验标定的 δ 值,而不是氢原子近似。