1. VQE变分本征求解器概述

变分量子本征求解器(Variational Quantum Eigensolver,VQE)是当前含噪声中等规模量子(NISQ)时代最具代表性的量子-经典混合算法之一。它要解决的问题非常明确:给定一个物理系统的哈密顿量 H,求出其最低本征值 E0 以及对应的本征态 |ψ0⟩。在量子化学中,这个最低本征值就是分子的基态能量;在材料模拟中,它对应体系基态、能隙与相稳定性的判断依据。VQE 的核心思想并不复杂:用一个带参数的量子线路制备试探态 |ψ(θ)⟩,在量子处理器上测量该态下哈密顿量的期望值 ⟨ψ(θ)|H|ψ(θ)⟩,再把这个期望值交给经典优化器调整参数 θ,循环往复,直到期望值收敛到最小值。

变分原理:VQE 的理论基石

VQE 之所以成立,依赖量子力学中的变分原理。对于任意归一化试探波函数 |ψ(θ)⟩,其能量期望值满足:

E(θ) = ⟨ψ(θ)|H|ψ(θ)⟩ ≥ E₀

其中 E0 是 H 的真实基态能量。等号仅在 |ψ(θ)⟩ 恰好等于基态时成立。这条不等式给出了一个极其重要的工程含义:期望值永远不可能低于真实基态能量,因此我们可以放心地把“最小化期望值”作为优化目标,而不会得到物理上错误的下界。同时,变分原理也揭示了 VQE 的精度上限——无论优化器多强,结果精度最终受限于试探态(即拟设,ansatz)所能张成的希尔伯特子空间是否包含真实基态。换句话说,拟设决定可达精度,优化器决定能否逼近该精度。

需要特别说明符号约定:本文中哈密顿量 H 的本征值按从小到大排列,E0 为基态能量,能量单位统一采用 Hartree(哈特里)或 eV,具体在示例中标注。测量得到的期望值始终是实数,优化目标即最小化该实数。

量子-经典混合架构

VQE 不是纯量子算法,也不是纯经典算法,而是一个闭环反馈系统。其运行结构可以概括为:经典侧准备哈密顿量的泡利分解与初始参数;量子侧按参数制备试探态并逐项测量泡利字符串的期望值;经典侧把加权求和得到的能量送入优化器,生成新参数;如此迭代。下图展示了这一混合架构中各环节的职责与数据流向。

VQE 量子-经典混合闭环架构 经典处理器 量子处理器 哈密顿量泡利分解 H = Σ cᵢ Pᵢ(如 H₂ 共 15 项) 经典优化器 COBYLA / SLSQP / 梯度下降 能量加权求和 E(θ) = Σ cᵢ ⟨Pᵢ⟩ 参数化拟设线路 U(θ) 制备 |ψ(θ)⟩ 逐项测量泡利字符串 多次采样估计 ⟨Pᵢ⟩ 测量结果统计 比特串计数 → 期望值 参数 θ ⟨Pᵢ⟩ 测量值 E(θ) 收敛?|ΔE| < 阈值

理解这张图的关键在于两点。第一,量子处理器只负责“制备态并测量”,它不承担优化计算;经典处理器只负责“分解、求和、更新参数”,它不直接接触量子态。第二,整个系统是迭代闭环,每一轮的能量评估都需要重新制备量子态并重复测量,因此测量采样次数(shots)直接决定能量估计的统计噪声,而噪声又会影响优化器对梯度的判断。这是 VQE 工程实现中最容易被低估的细节。

本征值问题的具体形态

在进入 VQE 之前,物理问题已经被转化为量子比特形式。以氢分子 H2 在 STO-3G 基组、键长 0.735 Å 附近为例,经过 Jordan-Wigner 或 Bravyi-Kitaev 变换后,电子哈密顿量可写成若干泡利字符串的线性组合,项数约为 15 项(4 个量子比特)。每一项形如 ci·Pi,其中 Pi ∈ {I, X, Y, Z}⊗n。VQE 的任务就是找到参数 θ,使 Σci⟨Pi⟩ 最小。这个例子说明:本征值问题在 VQE 语境下不是抽象矩阵对角化,而是“泡利字符串期望值的加权和最小化”。

VQE 在量子化学与材料模拟中的定位

VQE 的定位需要放在与经典方法的对比中理解。经典量子化学方法如 CCSD(T) 在弱关联体系上精度极高且成本可控,但在强关联体系(如过渡金属配合物、某些磁性材料、键断裂过程)中会因多参考特性而失效或成本爆炸。VQE 的优势在于拟设可以灵活设计以捕捉强关联,且线路深度相对浅,适配当前含噪声设备。其代价是:需要大量测量采样、优化可能陷入局部极小、拟设表达能力与线路深度之间存在权衡。因此 VQE 更适合作为“经典方法难以处理的中等规模强关联问题”的补充工具,而不是全面替代经典方法。

课程项目总览

本课程以一个可完整跑通的 H2 基态能量求解项目为主线,贯穿 30 个章节。项目最终交付物包括:一份哈密顿量泡利分解结果、一个参数化拟设线路、一套经典优化循环、一组不同键长下的能量曲线,以及与精确对角化结果的误差对比。建议在项目初期就固定以下约定,避免后续返工:量子比特数 n=4(H2,STO-3G);拟设采用硬件高效型或 UCCSD 简化版;优化器先用 COBYLA(无梯度、抗噪),后期再引入参数移位梯度;每轮能量评估 shots 设为 1024 或 4096,并记录统计误差。项目验收的判断依据是:在平衡键长附近,VQE 能量与精确基态能量之差(化学精度参考值约 1.6 mHartree)是否达到可接受范围,以及能量曲线是否呈现正确的解离行为。