上一节的 CNOT 造出了 Bell 态 ∣ 00 ⟩ + ∣ 11 ⟩ 2 \frac{\ket{00}+\ket{11}}{\sqrt2} 2 ∣ 00 ⟩ + ∣ 11 ⟩ ,
我们证明了它写不成乘积。但「写不成乘积」只是一句定性判决。
把参数稍微拧一下:
∣ ψ 1 ⟩ = ∣ 00 ⟩ + ∣ 11 ⟩ 2 , ∣ ψ 2 ⟩ = 0.99 ∣ 00 ⟩ + 1 − 0.99 2 ∣ 11 ⟩ (9.3.1) \ket{\psi_1}=\frac{\ket{00}+\ket{11}}{\sqrt2},\qquad
\ket{\psi_2}=0.99\ket{00}+\sqrt{1-0.99^2}\,\ket{11}\tag{9.3.1} ∣ ψ 1 ⟩ = 2 ∣ 00 ⟩ + ∣ 11 ⟩ , ∣ ψ 2 ⟩ = 0.99 ∣ 00 ⟩ + 1 − 0.9 9 2 ∣ 11 ⟩ ( 9.3.1 )
两个都写不成乘积,都是纠缠态。可 ∣ ψ 2 ⟩ \ket{\psi_2} ∣ ψ 2 ⟩ 几乎就是 ∣ 00 ⟩ \ket{00} ∣ 00 ⟩ ——
直觉说它的纠缠「淡得多」。直觉需要变成数字:
9.4 节的 Bell 实验、9.5 节的隐形传态都要消耗纠缠,
消耗多少、够不够用,必须可计算 。
先把地基打牢:可分离与纠缠
两比特纯态的一般形式是
∣ ψ ⟩ = α ∣ 00 ⟩ + β ∣ 01 ⟩ + γ ∣ 10 ⟩ + δ ∣ 11 ⟩ (9.3.2) \ket{\psi}=\alpha\ket{00}+\beta\ket{01}+\gamma\ket{10}+\delta\ket{11}\tag{9.3.2} ∣ ψ ⟩ = α ∣ 00 ⟩ + β ∣ 01 ⟩ + γ ∣ 10 ⟩ + δ ∣ 11 ⟩ ( 9.3.2 )
若它能写成 ∣ a ⟩ ⊗ ∣ b ⟩ \ket{a}\otimes\ket{b} ∣ a ⟩ ⊗ ∣ b ⟩ ,称为可分离态 (各过各的日子);
否则称为纠缠态 。判据其实一行就够:乘积态的系数矩阵
( α β γ δ ) \begin{pmatrix}\alpha&\beta\\\gamma&\delta\end{pmatrix} ( α γ β δ ) 秩为 1,即
α δ − β γ = 0 ⟺ 可分离 (9.3.3) \alpha\delta-\beta\gamma=0
\quad\Longleftrightarrow\quad \text{可分离}\tag{9.3.3} α δ − β γ = 0 ⟺ 可分离 ( 9.3.3 )
这个行列式不为零就是纠缠——它甚至暗示了「纠缠量」应该和 ∣ α δ − β γ ∣ |\alpha\delta-\beta\gamma| ∣ α δ − β γ ∣
的大小有关。我们最后会回到这一点。
旧工具哪里不够:波函数「局部读数」失灵
第 3 章的态矢量描述整个系统。但纠缠态里,单独一个 qubit 根本没有自己的态矢量 ——
你找不到任何 ∣ a ⟩ \ket{a} ∣ a ⟩ 能正确预言只测 A 粒子的所有实验。
这正是 3.10 节 引入密度算符的动机:
局部的一切都装在约化密度矩阵 里,
ρ A = T r B ∣ ψ ⟩ ⟨ ψ ∣ (9.3.4) \rho_A=\mathrm{Tr}_B\,\ket{\psi}\bra{\psi}\tag{9.3.4} ρ A = Tr B ∣ ψ ⟩ ⟨ ψ ∣ ( 9.3.4 )
(对 B 求偏迹 = 把 B 的所有可能结果加权平均掉。)关键的观察是:
Bell 态的约化密度矩阵:整体纯、局部最混 基础 ~6 min 展开 目标 :对 ∣ Φ + ⟩ = ∣ 00 ⟩ + ∣ 11 ⟩ 2 \ket{\Phi^+}=\frac{\ket{00}+\ket{11}}{\sqrt2} ∣ Φ + ⟩ = 2 ∣ 00 ⟩ + ∣ 11 ⟩ 求 ρ A \rho_A ρ A 。
第一步:写出整体密度算符。
∣ Φ + ⟩ ⟨ Φ + ∣ = 1 2 ( ∣ 00 ⟩ ⟨ 00 ∣ + ∣ 00 ⟩ ⟨ 11 ∣ + ∣ 11 ⟩ ⟨ 00 ∣ + ∣ 11 ⟩ ⟨ 11 ∣ ) (9.3.5) \ket{\Phi^+}\bra{\Phi^+}
=\frac12\Big(\ket{00}\bra{00}+\ket{00}\bra{11}+\ket{11}\bra{00}+\ket{11}\bra{11}\Big)\tag{9.3.5} Φ + ⟩ ⟨ Φ + = 2 1 ( ∣ 00 ⟩ ⟨ 00 ∣ + ∣ 00 ⟩ ⟨ 11 ∣ + ∣ 11 ⟩ ⟨ 00 ∣ + ∣ 11 ⟩ ⟨ 11 ∣ ) ( 9.3.5 ) 第二步:对 B 求偏迹。 规则:T r B \mathrm{Tr}_B Tr B 只保留 B 侧「左右标签相同」的项——
∣ i ⟩ A ∣ j ⟩ B ⟨ k ∣ A ⟨ l ∣ B \ket{i}_A\ket{j}_B\bra{k}_A\bra{l}_B ∣ i ⟩ A ∣ j ⟩ B ⟨ k ∣ A ⟨ l ∣ B 的偏迹是 δ j l ∣ i ⟩ A ⟨ k ∣ A \delta_{jl}\ket{i}_A\bra{k}_A δ j l ∣ i ⟩ A ⟨ k ∣ A 。逐项看:
∣ 00 ⟩ ⟨ 00 ∣ \ket{00}\bra{00} ∣ 00 ⟩ ⟨ 00 ∣ :B 侧是 ∣ 0 ⟩ ⟨ 0 ∣ \ket{0}\bra{0} ∣ 0 ⟩ ⟨ 0 ∣ ,迹为 1,留下 ∣ 0 ⟩ ⟨ 0 ∣ A \ket{0}\bra{0}_A ∣ 0 ⟩ ⟨ 0 ∣ A ;
∣ 00 ⟩ ⟨ 11 ∣ \ket{00}\bra{11} ∣ 00 ⟩ ⟨ 11 ∣ :B 侧是 ∣ 0 ⟩ ⟨ 1 ∣ \ket{0}\bra{1} ∣ 0 ⟩ ⟨ 1 ∣ ,迹为 0,整项消失 ;
∣ 11 ⟩ ⟨ 00 ∣ \ket{11}\bra{00} ∣ 11 ⟩ ⟨ 00 ∣ :同理消失;
∣ 11 ⟩ ⟨ 11 ∣ \ket{11}\bra{11} ∣ 11 ⟩ ⟨ 11 ∣ :留下 ∣ 1 ⟩ ⟨ 1 ∣ A \ket{1}\bra{1}_A ∣ 1 ⟩ ⟨ 1 ∣ A 。
ρ A = 1 2 ∣ 0 ⟩ ⟨ 0 ∣ + 1 2 ∣ 1 ⟩ ⟨ 1 ∣ = I 2 (9.3.6) \rho_A=\frac12\ket{0}\bra{0}+\frac12\ket{1}\bra{1}=\frac{I}{2}\tag{9.3.6} ρ A = 2 1 ∣ 0 ⟩ ⟨ 0 ∣ + 2 1 ∣ 1 ⟩ ⟨ 1 ∣ = 2 I ( 9.3.6 ) 第三步:读出物理。 ρ A = I / 2 \rho_A=I/2 ρ A = I /2 是最大混合态 ——上一节 Bloch 球的球心。
只看 A 粒子,沿任何轴测都是五五开,与「一枚被抛过但没看的硬币」在统计上无法区分。
交叉项(相干项)恰好在偏迹中被 B 侧的正交性杀掉了:
整体的相干搬进了关联里,局部只剩下噪声。
对比可分离态 ∣ a ⟩ ⊗ ∣ b ⟩ \ket{a}\otimes\ket{b} ∣ a ⟩ ⊗ ∣ b ⟩ :偏迹直接给 ρ A = ∣ a ⟩ ⟨ a ∣ \rho_A=\ket{a}\bra{a} ρ A = ∣ a ⟩ ⟨ a ∣ ,仍是纯态。
给「混」定量:Schmidt 分解与纠缠熵
任何两体纯态都能通过奇异值分解写成Schmidt 形式 :
∣ ψ ⟩ = ∑ k λ k ∣ u k ⟩ A ∣ v k ⟩ B , λ k ≥ 0 , ∑ k λ k = 1 (9.3.7) \ket{\psi}=\sum_k \sqrt{\lambda_k}\,\ket{u_k}_A\ket{v_k}_B,
\qquad \lambda_k\ge0,\ \sum_k\lambda_k=1\tag{9.3.7} ∣ ψ ⟩ = k ∑ λ k ∣ u k ⟩ A ∣ v k ⟩ B , λ k ≥ 0 , k ∑ λ k = 1 ( 9.3.7 )
其中 ∣ u k ⟩ \ket{u_k} ∣ u k ⟩ 、∣ v k ⟩ \ket{v_k} ∣ v k ⟩ 各是 A、B 侧的一组正交基。对两比特,k k k 至多取 2。
好处立竿见影:
ρ A = ∑ k λ k ∣ u k ⟩ ⟨ u k ∣ , ρ B = ∑ k λ k ∣ v k ⟩ ⟨ v k ∣ (9.3.8) \rho_A=\sum_k\lambda_k\ket{u_k}\bra{u_k},\qquad
\rho_B=\sum_k\lambda_k\ket{v_k}\bra{v_k}\tag{9.3.8} ρ A = k ∑ λ k ∣ u k ⟩ ⟨ u k ∣ , ρ B = k ∑ λ k ∣ v k ⟩ ⟨ v k ∣ ( 9.3.8 )
两侧约化密度矩阵的本征值相同 ,都是 Schmidt 系数 λ k \lambda_k λ k 。
只有一个非零 λ \lambda λ ⇔ 乘积态;两个非零 λ \lambda λ ⇔ 纠缠。
把混合程度用第 3 章的 von Neumann 熵来称重:
E ( ψ ) ≡ S ( ρ A ) = − T r ρ A log 2 ρ A = − ∑ k λ k log 2 λ k (9.3.9) E(\psi)\equiv S(\rho_A)=-\mathrm{Tr}\,\rho_A\log_2\rho_A
=-\sum_k\lambda_k\log_2\lambda_k\tag{9.3.9} E ( ψ ) ≡ S ( ρ A ) = − Tr ρ A log 2 ρ A = − k ∑ λ k log 2 λ k ( 9.3.9 )
这就是纠缠熵 (以 2 为底,单位是比特)。三个基准点:
态 Schmidt 系数 E E E 乘积态 ∣ a ⟩ ⊗ ∣ b ⟩ \ket{a}\otimes\ket{b} ∣ a ⟩ ⊗ ∣ b ⟩ 1 , 0 1,\ 0 1 , 0 0 0 0 cos χ ∣ 00 ⟩ + sin χ ∣ 11 ⟩ \cos\chi\ket{00}+\sin\chi\ket{11} cos χ ∣ 00 ⟩ + sin χ ∣ 11 ⟩ cos 2 χ , sin 2 χ \cos^2\chi,\ \sin^2\chi cos 2 χ , sin 2 χ 0 < E < 1 0<E<1 0 < E < 1 Bell 态 1 2 , 1 2 \tfrac12,\ \tfrac12 2 1 , 2 1 1 1 1 (最大)
纠缠熵为 1 的态叫最大纠缠态 ,说「消耗 1 个 ebit」就是指用掉一对这样的态——
9.5、9.6 节的协议记账时全用这个单位。
◑ 物理图像为什么熵是「对」的量尺? 因为它不是随手挑的:可以证明,
两地只靠本地操作加打电话(LOCC)造不出新纠缠、也提不高纠缠熵 ——
它像能量一样是笔守恒的资源账。
而且它可操作:n n n 对纠缠熵为 E E E 的态,可以(渐近地)蒸馏出约 n E nE n E 对 Bell 态,
反过来也能用 n E nE n E 对 Bell 态稀释制备它们。
「熵 = 可兑换的 Bell 对数目」,这才配叫度量。
∑ 数学形式对 ∣ ψ 2 ⟩ = 0.99 ∣ 00 ⟩ + 0.141 ∣ 11 ⟩ \ket{\psi_2}=0.99\ket{00}+0.141\ket{11} ∣ ψ 2 ⟩ = 0.99 ∣ 00 ⟩ + 0.141 ∣ 11 ⟩ :
λ 1 = 0.9801 , λ 2 = 0.0199 (9.3.10) \lambda_1=0.9801,\quad\lambda_2=0.0199\tag{9.3.10} λ 1 = 0.9801 , λ 2 = 0.0199 ( 9.3.10 ) E = − 0.9801 log 2 0.9801 − 0.0199 log 2 0.0199 ≈ 0.141 (9.3.11) E=-0.9801\log_2 0.9801-0.0199\log_2 0.0199\approx 0.141\tag{9.3.11} E = − 0.9801 log 2 0.9801 − 0.0199 log 2 0.0199 ≈ 0.141 ( 9.3.11 ) 不到 Bell 态的 15%——「淡纠缠」的直觉有了数字。
要蒸馏出 1 对 Bell 态,平均得消耗约 7 对这样的态。
快捷公式:concurrence
熵要先求偏迹再对角化,对两比特有条捷径。开头那个判据行列式
α δ − β γ \alpha\delta-\beta\gamma α δ − β γ 正好给出 concurrence (一致度):
C = 2 ∣ α δ − β γ ∣ , 0 ≤ C ≤ 1 (9.3.12) C=2\,|\alpha\delta-\beta\gamma|,\qquad 0\le C\le 1\tag{9.3.12} C = 2 ∣ α δ − β γ ∣ , 0 ≤ C ≤ 1 ( 9.3.12 )
C = 0 C=0 C = 0 可分离,C = 1 C=1 C = 1 最大纠缠。它与 Schmidt 系数的关系是
C = 2 λ 1 λ 2 C=2\sqrt{\lambda_1\lambda_2} C = 2 λ 1 λ 2 ,与熵单调对应:
E = h ( 1 + 1 − C 2 2 ) , h ( p ) ≡ − p log 2 p − ( 1 − p ) log 2 ( 1 − p ) (9.3.13) E=h\!\left(\frac{1+\sqrt{1-C^2}}{2}\right),\qquad
h(p)\equiv-p\log_2 p-(1-p)\log_2(1-p)\tag{9.3.13} E = h ( 2 1 + 1 − C 2 ) , h ( p ) ≡ − p log 2 p − ( 1 − p ) log 2 ( 1 − p ) ( 9.3.13 )
验证一下:∣ Φ + ⟩ \ket{\Phi^+} ∣ Φ + ⟩ 有 α = δ = 1 2 \alpha=\delta=\tfrac{1}{\sqrt2} α = δ = 2 1 、β = γ = 0 \beta=\gamma=0 β = γ = 0 ,
C = 2 ⋅ 1 2 = 1 C=2\cdot\tfrac12=1 C = 2 ⋅ 2 1 = 1 ;态 cos χ ∣ 00 ⟩ + sin χ ∣ 11 ⟩ \cos\chi\ket{00}+\sin\chi\ket{11} cos χ ∣ 00 ⟩ + sin χ ∣ 11 ⟩ 给 C = sin 2 χ C=\sin 2\chi C = sin 2 χ ——
χ \chi χ 从 0 转到 45°,纠缠从 0 平滑爬到 1。
! 两个常见误会
「纠缠 = 相关」?不止。 经典硬币对(同抛同看)也有完美相关,
但那对应可分离的混合态。纠缠的独特之处在于每个测量轴上都相关 :
Bell 态沿 z z z 测完全相关,沿 x x x 测还是完全相关——经典关联做不到双轴通吃,
这正是下一节 Bell 不等式要开采的矿脉。
混合态的纠缠没这么简单。 本节的熵和 concurrence 公式只对纯态 成立。
混合态先天带经典噪声,「混」不再等于「纠缠」,
需要 Wootters 公式等更重的工具(本书不展开)。
判断混合态是否纠缠在一般维度甚至是 NP 难问题。
要「看见」这些关联,可以去实验室模块 08:纠缠与 Bell 非定域性 ,
两个粒子各自读数完全随机、合起来却严格相关的景象就在那里。
∎ 本节关键公式可分离判据(纯态)
α δ − β γ = 0 ⇔ ∣ ψ ⟩ = ∣ a ⟩ ⊗ ∣ b ⟩ \alpha\delta-\beta\gamma=0\ \Leftrightarrow\ \ket{\psi}=\ket{a}\otimes\ket{b} α δ − β γ = 0 ⇔ ∣ ψ ⟩ = ∣ a ⟩ ⊗ ∣ b ⟩ 系数矩阵秩 1;不为零即纠缠
Schmidt 分解
∣ ψ ⟩ = ∑ k λ k ∣ u k ⟩ A ∣ v k ⟩ B \ket{\psi}=\sum_k\sqrt{\lambda_k}\ket{u_k}_A\ket{v_k}_B ∣ ψ ⟩ = k ∑ λ k ∣ u k ⟩ A ∣ v k ⟩ B 两侧 ρ 的本征谱相同,均为 {λ_k}
纠缠熵
E = S ( ρ A ) = − ∑ k λ k log 2 λ k E=S(\rho_A)=-\sum_k\lambda_k\log_2\lambda_k E = S ( ρ A ) = − k ∑ λ k log 2 λ k 乘积态 0,Bell 态 1(单位:ebit)
Concurrence
C = 2 ∣ α δ − β γ ∣ = 2 λ 1 λ 2 C=2|\alpha\delta-\beta\gamma|=2\sqrt{\lambda_1\lambda_2} C = 2∣ α δ − β γ ∣ = 2 λ 1 λ 2 两比特纯态专用捷径,与 E 单调对应
1. 态 (∣00⟩ + ∣01⟩ + ∣10⟩ + ∣11⟩)/2 的纠缠熵是多少?
提交 2. 对 Bell 态 ∣Φ⁺⟩ 只测 A 粒子(沿任意轴),能观察到什么?
提交
接下来
我们已经能给纠缠称重,但还欠一个更根本的交代。
爱因斯坦当年的质疑并没有被「熵等于 1」回答:
也许两个粒子在分开前就偷偷商量好了所有答案——所谓相关,
不过是各自揣着同一张小抄?这个「局域隐变量」假说听起来无懈可击,
却能被推上实验的审判席。下一节:Bell 不等式,
以及量子力学如何在 2015 年拿到无漏洞的终审判决。