3. 变分拟设设计基础
变分拟设(ansatz)是 VQE 中唯一由人设计的部分:它决定了参数化量子电路能表示哪些量子态,也决定了优化器要在多大的参数空间里搜索。同一个分子哈密顿量,换一个拟设,基态能量的精度、参数数量、电路深度和优化难度都会显著变化。因此,拟设设计不是“随便搭几层旋转门”,而是在表达能力、硬件可行性和参数规模之间做工程取舍。
一、硬件高效拟设:先能跑,再谈准
硬件高效拟设(Hardware-Efficient Ansatz, HEA)的核心思想是:不追求物理意义,而是直接使用目标量子设备最擅长执行的门——通常是单比特旋转门和相邻比特间的 CNOT/CZ。典型结构是“旋转层 + 纠缠层”交替堆叠:每一层对每个量子比特施加 RX(θ)、RY(θ) 或 RZ(θ),随后按硬件拓扑连接 CNOT 形成纠缠。重复 L 层,参数总数约为 L × n × k(n 为比特数,k 为每比特旋转门数)。
HEA 的优点是深度可控、天然适配硬件连通性,适合近期含噪设备。缺点是它没有编码任何化学结构信息,容易陷入“表达能力强但物理无关”的困境:参数多、优化地形复杂,且可能收敛到非物理的局部极小。实践中,HEA 常用于小分子快速验证或作为基准对照,而不是高精度基态计算的首选。
二、UCCSD 拟设:把化学知识写进电路
UCCSD(Unitary Coupled Cluster with Singles and Doubles)来自量子化学的耦合簇方法。它从 Hartree-Fock 参考态 |HF⟩ 出发,用幺正算符 e^(T−T†) 构造试探态,其中 T = T1 + T2 分别对应单激发和双激发。每个激发项对应一个参数,物理意义明确:单激发描述电子从占据轨道跃迁到虚拟轨道,双激发描述两个电子的关联跃迁。
在电路实现上,每个激发项通常用 Jordan-Wigner 或 Bravyi-Kitaev 映射转成若干 CNOT 和旋转门。以 4 比特的 H₂ 为例,若采用冻结核心近似,仅保留 2 个活性电子和 2 个活性轨道,UCCSD 只需 1 个双激发参数,电路深度远小于同等精度的 HEA。UCCSD 的优点是精度高、参数少、可解释;缺点是电路深度随系统规模增长较快,且对硬件连通性要求高,需要大量 SWAP 或长程 CNOT。
示例:H₂ 分子在 STO-3G 基组下的拟设对比(合理假设,非实测)
- 量子比特数:4(冻结 1s 核心后活性空间 2 轨道 → 4 自旋轨道)
- HEA:2 层,每层 4 个 RY + 4 个 CNOT,共 8 个参数,深度约 10
- UCCSD:仅 1 个双激发参数,电路约 6 个 CNOT,深度约 8
- 在理想模拟下,UCCSD 单参数即可达到化学精度(约 1.6 mHa);HEA 需优化 8 个参数才能接近,且可能落入局部极小
三、表达能力与参数数量的权衡
拟设的表达能力(expressibility)指它能在多大范围内覆盖希尔伯特空间中的量子态。表达能力越强,越有可能包含真实基态;但表达能力过强也会带来两个问题:一是参数数量增加,优化维度升高,训练成本变大;二是可能出现过参数化(over-parameterization),导致梯度消失或贫瘠高原(barren plateau),使优化器难以找到有效方向。
参数数量并非越多越好。对于固定比特数 n,UCCSD 的参数数由激发项数决定,通常远小于 HEA 的 L×n×k。但 UCCSD 的每一项都对应物理激发,因此“参数效率”更高。工程上常用以下判断依据:
- 化学精度目标:若只需 1.6 mHa 精度,小分子优先 UCCSD 或 UCCSD 的截断版本。
- 硬件深度限制:若设备相干时间短、连通性差,HEA 更易映射,但需接受精度损失。
- 参数规模:参数数超过约 10^2 时,经典优化器(如 BFGS、COBYLA)容易陷入局部极小,需考虑自适应拟设(ADAPT-VQE)或分层优化。
- 对称性:利用粒子数守恒、自旋对称性可裁剪拟设,减少无效参数。
四、从设计到验证:一个可操作的检查流程
实际设计拟设时,建议按以下顺序推进:先确定活性空间和映射方式,再选择拟设类型,然后估算参数数量和电路深度,最后在小规模上做数值验证。验证时不要只看最终能量,还要观察参数梯度范数和优化收敛曲线:若梯度在训练初期就接近零,说明拟设可能过于表达或初始化不当;若能量随迭代震荡不降,可能是纠缠层不足或学习率不合适。
# 伪代码:拟设参数规模快速估算
n_qubits = 4 # 活性自旋轨道数
n_electrons = 2 # 活性电子数
n_occ = n_electrons
n_vir = n_qubits - n_electrons
# UCCSD 激发项数(近似)
n_singles = n_occ * n_vir
n_doubles = n_occ * (n_occ - 1) / 2 * n_vir * (n_vir - 1) / 2
n_params_uccsd = n_singles + n_doubles
# HEA 参数数
L = 2 # 层数
k = 2 # 每比特旋转门数
n_params_hea = L * n_qubits * k
print(n_params_uccsd, n_params_hea) # 示例输出:1 16
上面的示例中,UCCSD 仅需 1 个参数,而 HEA 需要 16 个。这并不意味着 HEA 一定更差,而是说明:在化学精度要求明确、活性空间较小时,UCCSD 的参数效率明显更高;当系统变大、UCCSD 电路深度不可接受时,才考虑用 HEA 或自适应拟设做折中。
总结来看,拟设设计没有万能答案。HEA 胜在硬件友好和深度可控,UCCSD 胜在物理可解释和参数高效。实际工程中,应先明确目标精度和硬件约束,再选择拟设类型,并通过参数数量估算和梯度行为观察来验证设计是否合理。对于小分子高精度场景,UCCSD 及其截断版本通常是更稳妥的起点;对于近期含噪设备上的快速探索,HEA 可作为基准;当系统规模继续增大时,自适应拟设和对称性裁剪会成为必要的补充手段。