形式陈述
取 N = 2 n 维Hermitian矩阵 H = H † 、时间 t ∈ R 和精度 ε > 0 。量子Hamiltonian模拟的目标是构造量子电路 公理库 量子电路 Quantum circuit 用固定寄存器上的酉门、测量和经典控制表示有限量子操作序列,并明确基顺序、矩阵乘法与测后更新。 V ,使
(1) ‖ V − e − i H t ‖ ≤ ε . 这里采用谱算子范数 公理库 矩阵范数与诱导算子范数 Matrix norm · Induced matrix norm · Operator norm of a matrix 用诱导范数和常用可计算矩阵范数度量线性映射的放大能力,并区分算子范数、Frobenius 范数与谱半径。 ,因此误差对所有归一化输入统一成立;张量上任意外部参考的恒等算子后仍有同一界。若电路使用辅助位,还要说明辅助是否返回零态,或改用相应等距映射/通道误差,不能只比较某个后选择分支。
任务的输入不是一句“给定矩阵”就说清了。常见接口有两类:[1, Problems 1–2]
已知可模拟项 :H = ∑ j = 1 m H j ,并给出实现 e − i H j τ 的电路及其门成本,任意所需 τ 与精度均须计费
稀疏访问 :每行至多 d 个非零条目,给出可相干查询第 ℓ 个非零位置的oracle,以及查询条目数值的oracle;数值位宽、逆查询和支持索引约定都应写明
例如可把稀疏位置查询规定为
O F | i , ℓ , z ⟩ = | i , ℓ , z ⊕ f ( i , ℓ ) ⟩ , 条目查询为 O H | i , j , z ⟩ = | i , j , z ⊕ enc ( H i j ) ⟩ 。它们是对叠加输入有效的可逆接口。显式经典矩阵列表不自动提供这样的单位成本oracle;建立数据结构或编译访问电路的成本也属于问题。
直觉
H 描述瞬时变化,e − i H t 描述经历时间 t 后的整体演化。量子模拟并不要求把这个 N × N 矩阵的全部条目经典打印出来,而是让未知量子态受到相同作用。
维数可能指数大,但只有当 H 具有可用的局部结构或高效访问接口时,电路才可能比逐项处理整个矩阵更省。把一个庞大数据表藏在“调用一次H”里,会把真正昂贵的输入工作从账本中删掉。
图片加载失败
例子与边界
一个可直接复算的模拟:H等于X加Z
令
H = X + Z = ( 1 1 1 − 1 ) . 由Pauli 代数 公理库 Pauli 算子与 Pauli 群 Pauli operator · Pauli group · Pauli string · 泡利算子 · 泡利群 保留整体相位的 Pauli 张量积群,以二进制标签计算乘法、Hermitian 条件和对易符号。 的 X 2 = Z 2 = I 及 X Z + Z X = 0 ,有 H 2 = 2 I ,所以
(2) e − i t H = cos ( 2 t ) I − i sin ( 2 t ) 2 ( X + Z ) . 这条小矩阵恒等式只是用来检验电路,不是把一般量子模拟改成经典矩阵指数计算。
把复振幅视作实部与虚部组成的 2 N 维实向量,Schrödinger 方程就是常系数线性 ODE,其两个子流为相应酉演化。取整数 r ≥ 1 、Δ = t / r ,调用已有Lie–Trotter分裂 公理库 Lie–Trotter 与 Strang 分裂 Lie-Trotter splitting · Strang splitting 把可解子流按一阶或对称二阶次序复合,用非交换误差解释分裂的精度与结构保持。 :
V r = ( e − i Δ Z e − i Δ X ) r . 每片在时间上先执行 X 演化,再执行 Z 演化;矩阵乘法右侧先作用。两门分别是 R X ( 2 Δ ) 、R Z ( 2 Δ ) ,因此调用局部演化恰为 2 r 次。
对Hermitian两项,Duhamel积分中的演化均为酉,范数为1,单片误差可界为
‖ e − i Δ ( X + Z ) − e − i Δ Z e − i Δ X ‖ ≤ Δ 2 2 ‖ [ X , Z ] ‖ = Δ 2 . 再对 r 个酉片逐项替换,得到
(3) ‖ V r − e − i t H ‖ ≤ t 2 r . 这里沿用分裂方法的基本机制,只利用酉性把增长因子去掉;一般多项、高阶公式的对易子分析见文献 [2]。
数值误差与理论保证不是同一个数
当 t = 1 时,直接用式 (2) 计算可得:
r
局部演化次数
实际算子误差,约
式(3)上界
1
2
0.799214
1
2
4
0.362410
0.5
4
8
0.176261
0.25
10
20
0.069951
0.1
上界不是实际误差的等式。交换片内两门次序会改变近似电路;二者都收敛到同一目标,但有限 r 时不能直接认为矩阵相同。
如果要求总误差不超过 ε ,可先给分裂分配 ε / 2 ,选择 r ≥ 2 t 2 / ε 。若每个局部旋转又只有误差 η ,逐门替换额外带来至多 2 r η ;选 η ≤ ε / ( 4 r ) 才完成这份预算。固定门集合成 公理库 固定门集近似与测量误差预算 Finite gate approximation · Quantum gate synthesis error budget 从单qubit旋转的算子误差推到任意参考系统上的测量总变差,逐门分配合成预算,并处理受控门中的相对相位。 的实际门数需按该精度再计,不能把任意实角旋转永久当作精确免费门。
推论与应用
一个模拟保证对应哪些可观察结果
式 (1) 使任意输入的输出纯态向量差至多 ε ,进而迹距离至多 ε ;混合输入和带参考输入可通过纯化 公理库 量子态纯化 Quantum state purification · Purification of a quantum state 量子态纯化将一个混态表示成较大系统中纯态的约化态,最小辅助空间维数等于原态的秩。 得到同一测量分布保证。因此一个具体测量事件的概率差至多 ε 。
这不意味着一次运行能知道全部振幅。若最终需要某个可观察量的平均值,还要准备多份输入、重复模拟和测量;这些采样次数与模拟一次的门成本分别报告。Hadamard测试 公理库 Hadamard 测试 Hadamard test 逐振幅推导受控酉期望值的实部和虚部读出,固定S逆门的相位约定,并给出样本数、输入制备和黑盒控制访问的完整成本。 给出一个具体的期望值读出接口。
时间、尺度和黑盒幂都要计费
对常数 c > 0 ,( c H , t / c ) 给同一演化,所以自然参数包含 ‖ H ‖ | t | ,或某种访问构造的更大归一化尺度。稀疏度 d 、条目位宽和目标精度也会影响查询与门数。
相位估计 公理库 量子相位估计 Quantum phase estimation · QPE 用受控酉幂将特征相位写入控制寄存器,再以逆 Fourier 变换读出,推导精确情形、有限概率分布与实际查询成本。 若使用 U = e − i H τ ,其 U 2 j 对应时间 2 j τ ,不是因为写成一个幂就只算一次基本模拟。若接口真的额外提供高次幂,应把它列为更强的输入承诺。
一般黑盒Hamiltonian不能任意快进。不可快进定理 公理库 量子模拟的不可快进定理 No-fast-forwarding theorem 以范数一的加权路径Hamiltonian把隐藏输入奇偶传播到终点,证明每次稀疏访问只需常数次bit查询,从而得到最坏情形线性时间查询下界。 用范数有界的稀疏实例把输入奇偶编码到长时间传播中,证明最坏情形需要与时间成比例的查询。它不排除像已知对角矩阵这样的特殊易模拟族。
参考资料