5935 words
30 minutes
毕业论文|Chapter 2 — Bell 不等式与 CHSH 博弈

量子纠缠的验证#

Bell不等式#

要验证量子纠缠态的存在性,我们首先要证明隐变量理论的矛盾性,这样就证明了非局域性是量子力学的本质特征。前面EPR论文提出了量子力学的不完备性,而应该由额外的变量来补充的论据。这些变量试图去恢复理论中的局域性(一个系统上的测量结果不受过去与之相互作用且遥远系统的影响)和因果性。在之后的一段时间里,有许多试图去完善隐变量理论的工作,但是均已失败告终。

Bell通过将隐变量理论数学化并加入局域性假设,证明了其与量子力学的统计预测不相容(具有矛盾),即非局域性是量子力学的典型特征。

斯特恩-格拉赫实验#

(斯特恩-格拉赫实验可以展开说)

考虑一对自旋为12\frac{1}{2}的粒子,整个系统处于自旋单态(总自旋为0),两个粒子自由地朝反方向运动。自旋单态的表达式为:

ψ=12()\ket{\psi} = \frac{1}{\sqrt{2}}(\ket{\uparrow \downarrow} - \ket{\downarrow \uparrow})

式中\ket{\uparrow\downarrow}表示两粒子系统的状态,也可写成\ket{\uparrow} \otimes \ket{\downarrow}(左边为粒子1的自旋状态),\ket{\uparrow}表示单个粒子自旋向上,\ket{\downarrow}表示自旋向下。

定义粒子1和2的自旋算符σ1\vec{\sigma_1}σ2\vec{\sigma_2}(Pauli算符形式,自旋算符与Pauli算符的关系为S=2σ\vec{S} = \frac{\hbar}{2}\vec{\sigma})。通过斯特恩-格拉赫实验测量某个方向(例如a\vec{a}为粒子1的测量方向)的自旋分量。

测量粒子1的σ1a\vec{\sigma_1} \cdot \vec{a}可能得到+1+11-1。由于系统反对称性(可补充解释),此时测量粒子2的σ2a\vec{\sigma_2} \cdot \vec{a}结果必然与粒子1相反(1-1+1+1)。


隐变量理论的数学化#

Bell通过以下两个假设引出隐变量原理:

  1. 局域性假设:若两粒子自旋的测量在空间分离的方向进行,则一个磁铁的方向(测量粒子1的方向)不会影响另一个磁铁的测量结果。
  2. 因果性假设:通过测量σ1\sigma_1的某个方向分量,可预测对应方向粒子2的自旋分量结果。

假设1指定了隐变量理论的局域性;假设2说明测量结果可被提前预测,对应隐变量λ\lambda的存在性(λ\lambda可以是单变量、一组变量或函数)。测量结果满足:

A(a,λ)=±1,B(b,λ)=±1A(\vec{a}, \lambda) = \pm 1, \quad B(\vec{b}, \lambda) = \pm 1

隐变量理论的期望值为:

P(a,b)=dλρ(λ)A(a,λ)B(b,λ)P(\vec{a}, \vec{b}) = \int d\lambda \, \rho(\lambda) A(\vec{a}, \lambda) B(\vec{b}, \lambda)

而量子力学对自旋单态的期望值为:

σ1a  σ2b=ab\langle \vec{\sigma}_1 \cdot \vec{a} \; \vec{\sigma}_2 \cdot \vec{b} \rangle = -\vec{a} \cdot \vec{b}

主要证明#

Bell不等式的推导步骤如下:

  1. 归一化条件dλρ(λ)=1\int d\lambda \, \rho(\lambda) = 1
  2. 测量结果限制A(a,λ),B(b,λ)=±1A(\vec{a}, \lambda), B(\vec{b}, \lambda) = \pm 1
  3. 引入新方向c\vec{c},计算差值:
P(a,b)P(a,c)=dλρ(λ)A(a,λ)A(b,λ)[A(b,λ)A(c,λ)1]P(\vec{a}, \vec{b}) - P(\vec{a}, \vec{c}) = \int d\lambda \, \rho(\lambda) A(\vec{a}, \lambda) A(\vec{b}, \lambda) \left[ A(\vec{b}, \lambda) A(\vec{c}, \lambda) - 1 \right]
  1. 取绝对值并利用积分性质,最终得到Bell不等式:
P(a,b)P(a,c)1+P(b,c)|P(\vec{a}, \vec{b}) - P(\vec{a}, \vec{c})| \leq 1 + P(\vec{b}, \vec{c})

验证实例#

a=(0,0,1)\vec{a} = (0,0,1), b=(1,0,0)\vec{b} = (1,0,0), c=(1/2,0,1/2)\vec{c} = (1/\sqrt{2}, 0, 1/\sqrt{2})

  • P(a,b)=0P(\vec{a}, \vec{b}) = 0
  • P(a,c)=1/2P(\vec{a}, \vec{c}) = -1/\sqrt{2}
  • P(a,b)P(a,c)=1/20.707|P(\vec{a}, \vec{b}) - P(\vec{a}, \vec{c})| = 1/\sqrt{2} \approx 0.707
  • 1+P(b,c)=11/20.2931 + P(\vec{b}, \vec{c}) = 1 - 1/\sqrt{2} \approx 0.293

显然0.707>0.2930.707 > 0.293,违反Bell不等式,说明隐变量理论无法解释量子力学预测。

f(g(x))={0,g(x)<21,g(x)2f(g(x)) = \begin{cases} 0, & g(x) < 2 \\ 1, & g(x) \geq 2 \end{cases}

结论#

通过证明Bell不等式在量子力学中不成立,表明非局域性是量子力学的本质特征,从而验证了量子纠缠的合理性。

CHSH游戏:#

上面描述了两种可行的实验方案,但是其实还有一种更容易理解的思想实验(CHSH游戏),通过这个游戏的分析表明任何经典的局部隐变量理论都不能够解释量子纠缠的情况。由于该游戏确实在物理上可实现,因此就证明了经典物理从根本上无法解释某些量子现象,至少在“局部”层面上无法解释。

CHSH游戏的定义:#

游戏有三个人,一个Alice,一个Bob,一个裁判:

  1. 首先,裁判以均匀的概率随机地选定两个数 xy{0,1}x,y\in \{0,1\}
  2. 然后,裁判把 xx给Alice,yy给Bob
  3. 最后,Alice需要给出回应a{0,1}需要给出回应a \in \{0,1\},Bob也需要给出回应b{0,1}b\in\{0,1\}

如果x=y=1x = y = 1,则只有当Alice和Bob要给出不一样的回应才可以获胜;对于其他的情况,则需要Alice和Bob的回应一样才能获胜。

经典情况:#

在经典情况下,如果我们去找一个对于这个游戏的最优策略,那么只需要找到胜率最高的那种策略,根据下面的分析:

xxyy最佳的
00ab相同a和b相同
01ab相同a和b相同
10ab相同a和b相同
11ab不同a和b不同

我们可以从上面的表格看到,有三种情况都需要Alice和Bob选择一样的回应,所以对于经典的最优解就是:Alice和Bob两者的回应保持相同,则获胜的概率最大为 34\frac{3}{4}(对于这种策略,四种情况只有一种情况输)。

量子情况:#

对于量子情况,首先介绍一些对于量子比特和量子态的基础知识。

基态:#

s1sn{0,1}ns_1\dots s_n \in \{0,1\}^n是一个长度为n的二进制串,i{0,1,,2n1}i\in \{0,1,\dots,2^{n}-1\}是对应二进制串的值,比如n=2,那么此时二进制串为s1s2={00,01,10,11}s_1s_2=\{00,01,10,11\},对应这四个二进制串的值为0,1,2,3。然后根据这个我们可以定义 s1sn\ket{s_1\dots s_n}是一个长度为 2n2^n的向量,假设这个向量里面的二进制串所对应的值为 ii,则这个向量的第 i+1i+1个位置为1,其余的位置为0,这个向量s1sn\ket{s_1\dots s_n}我们就称之为n个量子比特的基态。对于n=2(两个量子比特)的情况,那么就有 00,01,10,11\ket{00},\ket{01},\ket{10,\ket{11}}这些基态,这些基态分别对应(下面的基态都是列向量):

00=[1,0,0,0]T,01=[0,1,0,0]T10=[0,0,1,0]T,11=[0,0,0,1]T\begin{aligned} \ket{00} &= [1,0,0,0]^T,\quad \ket{01} = [0,1,0,0]^T \\ \ket{10} &= [0,0,1,0]^T,\quad \ket{11} = [0,0,0,1]^T \end{aligned}
量子态:#

定义:n个量子比特的量子态是一个向量 x\vec{x}:

x=x0000++x1111=S{0,1}nxsS\begin{aligned} \vec{x} &= x_{0\dots0}\ket{0\dots0}+ \dots + x_{1\dots1}\ket{1\dots1} \\ &= \sum_{S\in\{0,1\}^n} x_s \ket{S} \end{aligned}

上式中 SS是遍历所有长度是n的二进制串,并且 x\vec{x}满足归一化条件:

x002++x112=1.|x_{0\ldots0}|^2 + \cdots + |x_{1\ldots1}|^2 = 1.

系数 x00,,x11x_{0\dots0},\dots,x_{1\dots1}可以是复数,但是我们只考虑它们是实数的情况。

量子门:#

就像函数 f(x)f(x)一样,有输出对应着输入,这个函数相当于对这个输入做了某种操作;那么对应某个系统的量子态,我们应该如何去操作它改变它呢?上面我们提到了量子态是一个向量 x\vec{x},那么量子门操作就是 f(x)f(\vec{x}),并且要遵循如下的性质:

  1. 归一化:对于所有满足归一化条件的 x\vec{x},在经过量子门操作之后的量子态 y=f(x)\vec{y}=f(\vec{x})也要满足归一化关系。

  2. 线性操作(可以在附录说一下为什么线性是必要的):对于所有的可能的输入的向量 ab\vec{a},\vec{b}以及所有实数 cc,我们有 f(a+b)=f(a)+f(b)f(\vec{a}+\vec{b})=f(\vec{a})+f({\vec{b}})f(ca)=cf(a)f(c\vec{a})=cf(\vec{a}).

    对于满足上面两个性质的量子门操作,我们称之为线性量子门。

量子系统,孤立系统:#

对于要如何去实现这样的二进制的量子态我们可以直接想到的是电子,我们可以依靠电子的自旋量子数是1/2的特性(这个特性使得电子再任意方向上的自旋都有两个本征态,自旋向上,自旋向下)去实现基态是 0,1\ket{0},\ket{1}的量子系统。顺着这个思路,我们可以得出量子系统其实就是多个有顺序的电子的集合表示为 E=(e1,,en)E =(e_1,\dots,e_n)

因为每个电子要么自旋向上(0)要么自旋向上(1),所以这个量子系统的状态一定是二进制串 S=s1sn{0,1}nS = s_1\dots s_n\in\{0,1\}^n中的一个,那么我们可以将测量这个量子系统中每一个电子的自旋的过程描述为 ESE \to S

在大部分情况下,我们只关心庞大世界中的某一个系统,那么此时这个孤立系统是怎么样的,和复合系统的区别是什么?一个孤立的量子系统 EE之所以孤立是因为他的测量事件和其他量子系统的测量是独立的,即:

这个孤立系统的测量事件:ESE \to S与所有其他的系统的测量 ESE’ \to S'是独立的,也就是这个系统的测量结果并不会影响其他系统的结果。

上面我们提到了测量事件和测量结果,那么对于一个处于某一个量子态的量子系统E,这些是如何定义的呢?

测量公理:#

假设 EE是一个由n个电子 (e1,,en)(e_1,\dots,e_n) 组成的孤立量子系统,根据前面的定义,我们可以用一个向量 x\vec{x}来表示这个系统的状态,如果我们去测量这个系统的状态,可能的结果是所有二进制串(长度为n) S{0,1}nS\in\{0,1\}^n中的一个,或者说这个系统的状态是不同基态的组合:

ψ=S{0,1}nxsS\ket{\psi} = \sum_{S\in\{0,1\}^n}x_s\ket{S}

那么此时我们知道当我们观测这个系统的时候,这个系统一定会处于基态中的一个,但是在观测前我们也许可以知道我们测量得到不同基态的概率 xs2|x_s|^2,测量公理定义了这样的概率是:

Pr[ES]=xs2Pr [E \to S] =|x_s|^2

以及将我们可能对这个处于 x\vec{x} 孤立系统E做出的物理操作定义为一个量子门 ff,这样就使得前面有关与量子门的操作和定义能够在实际物理世界中存在对应。

总而言之,这个公理将一个孤立的系统描述为某一个量子态,以及每一个物理上的操作(Physical process)描述为某一个量子门。

张量积(Tensor Product):#

回想一下我们CHSH游戏中指定了Alice和Bob共享一个EPR对 ψ=12(00+11)\ket{\psi}=\frac{1}{\sqrt{2}}(\ket{00}+\ket{11}),那么这个状态其实是Alice的电子和Bob的电子相组合成的系统(但是这个两个系统之间有纠缠)。前面我们知道对于一个系统的状态可以用一个向量 x=(x1,,xn)T\vec{x}=(x_1,\dots,x_n)^T(长度为n)来描述,如果此时还有另一个系统状态是 y=(y1,,ym)T\vec{y}=(y_1,\dots,y_m)^T(长度为m),那么两个系统组合起来的状态由向量 Φ\vec{\Phi} 表示,那么这个向量应该如何定义呢?

向量的张量积:#

我们这里引入向量之间的张量积,定义为:

Φ=xy\vec{\Phi} = \vec{x}\otimes\vec{y}

并且这个向量 Φ\vec{\Phi}的长度是nm,并且 \otimes被称之为张量积。向量x\vec{x}和向量y\vec{y}之间的张量积定义为:

xy=(a1b1,a1b2,,anbn1,anbn\vec{x}\otimes\vec{y} = (a_1b_1,a_1b_2,\dots,a_nb_{n-1},a_nb_n

比如,但两个系统(E1E2E_1,E_2)的状态分别为 a=(a1,a2)\vec{a}=(a_1,a_2)b=(b1,b2,b3)\vec{b}=(b_1,b_2,b_3),那么这两个系统组合的状态 ψ\vec{\psi}为:

ψ=ab=[a1(b1,b2,b3),a2(b1,b2,b3)]=(a1b1,a1b2,a1b3,a2b1,a2b2,a2b3)\begin{aligned} \vec{\psi} = \vec{a}\otimes \vec{b} &= [a_1(b_1,b_2,b_3),\,a_2(b_1,b_2,b_3)] \\ &= (a_1b_1,a_1b_2,a_1b_3,a_2b_1,a_2b_2,a_2b_3) \end{aligned}
不同孤立系统状态的张量积:#

上面对于两个系统组合起来的张量积表示有一个前提条件就是,它们两个都是孤立系统,即它们不能有相互作用或者纠缠(如EPR态)。它表述为一个定理:

如果两个系统 E1E_1E2E_2都是孤立的系统,并且两个分别由向量 x\vec{x}y\vec{y}表示,则它们的联合系统状态 E1E2E_1E_2 用向量 xy\vec{x}\otimes\vec{y}表示

作用在两个系统的量子门的张量积:#

我们已经定义了两个孤立系统的联合状态是怎么样的,那么我们想,如果我们想对系统(A和B)进行物理操作(相当于作用量子门),那么这些操作是如何影响整个系统的联合状态的呢?我们可以定义:

fVnVnf:V^n\to V^n表示对于系统A的一个线性操作(线性的量子门),g:VmVmg:V^m\to V^m表示对于系统B的一个线性操作,则对于这整个系统的联合状态 ab\vec{a}\otimes\vec{b}而言,我们对于整个系统的操作 fgf\otimes g作用效果为:

(fg)(ab)=f(a)g(b)(f\otimes g)(\vec{a}\otimes\vec{b}) = f(\vec{a})\otimes g(\vec{b})

这里对于整个系统AB的操作 fgf\otimes g 可以由一个公理指定(ff作用与A系统,gg作用于B系统)。(此处可以说明一下这里如果两个操作并不同时发生,这个定义还是有效的吗?,因为两个系统是孤立系统,所以互不影响。

CHSH游戏的一种解法:#

在CHSH游戏中,我们让Alice和Bob共享一个Bell态 12(00+11)\frac{1}{\sqrt{2}}(\ket{00}+\ket{11}),将第一个电子分配给Alice(左边的比特),将第二个电子分配给Bob(右边的比特),那么此时Alice和Bob都可以测量各自电子的自旋,并且得到的结果要么是+1要么是-1,最后将测量的结果回应给裁判(但是这里的结果是+1,-1,我们之后会处理将这个测量结果和0,1对应起来)但是此时两者是严格正相关的,也就是当Alice测量的+1的时候,整个系统的状态就塌缩到 00\ket{00},那么这时候Bob测量的结果就一定是+1,另一种情况也是一样。

这种情况下恰好就对应了经典情况下的最优策略(始终保持Alice和Bob的回应相同),这也就说明了为什么在量子情况要用Bell态,因为在Alice和Bob不对Bell态做出任何操作的情况下,就已经处于经典情况胜率的上界(75%)了。

所以如果要找到比经典情况最优策略胜率更高的策略,那么Alice和Bob就必须对这个Bell态进行物理操作了。

旋转门的引入:#

因为我们知道,Alice和Bob都是将各自电子自旋的测量结果(进行处理),最后返回给裁判是0和1,但是其实对于电子的自旋测量可以改变方向,那么就需要引入一个量子门操作 RθR_{\theta},它的几何意义其实就是将二维空间中的某一个向量顺时针旋转 θ\theta角度,由此可以得到他的数学定义是1

Rθ(a10+a21)=(a1cosθa2sinθ)0+(a1sinθ+a2cosθ)1\begin{aligned} R_{\theta}(a_1\ket{0}+a_2\ket{1}) &= (a_1\cos{\theta}-a_2\sin{\theta})\ket{0} \\ &\quad + (a_1\sin{\theta}+a_2\cos{\theta})\ket{1} \end{aligned}

有了旋转门之后,Alice和Bob就可以分别对于各自的电子进行旋转,之后再执行测量操作,我们可以将Alice和Bob的操作记为 RA0,1RB0,1R_{A_{0,1}}\otimes R_{B_{0,1}},这里的角标对应着,当Alice或者Bob收到x或者y为0和1的情况。

胜率的定义:#

我们这里将Alice和Bob的总胜率定义为 ww。在CHSH游戏中,Alice和Bob分别收到来自裁判的x,y,Alice和Bob分别给出回应a,b。对于每一对可能的(x,y)(x,y),他们获胜的条件可以简化为:

ab=xya \oplus b = x \cdot y

前面是逻辑上的异或操作(简而言之就是相同为0,不同为1),后面是普通的乘法操作(只有当 x=y=1x=y=1时,右边才等于1)2

假设对于每一种 (x,y)(x,y)的胜率记作 Pr[win(x,y)]Pr[win|(x,y)],由于每种 (x,y)(x,y)出现的概率相同都是1/4,所以总的胜率 ww为:

w=14(x,y)Pr[win(x,y)]w=\frac{1}{4}\sum_{(x,y)}Pr[win|(x,y)]
测量结果的处理:#

我们知道Alice和Bob对于任意方向电子自旋的测量结果是+1/-1,但是游戏规则中需要Alice和Bob做出回应0/1,那么我们就需要做出如下的处理:

a=1a2,b=1b2a = \frac{1 - a^\ast}{2}, \quad b = \frac{1 - b^\ast}{2}

这里的 Alice和Bob的测量结果a=±1,b=±1a^\ast=\pm1,b^\ast=\pm1,而 aabb就是Alice和Bob对于测量结果处理之后给裁判各自的回应,回应0就对应着测量结果1,回应1就对应测量结果是-1。

因此上面对于a,b的获胜条件,转化为测量结果aa^\astbb^\ast的获胜条件为:

ab=xy(1a2)(1b2)=xy\begin{aligned} a \oplus b &= x \cdot y \\ \left(\frac{1-a^\ast}{2}\right)\oplus\left(\frac{1-b^\ast }{2}\right) &= x \cdot y \end{aligned}

代入之后第二个式子可以用一个公式化简:

uv=u+v2uvu \oplus v = u+v-2uv

上面的 uvu和v都代表的是二进制数。所以得到:

1ab2=xyab=12xy\begin{aligned} \frac{1-a^\ast b^\ast}{2} &= x \cdot y \\ a^\ast b^\ast &= 1-2xy \end{aligned}

我们可以用一个表格看到12xy1-2xy其实可以等价为(1)xy(-1)^{xy}

xxyyxyxy12xy1-2xy(1)xy(-1)^{xy}
00011
01011
10011
111-1-1

由此可以得到获胜条件变成Alice和Bob的测量结果要满足:

ab=(1)xya^\ast b^\ast =(-1)^{xy}
用测量结果的乘积表示赢的概率:#

记任意一次测量 (a,b)(a^\ast,b^\ast)的概率分布为{P++,P+,P+,P}\{P_{++},P_{+-},P_{-+},P_{--}\}为对应结果的概率,角标的左边表示Alice的测量结果,右边表示Bob的测量结果,正负号代表测量结果为1和-1。所以双方测量结果乘积的期望是:

E[ab]=(+1)P+++(1)P++(1)P++(1)PE[a^\ast b^\ast]=(+1)P_{++}+(-1)P_{+-}+(-1)P_{-+}+(1)P_{--}

因为赢得条件是ab=(1)xya^\ast b^\ast =(-1)^{xy},所以对于每一种情况的(x,y)(x,y)而言,获胜的概率Pr[win(x,y)]Pr[win|(x,y)]就等价于 Pr[ab=(1)xy]Pr[a^\ast b^\ast=(-1)^{xy}]

对于任意一次测量(a,b)(a^\ast,b^\ast)的概率分布{P++,P+,P+,P}\{P_{++},P_{+-},P_{-+},P_{--}\},有:

Pr[ab=(1)xy]=ab=(1)xyPr[a,b]Pr[a^\ast b^\ast=(-1)^{xy}] = \sum_{a^\ast b^\ast=(-1)^{xy}}Pr[a^\ast,b^\ast]

式中Pr[a,b]Pr[a^\ast,b^\ast]代表两者测量的概率分布(离散的),所以对于等式的右侧,我们可以分两种情况:

ab=(1)xyPr[a,b]={P+++P,(1)xy=+1P++P+,(1)xy=1\sum_{a^\ast b^\ast=(-1)^{xy}} \mathrm{Pr}[a^\ast,b^\ast] = \begin{cases} P_{++}+P_{--}, & (-1)^{xy}=+1 \\ P_{+-}+P_{-+}, & (-1)^{xy}=-1 \end{cases}

可以看到:

P+++P=12[1+(P+++PP+P+)]=12(1+E[ab])\begin{aligned} P_{++}+P_{--} &= \frac{1}{2}\left[1+(P_{++}+P_{--}-P_{+-}-P_{-+})\right] \\ &= \frac{1}{2}(1+E[a^\ast b^\ast]) \end{aligned}

第二个等号是由归一化 P+++P++P++P=1P_{++}+P_{+-}+P_{-+}+P_{--}=1。同样的,对于P++P+P_{-+}+P_{+-}有:

P++P+=12(1E[ab])P_{-+}+P_{-+}=\frac{1}{2}(1-E[a^\ast b^\ast])

所以可以得到:

Pr[winx,y]=1+(1)xyE[ab]2Pr[win|x,y]=\frac{1+(-1)^{xy}E[a^\ast b^\ast]}{2}

在Alice和Bob收到(x,y)(x,y)后,他们对应的做出对于各自电子自旋测量方向的旋转操作AxA_xByB_y,对于他们共享的Bell态ψ\ket{\psi}而言就是量子门操作AxByA_x \otimes B_y,那么他们测量结果的乘积的期望E[ab]E[a^\ast b^\ast]就是:

E[ab]=ψAxByψE[a^\ast b^\ast] = \bra{\psi}A_x \otimes B_y\ket{\psi}

然后将这个期望值代入式23,再将式23代入总的获胜概率ww的定义,得到:

ω=14x,y1+(1)xyψAxByψ2=12+18x,y(1)xyψAxByψ.\begin{aligned} \omega&=\frac{1}{4}\sum_{x,y}\frac{1+(-1)^{xy}\,\langle\psi|A_x\otimes B_y|\psi\rangle}{2} \\ &=\frac{1}{2}+\frac{1}{8}\sum_{x,y}(-1)^{xy}\,\langle\psi|A_x\otimes B_y|\psi\rangle. \end{aligned}

我们可以定义总获胜概率为ww和失败概率为LL,以及两者之间的差值为D=wLD=w-L3,然后

D=w(1w)=2w1D = w-(1-w) = 2w-1

所以要得到总的获胜概率的上界等价于找到D的上界,由上可以得到:

D=2ω1=2[12+18x,y(1)xyψAxByψ]1=14x,y{0,1}(1)xyψAxByψ.\begin{aligned} D&=2\omega-1=2\left[\frac{1}{2}+\frac{1}{8}\sum_{x,y}(-1)^{xy}\langle\psi|A_x\otimes B_y|\psi\rangle\right]-1 \\ &=\frac{1}{4}\sum_{x,y\in\{0,1\}}(-1)^{xy}\,\langle\psi|A_x\otimes B_y|\psi\rangle. \end{aligned}

那么现在应该如何找到这个D的上界呢?

如果根据量子力学的期望的定义,D可以写成:

D=14(A0B0+A0B1+A1B0A1B1)D=\frac{1}{4}\left(\langle A_0B_0\rangle+\langle A_0B_1\rangle+\langle A_1B_0\rangle-\langle A_1B_1\rangle\right)

其中,A0B0\langle A_0B_0\rangle表示在量子态 ψ\ket{\psi}下算符 A0B0A_0 \otimes B_0的期望值,其余同理。

然后定义四个投影算符 P,P,Q,QP,P',Q,Q',我们可以看成是一个投影操作,那么这个算符的本征值只有两个0和14,我们作用到Alice和Bob各自电子上的算符仅仅改变了电子的自旋方向,但是所得到的测量值还是+1和-1,所以可以做出如下的变换:

A0=2P1B0=2Q1A1=2Q1B1=2P1\begin{aligned} A_0 &= 2P-1 \\ B_0 &= 2Q'-1 \\ A_1 &= 2Q-1 \\ B_1 &= 2P'-1 \end{aligned}

做出这个变换的理由是:我们对于电子的旋转的操作总可以等效成将电子的自旋向某个方向投影,并且保持了测量值的不变性。

然后定义R:

R=A0B0+A0B1+A1B0A1B1R =\langle A_0B_0\rangle+\langle A_0B_1\rangle+\langle A_1B_0\rangle-\langle A_1B_1\rangle

C=A0B0+A0B1+A1B0A1B1C =A_0B_0+A_0B_1+A_1B_0-A_1B_1 可以得到:

R=CR =\langle C\rangle

其中:

C=A0(B1+B0)+A1(B0B1)C = A_0(B_1+B_0)+A_1(B_0-B_1)

然后我们计算C2C^2,因为A0A1B0B1A_0,A_1,B_0,B_1的本征值都是 ±1\pm 1,所以有:

A02=A12=B02=B12=1A_0^2=A_1^2=B_0^2=B_1^2=1

5

由此可以推出:

C2=4+[A0,A1][B1,B0]C^2 = 4+[A_0,A_1][B_1,B_0]

其中[A,B][A,B]表示算符的对易子ABBAAB-BA,从这个等式可以看出,在经典情况下,Alice和Bob并不做出任何操作,所以等式右侧的第二项就是0,但是这一项正表现出了量子的非定域性特点。

然后我们将变换关系29,代入计算得到:

C2=4+16[P,Q][P,Q]C^2 =4+16[P,Q][P',Q']

由于P,Q都是投影算符,所以有:

P2=P,Q2=Q,P=P,Q=QP^2 = P ,\quad Q^2=Q ,\quad P^{\dagger} =P,\quad Q^{\dagger} = Q

因为A0=2P1A_0 = 2P-1,所以:

A02=(2P1)(2P1)=4P24P+I=4P24P2+I=I\begin{aligned} A_0^2 &= (2P-1)(2P-1) \\ &=4P^2-4P+I \\ &=4P^2-4P^2+I \\ &=I \end{aligned}

所以A0=1\|A_0\|=1,同理A12=IA_1^2 = IA1=1\|A_1\|=1

我们想得到[P,Q][P,Q]的范围,有A0A1=1\|A_0\|,\|A_1\|=1,以及:

[A0,A1]=(2PI)(2QI)(2QI)(2PI)=4(PQQP)=4[P,Q]\begin{aligned} [A_0,A_1] &= (2P-I)(2Q-I)-(2Q-I)(2P-I) \\ &=4(PQ-QP) \\ &=4[P,Q] \end{aligned}

所以我们有[P,Q]=14[A0,A1][P,Q] = \frac{1}{4}[A_0,A_1]。所以现在的想法是找一下[A0.A1][A_0.A_1]的范围。

考虑[A0,A1][A_0,A_1]的范数(大小),再根据三角不等式,可以得到:

[A0,A1]=A0A1A1A0A0A1+A1A02A0A1=2\begin{aligned} \|[A_0,A_1]\|&=\|A_0A_1-A_1A_0\| \\ &\leq \|A_0A_1\|+\|A_1A_0\| \\ &\leq 2\|A_0\|\|A_1\|=2 \end{aligned}

然后根据[P,Q]=14[A0,A1][P,Q] = \frac{1}{4}[A_0,A_1],得到:

[P.Q]12\|[P.Q]\| \leq \frac{1}{2}

对于[P,Q][P',Q']的范数也同理小于等于二分之一,所以对于式(35)而言有范数形式:

C2=4+16[P,Q][P,Q]\|C^2\| = 4+16\|[P,Q][P',Q']\|

又因为[P.Q][P,Q][P,Q][P,Q]\|[P.Q][P',Q']\|\leq \|[P,Q]\|\|[P',Q']\|(三角不等式),所以:

C24+16[P,Q][P,Q]4+16×148\begin{aligned} \|C^2\| &\leq 4+16\|[P,Q]\|\|[P',Q']\| \\ &\leq 4+16\times\frac{1}{4} \leq 8 \end{aligned}

第二个不等式代入了[P,Q]12\|[P,Q]\| \leq \frac{1}{2}[P,Q]12\|[P',Q']\| \leq \frac{1}{2}。所以得到了C22\|C\| \leq 2\sqrt{2} 这一结论。

因为R=CR = \langle C\rangle,所以对于R的绝对值有:

R=ψCψψCψ=Cψ|R| = |\bra{\psi}C\ket{\psi}| \leq \|\psi\|\|C\psi\| = \|C\psi\|

上式中不等号是由于施瓦茨不等式uvuv|uv|\leq \|u\|\|v\|,最后的等号是因为ψ\ket{\psi}是单位向量(所以它的模是1)。

然后通过算符范数的定义:C\|C\|是拉伸向量的最大的因子6。由这个定义我们可以得到:

C=supψ=1{Cψ}\|C\| = \sup_{\|\psi\|=1}\{\|C\psi\|\}

其中sup\sup是对于所有满足要求的ψ\psi找到目标函数的上界。于是对于式(19)有:

RCψsupψ=1{Cψ}=C|R| \leq \|C\psi\| \leq \sup_{\|\psi\|=1}\{\|C\psi\|\}=\|C\|

所以对于R2|R|^2有:

R2C2=C28|R|^2 \leq \|C\|^2 = \|C^2\| \leq 8

7

所以,对于R来说,有R22|R| \leq 2\sqrt{2}。然后根据定义D=14RD = \frac{1}{4}R,得到总胜率和总败率之间的偏差D的上界是22\frac{\sqrt{2}}{2},即:

D=14R22D = \frac{1}{4}R \leq \frac{\sqrt{2}}{2}

最后我们通过D=2w1D = 2w-1,得到量子情况下,总胜率ww的上界为:

w=D+1224+1285%w = \frac{D+1}{2} \leq \frac{\sqrt{2}}{4}+\frac{1}{2} \approx 85\%

所以可以通过CHSH游戏8表明,通过两个距离很远(足以排除以光速进行信息传递)的双方,各自可以访问(可以测量并且操作)一个量子纠缠态的一半,通过合适的测量操作,可以给出比经典最优策略更优的量子策略。因此对于这个游戏的分析表明没有经典的隐变量理论(局域性理论)可以解释量子纠缠所带来的一系列的相关性(一方的测量可以非局域性地影响另一方的测量结果)和后果(比经典更优的策略)。

附录#

相关概念#

  • [[Bell Inequality]] — 证明任何局域隐变量理论与量子力学不相容的数学不等式
  • [[CHSH Game]] — Bell 定理的博弈论表述:经典胜率 75% vs 量子胜率 85%
  • [[Quantum Entanglement]] — 量子非局域关联,Bell 不等式的物理基础
  • [[Hidden Variable Theory]] — 被 Bell 定理证伪的局域隐变量尝试

前情见 Chapter 1 — [[EPR Paradox]] 的起源。应用见 Chapter 3 — [[Shor Algorithm]]。

参考文献#

版权信息#

本文原载于 quantum51.top,遵循 CC BY-NC-SA 4.0 协议,复制请保留原文出处。

Footnotes#

  1. 原文注释:后面可以在附录添加这个的几何推导

  2. 原文注释:CHSH游戏的规则是:当 x=y=1x=y=1时,Alice和Bob的回应 aba和b要不一样(ab=1a\oplus b=1)才能获胜,剩下的情况需要 ab要一样a和b要一样ab=0a\oplus b=0)才能获胜

  3. 原文注释:这里总是认为我们的策略总是让Alice和Bob的获胜概率大于他们失败的概率,因为就算是不改变Bell态,获胜的概率都是75%

  4. 原文注释:why?

  5. 原文注释:解释放在附录里面

  6. 原文注释:详细定义见附录

  7. 原文注释:上式中C2=C2\|C\|^2 = \|C^2\|对于一般的normal算符是成立的(具体证明在附录)

  8. 原文注释:CHSH游戏是可以在实验上实现的

毕业论文|Chapter 2 — Bell 不等式与 CHSH 博弈
/posts/thesis-chapter-2/
Author
俞小金 / Eugene
Published at
2025-04-16

Some information may be outdated