选了课程 《量子计算与量子信息》,没学过量子力学的博主实在是听不懂啊 (ㄒoㄒ)
简略整理了下 可能大概也许 明白一二都没有 的课程最开始两节的内容,如有错误欢迎指出 ~ ~ ~
文章目录
- 矩阵论
- 复空间中的矩阵
- 矩阵上的运算
- 量子力学
- 量子态
- 基本假设
- 量子计算机的线路模型
- 单量子位门
- 双量子位门
- 三量子位门
- Bell 态的制备
- 一些有意思的应用
- 量子密码
- 量子计算
- 可逆计算
- 隐形传态
矩阵论
复空间中的矩阵
复空间
C
n
\mathbb C^n
Cn 上矩阵的厄米(Hermite):
A
†
:
=
(
A
∗
)
T
=
(
A
T
)
∗
A^\dagger := (A^*)^T = (A^T)^*
A†:=(A∗)T=(AT)∗,共轭转置,上角标
†
\dagger
† 读作 \dagger
。
一些特殊的矩阵,
-
厄米矩阵(Hermitian): H † = H H^\dagger = H H†=H,类比实空间中的对称阵 H T = H H^T=H HT=H,特征值都是实数,不同特征值对应的特征向量正交。
-
酉矩阵(Unitary): U † U = I U^\dagger U = I U†U=I,类比实空间中的单位正交阵 U T U = I U^TU=I UTU=I,特征值都是单位复数,不同特征值对应的特征向量正交。
-
正规矩阵(Normal): A A † = A † A AA^\dagger = A^\dagger A AA†=A†A,等价于它有 n n n 个标准正交的特征向量,也等价于它可以酉对角化(谱分解)。厄米矩阵、酉矩阵、对角阵,都是正规矩阵。
-
半正定矩阵:函数 f ( x ) : = x † A x f(x):=x^\dagger Ax f(x):=x†Ax 叫做二次型(Quadratic form),如果任意的 x ∈ C n x \in \mathbb C^n x∈Cn 都有 f ( x ) ≥ 0 f(x) \ge 0 f(x)≥0,则称 A A A 是半正定的(positive semi-definite),它等价于特征值为非负实数。
-
奇异值矩阵(Singular):对于可酉对角化的矩阵 A A A,定义 ∣ A ∣ : = A † A |A|:= \sqrt{A^\dagger A} ∣A∣:=A†A,它的特征值等于 A A A 的特征值绝对值。对称方阵 A † A A^\dagger A A†A 的特征值叫做 A A A 的奇异值,是 A A A 的特征值平方根。
一些算子,
- 投影算子:令 ∣ 1 ⟩ , ⋯ , ∣ d ⟩ |1\rang,\cdots,|d\rang ∣1⟩,⋯,∣d⟩ 是全空间的一组标准正交基,且 ∣ 1 ⟩ , ⋯ , ∣ k ⟩ |1\rang,\cdots,|k\rang ∣1⟩,⋯,∣k⟩ 是 k k k 维子空间的一组基,那么 P = ∑ i = 1 k ∣ i ⟩ ⟨ i ∣ P=\sum_{i=1}^k |i\rang\lang i| P=∑i=1k∣i⟩⟨i∣ 就是投影矩阵,它是个厄米矩阵。
- 投影算子的正交补: Q = I − P Q = I-P Q=I−P 是到 P P P 的正交补空间上的投影矩阵。
矩阵上的运算
一些矩阵运算,
- 矩阵相似:存在可逆阵 P P P,使得 P − 1 A P = B P^{-1}AP=B P−1AP=B
- 酉相似:存在酉阵 U U U,使得 U † A U = B U^\dagger AU = B U†AU=B
- 酉对角化:存在酉阵 U U U,使得 U † A U = D U^\dagger AU = D U†AU=D,其中 D D D 是由特征值组成对角线的对角阵, U U U 恰好是对应的标准正交特征向量。充要条件是 A A A 的正规阵,谱分解为 A = U D U † = ∑ i λ i ∣ i ⟩ ⟨ i ∣ A=UDU^\dagger=\sum_i\lambda_i|i\rang\lang i| A=UDU†=∑iλi∣i⟩⟨i∣
- 同时酉对角化:存在酉阵 U U U,使得 U † A U = D U^\dagger AU = D U†AU=D 且 U † B U = D ′ U^\dagger BU = D' U†BU=D′,其中 D , D ′ D,D' D,D′ 都是对角阵。充分条件是两个正规阵 A , B A,B A,B 可交换。
多种矩阵分解,
- LU分解:某矩阵 A ∈ C m × n A \in C^{m \times n} A∈Cm×n,如果可以写成 A = L U A=LU A=LU 的形式,其中 L L L 是下三角阵, U U U 是上三角阵,则称为三角分解(LU 分解、LR 分解)
- QR分解:某矩阵 A ∈ C m × n A \in C^{m \times n} A∈Cm×n,如果可以写成 A = Q R A=QR A=QR 的形式,其中 Q Q Q 是正交阵, R R R 是上三角阵,则称为 QR 分解
- 奇异值分解:秩 r r r 的 m × n m \times n m×n 矩阵 A A A,存在 m m m 阶酉阵 U U U 和 n n n 阶酉阵 V V V 以及秩 r r r 对角阵 Σ = d i a g ( s 1 , ⋯ , s r ) \Sigma=diag(s_1,\cdots,s_r) Σ=diag(s1,⋯,sr),其中 s 1 ≥ ⋯ ≥ s r > 0 s_1\ge\cdots\ge s_r>0 s1≥⋯≥sr>0,使得 A = U Σ V A=U \Sigma V A=UΣV
- 极式分解:任意线性算子 A A A,存在酉算子 U U U 以及唯一的半正定算子 J : = A † A , K : = A A † J:=\sqrt{A^\dagger A},K:=\sqrt{AA^\dagger} J:=A†A,K:=AA†,使得 A = U J = K U A=UJ=KU A=UJ=KU
- 谱分解:任意正规算子 M M M,存在酉阵 U U U 和对角阵 D D D,使得 M = U D U † M=UDU^\dagger M=UDU†。设 M M M 的不同特征值为 λ 1 , ⋯ , λ s \lambda_1,\cdots,\lambda_s λ1,⋯,λs,代数重数 k 1 , ⋯ , k s k_1,\cdots,k_s k1,⋯,ks,则存在唯一的幂等阵 E 1 , ⋯ , E s E_1,\cdots,E_s E1,⋯,Es(谱族),满足 ∑ i E i = I \sum_i E_i=I ∑iEi=I, r a n k ( E i ) = k i rank(E_i)=k_i rank(Ei)=ki, E i E j = O , i ≠ j E_iE_j=O,i\neq j EiEj=O,i=j,使得 M = ∑ i λ i E i M=\sum_i \lambda_i E_i M=∑iλiEi
量子力学
一些术语:
- 量子系统的状态:希尔伯特空间(完备的内积空间)中的一个向量
- 几率 = 概率,几率幅 = 概率的平方根
- 线性算子 = 矩阵
- 本征值 = 特征值,本征态 = 特征向量
量子态
量子力学中使用 Dirac 符号表示量子态:
- ⟨ a ∣ \langle a| ⟨a∣ 是左矢,叫做 “bra”,它是一个行向量
- ∣ b ⟩ |b\rangle ∣b⟩ 是右矢,叫做 “ket”,它是一个列向量
- 两者关系为 ∣ a ⟩ † = ⟨ a ∣ |a\rangle^\dagger = \langle a| ∣a⟩†=⟨a∣,而 ⟨ a ∣ b ⟩ \langle a|b\rangle ⟨a∣b⟩ 就是向量内积
单个量子比特(qubit)可表示为 ∣ ψ ⟩ = α ∣ 0 ⟩ + β ∣ 1 ⟩ |\psi\rangle = \alpha|0\rangle + \beta|1\rangle ∣ψ⟩=α∣0⟩+β∣1⟩,其中 α , β ∈ C \alpha,\beta \in \mathbb C α,β∈C 都是复数,需满足归一化条件 α α ∗ + β α ∗ = 1 \alpha\alpha^* + \beta\alpha^* = 1 αα∗+βα∗=1
于是
∣
ψ
⟩
=
α
∣
0
⟩
+
β
∣
1
⟩
|\psi\rangle = \alpha|0\rangle + \beta|1\rangle
∣ψ⟩=α∣0⟩+β∣1⟩ 也可写作:
∣
ψ
⟩
=
e
i
γ
(
cos
θ
2
∣
0
⟩
+
e
i
ϕ
sin
θ
2
∣
1
⟩
)
|\psi\rangle = e^{i\gamma}\left( \cos\dfrac{\theta}{2}|0\rangle + e^{i\phi}\sin\dfrac{\theta}{2}|1\rangle \right)
∣ψ⟩=eiγ(cos2θ∣0⟩+eiϕsin2θ∣1⟩)
其中
e
i
γ
e^{i\gamma}
eiγ 是整体相因子,我们认为它不影响 qubit 的值,可以省略,
∣
ψ
⟩
=
cos
θ
2
∣
0
⟩
+
e
i
ϕ
sin
θ
2
∣
1
⟩
|\psi\rangle = \cos\dfrac{\theta}{2}|0\rangle + e^{i\phi}\sin\dfrac{\theta}{2}|1\rangle
∣ψ⟩=cos2θ∣0⟩+eiϕsin2θ∣1⟩
那么 ∣ ψ ⟩ |\psi\rangle ∣ψ⟩ 就是由 θ ∈ [ − π , π ) , ϕ ∈ [ 0 , 2 π ) \theta \in [-\pi,\pi), \phi \in [0,2\pi) θ∈[−π,π),ϕ∈[0,2π) 定义的单位球面上的一个点。使用 Bloch Sphere 表示,其 Bloch 矢量为 ( sin θ cos ϕ , sin θ sin ϕ , cos θ ) (\sin\theta\cos\phi, \sin\theta\sin\phi,\cos\theta) (sinθcosϕ,sinθsinϕ,cosθ)
单量子比特系统的基态(单位正交基)为:
∣
0
⟩
=
[
1
0
]
∈
C
2
,
∣
1
⟩
=
[
0
1
]
∈
C
2
|0\rangle = \begin{bmatrix}1\\0\end{bmatrix} \in \mathbb C^2, |1\rangle = \begin{bmatrix}0\\1\end{bmatrix} \in \mathbb C^2
∣0⟩=[10]∈C2,∣1⟩=[01]∈C2
我们记
∣
a
b
⟩
:
=
∣
a
⟩
⊗
∣
b
⟩
|ab\rangle := |a\rangle \otimes |b\rangle
∣ab⟩:=∣a⟩⊗∣b⟩ 是两个量子比特串的直积(矩阵的 Kronecker 积),例如双量子比特系统的基态为:
∣
00
⟩
,
∣
01
⟩
,
∣
10
⟩
,
∣
11
⟩
|00\rangle, |01\rangle, |10\rangle, |11\rangle
∣00⟩,∣01⟩,∣10⟩,∣11⟩
或者选取 ”EPR 对“ 里的 Bell 态作为基态:
∣
ϕ
+
⟩
:
=
∣
00
⟩
+
∣
11
⟩
2
∣
ϕ
−
⟩
:
=
∣
00
⟩
−
∣
11
⟩
2
∣
ψ
+
⟩
:
=
∣
01
⟩
+
∣
10
⟩
2
∣
ψ
−
⟩
:
=
∣
01
⟩
−
∣
10
⟩
2
|\phi^+\rangle := \dfrac{|00\rangle+|11\rangle}{\sqrt2}\\ |\phi^-\rangle := \dfrac{|00\rangle-|11\rangle}{\sqrt2}\\ |\psi^+\rangle := \dfrac{|01\rangle+|10\rangle}{\sqrt2}\\ |\psi^-\rangle := \dfrac{|01\rangle-|10\rangle}{\sqrt2}\\
∣ϕ+⟩:=2∣00⟩+∣11⟩∣ϕ−⟩:=2∣00⟩−∣11⟩∣ψ+⟩:=2∣01⟩+∣10⟩∣ψ−⟩:=2∣01⟩−∣10⟩
以此类推, n n n 量子比特系统的基态有 2 n 2^n 2n 个(指数级)。
量子态叠加:一个量子态可以写成若干基态的线性组合,
∣
ψ
⟩
=
∑
i
p
i
∣
i
⟩
|\psi\rangle = \sum_{i} p_i|i\rangle
∣ψ⟩=i∑pi∣i⟩
其中 i ∈ { 0 , 1 } n i \in \{0,1\}^n i∈{0,1}n 是基态的索引,几率幅 p i p_i pi 满足 ∑ i p i 2 = 1 \sum_i p_i^2=1 ∑ipi2=1
量子态纠缠:如果 ∣ ψ ⟩ |\psi\rangle ∣ψ⟩ 不可以写成 n n n 个单量子比特的直积 ∣ ψ ⟩ = ∣ ϕ 1 ⋯ ϕ n ⟩ |\psi\rangle = |\phi_1\cdots\phi_n\rangle ∣ψ⟩=∣ϕ1⋯ϕn⟩ 的形式,那么这叫做纠缠态。例如 Bell 态,不同位置上的量子比特,因为纠缠而彼此影响(相关性极强)。
基本假设
- 波函数假设:系统状态的描述为 ψ ( r , t ) \psi(r,t) ψ(r,t),参数是位置和时间,也记作 ∣ ψ ⟩ |\psi\rangle ∣ψ⟩
- 算符假设:力学量可以用线性厄米算符表示(也就是厄米矩阵),共轭力学量算符不可对易(对易子 [ A , B ] : = A B − B A ≠ 0 [A,B]:=AB-BA\neq 0 [A,B]:=AB−BA=0)
- 测量假设:对力学可观测量的测量,使得系统以几率
∣
⟨
ϕ
m
∣
ψ
⟩
∣
2
|\langle\phi_m|\psi\rangle|^2
∣⟨ϕm∣ψ⟩∣2 随机落入该力学量的某一个本征态
∣
ϕ
m
⟩
|\phi_m\rangle
∣ϕm⟩(也就是
∣
ψ
⟩
=
∑
m
⟨
ϕ
m
∣
ψ
⟩
⋅
∣
ϕ
m
⟩
|\psi\rangle = \sum_m \langle\phi_m|\psi\rangle \cdot| \phi_m\rangle
∣ψ⟩=∑m⟨ϕm∣ψ⟩⋅∣ϕm⟩,其中的
⟨
ϕ
m
∣
ψ
⟩
\langle\phi_m|\psi\rangle
⟨ϕm∣ψ⟩ 是几率幅)
- 大量相同量子态的系统,构成了量子系综。
- 两力学量可同时观测,当仅当它们的力学量算符可对易 [ A , B ] = 0 [A,B]=0 [A,B]=0,测量后系统进入两力学量的某个共同本征态。
- 态演化假设:量子态遵循的演化方程为 Schrodinger 方程
- 全同性假设:全同粒子体系的波函数对于任意两个粒子的交换,要么是对称的(玻色子),要么是反对称的(费米子)
量子计算机的线路模型
量子运算需要保持量子态的归一化条件,因此每个逻辑门都是酉矩阵 U † U = I U^\dagger U=I U†U=I(可逆运算),具体为 ∣ ψ ⟩ ↦ U ∣ ψ ⟩ |\psi\rangle \mapsto U|\psi\rangle ∣ψ⟩↦U∣ψ⟩(左乘变换矩阵)。
Deutsch 定理:任意的 d d d 维酉变换,总是可以分解为 2 d 2 − d 2d^2-d 2d2−d 个二级酉变换( d d d 维酉阵,仅改变两个叠加分量上的几率幅)的乘积,并且可以用一系列的单量子位门( U ( α , ϕ ) U(\alpha,\phi) U(α,ϕ))以及双量子位门(CNOT 门)依次作用来实现。
单量子位门
经典信息论中,只有 “逻辑非” 这一个非凡逻辑门。而在量子信息理论中,由于定义在希尔伯特空间上,有特别多的非平凡的单量子逻辑门。
三个 Pauli 矩阵
σ
⃗
=
(
σ
x
,
σ
y
,
σ
z
)
\vec\sigma=(\sigma_x,\sigma_y,\sigma_z)
σ=(σx,σy,σz),
X
=
[
0
1
1
0
]
,
Y
=
[
0
−
i
i
0
]
,
Z
=
[
1
0
0
−
1
]
X = \begin{bmatrix} 0&1\\ 1&0 \end{bmatrix}, Y = \begin{bmatrix} 0&-i\\ i&0 \end{bmatrix}, Z = \begin{bmatrix} 1&0\\ 0&-1 \end{bmatrix}
X=[0110],Y=[0i−i0],Z=[100−1]
其中 X X X 是比特翻转(逻辑非),而 Z Z Z 是相位翻转。
阿达玛门(Hadamard),
H
=
1
2
[
1
1
1
−
1
]
H = \dfrac{1}{\sqrt 2}\begin{bmatrix} 1&1\\ 1&-1 \end{bmatrix}
H=21[111−1]
相位门(Phase),
S
=
[
1
0
0
i
]
S = \begin{bmatrix} 1&0\\ 0&i \end{bmatrix}
S=[100i]
T 门(
π
8
\dfrac{\pi}{8}
8π 门),
T
=
[
1
0
0
exp
(
π
4
i
)
]
T = \begin{bmatrix} 1&0\\ 0&\exp(\dfrac{\pi}{4}i) \end{bmatrix}
T=[100exp(4πi)]
易知 T 2 = S T^2=S T2=S, S 2 = Z S^2=Z S2=Z, Z 2 = Y 2 = X 2 = I Z^2=Y^2=X^2=I Z2=Y2=X2=I,
根据欧拉公式
e
i
A
θ
=
I
cos
(
θ
)
+
i
A
sin
(
θ
)
,
∀
A
2
=
I
e^{iA\theta}= I\cos(\theta) + iA\sin(\theta), \forall A^2=I
eiAθ=Icos(θ)+iAsin(θ),∀A2=I,绕轴
x
,
y
,
z
,
n
⃗
x,y,z,\vec n
x,y,z,n 旋转的旋转算子分别为:
R
x
(
θ
)
:
=
e
−
i
X
θ
/
2
=
[
cos
(
θ
/
2
)
−
i
sin
(
θ
/
2
)
−
i
sin
(
θ
/
2
)
cos
(
θ
/
2
)
]
R
y
(
θ
)
:
=
e
−
i
Y
θ
/
2
=
[
cos
(
θ
/
2
)
−
sin
(
θ
/
2
)
sin
(
θ
/
2
)
cos
(
θ
/
2
)
]
R
z
(
θ
)
:
=
e
−
i
Z
θ
/
2
=
[
e
−
i
θ
/
2
0
0
e
i
θ
/
2
]
R
n
⃗
(
θ
)
:
=
e
−
i
θ
n
⃗
σ
⃗
/
2
=
cos
(
θ
/
2
)
I
−
i
sin
(
θ
/
2
)
(
n
x
X
+
n
y
Y
+
n
z
Z
)
\begin{aligned} R_x(\theta) &:= e^{-iX\theta/2} = \begin{bmatrix} \cos(\theta/2) & -i\sin(\theta/2)\\ -i\sin(\theta/2)& \cos(\theta/2) \end{bmatrix}\\ R_y(\theta) &:= e^{-iY\theta/2} = \begin{bmatrix} \cos(\theta/2) & -\sin(\theta/2)\\ \sin(\theta/2)& \cos(\theta/2) \end{bmatrix}\\ R_z(\theta) &:= e^{-iZ\theta/2} = \begin{bmatrix} e^{-i\theta/2} & 0\\ 0 & e^{i\theta/2} \end{bmatrix}\\ R_{\vec n}(\theta) &:= e^{-i\theta\vec n\vec\sigma/2} = \cos(\theta/2)I - i\sin(\theta/2)(n_xX+n_yY+n_zZ)\\ \end{aligned}
Rx(θ)Ry(θ)Rz(θ)Rn(θ):=e−iXθ/2=[cos(θ/2)−isin(θ/2)−isin(θ/2)cos(θ/2)]:=e−iYθ/2=[cos(θ/2)sin(θ/2)−sin(θ/2)cos(θ/2)]:=e−iZθ/2=[e−iθ/200eiθ/2]:=e−iθnσ/2=cos(θ/2)I−isin(θ/2)(nxX+nyY+nzZ)
1 2 \dfrac{1}{2} 21 旋量:自旋两周(角度 θ = 4 π \theta=4\pi θ=4π)才回到本身(单位阵 I I I)。
定理:单量子比特的任意酉变换
U
U
U,都存在实参
α
,
β
,
γ
,
δ
\alpha,\beta,\gamma,\delta
α,β,γ,δ,使得:
U
=
e
i
α
R
x
(
β
)
R
y
(
γ
)
R
z
(
δ
)
U = e^{i\alpha} R_x(\beta)R_y(\gamma)R_z(\delta)
U=eiαRx(β)Ry(γ)Rz(δ)
双量子位门
受控非门(CNOT),运算为
∣
A
B
⟩
↦
∣
A
⟩
∣
A
⊕
B
⟩
|AB\rangle \mapsto |A\rangle|A \oplus B\rangle
∣AB⟩↦∣A⟩∣A⊕B⟩,其中
∣
A
⟩
|A\rangle
∣A⟩ 是控制位,而
∣
B
⟩
|B\rangle
∣B⟩ 是靶比特,
C
(
X
)
=
[
1
0
0
0
0
1
0
0
0
0
0
1
0
0
1
0
]
=
[
I
0
0
X
]
C(X) = \begin{bmatrix} 1&0&0&0\\ 0&1&0&0\\ 0&0&0&1\\ 0&0&1&0 \end{bmatrix} = \begin{bmatrix} I&0\\ 0&X \end{bmatrix}
C(X)=
1000010000010010
=[I00X]
交换门(Swap),运算为
∣
A
B
⟩
↦
∣
B
A
⟩
|AB\rangle \mapsto |BA\rangle
∣AB⟩↦∣BA⟩,则
[
1
0
0
0
0
0
1
0
0
1
0
0
0
0
0
1
]
=
C
C
′
C
\begin{bmatrix} 1&0&0&0\\ 0&0&1&0\\ 0&1&0&0\\ 0&0&0&1 \end{bmatrix} = CC'C
1000001001000001
=CC′C
其中 C ′ C' C′ 是以 ∣ B ⟩ |B\rangle ∣B⟩ 为控制位 ∣ A ⟩ |A\rangle ∣A⟩ 为靶比特的 CNOT 门。注意,Swap 过程中并没有复制 qubit,不受量子不可克隆性影响。
一般受控变换,运算为
∣
c
⟩
∣
t
⟩
↦
∣
c
⟩
U
c
∣
t
⟩
|c\rangle|t\rangle \mapsto |c\rangle U^c|t\rangle
∣c⟩∣t⟩↦∣c⟩Uc∣t⟩,其中
U
U
U 是任意的单量子比特酉变换,对应的酉阵为:
C
(
U
)
=
I
⊗
U
c
=
[
U
c
0
0
U
c
]
=
[
I
0
0
U
]
C(U) = I \otimes U^c = \begin{bmatrix} U^c & 0\\ 0 & U^c \end{bmatrix} = \begin{bmatrix} I & 0\\ 0 & U \end{bmatrix}
C(U)=I⊗Uc=[Uc00Uc]=[I00U]
而对于运算
∣
t
⟩
∣
c
⟩
↦
U
c
∣
t
⟩
∣
c
⟩
|t\rangle|c\rangle \mapsto U^c|t\rangle |c\rangle
∣t⟩∣c⟩↦Uc∣t⟩∣c⟩,则有(只需考虑四个基态):
U
=
[
a
b
c
d
]
,
U
c
⊗
I
=
[
1
a
b
1
c
d
]
U=\begin{bmatrix} a&b\\ c&d \end{bmatrix}, U^c \otimes I = \left[\begin{array}{cc|cc} 1\\ &a&&b\\ \hline &&1\\ &c&&d \end{array}\right]
U=[acbd],Uc⊗I=
1ac1bd
注意:受控酉变换作用之后,控制位也可能变化,而靶比特也可能不变!例如:受控相移,运算为
∣
c
⟩
∣
t
⟩
↦
∣
c
⟩
(
e
i
α
)
c
∣
t
⟩
|c\rangle|t\rangle \mapsto |c\rangle (e^{i\alpha})^c|t\rangle
∣c⟩∣t⟩↦∣c⟩(eiα)c∣t⟩,
I
⊗
[
e
i
α
0
0
e
i
α
]
c
=
[
1
0
0
e
i
α
]
⊗
I
I \otimes \begin{bmatrix} e^{i\alpha}&0\\ 0&e^{i\alpha} \end{bmatrix}^c = \begin{bmatrix} 1&0\\ 0&e^{i\alpha} \end{bmatrix} \otimes I
I⊗[eiα00eiα]c=[100eiα]⊗I
假设酉矩阵的 ABC 分解为:
U
=
e
i
α
A
B
C
U = e^{i\alpha} ABC
U=eiαABC(矩阵作用从右向左),那么就有:
C
(
U
)
=
(
[
1
0
0
e
i
α
]
⊗
I
)
⋅
(
I
⊗
A
)
⋅
C
(
X
)
⋅
(
I
⊗
B
)
⋅
C
(
X
)
⋅
(
I
⊗
C
)
C(U) = \left(\begin{bmatrix} 1&0\\ 0&e^{i\alpha} \end{bmatrix} \otimes I\right) \cdot (I \otimes A) \cdot C(X) \cdot (I \otimes B) \cdot C(X) \cdot (I \otimes C)
C(U)=([100eiα]⊗I)⋅(I⊗A)⋅C(X)⋅(I⊗B)⋅C(X)⋅(I⊗C)
如图所示:
三量子位门
双控非门(Toffoli),运算为
∣
a
b
c
⟩
↦
∣
a
b
⟩
∣
c
⊕
a
b
⟩
|abc\rangle \mapsto |ab\rangle|c\oplus ab\rangle
∣abc⟩↦∣ab⟩∣c⊕ab⟩,
C
2
(
X
)
=
[
I
I
I
X
]
C^2(X) = \begin{bmatrix} I\\ &I\\ &&I\\ &&&X \end{bmatrix}
C2(X)=
IIIX
受控交换门(Fredkin),运算为
∣
c
⟩
∣
a
b
⟩
↦
∣
c
⟩
∣
b
a
⟩
|c\rangle|ab\rangle \mapsto |c\rangle|ba\rangle
∣c⟩∣ab⟩↦∣c⟩∣ba⟩,
[
1
1
1
1
1
0
0
0
0
0
1
0
0
1
0
0
0
0
0
1
]
\begin{bmatrix} 1\\ &1\\ &&1\\ &&&1\\ &&&& 1&0&0&0\\ &&&& 0&0&1&0\\ &&&& 0&1&0&0\\ &&&& 0&0&0&1 \end{bmatrix}
11111000001001000001
Toffoli 门和 Fredkin 门,都是可逆计算的通用逻辑门。
一般双控门,运算为
∣
a
b
c
⟩
↦
∣
a
b
⟩
U
a
b
∣
c
⟩
|abc\rangle \mapsto |ab\rangle U^{ab}|c\rangle
∣abc⟩↦∣ab⟩Uab∣c⟩,
C
2
(
U
)
=
I
4
⊗
U
a
b
=
[
U
a
b
U
a
b
U
a
b
U
a
b
]
=
[
I
I
I
U
]
C^2(U) = I_4 \otimes U^{ab} = \begin{bmatrix} U^{ab}\\ &U^{ab}\\ &&U^{ab}\\ &&&U^{ab}\\ \end{bmatrix} = \begin{bmatrix} I\\ &I\\ &&I\\ &&&U \end{bmatrix}
C2(U)=I4⊗Uab=
UabUabUabUab
=
IIIU
Bell 态的制备
令
U
=
C
(
X
)
⋅
(
H
⊗
I
)
U = C(X) \cdot (H \otimes I)
U=C(X)⋅(H⊗I),那么容易验证
U
∣
00
⟩
=
∣
ϕ
+
⟩
U
∣
01
⟩
=
∣
ψ
+
⟩
U
∣
10
⟩
=
∣
ϕ
−
⟩
U
∣
11
⟩
=
∣
ψ
+
⟩
U|00\rangle = |\phi^+\rangle\\ U|01\rangle = |\psi^+\rangle\\ U|10\rangle = |\phi^-\rangle\\ U|11\rangle = |\psi^+\rangle\\
U∣00⟩=∣ϕ+⟩U∣01⟩=∣ψ+⟩U∣10⟩=∣ϕ−⟩U∣11⟩=∣ψ+⟩
一些有意思的应用
量子力学的神秘之源:态叠加原理,它导致了
- 量子态测量的不确定性:量子密码([BB84],信息论安全的一次一密)
- 量子态操作的并行性:量子计算机(经典信息就是量子信息的基态)
量子密码
由于量子态的叠加是线性的,因此运算 ∣ ϕ ⟩ ∣ 0 ⟩ ↦ ∣ ϕ ⟩ ∣ ϕ ⟩ |\phi\rangle|0\rangle \mapsto |\phi\rangle|\phi\rangle ∣ϕ⟩∣0⟩↦∣ϕ⟩∣ϕ⟩ 对于任意量子态并不总是存在,这直接导致了量子态不可克隆定理。根据测量假设,对量子态的测量结果是基态的真随机分布(遵照几率幅)。
量子秘钥分发(QKD, [BB84])协议如下:
量子计算
量子计算机的能力: B Q P BQP BQP 类(有界错误的概率量子图灵机可在 PPT 时间求解的问题的收集),
- 可以有效求解 P P P 类(经典计算机就是量子计算机的退化)
- 不能有效求解 P S P A C E PSPACE PSPACE 类以外的(指数级内存开销,无能为力)
- 尚不清楚它处于 P P P 和 P S P A C E PSPACE PSPACE 之间的什么位置(也不知道是否可以有效求解 N P NP NP 类)
- 量子可计算函数类,也满足 Church-Turing 题论(这与经典可计算函数类完全相同,仅仅是效率上的不同)
量子态可叠加原理,导致了量子计算机的并行性。然而,其计算结果是一个量子态(量子比特),一旦进行观测就会依照波函数几率幅,坍缩到某个基态(经典比特)上。波函数中包含的指数级信息不可以被全部访问,只能观测到一个全局信息。
Shor 算法巧妙地利用了 FFT 可将周期性信息集中到某个频率的振幅上的特点,可有效求解整数分解和离散对数问题。假设这两者不属于 P P P 类,那么量子计算机确实比经典计算机更强。但是,它们不属于 N P C NPC NPC 类,量子计算机能否有效求解 N P NP NP 类我们依然不知道。除了求解困难问题,量子计算机还可以有效模拟量子系统,但是如何提取期望信息也是个问题。
想要实现量子计算机,就需要在真实世界中寻找一些物理现象,使得它们表现为酉算子(或者仅仅是量子操作,与外部环境作用后移除)。在物理实现的多种方式中,离子阱(Ion Traps)似乎是目前表现最好的。
量子系统总会与外部环境作用(退相干),因此如何进行量子容错计算,也是个重要问题。一种方法是使用量子纠错码,例如 Shor 码、CSS 码、Steane 码。当然也有其他的不使用纠错码的容错技术。
可逆计算
Landauer 原理:在物理计算机中,每擦除一比特或者执行一个扇入运算,则至少增加环境热量 k B T ln 2 k_BT\ln2 kBTln2 焦耳, k B k_B kB 是玻尔兹曼常数。
信息是物理的,这解决了 Maxwell 佯谬(热力学的看门小精灵)。
如果执行可逆运算(例如酉算子),那么这就可以绕过上述的热力学限制,使得计算所消耗的能量大幅下降。
隐形传态
物体的信息包括经典信息与量子信息,
- EPR 制备中心产生处于 Bell 态 ∣ ψ − ⟩ A B |\psi^-\rangle_{AB} ∣ψ−⟩AB 的纠缠粒子对 A , B A,B A,B,通过信道分发给 Alice 和 Bob(光速)
- Alice 对粒子 A A A 和粒子 C C C 形成的系统利用 Bell 基进行测量(定域量子操作),这会发生纠缠交换(超光速的超距作用,使得 A , C A,C A,C 纠缠,而 A , B A,B A,B 不再纠缠),然后 Alice 将测量得到的经典信息通过信道传递给 Bob(光速)
- Bob 根据收到的经典信息,选择一个合适的定域酉变换( I , X , Z , Z X I, X, Z, ZX I,X,Z,ZX),对粒子 B B B 做变换使之称为粒子 C C C(的量子态)。