ENTCALC: Toolkit for calculating geometric entanglement in multipartite quantum systems

技术摘要：“ENTCALC：用于计算多体量子系统中几何纠缠的工具包”

问题陈述
量化多体系统中的量子纠缠是一项计算密集型任务。虽然存在多种纠缠度量（例如，纠缠熵、负度/negativity），但它们往往面临局限性：有些仅适用于纯态，另一些则无法检测某些纠缠态（如正部分转置/PPT 纠缠态），且判定可分性是一个 NP-hard 问题。具体而言，由于计算一般混合多体态的几何纠缠涉及对可分态集的优化，这一过程非常困难。现有方法通常将分析限制在纯态，或仅侧重于两体子系统（忽略全局多体结构），或者依赖于无法提供定量度量的纠缠见证（entanglement witnesses）。因此，需要一种能够为纯态和混合多体态提供高精度几何纠缠估计值的工具，特别是在解析解不可用的情况下。

方法论
作者提出了 ENTCALC，这是一个通过 Python (entcalcpy) 和 MATLAB (entcalc) 实现的软件工具包。该工具包通过计算几何纠缠 ($E_G$) 的上下界来估计其值，从而确定真实值所在的区间。

1. 几何纠缠定义：

   • 对于纯态 $|\psi\rangle$，$E_G$ 定义为 $1 - \max_{|\phi\rangle \in S} |\langle\phi|\psi\rangle|^2$，其中 $S$ 是全可分态的集合。
   • 对于混合态 $\rho$，它通过凸组合扩展（convex-roof extension）定义为：$E_G(\rho) = \min_{\{p_i, |\psi_i\rangle\}} \sum p_i E_G(|\psi_i\rangle)$。
   • 等价形式使用量子保真度（quantum fidelity）：$E_G(\rho) = 1 - \max_{\sigma \in S} F(\rho, \sigma)$。

2. 纯态算法：

   • SDP 方法： 该程序使用基于 PPT（正部分转置）松弛的半正定规划（SDP）公式。对于三体态，下界通过最大化与 PPT 态的重叠来推导。同时计算一个带有误差估计（$4\sqrt{\epsilon}$）的上界，其中 $\epsilon$ 与最优 PPT 态的最大特征值相关。
   • 梯度下降： 使用迭代算法寻找最接近的可分态，从而提供一个上界。该方法也被用作混合态的子程序。

3. 混合态算法：

   • 上界： 通过梯度下降算法计算，该算法通过迭代优化目标态的可分分解来最大化保真度。
   • 下界： 该程序实现了四种不同的下界方法，允许用户在准确度与计算成本之间进行权衡：
     • 下界 1 (SDP)： 基于在 PPT 态集上最大化保真度。效率高，但对于 PPT 纠缠态会产生零值。
     • 下界 2 (k-对称扩展)： 通过引入 $k$-对称扩展来提高准确度，构建了一个随着 $k \to \infty$ 向可分态集收敛的更紧密的层次结构。
     • 下界 3 (纯度-纠缠互补性)： 源于热力学与纠缠之间的关系，利用混合态的纯化（purification）并对辅助变量进行优化。
     • 下界 4 (全局纯度-纠缠互补性)： 与下界 3 类似，但在整个纯化系统上运行，以更高的计算维度提供更紧的下界。

核心贡献

• 双界估计： 不同于以往通常只提供单一界限或见证的方法，ENTCALC 同时计算几何纠缠的上下界。这使得识别出真实几何纠缠所在的一个狭窄区间成为可能，显著提高了估计的置信度。
• 处理混合态与 PPT 态： 该工具包能够量化包括 PPT 纠缠态（束缚纠缠）在内的混合态纠缠，而这类纠缠是基于负度（negativity）的度量无法检测到的。
• 新颖的理论扩展： 作者扩展了现有的理论框架，推导出了基于纯度-纠缠互补关系的、针对多体场景的新型下界。
• 软件实现： 发布了 Python 和 MATLAB 的开源软件包，利用了标准科学库（NumPy, QuTiP, CVXPY, SciPy, CVX, QETLAB），使高精度纠缠计算变得触手可及。

结果与应用
作者将 ENTCALC 应用于几个具有代表性的示例，观察到上下界之间的差距始终非常小（通常 $< 10^{-5}$），表明了极高的准确度：

1. PPT 纠缠态： 该程序成功量化了 $3 \otimes 3$ Horodecki PPT 纠缠态的几何纠缠，其界限差异约为 $5.15 \times 10^{-6}$。
2. 低秩态： 在 3、4 和 5 量子比特系统中，对 GHZ 态和 W 态的混合物进行了纠缠量化。
3. 自旋链中的热纠缠：
   • XX 模型： 即使在两体纠缠为零的情况下，该程序也检测到了三量子比特海森堡 XX 链热态中的多体纠缠。它识别出了铁磁耦合（$J>0$）和反铁磁耦合（$J<0$）下的束缚纠缠（PPT 纠缠）区域。
   • 带磁场的 XXX 模型： 工具包通过观察几何纠缠随磁场强度和逆温度变化的间断，捕捉到了量子相变的特征。
   • 六角形自旋链： 作者展示了通过调节外部磁场，可以激活非相邻位点（例如 6 量子比特环中的自旋 B、D 和 F）之间的纠缠，这些位点并不直接相互作用，而这一现象是两体分析无法捕捉到的。
4. 噪声通道： 该程序量化了 GHZ 态和 W 态对振幅阻尼（amplitude damping）和去极化噪声（depolarizing noise）的鲁棒性。结果证实，在高噪声水平下，W 态通常比 GHZ 态对振幅阻尼更具鲁棒性，而在低噪声水平下，GHZ 态最初表现出更高的纠缠度。

重要性与声明
论文声称，ENTCALC 是第一个允许对多体混合态进行高精度几何纠缠计算的软件包。通过提供上下界之间的紧密区间，该工具为解析解难以处理的理论研究提供了一种可靠的数值方法。作者强调，测试案例中极小的界限间隙证实了所采用界限的紧密性。他们将该程序定位为研究复杂现象（如自旋链中的束缚纠缠、量子相变以及噪声环境中的纠缠动力学）的宝贵资源，并指出它可以揭示两体度量所遗漏的特征（如非局域多体纠缠）。这项工作并不声称解决了通用的 NP-hard 可分性问题，而是为特定的、具有物理意义的态提供了一种实用的、高精度的近似方法。