机制回顾:为什么 QRAM 电路可以被经典模拟
ECS 操作的形式定义
Hann 等把”能被经典计算机高效处理”的操作形式化为高效可计算稀疏(efficiently computable sparse, ECS)操作 Hann 2021:设寄存器共 个量子比特,若输入态可写成
其中基态集合 与振幅 都能在多项式时间内被经典计算机枚举/求值,且操作本身在经典机上多项式时间可实现,则称该操作对这类态是 ECS 的。
置换结构:稀疏度只减不增
QRAM 查询电路只含 SWAP/CSWAP。对计算基,SWAP 是对换 ;CSWAP 在控制位为 0 时恒等、为 1 时对换——无论控制位取什么值,每个门在计算基上都是一个置换 。 层电路的酉算子是置换的复合 ,仍是一个置换:
作用于稀疏态,逐振幅只需更新基态标签:
每操作成本 ,与总希尔伯特空间维数 无关。更关键的是:置换不分裂振幅——一个基态永远只映射到一个基态,不会变成叠加。因此稀疏度沿电路单调不增:。这是 QRAM 与一般量子电路的本质分界:一般电路里的 Hadamard 一步就把稀疏度翻倍, 层后稀疏度可达 ,ECS 性质荡然无存。
输入侧的稀疏性来自查询语义:地址处于叠加 ( 个地址),bus 与输出寄存器处于基态。于是任意时刻全态为
其中 是该分支的数据基态, 是路由寄存器的乘积态——BB 路由中每条分支只激活自己路径上的 个节点,其余处于等待态(桶队架构)。振幅数 、每分支需显式存储的路由自由度 ,单次模拟(一条噪声轨迹)的成本为 。
比特数账:稠密模拟为何早早出局
BB 装置的总比特数:数据 + 路由器 个 qutrit(按 2 量子比特编码) + 地址 + bus,合计
稠密态矢量需 字节(复双精度),密度矩阵需 字节——后者是前者的平方,而含噪模拟恰恰绕不开混态。数值代入(本词条估算):
- ():,态矢量 勉强可行;密度矩阵 ,直接出局。
- ():,态矢量 ,出局。
- :,双指数 已无法书写。
结论:含噪 QRAM 的模拟必须二选一——要么对噪声做蒙特卡洛采样,把混态拆成纯态轨迹;要么换掉”稠密振幅数组”这一表示本身。以下三节分别对应三条技术路线。
蒙特卡洛轨迹:Hann 框架的复杂度推导
Kraus 采样与无偏保真度估计量
任何 CPTP 信道可作 Kraus 分解 ,其中 、,可解释为”以概率 施加 “。整条电路的噪声等价于在全部噪声位置抽取一条 Kraus 链 ,每条链把理想演化接续成纯态
混态输出 正是这些轨迹按概率的混合。对纯目标态 (理想输出)的保真度做展开:
即单轨迹量 的期望(本词条推导;采样式注入 Kraus 序列的做法出自 Hann 2021)。 取值有界于 ,由 Hoeffding 不等式, 条轨迹以置信 达到加性精度 只需
数值算例:、 时 条——与 无关。这是”采样换维度”的第一重红利:混态的 维表示被 个 成本的纯态轨迹替换, 只由统计精度决定。注意保真度是”对角型”观测量(只问输出与目标的交叠),天然适配轨迹平均;轨迹法不产出 的非对角块——这正是它省掉的东西,也是它的能力边界。
每轨迹成本与总量账
每轨迹成本 (稀疏度 个振幅 × 深度 层的置换更新),总量 。文献锚点:原始模拟器做到 Hann 2021。本词条估算 、全叠加查询():每轨迹 次振幅更新,乘 得 次基本操作,单机秒级——与文献可及规模自洽。该模拟器验证的标度结论( 下三能级 、双能级斜率 、amplitude damping 观测斜率 1.86)见一般噪声下的韧性;局域化的实验对照见误差局域化。
相干误差怎么办
韧性定理覆盖任意 CPTP 信道,包括相干错误;轨迹法对它的处理反而最简单。一个相干错误 本身就是单元素 Kraus 分解(),按链采样时它退化为一次确定性的置换之外的酉更新——但注意这正是 ECS 假设发挥作用的地方:模拟器需要能对稀疏态执行这个酉。对 QRAM 电路,路由门的相干版本(如带相位误差的 CSWAP)仍作用在 个活跃振幅上;对等待态组件的相干漂移,则按 Hann 的分解归入坏分支、由好分支的加权比例吸收(证明见一般噪声下的韧性)。数值上文献以去极化与 amplitude damping 为代表信道验证标度;相干误差下 会叠加出系统性相位偏差,轨迹平均同样无偏,只是不再有”每次错误独立洗牌”的方差结构,收敛行为依赖误差的具体取向(此句为本词条方法学注记,文献未单独给相干情形的模拟数据)。
Clifford 路线:网络化 QRAM 的事件驱动模拟
协议的 Clifford 性与 账
自旋–光子网络的传态式 QRAM(TD-QRAM,每层 GHZ 态 + Bell 测量写地址)的理想协议部分——GHZ 制备( + CNOT)、CNOT 链、Bell 测量、Pauli 帧修正——全部落在 Clifford 群内 Bugalho 2023。按戈特斯曼–克尼尔定理(Gottesman–Knill theorem),这类电路可被经典高效模拟:纯态由其稳定子群(stabilizer group)唯一刻画, 比特态对应 个独立生成元,排成 的表(tableau);Clifford 门是对生成元的共轭置换,单门更新 ,Pauli 测量 ;内存 比特。与 彻底无关,这就是 NetSquid 路线先以 稳定子模拟跑完整个路由协议的依据。本词条估算:(对应下述有效区”约 量子比特”),表占 比特 ;即便 个协议门,总更新 比特操作也在单机小时级。
噪声的解析重构:一阶展开及其有效区
非 Clifford 的部分只剩”等待期间的退相干”与”含错 CNOT”。事件驱动仿真记录每个量子比特的退相干时长与含错 CNOT 位置,事后对稳定子态解析地施加局部信道。推导(本词条推导,教科书级展开):设每处噪声为小参数信道 。保真度对信道是线性的:
对 逐项展开并归纳,得
其中 是协议理想部分的失真。等待时长 的退相干贡献 ,含错 CNOT 贡献 ,故
不保真度按”每次噪声事件”线性记账,这正是事件驱动记录时长后可事后解析重构的原因。一阶展开要求总不保真度 ;当 个组件各自贡献 时即 ——与文献声明的有效区一致。数值: 个组件、单组件 时 已在边界; 则 仍可用。
文献定量结论
与全密度矩阵对比(5 层 TD-QRAM)偏差小于 1%,近似在 (约 量子比特)内有效 Bugalho 2023。确定性核–电子 CNOT 主导不保真度,访问 128 量子比特 QRAM 需 CNOT 保真度约 99.9% 才能使访问保真度超过 90%。一致性检查(本词条推导):若平均每存储量子比特承担 次确定性 CNOT,由上节记账 在 时给出 ——正是 99.9% 门槛的量级来源。据此提出的 TS-QRAM 改进:以概率性光子 CNOT(约 10 ns 量级,快于确定性核–电子 CNOT 约四个数量级)按类中继器二叉树顺序链接,混合方案把少量确定性节点集中于链接树高层(如 6.2% 确定性节点),在 、 ms、CNOT 成功率约 70% 时,1024 个存储单元保真度仍高于 60%(器件背景见自旋–光子 QRAM)。
稀疏态与子树剪枝:误差局域化的算法化
分支分解表示
中科大团队的 QRAM-Simulator 把上节”逐分支乘积结构”直接做成数据结构:按分支(地址条件下的唯一路由路径)重组稀疏态,每分支只存 个活跃路径节点 Wang 2025。内存 。本词条估算全叠加查询():
- : 个记录 ;
- : 个记录 ;
- : 个记录 (单机出局)。
的 335 MB 与文献”20 层树(百万量级节点)在 <1 GB 内存内模拟”同量级(差额是实现开销)——这构成对表示复杂度 的独立一致性验证。
子树包含准则及其证明
剪枝基于子树包含准则:故障 ( 层节点 )只影响地址区间
内的分支;区间外分支共享同一无噪树态、可整体跳过 Wang 2025。证明:二叉树上 层节点由 比特前缀标记;从 出发可达的叶恰是”前 位等于 “的地址。把前缀 读作 位二进制整数,该集合就是以 为起点、连续 个整数的区间。区间外的分支不经节点 ,其态演化与该故障无关——这是 Hann 误差局域化证明的算法化(误差局域化),把”错误限制在指数少分支”从定性引理变成可执行的地址区间测试。数值算例:、、 得 ,恰 个地址。
采样增益:一个伸缩求和
剪枝把蒙特卡洛轨迹的期望工作量也压低了。设单个节点在一次查询中承受的噪声应用总概率为 (多时间步按 折算),则每轨迹的期望”分支更新”数为(本词条推导)
每层贡献恰为 ——浅层故障”稀有而全局”、深层故障”常见而局部”,乘积与层无关。方差由浅层主导: 时一个根附近故障波及 个分支,这解释了为何实现上需要子树态共享/写时复制,而不是简单地逐故障重算。经验对照:128 核、512 GB 基准上,剪枝模式把噪声主导区的运行时间降至全模拟的约 20%、内存约 60% Wang 2025。该模拟器随即被用于误差过滤的大规模研究:高噪声或大地址尺寸下的抑制异常、把后选择概率纳入后的近确定性判据,见误差过滤。
少地址查询:期望工作量与 N 无关
误差过滤研究恰恰用得上”少地址”情形:EF 数值只查询 个地址。此时伸缩求和还能再压一层:被查地址集 均匀散布时,节点 的子树与 的期望交叠为 ,于是
期望工作量变为 ——与 无关,只随地址数线性走(本词条推导)。直觉:剪枝把”装置有 个组件”这一事实从成本账里彻底删除了,剩下的只有”查询真正碰到的组件”。这正是同一个模拟器既能跑 全叠加基准、又能扫误差过滤的尺寸–噪声参数网格的算力来源;两类负载的瓶颈从此分家——前者受内存 约束,后者受轨迹数与每轨迹 更新约束。
算法侧与电路侧:SparQSim、Di Matteo 核算
SparQSim:理想化 oracle 的寄存器级模拟
SparQSim 走另一方向:把 QRAM 作为理想数据加载模块内置于寄存器级稀疏态模拟器—— 直接访问经典内存,无噪声、不展开门电路 Sun 2025。基准上,20 比特 QRAM 电路(稀疏度 1 分支)只需 <10⁻⁴ s / 4.8 MB(Qiskit 为 0.1577 s / 221.7 MB),并支撑量子线性系统求解器的全流程模拟。注意口径:这里”20 比特”是数据装载电路的总比特数,不是 个存储单元——SparQSim 模拟的是”以 QRAM 为 oracle 的算法”,而非 QRAM 装置本身。值得引用的是其自我批评:“以容错技术实现大规模 QRAM 可能并不现实,因为 T 门带来的指数级开销”——算法社区把 QRAM 当作单位成本 oracle 的惯性,正是成本评估的盲区(见oracle 假设谱系)。
稀疏表示在这里也有明确失效面:干涉型电路会让分支数爆炸。29 比特 QFT 达 分支(恰为 量级),退化为 。本词条换算其内存反超:稠密态矢量 ,而逐分支存储 ——分支记录的每分支 开销使内存反超通用模拟器约 29 倍。
Di Matteo 核算:不追踪态矢量的”模拟”
与态模拟正交的是电路级资源核算:Di Matteo、Gheorghiu、Mosca 把多族 QRAM 查询电路(串行/并行 bucket-brigade、大深度小宽度、大宽度小深度、k-分层混合)完整嵌入表面码容错框架 Di Matteo 2020。它”能覆盖 “的技巧在于:核算只数逻辑比特数、T-count/T-depth 与表面码周期,是解析计数,不需要态矢量;电路族”功能如同 qRAM”由构造与分析确立,功能验证只需小规模语义测试。两者解耦后,计数可外推到态模拟不可及的规模。成本度量取 ;T-depth 为 1 的 Toffoli 分解恒记 7 个 T。并行 BB 族需 对 compute–uncompute Toffoli(多项式编码 QRAM词条的核算),本词条换算 T-count 。锚点数字:4 KB 表()在码周期 200 ns 下需 0.348 ms(折合 1740 个码周期);8 GB()需 物理比特换 2.13 ms(完整账见容错资源红线与电路级优化)。这些天文数字正是态模拟社区转向”物理级含噪 BB”(Hann、Wang 路线)而非”容错电路级”的动机:前者的研究对象是不纠错的物理装置,后者在实用尺寸下已被红线排除。
数值算例:各方法在 n=10/20/36 的资源对照
下表汇总各路线的代表性资源数(标”文献”者直接引自对应论文;标”估算/换算”者为本词条推导,方法与中间数字见上文各节)。注意两套口径: 通常指地址比特数( 个存储单元);SparQSim 行的”20 比特”是电路总比特数。
| 方法 | 口径 | |||
|---|---|---|---|---|
| 稠密态矢量() | 振幅,出局 | 估算 | ||
| 稠密密度矩阵 | 已需 288 PB | 出局 | 出局 | 估算 |
| ECS 稀疏/剪枝(全叠加, 字节) | ;20 层 <1 GB(文献) | ,单机出局 | 估算+[[sources/QR043 | |
| 蒙特卡洛轨迹数(,) | 同左,与 无关 | 同左 | 推导 | |
| 稳定子表( 比特,网络化) | : | — | :( 已临边界) | 估算+[[sources/QR038 |
| 寄存器级稀疏(QRAM 为理想 oracle) | — | <10⁻⁴ s / 4.8 MB(Qiskit:0.1577 s / 221.7 MB) | — | 文献 [[sources/QR044 |
| BB 电路 T-count | 换算自 [[sources/QR009 | |||
| 容错时空核算 | — | — | 物理比特 / 2.13 ms | 文献 [[sources/QR009 |
三点读法。其一,表的上半部与下半部回答不同问题:上半部问”含噪物理装置的查询保真度几何”,下半部问”若强行容错要付多少”——红线词条的结论是后者在实用尺寸不可行,于是上半部(免纠错路线)成为唯一活路,这正是 一般噪声下的韧性 的定理与它的模拟验证之间的分工。其二,ECS 稀疏行的幂律增长()与稠密行的双指数之间隔着整整一层指数,这个差距把可模拟规模从 推到 。其三,轨迹数行与 无关不是巧合:蒙特卡洛买的是统计精度,维度灾难由稀疏表示单独消化,两条账互不耦合。
模拟验证了哪些理论预言
QRAM 模拟器的产出不是数字本身,而是对理论清单的逐条对账:
- Hann 韧性界: 下三能级 BB 的 、双能级斜率 、amplitude damping 斜率 1.86,均在 内验证 Hann 2021(定理与证明见一般噪声下的韧性)。
- 误差局域化:剪枝模拟器的子树包含准则本身就是局域化引理的算法化 Wang 2025;实验侧,三层 BB 错误注入给出被查询路径与最远分支约 50 倍的保真度斜率差(Shen 2026,见误差局域化)。
- EF 判据的精化:理想 抑制在高噪声/大地址下的失效区、后选择概率 的 无关下界 、基失真 时几乎达不到 抑制——全部由剪枝模拟器在大规模上首次刻画 Wang 2025(见误差过滤)。
- 资源态 QRAM 的标度指数:嵌套单热编码的数值模拟给出 ,(去极化)或 (Z 偏置)Cesa 2026——模拟先于器件给出可检验的标度预言。
局限
稀疏优势依赖电路结构:干涉型电路(如 29 比特 QFT 达 分支)退化为 ,内存反超通用模拟器;网络化方法的振幅阻尼截断近似只在 内有效;剪枝模拟器的噪声信道沿用逻辑级近似。ECS 性质对含任意门的一般量子电路并不成立——它是 QRAM 路由结构的专属红利。此外:轨迹法只产出保真度一类对角型观测量,拿不到输出混态的非对角块;蒙特卡洛的 账在要求百分位以下精度时迅速变贵;全叠加查询下剪枝模拟器的内存 ,把无辅助技巧的模拟 ceiling 钉在 一带—— 的世界目前只有解析核算(电路级)能够触达,而这恰好也是”QRAM 是否值得建”这一问题的回答尺度(见成本口径)。