Skip to main content
IBM Quantum Platform

量子状態の本質:隠れた変数とベルの不等式

このQiskit in Classroomsモジュールでは、以下のパッケージがインストールされた Python 環境が必要です:

  • qiskit v2.1.0 または新しい
  • qiskit-ibm-runtime v0.40.1 または新しい
  • qiskit-aer v0.17.0 または新しい
  • qiskit.visualization
  • numpy
  • pylatexenc

上記のパッケージをセットアップしてインストールするには、 Qiskitのインストールガイドをご覧ください。 実際の量子コンピュータでジョブを実行するには、 IBM Quantum® のアカウントを設定する必要があります。 IBM Cloud アカウントの設定ガイドの手順に従ってください。

このモジュールはテストされ、12秒のQPU時間を使用した。 これはあくまでも目安である。 実際の使用方法は異なる場合があります。

# Uncomment and modify this line as needed to install dependencies
#!pip install 'qiskit>=2.1.0' 'qiskit-ibm-runtime>=0.40.1' 'qiskit-aer>=0.17.0' 'numpy' 'pylatexenc'

ケイティ博士( McCormick )によるモジュールのウォークスルーを以下でご覧いただくか、 こちらをクリックして YouTube でご覧ください。



背景

量子力学の多くの計算では、系の状態を知ることから始める。 今日は、"測定前の粒子の状態について何が言えるのか?"という疑問に答えたい 測定することが許されないのであれば、どうやって知ることができるのか?

この疑問は量子力学の初期にまでさかのぼる。 この分野のパイオニアたちは、アインシュタインや他の多くの人たちと、粒子は単に測定前の未知の状態にあるだけだと言って、派閥に分かれてしまった。 また、マックス・ボルンや後のニールス・ボーアなどは、粒子の状態は、単に人間にとって未知であるだけでなく、測定以前に自然界で真に決定されていないものであると、より急進的な主張を行った。 そして、測定は確率的に粒子を明確な状態に崩壊させる。 この説明に不満だったアインシュタインが、"Gott würfelt nicht"(神はサイコロを振らない)と言ったのは有名な話だ

この見解の相違が表面化してから数十年間、多くの人がこの問題に答えが出ることはないだろう、あるいは見解の相違だと考えていた。 そして1964年、北アイルランドの物理学者ジョン・ベルが、この疑問に明確に答えることができるある実験の統計について探求した論文を書いた。 彼は、あるテストにおいて、定義された(しかし未知の)量子状態から得られる統計量と、自然界では未決定である量子状態から得られる統計量とが異なることを示した。

ベルの論文が発表された当時は、物理学の最前線にいる研究者以外には、統計の実験的検証は不可能であった。 しかし今日、 IBM クアンタムは、世界中の学生がクラウド経由で遠隔から、しかも無料で本物の量子デバイスを使用し、量子状態の本質を探求することを可能にした。 これが今日の仕事だ。

思考実験の設定:スピンの量子もつれ

スピンを持たない粒子が、それぞれスピンを持つ2つの粒子に崩壊する過程がある。 スピンは角運動量の一種であるため、角運動量保存の法則からすると、出てくる2つの粒子のスピンは正確に逆向きでなければならない。 実際、これは実験的に観察されている。

例えば、中性のπ中間子が陽電子と電子に崩壊することがある: π0e++e\pi^0\rightarrow e^+ + e^- これらの粒子が何であるか知らなくても心配はいらないし、この崩壊タイプが比較的あり得ないと知っているほどよく分かっていても心配はいらない。 ただ、出てくる粒子の一方がスピンアップしていれば、もう一方はスピンダウンしているはずで、その逆も然りである。 もちろん、"上 "と "下 "について特別なことは何もない。私たちがよく xx または yy と呼んでいるものに沿って測定が行われる場合、同じような反アライメントが観察される。この崩壊は、過去にどのような測定が行われたかについての疑問を回避することができるため、私たちにとって考慮すべき説得力のある文脈である。陽電子と電子は崩壊の瞬間まで存在すらしていなかったのだ。

私たちは、 π0\pi^0 中間子を崩壊させ、不均一磁場の影響下で、出てくる粒子の偏向を見ることができる。 スピンを偏向させるために使用される不均一場は、量子力学的スピンの存在の証拠を(偶然にも)集めるために最初に使用した研究者にちなんで、しばしばシュテルン・ゲルラッハ装置と呼ばれる。 電子と陽電子も荷電しているため(シュテルン・ゲルラッハ実験の銀原子とは異なる)、ここでの話は元の実験よりも複雑であることに注意されたい。 しかし、我々は荷電粒子が磁場の中でどのように動くかを知っているので、その影響を差し引くことができる。 以下では、計算で使われる偏向は電荷ではなく粒子のスピンによるものと仮定する。 その結果、我々の目的にとっては、どちらの観測者が陽電子を受け取り、どちらが電子を受け取るかは問題ではない。 実験のセットアップは次のようなものだ:

EPR実験のセットアップ図。 ソースは左と右に1つずつ粒子を送り出す。 それぞれシュテルン・ゲルラッハ装置、不均一磁場、測定スクリーンを通過する。

中間子が崩壊するとき、電子が一方向に、陽電子がもう一方向に追い出される。 これら2つの粒子はそれぞれ不均一な磁場中を移動し、磁場の方向に偏向するか、磁場と反対方向に偏向する。

多くの中間子の発生源があれば、その統計を取ることができる。 左側の観察者と右側の観察者(それぞれルーカスとリアーナと呼ぶ)が常に同じ軸に沿って測定する場合、これらの統計はあまり面白くないだろう:一方が測定値を上げるたびに、もう一方は測定値を下げる。 しかし、もし選手が好きな方向に沿って自由にスピンを測定することができれば、もっと面白いものが見つかるかもしれない。

粒子がスピン角運動量を持って飛び出し、それを2人の観測者が測定するという上記の実験は、当初アインシュタイン、ポドルスキー、ローゼン(EPR)がこの論文で提案したもので、これは「EPR実験」と呼ばれることもある。

選択肢

わかりやすくするために、この2つの歴史的視点をもう一度整理してみよう:

選択肢1 (アインシュタイン):2つのスピン(電子と陽電子)は、たとえそれが何であるかわからなくても、どのような軸に沿った測定結果も、自然によってあらかじめ決定されているという意味で、決定されている。 これは、スピンが空間において、私たちにはわからないが、実在する、明確に定義された方位を持っていると考えることができるだろう。 あるいは、 xxyyzz、またはその中間に沿った測定結果を決定する一連の情報や指示と考えることもできる。 陽電子のスピンを測定すると(例えばz方向に沿って)、陽電子はz方向または -z。 しかし、電子スピンは陽電子スピンと反対から始まったことが分かっているので、陽電子スピンが+zに沿って測定される場合、電子スピンは -z。 角運動量を保存する指示(スピンが反整列している)という初期条件以外には、2つのスピンの間には何のつながりもない。 このオプションは "隠れた変数 "と呼ばれることもある。つまり、異なる軸に沿った投影は決定されているが、私たちには隠されているということだ。

オプション2 (ボルン):スピンは両方とも初期状態では未定である...単に未知であるだけでなく、物理的に定義されていない。 陽電子スピンを測定することで、あらゆる可能性の空間が、+z軸または -z 軸に沿った単一の決定された状態に「折り畳まれる」。 この陽電子の測定によって、電子のスピンも、陽電子とは正反対の、zに沿った明確な投影に崩壊させられる。 この効果は、陽電子と電子の間の空間全体に広がって起こる。 これは "spooky action at a distance "と呼ばれているが、もっと大げさに "non-local physics "と呼ぶこともできるだろう。

理解度チェック

アインシュタインの説とボルンの説を実験的に区別できれば素晴らしいですね。 どちらの選択肢が正しかったとしても、同じ結果が得られるような実験にはどのようなものがあるでしょうか? この2つの選択肢で結果が異なるような実験を思いつきますか?

注: アインシュタインの案とボルンの案で異なる結果が得られるような実験を考案できれば、非常に印象的でしょう。人類がそのような実験を思いつくまでに、何十年もの歳月を要したのですから。

  • これまで説明した実験(つまり、陽電子と電子が反整列している状態で正味のスピンがない)にこだわると、 ±x\pm x±y\pm y±z\pm z に沿って両方のスピンを測定すると、角運動量保存のため、どの選択肢が正しいかとは無関係に、常に反対の符号が得られます。 ある粒子(例えば電子)のスピンをある方向(例えば +z+z )に沿って測定すると、もう1つの粒子(陽電子)のスピンは z-z に沿って測定されることになります。その代わりに、陽電子のスピンを xx の方向に沿って測定すると、 +x+xx-x のどちらかになる可能性が等しくなります。これは、隠された指示がそう言っているからか(アインシュタインのオプション1)、電子のスピンを測定した後に陽電子のスピンの確率分布が更新され、新しい確率分布が ±x\pm x (ボルンのオプション2)と半々で一致するからでしょう。 これらの点については、以下で詳しく説明する。

    spin-1 陽電子と電子のような)2つの粒子のスピンは反整列ではなく、整列していなければならない。 一方の粒子を +y+y に沿って測定した場合、もう一方の粒子を yy 軸に沿って測定すると、 +y+y のようになる。 以前と同様、これはどちらのオプションでも起こりうる。

    このレッスンの残りは、アインシュタインとボルンの選択肢を区別できる実験に費やされるので、ここではあまり詳しく説明しない。 しかし、トリックの一部は、2つの粒子を異なる方向( xxzz のような、あるいは従来のデカルト軸の間の方向)に沿って測定することである。 あとは、量子力学の予測と、隠れ変数のような古典的情報の予測から、異なる結果が得られる正確な確率を注意深く考えることから生まれる。

どちらのオプションでも、ルーカスとリアーナという2人の観測者が同じ軸に沿って測定した場合、どちらのオプションが正しいかにかかわらず、2人のスピンは反アラインになると予想される。 その理由については、以下の図をご覧いただきたい。

3つの図 まず、x軸とz軸の対角線上にある、反対方向を向いた2つのスピンベクトルの図がある。 次に、zに沿って測定が行われると、一方の友人は負のzに沿ってスピンを得、もう一方は正のzに沿ってスピンを得る。 この場合、一方は負のxに沿ってスピンを測定し、もう一方は正のxに沿ってスピンを測定する。

上の図は、アインシュタインのオプションを示したものである。 スピンの方向は逆に決まっている。 zz 軸に沿って測定すると、一方は +z+z に沿い、もう一方は z-z に沿う。陽電子は +z+z に沿い、電子は z-z に沿うと仮定する理由はない。画像は単にスピンが反対方向に測定されることを示している。 実際、アインシュタインのオプションの場合、あるスピンが最終的に測定される方向に沿ってスピンの成分を持つ必要はない。 アインシュタインのオプションの最も弱い声明は、スピンに格納されている命令のセットがあり、それがどの軸に沿って測定したときに測定結果がどうなるかを決定するというものである。 これらの命令が単純なベクトル(下図参照)の形をしていることは、説明するまでもないだろう。

3つの図 最初のものは、もはや明確に定義されたベクトルの方向を示していないが、代わりに測定結果を決定するルールのリストを持っている。 2つ目と3つ目の図は、前の写真と同じで、正負のzと正負のxに沿った測定値を示している。

下図はボルンのオプションを示したもので、陽電子と電子のスピンの方向は確率分布の中に散らばり、明確な方向はない。 分布の形を深読みしてはいけない。 それぞれのスピンは、互いに反対である限り、実際にはどの方向を向く確率もゼロではない。 ボルンのオプションの場合でも、角運動量は保存されなければならないことに注意してください。 つまり、スピンが +z+z を指すように確率の波の一方が "崩壊 "すれば、もう一方は z-z を指し、反対方向に偏向することになる。 オプションは同じに見える。

スピンベクトルの可能な方向がぼやけた図と、それに続く測定結果。

しかし、図4と図5に示すように、観測者LとRが、それぞれ120度離れた3つの軸のいずれかに沿って測定できる場合はどうなるか。 各観測者は、スピンを測定する軸(a、b、c)をランダムに決めることができる。 この2つは同じ軸に沿って測定する必要はない。 各観測者が測定するとき、自分の選んだ軸に正の投影を見つけるかもしれないし、負の投影を見つけるかもしれない。 例えば、ルーカスとリアーナは+aと -b、または+bと+cを測定するかもしれない。 たまたま同じ軸に沿って測定することを選択した場合、投影において反対の符号を得なけれ ばならないことに注意しよう:+aと -a、+bと -b、または+cと -c。 次のセクションでは、ルーカスとリアーナが測定した軸で同じ符号(++または--)になる確率と、反対の符号(+-)または(-+)になる確率を計算する方法を説明する。

2つの初期スピンベクトル、それに続く3つの線、a、b、cに沿った可能な測定値を示す図で、それぞれ正と負の方向を持つ。 これらの軸は互いに垂直ではないので、ルーカスとリアーナは、たまたま同じ軸に沿って測定することを選ばない限り、反対方向のスピンを測定することはない。

上の2つの図は、この新しい3軸測定シナリオで考えられる隠れ変数の解釈を示している。 つまり、スピンがベクトルとしてすでに決定されているか、あるいは、たとえ測定前には実験者が知ることができなくても、可能性のあるすべての測定の結果があらかじめ決定されているような物理的命令のセットが何らかの形でシステムに組み込まれているかのどちらかである。 代替案を以下に示す。 結果の確率分布が存在し、この分布からさまざまな測定結果の可能性についてある程度知ることができるが、測定に先立つ性質上、結果は未確定である。

ぼかした確率分布が、今度は3つの新しい軸a、b、cとともに、もう一度表示されている。 これには、さまざまな方向に沿って異なる測定確率が伴う。

私たちは、"2人のプレーヤーはどれくらいの頻度でスピンの投射の同じサインを見つけるべきか?"と自問することができる つまり、どの軸に沿って測定したかを記録しているわけではなく、単に同じ符号を見つけたか、違う符号を見つけたかを記録しているだけなのだ。 アインシュタインの選択肢とボルンの選択肢が、このより複雑な測定スキームで同じ結果をもたらすかどうかは明らかではない。 しかし、図4と図5から明らかなように、差があるのは possiblepossible。 アインシュタインのオプションで示されたケースの場合、 e+e+ スピンの軸 aa への投影を測定すると、 +a+a が確実に得られ、 ee- スピンの軸 bb への投影を測定すると、 b-b が(かろうじて)得られます。 しかし、ボルンのオプションでは、可能性は広く開かれている。 角運動量が保存されているのは事実だ。 しかし、2つの磁場は同じ軸を向いていないので、粒子は(磁場との相互作用によって)異なる軸に倒れなければならない状況に追い込まれる。 次のセクションでは、量子力学を使って、ルーカスとリアーナが測定された軸で同じ符号(++または--)を得る確率と、反対の符号(+-または-+)を得る確率を決定する。


予測

アインシュタインの立場(隠れた変数)は何を予測するのか?

もしアインシュタインのオプションが本当なら、 e+e+ee- の任意のペアは、それらのスピンにベクトル成分のセットを持つことになる。 例えば、電子は成分 (+a^,b^,+c^)(+\hat{a},-\hat{b}, +\hat{c}) を持つかもしれないが、その場合、陽電子は成分 (a^,+b^,c^)(-\hat{a},+\hat{b}, -\hat{c}) を持たなければならない。ここでは、各軸への投影の符号を指定しているだけで、大きさを指定しているわけではない。 このような崩壊が非常に多数( NN )起こるようにし、以下の表を作成するために測定値を集めたとしよう。

人口
パーティクル1
パーティクル2
N1N_1(+a^,+b^,+c^)(+\hat{a},+\hat{b},+\hat{c})(a^,b^,c^)(-\hat{a},-\hat{b},-\hat{c})
N2N_2(+a^,+b^,c^)(+\hat{a},+\hat{b},-\hat{c})(a^,b^,+c^)(-\hat{a},-\hat{b},+\hat{c})
N3N_3(+a^,b^,+c^)(+\hat{a},-\hat{b},+\hat{c})(a^,+b^,c^)(-\hat{a},+\hat{b},-\hat{c})
N4N_4(+a^,b^,c^)(+\hat{a},-\hat{b},-\hat{c})(a^,+b^,+c^)(-\hat{a},+\hat{b},+\hat{c})
N5N_5(a^,+b^,+c^)(-\hat{a},+\hat{b},+\hat{c})(+a^,b^,c^)(+\hat{a},-\hat{b},-\hat{c})
N6N_6(a^,+b^,c^)(-\hat{a},+\hat{b},-\hat{c})(+a^,b^,+c^)(+\hat{a},-\hat{b},+\hat{c})
N7N_7(a^,b^,+c^)(-\hat{a},-\hat{b},+\hat{c})(+a^,+b^,c^)(+\hat{a},+\hat{b},-\hat{c})
N8N_8(a^,b^,c^)(-\hat{a},-\hat{b},-\hat{c})(+a^,+b^,+c^)(+\hat{a},+\hat{b},+\hat{c})

上の表の各ケースについて、ルーカスとリアーナの軸には9つの選択肢が考えられる: aaaa ababacacbababbbbbcbccacacbcbcccc。この表から読むと、2人のオブザーバーが1行目と8行目で同じ符号を測定する確率はゼロである。 2行目から7行目まで、同じ符号を得るには4つの方法がある:

同じ符号: acac bcbc cacacbcb 反対のサイン aaaa abab, baba, bbbbcccc

つまり、アインシュタインのオプションが量子状態の正しい解釈であるならば、ルーカスとリアーナがランダムに選んだ軸に同じ符号のスピンを投影する確率は、すべての可能な母集団を合計したものになる: Psame=1iNi49(N2+N3+N4+N5+N6+N7)49 P_\text{same}=\frac{1}{\sum_i{N_i}} \frac{4}{9} (N_2+N_3+N_4+N_5+N_6+N_7)\leq \frac{4}{9} ここで等式が成り立つのは、 N1=N8=0N_1=N_8=0

理解度チェック

上のグラフの2行目では、ルーカスとリアーナの測定値が同じ符号になる可能性のあるすべての方法と、異なる符号になる可能性のあるすべての方法を列挙した。 これを3列目も繰り返す。

  • 同じ符号: abab baba bcbccbcb

    反対側のサイン aaaa acac, bbbb, cacacccc

上の表は "母集団 "を示している。つまり、もし隠れ変数の扱いが正しければ、自然が生み出す各タイプの命令の数はわからないということだ。 N1N_1 から N8N_8 の分布がどのようなものであっても、測定から同じ符号が得られる確率は常に4/9以下であることを示せ。

  • まず、 iNi=Ntot\sum_i{N_i} = N_{tot} が一定であるように、全測定試行回数が一定であると仮定することから始めよう。 N1=N8=0N_1=N_8=0 の場合、式は次のようになる

    Psame=1N2+N3+N4+N5+N6+N7×49×(N2+N3+N4+N5+N6+N7)=1Ntot×49×Ntot=49P_{same}=\frac{1}{N_2+N_3+N_4+N_5+N_6+N_7} \times \frac{4}{9} \times (N_2+N_3+N_4+N_5+N_6+N_7) = \frac{1}{N_{tot}} \times \frac{4}{9} \times N_{tot}= \frac{4}{9}

    ここで、 N10N_1 \neq 0 または N80N_8 \neq 0 のどちらかだとする

    Psame=1N1+N2+N3+N4+N5+N6+N7+N8×49×(N2+N3+N4+N5+N6+N7)=49P_{same}=\frac{1}{N'_1+N'_2+N'_3+N'_4+N'_5+N'_6+N'_7+N'_8} \times \frac{4}{9} \times (N'_2+N'_3+N'_4+N'_5+N'_6+N'_7) = \frac{4}{9}

    すべての試行の合計、 NtotN_tot は以前と同じである。 しかし、 N1N'_1 または N8N'_8 が0から増えているので、 N2N'_2 から N7N'_7 の合計は以前より少なくなっているはずである。 特に、 N2N'_2 から N7N'_7 までの合計は、 NtotN_{tot} よりも小さい

    Psame=1Ntot×49×(N2+N3+N4+N5+N6+N7)<49P_{same}=\frac{1}{N_{tot}} \times \frac{4}{9} \times (N'_2+N'_3+N'_4+N'_5+N'_6+N'_7) < \frac{4}{9}

    すべての可能なケースを組み合わせると、 Psame49P_{same} \leq \frac{4}{9} となる。

汎化

上記の処理では、特定の軸に沿った測定を考慮した。 もちろん、どの軸に沿った測定も可能だ。 つの粒子の2つのスピンベクトルを a\vec{a}b\vec{b} と呼ぶことにする。2粒子系の状態が lambdalambda のよく定義された値に対応するような隠れ変数を λ\lambda とします。 ρ(λ)\rho(\lambda)λ\lambda の確率密度とします。 最後に、スピンベクトルと隠れた変数が与えられたときに、どちらかの粒子(AまたはB)に対して行われる測定の結果があらかじめ決定されているものとして、記号 A(a,λ)A(\vec{a},\lambda)B(b,λ)B(\vec{b},\lambda) を選びます。 批判的に言えば、 AAb\vec{b} から独立しており、 BBa\vec{a} から独立していることに注意されたい。ここで、AとBの測定値間の相関に関連する質問をいくつでも投げかけることができる。 特に、次の式で与えられる期待値について尋ねることができる

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

これらの値について、 A(a,λ)1A(\vec{a},\lambda)\leq 1B(b,λ)1B(\vec{b},\lambda)\leq 1ρ(λ)\rho(\lambda) に対する正規化など、いくつかの標準的な仮定が与えられれば、2つの粒子間の相関が次の関係に従うことを示すことができる

E(a,b)E(a,d)+E(c,d)+E(c,b)2,|E(\vec{a},\vec{b})-E(\vec{a},\vec{d})|+|E(\vec{c},\vec{d})+E(\vec{c},\vec{b})|\leq 2,

ここで、 a\vec{a}b\vec{b} はシステムのスピン状態であり、 c\vec{c}d\vec{d} は参照スピン状態(システムの他の可能なスピン状態)です。 これは、現在「ベルの不等式」として知られている不等式の一種である。 ここではこの一般的な形は使わない。 その代わりに、ある特定の実験セットアップに焦点を当て、そのセットアップを量子回路にマッピングできるようにする。

ボーンの解釈(非決定論的量子力学)は何を予測するのか?

ルーカスはある軸を選び、ある粒子のスピンがプラスかマイナスのどちらかの方向にあることを見つける。 彼が何を得たにせよ、 zz の軸がその方向になるように、軸の向きを変えよう。 中間子の崩壊後、測定前の初期状態を次のように書くことができる

ψ=12(+LRL+R)|\psi \rangle =\frac{1}{\sqrt{2}}(|+\rangle_L|-\rangle_R-|-\rangle_L|+\rangle_R)

リアーナは自分の粒子のスピンを、ルーカスのスピンと相対的な角度( θ\theta )で別の方向に沿って測定する。 任意の方向 n^\hat{n} に沿ったスピン演算子は次式で与えられる

S^n^=2[cos(θ)sin(θ)eiϕsin(θ)eiϕcos(θ)]\hat{S}_{\hat{n}}=\frac{\hbar}{2}\begin{bmatrix} \cos(\theta) & \sin(\theta) e^{-i\phi} \\ \sin(\theta) e^{i\phi} & -\cos(\theta) \end{bmatrix}

この演算子の固有状態は

+n^=cos(θ/2)0+sin(θ/2)eiϕ1n^=sin(θ/2)0cos(θ/2)eiϕ1|+\rangle_{\hat{n}}=\cos(\theta/2)|0\rangle+\sin(\theta/2)e^{i\phi}|1\rangle \\ |-\rangle_{\hat{n}}=\sin(\theta/2)|0\rangle-\cos(\theta/2)e^{i\phi}|1\rangle

理解度チェック

+n^|+\rangle_{\hat{n}} が上記の演算子 S^n^\hat{S}_{\hat{n}} の固有状態であることを検証し、固有値を求めよ。

  • S^n^+n^=2[cos(θ)sin(θ)eiϕsin(θ)eiϕcos(θ)][cos(θ/2)sin(θ/2)eiϕ]\hat{S}_{\hat{n}}|+\rangle_{\hat{n}}=\frac{\hbar}{2}\begin{bmatrix} \cos(\theta) & \sin(\theta) e^{-i\phi} \\ \sin(\theta) e^{i\phi} & -\cos(\theta) \end{bmatrix} \begin{bmatrix} \cos(\theta/2) \\ \sin(\theta/2)e^{i\phi}\end{bmatrix}=2[cos(θ)cos(θ/2)+sin(θ)sin(θ/2)eiϕeiϕcos(θ/2)sin(θ)eiϕcos(θ)sin(θ/2)eiϕ]=\frac{\hbar}{2}\begin{bmatrix} \cos(\theta)\cos(\theta/2) + \sin(\theta)\sin(\theta/2)e^{i\phi} e^{-i\phi} \\ \cos(\theta/2)\sin(\theta) e^{i\phi} -\cos(\theta)\sin(\theta/2)e^{i\phi} \end{bmatrix}

    cos(θ)=cos2(θ/2)sin2(θ/2)\cos(\theta)=\cos^2(\theta/2)-\sin^2(\theta/2)sin(θ)=2cos(θ/2)sin(θ/2)\sin(\theta)=2\cos(\theta/2)\sin(\theta/2) を用いると、次のようになる

    =2[(cos(θ)+2sin2(θ/2))cos(θ/2)(2cos2(θ/2)cos2(θ/2)+sin2(θ/2))sin(θ/2)eiϕ]=\frac{\hbar}{2}\begin{bmatrix} \left(\cos(\theta) + 2\sin^2(\theta/2)\right) \cos(\theta/2) \\ \left(2\cos^2(\theta/2) -\cos^2(\theta/2)+\sin^2(\theta/2)\right)\sin(\theta/2)e^{i\phi} \end{bmatrix}=2[(cos2(θ/2)sin2(θ/2)+2sin2(θ/2))cos(θ/2)(2cos2(θ/2)cos2(θ/2)+sin2(θ/2))sin(θ/2)eiϕ]=\frac{\hbar}{2}\begin{bmatrix} \left(\cos^2(\theta/2)-\sin^2(\theta/2) + 2\sin^2(\theta/2)\right) \cos(\theta/2) \\ \left(2\cos^2(\theta/2) -\cos^2(\theta/2)+\sin^2(\theta/2)\right)\sin(\theta/2)e^{i\phi} \end{bmatrix}=2[cos(θ/2)sin(θ/2)eiϕ]=\frac{\hbar}{2}\begin{bmatrix} \cos(\theta/2) \\ \sin(\theta/2)e^{i\phi} \end{bmatrix}

    これは、 +n^|+\rangle_{\hat{n}} が固有状態であり、対応する固有値が 2\frac{\hbar}{2} であることを示している。

Lucasが選んだ軸に沿った正方向のスピンを測定する確率 +|+\rangle andand Rihannaも選んだ方向に沿った正方向のスピンを測定する確率 +n^|+\rangle_{\hat{n}} は次のようになります

P++=(L+R,n^+)ψ2P_{++}=\left|\left(_L\langle+|_{R,\hat{n}}\langle+|\right)|\psi\rangle\right|^2 P++=(L+R(cos(θ/2)++sin(θ/2)eiϕ))12(+LRL+R)2P_{++}=\left| \left(_L\langle+|_R\left(\cos(\theta/2)\langle+|+\sin(\theta/2)e^{-i\phi}\langle-|\right)\right) \frac{1}{\sqrt{2}}\left(|+\rangle_L|-\rangle_R-|-\rangle_L|+\rangle_R\right) \right|^2 P++=12(L+R(cos(θ/2)++sin(θ/2)eiϕ))(+LR)2P_{++}=\frac{1}{2}\left| \left(_L\langle+|_R\left(\cos(\theta/2)\langle+|+\sin(\theta/2)e^{-i\phi}\langle-|\right)\right) \left(|+\rangle_L|-\rangle_R\right) \right|^2 P++=12(sin(θ/2)eiϕpR)R2P_{++}=\frac{1}{2}\left| \left(\sin(\theta/2)e^{-i\phi}\vphantom{p}_R\langle-|\right) |-\rangle_R \right|^2 P++=12sin2(θ/2)P_{++}=\frac{1}{2}\sin^2(\theta/2)

理解度チェック

PP_{--} についても同様である。同じように 12sin2(θ).\frac{1}{2}\sin^2(\theta).

  • P=(LR,n^+)ψ2P_{--}=\left|\left(_L\langle-|_{R,\hat{-n}}\langle+|\right)|\psi\rangle\right|^2P=(LR(sin(θ/2)+cos(θ/2)eiϕ))12(+LRL+R)2P_{--}=\left| \left(_L\langle-|_R\left(\sin(\theta/2)\langle+|-\cos(\theta/2)e^{-i\phi}\langle-|\right)\right) \frac{1}{\sqrt{2}}\left(|+\rangle_L|-\rangle_R-|-\rangle_L|+\rangle_R\right) \right|^2P=12(LR(sin(θ/2)+cos(θ/2)eiϕ))(L+R)2P_{--}=\frac{1}{2}\left| \left(_L\langle-|_R\left(\sin(\theta/2)\langle+|-\cos(\theta/2)e^{-i\phi}\langle-|\right)\right) \left(-|-\rangle_L|+\rangle_R\right) \right|^2P=12(sin(θ/2)pR+)+R2P_{--}=\frac{1}{2}\left| \left(\sin(\theta/2) \vphantom{p}_R\langle+|\right) |+\rangle_R \right|^2P=12sin2(θ/2)P_{--}=\frac{1}{2}\sin^2(\theta/2)

これらの結果を加えると、測定された2つの軸の符号の確率は同じであることがわかる Psame=sin2(θ/2)P_{\text{same}}=\sin^2(\theta/2)

理解度チェック

この結果の計算をチェックするために何ができますか? はっきりさせておきたいのは、私たちは、それが自然界と一致しているかどうかをまだ検証してほしいとは頼んでいない。

  • (1) Pdiff=cos2(θ/2)P_{\text{diff}}=\cos^2(\theta/2) についても同様の計算を行い、確率の保存を検証する。

    (2) 既知のケースをチェックする。 θ=0\theta = 0 を挿入する。そうすると、 PsameP_{\text{same}}、2人の観測者がそれぞれ同じ軸に沿ってスピンを測定することになり、角運動量保存則に反することになる。 したがって、その確率はゼロになると予想され、実際、 θ=0\theta = 0 を挿入すると、 sin2(0/2)=0\sin^2(0/2) = 0 が得られる。

    (3) 別の既知のケースをチェックする。 θ=π\theta = \pi を試す。何を入手すべきか。 12\frac{1}{2} に気をつけて。

私たちは特に、軸が互いに対して 120deg120\deg にある場合をスケッチしていた。 ルーカスがどの方向( ±a\pm a, ±b\pm b, または ±c\pm c )を選んだとしても、それを zz と呼ぶ。そして、Rihannaは ±a\pm a, ±b\pm b, ±c\pm c のいずれかに沿って測定することをランダムに選択する。彼女の選択がLucasの選択(符号まで)と同じであれば、2人は zz に沿って測定していることになり、Rihannaが +z+z も測定する確率はゼロである。 リアーナの軸の選択はルーカスの選択とは無関係なので、これは1/3の確率で起こるはずだ。 それ以外の選択肢の場合、リアーナは zz (時間の1/3)から 120deg=2π/3120\deg = 2\pi/3 ラジアン、または zz (時間の1/3)から 240deg=4π/3240\deg = 4\pi/3 ラジアンのいずれかの軸に沿って測定することになる。 そしてもちろん、これらの軸のいずれかに沿って、スピンはプラス方向かマイナス方向かを測定することができる。 これで、ルーカスとリアーナが同じサインを得る確率の合計が出た:

Psame=13(0+sin2(π/3)+sin2(2π/3))=13(0+34+34)=12P_{\text{same}} = \frac{1}{3}\left( 0 + \sin^2(\pi/3) + \sin^2(2\pi/3) \right) = \frac{1}{3}\left( 0 + \frac{3}{4} + \frac{3}{4} \right) = \frac{1}{2}

WOW

我々はそれを示しただけだ

(Psame)max, Einstein<(Psame)max, Born.(P_\text{same})_\text{max, Einstein}<(P_\text{same})_\text{max, Born}.

一歩下がってみよう。

アインシュタインとボルンの選択肢は、測定の前に何が起こるかについての記述に違いがあるだけなので、常に同じ結果が得られるように思われた。 そして、ある軸に沿ったスピン測定の符号を事前に決定する指示があったと仮定すると、測定が同じ符号をもたらす確率の制約が得られました (Psame)Einstein49(P_{\text{same}})_\text{Einstein}\leq\frac{4}{9}。そして、量子力学のように確率分布を仮定しました。 を計算し、 (Psame)Born=12(P_{\text{same}})_\text{Born}=\frac{1}{2} に異なる値を求めた。量子力学からの予測は、隠れた変数の扱いによって許容される値よりも高い。 だから、実際に実験をして、量子力学的な状態が測定前に自然界で決定されているのか、それとも本当に起こりうる状態の確率的重ね合わせの中にあるのかを発見することができる。

この実験は、多くの場合光子であるさまざまな物理系を使って何度も行われてきた。 測定の偏り、測定のタイミング(同時性)、その他多くの微妙な考慮事項がある。 何十年もの間、こうした微妙な問題に対する懸念は着実に薄れてきた。 しかし、 IBM® 量子コンピュータを使用することで、ここで得られる答えが正しいということは、現在では広く合意されている。


実際の量子コンピューターでテストしよう!

上記の取り扱いに従い、ルーカスの測定方向を +z+z と定義しよう。これは代数的アプローチでも便利であったが、量子計算では特に便利である。一般的に測定されるのは、 zz に沿った量子ビットの射影だからである。 P++P_{++} の場合と同じ確率条件を与える量子回路を作りたい。平面の向きは自由で、 ϕ=0\phi=0、次のようになる。

P++=(L+R(cos(θ/2)++sin(θ/2)))12(+LRL+R)2P_{++}=\left| \left(_L\langle+|_R\left(\cos(\theta/2)\langle+|+\sin(\theta/2)\langle-|\right)\right) \frac{1}{\sqrt{2}}\left(|+\rangle_L|-\rangle_R-|-\rangle_L|+\rangle_R\right) \right|^2

議論の指針として、 IBM 量子コンピューターについて少し知っておく必要がある。 まず、量子ビットは 0=+z|0\rangle = |+\rangle_z の状態で初期化される。 前述したように、測定は zz。 つまり、上記の複雑な式を得るために、測定基底状態 00\langle 0|\langle 0| と量子ビットの初期状態 00|0\rangle |0\rangle の間にどのような演算子を挿入できるかを決定することが目的である。 そのためには、量子コンピューティングの基本的なゲートをおさらいする必要がある。

XX ゲート:NOT演算に相当する。 単一量子ビットゲート。

X0=1,X1=0X|0\rangle = |1\rangle,\\X|1\rangle=|0\rangle X=[0110]X=\begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix}

Qiskitでは、 XX ゲートで回路を作ると次のようになる:

from qiskit import QuantumCircuit

qc = QuantumCircuit(1)
qc.x(0)
qc.draw("mpl")

Output:

Output of the previous code cell

HH ハダマードゲート:重ね合わせ状態を作り出す。 単一量子ビットゲート。

H0=12(0+1),H|0\rangle = \frac{1}{2}\left(|0\rangle+|1\rangle\right), H1=12(01)H|1\rangle = \frac{1}{2}\left(|0\rangle-|1\rangle\right) H=12[1111]H=\frac{1}{2}\begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix}

ハダマードゲートを使った回路は次のように作られる:

qc = QuantumCircuit(1)
qc.h(0)
qc.draw("mpl")

Output:

Output of the previous code cell

CNOT Controlled-NOTゲート:このゲートは、コントロールとターゲットの2つの量子ビットを使用する。 変更されていない制御量子ビットの状態をチェックする。 しかし、もし制御量子ビットが 1|1\rangle の状態であれば、ゲートはターゲット量子ビットの状態を変化させる。もし制御量子ビットが 0|0\rangle の状態であれば、変化は全く起こらない。 以下の表記では、1番目の量子ビットをコントロール、2番目の量子ビットをターゲットとする。

CNOT00=00,CNOT01=01CNOT10=11CNOT11=10CNOT|00\rangle = |00\rangle, \\ CNOT|01\rangle = |01\rangle \\ CNOT|10\rangle = |11\rangle \\ CNOT|11\rangle = |10\rangle

CNOTゲートは、2つの量子ビットを必要とするので、回路的には少し違って見える。 これが実施方法だ:

qc = QuantumCircuit(2)
qc.cx(0, 1)
qc.draw("mpl")

Output:

Output of the previous code cell

qc.cx(0,1) に記載されている最初の量子ビットがコントロールで、2番目の量子ビットがターゲットであることに注意。 図式的に言えば、"+"マークか十字架があるのがターゲットだ。

Ry(θ)R_y(\theta) 回転Yゲート:Y軸を中心に状態を回転させる。 これは1量子ビットのゲートである。

Ry(θ)0=cos(θ/2)0+sin(θ/2)1,Ry(θ)0=sin(θ/2)0+cos(θ/2)1R_y(\theta)|0\rangle = \cos(\theta/2)|0\rangle+\sin(\theta/2)|1\rangle,\\R_y(\theta)|0\rangle = -\sin(\theta/2)|0\rangle+\cos(\theta/2)|1\rangle Ry(θ)=[cos(θ/2)sin(θ/2)sin(θ/2)cos(θ/2)]R_y(\theta)=\begin{bmatrix} \cos(\theta/2) & -\sin(\theta/2) \\ \sin(\theta/2) & \cos(\theta/2) \end{bmatrix}

最後に、回転ゲートは、ゲートのタイプ、回転量、ゲートが置かれる量子ビットを順に指定することで実装される:

import numpy as np

pi = np.pi

qc = QuantumCircuit(2)
qc.ry(pi / 2, 0)
qc.draw("mpl")

Output:

Output of the previous code cell

ゲート名 ry は、回転が発生する軸を指定する。 最初の引数 π/2\pi/2 は回転量を表し、2番目の引数はゲートが置かれる量子ビットを指定する。

理解度チェック

上記で紹介した構文、または新たに紹介した構文を使って、4つの異なるタイプの量子ゲートを含む任意の量子回路を作ってください。

  • もちろん、可能性は無限にある。 以下に例を 1 つ示します。

    
    qc=QuantumCircuit(2)
    qc.ry(pi/2,0)
    qc.cx(1,0)
    qc.x(1)
    qc.h(0)
    qc.cx(0,1)
    qc.draw("mpl")

物理実験から量子回路へ

これらのゲートの操作から、例えば、 P++P_{++}

12(+LRL+R)\frac{1}{\sqrt{2}}\left(|+\rangle_L|-\rangle_R-|-\rangle_L|+\rangle_R\right)

には、重ね合わせを得るためのハダマードゲートと、もつれを作り出すためのCNOTゲートが含まれる可能性が高い。

H、X、CNOTゲートを使って、 0L0R|0\rangle_L|0\rangle_R12(+LRL+R)\frac{1}{\sqrt{2}}\left(|+\rangle_L|-\rangle_R-|-\rangle_L|+\rangle_R\right) にする:

12(0L1R1L0R)\frac{1}{\sqrt{2}}\left(|0\rangle_L|1\rangle_R-|1\rangle_L|0\rangle_R\right) 12CNOTLR(0L1R1L1R)\frac{1}{\sqrt{2}}CNOT_{LR}\left(|0\rangle_L|1\rangle_R-|1\rangle_L|1\rangle_R\right)

ここで CNOTLRCNOT_{LR}、Lをコントロール、RをターゲットとするCNOTゲートを意味する。 これで、州のRの部分を因数分解することができる:

CNOTLR12(0L1L)1R\text{CNOT}_{LR}\frac{1}{\sqrt{2}}\left(|0\rangle_L-|1\rangle_L\right)|1\rangle_R CNOTLRHL1L1R\text{CNOT}_{LR} H_L|1\rangle_L|1\rangle_R CNOTLRHLXLXR0L0R\text{CNOT}_{LR} H_L X_L X_R|0\rangle_L|0\rangle_R

これで、量子ビットのデフォルトの開始状態で動作する量子ゲートとして、ケットを完全に記述したことになる。

ここで、 pL0R1\vphantom{p}_L\langle 0|_R\langle 1| に作用する Ry(θ)R_y(\theta) を使って、 P++P_{++} の式のブラジャーを求めることができる。

pL0R(cos(θ/2)0+sin(θ/2)1)\vphantom{p}_L\langle0|_R\left(\cos(\theta/2)\langle0|+\sin(\theta/2)\langle1|\right) pL0R(0cos(θ/2)+1sin(θ/2))\vphantom{p}_L\langle0|_R\left(|0\rangle \cos(\theta/2)+|1\rangle \sin(\theta/2)\right)^{\dagger} pL0(Ry,R(θ)0R)\vphantom{p}_L\langle0|\left(R_{y,R}(\theta)|0\rangle_R\right)^{\dagger} pL0R0Ry,R(θ)\vphantom{p}_L\langle0|_R\langle0|R_{y,R}(-\theta)

これらの結果を組み合わせると、確率 P++P_{++} を次のように書くことができる

p++=pL0R0Ry,R(θ)CNOTLRHLXLXR0L0R2p_{++}=\left|\vphantom{p}_L\langle0|_R\langle0|R_{y,R}(-\theta)\text{CNOT}_{LR} H_L X_L X_R|0\rangle_L|0\rangle_R\right|^2

これによって、量子回路の作り方が明示される。 ルーカスとリアーナが測定した粒子の量子状態を表す量子ビットに、X、H、CNOT、 RyR_y ゲートを適用し、測定を行って確率を求める。

IBM Quantumは、我々がQiskitパターンと呼ぶフレームワークを使って量子コンピューティングの問題に取り組むことを推奨している。 以下のステップで構成される。

  • ステップ1:問題を量子回路にマップする
  • ステップ2:実際の量子ハードウェア上で動作するように回路を最適化する
  • ステップ 3: ランタイム・プリミティブを使って、 IBM 量子コンピュータ上でジョブを実行する
  • ステップ4:結果の後処理

基本的に、上でやった作業はすべてステップ1だった。 出来上がった回路をQiskitを使って作ってみよう!

ステップ1: 結果を量子回路にマッピングする

# We'll begin by importing qiskit and a visualization module so that we can plot a
# histogram of our results.

from qiskit.visualization import plot_histogram

リアーナが選ぶ軸の1/3はルーカスの軸から 2π/32\pi/3 ラジアン、1/3はルーカスの軸から 4π/34 \pi/3 ラジアン、1/3は同じ軸を選ぶことを覚えておいてほしい。 だから、実際にはこれら3つのケースについて3つの量子回路を作り、その結果を合計する必要がある。 最初の1つについては丁寧に説明し、最後の2つについては簡単に述べる。

# We start by declaring our first quantum circuit, and giving it two qubits (the first "2") and two
# classical bits for storing outputs (the second "2")
# Define registers
from qiskit import ClassicalRegister, QuantumRegister

qr = QuantumRegister(2, "q")
cr = ClassicalRegister(2, "c")
qc1 = QuantumCircuit(qr, cr)

# We know from our analysis above that we need an X gate acting on each of the qubits (L and R)
qc1.x([0, 1])
# We need a Hadamard gate acting on Lucas's qubit, which we're calling the 0th qubit.
qc1.h(0)
# The controlled-NOT gate uses the 0th qubit (Lucas's) as the control and the 1st qubit (Rihanna's)
# as the target.
qc1.cx(0, 1)
# The rotation gate acts on the 1st qubit (Rihanna's) and has an argument of -2 pi/3
qc1.ry(-2 * pi / 3, 1)
# Finally, we want to measure all the qubits in the circuit to obtain measurement probabilities, and
# store the results in the classical bits.
qc1.measure([0, 1], [0, 1])
# Now we can draw the first of the three circuits that will check Bell's inequality for us.
qc1.draw(output="mpl")

Output:

Output of the previous code cell

以下のコードは、3つの回路をより合理的な方法で素早く構築する。 3つの回路の違いは、2つの量子ビットを yy

qcs = [QuantumCircuit(2, 2), QuantumCircuit(2, 2), QuantumCircuit(2, 2)]
for i in range(0, len(qcs)):
    qcs[i].x([0, 1])
    qcs[i].h(0)
    qcs[i].cx(0, 1)

qcs[0].ry(-2 * pi / 3, 1)
qcs[1].ry(-4 * pi / 3, 1)
qcs[2].ry(-2 * pi / 3, 1)
qcs[2].ry(-4 * pi / 3, 1)

for i in range(0, len(qcs)):
    qcs[i].barrier()
    qcs[i].measure([0, 1], [0, 1])

counts_list = [None] * len(qcs)
qcs[0].draw(output="mpl")

Output:

Output of the previous code cell

ここでは、「 StatevectorSampler 」と呼ばれるQiskitのリファレンスプリミティブを使用します。 サンプラーとは、システムのすべての可能な状態をサンプリングし、各状態が得られる確率(場合によっては準確率)を返すように設計されたプリミティブである。 「ショット」の数を指定し、各状態ごとの「カウント」を確認することができます。

from qiskit.primitives import StatevectorSampler

sampler = StatevectorSampler()

# Start a job that will return shots for all 100 parameter value sets.
for i in range(0, len(qcs)):
    pub = qcs[i]
    job = sampler.run([pub], shots=10000)
    # Extract the result for the 0th pub (this example only has one pub).
    result = job.result()
    data_pub = result[0].data
    counts = data_pub.c.get_counts()
    counts_list[i] = counts
#    plot_histogram(counts)

各サーキットのカウントを見ると、2つは基本的に同じで、3つ目はまったく異なっている。

plot_histogram(counts_list)

Output:

Output of the previous code cell

可能性のある結果のリストを作り、3つの回路それぞれからの各状態のカウントをすべて合計して、全体の確率を求めよう。

outcomes = ("00", "01", "10", "11")

# Here we convert "None"s into 0's so that we can sum.

for i in range(0, len(qcs)):
    for j in range(0, len(outcomes)):
        if counts_list[i].get(outcomes[j]) is None:
            counts_list[i].update({outcomes[j]: 0})

# Here we create a dictionary that holds all the outcomes and sums over their appearances
# in each of the circuits.

total_counts = {}
for i in range(0, len(outcomes)):
    total_counts[outcomes[i]] = sum(
        counts_list[j].get(outcomes[i]) for j in range(0, len(qcs))
    )

これで、各結果の合計カウントをプリントアウトし、ヒストグラムをプロットすることができる。

print(total_counts)
plot_histogram(total_counts)

Output:

{'00': 7493, '01': 7432, '10': 7605, '11': 7470}
Output of the previous code cell

理解度チェック

上記の図式は、隠れた変数や決定論によって予測される結果と矛盾しないのだろうか? それとも、確率論的な(そして非局所的な)量子力学と矛盾しないのだろうか?

  • これは確率的で非局所的な量子力学と矛盾しない。 隠れ変数処理は、同じ符号が得られる確率が4/9以下であることを予測した。 量子力学では確率は50%と予測されていた。 上記のヒストグラムは、00または11の確率を 49.97 %と表している。 これは確率論的量子力学の予測に非常に近いが、最も重要なのは、隠れ変数の扱いで許容される範囲よりも大きいということだ。

これは自然について何かを証明するものだろうか?

  • いいえ。 私たちはシミュレーターを使っていた! それは、確率論的量子力学の法則に従って動作するようにプログラムされたコンピューターである。 あるルールを提案し、そのルールに従うようにコンピュータをプログラムしたとしても、そのルールに従うコンピュータの能力は、そのルールが正しいことの証明にはならない! これを証明するには、本物の量子コンピューターを使うしかない!

ステップ2: 量子回路を実機で実行するために最適化する

当初はシミュレーターを使ってコードをデバッグしていたが、本当は実際のハードウェア上で走らせたい。 結局のところ、シミュレーターは上記の方程式に基づいて、量子力学的であるかのように装っているだけなのだ。 シミュレーターがその方程式が正しいと言ったとしても、私たちを納得させることはできないだろう。 何が起きているのかを教えてくれる本物の量子コンピューターが欲しい! そこで、使用する量子コンピューターを選択する。 必要な特性を持つ特定のデバイスを選ぶことが重要な場合もあるが、多くの場合、最も混雑していないデバイスを使いたい。

初回使用時に認証情報を保存するためのコードが以下にあります。 ノートブックを自分の環境に保存した後、必ずこの情報をノートブックから削除してください。そうすれば、ノートブックを共有するときにあなたの認証情報が誤って共有されることはありません。 詳しいガイダンスについては、 IBM Cloud アカウントの設定および信頼できない環境でのサービスの初期化を参照してください。

from qiskit_ibm_runtime import QiskitRuntimeService

# Syntax for first saving your token. Delete these lines after saving your credentials.
# QiskitRuntimeService.save_account(channel='ibm_quantum_platform',
# instance = '<YOUR_IBM_INSTANCE_CRN>', token='<YOUR-API_KEY>', overwrite=True, set_as_default=True)
# service = QiskitRuntimeService(channel='ibm_quantum_platform')

# Load saved credentials
service = QiskitRuntimeService()

backend = service.least_busy(
    operational=True, min_num_qubits=qcs[0].num_qubits, simulator=False
)
from qiskit.transpiler.preset_passmanagers import generate_preset_pass_manager

target = backend.target
pm = generate_preset_pass_manager(target=target, optimization_level=3)

qcs_isa = qcs

for i in range(0, len(qcs)):
    qcs_isa[i] = pm.run(qcs[i])
    qcs_isa[i].draw(output="mpl", idle_wires=False, style="iqp")
qcs_isa[2].draw(output="mpl", idle_wires=False, style="iqp")

Output:

Output of the previous code cell

ステップ3: ランタイムプリミティブを使用して、 IBM 量子コンピュータ上でジョブを実行する

実際の量子ハードウェア上で動作するように回路を最適化し、シミュレーターを使ってコードをデバッグした今、私たちは実際の量子コンピューターから統計データを収集し、アインシュタインとボルンの意見の相違に決着をつける準備が整った。

from qiskit_ibm_runtime import SamplerV2 as Sampler
# from qiskit_ibm_runtime import Session
# sampler.options.default_shots = 1000

# Start a job that will return shots for all 100 parameter value sets.
# The best practice is to use a session as shown below. This is available to Premium Plan, Flex
# Plan, and On-Prem (IBM Quantum Platform API) Plan users.
# result_list = [None] * len(qcs)
# real_counts_list = [None] * len(qcs)
# with Session(backend=backend) as session:
#    sampler = Sampler(mode=session)

#    for i in range(0, len(qcs)):
#        # Define the primitive unified bloc (pub)
#        pub = qcs[i]
#        job = sampler.run([pub], shots=10000)
#        # Extract the result for the 0th pub (this example only has one pub).
#        result_list[i] = job.result()
#        data_pub = result_list[i][0].data
#        counts = data_pub.c.get_counts()
#        real_counts_list[i] = counts
#    #   plot_histogram(counts)

# Open users can still carry out this experiment, but without reserving a session of use, meaning
# repeated queuing is possible.
from qiskit_ibm_runtime import Batch

batch = Batch(backend=backend)
sampler = Sampler(mode=batch)

result_list = [None] * len(qcs)
real_counts_list = [None] * len(qcs)

for i in range(0, len(qcs)):
    # Define the primitive unified bloc (pub)
    pub = qcs[i]
    job = sampler.run([pub], shots=10000)
    # Extract the result for the 0th pub (this example only has one pub).
    result_list[i] = job.result()
    data_pub = result_list[i][0].data
    counts = data_pub.c.get_counts()
    real_counts_list[i] = counts

# Close the batch because no context manager was used.
batch.close()
outcomes = ("00", "01", "10", "11")

# Here we convert "None"s into 0's so that we can sum.

for i in range(0, len(qcs)):
    for j in range(0, len(outcomes)):
        if real_counts_list[i].get(outcomes[j]) is None:
            real_counts_list[i].update({outcomes[j]: 0})

# Here we create a dictionary that holds all the outcomes and sums over their appearances
# in each of the circuits.

real_total_counts = {}
for i in range(0, len(outcomes)):
    real_total_counts[outcomes[i]] = sum(
        real_counts_list[j].get(outcomes[i]) for j in range(0, len(qcs))
    )

print(real_total_counts)
plot_histogram(real_total_counts)

Output:

{'00': 7542, '01': 7503, '10': 7304, '11': 7651}
Output of the previous code cell
# This syntax allows you to run the job on a simulator, in case you have exhausted your allotted
# time on real IBM quantum computers.
# But we strongly advise running this on real quantum computers, since this is meant to be a check
# of the behavior of real quantum systems.

# This uses a local simulator
# from qiskit_aer import AerSimulator

# This generates a simulator that mimics the real quantum system
# backend_sim = AerSimulator.from_backend(backend)

# Import an estimator, this time from qiskit (we import from Runtime for real hardware)
# from qiskit.primitives import BackendSamplerV2
# sampler = BackendSamplerV2(backend = backend_sim)

# result_list = [None] * len(qcs)
# counts_list = [None] * len(qcs)
# for i in range(0, len(qcs)):
# Define the primitive unified bloc (pub)
#    pub = qcs[i]
#    job = sampler.run([pub], shots=10000)
# Extract the result for the 0th pub (this example only has one pub).
#    result_list[i] = job.result()
#    data_pub = result_list[i][0].data
#    counts = data_pub.c.get_counts()
#    counts_list[i] = counts

# data_pubs = (result_list[0][0].data,result_list[1][0].data,result_list[2][0].data)
# outcomes = ("00", "01", "10", "11")

# Here we convert "None"s into 0's so that we can sum.

# for i in range(0, len(qcs)):
#    for j in range(0, len(outcomes)):
#        if counts_list[i].get(outcomes[j]) is None:
#            counts_list[i].update({outcomes[j]: 0})

# Here we create a dictionary that holds all the outcomes and sums over their appearances
# in each of the circuits.

# total_counts = {}
# for i in range(0, len(outcomes)):
#    total_counts[outcomes[i]] = sum(
#        counts_list[j].get(outcomes[i]) for j in range(0, len(qcs))
#    )

# print(total_counts)
# plot_histogram(total_counts)
counts_list

Output:

[None, None, None]

ステップ4:後処理と分析

一歩戻って復習してみよう:隠れ変数の扱いと3つのオフセット軸を使って、我々は測定値が同じ符号になる確率の制約を求めた Psame,hv=4/9P_{same,hv} =4/9。次に、量子力学と同じように確率分布を仮定し、その確率に異なる値を求めた。 Psame,gm=1/2P_{same,gm} = 1/2。量子力学からの予測は、隠れた変数の扱いによって許される予測よりも高い。 つまり、量子力学的な状態が測定前に自然界で決定されているのか、それとも本当に起こりうる状態の確率的重ね合わせの中にあるのかを実験的に知ることができる。

私たちは、ルーカスとリアーナがスピンの投影のいずれかの符号を測定した場合に対応する4つの可能な結果があるように量子回路を設計した:00、01、10、11である。 00と11の場合、ルーカスとリアーナは同じサインを測り、01と10の場合は反対のサインを測る。 非常によく近似すると、ルーカスとリアーナが同じサインを得る確率は約50%であり、 4/94/9 よりも確実に大きいことがわかる。これは、この確率分布を説明できる隠れ変数のセットが存在しないことを意味し、隠れ変数の扱いは実験と両立しない。

量子力学的な実験結果にはさまざまな解釈があり、実験セットアップには微妙な点が多く、たびたび見直されている。 しかし、これまでのところ、量子力学の原理と量子状態の確率論的解釈が結果を正確に説明している。 マックス・ボルンは正しかったようだ。

もう一度だけ、この意味を考えてみよう。 崩壊現象から2つの粒子が出現し、2つの粒子は異なる方向に、おそらく長い時間移動する。 スピンは明確に定義された状態にはなく、将来の測定結果を決定するための隠れ変数の指示も持っていない。 しかし、一方(例えば +z+z )の測定は、もう一方の粒子のスピンを zz の方向に沿って測定する実験の結果を必然的に決定する(それは z-z でなければならない)。このことは、ある粒子の物理学に関する何かが、もう一方の粒子(おそらくは遠くにある粒子)に対して行われることによって決定されることを意味する。 これが、人々が現実を「非局所的」と呼ぶに至った状況のひとつである。

これまで説明してきたような2つの粒子は、一方を測定することでもう一方に影響を与えることができるという意味で、何らかの形で「つながっている」。 私たちはこのような粒子を「もつれ」と呼んでいる。 エンタングルメントは単なる相関関係ではない。 例えば、磁石の北端を上にして一方に吐き出し、北端を下にしてもう一方に吐き出す古典的な機械を作ることができる。 このような磁石は完全に反相関している可能性がある。 しかし、一方を測定しても、他方には何の影響もない。 量子力学的なエンタングルメントでは、粒子Aは未定義の状態(あるいは多くの状態が混在した状態)にある可能性があり、全く別の粒子(例えばB)の測定によって、定義された状態に突き止めることができる。 クラシックの世界にはそんなものは存在しない。

これはしばしば、まったく新しい疑問や可能性をもたらす。 量子コンピュータのように、エンタングルメント(もつれ)を使って計算するようなアイデアだ! エンタングルメント(もつれ)を使って光よりも速く情報を送ろうとするなど、魅力的だが失敗に終わるものもある。 私たちは、あなたが思いつく限りの質問をし、他の人々がどのようにこれらの現象を調査したかを読むことを勧める。 量子力学の世界を探求するのは簡単ではないが、ここではいくつかのリソースを紹介しよう:

IBM 量子コース:

興味深い量子力学の論文

量子力学の教材

量子力学の教育研究


質問

指導者は、このノートがどのように使用されているかについての簡単なアンケートに答えることで、解答と一般的なカリキュラムにおける配置についてのガイダンスが付いたバージョンのノートを要求することができる。

重要な概念:

  • 量子力学は決定論的なのか確率論的なのか。
  • 隠れた変数、つまり局所的実在論は、量子力学の観測結果とは一致しない。 つまり、量子力学で観測される相関関係は、単に未知の変数によって説明できるものではない。
  • 量子力学は確率論的である。
  • エンタングルメントは現実であり、観測可能である。
  • エンタングルメントは相関関係だけではない。
  • 現実世界のシナリオを量子コンピューターにマッピングすることができる。
  • 隠れた変数とは、自然界では特定されているが、人間には未知の量のことで、この文脈では存在しない。

正誤問題:

  1. T/F アルベルト・アインシュタインは、量子力学は理論として不完全であると主張した。なぜなら、量子力学は結果の確率を記述しているだけで、その結果を決定する根本的なメカニズムを記述していないからである。
  2. T/F 「隠れた変数」とは、2つの量子力学的粒子がもつれ合う可能性があるという考え方のこと。
  3. T/F いかなる2つの相関系も量子力学的にもつれる。
  4. T/F 量子力学的なもつれは数学を正しく理解するために重要だが、実験ではそれを見ることはできない。
  5. T/F ほとんどの場合、量子力学では実験の正確な結果を知ることはできない。
  6. T/F 量子力学では、ある条件下では、粒子Aと粒子Bが接触しておらず、粒子を交換していなくても、粒子Aの状態が粒子Bの状態に影響されることがある。
  7. T/F 現実世界の実験を量子回路にマッピングすることができる。

MCの質問:

  1. spin-0 粒子が2つの spin-1/2 粒子AとBに崩壊するとする。 粒子 A に測定が行われ、そのスピンが +z+z に沿って投影されていることが明らかになる。粒子B

    • a. 間違いなく、スピンが投影されている z-z
    • b. 間違いなく、スピンが投影されている x-x
    • c. 間違いなく、スピンが投影されている y-y
    • d. どの軸を測っても、間違いなく負のスピンが投影されている。
  2. spin-0 粒子は、 spin-1/2 2つの粒子AとBに崩壊する。 粒子Aが +z+z に沿った投影を持つように測定された場合、粒子Bの測定で可能な投影はどれか 該当するものすべてに○をつけてください。

    • a. +x+x
    • b. x-x
    • c. +y+y
    • d. y-y
    • e. +z+z
    • f. z-z
  3. spin-0 粒子が2つの spin-1/2 粒子AとBに崩壊するとする。 測定前の粒子Aの状態を最もよく表すものは何か。

    • a. 粒子Aのスピンは +z+z に沿っている。
    • b. 粒子Aのスピンは z-z に沿っている。
    • c. 粒子Aのスピンは +x+x に沿っている。
    • d. 粒子Aのスピンはある方向に沿って定義されるが、他の方向には定義されない。
    • e. 粒子Aのスピンの向きは、測定前は自然界では未定である。
  4. ハダマルドゲートについて正しいものはどれか? 該当するものをすべて選択します。

    • a. H0=12(0+1)H|0\rangle = \frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)
    • b. H0=12(01)H|0\rangle = \frac{1}{\sqrt{2}}(|0\rangle-|1\rangle)
    • c. H12(01)=0H \frac{1}{\sqrt{2}}(|0\rangle-|1\rangle)=|0\rangle
    • d. H12(0+1)=0H \frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)=|0\rangle
  5. Xゲートについて正しいものはどれか? 該当するものをすべて選択します。

    • a. X0=1X|0\rangle = |1\rangle
    • b. X0=12(01)X|0\rangle = \frac{1}{\sqrt{2}}(|0\rangle-|1\rangle)
    • c. X0=0X|0\rangle = -|0\rangle
    • d. X1=0X|1\rangle = |0\rangle
    • e. X12(01)=X12(10)X\frac{1}{\sqrt{2}}(|0\rangle-|1\rangle)=X\frac{1}{\sqrt{2}}(|1\rangle-|0\rangle)
  6. 次のうち、2量子ビットゲートはどれか?

    • a. X
    • b. Ry(θ)R_y(\theta)
    • c. H
    • d. CNOT
  7. ある量子ビットが 0|0\rangle の状態にあるとする。 11\rangle の状態にあることを測定する確率は?

    • a. ノイズのないシミュレーターでは100%、実際の量子コンピューターではほぼ100%
    • b. ノイズのないシミュレータではほぼ100%、実際の量子コンピュータでは正確に100%
    • c. ノイズのないシミュレーターでは正確に0%、実際の量子コンピューターではほぼ0%
    • d. ノイズのないシミュレーターではほぼ0%、実際の量子コンピューターでは正確に0%

ディスカッションの質問:

  1. 友人A、B、Cが、ベルの不等式に関連するこの実験室の結果について話し合っている。 具体的には、彼らは、軸に沿って同じ符号を測定する量子力学的確率が、隠れ変数の扱いによって許容される確率よりも大きいことを示す画像を見ている。 (Psame)max,QM>(Psame)max,HV(P_\text{same})_\text{max,QM}>(P_\text{same})_\text{max,HV}。友人Aは "これは測定の前にスピンの状態を知らなかったことを意味する "と言っている 友人Bは「いや、それ以上だよ。 これは、測定前にスピンがすでに特定の方向を向いていないことを意味する。 でも、スピンの状態は何らかの方法で決定されたり、保存されたりするかもしれない」友人Cは言う。 つまり、未来のスピンの状態は、測定前に自然が決めたものでもないのだ 誰の意見に賛成ですか?

  2. エンタングルメントのような量子力学的現象が、現実が非局所的であることを示していることを説明しなさい。

  3. ここで得られた結果に納得するために、さらにどのような実験をしたいですか?

  4. ベルの不等式は、等間隔に配置された3つの軸 aabbcc でのみ探究できるだろうか?他の軸数でも可能か? これはどのようなものだろうか? それでも、隠れた変数が予測する確率と確率論的量子力学が予測する確率には違いがあるのだろうか?

このページは役に立ちましたか?
バグや誤字の報告、またはコンテンツの要求はGitHubで行ってください。