• 検索結果がありません。

ノイズがある場合の量子状態識別

N/A
N/A
Protected

Academic year: 2021

シェア "ノイズがある場合の量子状態識別"

Copied!
29
0
0

読み込み中.... (全文を見る)

全文

(1)

ノイズがある場合の量子状態識別

2010

2

福井大学 工学部 物理工学科

04380525

若林  隼吾

(2)

目 次

1章 序論 2

2章 量子力学の性質とBloch 3

2.1 重ね合わせの原理 . . . . 3

2.2 量子状態の測定 . . . . 3

2.3 Bloch . . . . 4

3章 ノイズがない場合の量子状態識別 6 3.1 qubit1個の場合 . . . . 6

3.2 qubitN個の場合 . . . . 7

3.2.1 個別測定(多数決法) . . . . 8

3.2.2 個別測定(適応法) . . . . 9

3.2.3 集団測定 . . . . 12

3.3 識別成功確率の比較 . . . . 13

4章 ノイズがある場合の量子状態識別 14 4.1 ノイズがある場合の量子状態 . . . . 14

4.2 qubit1個の場合 . . . . 15

4.3 qubitN個の場合 . . . . 16

4.3.1 個別測定(多数決法) . . . . 16

4.3.2 個別測定(適応法) . . . . 18

4.4 識別成功確率の比較 . . . . 21

5章 まとめ 25 5.1 結論 . . . . 25

5.2 今後の課題 . . . . 26

参考文献 27

謝辞 28

(3)

1

章 序論

量子情報の分野では量子コンピュータや量子暗号など の研究が進められ 、既存の技 術では成し得ない超並列処理が可能な計算機や絶対安全な通信の実用化が期待されて いる。ではそれらの技術の基盤となる、量子情報とはいったいど のようなものなのだ ろうか。

現在我々が日常的に使用しているコンピュータは「0」または「1」の値をとる「ビッ ト 」というものを操作することにより計算を行っている。情報処理分野ではこの古典 ビットを用いることによって様々な処理を行っているのである。しかし量子情報の分野 では古典ビットではなく量子ビット(qubit)というものを用いる。この量子ビットは量 子力学の性質である重ね合わせの原理を用いるため「0」と「1」状態だけでなく「0」

と「1」の重ね合わせの状態をとることができる。この量子力学の性質を利用したqubit を用いることによって量子コンピュータや量子暗号などの技術が可能になるのである。

そして、このように量子力学の性質を情報処理に積極的に使用するのが量子情報とい う分野である。

量子情報では量子力学の性質を利用するため重ね合わせの原理だけでなく不確定性 原理や測定後の状態が壊れるといった性質も絡んでくる。そのため、ある量子状態の

qubitに正規直交基底で測定を行うと測定は確率的で測定結果が確定することはなく、

測定後の状態は壊れてしまうのである。したがって同じqubitを何度も測定するという ことはできない。そのため最適な量子状態の識別方法というものが重要になってくる。

そこで本研究では、qubitがある確率で2つの状態のどちらかに準備されているとき、

この2つの状態を識別する最適な測定方法について考える。まず、qubit1個の場合 の最適な測定方法の確認からはじ まり、qubitの個数を増やしていったときに 、qubit を個々に測定し多数決で状態を識別する個別測定(多数決法)や前のqubitに対する測 定結果の情報から測定方法を変更する個別測定(適応法)、そしてqubit全体を一つの系 とみなし測定する集団測定といった測定方法での識別成功確率はど う変化するのかを 確認する。そして準備された状態にノイズかあった場合にそれぞれの測定方法に対す る影響を調べることを本研究の目的とし 、その結果から有効な測定方法を検討する。

(4)

2

章 量子力学の性質と

Bloch

この章では序論で触れた量子力学の性質である重ね合わせの原理や状態の測定と、量 子状態の図示に用いるBloch球について説明する。

2.1

重ね合わせの原理

量子情報で用いるqubit|0i|1iの状態しか取れない古典ビットとは異なり、2 つの基本的な状態|0i|1iの重ね合わせの状態をとることができる。その重ね合わ せの状態は次のように表される。

|φi=α|0i+β|1i (2.1)

ここで αβは複素数で 、|α|2 + |β|2 = 1で規格化されている。そして基本的状態

|0i,|1iは大きさ1がで互いに直交している正規直交基底である。

h0|0i=h1|1i= 1 (2.2)

h0|1i=h1|0i= 0 (2.3)

2.2

量子状態の測定

(2.1)式の状態|φiを正規直交基底{|0i,|1i}で測定すると図2.1のように測定後の 状態は|α|2|0iとなり、|β|2|1iとなる。すなわち測定は確率的で状態を壊す。

|0>

|1>

>=ƒ`|0>+ƒ´|1> or

2.1: 正規直交基底{|0i,|1i}での測定

(5)

2.3 Bloch

qubitの状態を視覚的にとらえるのに便利な方法としてBloch球を用いる方法がある。

まず、qubitの状態(2.1)式は長さ1に規格化されているので、次のように書くことが できる。

|φi=e

(

cosθ

2|0i+esinθ 2|1i

)

(2.4) ここでeは観測時の確率に影響を及ぼさないので無視することができる。

したがって(2.4)式は次のように表すことができる。

|φi= cosθ

2|0i+esinθ

2|1i (2.5)

θφ3次元単位球面上の点を定義し 、状態は図2.2に示す球面上の1点で表される。

x

y z

Ĭ

ƒ

n

|0>

|1>

2.2: Bloch

2.2の球面をBloch球という。状態|0iが球面上の北極に対応し 状態|1iが南極に 対応している。そしてBloch球上の単位ベクトルnBlochベクトルと呼ぶ。一般に

Blochベクトルで表記したqubitの状態は次のようになる。

|ni= cosθ

2|0i+esinθ

2|1i (2.6)

(6)

次にこのBlochベクトルを用いてqubitの内積公式を求める。まず|ni|n0iの内 積は

hn|n0i =

(

cosθ

2h0|+esinθ 2h1|

) (

cosθ0

2|0i+e0sinθ0 2|1i

)

=

(

cosθ 2cosθ0

2 + cos (φ−φ0) sinθ 2sinθ0

2

)

−i

(

sin (φ−φ0) sinθ 2sinθ0

2

)

(2.7) となる。したがって内積の2乗は次のように表される。

|hn|n0i|2 = ¯¯¯¯

¯ (

cosθ 2cosθ0

2 + cos (φ−φ0) sinθ 2sinθ0

2

)

−i

(

sin (φ−φ0) sin θ 2sinθ0

2

)¯¯¯¯¯

2

= 1

2(1 + cosθcosθ0+ sinθsinθ0cos (φ−φ0))

= 1 +n·n0

2 (2.8)

(7)

3

章 ノイズがない場合の量子状態 識別

ノイズがある場合の識別を考える前にノイズがない場合での状態の識別成功確率が ど のようなものになるのか見ていく必要がある。ここではqubit1個の場合に始ま り、qubitN個与えられたときの個別測定、集団測定での識別成功確率を求める。

3.1 qubit

1

個の場合

確率η1で状態1i=|n1iが、確率η2で状態2i=|n2iが与えられる場合を考え る。このときに正規直交基底{|mi,| −mi}でこのqubit1個に対して測定を行うと状 態は壊れるので測定後は|miまたは| −miとなる。このとき、測定後の状態が|mi なら1i| −miなら2iであると判断することとする。

or

1>=|n1>

2>=|n2>

or

|m>

|-m>

3.1: 正規直交基底{|mi,| −mi}での測定

確率η11iが与えられたときに|miを観測する確率は次のように表される。

η1|hm1i|2 =η11 +m·n1

2 (3.1)

一方、確率η22iが与えられたときに| −miを観測する確率は

η2|h −m2i|2 =η2

1 + (m·n2)

2 (3.2)

となる。したがってqubit1個の場合に識別が成功する確率は次式のように表される。

(8)

P(1) = η11 +m·n1

2 +η21 + (m·n2) 2

= 1 +m·1n1−η2n2)

2 (3.3)

ここで(3.3)式の識別成功確率を最大にするmは、mが単位ベクトルなのでシュワル

ツの不等式から次のように表すことができる。

m= η1n1−η2n2

1n1−η2n2| (3.4)

ここで|hφ12i|2

|hφ12i|2 = 1 +n1·n2

2 (3.5)

と表すことができるので、(3.5)式と(3.4)式を用いて(3.3)式を変形すると

P(1) = 1 +

11η2|hφ12i|2

2 (3.6)

となる。これがqubit1個与えられたときの最大識別成功確率となる。

3.2 qubit

N

個の場合

3.1ではqubit1個の場合での識別成功確率を求めた。では、与えられるqubit N個の場合、識別成功確率はど うなるのだろうか。

確率η1||φ1i|φ1{zi|φ1i · · ·}

N

, または確率η2||φ2i|φ2{zi|φ2i · · ·}

N

あたえられる場合を 考えてみる。まずこの場合には次2つの測定方法があげられる。

qubitを個別に測定し量子状態の識別を行う方法( 個別測定)

qubit全体を一つの系とみなして測定し量子状態の識別を行う方法( 集団測定)

この2つの測定方法では個別測定を含む測定方法である集団測定が優れている。しか し集団測定は実験的には難しい測定方法であるため個別測定での優れた識別方法が重 要となってくる。そして、その個別測定には次にあげる2つの方法がある。

(9)

qubitを個別に測定し 、多く出た方の状態であると推論する多数決を用いる方法

( 個別測定  多数決法)

qubitを個別に測定し、前の測定結果の情報からη12を変更し、あらたな{ |mi,| − mi }で測定する方法( 個別測定  適応法)

最適な識別方法を考えるためにはそれぞれの方法での識別成功確率を求める必要がある。

3.2.1

個別測定

(

多数決法

)

qubitN個与えられる場合、多数決で状態の識別を行うにはη11iが与えられ たときに識別が成功する確率、η22iが与えられたときに識別が成功する確率とそ れぞれの場合に識別が失敗する確率を求める必要がある。

まずη11iが与えられたときに識別が成功する確率は(3.1)式から次のように表 される。

η1|hm1i|2 = η11 +m·n1 2

= η1

1

2 + 12|hφ12i|2 2

14η1η2|hφ12i|2

(3.7)

一方で識別が失敗する確率は

η1|h −m1i|2 = η1

(

1 1 +m·n1 2

)

= η1

1

2 12|hφ12i|2 2

14η1η2|hφ12i|2

(3.8)

と表される。η22iが与えられる場合も同様に計算すると、識別が成功する確率は

η2|h −m2i|2 = η21 + (m·n2) 2

= η1

1

2 + 11|hφ12i|2 2

14η1η2|hφ12i|2

(3.9)

となり、識別が失敗する確率は次のようになる。

η2|hm2i|2 = η2

(

1 1 + (m·n2) 2

)

= η1

1

2 11|hφ12i|2 2

14η1η2|hφ12i|2

(3.10)

(10)

これらの式を用いてqubitの個数Nが奇数の場合を考える。

まずη11iN個与えられたときに状態の識別が成功するには測定結果で1i の個数が2iの個数より多く出る必要がある、したがって測定結果で1i2iより も多く出るときの全ての組み合わせでの識別成功確率を求めると次のように表される。

P1(N) =η1

N−1

2

k=0

NCNk

1

2 + 12|hφ12i|2 2

11η2|hφ12i|2

Nk

1

2 12|hφ12i|2 2

11η2|hφ12i|2

k

(3.11) 同様にη22iN個与えられた場合での識別成功確率を求めると

P2(N) =η2

N1

2

k=0

NCNk

1

2 + 11|hφ12i|2 2

11η2|hφ12i|2

Nk

1

2 11|hφ12i|2 2

11η2|hφ12i|2

k

(3.12) となる。したがってqubitN個の場合の識別成功確率は次式で表される。

P(N) =η1

N−1

2

k=0

NCNk

1

2 + 12|hφ12i|2 2

11η2|hφ12i|2

Nk

1

2 12|hφ12i|2 2

11η2|hφ12i|2

k

2

N1

2

k=0

NCNk

1

2 + 11|hφ12i|2 2

11η2|hφ12i|2

Nk

1

2 11|hφ12i|2 2

11η2|hφ12i|2

k

(3.13)

(3.13)式が個別測定(多数決法)による識別成功確率である。

ただし(3.13)式はη12が同確率のときに最適となる形となっている。そしてNが偶 数個の場合には測定結果で1i2iが同数出てくる組み合わせも考えられ 、多数 決で決めることはできないので確率 121i2iのど ちらかであると判断しなけ ればならない。そのため偶数個では N21まで求めていた項の数を N2 まで求め、最後の 項の組み合わせの数を 12NCN

2 とする必要がある。

3.2.2

個別測定

(

適応法

)

個別測定(適応法)での状態の識別では前の測定結果の情報を受けてη12を変更する

という方法を取る。具体的にN=2の時を見ていくと、1回目の測定結果が|miのとき 1iである確率をη10 とすると、その確率は次のように表される。

η10 = η1|hm1i|2

η1|hm1i|2+η2|hm2i|2 (3.14) ここで(3.14)式に(3.7),(3.10)式を代入し整理すると

(11)

η10 = η1

(√

11η2|hφ12i|22|hφ12i|2+ 1

)

11η2|hφ12i|2+η1−η2

(3.15)

となる。同様に1回目の測定結果が|miのときに2iである確率をη20 とし 、その確 率を求めると次式で表される。

η20 = η1|hm2i|2

η1|hm1i|2+η2|hm2i|2

= η2

(√

11η2|hφ12i|2+ 2η1|hφ12i|21

)

11η2|hφ12i|2+η1−η2

(3.16)

次に1回目の測定結果が| −miであったときを見ていくと、| −miのときに1i ある確率をη001とするとその確率は

η001 = η1|h −m1i|2

η1|h −m1i|2+η2|h −m2i|2

= η1

(√

11η2|hφ12i|2+ 2η2|hφ12i|21

)

11η2|hφ12i|2+η2−η1

(3.17)

となる。同様に1回目の測定結果が| −miのときに2iである確率をη200とし 、その 確率を求めると次式で表される。

η002 = η1|h −m2i|2

η1|h −m1i|2+η2|h −m2i|2

= η2

(√

11η2|hφ12i|21|hφ12i|2+ 1

)

11η2|hφ12i|2+η2−η1

(3.18)

次に1回目の測定結果が|miである確率をP1、1回目の測定結果が| −miである確 率をP2とするとN = 2のときの識別成功確率は次のようになる。

P(2) =P11 +

110η20 |hφ12i|2

2 +P21 +

1100η200|hφ12i|2

2 (3.19)

(12)

ここでP1,P2

P1 = η1|hm1i|2+η2|hm2i|2

= 1 2

1 + η1−η2

11η2|hφ12i|2

(3.20)

P2 = = η1|h −m1i|2+η2|h −m2i|2

= 1 2

1 + η2−η1

11η2|hφ12i|2

(3.21)

で表され 、η10η20 を計算すると

η10η02 =

2|hφ12i|2η1η2

(√

11η2|hφ12i|21−η2)1η2(1 +|hφ12i|2)+ 1

) (√

11η2|hφ12i|2+η1−η2

)2

=

2|hφ12i|2η1η2

(√

11η2|hφ12i|21−η2)1η2(1 +|hφ12i|2)+ 1

)

2

(√

11η2|hφ12i|21 −η2)1η2(1 +|hφ12i|2)+ 1

)

= η1η2|hφ12i|2 (3.22)

となる。同様にη100η200を求めると

η100η200 =

2|hφ12i|2η1η2

(√

11η2|hφ12i|22−η1)1η2(1 +|hφ12i|2)+ 1

) (√

11η2|hφ12i|2+η2−η1

)2

=

2|hφ12i|2η1η2

(√

11η2|hφ12i|22−η1)1η2(1 +|hφ12i|2)+ 1

)

2

(√

11η2|hφ12i|22−η1)1η2(1 +|hφ12i|2)+ 1

)

= η1η2|hφ12i|2 (3.23)

となり、同じ値となることが分かる。そして(3.20),(3.21),(3.22),(3.23)式を(3.19)式に 代入し整理すると次のようになる。

P(2) = 1 +

11η2|hφ12i|4

2 (3.24)

(3.24)式がN = 2のときの個別測定(適応法)による識別成功確率である。

(13)

N = 3以降も同様に求めることができるのでこの式をNで一般化すると次式で表さ れる。

P(N)= 1 +

11η2|hφ12i|2N

2 (3.25)

(3.25)式がqubitN個の場合の個別測定(適応法)による識別成功確率である。

3.2.3

集団測定

集団測定での状態の識別はqubit全体を一つの系とみなし 測定を行う。まず具体的 に個別測定(適応法)のときと同様にN = 2のときを見ていくと、2qubit系での状態は 4次元ベクトル空間で表されるが必要な成分は1i|φ1i, 2i|φ2iなので2次元空間 として考えることができる。したがってqubit1個のときと同じ形で求めることがで きる。1i = 1i|φ1i , 2i = 2i|φ2iとすると識別成功確率の式は次のように なる。

P(2) = 1 +

11η2|hψ12i|2

2 (3.26)

ここで|hψ12i|2

|hψ12i|2 = |(1|hφ1|) (2i|φ2i)|2

= |hφ12i|2|hφ12i|2

= |hφ12i|4 (3.27)

と表すことができるので、(3.26)式は

P(2) = 1 +

11η2|hφ12i|4

2 (3.28)

ノイズがある場合の識別に関する最近の研究には文献があるとなる。N = 3以降も同 様に考えることができるのでNで一般化した式は次のように表される。

P(N)= 1 +

11η2|hφ12i|2N

2 (3.29)

(3.29)式がqubitN個の場合の集団測定での識別成功確率である。そして、この式は

個別測定(適応法)での識別成功確率(3.25)式と一致することが分かる。

(14)

3.3

識別成功確率の比較

qubitN個のときの個別測定(多数決法)、個別測定(適応法)、集団測定それぞれの

測定での識別成功確率を計算により求めた結果、個別測定(適応法)と集団測定の識別 成功確率が一致することが分かった。

次にそれぞれの測定方法においてプログラム上で実際に数値を与え個別測定(多数決 法)と個別測定(適応法)の識別成功確率の比較を行うため、η1 =η2 = 0.5 ,|hφ12i|2 = 0.6とし 、Nを変化させたときの識別成功確率の値をプロットしたものを図3.2に示す。

0.75 0.8 0.85 0.9 0.95 1

0 5 10 15 20 25 30

N

adaptive, collective majority rule

success probability

3.2: qubitの個数を変化させたときの識別成功確率

3.2のグラフから、前述したように優れた測定方法である集団測定と個別測定(適 応法)の識別成功確率は一致し 、ノイズが無い場合の個別測定(多数決法)は個別測定

(適応法)よりも識別成功確率が劣るということが確認された。

(15)

4

章 ノイズがある場合の量子状態 識別

3章ではノイズがない場合の識別成功確率を求め、比較を行った。その結果、個

別測定(多数決法)が個別測定(適応法)よりも識別成功確率が劣るということが確認さ

れた。

ではノイズがある場合での識別成功確率はど うなるのだろうか。ノイズがある場合 の識別に関する最近の研究には文献[1]があるが 、本研究では準備された状態にノイズ があった場合にノイズがない場合のときと同じ 測定法を用いたらど のような影響がで るのかを調べる。したがって、この章ではノイズがない場合のときと同じ 測定法を用 いて、qubit1個の場合から考え、qubitN個のときの個別測定(多数決法)、個別

測定(適応法)での識別成功確率を求め、ノイズの影響を調べる。

4.1

ノイズがある場合の量子状態

まず、ノイズがある場合の状態|n0iは、図4.2に示すようにBloch球上で|niのま わりに等方的に分布していると仮定する。

y

x

|n>

Z

4.1: ノイズがない場合|ni

y

x

|n>

|n’>

|n’>

Z

4.2: ノイズがある場合|n0i したがってn0を平均すると方向はnと同じで長さがnよりも短くなる。よってn0の平 hn0i

hn0i=rn (0< r≤1) (4.1)

(16)

と表すことができる。r= 1はノイズが無い状態である。

4.2 qubit

1

個の場合

qubit1個の場合は確率η1で状態01i=|n01iが 、確率η2で状態02i=|n02i 与えられる場合を考える。

まず確率η101iが与えられたときに|miを観測する確率は

η1|hm01i|2 =η11 +m·n01

2 (4.2)

であり、次に確率η202iが与えられたときに| −miを観測する確率は

η2|h −m02i|2 =η21 + (m·n02)

2 (4.3)

となる。したがって識別成功確率は次のようになる。

P(1) = η1

1 +m·n01 2 +η2

1 + (m·n02) 2

= 1 +m·1n01−η2n02)

2 (4.4)

ここでmの式(3.4)を代入するとP(1)

P(1) = 1 2

1 + η11n1·n01−η2n2·n01) +η22n2·n02 −η1n1 ·n02)

11η2|hφ12i|2

(4.5)

で与えられる。ここでn01n02の平均をとると(4.1)式からhn01i=rn1,hn02i=rn2 置き換えることができる。したがって(4.5)式は

P(1) = 1 2

1 + 11n1·n1−η2n2·n1) +η22n2·n2−η1n1·n2)

11η2|hφ12i|2

= 1 +r

11η2|hφ12i|2

2 (4.6)

と表すことができる。以後の計算では置き換えた形を用いることとする。(4.6)式がノ イズが存在する場合でのqubit1個の識別成功確率である。式の形はノイズがない場合 の識別成功確率(3.6)式と似ているが分子にノイズの影響rが現れているためノイズが 大きくなる(rが小さくなる)と識別成功確率が低くなる。

(17)

4.3 qubit

N

個の場合

qubitN個与えられる場合ではノイズによる影響はど うなるのだろうか。

確率η1||φ01i|φ001{zi|φ0001 i · · ·}

N

, または確率η2||φ02i|φ002{zi|φ0002 i · · ·}

N

あたえられる場合を考 えてみる。

この場合ではノイズ影響下における個別測定(多数決法)と個別測定(適応法)の識別 成功確率を求め、比較を行う。

4.3.1

個別測定

(

多数決法

)

個別測定(多数決法)ではノイズがない場合と同様の手順で識別成功確率を求める。

まず、η101iが与えられたときに識別が成功する確率は(4.2)式から

η1|hm01i|2 = η11 +m·n01 2

= η1

1

2 + η1n1·n01−η2n2·n01 2

11η2|hφ12i|2

= η1

1

2 + r(12|hφ12i|2) 2

11η2|hφ12i|2

(4.7)

となる。この式から識別が失敗する確率を求めると次のようになる。

η1|h −m01i|2 = η1

(

1 1 +m·n01 2

)

= η1

1

2 r(12|hφ12i|2) 2

11η2|hφ12i|2

(4.8)

次にη202iが与えられたときに識別が成功する確率は(4.3)式から

η2|h −m02i|2 = η2

1 + (m·n02) 2

= η2

1

2 + η2n2·n02−η1n1·n02 2

11η2|hφ12i|2

= η2

1

2 + r(11|hφ12i|2) 2

11η2|hφ12i|2

(4.9)

である。一方、識別が失敗する確率は次式で表される。

参照

関連したドキュメント

そればかりか,チューリング機械の能力を超える現実的な計算の仕組は,今日に至るま

線遷移をおこすだけでなく、中性子を一つ放出する場合がある。この中性子が遅発中性子で ある。励起状態の Kr-87

事業セグメントごとの資本コスト(WACC)を算定するためには、BS を作成後、まず株

また適切な音量で音が聞 こえる音響設備を常設設 備として備えている なお、常設設備の効果が適 切に得られない場合、クラ

〔問4〕通勤経路が二以上ある場合

2 E-LOCA を仮定した場合でも,ECCS 系による注水流量では足りないほどの原子炉冷却材の流出が考

何日受付第何号の登記識別情報に関する証明の請求については,請求人は,請求人

つまり、p 型の語が p 型の語を修飾するという関係になっている。しかし、p 型の語同士の Merge