[论文解读] Towards Quantum Simulation of Chemical Dynamics with Prethreshold Superconducting Qubits

P. C. Stancil, Hao You+3 more|arXiv (Cornell University)|Jan 30, 2016

Quantum Information and Cryptography参考文献 1被引用 5

一句话总结

本论文提出在单激发子空间（SES）中使用预阈值超导量子比特来模拟低能化学动力学，将半经典的分子轨道耦合理论（SCMOCC）方法扩展至常速度近似之外。通过引入多通道势能平均化与埃伦费斯特定时对称化，实现了在低至1 eV/u能量下的精确模拟，借助经典计算机上的预处理自适应龙格-库塔传播，实现了220倍的加速，展示了未来对大规模分子碰撞体系进行量子模拟的可行性。

ABSTRACT

The single excitation subspace (SES) method for universal quantum simulation is investigated for a number of diatomic molecular collision complexes. Assuming a system of $n$ tunably-coupled, and fully-connected superconducting qubits, computations are performed in the $n$-dimensional SES which maps directly to an $n$-channel collision problem within a diabatic molecular wave function representation. Here we outline the approach on a classical computer to solve the time-dependent Schrödinger equation in an $n$-dimensional molecular basis - the so-called semiclassical molecular-orbital close-coupling (SCMOCC) method - and extend the treatment beyond the straight-line, constant-velocity approximation which is restricted to large kinetic energies ($\gtrsim 0.1$ keV/u). We explore various multichannel potential averaging schemes and an Ehrenfest symmetrization approach to allow for the application of the SCMOCC method to much lower collision energies (approaching 1 eV/u). In addition, a computational efficiency study for various propagators is performed to speed-up the calculations on classical computers. These computations are repeated for the simulation of the SES approach assuming typical parameters for realistic pretheshold superconducting quantum computing hardware. The feasibility of applying future SES processors to the quantum dynamics of large molecular collision systems is briefly discussed.

研究动机与目标

• 将化学动力学的量子模拟扩展至低碰撞能量（接近1 eV/u），传统方法因常速度近似而失效。
• 通过引入多通道势能平均化与埃伦费斯特定时对称化，克服半经典分子轨道耦合理论（SCMOCC）方法的局限性。
• 识别模拟SES动力学时最高效的经典数值传播器，以实现未来量子硬件的快速基准测试。
• 评估未来预阈值超导量子处理器在多达十个通道的大型分子碰撞体系中应用的可行性。
• 建立基于真实硬件参数的n-量子比特SES处理器，连接经典模拟与未来量子模拟的计算桥梁。

提出的方法

• 在n维单激发子空间（SES）中应用SCMOCC方法，以模拟n通道双原子碰撞动力学，将量子比特态映射至绝热分子通道。
• 引入多通道势能平均化方案，以处理非恒定相对速度，取代限制性强的直线恒速假设。
• 应用埃伦费斯特定时对称化，校正势能平均化导致的详细平衡破坏，提升低能区域的精度。
• 对Crank-Nicholson、切比雪夫与自适应龙格-库塔等数值传播器进行全面研究，预处理显著提升了效率。
• 利用盖尔绍林圆定理估算哈密顿量的特征值边界，以通过矩阵对角化实现高效的时域演化。
• 在真实预阈值超导量子比特硬件参数下重复经典模拟，以评估未来量子模拟的可行性。

实验结果

研究问题

• RQ1 能否通过多通道势能平均化，将SCMOCC方法扩展至低动能（≤1 eV/u）分子碰撞的模拟？
• RQ2 埃伦费斯特定时对称化在恢复平均势能模型中细节平衡并提升低能碰撞模拟精度方面效果如何？
• RQ3 在经典计算机上模拟SES动力学时，哪种数值传播器在精度与计算速度之间提供了最佳权衡？
• RQ4 在使用优化传播器的情况下，经典SCMOCC方法可高效模拟的最大系统规模（以通道数计）是多少？
• RQ5 具有n(n−1)/2个可调耦合器的预阈值超导量子比特硬件，在SES中模拟大规模分子碰撞动力学的潜力有多大？

主要发现

• 预处理自适应步长费尔伯格-龙格-库塔方法在n=5的Si³⁺+He问题中，相比标准Crank-Nicholson传播器实现了220倍加速。
• 自适应龙格-库塔方法将计算时间缩短至0.88秒，局部误差控制在10⁻⁴以内，显著优于固定步长与对角化方法。
• 埃伦费斯特定时对称化有效校正了平均势能模型中因势能平均化导致的详细平衡损失，实现了低能模拟的高精度。
• 结合多通道势能平均化的SCMOCC方法成功将模拟扩展至1 eV/u的低动能区域，克服了常速度近似的局限性。
• 对于n=5的Si³⁺+He体系，SCMOCC方法计算的电荷交换截面与使用相同绝热势能的先前QMOCC计算结果一致，验证了方法的有效性。
• 本研究证明，基于SES的量子模拟在当前预阈值超导量子比特硬件参数下，对多达十个通道的分子体系是可行的。