[论文解读] Quantum algorithm for time-dependent differential equations using Dyson series

Dominic W. Berry, Pedro C. S. Costa|arXiv (Cornell University)|Dec 7, 2022

Quantum Computing Algorithms and Architecture参考文献 17被引用 7

一句话总结

本文提出了一种利用Dyson级数求解时变线性微分方程的量子算法，通过块编码将解编码为量子态，并利用最优量子线性系统求解器进行求解。该算法在误差和导数依赖上实现对数尺度，复杂度在演化时间上线性增长，在精度上对数增长，相较于经典方法在高维情况下实现了指数加速。

ABSTRACT

Time-dependent linear differential equations are a common type of problem that needs to be solved in classical physics. Here we provide a quantum algorithm for solving time-dependent linear differential equations with logarithmic dependence of the complexity on the error and derivative. As usual, there is an exponential improvement over classical approaches in the scaling of the complexity with the dimension, with the caveat that the solution is encoded in the amplitudes of a quantum state. Our method is to encode the Dyson series in a system of linear equations, then solve via the optimal quantum linear equation solver. Our method also provides a simplified approach in the case of time-independent differential equations.

研究动机与目标

• 开发一种高效求解时变线性微分方程的量子算法，实现改进的复杂度标度。
• 解决使用量子方法求解此类方程的挑战，特别是当经典方法在维度增加时性能显著下降时。
• 利用Dyson级数实现时间有序演化，并将其编码为量子线性系统以实现高效求解。
• 实现对误差和导数的对数依赖，相较于先前方法简化了复杂度分析。
• 提供一个统一且简化的框架，适用于时变与时不变两种情形。

提出的方法

• 利用Dyson级数对系数矩阵A(t)的时间有序演化进行齐次解与特解的分解。
• 通过块编码将Dyson级数编码为块矩阵结构，类似于哈密顿量模拟中使用的方法。
• 利用Costa等人（2021）提出的最优量子线性系统求解器求解由Dyson级数导出的线性方程组。
• 使用酉操作U_A、U_b和U_x分别编码矩阵A(t)、非齐次项b(t)和初始态x₀。
• 应用幅度放大以高保真度制备最终的量子态，其调用次数在反误差上对数增长。
• 该方法避免了要求A(t)可对角化，并利用对数范数简化复杂度界。

实验结果

研究问题

• RQ1 Dyson级数能否被高效编码为量子线性系统，以实现对误差的对数尺度？
• RQ2 求解时变微分方程的复杂度如何随演化时间T、误差ε以及A(t)和b(t)的导数变化？
• RQ3 能否将哈密顿量模拟中使用的块编码方法适配到时变线性ODE中，以实现更优的复杂度分析？
• RQ4 该量子算法中对U_A、U_b和U_x的调用次数如何标度？与先前工作相比有何差异？
• RQ5 该方法能否简化时不变情形的复杂度分析，后者是时变形式的特例？

主要发现

• 算法对U_b和U_x的调用次数为O(ℛλT log(1/ε))，对U_A的调用次数为O(ℛλT log(1/ε) log(λT/ε))，额外门操作的复杂度为O(ℛλT log(1/ε) log²(λT/ε))。
• 误差依赖为1/ε的对数形式，这是量子算法的最优形式，相较于多项式或线性依赖有显著改进。
• 该方法在维度上相对于经典求解器实现了指数加速，同时保持对误差和A(t)、b(t)导数的对数依赖。
• 复杂度在总演化时间T和块编码参数λ上线性增长，在精度上对数增长。
• 通过避免块矩阵中的额外行，并采用类似哈密顿量模拟的基于Dyson级数的块编码，相较于先前工作简化了分析。
• 该方法还为时不变情形提供了简化的复杂度表达式，保持对数误差标度，并在系统参数依赖关系上实现更清晰的表达。