Skip to content

方法Method

固定门集近似与测量误差预算

Finite gate approximation · Quantum gate synthesis error budget

从单qubit旋转的算子误差推到任意参考系统上的测量总变差,逐门分配合成预算,并处理受控门中的相对相位。

形式陈述 ​

门的误差与实验的误差 ​

给定按时间顺序执行的理想量子门 U1,…,UL,编译后对应 V1,…,VL,均已嵌入同一寄存器空间。总算子为 U=UL⋯U1、V=VL⋯V1。使用 Euclidean 诱导的算子范数

‖A‖op=sup‖ψ‖2=1‖Aψ‖2,δj=‖Uj−Vj‖op.

本页完整证明:对任意输入,包括与不操作的参考系统纠缠的输入,任意共同最终 POVM 的输出分布 p,q 满足

(1)TV(p,q)≤D(ρU,ρV)≤‖U−V‖op≤∑j=1Lδj.

其中 D 是迹距离。因此总容差 ε 可以分配为 ∑jδj≤ε。这是最坏输入、最坏测量的保证,不是特定测量的实际误差等式。

固定有限字母表 ​

单 qubit 令 T=diag(1,eiπ/4),H 为 Hadamard,取

H―=iH,T―=e−iπ/8T=Rz(π/4),G={iH,−iH,T―,T―†}⊂SU(2).

这个集合有限且对逆封闭;特别是 (iH)†=−iH,不能因 H†=H 而漏掉逆。其生成群在 SU(2) 中稠密,这里引用 H,T 的通用性结果。[1, §§2–3] 对通常物理门,每个这样的字对应 H,T,T† 的字并只差整个字的标量相位;T†=T7 也可改成常数开销。

引用的 Solovay–Kitaev 结论是:对固定、逆封闭、稠密生成 SU(2) 的 G,可用

O([1+log⁡(1/δ)]c),c=log⁡5log⁡(3/2)≈3.9694

个字母逼近任意目标 SU(2) 门到算子误差 δ,0<δ<1。[2, Theorem 1、§3] 算法先固定一个常精度网表,再递归细化。这里不重证稠密性或交换子细化定理;本页证明的是旋转参数、已获得近似门及最终实验之间的全部误差转换。常数依赖 G 与预处理网表;逻辑门数不等于有限精度经典编译的位复杂度。

直觉

一次旋转多转了一点,不一定让当前测量结果明显改变:计算基测量甚至可能完全看不到 z 轴相位误差。但后续换基、干涉或参考纠缠可以显露差异,所以用算子范数先控制所有输入,再转成可观察概率。

多个小误差最坏可能朝同一方向积累。逐门替换、每次只比较一个差项,便得到总和预算,不需要假设误差独立或随机抵消。

例子与边界

用一个 T 门近似 Rz(π/3) ​

对实单位向量 n,令 A=nxX+nyY+nzZ。Pauli 反对易关系给出 A2=I,故指数级数分出奇偶项后有

Rn(θ)=e−iθA/2=cos⁡(θ/2)I−isin⁡(θ/2)A.

特别是

Rz(θ)=(e−iθ/200eiθ/2),Ry(θ)=(cos⁡(θ/2)−sin⁡(θ/2)sin⁡(θ/2)cos⁡(θ/2)).

同一轴上相差 Δ=θ−θ~ 时,酉不变性将误差化成 ‖Rn(Δ)−I‖。A 的特征值为 ±1,所以差矩阵的两个奇异值都是 2|sin⁡(Δ/4)|:

(2)‖Rn(θ)−Rn(θ~)‖op=2|sin⁡(Δ/4)|≤|Δ|/2.

这固定了矩阵相位;当 Δ=2π 时范数为 2,两酉信道却相同。

用物理 T 门时,比较的相位对齐代表为 Rz(π/4),故对目标 Rz(π/3),

δ=2sin⁡(π/48)≈0.1308062585.

输入 |+⟩,旋转后的 X 基测量满足 p(+∣θ)=(1+cos⁡θ)/2。按 +,− 次序,理想与近似分布是

p=(3/4,1/4),q=((2+2)/4,(2−2)/4).

于是

TV(p,q)=2−14≈0.1035533906,D(ρU,ρV)=sin⁡(π/24)≈0.1305261922<δ.

纯态重叠给出的平方保真度为 F=cos2⁡(π/24)≈0.9829629131。三种数值分别是当前测量偏差、最佳测量可区分度和统一算子上界,不能互换。

三个门的实际误差与预算 ​

连续三个 Rz(π/3) 合成 Rz(π);以三个 T 替代后,相位对齐的总门是 Rz(3π/4)。因此实际总范数差为

2sin⁡(π/16)≈0.3901806440≤3δ≈0.3924187754.

对 |+⟩ 再测 X,理想 + 概率为零,近似侧为 (2−2)/4≈0.1464466094。总预算虽然偏保守,却不需要预先知道输入和测量。

整体相位放入控制分支后不再整体 ​

独立操作可以按

dph(U,V)=minϕ‖U−eiϕV‖op

比较,因为标量相位不改变输出密度算子。但受控门满足

c(eiϕU)=(P(ϕ)⊗I)cU,P(ϕ)=diag(1,eiϕ),

多出的相位已落在控制线。最简单地,dph(I,−I)=0,但 cI=I、c(−I)=Z⊗I;对 |+⟩|0⟩,它们给出 |+⟩|0⟩ 与 |−⟩|0⟩,测控制线 X 的总变差是 1。

因此不能先丢弃目标门相位,再免费给它加控制。若先以相位固定的矩阵比较,则块对角结构恰有 ‖cU−cV‖=‖U−V‖,不会丢掉这一差别。

推论与应用

任意单qubit门的旋转参数 ​

任意 U∈U(2) 可选 e2iχ=det⁡U,使 e−iχU∈SU(2)。后者具有形式

V=(a−b―ba―),|a|2+|b|2=1.

写 a=|a|eiϕa、b=|b|eiϕb;系数为零时相应辐角任选。取

β=2atan2(|b|,|a|),α=ϕb−ϕa,γ=−ϕa−ϕb.

直接乘矩阵得到

Rz(α)Ry(β)Rz(γ)=(e−i(α+γ)/2cos⁡(β/2)−e−i(α−γ)/2sin⁡(β/2)ei(α−γ)/2sin⁡(β/2)ei(α+γ)/2cos⁡(β/2))=V.

因零系数所乘项消失,这个选法覆盖 a=0 或 b=0。若三角分别舍入 Δα,Δβ,Δγ,式 (2) 加上下文乘积界给出总误差至多 (|Δα|+|Δβ|+|Δγ|)/2。丢弃的 eiχ 仅在这个完整操作作为普通局部门时可忽略。

从算子误差到含参考系统的迹距离 ​

先证对任意有限维参考空间,‖(U−V)⊗IR‖=‖U−V‖。将向量按参考正交基写成 Ψ=∑rψr⊗|r⟩,平方范数逐项相加,给出“≤”;取达到单系统最大值的向量与任意参考单位向量的乘积,给出“≥”。

令 A=U⊗IR、B=V⊗IR,δ=‖U−V‖。对任意密度算子 ρ,

AρA†−BρB†=(A−B)ρA†+Bρ(A†−B†).

对 ρ=∑kpk|ψk⟩⟨ψk|,秩一算子的迹范数为两侧向量范数之积,从而

‖CρD‖1≤∑kpk‖Cψk‖‖D†ψk‖≤‖C‖‖D‖.

上面两个差项的迹范数因此各至多 δ,除以迹距离定义中的二,得到 D(ρU,ρV)≤δ。证明没有参考维数因子,也没有纯态假设。把 V 换成最优相位对齐代表,同样得到 dph 的版本。

对纯输出,D2=1−F,所以 δ≤1 时还得 F≥1−δ2。这一保真度式不能不加说明地当作任意两个混态的等式。

从量子态到完整测量分布 ​

对 POVM 效应 Ea,令 pa=Tr(Eaρ)、qa=Tr(Eaσ),并取 S={a:pa≥qa}。总概率相等,所以

TV(p,q)=∑a∈S(pa−qa)=Tr(ES(ρ−σ)),ES=∑a∈SEa.

0⪯ES⪯I。把 ρ−σ 分为正负谱部分,两者迹都等于 D(ρ,σ),于是最后一项不超过 D。故整个测量分布、任何单个结果或结果事件的概率误差都受同一个 δ 控制。共同后续信道的收缩性以及经典后处理也保持此界。

逐门替换与字长 ​

空乘积约定为恒等,以下恒等式的中间项逐项相消:

UL⋯U1−VL⋯V1=∑j=1LUL⋯Uj+1(Uj−Vj)Vj−1⋯V1.

两侧的酉乘积范数为一,用三角不等式即得 ‖U−V‖≤∑jδj,补齐式 (1)。普通局部门若各有相位对齐因子,它们相乘成为整个电路的一份标量相位,仍可忽略;这个说法没有授权忽略受控目标内部的相位。

若只有 M≥1 个单qubit门需要近似,其余门精确实现,可各分配 δ=ε/M。由前面明确引用的合成定理,总字长为

O(M[1+log⁡(M/ε)]c),0<ε<1.

这是使用固定编译字母表的逻辑门数。Dawson–Nielsen 的更小经典运行时间指数使用共享指针表示输出,逐门展开则至少支付输出长度;其分析还明确未计有限精度算术成本,[2, p. 7] 因而这里不把该字长界冒称为完整编译位复杂度。

QFT中的受控旋转如何落实 ​

对 P(θ)=diag(1,eiθ),一个相位安全的分解是

cP(θ)=(P(θ/2)⊗P(θ/2))CNOT(I⊗P(−θ/2))CNOT.

按矩阵从右到左作用,计算基 |a,b⟩ 的总相位为 θ(a+b−(a⊕b))/2=θab,末尾两次 CNOT 恢复基标签,故对所有叠加也成立。再将三个普通 P 门换成相应 Rz 代表,只改变这整个两qubit门的一份整体相位。

因此精确QFT电路的每个受控相位可由三次单qubit合成与两个精确 CNOT 实现。m≥2 时共有 M=3m(m−1)/2 个待合成旋转;逐个用 ε/M 精度,总测量误差至多 ε,门数为

O(m2[1+log⁡(m2/ε)]c).

此外 m 个 H、受控相位中的 m(m−1) 个 CNOT,以及位反转的 3⌊m/2⌋ 个 CNOT 都可精确计费。m=1 时只需 H,没有合成项。这完成从可变角度电路到固定有限门集的预算;物理路由、器件噪声与纠错资源需另建模型。

参考资料
  • [1] P. Oscar Boykin, Tal Mor, Matthew Pulver, Vwani Roychowdhury, Farrokh Vatan, On Universal and Fault-Tolerant Quantum Computing,v1,1999,§§2–3,pp. 3–5:旋转表示与 H,T 的稠密性。本文只引用门集通用性,不承接其物理容错结论。
  • [2] Christopher M. Dawson and Michael A. Nielsen, The Solovay–Kitaev algorithm,v2,2005,§2 Definition 1,pp. 2–3:字母表与范数;Theorem 1,p. 4;§3,pp. 5–7:预处理、字长递推、共享输出和算术精度边界。本页旋转误差、测量预算和数值例均直接展开,未重证交换子合成。
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系