变分量子本征求解器的实用方法与挑战是什么?

梳理VQE的编码、拟设、测量与噪声缓解流程,并以O2解离曲线对照CASCI基准。

直接答案

VQE用混合循环取代了通往基态能量的深度相位估计路线:在近中期量子硬件上制备参数化试探态,再由经典优化器最小化其能量。这篇论文梳理了从活性空间选择、费米子到量子比特映射,到拟设选择、测量分组和噪声感知误差缓解的实用工作流,并在O2解离曲线上表明,USCC-SDTQ拟设相对CASCI达到化学精度,而在VQE态之上施加NEVPT2修正后,与CASCI+NEVPT2的差异远低于该阈值。

1篇文献引用

本文由 WisPaper 驱动的搜索和论文分析生成。

VQE为何用经典循环换取电路深度

论文把VQE定位为对具体资源问题的回应。量子相位估计结合绝热演化是电子结构最早被提出的量子计算应用之一,但它需要大量受控哈密顿量操作,产生深电路,其电路深度、纠缠门数量和达到化学精度所需的量子比特数都远超当前含噪中等规模量子设备的能力。

VQE通过以下方式绕开这一限制:在量子设备上制备参数化试探态,并用经典优化程序最小化其能量。长时相干演化被更短的参数化电路取代,参数由经典方式更新,期望值测量取代相位估计。论文明确指出,这使算法在概念上简单,但在实践中各环节相互关联:编码与拟设选择决定电路结构、精度和资源开销,而优化、测量和硬件噪声又彼此制约。

活性空间与嵌入方法决定精度上限

在构建任何电路之前,论文把活性或关联子空间的选择视为决定性的精度杠杆。子空间大小决定精度与成本的权衡;化学惰性轨道被冻结,而所有强关联轨道都被纳入,直到目标可观测量收敛。论文指出,自然轨道占据数以及纠缠或互信息等高级选择标准可以提高精度。

随后,非活性环境在平均场层面通过量子嵌入方法被重新纳入,该方法为活性轨道构造有效哈密顿量,其中包含由Hartree-Fock或密度泛函理论计算得到的非活性轨道之间及其内部的相互作用。论文列举了基于电子密度、基于密度矩阵和基于格林函数的三类方法,包括密度泛函嵌入理论、密度矩阵嵌入理论、动力学平均场理论、自能嵌入理论和量子缺陷嵌入理论。剩余的活性空间模型显著更小,因而可处理。

费米子到量子比特映射决定泡利权重

论文描述了从费米子到量子比特的两步转换:先把费米子Fock态与量子比特态对应,再把产生和湮灭算符改写为泡利算符的张量积,使电子哈密顿量成为泡利字符串的线性组合。比较不同映射时最关键的量为泡利权重,即字符串中非恒等算符的数量,以及目标期望值中非对易泡利字符串的总数。

这些量相互关联,因为权重更低的字符串更容易对易,而只有非对易泡利字符串才需要单独执行电路,因此更大的对易组可以节省计算时间。Jordan-Wigner映射把轨道占据数局域地存储在各自量子比特上,但把宇称非局域地存储在所有前序量子比特上,因此作用于最后一个基函数的算符会产生作用于所有量子比特的泡利算符,即最大可能的泡利权重。宇称映射改为局域存储宇称而非占据数,泡利权重相近,但优势在于粒子数通常已知,且系统可以按自旋扇区拆分。论文还提醒,如果非恒等操作位于不相连的量子比特上,有限的硬件连通性会迫使引入swap操作,从而大幅增大电路。

拟设选择在表达力与对称性间权衡

论文把拟设构建呈现为表达力、对称性保持和电路深度之间的平衡。它讨论的一个具体起点是二次的、保持粒子数的哈密顿量,其解为费米高斯态,在该特定情形下是Slater行列式的线性组合。近期提出的高效制备算法最大深度随量子比特数N按O(N^2)增长,构造由Givens旋转组成的电路,这类参数化双量子比特门保持粒子数,可由一个Rxx与一个Ryy门的组合,或三个CX门配合单量子比特旋转实现。

论文指出这些解的使用方式存在权衡:既可以在量子计算机上制备Slater行列式,也可以把哈密顿量变换到单粒子解态,对Hartree-Fock而言即从原子基变换到分子轨道。在该基下基态是乘积态,由X门初始化,电路长度最优,但原子基下的哈密顿量可能比分子基下更简单,从而形成更短更简单电路与更少所需测量次数之间的权衡。在Hartree-Fock之外,CISD、CCSD、CASCI和选择性组态相互作用等方法给出的基态可表示为Slater行列式的叠加,可作为初始态。

测量成本与贫瘠高原主导整个流程

论文指出期望值估计往往主导整体计算成本,因为真实哈密顿量会分解为大量必须反复测量的泡利项。这正是所选映射的对易组结构在实践中重要的原因,而不只是理论上的意义。

经典优化必须在采样噪声和崎岖的优化景观中进行,论文补充说,它还可能进一步受到贫瘠高原的阻碍,此时梯度变得过小,无法高效训练。硬件噪声、有限相干时间和编译约束则要求误差缓解和硬件感知的电路设计。论文的总体主张是,成功的VQE实现需要在工作流的所有阶段仔细平衡各种权衡,而不是只优化单一组件。

泡利twirling把未知噪声变为可测信道

论文把噪声描述为既来自环境耦合,也来自不完美的控制操作,包括串扰和过旋转或欠旋转,在许多平台上,幺正误差还会被系统漂移进一步放大,需要定期重新校准控制脉冲。由于误差来源和类型繁多且难以表征,论文把通过随机编译实现的泡利twirling作为标准工具,将一般噪声转换为已知噪声信道。

对噪声过程在泡利群上做twirling,得到的信道其泡利转移矩阵只保留对角元,因此所有相干和非幺正贡献都被去除,剩下一个泡利信道。实践中,随机编译把Clifford门U替换为一族逻辑等价、由单量子比特泡利门修饰的电路,论文指出实验上只需选取数量合理且随机选择的一小组电路即可获得良好的twirling结果。随后用循环基准测试测量得到的泡利信道,论文称该方法对态制备和测量误差不敏感。论文也指出一个实际限制:测量n量子比特系统泡利转移矩阵的每一个元素并不可行,因此必须假设所关注的信道服从某种稀疏、可高效学习的模型,例如相邻量子比特上的2-局域泡利误差。

读出缓解修正计数但带来采样开销

针对测量误差,论文描述了一个转移矩阵,其元素给出含噪读出信道把结果a变为结果b的概率,可通过反复制备计算基态并记录结果获得。如果该矩阵已知且可逆,缓解后的概率可由其逆矩阵作用得到,但结果可能因采样噪声而出现负值,因此实践中常改为最小化平方残差。由于完整矩阵维度为2^n,直接处理只对少量量子比特可行;更大规模实例需要假设量子比特间噪声不相关,使矩阵成为单量子比特转移矩阵的张量积,或假设其稀疏。

论文还介绍了读出twirling,它通过对误差模型构造概率混合来给出期望值的无偏估计,即在每个量子比特上随机地在测量前施加X门并翻转对应比特。论文指出,一旦计算出全零态的twirled期望值,采样开销即由该量的平方倒数给出,使修正的代价明确而非免费。

O2解离曲线显示VQE达到化学精度

论文的定量基准是O2解离曲线。收敛性按酉选择耦合簇拟设的激发等级报告:从USCC-SD到USCC-SDTQ改进了波函数描述,使VQE结果更接近CASCI参考,误差图显示USCC-SDTQ在化学精度阈值内收敛,论文将该阈值定义为1.6×10^-3 Ha。

在VQE参考态之上施加的NEVPT2修正有效恢复了动态相关,并最小化最终解离能对具体拟设截断的依赖。在1.0至4.0埃的解离曲线上,VQE+NEVPT2与CASCI+NEVPT2之间的能量差远低于化学精度阈值,论文将此解读为NEVPT2无论建立在量子还是经典参考波函数上都表现同样良好的证据。相对Active Thermochemical Tables给出的实验解离能0.188035636 Ha,VQE+NEVPT2和CASCI+NEVPT2的误差都远在化学精度之内,而单参考CCSD(T)的误差显著更大,论文将其归因于多参考方法对具有强静态相关系统的重要性。

本文引用的文献

1

A Practical Introduction to VQE: Methods, Challenges, and Applications in Physics and Chemistry

论文梳理了VQE工作流,从活性空间选择和费米子到量子比特映射,到拟设构建、测量分组和噪声感知误差缓解,并报告在O2解离曲线上,USCC-SDTQ拟设相对CASCI在1.6×10^-3 Ha化学精度阈值内收敛,而VQE+NEVPT2与CASCI+NEVPT2的差异远低于该阈值,且两者相对实验解离能0.188035636 Ha都保持在化学精度之内。