量子计算Dirac Notation基本教学—从零基础到读懂量子信息论文（下）

第三部分：多量子比特系统

第5章 张量积与多 Qubit 态

5.1 教学目标

• 理解张量积的数学定义和物理意义
• 掌握多 qubit 态的 Dirac Notation 写法
• 学会区分可分离态与纠缠态

5.2 从单 Qubit 到多 Qubit：张量积

当系统包含多个 qubit 时，整体状态空间是各子系统空间的张量积。

数学定义 ：

设 ℋ_A 和 ℋ_B 是两个 Hilbert 空间，其张量积 ℋ_A ⊗ ℋ_B 由形式为 |a⟩ ⊗ |b⟩ 的向量张成，满足双线性性：

• (α|a₁⟩ + β|a₂⟩) ⊗ |b⟩ = α|a₁⟩⊗|b⟩ + β|a₂⟩⊗|b⟩
• |a⟩ ⊗ (γ|b₁⟩ + δ|b₂⟩) = γ|a⟩⊗|b₁⟩ + δ|a⟩⊗|b₂⟩

符号简化 ：

通常省略张量积符号 ⊗，直接并置：

|a⟩ ⊗ |b⟩ ≡ |a⟩|b⟩ ≡ |ab⟩

维度的指数增长：

• 1 qubit：2 维
• 2 qubits：2² = 4 维
• 3 qubits：2³ = 8 维
• n qubits：2ⁿ 维

这种指数增长是量子计算强大能力的来源，也是经典计算机难以模拟量子系统的根本原因。

5.3 双 Qubit 计算基

双 qubit 系统的计算基由四个向量组成：

|00⟩, |01⟩, |10⟩, |11⟩

列向量表示（按字典序排列）：

        (1)        (0)        (0)        (0)
|00⟩ = (0)  |01⟩ = (1)  |10⟩ = (0)  |11⟩ = (0)
        (0)        (0)        (1)        (0)
        (0)        (0)        (0)        (1)

一般双 Qubit 态：

|ψ⟩ = α₀₀|00⟩ + α₀₁|01⟩ + α₁₀|10⟩ + α₁₁|11⟩

归一化条件：

|α₀₀|² + |α₀₁|² + |α₁₀|² + |α₁₁|² = 1

张量积的具体计算 ：

给定两个单 qubit 态：

|a⟩ = a₀|0⟩ + a₁|1⟩
|b⟩ = b₀|0⟩ + b₁|1⟩

它们的张量积为：

|a⟩ ⊗ |b⟩ = (a₀|0⟩ + a₁|1⟩) ⊗ (b₀|0⟩ + b₁|1⟩)
          = a₀b₀|00⟩ + a₀b₁|01⟩ + a₁b₀|10⟩ + a₁b₁|11⟩

注意系数正好是乘积：α_{ij} = a_i b_j。

5.4 多 Qubit 态的 Dirac Notation 展开技巧

当处理多 qubit 态时，利用完备性关系进行展开是非常强大的技巧。

双 Qubit 完备性关系：

I ⊗ I = (|0⟩⟨0| + |1⟩⟨1|) ⊗ (|0⟩⟨0| + |1⟩⟨1|)
      = |00⟩⟨00| + |01⟩⟨01| + |10⟩⟨10| + |11⟩⟨11|

应用示例：将任意双 qubit 态 |ψ⟩ 展开到计算基：

|ψ⟩ = I⊗I |ψ⟩
    = Σ_{i,j} |ij⟩⟨ij|ψ⟩
    = Σ_{i,j} α_{ij} |ij⟩
其中 α_{ij} = ⟨ij|ψ⟩

单 Qubit 算符的张量积作用 ：

设 U 是单 qubit 算符，当它仅作用于第一个 qubit 时，整体算符为 U ⊗ I。作用于 |ab⟩：

(U ⊗ I)|ab⟩ = (U|a⟩) ⊗ (I|b⟩) = (U|a⟩) ⊗ |b⟩

线性展开法 ：

对于任意态 |ψ⟩ = Σ α_{ij}|ij⟩：

(U ⊗ I)|ψ⟩ = Σ α_{ij} (U|i⟩) ⊗ |j⟩

只需知道 U 对基态 |0⟩, |1⟩ 的作用即可计算。

5.5 可分离态与纠缠态的初步判断

定义（可分离态） ：

一个双 qubit 态 |ψ⟩ 是可分离的，当且仅当存在单 qubit 态 |a⟩, |b⟩ 使得：

|ψ⟩ = |a⟩ ⊗ |b⟩

定义（纠缠态） ：

不可分离的态称为纠缠态。

判断方法 1：系数条件

设 |ψ⟩ = α₀₀|00⟩ + α₀₁|01⟩ + α₁₀|10⟩ + α₁₁|11⟩

如果可分离，则存在 a₀,a₁,b₀,b₁ 使得：

α₀₀ = a₀b₀, α₀₁ = a₀b₁, α₁₀ = a₁b₀, α₁₁ = a₁b₁

这意味着：

α₀₀ α₁₁ = a₀b₀ a₁b₁ = a₀b₁ a₁b₀ = α₀₁ α₁₀

因此，可分离的必要条件是：

α₀₀ α₁₁ = α₀₁ α₁₀

若此式不成立，则必为纠缠态。

实例：

• |ψ₁⟩ = (|00⟩+|01⟩)/√2 = |0⟩⊗|+⟩：可分离（α₀₀=1/√2, α₀₁=1/√2, α₁₀=0, α₁₁=0, 0×0=0×1/√2 成立）
• |ψ₂⟩ = (|00⟩+|11⟩)/√2：纠缠（α₀₀=1/√2, α₁₁=1/√2, α₀₁=0, α₁₀=0, (1/√2)×(1/√2) ≠ 0×0）

5.6 多 Qubit 态的测量

完全测量 ：

在计算基下测量所有 qubit，得到结果 (i₁,i₂,...,i_n) 的概率为：

P(i₁,i₂,...,i_n) = |⟨i₁ i₂ ... i_n|ψ⟩|²

测量后态塌缩到 |i₁ i₂ ... i_n⟩。

部分测量 ：

只测量部分 qubit 时，需要使用投影算符和偏迹运算。

测量第一个 qubit 在计算基下 ：

投影算符：

• 测得 0：P₀^{(1)} = |0⟩⟨0| ⊗ I
• 测得 1：P₁^{(1)} = |1⟩⟨1| ⊗ I

概率：

P(0) = ⟨ψ| (|0⟩⟨0| ⊗ I) |ψ⟩
P(1) = ⟨ψ| (|1⟩⟨1| ⊗ I) |ψ⟩

测量后态：

若测得 0：

|ψ'⟩ = (|0⟩⟨0| ⊗ I) |ψ⟩ / √P(0)

实例：测量 Bell 态 |Φ⁺⟩ = (|00⟩+|11⟩)/√2 的第一个 qubit

P(0) = ⟨Φ⁺| (|0⟩⟨0| ⊗ I) |Φ⁺⟩ 
     = 1/2 ⟨00|0⟩⟨0|00⟩ + 1/2 ⟨11|0⟩⟨0|11⟩ 
     = 1/2 × 1 + 0 = 1/2
P(1) = 1/2

若测得 0，塌缩态：

|ψ'⟩ = (|0⟩⟨0| ⊗ I) (|00⟩+|11⟩)/√2 / √(1/2)
     = |00⟩ / √2 / (1/√2) = |00⟩

注意：测量破坏了纠缠，剩余 qubit 的状态变得确定。

5.7 练习与思考

练习 5.1

写出三 qubit 态 |010⟩ 的 8 维列向量表示。

解答 ：

按字典序 |000⟩, |001⟩, |010⟩, |011⟩, |100⟩, |101⟩, |110⟩, |111⟩

|010⟩ 是第三个基向量，所以在第 3 个位置为 1，其余为 0：

|010⟩ = (0, 0, 1, 0, 0, 0, 0, 0)^T

练习 5.2

判断态 |ψ⟩ = 1/2|00⟩ + 1/2|01⟩ + 1/2|10⟩ + 1/2|11⟩ 是否可分离。

解答 ：

系数：α₀₀=1/2, α₀₁=1/2, α₁₀=1/2, α₁₁=1/2

检查：α₀₀ α₁₁ = 1/4，α₀₁ α₁₀ = 1/4，相等。

实际上 |ψ⟩ = |+⟩ ⊗ |+⟩，是可分离态。

练习 5.3

计算 (X ⊗ Z)|01⟩。

解答：

(X ⊗ Z)|01⟩ = (X|0⟩) ⊗ (Z|1⟩) = |1⟩ ⊗ (-|1⟩) = -|11⟩

练习 5.4

对于态 |ψ⟩ = (|00⟩ + |01⟩ + |10⟩)/√3，计算测量第一个 qubit 得到 0 的概率，以及测量后的状态。

解答 ：

投影：P₀^{(1)} = |0⟩⟨0| ⊗ I

P₀^{(1)}|ψ⟩ = 1/√3 (|00⟩ + |01⟩)

概率：P(0) = ⟨ψ|P₀^{(1)}|ψ⟩ = 1/3 (1+1) = 2/3

测量后态：|ψ'⟩ = (|00⟩ + |01⟩)/√2 = |0⟩ ⊗ |+⟩

5.8 本章小结

概念	符号/公式	说明
张量积	|a⟩⊗|b⟩ ≡ |ab⟩	组合系统
双 qubit 基	|00⟩,|01⟩,|10⟩,|11⟩	4 维空间
一般态	Σ α_{ij}|ij⟩	4 个复系数
可分离条件	α₀₀α₁₁ = α₀₁α₁₀	必要但不充分（对纯态充分）
部分测量	P₀^{(1)} = |0⟩⟨0|⊗I	只测第一个 qubit

---

第6章 双 Qubit 量子门与纠缠生成

6.1 教学目标

• 掌握 CNOT 门的矩阵表示与 Dirac 表示
• 理解受控操作的通用形式
• 学会用 H + CNOT 生成 Bell 态
• 理解纠缠的本质特征

6.2 CNOT 门详解

CNOT（Controlled-NOT）是最重要的双 qubit 门，它与单 qubit 门一起构成通用门集合。

定义 ：

第一个 qubit 为控制位（control），第二个 qubit 为目标位（target）。

• 当控制位为 |0⟩ 时，目标位不变
• 当控制位为 |1⟩ 时，目标位翻转（施加 X 门）

真值表：

|00⟩ → |00⟩
|01⟩ → |01⟩
|10⟩ → |11⟩
|11⟩ → |10⟩

矩阵表示（控制位为第一 qubit，目标位为第二 qubit）：

CNOT = (1 0 0 0)
       (0 1 0 0)
       (0 0 0 1)
       (0 0 1 0)

Dirac 表示：

CNOT = |0⟩⟨0| ⊗ I + |1⟩⟨1| ⊗ X

这个形式完美体现了"条件作用"：

• 第一项：如果控制位是 |0⟩（通过投影 |0⟩⟨0| 提取），则对目标位施加 I（恒等）
• 第二项：如果控制位是 |1⟩（通过投影 |1⟩⟨1| 提取），则对目标位施加 X（翻转）

验证 Diract 表示的正确性 ：

作用于 |10⟩：

CNOT|10⟩ = (|0⟩⟨0|⊗I + |1⟩⟨1|⊗X) |10⟩
         = |0⟩⟨0|1⟩ ⊗ I|0⟩ + |1⟩⟨1|1⟩ ⊗ X|0⟩
         = 0 + |1⟩ ⊗ |1⟩ = |11⟩

正确！

CNOT 的变体：

• 控制位和目标位互换：CNOT_{2→1} = I ⊗ |0⟩⟨0| + X ⊗ |1⟩⟨1|
• 以 |0⟩ 为控制（而不是 |1⟩）：使用 X 门翻转控制位前后

6.3 受控操作的通用形式

任何单 qubit 门 U 都可以变成受控操作 Controlled-U。

定义：

Controlled-U = |0⟩⟨0| ⊗ I + |1⟩⟨1| ⊗ U

常见受控门：

• Controlled-Z (CZ)：

  CZ = |0⟩⟨0| ⊗ I + |1⟩⟨1| ⊗ Z
      = diag(1, 1, 1, -1)

  注意 CZ 是对称的：哪个是控制哪个是目标结果相同。

• Controlled-S：

  CS = |0⟩⟨0| ⊗ I + |1⟩⟨1| ⊗ S

• Controlled-Hadamard (CH)：

  CH = |0⟩⟨0| ⊗ I + |1⟩⟨1| ⊗ H

受控门与纠缠生成 ：

任何 Controlled-U 门当控制位处于叠加态时，都会在控制位和目标位之间产生纠缠（只要 U 不是恒等或相位门？实际上只要 U 不是单位阵的倍数，就会产生纠缠）。

6.4 生成 Bell 态的量子电路

Bell 态是最简单的最大纠缠态，通过 H + CNOT 电路生成。

电路：

|0⟩ ---- H ----●----
               |
|0⟩ ----------⊕----

⊕ 表示 CNOT 的目标位。

推导过程：

初始态：|00⟩

第一步：施加 H 于第一个 qubit

|ψ₁⟩ = (H ⊗ I)|00⟩ 
     = (H|0⟩) ⊗ |0⟩
     = (|0⟩+|1⟩)/√2 ⊗ |0⟩
     = (|00⟩ + |10⟩)/√2

第二步：施加 CNOT

|ψ₂⟩ = CNOT |ψ₁⟩
     = CNOT (|00⟩ + |10⟩)/√2
     = (CNOT|00⟩ + CNOT|10⟩)/√2
     = (|00⟩ + |11⟩)/√2

这就是 Bell 态 |Φ⁺⟩。

生成其他 Bell 态：

• 初始态改为 |01⟩ → 得到 |Ψ⁺⟩ = (|01⟩+|10⟩)/√2
• 初始态 |10⟩ → 先 X 再 H 再 CNOT → 得到 |Φ⁻⟩ = (|00⟩-|11⟩)/√2
• 初始态 |11⟩ → 得到 |Ψ⁻⟩ = (|01⟩-|10⟩)/√2

一般公式 ：

对于初始态 |xy⟩（x,y ∈ {0,1}），经过 H 再 CNOT 后得到：

|β_{xy}⟩ = (|0,y⟩ + (-1)^x |1, y⊕1⟩)/√2

这是四个 Bell 态的统一表达式。

6.5 纠缠态的深入分析

6.5.1 为什么 Bell 态不能分解

以 |Φ⁺⟩ = (|00⟩+|11⟩)/√2 为例，假设它可以写成张量积：

|Φ⁺⟩ = (a|0⟩ + b|1⟩) ⊗ (c|0⟩ + d|1⟩)
     = ac|00⟩ + ad|01⟩ + bc|10⟩ + bd|11⟩

比较系数：

ac = 1/√2
ad = 0
bc = 0
bd = 1/√2

从 ad=0 且 a≠0（因 ac≠0）得 d=0

从 bc=0 且 c≠0（因 ac≠0）得 b=0

但这样 bd=0，与 bd=1/√2 矛盾。

因此不可能分解。✓

6.5.2 纠缠的测量关联

对 |Φ⁺⟩ 在计算基下测量：

• 两个 qubit 都测 0：概率 1/2
• 两个 qubit 都测 1：概率 1/2
• 一个测 0 一个测 1：概率 0

这说明测量结果完全相关。但这不是经典相关，因为：

在不同基下测量 ：

如果在 Hadamard 基下测量两个 qubit：

将 |Φ⁺⟩ 用 |+⟩, |-⟩ 表示：

|Φ⁺⟩ = (|00⟩+|11⟩)/√2
     = 1/√2 [ (|+⟩+|-⟩)/√2 ⊗ (|+⟩+|-⟩)/√2 + (|+⟩-|-⟩)/√2 ⊗ (|+⟩-|-⟩)/√2 ]
     = 1/(2√2) [ |++⟩ + |+-⟩ + |-+⟩ + |--⟩ + |++⟩ - |+-⟩ - |-+⟩ + |--⟩ ]
     = 1/√2 (|++⟩ + |--⟩)

在 X 基下测量，结果也是完全相关（要么都是 +，要么都是 -）！

这种在多个互补基下都呈现完全相关是量子纠缠的独特特征，经典概率分布无法模拟。

6.5.3 纠缠与局域实在性------Bell 不等式

量子纠缠对经典直觉最深刻的挑战在于：它违反了 Bell 不等式。

CHSH 不等式 （经典关联的界限）：

对于任何局域实在理论，有：

|E(a,b) - E(a,b') + E(a',b) + E(a',b')| ≤ 2

其中 E 是关联函数。

量子力学预测 ：

对于 |Φ⁺⟩，选择适当的测量角度，可以得到：

|E(a,b) - E(a,b') + E(a',b) + E(a',b')| = 2√2 ≈ 2.828 > 2

这已经被实验反复证实，表明量子力学不能由任何局域隐变量理论复制。

6.6 其他重要双 Qubit 门

6.6.1 SWAP 门

作用：交换两个 qubit 的状态。

SWAP|ab⟩ = |ba⟩

矩阵表示：

SWAP = (1 0 0 0)
       (0 0 1 0)
       (0 1 0 0)
       (0 0 0 1)

Dirac 表示：

SWAP = |00⟩⟨00| + |01⟩⟨10| + |10⟩⟨01| + |11⟩⟨11|

CNOT 分解 ：

SWAP 可以用三个 CNOT 实现：

SWAP = CNOT₁₂ CNOT₂₁ CNOT₁₂

其中下标表示控制和目标。

6.6.2 CZ 门（Controlled-Z）

作用：当两个 qubit 都是 |1⟩ 时，施加 -1 相位因子。

CZ|11⟩ = -|11⟩
CZ|ab⟩ = |ab⟩ 对于其他情况

矩阵：diag(1, 1, 1, -1)

Dirac 表示：

CZ = |0⟩⟨0| ⊗ I + |1⟩⟨1| ⊗ Z
   = I ⊗ |0⟩⟨0| + Z ⊗ |1⟩⟨1|  （对称形式）

与 CNOT 的关系：

CZ = (I ⊗ H) CNOT (I ⊗ H)

即对目标 qubit 前后各加 H 门，CNOT 变成 CZ。

6.6.3 iSWAP 门

作用：交换并附加相位 i。

iSWAP = (1 0 0 0)
        (0 0 i 0)
        (0 i 0 0)
        (0 0 0 1)

重要性：在某些超导 qubit 系统中，iSWAP 是自然的两体相互作用。

6.7 多 Qubit 门与通用性

Toffoli 门（CCNOT） ：

三 qubit 门，两个控制位，一个目标位。当两个控制位都是 |1⟩ 时翻转目标位。

CCNOT = (|0⟩⟨0| ⊗ I - |11⟩⟨11|) ⊗ I + |11⟩⟨11| ⊗ X

重要性：Toffoli 门对于经典可逆计算是通用的；加上 H 门，对于量子计算也是通用的。

通用门集合 ：

一个门集合是通用的，如果任何幺正变换都可以用它以任意精度近似。常见的通用集合：

1. {H, T, CNOT}
2. {H, S, CNOT, T} （实际上 T 可以用 H 和 S 近似，但需要无限序列）
3. {R_x(θ), R_y(θ), R_z(θ), CNOT} （参数化门集合）

6.8 练习与思考

练习 6.1

计算 CNOT (H ⊗ I) |00⟩，并解释为什么结果不是 Bell 态。

解答：

(H ⊗ I)|00⟩ = (|0⟩+|1⟩)/√2 ⊗ |0⟩ = (|00⟩+|10⟩)/√2
CNOT (|00⟩+|10⟩)/√2 = (|00⟩+|11⟩)/√2

这确实是 Bell 态 |Φ⁺⟩。所以 H 后 CNOT 就是生成 Bell 态的标准电路。题目可能有误？

练习 6.2

证明 CNOT 是幺正且 Hermitian 的。

解答 ：

CNOT 矩阵：

(1 0 0 0)
(0 1 0 0)
(0 0 0 1)
(0 0 1 0)

显然 CNOT^T = CNOT，且只有实元素，所以 CNOT† = CNOT^T = CNOT，故 CNOT 是 Hermitian 的。

CNOT² = I（因为两次翻转等于不翻转），所以 CNOT† CNOT = CNOT² = I，故幺正。

练习 6.3

求一个电路，将可分离态 |+⟩ ⊗ |+⟩ 转换为 Bell 态 |Φ⁺⟩。

解答 ：

|+⟩ ⊗ |+⟩ = 1/2(|00⟩+|01⟩+|10⟩+|11⟩)

我们需要消除 |01⟩ 和 |10⟩ 项，这可以通过受控相位门实现？实际上，|+⟩⊗|+⟩ 不是可以通过局域操作变成 Bell 态的。因为纠缠不能通过局域操作和经典通信（LOCC）从可分离态产生。我们需要双 qubit 门。

电路：对 |+⟩⊗|+⟩ 施加 CZ 门？检查：

CZ|++⟩ = CZ 1/2(|00⟩+|01⟩+|10⟩+|11⟩)

= 1/2(|00⟩+|01⟩+|10⟩-|11⟩)

这不是 Bell 态。

其实 |+⟩⊗|+⟩ = H⊗H |00⟩，而生成 |Φ⁺⟩ 的电路是 CNOT (H⊗I)|00⟩。所以：

(H⊗H)|00⟩ → 需要某个操作变成 |Φ⁺⟩。

注意：|Φ⁺⟩ = CNOT (H⊗I) H⊗H |++⟩? 不，因为 H⊗H|00⟩ = |++⟩。

我们想从 |++⟩ 到 |Φ⁺⟩ = CNOT (H⊗I)|00⟩。但 |00⟩ = H⊗H |++⟩。

所以 |Φ⁺⟩ = CNOT (H⊗I) H⊗H |++⟩ = CNOT (I⊗H) |++⟩。

因此电路是：先施加 I⊗H（即对第二个 qubit 做 H），然后 CNOT。

验证：

(I⊗H)|++⟩ = |+⟩ ⊗ H|+⟩ = |+⟩ ⊗ |0⟩ = (|00⟩+|10⟩)/√2

CNOT 得到 (|00⟩+|11⟩)/√2 = |Φ⁺⟩。正确！

练习 6.4

写出 Controlled-Hadamard 门的 4×4 矩阵。

解答 ：

CH = |0⟩⟨0|⊗I + |1⟩⟨1|⊗H

矩阵为：

I    0
0    H

其中 0 是 2×2 零矩阵。所以：

CH = (1 0 0 0)
     (0 1 0 0)
     (0 0 1/√2 1/√2)
     (0 0 1/√2 -1/√2)

6.9 本章小结

门	矩阵	Dirac 表示	作用
CNOT	1 0 0 0; 0 1 0 0; 0 0 0 1; 0 0 1 0	|0⟩⟨0|⊗I + |1⟩⟨1|⊗X	条件翻转
CZ	diag(1,1,1,-1)	|0⟩⟨0|⊗I + |1⟩⟨1|⊗Z	条件相位
SWAP	1 0 0 0; 0 0 1 0; 0 1 0 0; 0 0 0 1	Σ|ij⟩⟨ji|	交换
Bell 态	|Φ⁺⟩ = (|00⟩+|11⟩)/√2 等	最大纠缠态	EPR 对

---

第四部分：密度矩阵与纠缠度量

第7章 密度算符形式

7.1 教学目标

• 理解纯态与混合态的区别
• 掌握密度矩阵的定义与性质
• 学会用密度矩阵计算测量概率和期望值
• 理解偏迹运算与约化密度矩阵

7.2 为什么需要密度矩阵？

到目前为止，我们一直用态向量 |ψ⟩ 描述量子系统。这种描述称为纯态描述，适用于系统处于完全确定的最大知识状态。

然而，在许多实际情况下：

• 系统处于热平衡态（温度不为零）
• 系统是更大纠缠系统的一部分（如 Bell 态中的一个 qubit）
• 我们对系统的制备存在经典不确定性（如以 50% 概率制备 |0⟩，50% 概率制备 |1⟩）

这些情况下，态向量描述不再适用，我们需要密度算符（density operator）ρ。

7.3 密度算符的定义

纯态的密度算符 ：

对于纯态 |ψ⟩，定义：

ρ = |ψ⟩⟨ψ|

混合态的密度算符 ：

若系统以概率 p_i 处于纯态 |ψ_i⟩，则密度算符为：

ρ = Σ_i p_i |ψ_i⟩⟨ψ_i|

其中 p_i ≥ 0，Σ_i p_i = 1。

密度算符的性质：

1. Hermitian：ρ† = ρ
2. 半正定：对所有 |ϕ⟩，⟨ϕ|ρ|ϕ⟩ ≥ 0
3. 迹为 1：Tr(ρ) = 1
4. 纯态判据：Tr(ρ²) = 1 当且仅当 ρ 是纯态；Tr(ρ²) < 1 对于混合态

7.4 密度矩阵的 Dirac Notation 计算

单 Qubit 示例：

纯态 |+⟩⟨+|：

|+⟩⟨+| = 1/2 (|0⟩+|1⟩)(⟨0|+⟨1|)
       = 1/2 (|0⟩⟨0| + |0⟩⟨1| + |1⟩⟨0| + |1⟩⟨1|)

矩阵形式：

ρ = 1/2 (1 1)
        (1 1)

验证 Tr(ρ) = 1/2 + 1/2 = 1，Tr(ρ²) = Tr(ρ) = 1（因为 ρ² = ρ）。

混合态 ρ = 1/2 |0⟩⟨0| + 1/2 |1⟩⟨1|：

ρ = 1/2 (1 0) + 1/2 (0 0) = (1/2 0)
        (0 0)       (0 1)   (0 1/2)

验证 Tr(ρ) = 1，Tr(ρ²) = 1/4 + 1/4 = 1/2 < 1，确实是混合态。

混合态与叠加态的本质区别：

• 叠加态 |+⟩ = (|0⟩+|1⟩)/√2：密度矩阵有非对角项（相干项），表示量子相干性
• 混合态 1/2(|0⟩⟨0|+|1⟩⟨1|)：密度矩阵只有对角项，表示经典概率混合

7.5 密度矩阵下的测量

投影测量的概率 ：

对于投影算符 P_m，测量结果 m 的概率为：

p(m) = Tr(P_m ρ)

测量后状态 ：

若测得结果 m，系统塌缩到：

ρ' = (P_m ρ P_m) / p(m)

期望值 ：

可观测量 A 的期望值为：

⟨A⟩ = Tr(A ρ)

实例：对 ρ = 1/2 |0⟩⟨0| + 1/2 |1⟩⟨1| 在计算基下测量

P(0) = Tr(|0⟩⟨0| ρ) = ⟨0|ρ|0⟩ = 1/2
P(1) = Tr(|1⟩⟨1| ρ) = ⟨1|ρ|1⟩ = 1/2

若测得 0，ρ' = (|0⟩⟨0| ρ |0⟩⟨0|) / (1/2) = |0⟩⟨0|。

7.6 偏迹与约化密度矩阵

当复合系统处于态 ρ_AB，我们想要描述子系统 A 的局域状态时，使用偏迹（partial trace）运算。

定义 ：

ρ_A = Tr_B(ρ_AB)

其中 Tr_B 表示对子系统 B 的基求迹：

Tr_B(|a₁⟩⟨a₂| ⊗ |b₁⟩⟨b₂|) = |a₁⟩⟨a₂| · Tr(|b₁⟩⟨b₂|) = |a₁⟩⟨a₂| · ⟨b₂|b₁⟩

线性扩展 ：

对于一般算符 O_AB = Σ_i A_i ⊗ B_i，

Tr_B(O_AB) = Σ_i A_i · Tr(B_i)

重要例子：Bell 态的约化密度矩阵

对于 |Φ⁺⟩ = (|00⟩+|11⟩)/√2，密度矩阵为：

ρ_AB = |Φ⁺⟩⟨Φ⁺| 
     = 1/2 (|00⟩⟨00| + |00⟩⟨11| + |11⟩⟨00| + |11⟩⟨11|)

对 B 求偏迹：

ρ_A = Tr_B(ρ_AB)
    = 1/2 (|0⟩⟨0|·⟨0|0⟩ + |0⟩⟨1|·⟨1|0⟩ + |1⟩⟨0|·⟨0|1⟩ + |1⟩⟨1|·⟨1|1⟩)
    = 1/2 (|0⟩⟨0| + |1⟩⟨1|)
    = I/2

这正是最大混合态！这揭示了一个深刻的事实：

在最大纠缠态中，单个 qubit 的状态是完全不确定的（最大混合态），但整体状态却是完全确定的（纯态）。

这是纠缠的典型特征：整体纯，局部混。

7.7 纠缠的密度矩阵判据

纯态可分离判据 ：

纯态 |ψ⟩_AB 可分离 ⇔ ρ_A = Tr_B(|ψ⟩⟨ψ|) 是纯态（即 Tr(ρ_A²)=1）

PPT 判据 （Positive Partial Transpose，部分转置正定性）：

对于两体系统，对其中一个子系统做部分转置，若结果矩阵不是半正定的，则原态必为纠缠态。

对于双 qubit 系统，PPT 判据是充要条件。

部分转置操作 ：

对子系统 B 的基取转置：

ρ^{T_B} = (I ⊗ T)(ρ)

其中 T 是转置运算（在固定基下）。

示例：对 |Φ⁺⟩⟨Φ⁺| 做部分转置

ρ = 1/2 (|00⟩⟨00| + |00⟩⟨11| + |11⟩⟨00| + |11⟩⟨11|)
ρ^{T_B} = 1/2 (|00⟩⟨00| + |01⟩⟨10| + |10⟩⟨01| + |11⟩⟨11|)

求本征值：1/2, 1/2, 1/2, -1/2。有负本征值，故非半正定，所以是纠缠态。

7.8 纠缠熵与纠缠度量

纠缠熵 （Entanglement Entropy）：

对于纯态 |ψ⟩_AB，定义纠缠熵为 von Neumann 熵：

S(ρ_A) = -Tr(ρ_A log₂ ρ_A)

其中 ρ_A 是约化密度矩阵。

性质：

• 可分离态：S = 0
• Bell 态：ρ_A = I/2，本征值为 {1/2, 1/2}，S = -1/2 log₂(1/2) - 1/2 log₂(1/2) = 1
• 最大纠缠态达到最大值 log₂(d)，d 是子系统维度（对 qubit 为 1）

并发度 （Concurrence）：

对于双 qubit 混合态，一个常用的纠缠度量是并发度 C：

C(ρ) = max(0, λ₁ - λ₂ - λ₃ - λ₄)

其中 λ_i 是矩阵 √(√ρ ρ̃ √ρ) 的本征值按降序排列，ρ̃ = (Y⊗Y)ρ*(Y⊗Y)。

7.9 练习与思考

练习 7.1

计算纯态 |ψ⟩ = (|0⟩ + i|1⟩)/√2 的密度矩阵，并验证 Tr(ρ²) = 1。

解答：

ρ = |ψ⟩⟨ψ| = 1/2 (1) (1, -i) = 1/2 (1 -i)
         (i)               (i 1)

注意 ⟨ψ| = (1/√2)(⟨0| - i⟨1|) = (1/√2, -i/√2)

所以外积为 1/2 (1)(1, -i) 和 1/2 (i)(1, -i) 等。矩阵是：

ρ = 1/2 (1 -i)
        (i 1)

Tr(ρ) = 1/2 + 1/2 = 1

ρ² = 1/4 (1 -i)(1 -i) = 1/4 (1 -i+i+1?) 仔细算：

(i 1)(i 1) (i-1? 等)

ρ² = 1/4 (1×1 + (-i)×i, 1×(-i)+(-i)×1) = 1/4 (1+1, -i-i) = 1/4 (2 -2i) = 1/2 (1 -i) = ρ

(i×1 + 1×i, i×(-i)+1×1) (i+i, 1+1) (2i 2) (i 1)

所以 ρ² = ρ，Tr(ρ²) = Tr(ρ) = 1。✓

练习 7.2

对于混合态 ρ = 3/4 |0⟩⟨0| + 1/4 |1⟩⟨1|，求在 X 基下测量得到 + 的概率。

解答 ：

P(+) = Tr(|+⟩⟨+| ρ)

|+⟩⟨+| = 1/2 (1 1)

(1 1)

ρ = (3/4 0)

(0 1/4)

P(+) = Tr(1/2 (1 1)(3/4 0)) = 1/2 Tr(3/4 1/4) = 1/2 (3/4 + 1/4) = 1/2

(1 1)(0 1/4) (3/4 1/4)

实际上任何对角基的混合在互补基下都是等概率的。

练习 7.3

对于态 |Ψ⁻⟩ = (|01⟩ - |10⟩)/√2，计算约化密度矩阵 ρ_A，并求其 von Neumann 熵。

解答 ：

ρ_AB = |Ψ⁻⟩⟨Ψ⁻| = 1/2 (|01⟩⟨01| - |01⟩⟨10| - |10⟩⟨01| + |10⟩⟨10|)

ρ_A = Tr_B(ρ_AB) = 1/2 (|0⟩⟨0|·⟨1|1⟩ + |1⟩⟨1|·⟨0|0⟩) = 1/2 (|0⟩⟨0| + |1⟩⟨1|) = I/2

S(ρ_A) = -1/2 log₂(1/2) - 1/2 log₂(1/2) = 1

练习 7.4

判断以下密度矩阵是否代表纯态：

ρ = (1/2 1/4)

(1/4 1/2)

解答 ：

计算 Tr(ρ²)：

ρ² = (1/2 1/4)(1/2 1/4) = (1/4+1/16 1/8+1/8) = (5/16 1/4)

(1/4 1/2)(1/4 1/2) (1/8+1/8 1/16+1/4) (1/4 5/16)

Tr(ρ²) = 5/16 + 5/16 = 5/8 = 0.625 < 1，所以是混合态。

7.10 本章小结

概念	公式	性质
纯态密度矩阵	ρ = |ψ⟩⟨ψ|	Tr(ρ²)=1
混合态密度矩阵	ρ = Σ p_i |ψ_i⟩⟨ψ_i|	Tr(ρ²)<1
测量概率	p(m) = Tr(P_m ρ)	Born 规则的推广
期望值	⟨A⟩ = Tr(A ρ)	可观测量平均
偏迹	ρ_A = Tr_B(ρ_AB)	子系统状态
von Neumann 熵	S(ρ) = -Tr(ρ log ρ)	纠缠度量（纯态）

---

第五部分：量子算法基础与论文阅读

第8章 Dirac Notation 在量子算法中的应用

8.1 教学目标

• 理解量子并行性的数学表述
• 掌握 Deutsch-Jozsa 算法的 Dirac 推导
• 学会用量子电路图与 Dirac 公式相互转换
• 能够阅读简单量子算法的论文片段

8.2 量子并行性

量子计算的一个核心优势是量子并行性：通过对叠加态应用线性算符，我们同时计算了函数在所有可能输入上的值。

数学表述 ：

设有函数 f: {0,1} → {0,1}，我们可以构造一个幺正算符 U_f，使得：

U_f |x⟩|y⟩ = |x⟩|y ⊕ f(x)⟩

若输入寄存器处于叠加态 (|0⟩+|1⟩)/√2，输出寄存器初始化为 |0⟩：

U_f ( (|0⟩+|1⟩)/√2 ⊗ |0⟩ ) 
= 1/√2 ( U_f|00⟩ + U_f|10⟩ )
= 1/√2 ( |0,f(0)⟩ + |1,f(1)⟩ )

这个态同时包含了 f(0) 和 f(1) 的信息！

Dirac 推导：

U_f (H⊗I) |00⟩ = U_f 1/√2 (|00⟩ + |10⟩)
                = 1/√2 (|0, f(0)⟩ + |1, f(1)⟩)

注意：虽然态同时包含了多个函数值，但直接测量只能得到一个结果。量子算法的艺术在于通过干涉将需要的全局信息"集中"到可测量的输出上。

8.3 Deutsch-Jozsa 算法

问题：给定一个函数 f: {0,1}^n → {0,1}，已知它要么是常数函数（对所有输入输出相同），要么是平衡函数（对一半输入输出 0，一半输出 1）。判断 f 是哪一类。

经典复杂度：最坏情况需要查询 2^{n-1}+1 次。

量子算法：只需 1 次查询！

电路：

|0⟩^⊗n ---- H^⊗n ---- U_f ---- H^⊗n ---- 测量
|1⟩ ---------- H --------⊕------------------ 丢弃

推导（n=1 情况）：

初始态：|ψ₀⟩ = |01⟩

第一层 H 门：

|ψ₁⟩ = (H ⊗ H) |01⟩
     = (H|0⟩) ⊗ (H|1⟩)
     = |+⟩ ⊗ |-⟩
     = 1/2 (|0⟩+|1⟩) ⊗ (|0⟩-|1⟩)

应用 U_f：

回忆 U_f |x⟩|y⟩ = |x⟩|y ⊕ f(x)⟩。

对于 |x⟩ ⊗ (|0⟩-|1⟩)/√2：

U_f |x⟩ (|0⟩-|1⟩)/√2 
= 1/√2 (|x⟩|f(x)⟩ - |x⟩|1⊕f(x)⟩)

若 f(x)=0，得 |x⟩(|0⟩-|1⟩)/√2 = (-1)^0 |x⟩(|0⟩-|1⟩)/√2

若 f(x)=1，得 |x⟩(|1⟩-|0⟩)/√2 = (-1)^1 |x⟩(|0⟩-|1⟩)/√2

所以一般地：

U_f |x⟩ (|0⟩-|1⟩)/√2 = (-1)^{f(x)} |x⟩ (|0⟩-|1⟩)/√2

这个技巧称为相位反冲（phase kickback）。

因此：

|ψ₂⟩ = U_f |ψ₁⟩ 
     = 1/2 [ (-1)^{f(0)}|0⟩ + (-1)^{f(1)}|1⟩ ] ⊗ (|0⟩-|1⟩)

最后对第一个 qubit 施加 H：

|ψ₃⟩ = (H ⊗ I) |ψ₂⟩

计算 H 作用：

若 f(0)=f(1)（常数函数）：

(-1)^{f(0)} [ H|0⟩ + H|1⟩ ] /2 ⊗ (|0⟩-|1⟩)/√2
= (-1)^{f(0)} (√2 |0⟩) /2 ⊗ (|0⟩-|1⟩)/√2
= (-1)^{f(0)} |0⟩ ⊗ (|0⟩-|1⟩)/√2

测量第一个 qubit 必得 0。

若 f(0)≠f(1)（平衡函数）：

不失一般性设 f(0)=0, f(1)=1：

1/2 [ |0⟩ - |1⟩ ] ⊗ (|0⟩-|1⟩) （忽略归一化）
H(|0⟩-|1⟩)/2 = √2 |1⟩ /2

测量第一个 qubit 必得 1。

因此，单次查询即可确定 f 是常数还是平衡！

n 比特推广 ：

结果为：测量得到全 0 当且仅当 f 是常数函数，否则得到其他结果。

8.4 常见量子算法的 Dirac Notation 片段

8.4.1 Grover 算法的核心迭代

Grover 迭代由两个反射组成：

G = (2|ψ⟩⟨ψ| - I) (I - 2|t⟩⟨t|)

其中 |ψ⟩ = H^⊗n|0⟩⊗n 是均匀叠加态，|t⟩ 是目标态。

在 Dirac Notation 下，Oracle 的作用写为：

O = I - 2|t⟩⟨t|

它给目标态加上 -1 相位，其他态不变。

扩散算符：

D = 2|ψ⟩⟨ψ| - I

迭代 k 次后，状态振幅可由几何旋转计算，最终以高概率测得目标态。

8.4.2 量子 Fourier 变换

QFT 在计算基下的作用：

QFT |j⟩ = 1/√N Σ_{k=0}^{N-1} e^{2π i j k / N} |k⟩

其中 N = 2^n。

对于单 qubit，QFT 就是 H 门：

H|0⟩ = (|0⟩+|1⟩)/√2
H|1⟩ = (|0⟩-|1⟩)/√2

Shor 算法中的周期查找 ：

利用 QFT 将周期性信息转化为可测量的峰值，从而找出函数的周期，进而分解整数。

8.5 阅读量子信息论文的基础能力

现在你已经具备了阅读量子计算基础论文所需的核心符号体系。让我们看几个典型的论文片段。

片段 1（态的定义）：

We initialize the system in the state |ψ₀⟩ = |0⟩^⊗n |1⟩.

解读：n 个 qubit 初始化为 |0⟩，最后一个 qubit 初始化为 |1⟩。

片段 2（纠缠生成）：

Applying a Hadamard on the first qubit followed by a CNOT yields the Bell state |Φ⁺⟩ = (|00⟩+|11⟩)/√2.

解读：这正是我们学过的标准 Bell 态生成电路。

片段 3（测量）：

The probability of outcome x is given by p(x) = ⟨ψ| P_x |ψ⟩, where P_x = |x⟩⟨x| is the projector onto the computational basis state |x⟩.

解读：Born 规则的投影算符形式。

片段 4（纠缠度量）：

For the pure state |ψ⟩_AB, the entanglement entropy is E(ψ) = S(ρ_A), where ρ_A = Tr_B(|ψ⟩⟨ψ|) is the reduced density matrix of subsystem A.

解读：偏迹和 von Neumann 熵用于度量纠缠。

片段 5（门序列）：

The circuit implements the unitary U = (H ⊗ I) CNOT (T ⊗ S).

解读：门从右向左应用：先 T 和 S 分别作用，然后 CNOT，最后 H 作用在第一个 qubit。

8.6 从公式到电路：相互转换

给定电路，写出门序列：

|0⟩ ---- H ----●---- H ----
               |
|0⟩ ----------⊕------------

门序列：U = (H ⊗ I) CNOT (H ⊗ I)

给定公式，画出电路 ：

U = (H ⊗ H) CZ (I ⊗ H)

电路：

|ψ₁⟩ ---- H ----●----
                |
|ψ₂⟩ ---- H ----Z---- H ----

（注意 CZ 可以用 ● 和 Z 表示，或者用 ● 表示控制，目标为 Z）

8.7 常见计算技巧汇总

在推导中反复使用的技巧：

技巧 1：插入完备性关系

|ψ⟩ = Σ_i |i⟩⟨i|ψ⟩

技巧 2：利用线性性

U Σ_i c_i |ψ_i⟩ = Σ_i c_i U|ψ_i⟩

技巧 3：相位反冲

U_f |x⟩(|0⟩-|1⟩) = (-1)^{f(x)} |x⟩(|0⟩-|1⟩)

技巧 4：Hadamard 作用在计算基

H^⊗n |x⟩ = 1/√(2^n) Σ_{y} (-1)^{x·y} |y⟩

其中 x·y 是逐位内积模 2。

技巧 5：Pauli 算符的乘积规则

X Z = -Z X
X Y = i Z (等)

技巧 6：偏迹的"掉入"技巧

对于可分离算符 A⊗B，Tr_B(A⊗B) = A·Tr(B)

8.8 练习与思考

练习 8.1

用量子并行性计算 f(x)=x 对于叠加态 (|0⟩+|1⟩)/√2 的结果（使用 U_f 作用于 |x⟩|0⟩）。

解答：

U_f ( (|0⟩+|1⟩)/√2 ⊗ |0⟩ ) 
= 1/√2 (U_f|00⟩ + U_f|10⟩)
= 1/√2 (|0, 0⊕0⟩ + |1, 0⊕1⟩)
= 1/√2 (|00⟩ + |11⟩)

这正是 Bell 态 |Φ⁺⟩。

练习 8.2

推导：对于任意 n，H^⊗n |0⟩^⊗n = 1/√(2^n) Σ_{x=0}^{2n-1} |x⟩。

解答 ：

单 qubit：H|0⟩ = 1/√2 (|0⟩+|1⟩) = 1/√2 Σ_{x=0}^{1} |x⟩

两个 qubit：H^⊗2 |00⟩ = (H|0⟩)⊗(H|0⟩) = 1/2 (|0⟩+|1⟩)⊗(|0⟩+|1⟩) = 1/2 Σ_{x,y} |xy⟩

一般地，n 重张量积产生 2^n 项，每项系数为 1/√(2^n)。

练习 8.3

证明：CNOT (X ⊗ I) CNOT = X ⊗ X。

解答 ：

利用 CNOT 的 Dirac 形式和 X 的比特翻转性质。

直观：CNOT 将第一个 qubit 的 X 错误传播到第二个 qubit？实际上这是纠错码中的一个恒等式。可以用矩阵乘法验证，或利用 CNOT 的性质：

CNOT (|1⟩⟨0|⊗I) CNOT = ...

这里从略，但结果成立。

练习 8.4

阅读以下论文片段并解释：

"The initial state is prepared as |ψ⟩ = ⊗{i=1}^n |+i⟩. After applying the unitary U_φ = Π{i<j} CZ{ij}, the state becomes a graph state |G⟩."

解答 ：

初始态是所有 qubit 的 |+⟩ 态的张量积。然后施加一系列 CZ 门（每对 qubit 之间根据图的边连接），结果得到图态 |G⟩。图态是一类重要的多体纠缠态，用于测量基量子计算。

8.9 本章小结

算法/概念	关键 Dirac 表述	核心思想
量子并行性	U_f Σ|x⟩|0⟩ = Σ|x⟩|f(x)⟩	同时计算所有输入
相位反冲	U_f |x⟩|-⟩ = (-1)^{f(x)}|x⟩|-⟩	函数值编码到相位
Deutsch-Jozsa	测量 H^⊗n 的结果判函数类型	单次查询解决指数经典问题
Grover	G = (2|ψ⟩⟨ψ|-I)(I-2|t⟩⟨t|)	振幅放大
QFT	QFT|j⟩ = 1/√N Σ e^{2πijk/N}|k⟩	周期查找的基础

---

第六部分：综合应用与进阶主题

第9章 多体纠缠与量子纠错入门

9.1 教学目标

• 理解 GHZ 态和 W 态的区别
• 了解图态与 stabilizer 形式
• 初步掌握量子纠错的基本思想
• 能够用 Dirac Notation 描述简单的纠错码

9.2 三 Qubit 纠缠态：GHZ 与 W

对于三 qubit，存在两种不等价的最大纠缠态。

GHZ 态（Greenberger-Horne-Zeilinger）：

|GHZ⟩ = (|000⟩ + |111⟩)/√2

性质：

• 测量任意一个 qubit 在计算基下，若结果为 0，则其他两个 qubit 塌缩到 |00⟩；若为 1，则塌缩到 |11⟩
• 在 X 基下测量，三 qubit 结果之间存在特定的奇偶关联
• 丢失任意一个 qubit，剩余两个 qubit 的状态变为完全可分离的混合态

W 态：

|W⟩ = (|001⟩ + |010⟩ + |100⟩)/√3

性质：

• 测量任意一个 qubit，以 2/3 概率得到 0，1/3 概率得到 1
• 若测得 0，剩余两 qubit 处于纠缠态 (|01⟩+|10⟩)/√2
• 若测得 1，剩余两 qubit 处于可分离态 |00⟩
• 丢失任意一个 qubit，剩余两个 qubit 仍然保持一定的纠缠

GHZ 与 W 的不可互转性 ：

在 LOCC（局域操作和经典通信）下，GHZ 态和 W 态不能互相转换。它们代表了两种不同类型的多体纠缠结构。

9.3 Stabilizer 形式简介

对于多 qubit 态，直接用振幅描述变得指数困难。Stabilizer 形式提供了一种高效的描述方法，适用于一大类重要的量子态（stabilizer states）。

定义 ：

一个 stabilizer 群 S 是一组相互对易的 Pauli 算符乘积（且不包含 -I），满足：

• 对所有 s ∈ S，s|ψ⟩ = |ψ⟩
  即 |ψ⟩ 是所有 s 的本征值为 +1 的本征态。

示例：Bell 态 |Φ⁺⟩ 的 stabilizer 为：

⟨ X⊗X, Z⊗Z ⟩

验证：

X⊗X (|00⟩+|11⟩)/√2 = (|11⟩+|00⟩)/√2 = |Φ⁺⟩

Z⊗Z (|00⟩+|11⟩)/√2 = (|00⟩+(-1)(-1)|11⟩)/√2 = |Φ⁺⟩

GHZ 态的 stabilizer：

⟨ X⊗X⊗X, Z⊗Z⊗I, Z⊗I⊗Z ⟩

（实际上只需 n 个独立生成元）

Stabilizer 形式的优势：

• 态的描述只需 O(n) 个生成元，而不是 O(2ⁿ) 个振幅
• 幺正演化在 stabilizer 下对应共轭作用，计算效率高
• Gottesman-Knill 定理：stabilizer 电路可在经典计算机上高效模拟

9.4 量子纠错基础

量子态极易受环境噪声干扰。量子纠错码通过冗余编码保护量子信息。

经典重复码的量子对应：三 qubit 比特翻转码

编码 ：

逻辑 |0⟩ 编码为物理 |000⟩

逻辑 |1⟩ 编码为物理 |111⟩

即：|ψ_L⟩ = α|000⟩ + β|111⟩

错误模型：每个 qubit 独立地以概率 p 发生 X 错误（比特翻转）。

错误症状测量 ：

测量 stabilizer 生成元 Z₁Z₂ 和 Z₂Z₃。

• 若无错误：Z₁Z₂ = +1, Z₂Z₃ = +1
• 若第 1 qubit 翻转：Z₁Z₂ = -1, Z₂Z₃ = +1
• 若第 2 qubit 翻转：Z₁Z₂ = -1, Z₂Z₃ = -1
• 若第 3 qubit 翻转：Z₁Z₂ = +1, Z₂Z₃ = -1

根据症状，施加相应的纠正操作（X 门）即可恢复原态。

为什么测量 stabilizer 不破坏叠加？

因为逻辑态是 stabilizer 的本征态，测量 stabilizer 不会区分 |0_L⟩ 和 |1_L⟩（它们本征值相同），从而保护了叠加。

相位翻转码 ：

通过基变换，将相位翻转（Z 错误）转化为比特翻转：

|+⟩ = (|0⟩+|1⟩)/√2，|-⟩ = (|0⟩-|1⟩)/√2

编码：|0_L⟩ = |+++⟩, |1_L⟩ = |---⟩

Z 错误在 |+⟩, |-⟩ 基下表现为 X 错误，因此可以用同样的症状测量纠正。

Shor 码 ：

结合比特翻转码和相位翻转码，可以纠正任意单 qubit 错误。它用 9 个物理 qubit 编码 1 个逻辑 qubit。

9.5 容错量子计算概念

量子纠错的目标不仅是存储，还要在编码态上进行计算。

横向操作 ：

某些逻辑门可以通过对每个物理 qubit 独立施加相同的门来实现。例如：

• 横向 H：H_L = H ⊗ H ⊗ ... ⊗ H
• 横向 CNOT：CNOT_L = CNOT ⊗ CNOT ⊗ ... ⊗ CNOT

非横向操作 ：

某些门（如 T 门）不能横向实现，需要借助魔法态蒸馏等技巧。

容错阈值定理 ：

如果每个物理门的错误率低于某个阈值（约 10⁻⁴ ~ 10⁻²），则可以通过级联编码任意延长量子计算时间。

9.6 练习与思考

练习 9.1

验证 GHZ 态是 stabilizer 生成元 X⊗X⊗X 和 Z⊗Z⊗I, Z⊗I⊗Z 的 +1 本征态。

解答 ：

|GHZ⟩ = (|000⟩+|111⟩)/√2

X⊗X⊗X (|000⟩+|111⟩)/√2 = (|111⟩+|000⟩)/√2 = |GHZ⟩，所以本征值 +1。

Z⊗Z⊗I|000⟩ = |000⟩，Z⊗Z⊗I|111⟩ = (-1)(-1)|111⟩ = |111⟩，所以也是 +1 本征值。

同理 Z⊗I⊗Z。

练习 9.2

写出三 qubit 比特翻转码的编码电路。

解答：

|ψ⟩ ----●----●----
        |    |
|0⟩ ----⊕----|----
             |
|0⟩ ---------⊕----

即两个 CNOT，控制位为待编码 qubit，目标位为两个辅助 qubit 初始化为 |0⟩。

练习 9.3

若三 qubit 编码态 α|000⟩+β|111⟩ 发生第二个 qubit 翻转，写出错误态和症状。

解答 ：

错误态：X₂ (α|000⟩+β|111⟩) = α|010⟩+β|101⟩

症状：Z₁Z₂ = ⟨010|Z₁Z₂|010⟩? 直接计算：Z₁Z₂|010⟩ = (1)(-1)|010⟩ = -|010⟩，Z₁Z₂|101⟩ = (-1)(1)|101⟩ = -|101⟩，所以本征值为 -1。

Z₂Z₃|010⟩ = (-1)(1)|010⟩ = -|010⟩，Z₂Z₃|101⟩ = (1)(-1)|101⟩ = -|101⟩，所以本征值为 -1。

症状为 (-1, -1)，指示第 2 qubit 错误。

练习 9.4

解释为什么 stabilizer 测量不会导致波函数塌缩到编码基的某个特定态。

解答 ：

因为逻辑基态 |0_L⟩ 和 |1_L⟩ 都是所测 stabilizer 的同一本征值 +1 的本征态。测量 stabilizer 只会将态投影到该本征空间，而不会区分 |0_L⟩ 和 |1_L⟩。因此叠加态 α|0_L⟩+β|1_L⟩ 在测量后仍然保持叠加，只是可能被投影到错误子空间（如果发生了错误）。纠正后恢复原叠加。

9.7 本章小结

概念	关键公式/描述	重要性
GHZ 态	(|000⟩+|111⟩)/√2	最大纠缠，但脆弱
W 态	(|001⟩+|010⟩+|100⟩)/√3	对丢失鲁棒
Stabilizer	S |ψ⟩ = |ψ⟩	高效描述多体态
三 qubit 码	|0_L⟩=|000⟩, |1_L⟩=|111⟩	纠正一个 X 错误
症状测量	测量 stabilizer 生成元	检测错误不破坏态
横向操作	U_L = U^⊗n	容错实现逻辑门

---

第七部分：总结与展望

第10章 Dirac Notation 掌握清单与进阶路线

10.1 核心知识体系回顾

经过以上九个章节的学习，你现在应该已经掌握了以下核心能力：

第一层：符号识别

• 能够区分 |ψ⟩（ket）、⟨ψ|（bra）、⟨ϕ|ψ⟩（内积）、|ϕ⟩⟨ψ|（外积）
• 理解这些符号对应的线性代数对象：列向量、行向量、标量、矩阵

第二层：基展开与计算

• 熟练写出单 qubit 态在计算基和 Hadamard 基下的展开
• 能够计算内积、外积、张量积
• 掌握 Born 规则的概率计算

第三层：量子门作用

• 知道 Pauli 门、H 门、S 门、T 门对基态的作用
• 能够利用线性性推导门对任意态的作用
• 理解 CNOT 的 Dirac 表示及其条件作用本质

第四层：纠缠分析

• 能够判断一个双 qubit 纯态是否可分离
• 理解 Bell 态的生成电路和性质
• 知道纠缠态的整体-局域关系（约化密度矩阵）

第五层：算法阅读

• 能够用量子并行性解释算法的加速
• 理解 Deutsch-Jozsa 算法的推导
• 能够阅读包含 Dirac Notation 的简单论文片段

第六层：进阶概念

• 理解密度矩阵和偏迹
• 知道 stabilizer 形式的基本思想
• 了解量子纠错的编码-症状-纠正流程

10.2 最重要的三条计算规则（重温）

无论推导多么复杂，最终都归结为以下三条：

规则 1：线性性

U (α|ψ⟩ + β|ϕ⟩) = α U|ψ⟩ + β U|ϕ⟩

规则 2：正交归一消元

⟨i|j⟩ = δ_{ij}

规则 3：张量积分步作用

(A ⊗ B)(|ψ⟩ ⊗ |ϕ⟩) = (A|ψ⟩) ⊗ (B|ϕ⟩)

10.3 常见误区与避免方法

误区	正确理解
叠加态 = 概率混合	叠加态有相干项，可以干涉；混合态没有
全局相位可忽略，相对相位也可忽略	相对相位至关重要，影响干涉结果
纠缠 = 经典相关	纠缠在互补基下也呈现关联，违反 Bell 不等式
忘记归一化	计算概率或期望值时务必保持归一化
Bra 与 Ket 的内积顺序随意	⟨ϕ|ψ⟩ = ⟨ψ|ϕ⟩*，交换需取复共轭

10.4 进一步学习路线图

初级→中级：

1. 完成 Nielsen & Chuang 前四章的所有习题
2. 手算 Deutsch-Jozsa、Grover（n=2）、量子隐形传态的完整推导
3. 用 Python + Qiskit 实现简单电路，验证理论计算

中级→高级：

1. 学习 stabilizer 形式与 Gottesman-Knill 定理
2. 理解 Shor 算法的数论部分和周期查找量子部分
3. 掌握量子纠错码的构造（Steane 码、Shor 码）
4. 阅读原始论文：Shor 1994, Grover 1996, Bennett & Brassard 1984

高级→研究级：

1. 深入学习容错量子计算与阈值定理
2. 研究特定的量子硬件平台（超导、离子阱、光量子）的噪声模型
3. 学习张量网络方法（MPS、PEPS）处理多体量子系统
4. 跟踪量子优势实验的最新进展

10.5 推荐资源

教材：

• Nielsen & Chuang, Quantum Computation and Quantum Information（圣经）
• Kaye, Laflamme & Mosca, An Introduction to Quantum Computing
• Mermin, Quantum Computer Science（侧重计算机科学视角）

在线课程：

• Qiskit Textbook (qiskit.org/textbook)
• MIT 8.370 Quantum Computation (EdX)
• Perimeter Institute 量子信息课程 (YouTube)

论文数据库：

• arXiv.org quant-ph
• Quantum (quantum-journal.org)

10.6 最后的话

Dirac Notation 是量子信息科学的通用语言。掌握它就像掌握了阅读量子世界地图的钥匙。

从最初困惑于 |ψ⟩ 是什么，到现在能够流畅地写出 Bell 态的生成、理解量子算法的推导、判断纠缠的存在，你已经走过了量子计算入门最艰难也最重要的一段路。

记住：量子计算的反直觉性不会消失，但数学会给你坚实的立足点。当直觉失效时，相信 Dirac Notation，相信线性代数。

Ket 表示态，Bra 表示对偶，内积给振幅，外积给算符，门是线性变换，多 qubit 用张量积，纠缠就是无法拆成张量积的联合态。

继续前进，量子世界等待你的探索！⚛️

---

附录 A：常用符号速查表

符号	含义	示例
|ψ⟩	量子态 ket	|0⟩, |+⟩, |ψ⟩
⟨ψ|	量子态 bra	⟨0|, ⟨+|
⟨ϕ|ψ⟩	内积（振幅）	⟨0|1⟩ = 0
|ϕ⟩⟨ψ|	外积（算符）	|0⟩⟨1|
⊗	张量积	|0⟩⊗|1⟩ = |01⟩
U	幺正算符	H, X, CNOT
U†	共轭转置	H† = H
Tr	迹	Tr(ρ) = 1
ρ	密度算符	ρ = |ψ⟩⟨ψ|
I	恒等算符	I|ψ⟩ = |ψ⟩
H	Hadamard 门	H|0⟩ = |+⟩
X, Y, Z	Pauli 门	X|0⟩ = |1⟩
S, T	相位门	S|1⟩ = i|1⟩
CNOT	受控非门	CNOT|10⟩ = |11⟩
CZ	受控 Z 门	CZ|11⟩ = -|11⟩

---

附录 B：标准量子态表

态名称	Dirac 表示	性质
|0⟩	(1,0)^T	计算基
|1⟩	(0,1)^T	计算基
|+⟩	(|0⟩+|1⟩)/√2	X 基
|-⟩	(|0⟩-|1⟩)/√2	X 基
|+i⟩	(|0⟩+i|1⟩)/√2	Y 基
|-i⟩	(|0⟩-i|1⟩)/√2	Y 基
|Φ⁺⟩	(|00⟩+|11⟩)/√2	Bell 态
|Φ⁻⟩	(|00⟩-|11⟩)/√2	Bell 态
|Ψ⁺⟩	(|01⟩+|10⟩)/√2	Bell 态
|Ψ⁻⟩	(|01⟩-|10⟩)/√2	Bell 态
|GHZ⟩	(|000⟩+|111⟩)/√2	三体纠缠
|W⟩	(|001⟩+|010⟩+|100⟩)/√3	三体纠缠

---

附录 C：练习题答案提示

因篇幅所限，各章练习题答案已嵌入相应章节的"练习与思考"部分。建议学习者先独立尝试，再对照答案验证。

---

文档信息

• 版本：1.0
• 适用对象：量子计算初学者至中级学习者
• 前置要求：线性代数基础（向量、矩阵、内积）、复数运算

---