里德堡相互作用与 Förster 共振

两个里德堡原子之间的相互作用就是两个强偶极子之间的静电多极能：共振时表现为 $1/R^3$ 的偶极-偶极”翻转-翻转”交换，非共振时退化为 $1/R^6$ 的范德瓦尔斯能移；连接这两种极限的旋钮是 Förster 缺陷——由 $n^{-3}$ 能级间隔决定的邻近对态能量失配，它可以用直流 Stark 场或微波 dressing 人为调零、调负或调正，从而把相互作用”工程化”。

物理图像

处于基态的中性原子之间几乎没有相互作用；但两个里德堡原子各自携带一个 $\sim n^{2}e{a}_0$ 的巨大电偶极子，它们之间的静电能即使在几微米、十几微米的间距上仍达 MHz–GHz 量级。直观图像有两条 Henriet 2020：共振交换——若两个对态（如 $|nS,n^{\prime }P⟩$ 与 $|n^{\prime }P,nS⟩$）能量简并，偶极-偶极作用使两个原子以”翻转-翻转”（flip-flop）方式相干交换激发，相互作用能 $\propto 1/R^3$；范德瓦尔斯——若最近的对态失谐很大，二阶微扰给出 $\propto 1/R^6$ 的纯能移。作为经验规则，碱金属 $nS$ 态之间的范德瓦尔斯相互作用近似各向同性且为排斥（$C_6<0$，按”能量移正为排斥”的约定），而 $nS$–$nP$ 之间则呈 $\propto 1-3{\cos }^2\theta$ 的简单偶极形式 Morgado 2021。

相互作用的强度与各向异性可以在很大的范围内设计：选态（$nS$、$nD$、高 $\ell$ 态）、直流场 Stark 调谐、微波 dressing 都能改变有效相互作用——这把里德堡阵列变成一个可编程的相互作用资源，支撑里德堡阻塞门与模拟优化的整个相图 Morgado 2021。

里德堡-里德堡相互作用物理起源示意：一对偶极耦合态 |a>,|c> 与 |b>,|d>、Förster 缺陷 ΔF，以及共振 1/R³ 与范德瓦尔斯 1/R⁶ 两种极限（图源：Morgado 2021 Fig. 6）

理论模型

多极展开与偶极-偶极哈密顿量

取两原子间距 $R$ 大于 Le Roy 半径 $R_{LR}=2\sqrt{⟨r_1^2⟩}+2\sqrt{⟨r_2^2⟩}$（电子云不重叠、随 $n^2$ 增大），相互作用可用多极展开 Weber 2017

$${\hat{H}}_{\mathrm{int}}(R)=\sum _{{\kappa }_1,{\kappa }_2\ge 1}\frac{{\hat{V}}_{{\kappa }_1{\kappa }_2}}{4\pi {ϵ}_0{R}^{{\kappa }_1+{\kappa }_2+1}},$$

展开阶 $\varrho ={\kappa }_1+{\kappa }_2+1\ge 3$，最低阶即偶极-偶极作用（$\varrho =3$），次高阶为偶极-四极（$\varrho =4$）等；截断阶数随间距减小而提高。偶极-偶极项的矢量形式为 Morgado 2021

$${\hat{V}}_{dd}=-\frac{1}{4\pi {ϵ}_0{R}^3}\left[{\hat{\mathbf{\mu }}}_j\cdot {\hat{\mathbf{\mu }}}_k-3({\hat{\mathbf{\mu }}}_j\cdot \hat{\mathbf{n}})({\hat{\mathbf{\mu }}}_k\cdot \hat{\mathbf{n}})\right],$$

其中 $\hat{\mathbf{n}}$ 是沿两原子连线的单位矢量。其矩阵元由径向偶极矩（$\propto n^2$）与角度因子（Clebsch–Gordan 系数与球谐函数）共同决定。

两通道模型：Förster 缺陷与两种极限

保留一对近共振对态 $|ab⟩$、$|cd⟩$，两原子哈密顿量为 Morgado 2021

$$\hat{H}=\frac{{\mu }_{ac}{\mu }_{bd}}{R_{j,k}^3}\left(|ab⟩⟨cd|+h.c.\right)+{\Delta }_F|cd⟩⟨cd|,$$

其中 ${\mu }_{\alpha \beta }=⟨\alpha |\hat{\mu }|\beta ⟩/\sqrt{4\pi {ϵ}_0}$，而

$${\Delta }_F=(E_c+E_d)-(E_a+E_b)$$

称为 Förster 缺陷（Förster defect）。对角化这个二能级问题，两条绝热势能曲线为

$$\Delta {\phi }_{\pm }(R)=\frac{{\Delta }_F}{2}\pm \sqrt{{\left(\frac{{\Delta }_F}{2}\right)}^2+\frac{{\mu }_{ac}^2{\mu }_{bd}^2}{R^6}}.$$

两种极限一目了然：

• 共振区（${\Delta }_F\approx 0$）：两支在 $R$ 处劈裂 $2{\mu }_{ac}{\mu }_{bd}/R^3$，布居以频率 $\propto 1/R^3$ 在两对态间相干往返（翻转-翻转交换），相互作用强度 $V_{j,k}=2{\mu }_{ac}{\mu }_{bd}/R_{j,k}^3$；
• 非共振区（$|{\Delta }_F|\gg |{\mu }_{ac}{\mu }_{bd}|/R^3$）：二阶微扰给出范德瓦尔斯能移

$$V_{j,k}=-\frac{C_6}{R_{j,k}^6},C_6=\frac{|{\mu }_{ac}{|}^2|{\mu }_{bd}{|}^2}{{\Delta }_F}.$$

两区的分界由交叉距离 $R_c$（定义 $|{\Delta }_F|=C_3\sqrt{D_{\phi }}/R_c^3$，$D_{\phi }$ 为角度本征值）刻画；一般情形（多通道、含角度）下 ${\hat{V}}_{dd}{\hat{V}}_{dd}^{†}$ 的本征值 $C_3^2{D}_{\phi }/R^6$ 给出一族势能曲线 $\Delta \phi (R)=\delta /2+\sqrt{{\delta }^2/4+C_3^2{D}_{\phi }/R^6}$，其中 $\delta$ 为广义 Förster 缺陷 Saffman 2010。

$n$ 标度

把标度律代回：$\mu \propto n^{2}e{a}_0$ 给出共振偶极系数 $C_3\propto n^4$；而 $C_6={\mu }^4/{\Delta }_F$ 中 ${\Delta }_F\propto n^{-3}$（或等价地 $\alpha \propto n^7$、$\Delta E\propto n^{-3}$），故

$$C_6\propto n^{11},$$

这正是里德堡门倾向于使用高 $n$ 态的原因 Adams 2019 Saffman 2010。注意 $C_6$ 的符号与大小由 ${\Delta }_F$ 决定：$\delta \to 0$ 处 $C_6$ 发散并变号，势能曲线在 $\delta \approx 0$ 的态对附近出现共振增强或” Förster 零”（近似平坦的 $C_6\approx 0$ 通道）Saffman 2010 Young 2021。

参数与量级

量	典型值	来源
偶极-偶极标度	$C_3\propto n^4$，$V\propto 1/R^3$	Adams 2019
范德瓦尔斯标度	$C_6\propto n^{11}$，$V\propto 1/R^6$	Saffman 2010
Rb $60s+60s$ Förster 缺陷 $\delta$	$1.7$ GHz	Saffman 2010
Rb $43d_{5/2}+43d_{5/2}\to 41f$ 缺陷	$-6.0/-8.3$ MHz（近” Förster 零”）	Saffman 2010
Rb $59d$ 对零场缺陷	$8.69$ MHz，$34.3$ mV/cm 调至共振	Weber 2017
偶极交换振荡	$9.2$ MHz（$R=9.1\mu$m，$\theta =0^{\circ }$）	Weber 2017
${}^{88}$Sr $61s{}^3{S}_1$ 的 $C_6/h$	$\approx -150\mathrm{GHz}\cdot \mu {\mathrm{m}}^6$	Kazemi 2025
Stark 场稳定度要求	$<1$ MHz 移动需 $\sim 0.01(100/n{)}^{7/2}$ V/cm	Saffman 2010

方法与实现

Stark 调谐 Förster 共振。 由于不同态极化率不同，直流电场可以调节 ${\Delta }_F$。标志性实验：两光镊中的 Rb 原子被制备在 $|59d_{3/2},m_j=3/2⟩$，零场下与 $|61p_{1/2};57f_{5/2}⟩$ 对态失配 $8.69$ MHz，加 $34.3$ mV/cm 电场调至简并后，$R=9.1\mu$m、$\theta =0^{\circ }$ 处观察到 $9.2$ MHz 的无阻尼偶极交换振荡；$\theta =14^{\circ }$ 时多通道耦合导致退相干——完整的角度花样与理论计算定量吻合 Weber 2017。

基于 Förster 通道的门。 在 Stark 调谐的 Förster 共振上做双绝热快速通道（double adiabatic passage）可以积累确定性几何相位，理论方案以 Cs $90S+96S\to 90P+95P$ 共振为例给出 CZ/CNOT 门，对原子位置涨落（$\sim 1/R^3$ 的敏感性）显著更鲁棒；射频辅助的 Förster 共振可进一步扩展通道 Beterov 2016 Beterov 2018。这与ARP 门同属绝热通道门家族。

态选择与”零点”工程。 $s$ 态不存在 Förster 零、范德瓦尔斯近似各向同性且排斥，是阻塞实验的默认选择；$d$ 态则常出现 $C_6\approx 0$ 的”零”通道与强各向异性，可用于需要在某些方向”关掉”相互作用的方案 Saffman 2010。把不同里德堡态的相互作用强度比放大（如控制–靶强、控制–控制弱），可实现非对称阻塞与多比特门，通常借助微波 dressing 混合邻近里德堡态完成 Young 2021——这正是多比特阻塞门与奇偶相位门的相互作用侧基础。

直流场敏感性的抑制。 $\alpha \propto n^7$ 的巨大极化率使里德堡门对杂散静电场极其敏感（补偿到 MHz 级需要 $0.01(100/n{)}^{7/2}$ V/cm 的场稳定度 Saffman 2010）。抑制手段包括：非共振微波 dressing 抵消态间偶极矩差——实验上用 $38.465$ GHz 场在 $\sim 1$ V/cm 直流场中消除了 Rb $49s/48s$ 的静态偶极矩差 Jones 2013；双频 dressing 同时压低两个态的场敏感性，对 Cs $90S_{1/2}$ 与 $90P_{3/2}$ 分别给出 $95$ 倍与 $1600$ 倍的抑制 Booth 2018。这类 dressing 属于量子最优控制框架中”用控制通道改造漂移哈密顿量”的典型例子。

计算工具。 精确的相互作用势（含高阶多极、外场、自旋-轨道耦合）可用开源程序 pairinteraction 与 ARC 直接对角化计算 Weber 2017 Robertson 2021；它们是解读与设计上述实验的标准工具。

实验特征与测量

Förster 共振的最直接信号是相干交换振荡：将一对原子激发到调谐共振的里德堡态，随保持时间测量布居在 $|dd⟩$ 与交换对态之间的振荡，振荡频率即 $2\mu \mu /R^3$；改变原子间夹角 $\theta$ 可以测绘相互作用的偶极角度花样（如 $1-3{\cos }^2\theta$）Weber 2017。范德瓦尔斯区则通过阻塞抑制、 blockade 半径测量或阵列上的 $C_6$ 标定来表征，这些属于里德堡阻塞与里德堡模拟优化词条的范畴。对多原子奇偶相位门，$C_6<0$ 的排斥性范德瓦尔斯相互作用（如 ${}^{88}$Sr 的 $\approx -150\mathrm{GHz}\cdot \mu {\mathrm{m}}^6$，间距 $\ge 3\mu$m）直接给出可用的多体相位速率 Kazemi 2025。

与其他概念的关系

• 相互作用的第一个应用是里德堡阻塞：$V(R)$ 与阻塞半径 $R_b$ 的定义都源自 $C_6/R^6$ 或 $C_3/R^3$。
• 阻塞机制上的两比特门见阻塞 CZ 门；跨 Förster 共振的绝热通道方案见ARP 门。
• 相互作用强度比的工程（非对称阻塞）支撑多比特阻塞门、Toffoli/CCZ 门与奇偶相位门。
• $C_6$ 的符号、各向异性与 $R_b/a$ 之比决定里德堡模拟优化的哈密顿量与相图。
• 杂散电场引起的相互作用噪声属于误差预算与真实建模；其抑制（dressing）依赖量子最优控制框架。
• 相互作用载体的微观参数（量子亏损、偶极矩、寿命）见里德堡原子与里德堡态。

参考文献

• M. Saffman, T. G. Walker, K. Mølmer, “Quantum information with Rydberg atoms”, Rev. Mod. Phys. 82, 2313 (2010). NA001 · arXiv:0909.4777 · DOI 10.1103/RevModPhys.82.2313
• C. S. Adams, J. D. Pritchard, J. P. Shaffer, “Rydberg atom quantum technologies”, J. Phys. B: At. Mol. Opt. Phys. (2019). NA003 · arXiv:1907.09231 · DOI 10.1088/1361-6455/ab52ef
• M. Morgado, S. Whitlock, “Quantum simulation and computing with Rydberg-interacting qubits”, AVS Quantum Sci. 3, 023501 (2021). NA005 · arXiv:2011.03031 · DOI 10.1116/5.0036562
• L. Henriet et al., “Quantum computing with neutral atoms”, Quantum 4, 327 (2020). NA004 · arXiv:2006.12326 · DOI 10.22331/q-2020-09-21-327
• S. Weber, C. Tresp, H. Menke, A. Urvoy, O. Firstenberg, H. P. Büchler, S. Hofferberth, “Calculation of Rydberg interaction potentials”, J. Phys. B: At. Mol. Opt. Phys. 50, 133001 (2017). Weber 2017 · arXiv:1612.08053 · DOI 10.1088/1361-6455/aa743a
• I. I. Beterov, M. Saffman, E. A. Yakshina, D. B. Tretyakov, V. M. Entin, S. Bergamini, E. A. Kuznetsova, I. I. Ryabtsev, “Two-qubit gates using adiabatic passage of the Stark-tuned Förster resonances in Rydberg atoms”, Phys. Rev. A 94, 062307 (2016). Beterov 2016 · DOI 10.1103/PhysRevA.94.062307
• I. I. Beterov et al., “Adiabatic passage of radio-frequency-assisted Förster resonances in Rydberg atoms for two-qubit gates and the generation of Bell states”, Phys. Rev. A 97, 032701 (2018). Beterov 2018 · DOI 10.1103/PhysRevA.97.032701
• J. T. Young, P. Bienias, R. Belyansky, A. M. Kaufman, A. V. Gorshkov, “Asymmetric Blockade and Multiqubit Gates via Dipole-Dipole Interactions”, Phys. Rev. Lett. 127, 120501 (2021). NA069 · arXiv:2006.02486 · DOI 10.1103/PhysRevLett.127.120501
• L. A. Jones, J. D. Carter, J. D. D. Martin, “Rydberg atoms with a reduced sensitivity to dc and low-frequency electric fields”, Phys. Rev. A 87, 023423 (2013). Jones 2013 · DOI 10.1103/PhysRevA.87.023423
• D. W. Booth, J. Isaacs, M. Saffman, “Reducing the sensitivity of Rydberg atoms to dc electric fields using two-frequency ac field dressing”, Phys. Rev. A 97, 012515 (2018). Booth 2018 · DOI 10.1103/PhysRevA.97.012515
• J. Kazemi, M. Schuler, C. Ertler, “Multiqubit parity gates for Rydberg atoms in various configurations”, Phys. Rev. Research (2025). NA084 · DOI 10.1103/56qk-rmsz