当前位置：首页 > news >正文

量子计算中的对称保护拓扑序：理论与硬件实现

news 2026/5/8 21:27:19

1. 量子计算中的对称保护拓扑序：从理论到硬件实现

在量子多体物理领域，对称保护拓扑序（Symmetry-Protected Topological order，SPT）代表了一类超越传统朗道范式的量子物质新分类。这类物质状态的特殊之处在于：它们无法用任何局域序参量来区分，却展现出非平凡的全局纠缠特性和受对称性保护的边缘态。理解并操控这类量子态，不仅对基础物理研究至关重要，也为拓扑量子计算等前沿应用提供了可能。

1.1 SPT序的核心物理特征

传统相变理论（朗道范式）认为，物质的不同相可以通过对称性破缺和相应的局域序参量来区分。然而，SPT相的发现彻底改变了这一认知——这些相保持所有对称性，却仍能通过全局的拓扑不变量来区分。其核心特征包括：

非局域弦序参数：在自旋链中表现为长程的弦关联函数，例如对于偶数长度的弦，测量算符为 $S^E_l = (-1)^{l/2}\langle \prod_{i=1}^l Z_i \rangle$。这种关联在SPT相中趋于非零常数，而在普通相中则指数衰减到零。
纠缠谱简并：当系统被分割时，约化密度矩阵的本征值呈现特定的简并模式。例如在偶数Haldane相中，切割J0键产生的纠缠谱呈现偶数重简并，而切割J1键则无此特征。
对称性保护的边缘态：在开边界条件下，系统边界会出现受对称性保护的局域化边缘态。这些态对保持对称性的微扰具有鲁棒性，是SPT相的"指纹"。

这些特性使得SPT材料在量子信息领域具有独特价值，例如可用于实现低耗散的量子互连和拓扑保护的量子记忆单元。

1.2 量子模拟的挑战与机遇

尽管冷原子、里德堡模拟器等平台已成功演示了小规模SPT系统，但这些方法受限于物理体系固有的耦合方式和测量手段。数字量子计算机提供了全新的可能性：

哈密顿量的自由编程：可以模拟任意设计的相互作用，而不仅限于特定材料中自然存在的耦合。
非局域测量的灵活性：能够直接测量弦序参数等传统实验难以获取的非局域观测量。
动态过程的精确控制：可研究量子淬火等非平衡过程，探索SPT相的动力学稳定性。

然而，实现这些优势面临重大技术挑战。量子硬件有限的相干时间和门保真度，要求我们必须找到制备SPT态的高效方法——既要保证状态质量，又要尽可能减少所需的量子门数量。这正是本研究突破的关键所在。

2. 张量网络与近似量子编译技术

2.1 从DMRG到量子电路的桥梁

密度矩阵重整化群（DMRG）是研究一维量子多体系统的黄金标准，能够高效计算系统的基态矩阵乘积态（MPS）表示。然而，直接将MPS转化为量子电路面临两个主要问题：

电路深度爆炸：传统方法需要与系统尺寸成比例的电路深度，远超当前量子硬件的容错能力。
纠缠捕获不足：简单电路结构难以有效表达SPT态中的非局域纠缠特征。

本研究采用的解决方案是**基于张量网络的近似量子编译（AQC）**技术。其核心思想是将MPS制备问题转化为一个优化问题：寻找一组量子门参数，使得制备出的态与目标MPS的保真度最大化。

2.2 AQC-Tensor算法的实现细节

具体实现流程如下：

MPS压缩预处理：首先使用DMRG计算目标哈密顿量的基态MPS，典型键维数在19-40之间。然后通过变分压缩将键维数降至5-8，同时保持99.9%以上的保真度。这一步骤大幅降低了后续优化的计算复杂度。
砖块电路ansatz设计：采用L层最近邻砖块结构作为变分ansatz。每层由作用于所有偶数对(i,i+1)的两比特门和随后作用于奇数对的门组成。这种结构在保持浅深度的同时，能有效捕获短程纠缠。
张量网络辅助优化：利用AQC-Tensor算法高效计算成本函数（1-保真度）及其梯度。通过将量子电路和张量网络表示相结合，避免了直接模拟大规模量子系统的计算开销。
电路深度优化：对不同层数L进行扫描，寻找在给定深度下能达到的最高保真度。对于100位点系统，最优电路深度在18-39个CNOT层之间，对应总CNOT门数在891-1932个。

表1比较了四个典型参数点的编译结果：

状态标记	J0	J1	原始键维数	压缩后键维数	编译保真度	CNOT深度
O½	0.5	1	22	5	98.9%	18
E½	1	0.5	19	5	98.9%	21
E-1	1	-1	26	8	99.0%	21
E-2	1	-2	40	8	97.9%	39

技术细节：对于最难编译的E-2态（J0=1, J1=-2），其较大的关联长度（约3个晶格常数）导致需要更深的电路来捕获纠缠结构。这与理论预期一致——制备一般平移不变MPS所需电路深度与关联长度成正比。

3. 量子硬件实验与SPT特征验证

3.1 实验设置与误差缓解

实验在IBM Pittsburgh超导量子处理器上执行，主要技术挑战包括：

串扰抑制：通过精心安排量子门时序和频率调谐，减少邻近比特间的非预期耦合。
误差缓解技术组合：
- Pauli-twirling：平均化系统误差
- TREX（Twirled Readout Error Extinction）：校正测量误差
- ZNE（Zero Noise Extrapolation）：通过人为引入额外噪声并外推至零噪声极限，估计理想结果

特别对于弦序测量，ZNE表现出色。如图3所示，通过在不同噪声强度下测量并线性外推，有效抑制了系统误差，使20位点弦序参数的测量精度达到理论值的90%以上。

3.2 SPT特征的系统性验证

3.2.1 弦序参数的测量

对于E-2态（偶数Haldane相），实验观测到：

偶数弦序SE：在l=20时仍保持约0.4的值，表现出典型的长程序
奇数弦序SO：随长度增加快速衰减至零

这与理论预期完美吻合——在偶数Haldane相中，只有偶数弦序参数才是真正的序参量。通过平均链上5个不同区间的测量结果（s=20,30,...,60），进一步抑制了统计波动。

3.2.2 边缘态的直接观测

在O½态（奇数Haldane相）中，观测到清晰的边缘磁化现象：

链两端约5个位点范围内存在显著的净磁化（≈0.25）
磁化强度向体内指数衰减，关联长度ξ≈1.8（与DMRG计算的1.97吻合）
体内磁化趋于零，符合SPT相无局域序的特征

通过拟合Sz = (⟨Z₂ᵢ⟩+⟨Z₂ᵢ₊₁⟩)/2 ∼ e^(-x/ξ)，可提取精确的关联长度。这一测量直接证实了对称性保护的边缘态存在。

3.2.3 纠缠谱诊断

通过量子态层析技术，重构了6位点区块的约化密度矩阵，并分析其本征值分布。关键发现：

E-2态：
- 切割J0键：前两个本征值接近简并（λ₁≈λ₂≈0.5）
- 切割J1键：仅一个主导本征值（λ₁≈1）
O½态：
- 表现出互补模式——J1键切割产生简并谱
- 随着区块增大，简并度逐渐显现

这种交替简并模式是SPT序最本质的特征，无法用任何局域扰动消除。尽管硬件噪声导致本征值不完全简并，但整体模式已足够诊断SPT序的存在。

4. 技术挑战与解决方案实录

4.1 编译过程中的典型问题

在实际操作中，我们遇到了几个关键挑战：

梯度消失问题：
- 现象：优化过程中成本函数梯度变得极小，导致收敛停滞
- 诊断：参数空间存在平坦区域，特别是当初始猜测远离解时
- 解决方案：采用"热启动"策略——先用较小系统尺寸（如20位点）优化，再将参数作为大系统的初始值
局部最优陷阱：
- 案例：E-2态在L=6层时多次收敛到95%保真度
- 突破：引入随机扰动跳出局部最优，最终达到97.9%
- 经验：对高关联长度态，需要更多优化尝试和耐心
测量误差放大：
- 问题：弦序参数涉及多个Pauli-Z乘积，测量方差随长度增加
- 技巧：采用分组测量策略，将长弦拆分为重叠短段进行联合估计