中性原子 QRAM

当 QRAM 的问题从"能不能路由"转向"能不能容错",中性原子凭借原子重排与擦除探测能力,成为资源态 QRAM 蓝图选择的宿主平台。

三大平台禀赋

其一,相干动态重排:光镊可在保持量子态的前提下重排原子布图,使”工厂区”与”查询区”之间的时空调度成为可编程资源。其二,擦除转换:Scholl 等 2023(Nature 622, 273)证明中性原子处理器可把双比特门错误转化为可检测的擦除/相位错误,天然对接 Z 偏置的非对称纠错。其三,逻辑块搬运:Bluvstein 等 2024(Nature 626, 58)的逻辑原子处理器展示了逻辑块之间以原子搬运直接接口的可行性。三者叠加,恰好覆盖 QRAM 从制备、纠错到互连的全部关键环节。

资源态 QRAM 蓝图

Cesa 等 2026(Nat. Commun.)是迄今最完整的中性原子 QRAM 方案。要点:查询的全部在线操作被 Clifford 化,非 Clifford 成本整体转移进可离线批量制备的资源态。“QRAM 工厂”先制备 $O(N)$ 个 Bell 对,经光镊相干重排完成布图重构,再以 $O(N\log N)$ 个 Toffoli(电路深度 $2\log N$)组装资源态;工厂周期 $T_{\Phi }\lesssim 3\sqrt{6}T{N}^{1/4}+2\tau \log N$,其中 $T\approx 33\mu \mathrm{s}$ 为最小陷阱间距搬运时间、$\tau \approx 500\mu \mathrm{s}$ 为中途测量时间。对 1 kB 存储,查询时间约 13 ms、工厂速率约 0.1 kHz,且查询与下一份资源态的制备可流水并行;重排耗时在 $N\lesssim 2^{24}$(约 2 MB)内可忽略。

噪声与容错接口

数值模拟给出查询不保真度 $\sim {\log }^{\alpha }N$:去极化噪声下 $\alpha \approx 3.8$,Z 偏置(经擦除转换)下 $\alpha \approx 1.9$;门错误率 $\epsilon \sim 10^{-4}$ 时 $N\gtrsim 512$ 保真度可过半,Z 偏置情形 $N\approx 512$ 即可达 99%。由于在线全为 Clifford 操作,查询可在表面码/CSS 码上横向容错执行,并可自然融入测量式纠错。综述的判断是:中性原子正是凭此在”容错接口”这一最艰难的维度上后来居上——这也是资源态 QRAM 选择该平台的原因;该思想亦可迁移至囚禁离子、超导、光子等具备搬运能力的平台。中期规模($2^{10}$–$2^{24}$,重排开销可忽略、工厂与查询流水)被视为该蓝图的甜蜜区。

以上三节是平台定位;以下四节推导回答一个更物理的问题:这套蓝图凭什么落在中性原子上——门错误从哪里来、温度与位置涨落如何进入预算、$2N-1$ 个路由器如何铺进二维阵列。除注明外,器件参数取 ${}^{87}\mathrm{Rb}$ 典型值,推导为教科书级、数值属本词条估算。

推导一:光镊囚禁——原子比特的运动学基座

QRAM 对存储阵列/路由节点的最基本要求:每个原子长时间保持量子态、位置确定到远小于相互作用与寻址尺度。这条要求由光镊(optical tweezer,即光偶极阱 optical dipole trap)承担。本节把运动学参数从器件标量(阱深、腰斑)一路推到 QRAM 关心的量(波包尺寸、速度弥散、加热率、搬运下限)。

偶极阱:从光强到谐振子频率

红失谐激光对基态原子的偶极势(取极化率实部):

$U(\mathbf{r})=-\frac{\mathrm{Re}\alpha (\omega )}{2{\epsilon }_0}I(\mathbf{r}),I(\mathbf{r})=I_0{e}^{-2r^2/w_0^2},$

高斯束在焦点附近展开为抛物线 $U(\mathbf{r})\approx -U_0+2U_0{r}^2/w_0^2$,与 $\frac{1}{2}m{\omega }^2{r}^2$ 逐项比对,得横向陷阱频率

${\omega }_{\perp }=\frac{2}{w_0}\sqrt{\frac{U_0}{m}}.$

代入 $U_0/k_B=1\mathrm{mK}$、$m=1.443\times 10^{-25}\mathrm{kg}$、$w_0=0.7\mu \mathrm{m}$:$\sqrt{U_0/m}=0.309m/s$,故 ${\omega }_{\perp }/2\pi \approx 141\mathrm{kHz}$。轴向曲率由瑞利区给出:$z_R=\pi w_0^2/{\lambda }_t\approx 1.8\mu \mathrm{m}$(${\lambda }_t=850\mathrm{nm}$),$U(z)\approx -U_0(1-z^2/z_R^2)$ 给出 ${\omega }_z=\sqrt{2U_0/(m{z}_R^2)}\approx 2\pi \times 38\mathrm{kHz}$。关键量纲关系是 $\hslash {\omega }_{\perp }/k_B\approx 6.8\mu \mathrm{K}$:10 μK 的原子平均只占 $\bar{n}=T/(\hslash {\omega }_{\perp }/k_B)-\frac{1}{2}\approx 1$ 个振动量子——接近基态,这是进入兰姆–迪克区的前提。

Lamb–Dicke 区与热分布

谐振子基态尺寸、热速度弥散与热波包尺寸:

$$x_0=\sqrt{\frac{\hslash }{2m{\omega }_{\perp }}}\approx 20\mathrm{nm},{\sigma }_v=\sqrt{\frac{k_{B}T}{m}}\approx 3.1cm/s(T=10\mu \mathrm{K}),x_{\mathrm{rms}}=\frac{{\sigma }_v}{{\omega }_{\perp }}\approx 35\mathrm{nm}.$$

兰姆–迪克参数(Lamb–Dicke parameter)$\eta =k_L{x}_0=2\pi x_0/780\mathrm{nm}\approx 0.16\ll 1$:激光”看”到的是一个位置固定的质点,反冲与多普勒效应都进入微扰级;热波包 $35\mathrm{nm}$ 只占腰斑的 5%、阻塞半径($\approx 6.6\mu \mathrm{m}$,推导二)的 0.5%。原子位置的不确定性远小于一切相互作用与寻址尺度——这是原子阵列能做确定性路由的运动学根据。

散射、加热与搬运的绝热下限

囚禁光散射按双能级近似 ${\Gamma }_{\mathrm{sc}}\approx (\Gamma /\Delta )(U_0/\hslash )$:$\Gamma /2\pi =6.07\mathrm{MHz}$(D2 线宽)、$\Delta /2\pi \approx -32\mathrm{THz}$,而 1 mK 阱深对应光移 $U_0/h\approx 20.8\mathrm{MHz}$,故 ${\Gamma }_{\mathrm{sc}}\sim 4{\mathrm{s}}^{-1}$ 每毫开尔文阱深(850 nm;1064 nm 失谐大 3 倍,降到 $\sim 1{\mathrm{s}}^{-1}$)。每次散射注入约 $2E_R\approx 0.31\mu \mathrm{K}$ 反冲($E_R={\hslash }^2{k}_t^2/2m$),即 $\sim 1\mu \mathrm{K}/s$ 量级的加热——分钟级连续运转的 QRAM 工厂必须更远离共振或用魔波长,这是估算层面就能看清的工程约束。位置涨落进入门预算的定量推导见推导三。

搬运(恒加速度轨迹 $T(d)=T_0\sqrt{d/d_0}$ 的推导见 QRAM 工厂)另有独立下限:镊子轨迹的频率成分须远低于陷阱频率,否则惯性力 $-m{\ddot{x}}_t(t)$ 向运动态泵振动量子。最小间距搬运时间 $T\approx 33\mu \mathrm{s}$(QR004 取值)与 $1/{\omega }_{\perp }=1.1\mu \mathrm{s}$ 之比约 30——33 μs 不是任意常数,而是”固定绝热度”约束下 $\sqrt{d}$ 律在最小间距处的实现;再快就要付非绝热加热,直接吃推导三的温度预算。

推导二:里德堡阻塞门——机制与保真度账

阻塞机制:单激发流形与集体增强

两原子各以 Rabi 频率 $\Omega$ 耦合 $|g⟩\leftrightarrow |r⟩$(里德堡态),相互作用 $V|rr⟩⟨rr|$(范德瓦尔斯势 $V=C_6/r^6$):

$H=\frac{\hslash \Omega }{2}{\sum }_{i=1}^2(|g_i⟩⟨r_i|+h.c.)+V|rr⟩⟨rr|.$

当 $V\gg \hslash \Omega$ 时 $|rr⟩$ 失谐出共振、可绝热剔除,动力学限制在 $\{|gg⟩,|W⟩\equiv (|rg⟩+|gr⟩)/\sqrt{2}\}$ 单激发流形上。耦合矩阵元逐项算:

$⟨W|H|gg⟩=\frac{1}{\sqrt{2}}(⟨rg|+⟨gr|)\frac{\hslash \Omega }{2}(|rg⟩+|gr⟩)=\frac{\hslash \Omega }{\sqrt{2}}=\hslash \cdot \frac{\sqrt{2}\Omega }{2},$

即 $|gg⟩\leftrightarrow |W⟩$ 以 $\sqrt{2}\Omega$ 振荡——两条激发路径相干叠加的集体增强(collective enhancement)。阻塞把”两个原子”变成共享一个激发的复合体,这是原子阵列能把门做快、把路由做并行的微观根源。

CZ 序列与门时间

控制相位门用三脉冲序列 $\pi$(控制)-$2\pi$(目标)-$\pi$(控制):仅当控制比特处于 $|1⟩$ 时,其里德堡激发把目标失谐出共振、阻塞目标的 $2\pi$ 脉冲,于是 $|11⟩$ 不吃动力学相位而其余基矢吃满 $2\pi$——配适当单比特相位修正即 CZ。总脉冲面积 $4\pi$,即里德堡暴露时间

$t_R\approx \frac{4\pi }{\Omega }=200\mathrm{ns}(\Omega /2\pi =10\mathrm{MHz}),$

这是全部误差预算的时间单位。

误差预算:四项分解

阻塞泄漏。 有限 $V$ 使 $|rr⟩$ 以幅度 $\sim \Omega /2V$ 混入,泄漏概率 ${\epsilon }_{\mathrm{bl}}\sim (\Omega /V{)}^2$。取 $C_6/h\approx 862\mathrm{GHz}\cdot \mu {\mathrm{m}}^3$(${}^{87}\mathrm{Rb}$ $n\approx 70$ S 态文献典型值,本词条估算):

$V(d)/h=\frac{C_6/h}{d^3}:d=5\mu \mathrm{m}\Rightarrow 6.9\mathrm{GHz},d=3\mu \mathrm{m}\Rightarrow 31.9\mathrm{GHz},$

代入得 ${\epsilon }_{\mathrm{bl}}\sim 2\times 10^{-6}$ 与 $1\times 10^{-7}$——3–5 μm 间距下阻塞泄漏根本不是主导错误。阻塞半径

$r_b=(\frac{C_6}{\hslash \Omega }{)}^{1/6}\approx 6.6\mu \mathrm{m}(\Omega /2\pi =10\mathrm{MHz})$

的作用不是把门做准,而是划定并行边界(推导四)。

里德堡衰变。 ${\epsilon }_{\mathrm{dec}}\approx {\Gamma }_n{t}_R$:$n\approx 70$ 的 Rb S 态室温黑体辐射下寿命 ${\tau }_n\approx 150\mu \mathrm{s}$(教科书值),代入 $t_R=200\mathrm{ns}$ 得 $1.3\times 10^{-3}$;计入黑体跃迁逃出门子空间与两原子因子,取 $1$–$3\times 10^{-3}$。

多普勒失谐。 推导三给 ${\epsilon }_{\mathrm{D}}\approx {\sigma }_{\phi }^2/2$,10 μK、对向传播两光子方案、200 ns 门:$\approx 4.9\times 10^{-4}$。

其余(激光相位/强度噪声、寻址、运动残余):合并记 $\sim 10^{-3}$。合计

$\epsilon \approx \underset{\text{阻塞}}{\underbrace{2\times 10^{-6}}}+\underset{\text{衰变}}{\underbrace{(1\text{–}3)\times 10^{-3}}}+\underset{\text{多普勒}}{\underbrace{4.9\times 10^{-4}}}+\underset{\text{其余}}{\underbrace{\sim 10^{-3}}}\approx 5\times 10^{-3}.$

这与 Cesa 等 2026 补充材料”当前实验场景 $\epsilon \simeq 5\times 10^{-3}$(约 500 原子)“的取值同量级:分解本身是本词条估算,但合计表明 $5\times 10^{-3}$ 不是神秘的经验地板,而是”$150\mu \mathrm{s}$ 寿命 $\times$ $200\mathrm{ns}$ 门 + $10\mu \mathrm{K}$ 多普勒”的直接算术。

通向 $10^{-4}$ 的逐项路线

通道	现状贡献	路线(教科书级标度)	预期
衰变	$1$–$3\times 10^{-3}$	低温抑黑体(77 K 约 3 倍;4 K 一个量级)或圆态($\sim 30\mathrm{ms}$)	$\lesssim 3\times 10^{-4}$
多普勒	$4.9\times 10^{-4}$	Raman 边带冷却至 $\sim 1\mu \mathrm{K}$	$\sim 5\times 10^{-5}$
其余	$\sim 10^{-3}$	激光锁定、脉冲整形、光移补偿	$\lesssim 10^{-4}$
阻塞	$2\times 10^{-6}$	$C_6\propto n^{11}$:$n:70\to 90$ 再降约 16 倍	忽略

四项同时到位即达 $\epsilon \sim 10^{-4}$——Cesa 等 2026 容错账目要求的工作点。表中倍数为标度推算,具体器件数字属本词条估算。

推导三:温度与位置涨落的退相干预算

高斯相位扩散:平均保真度的显式推导

设门期间原子累积随机相位 $\phi$(多普勒失谐与光移噪声的净效应),通道为 $\rho \to {\mathbb{E}}_{\phi }[e^{i\phi Z}\rho e^{-i\phi Z}]$。第一步,对任意输入 $|\psi ⟩=a|0⟩+b|1⟩$,旋转态与理想态的内积

$⟨{\psi }_{\phi }|\psi ⟩=|a{|}^2{e}^{i\phi }+|b{|}^2{e}^{-i\phi }=\cos \phi +i(|a{|}^2-|b{|}^2)\sin \phi ,$

对 $\phi$ 平均后 $\mathbb{E}[\sin \phi ]=0$(零均值),实部只剩 $\mathbb{E}[\cos \phi ]$:平均门保真度 $\bar{F}={\mathbb{E}}_{\phi }[\cos \phi ]$,与输入态无关。第二步,对高斯分布 $\phi \sim \mathcal{N}(0,{\sigma }_{\phi }^2)$ 用特征函数:

$\bar{F}=\frac{1}{\sqrt{2\pi }{\sigma }_{\phi }}\int e^{-{\phi }^2/2{\sigma }_{\phi }^2}\cos \phi d\phi =e^{-{\sigma }_{\phi }^2/2}.$

平均密度矩阵的相干项乘同一因子:

$$\bar{\rho }=\frac{1}{2}\left(\begin{array}{cc}1& e^{-{\sigma }_{\phi }^2/2}\\ e^{-{\sigma }_{\phi }^2/2}& 1\end{array}\right),1-\bar{F}\approx \frac{{\sigma }_{\phi }^2}{2}({\sigma }_{\phi }\ll 1).$$

多普勒项:$\epsilon \propto T{t}_R^2$ 与适用判据

失谐是速度的线性投影 $\delta ={\mathbf{k}}_{\mathrm{eff}}\cdot \mathbf{v}$。780 nm + 480 nm 对向传播的两光子里德堡激发,有效波矢

$k_{\mathrm{eff}}=k_{480}-k_{780}=2\pi (\frac{1}{480\mathrm{nm}}-\frac{1}{780\mathrm{nm}})=5.0\times 10^6{\mathrm{m}}^{-1},$

热速度弥散 ${\sigma }_v=3.1cm/s$(10 μK),故失谐弥散 ${\sigma }_{\delta }=k_{\mathrm{eff}}{\sigma }_v=2\pi \times 24.8\mathrm{kHz}$。适用判据先行:门时间远短于横向振动周期($t_R=200\mathrm{ns}\ll 2\pi /{\omega }_{\perp }=7.1\mu \mathrm{s}$)时,速度在门内冻结、在门间按麦克斯韦分布重新抽取,高斯扩散模型自洽;反之 $t_R\gtrsim$ 振动周期时速度门内自平均,退相干被运动平均化(motional averaging)压制——QRAM 的快门恰落在扩散区。于是

${\sigma }_{\phi }=k_{\mathrm{eff}}{\sigma }_v{t}_R=0.031⟹{\epsilon }_{\mathrm{D}}\approx \frac{{\sigma }_{\phi }^2}{2}=4.9\times 10^{-4}(10\mu \mathrm{K},200\mathrm{ns}).$

由 ${\sigma }_v\propto \sqrt{T}$ 得核心标度 ${\epsilon }_{\mathrm{D}}\propto T{t}_R^2$:

$T\backslash t_R$	0.1 μs	0.2 μs	0.5 μs
1 μK	$1.2\times 10^{-5}$	$4.8\times 10^{-5}$	$3.0\times 10^{-4}$
10 μK	$1.2\times 10^{-4}$	$4.9\times 10^{-4}$	$3.0\times 10^{-3}$
30 μK	$3.6\times 10^{-4}$	$1.5\times 10^{-3}$	$9.1\times 10^{-3}$

读表要点:未做亚多普勒冷却的 30 μK 原子在 0.5 μs 慢门下,多普勒一项就吃掉整个预算;1 μK + 200 ns 已压到 $5\times 10^{-5}$,不再是大头。

有限囚禁深度:关阱循环与光移账

里德堡态在近红外囚禁光下的极化率与基态反号且大得多,光移典型达 $2\pi \times$ 几十 MHz(依赖 $n$ 与波长;本词条估算)——这是多数实验在里德堡脉冲期间关阱的原因。关阱的代价可显式算:阱突然关闭,谐振子基态自由演化 $T_{\mathrm{off}}$ 后与初始高斯波包的重叠为

$|⟨0|e^{-i{H}_{\mathrm{free}}{T}_{\mathrm{off}}/\hslash }|0⟩{|}^2=\frac{1}{\sqrt{1+({\omega }_{\perp }{T}_{\mathrm{off}}/2{)}^2}},$

代入 $T_{\mathrm{off}}=200\mathrm{ns}$:$({\omega }_{\perp }{T}_{\mathrm{off}}/2{)}^2=7.8\times 10^{-3}$,每轮关-开循环损失 $\approx 0.4\%$ 运动重叠——与门错误同量级,且 $\propto T_{\mathrm{off}}^2$。替代路线是阱全程开、把里德堡光移当可校准失谐:未校准的旋转误差 $\sim ({\delta }_L/\Omega {)}^2$(${\delta }_L=0.03\times 20.8\mathrm{MHz}\approx 0.6\mathrm{MHz}$ 时为 $3.8\times 10^{-3}$);校准后的残余来自阱深噪声(强度噪声 $1\%$ → 失谐噪声 $\sim 6\mathrm{kHz}$ → 每门 ${\sigma }_{\phi }\approx 7.8\times 10^{-3}$ → $\epsilon \approx 3\times 10^{-5}$)与波包内光移梯度($2(x_{\mathrm{rms}}/w_0{)}^2{\delta }_L\sim$ kHz,更低)。阱开与阱关是把”光移账”换”运动学账”的权衡,两边账目都给全了。

空闲退相干:按测量轮记账

资源态蓝图的原子在两轮 SQPM 之间空闲,空闲退相干按测量轮 $\tau =500\mu \mathrm{s}$(Cesa 等 2026)逐轮记账:${\epsilon }_{\mathrm{idle}}\approx \tau /T_2$。超精细时钟态的 $T_2$ 经自旋回波可达秒级:$T_2=1\mathrm{s}$ 给 $5\times 10^{-4}$/轮(已超 $\epsilon =10^{-4}$ 预算),$T_2=10\mathrm{s}$ 给 $5\times 10^{-5}$/轮(达标)。未被查询路径经过的原子不进此账——桶队型误差局域化让空闲错误大部分与输出解纠缠(一般 CPTP 界见 一般噪声下的韧性),要紧的只是路径上的 ${\log }_{2}N$ 个节点。

位置涨落:相互作用噪声与串扰

镊子指向噪声 $\delta d$ 改变原子间距,相互作用按 $V\propto d^{-3}$ 缩放:$\delta V/V=-3\delta d/d$。对阻塞门不敏感($V$ 只当开关);对把 $V$ 或 $C_3/r^3$ 手性耦合当相位资源用的节点则直接进预算:手性路由器边长 $\sim 9.4\mu \mathrm{m}$(Palaiodimopoulos 等 2024),$\delta d=50\mathrm{nm}$ 即 $\delta J/J\approx 1.6\%$。光学串扰反而不成问题:聚焦束 $w\approx 1\mu \mathrm{m}$ 对 3 μm 外邻居的激发比 $e^{-2d^2/w^2}\sim 10^{-8}$;真正的并行约束是阻塞半径:相邻门若都落在 $r_b\approx 6.6\mu \mathrm{m}$ 内会互相阻塞,须全局面调度。

推导四:原子阵列到 BB 树的物理映射

节点=原子:“传桶”的原子实现

桶队架构(桶队架构)的 $2N-1$ 个路由器逐个映射为原子:$|0⟩,|1⟩$ 存地址选择(超精细态)、$|r⟩$ 供里德堡切换。地址信息以原子内态形式逐级”传桶”:父节点的超精细态经阻塞受控跃迁复制到左/右子节点之一,每条地址分支每层激活且只激活一个节点(激活门 $={\log }_{2}N$ 的 GLM 账见 评价指标);数据留在叶子原子的时钟态里,叶子被激活时施加 $Z^{D_{\mathbf{x}}}$ 相位读出——与资源态框架的 $W_D={⨂}_{\mathbf{x}}{Z}_{\mathbf{x}}^{D_{\mathbf{x}}}$ 同一语义。两能级/三能级路由器的电路差异见 两能级路由器。

平面 packing 账

链式布局(叶子排一行)线性尺度 $=$ 叶子数 $\times$ 间距:1 kB($N=2^{13}$、间距 3 μm)需 $8192\times 3\mu \mathrm{m}\approx 24.6\mathrm{mm}$——单区不可行。方形 packing:$2N-1=16383\approx 128^2$ 个节点恰可塞进 $\sqrt{2N-1}\times$ 间距 $\approx 384\mu \mathrm{m}$ 见方的区域(完全二叉树在 $O(\sqrt{N})$ 边长网格上有标准 pyramid 嵌入;存在性为本词条估算)。物理尺寸在已演示阵列量级内,但 $O(N)$ 个原子的制备、校准与并行门调度是从演示到器件的真实距离。

zoned 分区与搬运接口

Bluvstein 等 2024 的逻辑原子处理器已演示存储区/纠缠区/读出区的分区(zoned architecture)运转:原子相干搬运进纠缠区逐层做阻塞门、读出区做非破坏成像。这与 QRAM 工厂 的重排项($3\sqrt{6}T{N}^{1/4}$、每周期 6 次长程搬运)对接:树不必常驻纠缠区,而是”工厂区制备 → 重排 → 查询区消费”的时空流水;搬运期间运动学一致性由推导一的绝热条件保障。

阻塞半径 vs 间距:并行的几何约束

$r_b\approx 6.6\mu \mathrm{m}$ 大于 3–5 μm 间距:相邻查询路径物理上不可同时激活。出路有三:全局门按阻塞结构成批驱动(并行性编码进脉冲序列)、频率复用、时分交错。注意单一旋钮 $\Omega$ 连锁三者——降低 $\Omega$ 可缩 $r_b$,但线性拉长 $t_R$、推高衰变(${\epsilon }_{\mathrm{dec}}\propto 1/\Omega$)与多普勒(${\epsilon }_{\mathrm{D}}\propto 1/{\Omega }^2$):门速度、阻塞半径与并行度不可分账。

另两类原子节点:手性路由器与单里德堡记忆

手性路由器(Palaiodimopoulos 等 2024;机制详见 手性量子路由器)用三自旋磁通三角让控制比特的量子态决定流向,是”由量子态决定流向”的里德堡实现。其耦合 $|J|\approx 2\pi \times 200\mathrm{kHz}$ 给出节点开关时间 $\sim 1/(2|J|)\approx 2.5\mu \mathrm{s}$,文献报的退相干率 ${\Gamma }_{\mathrm{tot}}\approx 1/62,1/45,1/24\mu {\mathrm{s}}^{-1}$(0/77/300 K)给出相干窗口 24–62 μs。逐层串联的树深上限由比值决定:$2.5n\mu \mathrm{s}\lesssim 24\mu \mathrm{s}$(300 K)→ $n\lesssim 9$(512 单元);77 K 窗口 45 μs → $n\lesssim 18$。节点开关时间 × 深度 vs 相干窗口——这是把 QR039 模拟读数换算成 QRAM 尺度上限的可行性账(倍率属本词条估算,读数取自文献);代价还有原子数:$3\times (2N-1)$(三角占三原子),原子成本换被动手性。

单里德堡原子记忆(Patton & Fischer 2013)用 BEC 双组分宏观波函数加单原子里德堡两光子过程,实现 $\sim 10\mathrm{ns}$、宣称 97% 保真度的单比特写入/读出——标定里德堡介导接口的速度上限,但不提供树的规模。系综路线(Jiang 等 2019)在宏观原子系综中实现 105 qubit/210 存储单元、任意次序可编程读出且超过经典界——有存储、无逐单元逻辑,与阵列路线货币不同(系综的货币是光学深度/协同度,见 里德堡 EIT QRAM)。三类读数放在一起才能看清:“中性原子 QRAM”不是一条路线,而是一个谱系。

数值算例:两个工作点的可行性账

当前参数($\epsilon \approx 5\times 10^{-3}$)

纯级联桶队的 union bound 给 $1-F\lesssim n\epsilon$;$N=32$($n=5$)为 2.5%。但 Cesa 等 2026 模拟(含制备/测量/存储全链路错误)给 $F\approx 52\%$:两者之差是”全链路”与”只数激活门”的差,提醒可行性由最弱环节决定,不是路由门。同参数 $N=4,8,16,64$ 对应 $F\approx 91\%,80\%,66\%,35\%$;表面码 $d=3,5,7$ 下单次纠错查询对应 $N\approx 16,8,4$。硬件现状:约 500 原子级处理器;Bluvstein 等 2024 已达 280 物理比特的逻辑处理器。

1 kB 目标($\epsilon \to 10^{-4}$)

原子数 $2N-1=16383$;物理级联不保真度 $\epsilon n\approx 13\times 10^{-4}$、级联墙钟 $\sim n{t}_R\approx 2.6\mu \mathrm{s}$。对照资源态蓝图:在线 $T_{\mathrm{query}}\approx 13\mathrm{ms}$、不保真度 ${\log }^{\alpha }N$——$n=13$ 时 $\alpha =1.9$ 与 $\alpha =3.8$ 之差是 $13^{1.9}\approx 1.3\times 10^2$ vs $13^{3.8}\approx 1.7\times 10^4$,两个量级(偏置噪声的价值,口径同 评价指标)。张力很清楚:物理级联又快又准,却既无擦除预算也无容错接口;资源态路线慢三个量级,换来 ${\log }^{1.9}$ 标度与全 Clifford 在线。擦除转换(Scholl 等 2023)是粘合剂:把双比特门错误转成可检测擦除/相位错误,同时改善物理级联的可观测性与资源态模拟的 $\alpha$。

预算汇总

预算项	表达式	10 μK / 200 ns 现状	1 μK / 快门目标
多普勒	$\frac{1}{2}(k_{\mathrm{eff}}{\sigma }_v{t}_R{)}^2$	$4.9\times 10^{-4}$	$\sim 5\times 10^{-5}$
里德堡衰变	${\Gamma }_n{t}_R$	$1$–$3\times 10^{-3}$	$\lesssim 3\times 10^{-4}$
阻塞泄漏	$(\Omega /V{)}^2$	$2\times 10^{-6}$	忽略
关阱/光移	$({\omega }_{\perp }{T}_{\mathrm{off}}/2{)}^2$ 或 $({\delta }_L/\Omega {)}^2$	$10^{-4}$–$10^{-3}$	$\lesssim 10^{-4}$
空闲(每测量轮)	$\tau /T_2$	$5\times 10^{-4}$($T_2=1\mathrm{s}$)	$5\times 10^{-5}$($T_2=10\mathrm{s}$)
合计	—	$\approx 5\times 10^{-3}$	$\sim 10^{-4}$

讨论:边界与开放问题

其一,速度与保护的互换是结构性的:物理级联的微秒级查询与资源态的 13 ms 查询不是同一商品的两种价格——前者卖裸查询保真度,后者卖可纠错性与擦除预算。QRAM 作为算法子系统使用时必须连同纠错层一起记账(成本口径)。

其二,里德堡衰变的双重身份:它既是门预算里最大的一项,又经擦除转换成为 Z 偏置噪声的来源——把 $\alpha$ 从 3.8 压到 1.9 的正是它(Scholl 等 2023、Cesa 等 2026)。同一物理过程在两个账本上符号相反,这是硬件词条区别于纯电路词条之处。

其三,产率与布图仍是空白:万级原子的并行门调度、重排路径规划、擦除探测(非破坏成像,单轮 $\tau \approx 500\mu \mathrm{s}$)对工厂周期的挤占,均无实验读数;QRAM 工厂 的吞吐推导在这些环节取的是文献标称值。平台对比的总体定位见 硬件平台全景。

其四,平台迁移的边界:资源态蓝图的在线协议对”能搬运”的平台普遍可用(Cesa 等 2026 明确列出囚禁离子、超导、光子),但推导一、三的运动学账(陷阱频率、多普勒、重排绝热性、散射加热)是中性原子专属的——迁移须整体重算,而非只换门时间参数。

延伸阅读

• 资源态 QRAM
• QRAM 工厂
• Clifford 容错接口
• 硬件平台全景
• 手性量子路由器
• 里德堡 EIT QRAM
• 桶队架构
• 两能级路由器
• 一般噪声下的韧性
• 评价指标