Ground state energy and phase transitions of Long-range XXZ using VQE

本文提出了一种带有约束拟设电路的变分量子本征求解器（VQE）方法，通过分析不同相中基态能量误差的敏感性来探测长程 XXZ 链中的无穷阶相变，同时通过与精确对角化结果进行对比来验证该方法。

要点

想象一下，你正试图在一片广袤、大雾弥漫的景观中寻找最低点。在物理学世界中，这个“最低点”被称为基态能量（ground state energy），它告诉我们一个系统（比如磁铁中的原子）在完全平静时想要趋向于何种状态。

通常情况下，为复杂的系统寻找这个最低点，就像是在解一个规模大到人类大脑甚至世界上最强大的超级计算机都无法处理的谜题。这正是该论文作者引入一种新工具的地方：变分量子特征值求解器（Variational Quantum Eigensolver, VQE）。你可以把 VQE 想象成一个聪明的混合动力机器人，它利用量子计算机做出一个关于最低点的“猜测”，再利用经典计算机不断优化这个猜测，直到它尽可能接近真相。

挑战：两种类型的“边界”

研究人员研究了一个特定的模型，叫做长程 XXZ 链（Long-Range XXZ chain）。想象一排微小的磁铁（自旋），它们可以彼此交流。通常情况下，磁铁只与相邻的邻居交流，但在这种模型中，它们可以向整条线进行“大声呼喊”（长程相互作用）。

团队想要寻找这些磁铁行为发生彻底改变的“边界”。这些被称为相变（phase transitions）。他们发现了两种类型的边界：

1. “悬崖”（一级相变/First-Order Transition）： 这就像从陡峭的悬崖边走下。能量会突然且剧烈地变化。它很容易被发现，因为地面直接掉落了。
2. “斜坡”（无穷阶相变/Infinite-Order Transition）： 这要棘手得多。它就像是在走一段非常平缓、平滑的小山坡。没有突然的坠落或悬崖；变化发生得如此缓慢，以至于标准工具根本看不见这个边界。通常，科学家需要一个特殊的“全局地图”（复杂的序参数）来找到它，而这很难计算。

秘密武器：“糟糕的猜测”策略

这是论文中最巧妙的部分。通常，科学家使用 VQE 只是为了得到最低能量的精确数值。但作者意识到了一件有趣的事：猜测的质量取决于你所处的位置。

他们设计了一个特定的规则来构建他们的量子机器人（即“拟设电路/ansatz circuit”）：它必须保持总自旋（磁化强度）恒定。

• 在“正确的”邻里间： 如果机器人在一个磁铁自然倾向于保持该特定恒定自旋状态的相中，机器人就会做出一个很棒的猜测。误差（机器人的猜测与真实答案之间的差异）非常微小。
• 在“错误的”邻里间： 如果机器人在一个磁铁并不想处于那种状态的相中，机器人就会陷入挣扎。它试图强迫磁铁呈现错误的形状，导致误差变得巨大。

“指南针”类比

为了寻找那些隐形的边界，作者不仅仅观察误差的大小，他们还观察了误差的方向。

想象你在森林中行走，每走一步都会丢下一个指南针。

• 在森林的一个部分（相 A），指南针的指针指向随机的方向，疯狂旋转。
• 在另一个部分（相 B），指南针的指针全部整齐地指向同一个方向。

作者使用了一种称为**方向相干性（Directional Coherence）**的技术来测量这一点。他们计算了数千个点的“误差”，并观察变化的趋势。

• 当指南针的指针表现得混乱时，他们知道自己处于一个相中。
• 当指针突然变得整齐划一时，他们知道自己跨越了边界。

通过观察机器人的“踉跄”方式，他们成功识别出了既容易发现的“悬崖”边界，也包括隐藏的“斜坡”边界。他们不需要新的、复杂的地图；他们只需要观察机器人的错误表现。

研究结果

• 对于容易发现的边界（一级相变）： 他们看到了误差像悬崖一样突然跳升。
• 对于困难的边界（无穷阶相变）： 他们看到了“指南针指针”（误差梯度）从混沌变得有序。这揭示了标准方法所遗漏的边界。
• 准确度： 通过让机器人的“大脑”变得稍微“深”一点（增加电路层数），他们甚至能以极高的精度计算系统的精确能量（误差在 3-7% 以内），即使是在困难的情况下也是如此。

核心结论

这篇论文声称，你并不总是需要完美的计算才能找到系统的变化点。有时，研究你的计算如何失败（即误差分析）实际上比计算本身更能告诉你关于系统结构的信息。他们成功地利用这种“误差分析”方法，绘制出了长程磁性链中简单及复杂相变的图谱，证明了量子算法不仅可以用来解决问题，还可以用来发现物质行为的隐藏规则。

技术摘要

技术摘要：利用变分量子特征值求解器研究长程 XXZ 链的基态能量与相变

问题陈述
本文探讨了识别量子多体系统（特别是长程 XXZ 链）中相变的挑战。虽然变分量子特征值求解器（VQE）在寻找缺乏解析解的哈密顿量基态能量方面已得到广泛应用，但其在检测相边界（尤其是无穷阶相变）方面的效用仍有待深入探索。传统的相变检测方法依赖于基态能量或二点相关函数的导数奇异性。这些技术对于有限阶相变非常有效，但对于通常需要全局范围序参量或复杂方法（如量子蒙特卡洛）才能处理的无穷阶相变往往失效。作者旨在证明，传统上仅用于能量最小化的 VQE，可以通过分析能量误差相对于系统相位的行为，被改编用于识别一阶和无穷阶相变边界。

方法论
本研究对由以下哈密顿量定义的一维长程 XXZ 链采用了 VQE 算法：
$$H = -J \sum_{i \neq j} \frac{S^x_i S^x_j + S^y_i S^y_j + \Delta S^z_i S^z_j}{|i - j|^\alpha}$$
其中 $J$ 是耦合常数，$\Delta$ 是各向异性，$\alpha$ 定义了相互作用范围。

1. Ansatz 电路设计： 作者设计了一种受哈密顿量变分 Ansatz（HVA）启发但经过效率优化的特定变分 Ansatz。

   • 初始化： 系统从 Néel 态（$|0101\dots\rangle$）开始，通过对 $|000\dots\rangle$ 态施加交替的 Pauli-X 门来准备。
   • 约束： 电路构造旨在在整个优化过程中保持恒定的净自旋（z 方向磁化强度为零）。
   • 结构： 电路将哈密顿量分解为偶键项和奇键项。与标准 HVA 中交换参数的对易项不同，作者对每一个单独的门进行了参数化。这使得电路可以在保持总计 $4(N-1)$ 个参数的同时，实现较低的电路深度（$p=2$）。
   • 门修改： 使用标准 $R_Z$ 门的初始实现被发现效果不佳；这些门被替换为受控-$R_Z$（CRZ）门，以提高基态收敛性。

2. 相变检测策略：

   • 核心假设是，VQE 基态能量的准确性取决于系统相对于初始 ansatz 态所处的相位。
   • 作者计算了 VQE 结果与精确对角化（ED）结果之间的能量差（误差），即 $E_d = E_{exact} - E_{VQE}$。
   • 方向相干性： 为了绘制相边界，作者分析了跨参数空间（$\Delta, \alpha$）的能量差（$E_d$）梯度向量。他们计算了归一化梯度分量，并定义了“方向相干性”，即在移动窗口内梯度角度的平均正弦和余弦值的模。
   • 逻辑： 在 ansatz 适用性良好的相位中（例如给定 Néel 初始化下的反铁磁 AFM 相），误差较低且梯度向量表现出系统性的模式。在 ansatz 失效的相位中（例如铁磁 FM 或无序顺磁 PM 相），误差较高且梯度向量变得随机。这种行为之间的边界即标志着相变。

主要结果
研究针对 $N=12$ 个位点且具有开边界条件的系统进行，比较了 $J=1$（深度 1 和 2）和 $J=-1$（深度 2）的情况。

1. 一阶相变（PM 到 FM）：

   • 对于 $\Delta = 1$ 时从顺磁（PM）到铁磁（FM）的转变，该方法成功识别了边界。
   • 由于初始 Néel 态具有零净磁化强度，VQE ansatz 无法访问 FM 基态。因此，当系统进入 FM 相时，$E_d$ 会剧烈增加。
   • 方向相干性图显示，在 $\Delta = 1$ 处梯度向量的方向发生了剧烈变化，清晰地勾勒出了相变，这与平均场理论的预测一致。

2. 无穷阶相变（PM 到 AFM）：

   • 从 PM 相到反铁磁（AFM）相的转变是一个无穷阶相变，仅靠基态能量导数很难检测。
   • 作者观察到，虽然 $E_d$ 在此转变过程中逐渐增加，但梯度向量的“方向相干性”提供了清晰的区别。
   • 在 AFM 区域，梯度向量排列随机；在 PM 区域，梯度向量具有相干性。一个 0.7 的方向相干性阈值成功分离了这两个区域。
   • 识别出的临界点在 $\Delta = -3$ 时为 $\alpha_c \approx 1.75$，这与通过几何纠缠得到的结果（$\alpha \approx 1.78$）高度吻合。
   • 至关重要的是，作者指出，这种转变在标准的基态能量梯度或能量间隙图中是无法识别的，这验证了其基于误差方法的必要性。

3. 基态能量精度：

   • 使用深度为 2 的电路，作者实现了高精度的基态能量计算。
   • 对于 $J < 0$，最大相对误差约为 3%。
   • 对于 $J > 0$，最大相对误差为 7%。作者将其归因于在没有额外修改（如应用 $R_y$ 门以扩展希尔伯特空间）的情况下，优化器难以到达基态。尽管如此，该方法仍然有效。

意义与主张
本文声称是首个专门利用 VQE 来识别无穷阶相变边界的研究。其主要贡献在于提出了一种技术，该技术利用 ansatz 在某些相位中的“失败”来绘制相边界，而不是仅仅依赖于能量最小化的“成功”。

• 新颖性： 该方法证明，通过从误差梯度和方向相干性的角度分析基态能量数据，可以探测到传统能量诊断手段难以触及的无穷阶相变。
• 效率： 所提出的 ansatz 能以较低的电路深度（$p=2$）和可控的参数数量实现精确结果，适用于近期的量子设备。
• 普适性： 作者认为，这种方法并不局限于 LRXXZ 模型，而是为识别其他系统中基态对称性随边界变化的相变提供了一种新范式。该技术依赖于设计一个对正在评估的特定相位敏感的 ansatz。

作者保持了谦逊的语气，承认未考虑噪声效应，并强调方向相干性指标是一个分析误差行为的工具，而非基本的物理定律。他们强调，该方法的成功取决于 ansatz 的特定设计以及初始态的选择。