利用变分量子本征求解器研究横场伊辛模型的纠缠与拟设表达能力

A Study of Entanglement and Ansatz Expressivity for the Transverse-Field Ising Model using Variational Quantum Eigensolver

arXiv: 2602.17662v1

论文信息

标题: A Study of Entanglement and Ansatz Expressivity for the Transverse-Field Ising Model using Variational Quantum Eigensolver

作者: Ashutosh P. Tripathi, Nilmani Mathur, Vikram Tripathi

发布日期: 2026-02-19

arXiv ID: 2602.17662v1

PDF 链接: 下载 PDF

3 分钟速览

  • 研究问题：这篇论文要解决什么 研究变分量子本征求解器（VQE）在模拟横场伊辛模型（TFIM）时的表现，重点关注不同量子电路结构（ansatz）在强纠缠和基态简并区域的准确性和稳定性。

  • 核心方法：用什么办法解决 采用精确态矢量模拟，对比三种 ansatz（硬件高效型 HEA、哈密顿量启发型 HVA、加入对称破缺的 HVA-SB），在 1D、2D、3D 系统中基准测试其能量方差、纠缠熵和自旋关联。

  • 关键结果：最重要的一个结论或数字 HEA 表达力最强但会低估低场区域的纠缠熵（见论文图 4b），而 HVA 虽然参数少却能在哈密顿量子空间内实现更稳定的能量估计，能量方差显著低于 HEA（见论文图 3）。

  • 主要局限：作者自己承认的、或方法本身固有的限制 浅层电路难以捕捉强关联态导致纠缠熵系统性地被低估，HVA 在低纠缠区域表现较差，3D 模拟中优化难度急剧增加，且有限尺寸效应使得对称破缺态的解释需要向热力学极限外推。

  • 适合读者：什么背景的人值得读 适合从事量子计算、凝聚态物理变分模拟、或 NISQ 算法设计的研究者和工程师，尤其是关注 ansatz 选择和优化策略对模拟精度影响的读者。

论文背景和研究动机

变分量子本征求解器（VQE）被广泛认为是含噪中等规模量子（NISQ）时代最有前景的混合量子-经典算法之一。它的核心思想是将量子系统的试探波函数用参数化量子电路（PQC）表示，通过经典优化器调节参数来逼近基态能量。这个框架在量子化学、凝聚态物理等领域已有不少应用案例。

然而 VQE 的有效性高度依赖两个假设：其一，所选用的 PQC 能够表示真实的基态；其二，经典优化器能够收敛到正确的参数配置。当系统处于强纠缠或基态近简并状态时，这两个假设都可能失效。横场伊辛模型（TFIM）恰好为检验这些问题提供了理想的测试平台——它在横场较小时基态高度纠缠且近简并，这正是 VQE 容易出错的场景。

已有的 TFIM 与 VQE 相关研究未能系统性地将 ansatz 表达力、纠缠特性和优化器行为放在同一框架下比较。这项工作的动机正是填补这一空白：通过基准测试不同 ansatz 在 TFIM 各维度和参数区间的表现，揭示电路结构与模拟精度之间的本质联系，为后续研究提供选型依据。

核心方法和技术细节

系统模型与 Hamiltonian

论文研究带周期性边界条件的铁磁横场伊辛模型，Hamiltonian 为：

H=Jz∑⟨ij⟩σizσjz+hx∑iσixH = J_z \sum_{\langle ij \rangle} \sigma_i^z \sigma_j^z + h_x \sum_i \sigma_i^x

其中 Jz=−1.0J_z = -1.0，hxh_x 为可调的横场强度。系统在一维达到 15 个量子比特，二维为 4×44 \times 4 格点，三维采用最大 27 个量子比特的晶格。

该模型基态具有清晰的物理图像：当 hx≪hch_x \ll h_c（hch_c 为量子相变临界场）时，基态近似为 ∣00…00⟩+∣11…11⟩\vert 00\ldots 00 \rangle + \vert 11\ldots 11 \rangle 的叠加，呈现高纠缠的 Schrödinger 猫态特征，且与第一激发态近简并；当 hx≫hch_x \gg h_c 时趋向乘积态 ∣++…++⟩\vert ++\ldots ++ \rangle，纠缠较低。

三种 ansatz 电路

论文比较了三类参数化量子电路：

  1. 硬件高效 ansatz（HEA，即 Qiskit 的 EfficientSU2）：由单比特旋转门和 CNOT 纠缠层交替堆叠而成，参数数量随系统规模线性增长，表达力最强（见论文图 2）。

  2. 哈密顿变分 ansatz（HVA）：通过 Trotter 化 TFIM Hamiltonian 直接构造电路——RZZR_{ZZ} 门实现最近邻耦合项，RxR_x 门来自横场项。参数量少但搜索空间被严格限制在物理相关的 Hilbert 子空间内。

  3. 对称破缺 HVA（HVA-SB）：在 HVA 基础上附加 RzR_z 层来打破 Z2Z_2 宇称对称性，赋予电路访问近简并激发态的能力，在有限尺寸下获得更接近热力学极限的对称破缺态描述。

表达力与观察量定义

表达力通过帧势（frame potential）量化：Ft=Eψ,ϕ[∣⟨ψ∣ϕ⟩∣2t]F_t = \mathbb{E}{\psi,\phi}[\vert \langle \psi \vert \phi \rangle \vert^{2t}]，其中态从随机采样的电路参数生成。帧势越小表示产生的态在 Hilbert 空间中越均匀，表达力越强。图 2 的直方图显示 HEA 的重叠分布最窄，表达力最高。

论文定义了以下基准指标：

  • 每自旋磁化强度：M=1N∑i⟨σiz⟩M = \frac{1}{N}\sum_i \langle \sigma_i^z \rangle，由于基态对称性其精确值为零，非零值反映基态与激发态的混合程度。
  • 自旋关联函数：MCorr=1N/2∑i⟨σizσi+N/2z⟩M_{\mathrm{Corr}} = \frac{1}{N/2}\sum_i \langle \sigma_i^z \sigma_{i+N/2}^z \rangle，对近简并有限系统更为可靠。
  • 能量方差：Var(E)=(⟨H2⟩−⟨H⟩2)/⟨H⟩2\mathrm{Var}(E) = (\langle H^2 \rangle - \langle H \rangle^2)/\langle H \rangle^2，衡量态偏离精确本征态的程度，已有工作证明变分能量误差与方差呈线性关系。
  • 冯·诺依曼纠缠熵：SEE=−∑iλilog⁡λiS_{EE} = -\sum_i \lambda_i \log \lambda_i，通过单比特约化密度矩阵或 Schmidt 分解计算。

模拟环境与优化策略

所有模拟在 NVIDIA CUDA-Q 框架下进行精确态矢量仿真，使用 L-BFGS 优化器处理 HEA（因其优化景观较平滑），用无梯度 COBYLA 算法处理 HVA 和 HVA-SB（因其优化地形崎岖）。电路深度测试了 4、8、10、15 层，部分 3D 模拟采用去掉 RzR_z 的实振幅变体以降低参数数量。

创新点和贡献

这项工作的贡献不在于提出新算法，而在于提供了目前较为系统的 ansatz 对比基准研究，具有以下创新性：

第一，多维度的表达力-精度权衡分析。 论文同时考察了能量精度（能量方差）、纠缠表征能力（冯·诺依曼熵）、以及相变信号捕捉能力（自旋关联函数），揭示了 HEA 与 HVA 在不同物理区间的互补性：HEA 表达力强但在低场高纠缠区倾向于选择对称破缺态导致纠缠被低估，HVA 在哈密顿子空间内稳定但低纠缠区域表现欠佳（见论文第 3 节，图 4、图 5）。

第二，HVA 的能量方差跳变现象。 论文观察到随电路层数增加，HVA 和 HVA-SB 的能量方差并非逐步改善，而是在某个临界层数出现从 “差” 到 “准” 的突然转变（见论文图 3 及相关讨论）。这一现象揭示了限制在物理子空间内的 ansatz 具有门槛效应，对电路深度选择具有指导意义。

第三，对称破缺层的实际价值验证。 将 HVA-SB 引入 TFIM 并非首创，但论文通过对比对称破缺与否对低场态磁化和纠缠的影响，定量说明了该设计如何在有限系统中模拟热力学极限的自发对称破缺行为，同时不丢失临界点附近的定性特征。

实验结果分析

一维系统（10 自旋）

能量方差随 hxh_x 的变化（图 3）清晰展示了三类 ansatz 的差异：HVA 在整个横场范围内保持最低的方差，表明其虽参数少但搜索效率高；HEA 方差波动较大，优化不稳定；HVA-SB 介于两者之间。值得注意的是，增加层数时 HEA 收敛缓慢，而 HVA 呈现阈值效应。

自旋关联（图 4a）在临界点附近出现明显转折，三种 ansatz 均能定性反映相变。纠缠熵（图 4b）则暴露了关键问题：HVA 在低场区给出每自旋 ln⁡2\ln 2 的最大纠缠熵（对应猫态），而 HEA 和 HVA-SB 给出接近零的熵值（对应对称破缺的乘积态）。这解释了为什么 HEA 能获得合理能量却低估纠缠——它找到的是有限系统的近似对称破缺态而非真正的基态叠加。

二维系统（4×44\times4 格点）

二维情况下优化难度显著上升。HVA 在低纠缠（大 hxh_x）区域表现不佳，HVA-SB 的能量方差甚至可超过 HEA（图 5 相关分析）。自旋关联和纠缠熵仍能识别临界点（约 hx≈3.0h_x \approx 3.0），但不同初始参数会落入不同局部极小值，增加了定量比较的复杂性。

三维系统与规模扩展

在三维格点中，论文采用简化的实振幅 ansatz 才获得可接受的优化效果。基态平均能量几乎不随格点尺寸变化（图 6a），纠缠熵仍保留临界行为特征（图 6b），但绝对数值低于精确值，进一步证实了浅层电路在纠缠表征上的局限性。

实践建议

对于正在或计划使用 VQE 进行凝聚态模拟的实践者，这项工作提供了几条可操作的指导：

优先根据物理特性选择 ansatz。 HEA 的通用性虽强，但在有限系统中容易落入对称破缺的局部极值，适合作为能量参考值生成器，而非精确纠缠测量的工具。若目标是测量纠缠熵这类对基态质量敏感的观察量，HVA 或其对称破缺变体更为可靠——至少它们的偏差方向是可预期的。

合理分配电路深度资源。 HVA 类的层数并非 “越多越好”——实验表明存在性能突变的层数阈值，超过该阈值后继续增加深度收益递减。在模拟前通过小规模基准找到该阈值（例如监控能量方差的骤降点）可以节省大量计算资源。

优化器与 ansatz 需要配对考虑。 论文的经验是 HVA 的崎岖优化景观需要 COBYLA 等无梯度方法，而平滑的 HEA 可从 L-BFGS 中获益。如果混合使用 ansatz，应为每段电路单独配置优化策略。

参数迁移加速收敛。 3D 模拟中采用的 “参数重用” 策略——先优化某个 hxh_x 值，再将其参数作为邻近值的初始点——能显著减少迭代次数，这一技巧在需要扫描相图时尤其有效。