形式陈述
门的误差与实验的误差
给定按时间顺序执行的理想量子门 公理库 量子电路 Quantum circuit 用固定寄存器上的酉门、测量和经典控制表示有限量子操作序列,并明确基顺序、矩阵乘法与测后更新。 U 1 , … , U L ,编译后对应 V 1 , … , V L ,均已嵌入同一寄存器空间。总算子为 U = U L ⋯ U 1 、V = V L ⋯ V 1 。使用 Euclidean 诱导的算子范数 公理库 矩阵范数与诱导算子范数 Matrix norm · Induced matrix norm · Operator norm of a matrix 用诱导范数和常用可计算矩阵范数度量线性映射的放大能力,并区分算子范数、Frobenius 范数与谱半径。
‖ A ‖ o p = sup ‖ ψ ‖ 2 = 1 ‖ A ψ ‖ 2 , δ j = ‖ U j − V j ‖ o p . 本页完整证明:对任意输入,包括与不操作的参考系统纠缠的输入,任意共同最终 POVM 的输出分布 p , q 满足
(1) TV ( p , q ) ≤ D ( ρ U , ρ V ) ≤ ‖ U − V ‖ o p ≤ ∑ j = 1 L δ j . 其中 D 是迹距离 公理库 迹距离 Trace distance · Quantum trace distance 迹距离是两个密度算子之差的迹范数的一半,恰好刻画单份量子态在最优测量下的可区分程度。 。因此总容差 ε 可以分配为 ∑ j δ j ≤ ε 。这是最坏输入、最坏测量的保证,不是特定测量的实际误差等式。
固定有限字母表
单 qubit 令 T = diag ( 1 , e i π / 4 ) ,H 为 Hadamard,取
H ― = i H , T ― = e − i π / 8 T = R z ( π / 4 ) , G = { i H , − i H , T ― , T ― † } ⊂ S U ( 2 ) . 这个集合有限且对逆封闭;特别是 ( i H ) † = − i H ,不能因 H † = H 而漏掉逆。其生成群在 S U ( 2 ) 中稠密,这里引用 H , T 的通用性结果。[1, §§2–3] 对通常物理门,每个这样的字对应 H , T , T † 的字并只差整个字的标量相位;T † = T 7 也可改成常数开销。
引用的 Solovay–Kitaev 结论是:对固定、逆封闭、稠密生成 S U ( 2 ) 的 G ,可用
O ( [ 1 + log ( 1 / δ ) ] c ) , c = log 5 log ( 3 / 2 ) ≈ 3.9694 个字母逼近任意目标 S U ( 2 ) 门到算子误差 δ ,0 < δ < 1 。[2, Theorem 1、§3] 算法先固定一个常精度网表,再递归细化。这里不重证稠密性或交换子细化定理;本页证明的是旋转参数、已获得近似门及最终实验之间的全部误差转换。常数依赖 G 与预处理网表;逻辑门数不等于有限精度经典编译的位复杂度。
直觉
一次旋转多转了一点,不一定让当前测量结果明显改变:计算基测量甚至可能完全看不到 z 轴相位误差。但后续换基、干涉或参考纠缠可以显露差异,所以用算子范数先控制所有输入,再转成可观察概率。
多个小误差最坏可能朝同一方向积累。逐门替换、每次只比较一个差项,便得到总和预算,不需要假设误差独立或随机抵消。
例子与边界
用一个 T 门近似 R z ( π / 3 )
对实单位向量 n ,令 A = n x X + n y Y + n z Z 。Pauli 反对易关系给出 A 2 = I ,故指数级数分出奇偶项后有
R n ( θ ) = e − i θ A / 2 = cos ( θ / 2 ) I − i sin ( θ / 2 ) A . 特别是
R z ( θ ) = ( e − i θ / 2 0 0 e i θ / 2 ) , R y ( θ ) = ( cos ( θ / 2 ) − sin ( θ / 2 ) sin ( θ / 2 ) cos ( θ / 2 ) ) . 同一轴上相差 Δ = θ − θ ~ 时,酉不变性将误差化成 ‖ R n ( Δ ) − I ‖ 。A 的特征值为 ± 1 ,所以差矩阵的两个奇异值都是 2 | sin ( Δ / 4 ) | :
(2) ‖ R n ( θ ) − R n ( θ ~ ) ‖ o p = 2 | sin ( Δ / 4 ) | ≤ | Δ | / 2. 这固定了矩阵相位;当 Δ = 2 π 时范数为 2 ,两酉信道却相同。
用物理 T 门时,比较的相位对齐代表为 R z ( π / 4 ) ,故对目标 R z ( π / 3 ) ,
δ = 2 sin ( π / 48 ) ≈ 0.1308062585 . 输入 | + ⟩ ,旋转后的 X 基测量满足 p ( + ∣ θ ) = ( 1 + cos θ ) / 2 。按 + , − 次序,理想与近似分布是
p = ( 3 / 4 , 1 / 4 ) , q = ( ( 2 + 2 ) / 4 , ( 2 − 2 ) / 4 ) . 于是
TV ( p , q ) = 2 − 1 4 ≈ 0.1035533906 , D ( ρ U , ρ V ) = sin ( π / 24 ) ≈ 0.1305261922 < δ . 纯态重叠给出的平方保真度 公理库 量子态保真度与最优纯化重叠 Quantum fidelity · Uhlmann fidelity · 量子态忠实度 · Uhlmann 定理 用根保真度比较混态,并通过环境酉优化证明它恰好等于固定一方纯化后的最大重叠,完整计算一对非对易混态。 为 F = cos 2 ( π / 24 ) ≈ 0.9829629131 。三种数值分别是当前测量偏差、最佳测量可区分度和统一算子上界,不能互换。
三个门的实际误差与预算
连续三个 R z ( π / 3 ) 合成 R z ( π ) ;以三个 T 替代后,相位对齐的总门是 R z ( 3 π / 4 ) 。因此实际总范数差为
2 sin ( π / 16 ) ≈ 0.3901806440 ≤ 3 δ ≈ 0.3924187754 . 对 | + ⟩ 再测 X ,理想 + 概率为零,近似侧为 ( 2 − 2 ) / 4 ≈ 0.1464466094 。总预算虽然偏保守,却不需要预先知道输入和测量。
整体相位放入控制分支后不再整体
独立操作可以按
d p h ( U , V ) = min ϕ ‖ U − e i ϕ V ‖ o p 比较,因为标量相位不改变输出密度算子。但受控门满足
c ( e i ϕ U ) = ( P ( ϕ ) ⊗ I ) c U , P ( ϕ ) = diag ( 1 , e i ϕ ) , 多出的相位已落在控制线。最简单地,d p h ( I , − I ) = 0 ,但 c I = I 、c ( − I ) = Z ⊗ I ;对 | + ⟩ | 0 ⟩ ,它们给出 | + ⟩ | 0 ⟩ 与 | − ⟩ | 0 ⟩ ,测控制线 X 的总变差是 1 。
因此不能先丢弃目标门相位,再免费给它加控制。若先以相位固定的矩阵比较,则块对角结构恰有 ‖ c U − c V ‖ = ‖ U − V ‖ ,不会丢掉这一差别。
推论与应用
任意单qubit门的旋转参数
任意 U ∈ U ( 2 ) 可选 e 2 i χ = det U ,使 e − i χ U ∈ S U ( 2 ) 。后者具有形式
V = ( a − b ― b a ― ) , | a | 2 + | b | 2 = 1. 写 a = | a | e i ϕ a 、b = | b | e i ϕ b ;系数为零时相应辐角任选。取
β = 2 atan 2 ( | b | , | a | ) , α = ϕ b − ϕ a , γ = − ϕ a − ϕ b . 直接乘矩阵得到
R z ( α ) R y ( β ) R z ( γ ) = ( e − i ( α + γ ) / 2 cos ( β / 2 ) − e − i ( α − γ ) / 2 sin ( β / 2 ) e i ( α − γ ) / 2 sin ( β / 2 ) e i ( α + γ ) / 2 cos ( β / 2 ) ) = V . 因零系数所乘项消失,这个选法覆盖 a = 0 或 b = 0 。若三角分别舍入 Δ α , Δ β , Δ γ ,式 (2) 加上下文乘积界给出总误差至多 ( | Δ α | + | Δ β | + | Δ γ | ) / 2 。丢弃的 e i χ 仅在这个完整操作作为普通局部门时可忽略。
从算子误差到含参考系统的迹距离
先证对任意有限维参考空间,‖ ( U − V ) ⊗ I R ‖ = ‖ U − V ‖ 。将向量按参考正交基写成 Ψ = ∑ r ψ r ⊗ | r ⟩ ,平方范数逐项相加,给出“≤ ”;取达到单系统最大值的向量与任意参考单位向量的乘积,给出“≥ ”。
令 A = U ⊗ I R 、B = V ⊗ I R ,δ = ‖ U − V ‖ 。对任意密度算子 ρ ,
A ρ A † − B ρ B † = ( A − B ) ρ A † + B ρ ( A † − B † ) . 对 ρ = ∑ k p k | ψ k ⟩ ⟨ ψ k | ,秩一算子的迹范数为两侧向量范数之积,从而
‖ C ρ D ‖ 1 ≤ ∑ k p k ‖ C ψ k ‖ ‖ D † ψ k ‖ ≤ ‖ C ‖ ‖ D ‖ . 上面两个差项的迹范数因此各至多 δ ,除以迹距离定义中的二,得到 D ( ρ U , ρ V ) ≤ δ 。证明没有参考维数因子,也没有纯态假设。把 V 换成最优相位对齐代表,同样得到 d p h 的版本。
对纯输出,D 2 = 1 − F ,所以 δ ≤ 1 时还得 F ≥ 1 − δ 2 。这一保真度式不能不加说明地当作任意两个混态的等式。
从量子态到完整测量分布
对 POVM 效应 E a ,令 p a = Tr ( E a ρ ) 、q a = Tr ( E a σ ) ,并取 S = { a : p a ≥ q a } 。总概率相等,所以
TV ( p , q ) = ∑ a ∈ S ( p a − q a ) = Tr ( E S ( ρ − σ ) ) , E S = ∑ a ∈ S E a . 0 ⪯ E S ⪯ I 。把 ρ − σ 分为正负谱部分,两者迹都等于 D ( ρ , σ ) ,于是最后一项不超过 D 。故整个测量分布、任何单个结果或结果事件的概率误差都受同一个 δ 控制。共同后续信道的收缩性以及经典后处理也保持此界。
逐门替换与字长
空乘积约定为恒等,以下恒等式的中间项逐项相消:
U L ⋯ U 1 − V L ⋯ V 1 = ∑ j = 1 L U L ⋯ U j + 1 ( U j − V j ) V j − 1 ⋯ V 1 . 两侧的酉乘积范数为一,用三角不等式即得 ‖ 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 , e i θ ) ,一个相位安全的分解是
c P ( θ ) = ( P ( θ / 2 ) ⊗ P ( θ / 2 ) ) CNOT ( I ⊗ P ( − θ / 2 ) ) CNOT . 按矩阵从右到左作用,计算基 | a , b ⟩ 的总相位为 θ ( a + b − ( a ⊕ b ) ) / 2 = θ a b ,末尾两次 CNOT 恢复基标签,故对所有叠加也成立。再将三个普通 P 门换成相应 R z 代表,只改变这整个两qubit门的一份整体相位。
因此精确QFT电路 公理库 量子 Fourier 变换 Quantum Fourier transform · QFT 从有限 Fourier 矩阵的酉性与二进制分解推导 Hadamard、受控相位门和位反转电路,并区分振幅变换与经典输出。 的每个受控相位可由三次单qubit合成与两个精确 CNOT 实现。m ≥ 2 时共有 M = 3 m ( m − 1 ) / 2 个待合成旋转;逐个用 ε / M 精度,总测量误差至多 ε ,门数为
O ( m 2 [ 1 + log ( m 2 / ε ) ] 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:预处理、字长递推、共享输出和算术精度边界。本页旋转误差、测量预算和数值例均直接展开,未重证交换子合成。