3. 纠缠的量化与判据
前两节我们学会了如何用门电路把多个量子比特制备到纠缠态。但“制备出来了”并不等于“知道它有多纠缠”。在实验和算法中，我们往往需要回答三个层次的问题：这个态到底纠缠不纠缠？纠缠有多强？它属于哪一类多体纠缠？本节围绕这三问，梳理主流的量化度量与判据，并给出可操作的判断路径。
一、两比特纠缠的度量：熵与 concurrence
对纯态而言，最自然的度量是纠缠熵。把两比特系统 AB 分成 A 和 B 两部分，对 B 求偏迹得到约化密度矩阵 ρA，则纠缠熵定义为 S = −Tr(ρA log₂ρA)。对两比特纯态，ρA 的两个本征值可写成 (1±√(1−C²))/2，其中 C 就是 concurrence。于是 S 与 C 一一对应：C=0 时 S=0（直积态），C=1 时 S=1（最大纠缠，如 Bell 态）。这说明对两比特纯态，熵和 concurrence 是等价的度量，只是后者计算更直接。
concurrence 的定义为 C = |⟨ψ*|σy⊗σy|ψ⟩|，对纯态就是两比特系数构成的行列式的模的两倍。例如 |ψ⟩ = α|00⟩ + β|11⟩，则 C = 2|αβ|。当 α=β=1/√2 时 C=1；当 β=0 时 C=0。对混态，concurrence 需借助“自旋翻转”矩阵 ρ̃ = (σy⊗σy)ρ*(σy⊗σy) 计算：先求 ρρ̃ 的四个本征值 λi（降序），则 C = max(0, √λ₁−√λ₂−√λ₃−√λ₄)。这是实验数据处理中最常用的公式之一。
另一类度量是形成纠缠 EF，它定义为该态可分解为纯态系综时平均纠缠熵的最小值。对两比特态，EF 是 concurrence 的单调函数：EF = h((1+√(1−C²))/2)，其中 h(x)=−x log₂x−(1−x)log₂(1−x)。因此实际工作中常先算 C，再换算 EF，二者对两比特纠缠的排序一致。
二、PPT 判据与纠缠目击者
当系统维度升高或处于混态时，concurrence 不再够用。此时最常用的必要判据是 PPT 判据（正部分转置判据）：若 ρ 可分，则其部分转置 ρTB 必为半正定，即所有本征值非负。反之，若部分转置后出现负本征值，则 ρ 一定纠缠。对 2×2 和 2×3 系统，PPT 还是充分条件；对更高维或更多体，PPT 只是必要条件，存在“PPT 纠缠态”（束缚纠缠），需借助其他判据。
PPT 判据需要知道完整密度矩阵，实验上往往只能测有限个观测量。这时用纠缠目击者（entanglement witness）：构造一个厄米算符 W，使得对所有可分态都有 Tr(Wρsep) ≥ 0，而存在某个纠缠态使 Tr(Wρ) < 0。测量到负期望值即“目击”到纠缠。最经典的例子是 Bell 态 |Φ⁺⟩=(|00⟩+|11⟩)/√2 的目击者 W = 1/2·I − |Φ⁺⟩⟨Φ⁺|，只需测 XX、YY、ZZ 三个关联量即可构造：⟨W⟩ = 1/2 − (⟨XX⟩−⟨YY⟩+⟨ZZ⟩)/4。若测得 ⟨W⟩ < 0，就证明存在纠缠，无需层析。
目击者的优点是测量量少、可在实验上直接执行；缺点是“漏检”——某个目击者只能捕捉它针对的那一类纠缠，未目击不等于可分。因此实践中常组合多个目击者，或先用目击者快速筛查，再用层析加 PPT 做确认。
三、多体纠缠的分类：GHZ、W 与簇态
三比特及以上，纠缠不再是“有或无”的单一量，而出现 inequivalent 的纠缠类——不同类之间无法通过局域操作和经典通信（LOCC）相互转换。最典型的两类是 GHZ 态与 W 态。
GHZ 态 |GHZ⟩=(|000⟩+|111⟩)/√2 是“全有或全无”型纠缠：任意一个比特被测量后，其余两比特塌缩到直积态，即单比特约化密度矩阵完全混态（S=1），但任意两比特的约化态却是可分的。它的纠缠“脆弱”但全局性强，是量子纠错、GHZ 悖论检验和计量学的核心资源。W 态 |W⟩=(|001⟩+|010⟩+|100⟩)/√3 则是“一激发共享”型：丢失一个比特后，剩余两比特仍保留部分纠缠，因此对粒子丢失更鲁棒。实验上常通过比较“去掉一个比特后剩余态的纠缠是否为零”来区分二者。
簇态（cluster state）是另一类重要的多体纠缠，典型为线性四比特簇态，由 Hadamard 与 CZ 门作用在 |+⟩⊗n 上生成。它满足一组本征方程 Ki|C⟩ = |C⟩，其中 Ki = Xi∏j∈N(i)Zj 为稳定子。簇态是测量型量子计算（MBQC）的通用资源：只需对部分比特做单比特测量，就能在剩余比特上“驱动”出任意量子线路。与 GHZ 的全局纠缠不同，簇态的纠缠分布在相邻比特之间，具有局域结构，因而更适合作为可扩展的计算资源。
| 纠缠类型 | 典型形式 | 关键特征 | 典型用途 |
|---|---|---|---|
| GHZ 态 | (|000⟩+|111⟩)/√2 | 单比特约化完全混，任意两体可分 | 计量、纠错、基础检验 |
| W 态 | (|001⟩+|010⟩+|100⟩)/√3 | 丢失一比特后仍保留纠缠 | 鲁棒纠缠传输 |
| 簇态 | H⊗n 后接 CZ 链 | 稳定子结构，局域纠缠 | 测量型量子计算 |
四、实操：一次完整的纠缠判定流程
假设我们在超导平台上制备了一个两比特态，想判断它是否纠缠并估计纠缠度。可执行如下流程（示例参数为合理假设，非实测）：
- 测量关联量：在 XX、YY、ZZ 三组基下各做约 2000 次采样，得到 ⟨XX⟩≈0.62、⟨YY⟩≈−0.58、⟨ZZ⟩≈0.60（示例值）。
- 构造目击者：⟨W⟩ = 1/2 − (⟨XX⟩−⟨YY⟩+⟨ZZ⟩)/4 = 0.5 − (0.62+0.58+0.60)/4 = 0.5 − 0.45 = 0.05 > 0，未目击。注意这并不证明可分，只是该目击者未捕捉到。
- 层析重建：做 9 组基（3×3）测量，用最大似然估计重建 ρ，再计算部分转置本征值。若最小本征值为 −0.03（示例），则 PPT 判据判定纠缠。
- 量化：由 ρ 算 concurrence，若得 C≈0.35，则 EF = h((1+√(1−0.1225))/2) ≈ h(0.968) ≈ 0.20 ebit（示例换算）。
这个流程体现了“快速筛查—严格判定—定量估计”的分层思路：目击者便宜但会漏检，PPT 严格但需要完整层析，concurrence 与形成纠缠给出可比较的数值。对多体系统，则先通过“去掉一个比特后剩余态是否仍纠缠”区分 GHZ 与 W，再用稳定子测量验证簇态结构。