量子计算与量子信息读书笔记——第1章简介与概述(下)

注:这是全文的后半部分,前文传送门:
https://blog.csdn.net/qq_37638320/article/details/129971731

量子算法

经典计算

原理:任何经典电路都可以用等价的仅含可逆原件的,由可逆门Toffoli门构成的电路代替

量子电路不能直接用于模拟经典电路:酉量子逻辑门本质上是可逆的,而许多经典逻辑门本质上是不可逆的

Toffoli门

  • 概念:控制比特、目标比特
  • 作用:当两个控制比特都置为1时目标比特翻转,否则不变: ( a , b , c ) → ( a , b , c ⊕ a b ) (a,b,c)\to(a,b,c\oplus ab) (a,b,c)(a,b,cab)
  • 矩阵表示: [ 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 ] \begin{bmatrix} 1& 0& 0& 0& 0& 0& 0&0 \\ 0& 1& 0& 0& 0& 0& 0&0 \\ 0& 0& 1& 0& 0& 0& 0&0 \\ 0& 0& 0& 1& 0& 0& 0&0 \\ 0& 0& 0& 0& 1& 0& 0&0 \\ 0& 0& 0& 0& 0& 1& 0&0 \\ 0& 0& 0& 0& 0& 0& 0&1 \\ 0& 0& 0& 0& 0& 0& 1&0 \end{bmatrix} 1000000001000000001000000001000000001000000001000000000100000010
  • 电路表示:image-20230403222401773
  • 性质:
    1. 可逆门,逆为自身
    2. 可用于模拟不可逆经典逻辑门,并且保证量子计算机可以进行任何经典(确定性)计算机能够完成的计算
NAND扇出
image-20230403223222745image-20230403223521779
上方两比特表示NAND门的输入,第3个比特设置为标准状态1,称为辅助态第2个比特是输入(0/1),其它两个比特是标准辅助态,扇出的输出在第2个和第3个比特

有了这两个操作就可以模拟电路中的其它元件,所以任意经典电路都可以被等价的可逆电路模拟

非确定经典计算机

image-20230403230902442

∣ 0 ⟩ \ket 0 0通过一个阿达玛门得到 ∣ 0 ⟩ + ∣ 1 ⟩ 2 \dfrac{\ket0+\ket1}{\sqrt2} 2 0+1,再测量这个态,结果是 ∣ 0 ⟩ , ∣ 1 ⟩ \ket0,\ket1 0,1 50 % 50\% 50%的概率。为量子计算机提供了有效模拟不确定的经典计算机的能力

量子并行性

内涵:量子并行性允许量子计算机同时计算在不同 x x x值下的函数值 f ( x ) f(x) f(x)

原理:设 f ( x ) : { 0 , 1 } → { 0 , 1 } f(x):\{0,1\}\to\{0,1\} f(x):{0,1}{0,1},定义 U f : ∣ x , y ⟩ → ∣ x , y ⊕ f ( x ) ⟩ U_f:\ket{x,y}\to\ket{x,y\oplus f(x)} Uf:x,yx,yf(x)

  1. 第一个寄存器被称为数据寄存器,第二个寄存器被称为目标寄存器
  2. 酉变换
  3. y = 0时,第二个比特的终态是值 f ( x ) f(x) f(x)
  4. 暂且将 U f U_f Uf看作一个黑盒

image-20230403232032247

在数据寄存器中制备叠加态 ∣ 0 ⟩ + ∣ 1 ⟩ 2 \dfrac{\ket0+\ket1}{\sqrt2} 2 0+1,然后作用 U f U_f Uf
∣ 0 ⟩ + ∣ 1 ⟩ 2 ∣ 0 ⟩ = 1 2 ∣ 00 ⟩ + 1 2 ∣ 10 ⟩ → 1 2 ∣ 0 , f ( 0 ) ⟩ + 1 2 ∣ 1 , f ( 1 ) ⟩ \dfrac{\ket0+\ket1}{\sqrt2}\ket0= \dfrac{1}{\sqrt2}\ket{00}+\dfrac{1}{\sqrt2}\ket{10}\to\dfrac{1}{\sqrt2}\ket{0,f(0)}+\dfrac{1}{\sqrt2}\ket{1,f(1)} 2 0+10=2 100+2 1102 10,f(0)+2 11,f(1)
得到态 ∣ 0 , f ( 0 ) ⟩ + ∣ 1 , f ( 1 ) ⟩ 2 \dfrac{\ket{0,f(0)}+\ket{1,f(1)}}{\sqrt2} 2 0,f(0)+1,f(1):同时计算了 f ( 0 ) , f ( 1 ) f(0),f(1) f(0),f(1)的值,与经典的并行性用多个电路同时计算 f ( x ) f(x) f(x)不同,利用量子计算机处于不同状态的叠加态的能力,单个 f ( x ) f(x) f(x)电路用来同时计算多个 x x x的函数值

任意比特函数的推广

image-20230403233213226

初态制备为 ∣ 0 ⟩ \ket0 0,将2个阿达玛门同时作用在n比特上(记作 H ⊗ 2 H^{\otimes2} H2),输出为
( ∣ 0 ⟩ + ∣ 1 ⟩ 2 ) ( ∣ 0 ⟩ + ∣ 1 ⟩ 2 ) = ∣ 00 ⟩ + ∣ 01 ⟩ + ∣ 10 ⟩ + ∣ 11 ⟩ 2 (\dfrac{\ket0+\ket1}{\sqrt2})(\dfrac{\ket0+\ket1}{\sqrt2})=\dfrac{\ket{00}+\ket{01}+\ket{10}+\ket{11}}{2} (2 0+1)(2 0+1)=200+01+10+11
一般地,当初态均为 ∣ 0 ⟩ \ket0 0时, n n n比特上均作用阿达玛门(记作 H ⊗ n H^{\otimes n} Hn)后得到
1 2 n ∑ x ∣ x ⟩ \dfrac{1}{\sqrt{2^n}}\sum\limits_x\ket x 2n 1xx
n比特输入x和单比特输出 f ( x ) f(x) f(x)函数的量子并行性计算:制备 n + 1 n+1 n+1比特的量子态 ∣ 0 ⟩ ⊗ n ∣ 0 ⟩ \ket0^{\otimes n}\ket0 0n0,然后将阿达玛变换作用在前 n n n个量子比特上,再通过量子电路 U f U_f Uf,得到态 1 2 n ∑ x ∣ x ⟩ ∣ f ( x ) ⟩ \dfrac{1}{\sqrt{2^n}}\sum\limits_x\ket x\ket{f(x)} 2n 1xxf(x)(即 ∣ 0 ⟩ ⊗ n ∣ 0 ⟩ ⟶ H ⊗ n ⟶ U f 1 2 n ∑ x ∣ x ⟩ ∣ f ( x ) ⟩ \ket0^{\otimes n}\ket0\overset{H^{\otimes n}}{\longrightarrow}\overset{U_f}{\longrightarrow}\dfrac{1}{\sqrt{2^n}}\sum\limits_x\ket x\ket{f(x)} 0n0HnUf2n 1xxf(x),计算f一次计算f的所有值)

局限性:测量 ∑ x ∣ x , f ( x ) ⟩ \sum_x\ket{x,f(x)} xx,f(x)只能得到一个 x x x的函数值 f ( x ) f(x) f(x)(需要更高的信息抽取能力)

Deutsch算法

原理:将量子并行性和量子力学中的重要性质相干性结合起来

image-20230403234818416

输入态 ∣ ψ 0 ⟩ = ∣ 01 ⟩ \ket{\psi_0}=\ket{01} ψ0=01

通过两个阿达玛门后得到 ∣ ψ 1 ⟩ = [ ∣ 0 ⟩ + ∣ 1 ⟩ 2 ] [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] \ket{\psi_1}=\left[\dfrac{\ket0+\ket1}{\sqrt2}\right]\left[\dfrac{\ket0-\ket1}{\sqrt2}\right] ψ1=[2 0+1][2 01]

易知, ∣ x ⟩ ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 ⟶ U f ∣ x ⟩ ( ∣ 0 ⊕ f ( x ) ⟩ − ∣ 1 ⊕ f ( x ) ⟩ ) / 2 = { ∣ x ⟩ ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 , f ( x ) = 0 , − ∣ x ⟩ ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 , f ( x ) = 1. = ( − 1 ) f ( x ) ∣ x ⟩ ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 \ket x(\ket0-\ket1)/\sqrt2\overset{U_f}{\longrightarrow}\ket x(\ket{0\oplus f(x)}-\ket{1\oplus f(x)})/\sqrt2=\left\{\begin{array}{ll}\ket x(\ket0-\ket1)/\sqrt2,&f(x)=0,\\-\ket x(\ket0-\ket1)/\sqrt2,&f(x)=1.\end{array}\right.=(-1)^{f(x)}\ket x(\ket0-\ket1)/\sqrt2 x(01)/2 Ufx(0f(x)1f(x))/2 ={x(01)/2 ,x(01)/2 ,f(x)=0,f(x)=1.=(1)f(x)x(01)/2

所以
∣ ψ 2 ⟩ = U f ( ∣ 0 ⟩ 2 ∣ 0 ⟩ − ∣ 1 ⟩ 2 + ∣ 1 ⟩ 2 ∣ 0 ⟩ − ∣ 1 ⟩ 2 ) = 1 2 U f ( ∣ 0 ⟩ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ) + 1 2 U f ( ∣ 1 ⟩ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ) = 1 2 ( − 1 ) f ( 0 ) ∣ 0 ⟩ ∣ 0 ⟩ − ∣ 1 ⟩ 2 + 1 2 ( − 1 ) f ( 1 ) ∣ 1 ⟩ ∣ 0 ⟩ − ∣ 1 ⟩ 2 = [ ( − 1 ) f ( 0 ) ∣ 0 ⟩ + ( − 1 ) f ( 1 ) ∣ 1 ⟩ 2 ] [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] = { ± [ ∣ 0 ⟩ + ∣ 1 ⟩ 2 ] [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] , f ( 0 ) = f ( 1 ) , ± [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] , f ( 0 ) ≠ f ( 1 ) . \begin{array}{ll}\ket{\psi_2}&=U_f(\dfrac{\ket0}{\sqrt2}\dfrac{\ket0-\ket1}{\sqrt2}+\dfrac{\ket1}{\sqrt2}\dfrac{\ket0-\ket1}{\sqrt2})\\ &=\dfrac{1}{\sqrt2}U_f(\ket0\dfrac{\ket0-\ket1}{\sqrt2})+\dfrac{1}{\sqrt2}U_f(\ket1\dfrac{\ket0-\ket1}{\sqrt2})\\ &=\dfrac{1}{\sqrt2}(-1)^{f(0)}\ket0\dfrac{\ket0-\ket1}{\sqrt2}+\dfrac{1}{\sqrt2}(-1)^{f(1)}\ket1\dfrac{\ket0-\ket1}{\sqrt2}\\ &=\left[\dfrac{(-1)^{f(0)}\ket0+(-1)^{f(1)}\ket1}{\sqrt2}\right]\left[\dfrac{\ket0-\ket1}{\sqrt2}\right]\\ &=\left\{ \begin{array}{ll} \pm\left[\dfrac{\ket0+\ket1}{\sqrt2}\right]\left[\dfrac{\ket0-\ket1}{\sqrt2}\right],&f(0)=f(1),\\ \pm\left[\dfrac{\ket0-\ket1}{\sqrt2}\right]\left[\dfrac{\ket0-\ket1}{\sqrt2}\right],&f(0)\not=f(1). \end{array} \right. \end{array} ψ2=Uf(2 02 01+2 12 01)=2 1Uf(02 01)+2 1Uf(12 01)=2 1(1)f(0)02 01+2 1(1)f(1)12 01=[2 (1)f(0)0+(1)f(1)1][2 01]= ±[2 0+1][2 01],±[2 01][2 01],f(0)=f(1),f(0)=f(1).
最后阿达玛门作用在第一个量子比特得
∣ ψ 3 ⟩ = { ± ∣ 0 ⟩ [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] , f ( 0 ) = f ( 1 ) , ± ∣ 1 ⟩ [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] , f ( 0 ) ≠ f ( 1 ) . \ket{\psi_3}=\left\{ \begin{array}{ll} \pm\ket0\left[\dfrac{\ket0-\ket1}{\sqrt2}\right],&f(0)=f(1),\\ \pm\ket1\left[\dfrac{\ket0-\ket1}{\sqrt2}\right],&f(0)\not=f(1). \end{array} \right. ψ3= ±0[2 01],±1[2 01],f(0)=f(1),f(0)=f(1).
{ f ( 0 ) ⊕ f ( 1 ) = 0 , f ( 0 ) = f ( 1 ) , f ( 0 ) ⊕ f ( 1 ) = 1 , f ( 0 ) ≠ f ( 1 ) . \left\{\begin{array}{ll}f(0)\oplus f(1)=0,&f(0)=f(1),\\f(0)\oplus f(1)=1,&f(0)\not= f(1).\end{array}\right. {f(0)f(1)=0,f(0)f(1)=1,f(0)=f(1),f(0)=f(1).

∣ ψ 3 ⟩ = ± ∣ f ( 0 ) ⊕ f ( 1 ) ⟩ [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] \ket{\psi_3}=\pm\ket{f(0)\oplus f(1)}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right] ψ3=±f(0)f(1)[2 01]

通过测量第一个量子比特,我们可以确定 f ( 0 ) ⊕ f ( 1 ) f(0)\oplus f(1) f(0)f(1).

结论:量子电路有能力通过计算一次 f ( x ) f(x) f(x)来确定 f ( x ) f(x) f(x)全局性质

众多量子算法其设计的本质都是在精心选择函数和最终变换,使得能高效地确定有用的全局信息

Deutsch-Jozsa算法

Deutsch问题:Alice从 0 ∼ 2 n − 1 0\sim2^n-1 02n1中选择一个数 x x x寄给Bob,Bob计算出 f ( x ) f(x) f(x)的值寄回Alice.其中 f ( x ) f(x) f(x)值非零即一,且 f ( x ) f(x) f(x)要么对所有 x x x均为常数,要么为平衡函数,即对于所有可能的 x x x一半取 1 1 1,另一半取 0 0 0。问Alice最少通信多少次能成功?

确定性经典算法: 2 n − 1 + 1 2^{n-1}+1 2n1+1

量子计算:交换量子比特,使用一个酉变换 U f U_f Uf计算 f ( x ) f(x) f(x),1次

image-20230404201635672.png

前n个量子比特为查询寄存器,最后一个量子比特为答案寄存器。

输入态 ∣ ψ 0 ⟩ = ∣ 0 ⟩ ⊗ n ∣ 1 ⟩ \ket{\psi_0}=\ket0^{\otimes n}\ket1 ψ0=0n1

经过阿达玛门作用,有 ∣ ψ 1 ⟩ = ∑ x ∈ { 0 , 1 } n ∣ x ⟩ 2 n [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] \ket{\psi_1}=\sum\limits_{x\in\{0,1\}^n}\dfrac{\ket x}{\sqrt{2^n}}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right] ψ1=x{0,1}n2n x[2 01]

由上知, ∣ x ⟩ ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 ⟶ U f ( − 1 ) f ( x ) ∣ x ⟩ ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 \ket x(\ket0-\ket1)/\sqrt2\overset{U_f}{\longrightarrow}(-1)^{f(x)}\ket x(\ket0-\ket1)/\sqrt2 x(01)/2 Uf(1)f(x)x(01)/2

故使用 U f : ∣ x , y ⟩ → ∣ x , y ⊕ f ( x ) ⟩ U_f:\ket{x,y}\to\ket{x,y\oplus f(x)} Uf:x,yx,yf(x)
∣ ψ 2 ⟩ = U f ( ∑ x ∣ x ⟩ 2 n [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] ) = ∑ x 1 2 n U f ( ∣ x ⟩ [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] ) = ∑ x ( − 1 ) f ( x ) ∣ x ⟩ 2 n [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] \begin{array}{ll} \ket{\psi_2}&=U_f(\sum\limits_{x}\dfrac{\ket x}{\sqrt{2^n}}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right])\\ &=\sum\limits_{x}\dfrac{1}{\sqrt{2^n}}U_f(\ket x\left[\dfrac{\ket0-\ket1}{\sqrt2}\right])\\ &=\sum\limits_{x}\dfrac{(-1)^{f(x)}\ket x}{\sqrt{2^n}}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right] \end{array} ψ2=Uf(x2n x[2 01])=x2n 1Uf(x[2 01])=x2n (1)f(x)x[2 01]
H ∣ 0 ⟩ = ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 = ( ( − 1 ) 0 ∗ 0 ∣ 0 ⟩ + ( − 1 ) 0 ∗ 1 ∣ 1 ⟩ ) / 2 , H ∣ 1 ⟩ = ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 = ( ( − 1 ) 1 ∗ 0 ∣ 0 ⟩ + ( − 1 ) 1 ∗ 1 ∣ 1 ⟩ ) / 2 H\ket0=(\ket0+\ket1)/\sqrt2=((-1)^{0*0}\ket0+(-1)^{0*1}\ket1)/\sqrt2,H\ket1=(\ket0-\ket1)/\sqrt2=((-1)^{1*0}\ket0+(-1)^{1*1}\ket1)/\sqrt2 H0=(0+1)/2 =((1)000+(1)011)/2 ,H1=(01)/2 =((1)100+(1)111)/2

所以 H ∣ x ⟩ = ∑ z ∈ { 0 , 1 } ( − 1 ) x z ∣ z ⟩ / 2 H\ket x=\sum_{z\in\{0,1\}}(-1)^{xz}\ket z/\sqrt2 Hx=z{0,1}(1)xzz/2

推广到多量子比特,有 H ⊗ n ∣ x 1 , ⋯   , x n ⟩ = ∏ i = 1 n ∑ z i ∈ { 0 , 1 } ( − 1 ) x i z i ∣ z i ⟩ 2 = ∑ z 1 , ⋯   , z n ( − 1 ) x 1 z 1 + ⋯ + x n z n ∣ z 1 , ⋯   , z n ⟩ 2 n H^{\otimes n}\ket{x_1,\cdots,x_n}=\prod\limits_{i=1}^n\sum\limits_{z_i\in\{0,1\}}\dfrac{(-1)^{x_iz_i}\ket{z_i}}{\sqrt2}=\dfrac{\sum\limits_{z_1,\cdots,z_n}(-1)^{x_1z_1+\cdots+x_nz_n}\ket{z_1,\cdots,z_n}}{\sqrt{2^n}} Hnx1,,xn=i=1nzi{0,1}2 (1)xizizi=2n z1,,zn(1)x1z1++xnznz1,,zn

x ⋅ z x\cdot z xz x , z x,z x,z二进制表示下按位内积,则 H ⊗ n ∣ x ⟩ = ∑ z ( − 1 ) x ⋅ z ∣ z ⟩ 2 n H^{\otimes n}\ket x=\dfrac{\sum_z(-1)^{x\cdot z}\ket z}{\sqrt{2^n}} Hnx=2n z(1)xzz

利用该结果计算 ∣ ψ 3 ⟩ \ket{\psi_3} ψ3
∣ ψ 3 ⟩ = H ⊗ n ∣ ψ 2 ⟩ = ∑ x ( − 1 ) f ( x ) 2 n H ⊗ n ∣ x ⟩ [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] = ∑ x ( − 1 ) f ( x ) 2 n ∑ z ( − 1 ) x ⋅ z ∣ z ⟩ 2 n [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] = ∑ x ∑ z ( − 1 ) x ⋅ z + f ( x ) ∣ z ⟩ 2 n [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] = ∑ z ∑ x ( − 1 ) x ⋅ z + f ( x ) ∣ z ⟩ 2 n [ ∣ 0 ⟩ − ∣ 1 ⟩ 2 ] \begin{array}{ll} \ket{\psi_3}&=H^{\otimes n}\ket{\psi_2}\\ &=\sum\limits_x\dfrac{(-1)^{f(x)}}{\sqrt{2^n}}H^{\otimes n}\ket x\left[\dfrac{\ket0-\ket1}{\sqrt2}\right]\\ &=\sum\limits_x\dfrac{(-1)^{f(x)}}{\sqrt{2^n}}\sum\limits_z\dfrac{(-1)^{x\cdot z}\ket z}{\sqrt{2^n}}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right]\\ &=\sum\limits_x\sum\limits_z\dfrac{(-1)^{x\cdot z+f(x)}\ket z}{2^n}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right]\\ &=\sum\limits_z\sum\limits_x\dfrac{(-1)^{x\cdot z+f(x)}\ket z}{2^n}\left[\dfrac{\ket0-\ket1}{\sqrt2}\right] \end{array} ψ3=Hnψ2=x2n (1)f(x)Hnx[2 01]=x2n (1)f(x)z2n (1)xzz[2 01]=xz2n(1)xz+f(x)z[2 01]=zx2n(1)xz+f(x)z[2 01]
现对查询寄存器进行观察:

f ( x ) f(x) f(x)为常函数时,态 ∣ 0 ⟩ ⊗ n \ket0^{\otimes n} 0n对应的振幅为 ∑ x ( − 1 ) x ⋅ 00 ⋯ 0 ‾ + f ( x ) / 2 n = ∑ x ( − 1 ) f ( x ) / 2 n = ( − 1 ) f ( x ) = ± 1 \sum_x(-1)^{x\cdot \overline{00\cdots0}+f(x)}/2^n=\sum_x(-1)^{f(x)}/2^n=(-1)^{f(x)}=\pm1 x(1)x000+f(x)/2n=x(1)f(x)/2n=(1)f(x)=±1

由于 ∣ ψ 3 ⟩ \ket{\psi_3} ψ3为单位长度,所以其它振幅均为0,观测会使得查询寄存器中所有量子比特均为0;

f ( x ) f(x) f(x)为平衡函数时,态 ∣ 0 ⟩ ⊗ n \ket0^{\otimes n} 0n的振幅为 ∑ x ( − 1 ) f ( x ) / 2 n = 0 \sum_x(-1)^{f(x)}/2^n=0 x(1)f(x)/2n=0,测量会使得查询寄存器至少有一个量子比特不是0.

局限性:

  1. 尚未找到应用场景
  2. 与经典算法可比性差,计算函数的方式很不同
  3. 概率经典计算机同样可以优于确定性模型
量子算法总结
基于傅里叶变换的量子算法

离散傅里叶变换: N N N元复数集 x 0 , ⋯   , x N − 1 x_0,\cdots,x_{N-1} x0,,xN1到复数集 y 0 , ⋯   , y N − 1 y_0,\cdots,y_{N-1} y0,,yN1的变换,其中
y k ≡ 1 N ∑ j = 0 N − 1 e 2 π i j k / N x j y_k\equiv\dfrac{1}{\sqrt N}\sum\limits_{j=0}^{N-1}e^{2\pi ijk/N}x_j ykN 1j=0N1e2πijk/Nxj
量子力学形式:在计算基 ∣ j ⟩ \ket j j上定义 n n n量子比特的线性变换 U ( 0 ⩽ j ⩽ 2 n − 1 ) U(0\leqslant j\leqslant 2^n-1) U(0j2n1)
∣ j ⟩ → 1 2 n ∑ k = 0 2 n − 1 e 2 π i j k / 2 n ∣ k ⟩ \ket j\to\dfrac{1}{\sqrt{2^n}}\sum\limits_{k=0}^{2^n-1}e^{2\pi ijk/2^n}\ket k j2n 1k=02n1e2πijk/2nk
易验证其为酉变换。

作用于叠加态:
∑ j = 0 2 n − 1 x j ∣ j ⟩ → 1 2 n ∑ k = 0 2 n − 1 [ ∑ j = 0 2 n − 1 e 2 π i j k / 2 n x j ] ∣ k ⟩ = ∑ k = 0 2 n − 1 y k ∣ k ⟩ \sum\limits_{j=0}^{2^n-1}x_j\ket j\to\dfrac{1}{\sqrt{2^n}}\sum\limits_{k=0}^{2^n-1}\left[\sum\limits_{j=0}^{2^n-1}e^{2\pi ijk/2^n}x_j\right]\ket k=\sum\limits_{k=0}^{2^n-1}y_k\ket k j=02n1xjj2n 1k=02n1[j=02n1e2πijk/2nxj]k=k=02n1ykk
效率:经典 O ( N log ⁡ N ) = n 2 n O(N\log N)=n2^n O(NlogN)=n2n,量子 O ( log ⁡ 2 N ) = n 2 O(\log^2N)=n^2 O(log2N)=n2

局限:无法直接得到变换后的结果 y k y_k yk

实例:Deutsch-Jozsa算法、Shor算法

量子搜索算法

问题:给定一个大小为 N N N的状态空间,没有关于其结构的先验知识,我们想要在搜索空间中找到一个满足已知性质的元素

效率:经典 O ( N ) O(N) O(N),量子 O ( N ) O(\sqrt N) O(N )

量子模拟

概念:模拟自然发生的量子系统( k n kn kn个量子比特)

局限:有效模拟并不意味着得到量子系统中所期望得到的信息(在测量后,一个 k n kn kn个量子比特的模拟将会坍缩到一个确定的态,只留下 k n kn kn比特的信息;在波函数中 c n c^n cn比特的“隐藏信息”不能全部访问)

问题:研究能够有效抽取期望得到的答案的方法

作用:获取其它量子算法灵感的一般途径

摩尔定律的量子推论:如果能够每两年向量子计算机增加一个量子比特,量子计算机将保持与经典计算机相同的步伐

量子计算的能力

计算复杂性理论

  • 复杂性类:一系列计算问题的集合,所有这些问题在求解所需的计算资源上具有相同的性质

    复杂性类说明
    P在经典计算机中能够快速求解
    NP在经典计算机中能够快速验证解
    NP-complete一个算法如果能够求解一个特殊的NP完全问题,那么稍微增加一点计算资源就可以被用于求解其它NP问题
    PSPACE能用有限空间资源求解的问题
    BPP如果允许有界的错误概率,能够用随机算法在多项式时间内求解
    BQP能够被量子计算机有效求解的计算问题(允许有界的错误概率)
  • 经典与量子复杂性类的关系:

    image-20230404224814559

现有结论:

  1. 尽管量子计算机可以求解素因子分解问题(普遍相信在NP中),但由于不确定该问题是否为NP完全问题,无法知道量子计算机是否可以求解所有NP问题;且基于搜索的方法并不能有效的给出所有NP问题的解
  2. 量子计算机能有效地求解P类问题,但PSPACE以外的问题不能被有效求解
  3. 量子计算机计算能力严格强于经典计算机这一命题的证明很可能是不平凡的

实验量子信息处理

Stern-Gerlach实验

目的:表明量子比特结构是两能级量子系统

原理:氢原子具有一个质子和一个环绕的电子,可以看作质子周围的一点“电流”,该电流使原子具有磁场,每个原子都有“磁偶极矩”,使得每个原子表现得像一块小的条形磁铁,其轴线对应于电子旋转的轴。把小磁铁扔过一个磁场会导致磁铁被磁场偏转,我们期望在实验中看到类似的原子偏转。

注:受当时技术水平限制,最初的实验是用Ag完成的实验,Ag具有复杂的结构,解释复杂效应的难度更大

内容:热原子从烤箱中“射出”,通过一个导致原子偏转的磁场,然后记录每个原子的位置。

image-20230405103134682

原子如何偏转取决于原子的磁偶极矩(即电子旋转的轴)以及Stern-Gerlach装置产生的磁场。通过适当地构造Stern-Gerlach器件,我们可以使原子偏转一定量,而该量取决于原子磁偶极矩的 z z z分量,其中 z z z是一些固定的外轴

现象与解释:

  1. 离开炉子的热原子从一组离散的角度出现:磁偶极矩是量子化的
  2. 氢原子具有零磁矩,但观测到了两个光束,一个被磁场向上 ( ∣ + Z ⟩ ) (\ket{+Z}) (+Z)偏转,另一个则向下 ( ∣ − Z ⟩ ) (\ket{-Z}) (Z)偏转:氢原子中的电子有一个叫做自旋的量,电子的自旋对氢原子的磁偶极矩产生额外的贡献

进一步实验以描述电子自旋:

image-20230405105202980

将两个Stern-Gerlach设备级联在一起,将第二个设备侧向倾斜,使得磁场沿 x ^ \hat x x^偏转原子;阻止第一台设备的 ∣ − Z ⟩ \ket{-Z} Z输出,而将 ∣ + Z ⟩ \ket{+Z} +Z输出发送至沿着 x ^ \hat x x^轴的第二设备,将检测器放置在最终处以测量沿 x ^ \hat x x^轴的原子分布

现象与解释:

  1. 指向 + z ^ +\hat z +z^方向的经典磁偶极子在 x x x方向没有净磁矩,但实验上观察到两个强度相等的峰值:通过第二个设备的每个原子可以被描述成 ∣ + Z ⟩ ∣ + X ⟩ \ket{+Z}\ket{+X} +Z+X ∣ + Z ⟩ ∣ − X ⟩ \ket{+Z}\ket{-X} +ZX,来指代可能被观测到的自旋的两个值

image-20230405105830259

通过发送之前的输出使之通过第二个 z ^ \hat z z^导向的Stern-Gerlach设备

现象与解释:

  1. 在最终输出处再次观察到两个强度相等的光束:一个 ∣ + Z ⟩ \ket{+Z} +Z态由相等比例的 ∣ + X ⟩ \ket{+X} +X ∣ − X ⟩ \ket{-X} X态组成,一个 ∣ + X ⟩ \ket{+X} +X态由相等比例的 ∣ + Z ⟩ \ket{+Z} +Z ∣ − Z ⟩ \ket{-Z} Z态组成

量子比特模型解释:

∣ 0 ⟩ , ∣ 1 ⟩ \ket0,\ket1 0,1表示量子比特的状态,并分配 ∣ + Z ⟩ ← ∣ 0 ⟩ , ∣ − Z ⟩ ← ∣ 1 ⟩ , ∣ + X ⟩ ← ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 , ∣ − X ⟩ ← ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 \ket{+Z}\leftarrow\ket0,\ket{-Z}\leftarrow\ket1,\ket{+X}\leftarrow(\ket0+\ket1)/\sqrt2,\ket{-X}\leftarrow(\ket0-\ket1)/\sqrt2 +Z0,Z1,+X(0+1)/2 ,X(01)/2 ,假设 z ^ \hat z z^Stern-Gerlach装置在计算基 ∣ 0 ⟩ , ∣ 1 ⟩ \ket0,\ket1 0,1下测量自旋(即量子比特), x ^ \hat x x^Stern-Gerlach实验装置在基 ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 , ( ∣ 0 ⟩ − ∣ 1 ⟩ ) / 2 (\ket0+\ket1)/\sqrt2,(\ket0-\ket1)/\sqrt2 (0+1)/2 ,(01)/2 下测量自旋。

在级联 z ^ − x ^ − z ^ \hat z-\hat x-\hat z z^x^z^实验中,从第一个实验射出的自旋处于状态 ∣ + Z ⟩ = ∣ 0 ⟩ = ( ∣ + X ⟩ + ∣ − X ⟩ ) / 2 \ket{+Z}=\ket0=(\ket{+X}+\ket{-X})/\sqrt2 +Z=0=(+X+X)/2 ,那么第二个装置得到 ∣ + X ⟩ \ket{+X} +X的概率是 1 2 \dfrac{1}{2} 21;相似地,第三个装置得到 ∣ + Z ⟩ \ket{+Z} +Z的概率是 1 2 \dfrac{1}{2} 21.因此,量子比特模型正确预言了级联Stern-Gerlach实验的结果。

问题:没有证明量子比特是毫无疑问的理解电子自旋的正确方式(由于此类实验很多,现在相信量子比特是描述电子自旋的最好模型)

更进一步,我们相信量子比特模型以及它向更高维度的推广(量子力学)能够描述每个物理系统

实用量子信息处理的前景

问题:是否存在某种原理禁止我们进行一种或多种形式的量子信息处理?

  1. 噪声可能对有用的量子信息处理构成根本性障碍

    量子纠错码理论:量子噪声需要解决,但不是根本性的原理问题

    量子计算的阈值定理:如果量子计算机中的噪声水平可以降低到某个常数“阈值”以下,那么就可以使用量子纠错码来进一步地降低噪声,只需要很小的计算复杂性开销,基本上可以降低到任意小

  2. 量子力学可能是不正确的

应用:既然构建量子信息处理设备没有根本性的障碍,为什么我们要投入大量的时间和金钱这样做?

  1. 【小规模】量子态层析:确定系统的量子状态方法,通过重复制备同一个量子态,然后以不同的方式测量,以建立量子态的完整描述

    量子过程层析成像:完全表征量子系统的动态学

  2. 小规模通信原语:分发少量需要高度安全的关键材料

  3. 【中等规模】量子系统模拟:确定材料的基本性质;实验室设计和测试新分子性质的工具(狭义第一性原理计算)

  4. 【大规模】大整数素因子分解,计算离散对数,量子搜索

实现:

  1. 光学技术(电磁辐射)
  2. 囚禁不同原子(离子阱,中性原子阱,电磁辐射)——精湛的状态制备和量子测量
  3. 核磁共振NMR——极好的动态演化

量子信息处理不仅仅是另一种信息处理技术:量子计算是信息处理的一个抽象范式,可能在技术上有许多不同的实现。人们可以比较量子计算的两个不同方案的技术优点,即使量子计算机的一个非常糟糕的方案,它与精湛设计的经典计算机也具有定性的本质不同

量子信息

内涵

  • 广义:与使用量子力学的信息处理相关的所有操作
  • 具体:对基本量子信息处理任务的研究(量子信息理论)

基本目标

  1. 确定量子力学中的基本静态资源类:例如量子比特、经典比特、贝尔态
  2. 确定量子力学中动力学过程的基本类:例如内存(在一段时间内存储量子态的能力)、量子信息传输、复制(或试图复制)量子态、保护量子信息处理免受噪声影响
  3. 量化执行基本动态过程所需的资源折衷

核心问题

  1. 是什么使量子信息处理成为可能?
  2. 是什么分离了量子世界与经典世界?
  3. 量子计算中正在利用哪些经典世界中无法获得的资源?
量子信息理论问题
使用量子信道传输经典信息

基本结果:香农无噪声信道编码定理,有噪声信道编码定理

信息源:由一组概率 p j , j = 1 , 2 , ⋯   , d p_j,j=1,2,\cdots,d pj,j=1,2,,d描述,信源每次随机地以概率 p j p_j pj产生“字母” j j j,信源每次使用都相互独立

香农无噪声信道编码定理

  • 香农熵:由概率 p j p_j pj描述的经典信源可以被压缩,以使得平均每次使用信源可用 H ( p j ) H(p_j) H(pj)位信息来表示,其中 H ( p j ) = − ∑ j p j log ⁡ ( p j ) H(p_j)=-\sum_jp_j\log(p_j) H(pj)=jpjlog(pj)是信源概率分布的函数,即香农熵(如果尝试比这更少的比特来对信源进行压缩时会导致在信息解压缩时出错的概率很高)

  • 基本目标的实现:

    目标结果
    目标1确定了两个静态资源:比特和信息源
    目标2确定了两阶段动态过程:压缩信息源,解压缩以恢复信息源
    目标3找到最优数据压缩方案确定消耗资源的量化标准

香农有噪声信道编码定理

  • 量化了可以通过有噪声信道可靠传输的信息量(空间/时间(存在噪声时信息的存储问题))

  • 纠错码:在信道发送的信息中引入足够多的冗余,这样即使在某些信息被破坏之后,仍然可以恢复原始消息(提供了计算任意噪声信道容量的通用过程)

  • 基本目标的实现:

    目标结果
    目标1确定了两个静态资源:比特和信息源
    目标2确定了三个动态过程:信道中的噪声,使用纠错码对状态执行编码和解码
    目标3对于固定噪声模型,香农定理告诉我们如果要实现可靠的信息传输,最佳纠错方案必须引入多少冗余

量子信息理论问题

  • 问题1:使用量子态作为介质传输经典信息会发生什么?

    结论:在无噪声信道中传输信息时,量子比特不会导致所需的通信量有任何显著的节省

  • 问题2:通过带有噪声的量子信道传输经典信息会发生什么?

    HSW定理:提供带有噪声的量子信道容量的下界

  • 问题3:是否可以使用纠缠态编码来提高容量?

    迄今为止的所有证据表明,使用纠缠态编码无助于提高容量。但证明或证否仍悬而未决

通过量子信道传输量子信息

量子信息源:由一组概率 p j p_j pj和相应的量子态 ∣ ψ j ⟩ \ket{\psi_j} ψj描述,信源的每次使用都有概率 p j p_j pj产生状态 ∣ ψ j ⟩ \ket{\psi_j} ψj,信源每次使用相互独立

压缩信源输出

  • 经典数据压缩标准技术:考虑单比特量子信源,以概率 p p p输出状态 ∣ 0 ⟩ \ket0 0,以概率 1 − p 1-p 1p输出状态 ∣ 1 ⟩ \ket1 1,与发射单比特的经典源相同,使用与经典类似的技术来压缩信息源,只需要 H ( p , 1 − p ) H(p,1-p) H(p,1p)量子比特来存储压缩信息;若信源以概率 p p p输出状态 ∣ 0 ⟩ \ket0 0,以概率 1 − p 1-p 1p输出状态 ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 (\ket0+\ket1)/\sqrt2 (0+1)/2 将不再适用,由于一般来说无法区分状态 ∣ 0 ⟩ , ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 \ket0,(\ket0+\ket1)/\sqrt2 0,(0+1)/2

  • Schumacher无噪声信道编码定理:

    • 保真度:信源产生的量子态可能会因为压缩——解压缩过程略微失真,为了量化失真,为压缩方案引入保真度,度量由压缩方案引入的平均失真

    • Schumacher无噪声信道编码定理量化了在接近1的保真度恢复信息源的限制下,进行量子数据压缩所需的资源

      信源产生正交量子态:定理退化为信源可以被压缩但不超过经典极限 H ( p j ) H(p_j) H(pj)

      信源产生非正交量子态:量子信源可以被压缩到冯·诺依曼熵,当且仅当状态 ∣ ψ j ⟩ \ket{\psi_j} ψj彼此正交时,冯·诺依曼熵与香农熵一致;否则信源 ( p j , ∣ ψ j ⟩ ) (p_j,\ket{\psi_j}) (pj,ψj)的冯·诺依曼熵通常严格小于香农熵 H ( p j ) H(p_j) H(pj)

    • 实例:对以概率 p p p输出状态 ∣ 0 ⟩ \ket0 0,以概率 1 − p 1-p 1p输出状态 ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 (\ket0+\ket1)/\sqrt2 (0+1)/2 的信源进行压缩

      信源使用 n n n次( n n n为很大的数),由大数定理,信源以高概率发送 n p np np ∣ 0 ⟩ \ket0 0 n ( 1 − p ) n(1-p) n(1p) ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 (\ket0+\ket1)/\sqrt2 (0+1)/2 ,在重新将系统排序的意义下,有
      ∣ 0 ⟩ ⊗ n p ( ∣ 0 ⟩ + ∣ 1 ⟩ 2 ) ⊗ n ( 1 − p ) \ket0^{\otimes np}\left(\dfrac{\ket0+\ket1}{\sqrt2}\right)^{\otimes n(1-p)} 0np(2 0+1)n(1p)
      假设将右侧的 ∣ 0 ⟩ + ∣ 1 ⟩ \ket0+\ket1 0+1展开,由于 n ( 1 − p ) n(1-p) n(1p)很大,由大数定理推断乘积项中大约一半是0,一半是1,可以用如下形式的状态的叠加来近似:
      ∣ 0 ⟩ ⊗ n ( 1 − p ) / 2 ∣ 1 ⟩ ⊗ n ( 1 − p ) / 2 \ket0^{\otimes n(1-p)/2}\ket1^{\otimes n(1-p)/2} 0n(1p)/21n(1p)/2
      因此信源发送的态可由以下形式的叠加态来估计
      ∣ 0 ⟩ ⊗ n ( 1 + p ) / 2 ∣ 1 ⟩ ⊗ n ( 1 − p ) / 2     ( ∗ ) \ket0^{\otimes n(1+p)/2}\ket1^{\otimes n(1-p)/2}\ \ \ (*) 0n(1+p)/21n(1p)/2   ()
      这种形式的态一共有 C n n ( 1 + p ) / 2 C_n^{n(1+p)/2} Cnn(1+p)/2,由斯特林公式近似为 N ≡ C n n ( 1 + p ) / 2 = 2 n H [ ( 1 + p ) / 2 , ( 1 − p ) / 2 ] N\equiv C_n^{n(1+p)/2}=2^{nH[(1+p)/2,(1-p)/2]} NCnn(1+p)/2=2nH[(1+p)/2,(1p)/2]

      压缩方法:将 ( ∗ ) (*) ()的所有态标记为 ∣ c 1 ⟩ ∼ ∣ c N ⟩ \ket{c_1}\sim\ket{c_N} c1cN

      压缩: ∀ 1 ⩽ j ⩽ N \forall 1\leqslant j\leqslant N ∀1jN j j j表示 n H [ ( 1 + p ) / 2 , ( 1 − p ) / 2 ] nH[(1+p)/2,(1-p)/2] nH[(1+p)/2,(1p)/2]位二进制数。在信源发出的 n n n量子比特上执行酉变换 U : ∣ c j ⟩ ↦ ∣ j ⟩ ∣ 0 ⟩ n − n H [ ( 1 + p ) / 2 , ( 1 − p ) / 2 ] U:\ket{c_j}\mapsto\ket j\ket0^{n-nH[(1+p)/2,(1-p)/2]} U:cjj0nnH[(1+p)/2,(1p)/2],然后将最后 n − n H [ ( 1 + p ) / 2 , ( 1 − p ) / 2 ] n-nH[(1+p)/2,(1-p)/2] nnH[(1+p)/2,(1p)/2]个量子比特丢弃,剩下 n H [ ( 1 + p ) / 2 , ( 1 − p ) / 2 ] nH[(1+p)/2,(1-p)/2] nH[(1+p)/2,(1p)/2]量子比特的压缩态 ∣ j ⟩ \ket j j

      解压缩:在压缩态的末尾添上 ∣ 0 ⟩ n − n H [ ( 1 + p ) / 2 , ( 1 − p ) / 2 ] \ket0^{n-nH[(1+p)/2,(1-p)/2]} 0nnH[(1+p)/2,(1p)/2],然后执行 U − 1 U^{-1} U1

    • 理解:利用状态不正交的事实,导致它们比正交状态有更多的物理相似性,利用冗余可以实现数据压缩

量子可区分性

量子力学并不总是可以区分任何状态:例如,量子力学不允许任何过程可靠地区分状态 ∣ 0 ⟩ , ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 \ket0,(\ket0+\ket1)/\sqrt2 0,(0+1)/2 .假设试图通过计算基上测量来区分时,虽然测量结果为1意味着状态必须是 ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 (\ket0+\ket1)/\sqrt2 (0+1)/2 ,但我们无法从测量结果为0推断出有关量子态身份的任何信息

非正交量子态的不可区分性:量子态包含无法通过测量获得的隐藏信息(量子信息理论的核心之一是发展度量来量化非正交量子态的可区分度)

  1. 若我们可以区分任意量子态,则使用纠缠能够比光速更快地进行通信(Bob能判断出与之共享量子比特的Alice在哪组基下进行了测量)
  2. 量子货币(在钞票中嵌入量子序列,由于伪造者无法确定的得到其中的量子比特状态而不破坏它们;通过在 ( ∣ 0 ⟩ , ∣ 1 ⟩ ) (\ket0,\ket1) (0,1) ( ∣ + ⟩ , ∣ − ⟩ ) (\ket+,\ket-) (+,)上测量来验证真伪)
纠缠的产生与转化

相关问题:

  1. 产生纠缠:如果双方要创建一个在他们之间共享的特定纠缠态,假设它们之间没有共享纠缠,那么双方必须交换多少量子比特?

    分布式量子计算可以简单视为用于在两方或者更多方之间产生纠缠的方法;实现这样的分布式量子计算的通信量下界也是实现相应的纠缠态所需要的通信量的下界

  2. 转化纠缠:假设Alice和Bob之间共享了一个贝尔态,他们希望将其转换为其它类型的纠缠状态,他们需要什么资源才能完成这项任务?

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值