混合声学 QRAM

现代超导 QRAM 硬件路线的奠基之作:量子信息住进高 Q 声子模式,transmon 以失谐"虚拟激发"工程声子间门,其退相干被逆 Purcell 效应压低。本词条从 cQAD 哈密顿量出发,把虚拟耦合速率、模式间距工程、门保真度优化与查询级误差预算逐级推出。

机制与器件回顾

Hann 等 2019(PRL 123, 250501)提出多模电路量子声动力学(circuit quantum acoustodynamics, cQAD)架构:量子信息存于高相干声子模式,通过向 transmon 施加失谐驱动、以其 Kerr 非线性作四波混频器(four-wave mixer),以”虚拟激发”方式工程声子–声子耦合。声速与光速相差五个数量级,微波频段声子波长在亚微米量级,单个芯片上”数百声子模式同时耦合一个 transmon”已被实验演示;GHz 频段实测质量因子为 SAW $\sim 10^5$、BAW $\sim 10^8$、声子晶体 $\sim 10^{10}$。双音驱动得到束分裂器型耦合(演化 $\pi /2g_v$ 为 SWAP、$\pi /4g_v$ 为 50:50 束分裂器);单音驱动借助辅助模得到三模耦合,演化 $\pi /g_v$ 实现 CZ。束分裂器、SWAP 与 CZ 构成双轨编码下的通用门集。关键优势在于:虚拟门中 transmon 退相干经逆 Purcell 效应被 $(g/\delta {)}^2$ 压低,在长寿命声子区间门保真度可超 99% 并优于直接激发 transmon 的方案。代价是慢——现实虚拟耦合速率仅 10–100 kHz。工程上,近似均匀的模式间距会使共振条件简并,必须主动工程化非均匀模式间距(如与微波谐振器杂化产生约 1 MHz 频移);模块化分析给出每模块最优模式数约 12–15,单芯片数百声子模式可同时耦合一个 transmon。

以下推导依次回答:虚拟耦合从哪来、模式为何必须错频、transmon 退相干如何被压低而门保真度最优值多大(推导一至三)、行程声子对照的传播账(推导四)、整次查询的操作数与误差预算(推导五)。除标注”文献主张”外均为教科书级推导或本词条估算;速率符号均指角频率(除以 $2\pi$ 得 Hz)。

推导一:从 cQAD 哈密顿量到虚拟声子门

多模哈密顿量与色散消元

transmon 建模为带 Kerr 非线性的非谐振子,与 $M$ 个声子模式压电耦合 Hann 等 2019:

$$H={\omega }_q{q}^{†}q-\frac{\alpha }{2}{q}^{†}{q}^{†}qq+\sum _k\left({\omega }_k{m}_k^{†}{m}_k+g_k{q}^{†}{m}_k+g_k^{\ast }q{m}_k^{†}\right)+H_d,H_d=\sum _j\left({\Omega }_j{q}^{†}{e}^{-i{\omega }_{j}t}+H.c.\right),$$

其中 $q$、$m_k$ 为 transmon 与第 $k$ 个声子模式的湮灭算符,$\alpha$ 为非谐性(典型 $\alpha /2\pi \approx 150\mathrm{MHz}$),$g_k$ 为压电耦合(典型几 MHz),${\Omega }_j,{\omega }_j$ 为第 $j$ 个驱动音的幅度与频率。进入以 $H_0={\omega }_q{q}^{†}q+{\sum }_k{\omega }_k{m}_k^{†}{m}_k$ 定义的旋转系,记失谐 ${\delta }_k={\omega }_k-{\omega }_q$、${\delta }_j={\omega }_j-{\omega }_q$:

$$H=\sum _j\left({\Omega }_j{q}^{†}{e}^{-i{\delta }_{j}t}+H.c.\right)+\sum _k\left(g_k{m}_k{q}^{†}{e}^{-i{\delta }_{k}t}+H.c.\right)-\frac{\alpha }{2}{q}^{†2}{q}^2.$$

色散区 $|{\delta }_k|\gg g_k$ 下做逐阶消元(Nigg 型变换):引入无量纲量 ${\lambda }_k=g_k/{\delta }_k\ll 1$、${\xi }_j={\Omega }_j/{\delta }_j\ll 1$,用幺正变换把耦合项与驱动项吸收进 transmon 算符,等效于代换

$$q⟶Q=q+\sum _j{\xi }_j{e}^{-i{\delta }_{j}t}+\sum _k{\lambda }_k{m}_k{e}^{-i{\delta }_{k}t},$$

物理图像:transmon 算符 = 裸比特 + 各驱动的相干振荡 + 各声子模式经 ${\lambda }_k$ 的虚跃迁尾巴。代换后唯一剩下的相互作用来自 Kerr 项:

$$H=-\frac{\alpha }{2}{Q}^{†}{Q}^{†}QQ.$$

双音驱动:束分裂器与 SWAP

展开 $Q^{†2}{Q}^2$ 并挑出双共振项。取两个驱动音 ${\omega }_{1,2}$ 与两个模式 $A,B$,含四个不同算符的项为

$$-\frac{\alpha }{2}{Q}^{†2}{Q}^2\supset -2\alpha {\xi }_1^{\ast }{\xi }_2{\lambda }_A{\lambda }_B^{\ast }{m}_A{m}_B^{†}{e}^{i({\delta }_1+{\delta }_B-{\delta }_2-{\delta }_A)t}+H.c.,$$

系数来源:${\xi }_1^{\ast },{\lambda }_B^{\ast }$ 取自两个 $Q^{†}$、${\xi }_2,{\lambda }_A$ 取自两个 $Q$,排列数 $2\times 2=4$,乘 $-\alpha /2$ 得 $-2\alpha$。当驱动差频匹配模式差频,

$${\delta }_2-{\delta }_1={\delta }_B-{\delta }_A⟺{\omega }_2-{\omega }_1={\omega }_B-{\omega }_A,$$

相位因子化为常数,RWA 后留下束分裂器型有效哈密顿量 $H=g_v^{(1)}{m}_A{m}_B^{†}+H.c.$,速率

$$g_v^{(1)}=-2\alpha \frac{g_A}{{\delta }_A}\frac{g_B^{\ast }}{{\delta }_B}\frac{{\Omega }_1^{\ast }}{{\delta }_1}\frac{{\Omega }_2}{{\delta }_2}(1-{\beta }^{(1)}).$$

这是四波混频的量子版本:一个声子在 $A$ 湮灭、于 $B$ 再生,能量差由双音偿付,transmon 只提供非线性、全程留在基态——“虚拟激发”由此得名。数值代入基线参数($g/2\pi =10\mathrm{MHz}$、$\delta /2\pi =100\mathrm{MHz}$ 即 $\lambda =0.1$、$\xi =0.17$、$\alpha /2\pi =150\mathrm{MHz}$):$g_v^{(1)}=2\alpha {\lambda }^2{\xi }^2=5.4\times 10^{5}rad/s$,即 $g_v^{(1)}/2\pi \approx 87\mathrm{kHz}$,与文献数值 Floquet 计算给出的上限 $|g_v^{(1)}|/2\pi <100\mathrm{kHz}$ 一致;对应 $t_{\mathrm{SWAP}}=\pi /2g_v\approx 5.8\mu \mathrm{s}$、$t_{\mathrm{BS}}=\pi /4g_v\approx 2.9\mu \mathrm{s}$。

单音驱动:三模耦合与 CZ

单个驱动音 ${\omega }_1$ 与三个模式 $A,B,C$ 的共振项同理挑出($Q^{†}$ 贡献 ${\xi }_1^{\ast }$,$Q$ 贡献 ${\lambda }_A{m}_A$、${\lambda }_C{m}_C$,$Q^{†}$ 贡献 ${\lambda }_B^{\ast }{m}_B^{†}$,排列数仍为 4):

$$H=g_v^{(2)}{m}_A{m}_B{m}_C^{†}+H.c.,g_v^{(2)}=-2\alpha \frac{g_A}{{\delta }_A}\frac{g_B}{{\delta }_B}\frac{g_C^{\ast }}{{\delta }_C}\frac{{\Omega }_1^{\ast }}{{\delta }_1}(1-{\beta }^{(2)}),$$

共振条件 ${\omega }_1={\omega }_A+{\omega }_B-{\omega }_C$:两个声子合成一个、或一个声子劈成两个,能量由单音偿付。以 $C$ 为辅助模并初始化在 $|0{⟩}_C$,该耦合对 $|0,1⟩$ Fock 编码实现 CZ:在二维子空间 $\{|110⟩,|011⟩{\}}_{ABC}$ 中,$H_{\mathrm{eff}}=g_v^{(2)}|011⟩⟨110|+H.c.$ 是标准 Rabi 对,

$$U(t)|110⟩=\cos (g_v^{(2)}t)|110⟩-i\sin (g_v^{(2)}t)|011⟩,$$

演化 $\pi /g_v^{(2)}$ 给出 $|110⟩\to |001⟩\to -|110⟩$;其余基矢(如 $|100⟩$ 因 $m_B$ 打空、$|111⟩$ 因 $m_C^{†}$ 打满)不被触及,于是 $|11{⟩}_{AB}$ 获得相对几何相位,即 CZ。这类三模耦合不守恒声子数,但门只需其周期演化,一个周期内总声子数回归。

修正项:β 因子、AC Stark 移与并行操作

展开 $Q^{†2}{Q}^2$ 时其余项并非全部快变:最低阶修正给出 ${\beta }^{(1)}={\beta }^{(2)}=\frac{\alpha }{{\delta }_B+{\delta }_1+\alpha }$。SWAP 的驱动远失谐,该修正可忽略;CZ 中 ${\delta }_1,{\delta }_B$ 可与 $\alpha$ 同量级——当 ${\delta }_B+{\delta }_1\to -\alpha$ 时分母趋零,微扰失效(两光子共振打到 transmon 的 $|f⟩$ 能级,须避开),一侧压低、一侧反常增强,$|g_v^{(2)}|/2\pi <25\mathrm{kHz}$ 的天花板一半来自这里。强驱动还通过 AC Stark 移把 transmon 推离声子:$S=-2\alpha {\sum }_j|{\Omega }_j{|}^2/[{\delta }_j({\delta }_j+\alpha )]$,所有失谐须换为 Stark 移后的 $\tilde{\delta }$,共振条件改写为 ${\omega }_2-{\omega }_1=({\omega }_B+S_B)-({\omega }_A+S_A)$;Stark 移使 $g_v(\xi )$ 非单调,有效驱动被限制在 $\xi \le {\xi }_{\mathrm{crit}}$($\xi \lesssim 0.4$ 内与 Floquet 数值一致)。并行操作是虚拟门的独特优势:三路驱动可同时共振两组 SWAP(如 ${\omega }_2-{\omega }_1={\omega }_B-{\omega }_A$、${\omega }_3-{\omega }_2={\omega }_D-{\omega }_C$),但 $d$ 路驱动须核查全部 $d(d-1)/2$ 个差频组合以防误共振,这是并行度的实际上限之一。

推导二:声速、模式间距与频率拥挤

从声速到自由谱距

声学器件的色散关系是后续一切账目的几何基础。半无限基底上的瑞利表面声波(Rayleigh SAW)无色散:

$$\omega =v_R|k|,v_g=\frac{d\omega }{dk}=v_R\approx 0.9v_t,$$

群速度等于相速度,典型 $v_R\sim 3km/s$(教科书数值):GHz 频率下波长 $\lambda =v_R/f\sim 0.5\mu \mathrm{m}$,“单芯片数百模式”的密度来源正在于此。BAW/声子晶体谐振器则是驻波腔:长 $L$ 的法布里–珀罗声腔的模式间距即自由谱距(free spectral range, FSR)

$$\nu =\frac{\pi v}{L}\left(\text{两端反谐振边界}\right)⟹L=\frac{\pi v}{\nu };$$

本词条估算:取纵波 $v\approx 5km/s$、$\nu /2\pi =10\mathrm{MHz}$,得 $L\approx 250\mu \mathrm{m}$,与文献中数百微米量级的 BAW 衬底一致;$1\mathrm{GHz}$ 带宽内约容纳 100 个 FSR 间隔的模式,“数百声子模式”由此而来。声子晶体则靠带隙内缺陷模,模式频率由单个谐振器几何逐一设定,天然免于均匀间距问题。

非均匀度 Δν 与三种工程方案

均匀间距 ${\omega }_{j+1}-{\omega }_j=\nu$ 会使选择性失效:双音共振条件 ${\omega }_2-{\omega }_1={\omega }_B-{\omega }_A$ 对所有等差模式对同时成立,驱动无法只选中目标对。文献为此定义非均匀度 $\Delta \nu$(大致为模式间距变化的尺度):给定存量子信息的模式集 $\mathcal{S}$ 与拟耦合模式对集 $\mathcal{P}$,

$$\Delta \nu =\underset{\{i,j\}\in \mathcal{P}}{\min }\underset{\{k\in \mathcal{S},\ell \}\ne \{i,j\}}{\min }|{\nu }_{ij}-{\nu }_{k\ell }|,{\nu }_{ij}=|{\omega }_i-{\omega }_j|,$$

即最坏情况下目标共振条件与其他所有共振条件的最小间隔;高选择性要求 $g_v/\Delta \nu \ll 1$。工程化方案有三:(i) 外部模式杂化(点缺陷型)——声子与微波谐振器模式杂化,带宽 $\mathcal{D}$ 内模式被推开约 1 MHz(实验已演示),信息可经 SWAP 进出 $\mathcal{D}$;(ii) 双模式族(周期型)——transmon 同时耦合两组 FSR 不同的模式族,跨族耦合在带宽 $\approx \nu /\Delta \nu$ 的窗口 $\mathcal{S}$ 内全可选择性实现,$\Delta \nu =|{\nu }_1-{\nu }_2|$,三模耦合还需外加高相干外部模式;(iii) 复合谐振器——界面部分反射调制 FSR。注意 $\Delta \nu$ 依赖 $\mathcal{S},\mathcal{P}$ 的选取:选择性与有效规模/连通度存在权衡,本词条估算给出其操作含义——每模块可选地址模式数 $|\mathcal{S}|\approx \nu /\Delta \nu$,取 $\nu /2\pi =10\mathrm{MHz}$、$\Delta \nu /2\pi =0.85\mathrm{MHz}$ 得 $|\mathcal{S}|\approx 12$,这正是文献量子体积分析给出”每模块最优模式数 12–15”($M\le 12$ 于其 BAW 例、$M\approx 15$ 于声子晶体例)的另一面。

推导三:逆 Purcell 效应与门保真度优化

穿衣衰减率

色散消元后的”声子”是穿衣(dressed)模式。对单模作二阶微扰,声子型本征态为

$$|{\tilde{1}}_k⟩=|1_k,0_q⟩+{\lambda }_k|0_k,1_q⟩+O({\lambda }_k^2),$$

transmon 分量幅度 ${\lambda }_k=g_k/{\delta }_k$;transmon 以速率 $\gamma$ 衰减时,穿衣声子按其 qubit 分量权重获得附加衰减率 $\gamma {\lambda }_k^2$——这就是逆 Purcell 效应(inverse Purcell effect):cQED 中 qubit 因身居腔内而被加速,此处声子则因拖着 qubit 尾巴而染上其耗散。加上声子固有衰减 $\kappa$,第 $j$ 模穿衣衰减率(驱动把失谐改为 $\tilde{\delta }$)

$${\kappa }_{\gamma }^j=\kappa +\gamma (g_j/{\tilde{\delta }}_j{)}^2.$$

数值感受:transmon $T_1=100\mu \mathrm{s}$ 给 $\gamma =10^4{\mathrm{s}}^{-1}$,$g/\delta =0.1$ 时逆 Purcell 项 $=100{\mathrm{s}}^{-1}$($/2\pi \approx 16\mathrm{Hz}$),已可与 Q=10⁹ BAW 的固有 $\kappa /2\pi \approx 6\mathrm{Hz}$($f=6\mathrm{GHz}$)同量级甚至反超——存储器的天花板由 qubit 尾巴决定。消掉中间量还能得到一个设计恒等式(本词条推导,双模同 $\lambda$ 同 $\xi$ 特例):由 $g_v=2\alpha (g/\delta {)}^2{\xi }^2$ 反解 $(g/\delta {)}^2=g_v/(2\alpha {\xi }^2)$,故逆 Purcell 项 $=\gamma g_v/(2\alpha {\xi }^2)$,即在驱动幅度受限时,逆 Purcell 尾巴正比于目标门速,出路只有更大的非线性 $\alpha$、更满的驱动 $\xi$ 或更慢的门。

不保真度公式与 2/3 幂最优

虚拟门不保真度 = 退相干项 + 频率拥挤(spectral crowding)项 Hann 等 2019:

$$1-{\mathcal{F}}_v\approx {\bar{\kappa }}_{\gamma }\frac{c_v\pi }{2g_v}+(\frac{g_v}{\Delta \nu }{)}^2,1-{\mathcal{F}}_d\approx (\kappa +\gamma )\frac{c_d\pi }{2g_d}+(\frac{g_d}{\nu }{)}^2,$$

第一项 = 总退相干率 × 门时长($t_v=c_v\pi /2g_v$:$c_v=1$ 于 SWAP、$c_v=2$ 于 CZ;三模时 ${\bar{\kappa }}_{\gamma }=({\kappa }_{\gamma }^A+{\kappa }_{\gamma }^B+{\kappa }_{\gamma }^C)/2$,因子 1/2 来自门过程中两个声子合为一个的时间平均)。第二项:门驱动对其他模式对是失谐驱动,误激发概率 $\sim {\sum }_n(g_v/{\delta }_n^{\prime }{)}^2$,多余耦合依次错开 $\Delta \nu$ 的整数倍,级数 ${\sum }_{m\ge 1}1/m^2={\pi }^2/6$ 收敛,量级即 $(g_v/\Delta \nu {)}^2$——与系统模式总数无关,这是该架构可扩展性的关键。对 $g_v$ 优化 $f(g)=\frac{{\bar{\kappa }}_{\gamma }{c}_v\pi }{2g}+\frac{g^2}{\Delta {\nu }^2}$:令 $f^{\prime }(g)=0$ 得 $g^3={\bar{\kappa }}_{\gamma }{c}_v\pi \Delta {\nu }^2/4$;在最优点 $a/g=2b{g}^2$(两项之比 2:1),代入得幂律

$$g_{v,\mathrm{opt}}=(\frac{{\bar{\kappa }}_{\gamma }{c}_v\pi \Delta {\nu }^2}{4}{)}^{1/3},1-{\mathcal{F}}_{v,\mathrm{opt}}\approx \frac{3}{2}[\frac{c_v\pi {\bar{\kappa }}_{\gamma }}{\sqrt{2}\Delta \nu }{]}^{2/3},$$

不保真度只按 ${\bar{\kappa }}_{\gamma }^{2/3}$ 变差——2/3 幂把指数级严苛软化成温和的幂律。数值验算(SWAP,基线参数):取 ${\kappa }_{\gamma }^j\approx 138{\mathrm{s}}^{-1}$(上节),${\bar{\kappa }}_{\gamma }=2{\kappa }_{\gamma }\approx 276{\mathrm{s}}^{-1}$、$\Delta \nu =2\pi \times 1\mathrm{MHz}$:$g_{v,\mathrm{opt}}=(276\times \pi \times (6.28\times 10^6{)}^2/4{)}^{1/3}=2.0\times 10^5{\mathrm{s}}^{-1}$,即 $/2\pi \approx 33\mathrm{kHz}$,恰好落在可实现的 10–100 kHz 窗口内;此时退相干项 $=276\times 7.7\mu \mathrm{s}\approx 2.1\times 10^{-3}$、拥挤项 $=(0.033{)}^2\approx 1.1\times 10^{-3}$,合计 $1-\mathcal{F}\approx 0.3\%$。CZ($c_v=2$、${\bar{\kappa }}_{\gamma }\approx 207{\mathrm{s}}^{-1}$)最优点 37 kHz 超出 25 kHz 上限,按 25 kHz 工作得 $1-\mathcal{F}\approx 0.5\%$、$t_{\mathrm{CZ}}=20\mu \mathrm{s}$。敏感性:若 ${\bar{\kappa }}_{\gamma }$ 恶化 10 倍(Q 跌至 10⁸ BAW 水平),失真只涨 $10^{2/3}\approx 4.6$ 倍(SWAP $\sim 1.5\%$),仍是可用门。

虚拟 vs 直接:判据与数值

直接门把信息共振搬进 transmon($t\sim 1/g_d$ 量级内声子与 qubit 各占一半时间),退相干率为 $\kappa +\gamma$,transmon 与声子等权贡献;两式在各自最优点给出判据

$${\mathcal{F}}_v>{\mathcal{F}}_d⟺{\kappa }_{\gamma }\lesssim (\kappa +\gamma )\frac{\Delta \nu }{\nu }$$

(推导:在最优点直接比较两个 2/3 幂表达式,$c$ 因子吸收进 $\lesssim$)。基线参数下右端 $(37.7+10^4)\times 0.1\approx 10^3{\mathrm{s}}^{-1}\gg 138{\mathrm{s}}^{-1}$,虚拟门全胜;只有当声子相干不再比 transmon 长得多时判据才反转。文献主张:长寿命声子区间门保真度可超 99% 并优于直接激发 transmon 的方案,其基准参数域的数值计算支持这一结论。此外虚拟门兼容多声子编码与部分玻色纠错码(直接门只能 $|0,1⟩$ 编码),且可并行。

驱动上限:ξcrit 与驱动诱导加热

虚拟门的”慢”有三重来源:色散压低(${\lambda }^2$)、四波混频再压低(${\xi }^2$)、以及两重上限——解析式仅在 $\xi \lesssim 0.4$ 可信,且 AC Stark 移造成 $g_v(\xi )$ 非单调,超过 ${\xi }_{\mathrm{crit}}$ 后门速不升反降。更隐藏的一项是驱动诱导加热:驱动通过 transmon 自身耗散与高频退相干噪声把 qubit 泵热,双音情形加热率正比于 $\gamma |\alpha {\xi }^2/(2\delta +\alpha ){|}^2$ 与 $2{\gamma }_{\varphi }^{(\mathrm{hf})}|\xi {|}^2$ 之和(SM 式 S42),弱驱动下相对逆 Purcell 项可忽略;并行门增多时该加热线性累加,可能成为保真度的下一堵墙。这也是”超导电路耦合是否从根本上限制声子相干”这一公开问题在门层面的投影。

推导四:行程声子的群速度与传播损耗账

同一声学谱系的另一支不走驻波而走行程(itinerant)声子 Wang 等 2025(详见 声子路由 QRAM),其误差账结构不同,值得对照。无色散关系 $\omega =v_R|k|$ 下,波包以群速度 $v_R\approx 3km/s$ 直线传播;若行程声子寿命为 ${\tau }_{ph}$(该方案器件参数取 $2\mu \mathrm{s}$),定义衰减长度

$$\ell =v_R{\tau }_{ph}\approx 6\mathrm{mm},$$

则传播距离 $L$ 的损耗概率为 $1-e^{-L/\ell }\approx L/\ell$:片上 $L=1\mathrm{mm}$ 即损耗 $\sim 16\%$——行程方案必须靠马赫–曾德尔构型把传播压到最短、并靠混合双轨编码把损耗探测出来(后选择)。对照驻波账目:驻波模式损耗 $={\kappa }_{\gamma }\times$存储时长,与几何无关、只与时间成正比;QRAM 查询中驻波方案的有效占据时间为 $O(n)T_q$,行程方案则为 $O(n)$ 次穿越、每次 $L/\ell$。本词条估算:两者同为 $O(n)$ 的线性账,常数不同——驻波案的常数是 ${\kappa }_{\gamma }{T}_q$(慢门的代价),行程案的常数是 $nL/\ell$(快门但传播贵);哪边划算取决于 ${\kappa }_{\gamma }{T}_q$ 与 $L/\ell$ 之比,这正是后续 3D 腔 Weiss 等 2024(静态高 Q,门慢)与行程声波(动态低 Q,门快)两条路线分野的定量根源。

推导五:QRAM 查询的操作数账与误差预算

路由器的门合成:BS–CZ–BS

QRAM 的基本砖块是量子路由器:四个声子模式(顶部输入、路由模式 $\varphi$、左右输出),四阶段操作——初始化(顶模式与 $\varphi$ 换态,存下地址比特)、下行路由、上行路由、提取;下行路由 = 1 个 SWAP + 1 个 CSWAP(当 $\varphi =|0⟩$ 时 CSWAP 平凡、当 $\varphi =|1⟩$ 时 SWAP 因两目标模皆真空而平凡),上行与提取按时间反演执行(SM III)。CSWAP 由 CZ 加束分裂器合成,可在单声子子空间内显式验证(本词条推导):束分裂器 $U(\theta )=e^{-i\theta (m_a^{†}{m}_b+m_a{m}_b^{†})}$ 在 $\{|10⟩,|01⟩\}$ 子空间是旋转,

$$U(\theta )\doteq \left(\begin{array}{cc}\cos \theta & -i\sin \theta \\ -i\sin \theta & \cos \theta \end{array}\right),$$

控制模 $c=1$ 分支上 CZ(控制 $c$、目标 $b$)对 $b$ 模为 $D=\mathrm{diag}(1,-1)$;复合 $W=U({\theta }_1)DU({\theta }_2)$ 的矩阵元只依赖差 $\Delta ={\theta }_1-{\theta }_2$:

$$W\doteq \left(\begin{array}{cc}\cos \Delta & i\sin \Delta \\ i\sin \Delta & -\cos \Delta \end{array}\right)\stackrel{\Delta =\pi /2}{\to }i\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right),$$

即 50:50 束分裂器–CZ–束分裂器在单声子子空间内给出受限 SWAP(差一个相对相位,可并入驱动相位账目;路由逻辑只依赖声子占据哪条路径,不受影响)。当两目标模总声子数 $\ge 2$ 时该合成失效——这正是文献图 4(b)“仅在总声子数 <2 的相关子空间内实现该门”告诫的内容。

操作数账与查询时间

地址比特进树、总线取数、再全部退树,与 bucket-brigade 架构 的接力结构逐层相同(查询等距映射的一般语义见 QRAM 总览)。本词条按上述门表清点原语(估算;SWAP/BS 按 $2\pi \times 100\mathrm{kHz}$、CZ 按 $2\pi \times 25\mathrm{kHz}$ 上限):地址比特 $j$ 下行穿过 $j-1$ 个路由器、上行再穿过 $j-1$ 个,每次过路 = 1 全块(SWAP+CSWAP = 3 个二模门 + 1 个三模门);初始化与提取各 1 个独立 SWAP;总线往返 $2n$ 次过路。合计全块数 $P=n^2+n$、独立 SWAP $2n$,即

$$N_{2\mathrm{mode}}=3(n^2+n)+2n=3n^2+5n,N_{3\mathrm{mode}}=n^2+n.$$

$n$	$N=2^n$	全块 $P$	二模门	三模门	串行时长	并行时长($4n$ 层)
3	8	12	42	12	0.39 ms	0.36 ms
7	128	56	182	56	1.75 ms	0.84 ms
10	1024	110	350	110	3.4 ms	1.2 ms

虚拟门可跨路由器并行(同层路由器同时驱动),把深度从 $\Theta (n^2)$ 压回 $O(n)$ 层,但须按推导一核查全部差频组合;硬件上,全树需 $3N-2$ 个声子模式(每路由器顶/路由/左右端口,相邻共享;本词条估算),$n=4$ 时 46 个——已逼近”单谐振器数百模式”的容量边界,$n=10$ 时约 3070 个,必须模块化(每模块 12–15 模式、各带一个 transmon,分域独立控制)。经典数据走”总线 $|+⟩$ + 底层 CZ”的免拷贝读出(相位对 $|0⟩$ 平凡,不产生多余拷贝);只读经典库可进一步把数据完全存于经典阵列、以与辅助模接连两次 SWAP 的驱动相位差完成写总线——底层 $N$ 个数据模式整个从量子硬件中消失,这是”修改电路变体”最实惠的操作数减免。量子数据则由 $|1⟩$ 指针比特随路径下行、经受控 SWAP 抽出,满足不可克隆定理的限制;无误查询中指针全程保持 $|1⟩$ 且与系统解纠缠,查询后测量指针即可 herald 错误。

声子损耗账与查询保真度

声子损耗下 $1-F\sim \epsilon \log N$,论文并断言对任意独立非相干错误均为对数标度——这一断言两年后被 Hann 等 2021 在一般噪声模型下严格证明(见 噪声韧性的严格证明;union bound 等通用账目工具见 评价指标,逐层相消的微观机制见 bucket-brigade 架构)。现在把 $\epsilon$ 落到本平台参数上(本词条估算):损耗概率 $\epsilon \approx {\kappa }_{\gamma }{T}_q$,$T_q$ 为查询时长;有效占据–时间预算约 $2n{T}_q$($n$ 个地址声子全程驻留,总线/指针往返折合 $n$ 个全程),故每次查询的期望损耗事件数

$$E_{\mathrm{loss}}\approx 2n{\kappa }_{\gamma }{T}_q,1-F_{\mathrm{loss}}\lesssim c{E}_{\mathrm{loss}}(c\text{为}O(1)\text{常数}),$$

对数标度藏在 $T_q=O(\log N)$ 里:并行深度 $4n$ 层。代入三档存储器参数($T_q=1.2\mathrm{ms}$,$n=10$):

存储器	$\kappa /2\pi$	逆 Purcell 项	${\kappa }_{\gamma }$	$\tau =1/{\kappa }_{\gamma }$	$E_{\mathrm{loss}}$
BAW Q=10⁹($\delta /2\pi =100\mathrm{MHz}$)	6 Hz	$100{\mathrm{s}}^{-1}$	$\approx 138{\mathrm{s}}^{-1}$	7 ms	$\approx 3.3$(不可用)
声子晶体 Q=10¹⁰($\delta /2\pi =1\mathrm{GHz}$)	1 Hz	$1{\mathrm{s}}^{-1}$	$\approx 7.3{\mathrm{s}}^{-1}$	0.14 s	$\approx 0.18$
路线图:Q=10¹⁰ 且 $T_1=10\mathrm{ms}$	0.1 Hz	$10^{-2}{\mathrm{s}}^{-1}$	$\approx 0.6{\mathrm{s}}^{-1}$	1.6 s	$\approx 0.015$

第一行说明基线参数撑不住整次查询:7 ms 存储寿命对 1.2 ms 查询只是杯水车薪;这正是论文”可行性要求 $1/\kappa \gg 1/g_v$“的定量含义:门级只要求寿命 $\gg 10\mu \mathrm{s}$,查询级要求 ${\kappa }_{\gamma }\lesssim {\delta }_F/(2n{T}_q)$(失真预算 ${\delta }_F=1\%$、$n=10$ 时 ${\kappa }_{\gamma }\lesssim 0.4{\mathrm{s}}^{-1}$),两者相差三个数量级以上。门错误预算同样入账:单门 $0.3\%$–$0.5\%$ 时,Hann 定理的非空洞区要求 $\epsilon T\log N\le 1/4$,即 $\epsilon \lesssim 1/(4n^2)\approx 2.5\times 10^{-3}$($n=10$)——基线 $\Delta \nu /2\pi =1\mathrm{MHz}$ 下恰好压线;把 $\Delta \nu$ 提到 $3\mathrm{MHz}$ 可使门失真按 $\Delta {\nu }^{-2/3}$ 降到 $\sim 0.2\%$ 以下、进入非空洞区,但选择性窗口 $\nu /\Delta \nu$ 随之缩到 3–4 模式,模块被迫变小——$\Delta \nu$ 同时钉住门误差与模块规模,是全架构最核心的单一旋钮。

可行性判定

综合五本账(本词条估算,参数取当代文献值):(i) 门级——99%–99.7% 的 SWAP/CZ 在基线参数下即可达成,2/3 幂律对存储器劣化高度宽容;(ii) 查询级——$n=3$–$4$(N=8–16)的单谐振器演示在路线图级声子晶体存储(${\kappa }_{\gamma }\sim 1{\mathrm{s}}^{-1}$)下损耗贡献 1%–3%,加上门错误可期个位数百分比总失真,是最近期的实验目标;$n=10$ 则三墙并立——存储寿命须进秒级(Q≥10¹⁰ 且逆 Purcell 项被 GHz 级失谐或毫秒级 $T_1$ 压到 $1{\mathrm{s}}^{-1}$ 以下)、门错误须进 Hann 非空洞区($\epsilon \lesssim 2.5\times 10^{-3}$)、驱动加热与误共振随并行度累加。三条出路都在文献框架内:更高的声子 Q、更大的 $g_k$(换得同 $g_v$ 下更小的 $\lambda$)、更大的 $\alpha$ 或相干性更高的辅助量子比特(恒等式 $\gamma g_v/2\alpha {\xi }^2$)。此外多光子共振带须逐台器件标定避开,Stark 移须实时补偿。总体判定:机制全部成立、门级已达标,查询级可行性押注在存储器相干再进两个数量级上——这正是该提案把”量子信息住进声子、transmon 只当媒人”作为设计纲领的全部动机。

边界澄清与后续

该文同时给出 RAQM ≠ QRAM 的权威界定:只接收单一经典地址、返回对应量子比特的”random access quantum memory”一般不能执行叠加寻址操作,不应与 QRAM 混淆。局限在于:需要远长于门时间的声子相干(超导电路耦合是否从根本上限制声子相干是公开问题)、虚拟门慢、纯提案无整机实验。在五路硬件谱系(见 硬件平台全景)中,后续 Weiss 等 2024 把”高相干存储 + transmon 仅作辅助”的思想升级到 3D 腔,Wang 等 2025 则改用行程声波路由,分别对应推导四中驻波/行程两种损耗账的极化。

延伸阅读

• 超导腔 QRAM
• 声子路由 QRAM
• 噪声韧性的严格证明
• 边界:RAQM 与量子存储
• 评价指标
• bucket-brigade 架构