1. 课程导入与框架建立
量子纠缠熵与关联度量计算，研究的是一件很具体的事：给定一个量子态，如何用可计算的数值刻画它内部的量子关联，并解释这些数值的物理含义。这句话里包含三个必须分开对待的对象——量子态本身、刻画关联的度量、以及把度量算出来的流程。初学者最常见的困惑，往往不是公式难，而是把“纠缠熵”“互信息”“关联函数”混成一类东西，导致拿到一个态不知道该算什么、算出来的数又代表什么。本节先把这三者的定位摆清楚，再给出贯穿全课程的计算框架。
三个核心对象的分工
先看纠缠熵。它针对的是“一个子系统”，而不是整个态。设整个系统分为 A、B 两部分，对总态取偏迹得到 A 的约化密度矩阵 ρA，纠缠熵就是 ρA 的冯·诺依曼熵：
S(ρ_A) = -Tr(ρ_A log₂ ρ_A) = -Σ_i λ_i log₂ λ_i其中 λi 是 ρA 的本征值。它的关键性质是：当且仅当 A 与 B 之间没有纠缠（总态为直积态）时 S=0；纯态下 S(ρA)=S(ρB)，所以“A 的纠缠熵”同时也是“A 与 B 之间纠缠量”的度量。注意它只对纯态全局态才严格等于纠缠度量，混合态的纠缠熵还混入了经典不确定性，这一点后面会专门处理。
再看互信息。它衡量的是 A 与 B 之间“总共”的关联，既包括量子纠缠，也包括经典关联：
I(A:B) = S(ρ_A) + S(ρ_B) - S(ρ_AB)互信息恒非负，且对直积态为零。它比纠缠熵“宽”：一个只有经典关联的可分离态，纠缠熵为零但互信息可以大于零。因此在判断“有没有量子纠缠”时，互信息只能作为上界式的粗筛，不能单独作为纠缠判据。
第三是关联函数。它不直接对密度矩阵求熵，而是选取一组可观测量（如泡利算符 σx、σz），测量它们在不同子系统上的联合期望值：
C_ij = ⟨O_i^A ⊗ O_j^B⟩ - ⟨O_i^A⟩⟨O_j^B⟩关联函数是实验上最容易直接测的量，也是判断纠缠（如通过贝尔不等式、纠缠判据）的入口。它和熵的关系是：熵是态的整体信息量，关联函数是态在特定测量方向上的投影。二者互补，不能互相替代。
为什么必须区分这三者
举一个具体的教学示例（以下数值为按定义手算的示例，非实测）。取两量子比特的 Werner 态 ρ = p|Φ⁺⟩⟨Φ⁺| + (1−p)I/4，其中 |Φ⁺⟩=(|00⟩+|11⟩)/√2，p∈[0,1]。当 p=1 时是纯的贝尔态，对 A 取偏迹得 ρA=I/2，本征值均为 1/2，于是 S(ρA)=−2×(1/2)log₂(1/2)=1 bit，这就是最大纠缠。当 p=0.5 时，ρA 仍为 I/2，纠缠熵依旧是 1 bit，但此时态已是混合态，1 bit 并不等于“纠缠量为 1”。若直接套用纯态公式就会误判。正确做法是改用可分离判据或形成纠缠度等混合态度量。这个例子说明：同一个熵值在不同态类型下含义不同，必须先判断态的类型再选度量。
再看互信息。对上述 p=0.5 的态，S(ρA)=S(ρB)=1，而 S(ρAB) 需由四维密度矩阵本征值计算，结果小于 2，因此 I(A:B)>0。但该态在 p≤1/3 时可分离，互信息却不为零——这正是“互信息包含经典关联”的直接体现。所以互信息适合做关联总量的粗描述，不适合单独下“有纠缠”的结论。
本课程的计算流程
把上面三条路径统一起来，实际计算按以下顺序推进，每一步都有明确的输入输出：
- 写出密度矩阵：明确系统维数、基矢顺序（如 |00⟩,|01⟩,|10⟩,|11⟩），把态写成矩阵形式。基矢顺序一旦确定，全文与图示保持一致。
- 确定二分方式：是 A 与 B 两体，还是多体中取某一子区域。二分不同，熵值不同，必须先声明。
- 求约化密度矩阵：对不需要的子系统取偏迹，得到 ρA。
- 求本征值并算熵：对 ρA 做本征分解，代入 S=−Σλilog₂λi；本征值为零的项按 0log0=0 处理。
- 算互信息或关联函数：按需选择路径二或路径三，得到交叉验证用的第二个数。
- 与判据比对并解释：用 S=0、I≥S、可分离判据等条件检查结果是否自洽，再回到物理含义。
学习目标与后续路线
完成本课程后，你应当能够：独立完成两体及多体子区域的纠缠熵手算；正确区分纯态与混合态下熵的物理含义；用互信息和关联函数对熵结果做交叉验证；并能把计算结果与可分离判据、贝尔型判据对应起来解释。后续章节将依次展开密度矩阵与偏迹的具体运算、本征值求解与数值稳定性处理、混合态纠缠度量的选择、多体子区域熵的面积律与标度行为，以及关联函数在实验数据上的拟合与判读。本节建立的地图，就是这些内容共同的坐标系。