量子计算系列:
唠唠闲话
量子系列第二篇,内容如下:
本篇介绍 Bloch 球及其性质,在 Nakahara 第三章的基础上,补充命题和推导。
纯态与 Bloch 球面
量子比特
经典比特的取值为 1 或 0,由二极管的通电断电实现。
量子比特的取值是二维复空间 C2 上的单位向量,以 ∣0⟩=(10) 和 ∣1⟩=(01) 为一组基。
物理上,实现基元 ∣0⟩ 和 ∣1⟩ 有多种方案:
-
自旋 为 21 的粒子,比如电子或质子
- ∣0⟩=∣↑⟩ 为上旋
- ∣1⟩=∣↓⟩ 为下旋
-
偏振光子(polarized photon)
- ∣0⟩ 取竖直偏振的光子 ∣↕⟩
- ∣1⟩ 取水平偏振的光子 ∣↔⟩
注记:
- 量子比特用复向量定义,因而属于纯态;
- 自旋是粒子的内禀属性,自旋度为 S 的粒子能提供 d=2S+1 种正交的自旋状态,用于构造 d 种状态叠加的粒子,但这类粒子的重要性还有待探索。
- 量子计算与经典计算在硬件实现上有本质的不同,后者处理离散的0-1,而前者处理向量空间。
Bloch 球面
回顾上篇,纯态 ∣ψ⟩ 代表集合等价类
{eiα∣ψ⟩∣α∈R}
我们给出代表元的标准形式,用以建立纯态与球面的对应。引理: C2 中的纯态有如下标准形式
∣ψ(θ,ϕ)⟩=cos2θ∣0⟩+eiϕsin2θ∣1⟩,θ∈[0,π],ϕ∈[0,2π]
证明:
将代表元 a∣0⟩+b∣1⟩∈C2 标准化:
- 不妨设系数 a 为实数,否则调整幅角使 eiαa 为实数;将系数 b 拆成指数形式,代表元化
a∣0⟩+beiα∣1⟩, a,b,α∈R
- 由归一化条件 a2+b2=1,代表元化
cos2θ∣0⟩+eiϕsin2θ∣1⟩,θ∈[0,4π],ϕ∈[0,2π]
夹角取 2θ 方便与球面建立对应。
- 将 θ 取值范围缩小为 [0,2π],当 θ>2π 时,如下变量代换
cos2θ∣0⟩+eiϕsin2θ∣1⟩==whereθ′=θ−2π∈[0,2π]−(cos2θ−2π∣0⟩+eiϕsin2θ−2π∣1⟩)−(cos2θ′∣0⟩+eiϕsin2θ′∣1⟩)
- 将 θ 取值范围缩小为 [0,π],当 θ>π 时,如下变量代换
cos2θ+π∣0⟩+eiϕsin2θ+π∣1⟩==whereθ′=π−θ∈[0,π]−cos(2π−2θ)∣0⟩+eiϕsin(2π−2θ)∣1⟩−(cos2θ′∣0⟩+ei(ϕ+π)sin2θ′∣1⟩)
现在,对比二维球面 S2 与纯态标准形式
S2={n(θ,ϕ)=(sinθcosϕ,sinθsinϕ,cosθ)t:θ∈[0,π],ϕ∈[0,2π]}∣ψ(θ,ϕ)⟩=cos2θ∣0⟩+eiϕsin2θ∣1⟩,θ∈[0,π],ϕ∈[0,2π]
二者参数 θ 和 ϕ 的取值范围相同,且易证
n(θ1,ϕ1)=n(θ2,ϕ2)⇔∣ψ(θ1,ϕ1)⟩=∣ψ(θ2,ϕ2)⟩
于是 C2 纯态与二维球面 S2 借由 (θ,ϕ) 建立一一对应。
我们称 S2 为 Bloch 球面,称向量 n(θ,ϕ) 为量子态 ∣ψ(θ,ϕ)⟩ 的 Bloch 向量。
C2 纯态构成二维曲面的一种理解:
- 从实数域看, ∣ψ⟩∈C2 是四维空间 R4 上的向量,有 4 个自由变量;
- 由归一化条件 ⟨ψ∣ψ⟩=1,ψ 的自由变量剩 3 个;
- 再由整体相位 eiα,自由变量剩 2 个,得到二维球面 S2。
C2 是维数最低的量子系统,恰好能用 Bloch 球理解其量子态,但更高维的量子系统就很难再找到这种对应了。
几何性质
量子态 ∣ψ⟩ 是复空间 C2 中的向量,从实数域上看,C2 是个四维空间,我们很难想象。但借助 C2 纯态与二维球面 S2 的对应,我们可以用 Bloch 球面的几何性质理解 C2 中纯态的几何性质。
下边探讨二者几何性质的联系,并导出泡利矩阵。
坐标轴
如下建立直角坐标系
坐标轴的量子态与狄拉克记号::
狄拉克符号 |
坐标轴方向 |
S2 坐标 |
(θ,ϕ) |
量子态 |
∣0⟩ |
z+ |
(0,0,1) |
(0,0) |
∣0⟩=(10) |
∣1⟩ |
z− |
(0,0,−1) |
(π,0) |
∣1⟩=(01) |
∣+⟩ |
x+ |
(1,0,0) |
(2π,0) |
2∣0⟩+∣1⟩=21(11) |
∣−⟩ |
x− |
(−1,0,0) |
(2π,π) |
2∣0⟩−∣1⟩=21(1−1) |
∣+i⟩ |
y+ |
(0,1,0) |
(2π,2π) |
2∣0⟩+i∣1⟩=21(1i) |
∣−i⟩ |
y− |
(0,−1,0) |
(2π,23π) |
2∣0⟩−i∣1⟩=21(1−i) |
正交与反向
正交的量子态对应方向相反的 Bloch 向量:
⟨ψ(θ2,ϕ2)∣ψ(θ1,ϕ1)⟩=0⇔n(θ1,ϕ1)=−n(θ2,ϕ2)
证明:
令 α+=n(θ,ϕ),α−=n(π−θ,π+ϕ) 为方向相反的两个向量,计算內积得
⟨α+∣α−⟩=⟨ψ(θ,ϕ)∣ψ(π−θ,π+ϕ)⟩=(cos2θe−iϕsin2θ)(cos2π−θei(ϕ+π)sin2π−θ)=(cos2θe−iϕsin2θ)(sin2θ−eiϕcos2θ)=0
推论:量子态 ∣α+⟩ 和 ∣α−⟩ 作成 C2 的一组正交基。
球投影算子
对向量 α+=n(θ,ϕ)∈S2 ,定义球投影算子
σα+=∣α+⟩⟨α+∣−∣α−⟩⟨α−∣=(cos2θeiϕsin2θ)(cos2θe−iϕsin2θ)−(sin2θ−eiϕcos2θ)(sin2θ−e−iϕcos2θ)=(cosθeiϕsinθe−iϕsinθ−cosθ)
由定义知,∣α±⟩ 分别为 σα+ 的特征值 ±1 的特征向量。测量上,算子 σα+ 对纯态 ∣ψ′⟩ 的测量期望为 ∣ψ′⟩ 的 Bloch 向量 n′ 在 α+=n 方向的投影,即:
⟨σα+⟩=⟨ψ(θ′,ϕ′)∣σα+∣ψ(θ′,ϕ′)⟩=(cos2θ′e−iϕ′sin2θ′)(cosθeiϕsinθe−iϕsinθ−cosθ)(cos2θ′eiϕ′sin2θ′)=cosθcosθ′+sinθsinϕsinθ′sinϕ′+sinθcosϕsinθ′cosϕ′=n(θ′,ϕ′)⋅n(θ,ϕ)
Ps:为不引起歧义,我们称 σα+ 为球投影算子,相对 Bloch 球的几何而言,注意与量子态的投影算子 Pα+ 区分:
σα+Pα+=∣α+⟩⟨α+∣−∣α−⟩⟨α−∣=∣α+⟩⟨α+∣
泡利矩阵
特别地,将坐标轴正方向的球投影算子 σx,σy,σz 称为泡利矩阵,由定义得,泡利矩阵 σz
σz=∣0⟩⟨0∣−∣1⟩⟨1∣=(10)(10)−(01)(01)=(100−1)
同理可求得 σx 和 σy 为
σx=(0110), σy=(0i−i0)
纯态 ∣ψ⟩ 的 Bloch 向量可由泡利矩阵量得到
n=(⟨σx⟩,⟨σy⟩,⟨σz⟩)
其中 ⟨σx⟩,⟨σy⟩,⟨σz⟩ 为 ∣ψ⟩ 的测量期望。
直觉上,泡利矩阵提供物理量子系统 C2 的坐标系。
实际上,记 σ=(σx,σy,σz) ,则有
σα+=(cosθeiϕsinθe−iϕsinθ−cosθ)=(cosθ(cosϕ+isinϕ)sinθ(cosϕ−isinϕ)sinθ−cosθ)=sinθcosϕ⋅σx+sinθsinϕ⋅σy+cosθ⋅σz=(sinθcosϕ,sinθsinϕ,cosθ)⋅(σx,σy,σz)=n(θ,ϕ)⋅σ
向量 n 与 σ 点积,得到 n 方向的球投影算子 σα+ 。
混合态与 Bloch 球
C2 的量子态与 Bloch 球的向量一一对应,其中纯态对应 Bloch 球面,混合态对应 Bloch 球内点,下边通过“坐标分解”建立这一对应。
坐标分解
回顾上篇:
- 混合态用密度矩阵 ρ 表示
- H 上的密度矩阵全体记为 S(H)
- S(H) 与量子态一一对应,其中秩为 1 的部分对应纯态
- 密度矩阵 ρ 由下边三条性质等价刻画:
- ρ†=ρ
- ρ 为半正定矩阵
- Tr(ρ)=1
借助泡利矩阵,我们给出密度矩阵 ρ∈S(C2) 的坐标分解:
- 记 I=(1001) 为单位矩阵,则 {I,σx,σy,σz} 作成 EndC(H) 的一组基元,设:
ρ=bxσx+byσy+bzσz+b0I,bx,by,bz,b0∈C
- 断言 bx,by,bz,b0 为实数,由于
⇒⇒⇒ρ†=ρ bxσx+byσy+bzσz+b0I=bx∗σx+by∗σy+bz∗σz+b0∗I (bx−bx∗)σx+(by−by∗)σy+(bz−bz∗)σz+(b0−b0∗)I=0 bx−bx∗=by−by∗=bz−bz∗=b0−b0∗=0
- 由 Tr(ρ)=1 得 b0=21,将 ρ 写为如下形式
ρ=21(axσx+ayσy+azσz+I)=21((ax,ay,az)⋅σ+I)=21(n⋅σ+I),n∈R3
- 易见 n⋅σ 的特征值为 ±∥n∥ ,进而 ρ 的特征值为 2±∥n∥+1,再由 ρ 的半正定性,得
2±∥n∥+1≥0, i.e. ∥n∥≤1
综上,密度矩阵标准形式为
ρ=21(n⋅σ+I),n∈B3
其中 B3 为三维空间的单位球。
反之,易见上述形式给出的 ρ 为密度矩阵,于是建立对应关系:
S(C2)ρ=21(n⋅σ+I)⟶∼B3⟼n
定义相容性
前边将纯态 ∣ψ⟩ 与 Bloch 球面建立对应;注意到纯态也可以通过密度矩阵建立和 Bloch 球的对应,下边式子说明两种方式得到同一对应。
∣ψ(θ,ϕ)⟩=(cos2θeiϕsin2θ)ρ=∣ψ⟩⟨ψ∣=(cos2θeiϕsin2θ)(cos2θe−iϕsin2θ) =(2cosθ+1eiϕ2sinθe−iϕ2sinθ2cosθ+1) =21(n(θ,ϕ)⋅σ+I)
两种途径定义的 Bloch 向量都是
n(θ,ϕ)=(sinθcosϕ,sinθsinϕ,cosθ)
此外,球投影算子 σα+ 对混合态的测量期望仍是 Bloch 向量在 α+ 方向上的投影。以 σx 为例,由密度矩阵的测量公式,得
⟨σx⟩=Tr(σxρ)=Tr(21σx(n⋅σ+I))=Tr(21σx(axσx+ayσy+azσz+I))=Tr(21(axI+iayσz−iazσy+σx))=Tr(2axI)=ax
ax 为 Bloch 向量 n 的 x 轴坐标,也即 n 在 x 轴方向上的投影。
谱分解与完全混合态
下边从量子态和 Bloch 球两个角度,理解密度矩阵的谱分解和完全混合态。
从 S(H) 上看
-
凸线性组合
称 x=i=1∑rλixi 为向量集 {xi}i=1r 的一个凸线性组合,若 λi≥0 且 i=1∑rλi=1 。
易见密度矩阵的凸线性组合仍是密度矩阵。
-
密度矩阵为正规矩阵,做谱分解
ρ=i=1∑rλi∣λi⟩⟨λi∣=i=1∑rλiρi
其中 ρi=∣λi⟩⟨λi∣ 为纯态的密度矩阵,ρ 为纯态 {ρi}i=1r 的凸线性组合。 一般来说,纯态 ρi 的取法不唯一。
-
考虑 ρ∈S(C2),此时 ρ=λ1ρ1+λ2ρ2,由“不同特征值的特征向量线性无关”知:
- 当 ρ=21I 时,ρ1 和 ρ2 唯一确定
- 当 ρ=21I 时,必有 λ1=λ2=21 且 ρ1 和 ρ2 可取任意两个正交态,此时称 ρ 为完全混合态
从 Bloch 球上看
-
密度矩阵的凸线性组合,对应 Bloch 向量的凸线性组合
ρ1ρ2⇒ρ⇒n=21(n1⋅σ+I)=21(n2⋅σ+I)=λ1ρ1+λ2ρ2=21((λ1n1+λ2n2)⋅σ+I)=λ1n1+λ2n2
-
从 Bloch 球 B3 上看:
- 谱分解将 Bloch 向量 n 分解为方向相反的两个球面向量 n−,n+(正交=反向,纯态=球面)
- 当 n=0 时, n−,n+ 唯一确定,为 n 所在直线与 Bloch 球面的两个交点
- 当 n=0 为 Bloch 球原点时,n 可分解为任意两个方向相反的球面向量,对应完全混合态
酉变换
回顾上篇,演化算子 U 是酉算子,作用量子态上:
- 向量形式 ∣ψ⟩→U∣ψ⟩
- 密度矩阵形式 ρ→UρU†
注意到球投影算子 σα+ 不仅满足 Hermitian 性质,还是个酉算子,因而可作用在量子态上。
从 Bloch 球上看,σα+ 的作用效果是将空间绕 x 轴旋转 180°。以 σx 为例:
let ρσxρσx†(nx,ny,nz)=21(n⋅σ+I)=21σx(n⋅σ+I)σx=21(n⋅σxσσx+I)=21((nx,ny,nz)⋅(σx,−σy,−σz)+I)=21((nx,−ny,−nz)⋅σ+I)⟶rotate 180°(nx,−ny,−nz)
思考题:C2 上既能当成测量,又能当成演化的算子有哪些?或者判断正误:
{σα+∣α+∈S2}=?{σ∈End(C2)∣σ=σ†=σ−1}
Ps:本节提到了三类算子,注意区分
- 密度矩阵:Hermitian + 半正定 + 迹为 1
- 测量算子:Hermitian
- 演化算子:Unitary
- σα+, 视语境可能是测量算子,也可能是演化算子
- Pα=∣α+⟩⟨α+∣,视语境可能是投影算子,也可能是纯态的密度矩阵
总结
本篇建立 C2 量子态与 Bloch 球的对应,通过 Bloch 球的几何性质理解量子态,主要内容:
- 量子比特 ∣0⟩ 和 ∣1⟩
- 纯态 ∣ψ⟩ 标准化,与 Bloch 球面建立对应
- 正交的量子态对应 Bloch 球面方向相反的向量,导出“球投影算子”
- 球投影算子导出泡利矩阵,提供量子系统的“坐标系”
- 借助泡利矩阵,将密度矩阵做“坐标分解”,建立与 Bloch 球的对应
- 进一步讨论:谱分解,完全混合态,测量 + 酉变换
关于 Bloch 球
- Bloch 球只能表示 C2 的量子态,更高维数的量子系统很难找到类似 Bloch 球这样精巧的对应。
- C2n 的分离态可用多个 Bloch 球表示。因此一些特殊的量子门,比如 Walsh-Hadamard 门和量子傅里叶变换,其作用在分离态后仍得到分离态,可用 Bloch 球演示变化过程。
- 对一般的量子门,如果将分离态变换为纠缠态,就无法用 Bloch 球演示了。
拓展-MUB 问题
Bloch 球的几何直观可用于解决一些难题,比如清华视频课里提到的 MUB 问题(Mutually unbiased bases):
称希尔伯特空间 Hd=Cd 的两组标准正交基 {∣ψi⟩}i=1d,{∣φi⟩}i=1d 为“两两平等基”,若
∣⟨φj∣ψi⟩∣2=d1 ,∀ i,j
d = 2 时,由 Bloch 球性质易证,存在三组基元,两两作成 MUB,比如取坐标轴向量
{∣0⟩,∣1⟩},{∣+⟩,∣−⟩},{∣+i⟩,∣−i⟩}
设 Hd 上这样的基元有 M(d) 组,以下问题还是开问题:
- 是否有 M(d)≥2, ∀ d
- 计算 M(d)
已有结论:
M(N)=N+1if N=pr, where p is prime
比如 M(2)=3,M(3)=4,M(4)=5,但 M(6) 还未知。
参考维基:Mutually_unbiased_bases
qisikit