トフォリゲート(2021年春 Lab1)

トフォリゲート(2021年春 Lab1)#

Attention

このNotebookはQiskit v0.44の仕様に合わせてコードを改変しています。

歴史的背景#

40年前、50人の思想家からなるグループがMITのエンディコット・ハウスの芝生で写真を撮るために並びました。1981年に開催された、MITとIBMが共同で開催した計算物理に関する会議(Physics of Computation Conference)において、歴史を作っていると思った人はその時、ほとんどいませんでした。そこは、コンピューティングの物理、特に今、急成長中の分野で、教科書や大学のコースにふさわしい重要な科目である、量子コンピューティングが間違いなく誕生した場所であったのです。

この会議で、ファインマンがあの有名な言葉を述べました：「自然は古典的ではない。自然のシミュレーションをしたいなら、量子力学的にしなければならない。なんと、それは素晴らしい問題だ。なぜなら、そう簡単な問題には見えないからだ。」[1] 今月初めに、我々は、この重要な会議の40周年を祝いました。こちらでその詳細を読むことができます。

この会議で議論されたテーマのうちの一つに、可逆計算がありました。それは、MITのトマソ・トフォリとエドワード・フレドキンがその過去数年にわたりすでに考えていたことでした。[2-3] トフォリは、AND/NANDゲートの可逆バージョンを思いつきました(これは現在トフォリゲートまたは制御制御NOTゲートと呼ばれています)。NANDゲートは古典コンピューターでは普遍的であるので、トフォリゲートは、普遍的で可逆な論理ゲートになります。量子コンピューティングは、可逆コンピューティングの特別な形です：量子コンピューター上ではどんな可逆なゲートも実装でき、よってトフォリゲートも量子論理ゲートとなります。しかし、トフォリゲートのみでは、量子コンューティングにおける普遍的なゲートとはなりません。

この演習では、トフォリゲート、また、量子コンピューターにおける普遍的なゲートセットについて探究します。

参考文献#

  1. Feynman, Richard P. “Simulating physics with computers.” Int. J. Theor. Phys 21.6/7 (1982).

  2. Toffoli, Tommaso. “Reversible computing.” International colloquium on automata, languages, and programming. Springer, Berlin, Heidelberg, 1980.

  3. Fredkin, Edward, and Tommaso Toffoli. “Conservative logic.” International Journal of theoretical physics 21.3 (1982): 219-253.

古典的な論理ゲート#

古典的な計算でよく使用されるモデルの1つは、ブール論理ゲートまたは古典論理ゲートです。このゲートは、ブール関数、つまりバイナリー（0,1）の入力と出力のみを持つ関数を表します。 ブール論理の興味深い側面の1つは、少数の異なる論理ゲートの組み合わせを使用するだけで、すべての可能なバイナリー関数を形成できることです。このようなセットは、機能的に完全なセットと呼ばれます。このような有名なセットの1つは、ANDとNOTです。この2つのゲートは、すべての可能な機能を表現するのに十分です。ORとNOTについても同じことが言えます。NANDやNORなどは、単体で普遍性をもつ小さなセットですが、関数AND、NOT、ORは、古典的な計算の基本ブロックとして見なされることがよくあります。

目標

IBMの量子システムの基本ゲートセット(CX, RZ, SX, Xゲート)を使ってトフォリゲートを構築します。s! ``iv>

目的

この演習では、量子ゲートの基本的な概念を学び、量子回路の構築方法を学ぶことが目的です。

  1. Circuit Composerウィジェットを使って視覚的に学びます。

  2. Qiskitを使ってプログラミングを学びます。

すでに量子ゲートとQiskitに慣れている場合は、直接、問題 に飛ぶことができます。

# 不要な警告を取り除きます
import warnings
from matplotlib.cbook import MatplotlibDeprecationWarning
warnings.filterwarnings('ignore', category=MatplotlibDeprecationWarning)

# Qiskitの標準的なライブラリーをインポートします
from qiskit import QuantumCircuit, execute, Aer, IBMQ, QuantumRegister, ClassicalRegister
from qiskit.compiler import transpile, assemble
from qiskit.visualization import *

# piを設定しておくと便利です
import math
pi=math.pi

/tmp/ipykernel_2072/3800019171.py:3: MatplotlibDeprecationWarning: MatplotlibDeprecationWarning was deprecated in Matplotlib 3.6 and will be removed two minor releases later. Use matplotlib.MatplotlibDeprecationWarning instead.
  from matplotlib.cbook import MatplotlibDeprecationWarning

量子回路とは何でしょうか?#

量子回路は、計算が一連の量子ゲートによって行われる、量子計算におけるモデルです。量子ゲートは、しばしば、ブロッホ球における回転として表現されます。よく知られている量子ゲートをいくつか見てみましょう。

Xゲート#

Xゲートは、パウリのX行列によって表されます：

\(X = \begin{pmatrix} 0 & 1 \\ 1 & 0 \\ \end{pmatrix}\)

Xゲートは、ブロッホ球でのX軸まわりの\(\pi\)ラジアンの回転と等しいです。\(|0\rangle\)を\(|1\rangle\)へ、\(|1\rangle\)を\(|0\rangle\)へマップします。古典コンピューターにおけるNOTゲートの量子版であり、ビットフリップとも呼ばれます。

x_gate=QuantumCircuit(1) # 量子ビット1で量子回路を作ります
x_gate.x(0)
x_gate.draw(output='mpl')
_images/49df365c740768b289384c75f60085bfbe96e86c1ff52b8a3fe534acb5dac56e.png 
backend = Aer.get_backend('statevector_simulator')
result = execute(x_gate, backend).result().get_statevector()
plot_bloch_multivector(result)
_images/c2b83fc65202d73b0dd4331db44d1307b3db1d52da168d90a6b53b9aae601e70.png

SXゲート#

SXゲートは、ブロッホ球のX軸を中心とした𝜋 / 2の回転に相当します。Xゲートの平方根であることを示すためにSXゲートと呼ばれます。このゲートを2回続けて適用すると、標準のパウリXゲートが生成されます。SXの逆はSXダガーで、反対方向に\(\pi/2\)回転します。

\(SX = \frac{1}{2}\begin{pmatrix} 1+i & 1-i \\ 1-i & 1+i \\ \end{pmatrix}\)

sx_gate = QuantumCircuit(1)
sx_gate.sx(0)  
sx_gate.draw(output='mpl')
_images/6829a903bb7453760e5981864ea5e5488b27dff28c0f79f3a191870b66090875.png 
backend = Aer.get_backend('statevector_simulator')
result = execute(sx_gate, backend).result().get_statevector()
plot_bloch_multivector(result)
_images/6526500fcbd3e6fba60ddd5541a690b9413248a40700d2507cf25855e855773a.png

RZゲート#

Rzゲートは、Z軸のまわりに\(\phi\) 回転します(ここで、\(\phi\) は実数です)。行列は以下のようになります：

\(RZ = \begin{pmatrix} 1 & 0 \\ 0 & e ^{i \phi } \\ \end{pmatrix}\)

rz_gate = QuantumCircuit(1)
rz_gate.rz(pi/2, 0)
rz_gate.draw(output='mpl')
_images/1256449adf021cbf29d6e6e5aaa7467b441dcc9b993897925b676c6a91c881fe.png 
backend = Aer.get_backend('statevector_simulator')
result = execute(rz_gate, backend).result().get_statevector()
plot_bloch_multivector(result)
_images/ecd24fbd2154895add14ca8b30131419e6bf77ec76c238acea3d3853aec4f0e1.png

Z軸まわりの回転のため、デフォルトの状態\(|0\rangle\)に適用しても違いがわからないので、SXゲートを適用して生成された状態を代わりに使用し、そこにRZを適用します。

rz_gate.sx(0)
rz_gate.rz(pi/2, 0)
rz_gate.draw(output='mpl')
_images/97dd7adc85d71ef01461510acf4c47b4dcce37736c4f47550c4aa5c24eb0fde2.png 
backend = Aer.get_backend('statevector_simulator')
result = execute(rz_gate, backend).result().get_statevector()
plot_bloch_multivector(result)
_images/64c1a137db142b45bf29834a497321a22e10cab08659d9a984262bde1f3ccb64.png

Hゲート#

アダマールゲートは、\(X\)軸と\(Z\)軸の間の軸の周りに\(\pi\)回転します。 基本状態\(|0\rangle\)を\(\frac{|0\rangle + |1\rangle}{\sqrt{2}}\)にマップします。これは、測定値が1または0になる確率が等しくなり、つまり状態の「重ね合わせ」が作られることを意味します。 この状態は\(|+\rangle\)とも書かれます。 アダマールは、\(|0\rangle\) \(|1\rangle\)基底と\(|+\rangle\) \(|-\rangle\)基底を変換します。

\(H = \frac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ 1 & -1 \\ \end{pmatrix}\)

# |0>の量子ビットにHゲートをかけてみます
h_gate = QuantumCircuit(1)
h_gate.h(0)
h_gate.draw(output='mpl')
_images/7335203eb155fbece57ddb18d08036342cf6eead9099aa53b8769642c41f37b6.png 
# 結果を見てみましょう
backend = Aer.get_backend('statevector_simulator')
result = execute(h_gate, backend).result().get_statevector()
plot_bloch_multivector(result)
_images/64c1a137db142b45bf29834a497321a22e10cab08659d9a984262bde1f3ccb64.png

CXゲート(CNOTゲート)#

制御NOT(またはCNOT、またはCX)ゲートは、2量子ビットに作用します。 最初の量子ビットが\(|1\rangle\)の時のみ、2個目の量子ビットにNOT演算(Xゲートをかけるのと等しい)を適用し、それ以外の場合は何もしません。

注意：Qiskitでは、文字列のビットに、右から左へ向かって番号を付けます。

\(CX = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \\ \end{pmatrix}\)

cx_gate = QuantumCircuit(2)
cx_gate.cx(0,1)
cx_gate.draw(output='mpl')
_images/286bdfb73552dfb6de93d5162d546ed37881f8141271fddfbfa53deccebf8e1a.png

CCXゲート (トフォリゲート)#

CCXゲート(制御制御Xゲート)は、トフォリゲートとも呼ばれます。CCXゲートは3量子ビットゲートで、2個の制御ビットと1個の目標ビットがそれぞれ入力と出力として使われます。 最初の2ビットが\(|1\rangle\)の状態にある時、パウリX(またはNOT)を3番目のビットに適用します。それ以外の場合は、何もしません。

注意：Qiskitでは、文字列のビットに、右から左へ向かって番号を付けます。

\(CCX = \begin{pmatrix} 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{pmatrix}\)

ccx_gate = QuantumCircuit(3)
ccx_gate.ccx(0,1,2)
ccx_gate.draw(output='mpl')
_images/240ede2ce556deb984c5a6d688a450c1ead239d415a3913e45cc0e8e18acbe50.png

量子ゲートで古典論理ゲートを作る#

NOTゲート：#

NOTゲートはビットの値を反転し、以前述べたようにXゲートはNOTゲートとみなすことができます。NOTゲートの真理値表は以下のようになります：

入力

出力

1

0

0

1

 
not_gate=QuantumCircuit(1,1) # １量子ビットと1古典ビットで量子回路を作ります
not_gate.x(0)
not_gate.measure(0,0)
not_gate.draw(output='mpl')
_images/fda9c11e42b63caa7fee81de850a1dd3068787ce4d8b6252560fe91e3d6aa0ea.png

ANDゲート：#

ANDの出力は、入力が両方とも真であった時のみ真になります。ANDゲートの真理値表は以下のようになります：

A (入力)

B (入力)

出力

0

0

0

0

1

0

1

0

0

1

1

1

トフォリゲートを使って、2つの制御ビットを入力ビット、目標ビットを出力ビットとみなすことで、ANDゲートの解を得ることができます。

and_gate=QuantumCircuit(3,1) # 3量子ビットと1古典ビットで量子回路を作ります
and_gate.ccx(0,1,2)
and_gate.measure(2,0)
and_gate.draw(output='mpl')
_images/80598feb96971ef98a49c7ef693503777ca02aef5e3e58f10dbb1e2945cbb029.png

ORゲート：#

ORゲートは、少なくとも一つの入力が真であった時に真を返すゲートです。 真理値表は、以下のようになります：

A (入力)

B (入力)

出力

0

0

0

0

1

1

1

0

1

1

1

1

 
or_gate=QuantumCircuit(3,1) # 3量子ビットと1古典ビットで量子回路を作ります
or_gate.cx(1,2)
or_gate.cx(0,2)
or_gate.ccx(0,1,2)
or_gate.measure(2,0)
or_gate.draw(output='mpl')
_images/50c8e70f6bd3aedd14b407642ed13ce57c74ac887f6448b7ecef8707ca549bff.png

Circuit Composerウィジェットの使用#

省略

量子ゲートを組み合わせそのコストを出す#

実際の量子コンピューターは通常、すべてのゲートが物理的に実装されていることはありません。代わりに、少ないゲートでユニバーサル・ゲートセットを構成する基本ゲートのセットを使います。これは古典の場合に似て、すべての可能な演算を実装することができる命令のセットです。

このため、Qiskitは、基本ゲートセットのみを使った回路に回路を分解しなければなりません。これは、量子回路がIBM Quantumのシステムに送られる際に、普通は、Qiskit transpilerで自動的に行われます。しかし、学習のために、手を使って、基本ゲートで回路を構築することを期待されています。IBM Quantumシステムの基本ゲートは、通常、CX, ID, RZ, SX, Xゲートです。例としてibmq_mumbaiシステムを参考にすることができます。

では、以下の回路をみてみましょう：

qc = QuantumCircuit(2)
qc.sxdg(0)
qc.z(1)
qc.draw(output='mpl')
_images/20a668aff1fcb2f2485b6a3e6129da3e5e0345a06386da6d0bc5ff3c5bfeb5e0.png

次に、基本ゲートのみを使って、上記の回路を量子コンピューターのためにどのように分解するかみてみましょう。

qc = QuantumCircuit(2)
qc.sx(0)
qc.sx(0)
qc.sx(0)
qc.rz(pi,1)
qc.draw(output='mpl')
_images/e7f7244b41962daea88bc38b6776628687c0a7c6cc1f1ffc31f9e42b2b06602a.png

ご覧の通り、基本ゲートのみを使いましたが、そのためにより多くのゲートが必要になりました。

ご想像の通り、より多くのゲートを使った回路は、実行時により複雑になっていきます。そのため、回路のコストを計算したい場合、使われているゲートの数を考慮します。

しかし、すべてのゲートが等しい価値とはみなされないため、回路のコストを計算する時、以下の式を使います：

*コスト =(CXゲートの数)10 + (その他の基本ゲートの数)

アダマール#

上で述べたように、全ての演算は基本ゲートのみを用いて表現できます。ここで、例として、アダマールゲートを基本ゲートセットのみを用いて構築してみます。 \(X\)軸と\(Z\)軸の中間の軸にそって回転する基本ゲートはないため、代わりに\(X\)軸のまわりの回転と\(Z\)軸のまわりの回転を使って同じ結果を得るようにします。

どのような回転が必要だと思いますか？

q=QuantumRegister(1)
c=ClassicalRegister(1)
qc=QuantumCircuit(q,c)
qc.rz(pi/2, 0)
qc.sx(0)
qc.rz(pi/2, 0)
qc.draw(output='mpl')
_images/7030ca3eb26f118303c3ba64d473d050f979e9841bc312050098a0a400ea9c0a.png

覚えているかもしれませんが、これはさきほどRZゲートの回転を表示したときの回路です。 さきほど、\(|0\rangle\)または\(|1\rangle\)の状態にあるとき、最初のRZは何もしないことがわかりました。そのため、少し役に立たないと感じるかもしれません。 しかし、\(|+\rangle\)と\(|-\rangle\)の状態では、最初の回転は効果があります。逆のシナリオは、SXゲートを適用した後、再び\(|0\rangle\)または\(|1\rangle\)の状態になり、2回目のRZは効果がないというものです。

制御回転#

先ほど制御NOTの働きを見ましたが、\(Y\)軸まわりの制御回転をどのように構築するかという例を示しましょう。回転\(\theta\)はどの角度でもよく、\(\pi\)である必要はありません。以下は例です。

qc = QuantumCircuit(2)
theta = math.pi # thetaは何度でもよいです (任意の値としてpiを選択)
qc.ry(theta/2,1)
qc.cx(0,1)
qc.ry(-theta/2,1)
qc.cx(0,1)
qc.draw(output='mpl')
_images/db3b4ef66192379fce81d3ff68c79b487a8f92235a9adf8316702459e73821a7.png

この回路を見てみると、最初の量子ビットが0の場合、2つの回転がキャンセルされ何も起こらないことが分かります。

一方、最初の量子ビットが1の場合、\(\theta / 2\)回転を２回適用し、\(\theta\)回転となる状態を得ます。これは、\(X\)軸と\(Y\)軸が直行しているために起こります。

その他の軸まわりの回転については、他のトリックを使う必要があるでしょう。

制御制御回転#

上記の例で\(Y\)軸まわりの制御回転を作る例を確認しました。次に、(ある軸まわりに)制御回転できると仮定し、そこから制御制御回転(CCXゲートのように二つの制御量子ビットが1のときのみ作用する回転)を構築したいと思います。

qc = QuantumCircuit(3)
theta = math.pi # thetaは何でも良いです(任意にpiを選択)
qc.cp(theta/2,1,2)
qc.cx(0,1)
qc.cp(-theta/2,1,2)
qc.cx(0,1)
qc.cp(theta/2,0,2)
qc.draw()

                                           
q_0: ───────────■──────────────■───■───────
              ┌─┴─┐          ┌─┴─┐ │       
q_1: ─■───────┤ X ├─■────────┤ X ├─┼───────
      │P(π/2) └───┘ │P(-π/2) └───┘ │P(π/2) 
q_2: ─■─────────────■──────────────■───────

この回路では、1個目と2個目の量子ビットが0の場合は、何も起こりません。 2個目の量子ビットのみが1の場合は、最初に\(\pi/2\)回転が適用され、その後\(-\pi/2\)回転が適用されるのでキャンセルされます。1個目の量子ビットのみが1の場合には、2個目の量子ビットが最初のCXで1になり、\(-\pi/2\)の回転が適用され、その後\(\pi/2\)の回転が適用されるので、この2つの回転はまたお互いにキャンセルされます。

1個目の量子ビットと2個目の量子ビットが1の場合は、最初に\(\pi/2\)回転が起こり、その後2個目の量子ビットは0になるので、次の回転は適用されず、その後2個目の量子ビットは1に戻ります。その後、1個目の量子ビットが1のため、別の\(\pi/2\)回転が適用されます。よって、２回の\(\pi/2\)回転が起こり、\(\pi\)回転になります。

演習 1#

トフォリゲートの構築

先ほど、我々の基本ゲートセットでアダマールゲートをどのように構築するか見ましたが、今度は、トフォリゲートを同じように構築したいと思います。なぜトフォリゲートなのでしょうか？上記で述べたように、トフォリゲートは、古典コンピューターにおいてNANDゲートと同じように、ユニバーサル・ゲートでありますが、可逆的です。さらに、アダマールゲートとともに量子コンピューターではユニバーサル・ゲートを構築します。

基本ゲートを使ってより複雑なゲートを表現するいくつかの例を見てきました。今回は、得られた知識を使って、我々の基本ゲートのみを使ってトフォリゲートを構築します。このために、先ほどの例である、制御回転を構築し使う方法が重宝するでしょう。大きなチャレンジは、制御回転を構築するところです。

Warning

IBM Quantumシステムの基本ゲートはCX、RZ、SX、Xであり、そのほかのゲートは使えないことを覚えておいてください。

もちろん、コストの最小化にもトライしてください。(コストは CNOT の数の 10倍 プラス その他のゲート の数です。)

\[ Cost = 10 N_{CNOT} + N_{other} \]
from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit
from qiskit import IBMQ, Aer, execute


circuit = QuantumCircuit(3)

# この下にコードを書いてくだ





# この上にコードを書いてください

circuit.measure_all() # すべて量子ビットを測定します

解答例#

from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit
from qiskit import IBMQ, Aer, execute

circuit = QuantumCircuit(3)

theta = pi # 任意の値

# circuit.h(2)の置換
circuit.rz(pi/2, 2)
circuit.sx(2)
circuit.rz(pi/2, 2)
# 既に基本ゲート
circuit.cx(1,2)
# circuit.tdg(2)の置換
circuit.rz(7*pi/4, 2)
# 既に基本ゲート
circuit.cx(0,2)
# circuit.t(2)の置換
circuit.rz(pi/4, 2)
# 既に基本ゲート
circuit.cx(1,2)
# circuit.tdg(2)の置換
circuit.rz(7*pi/4, 2)
# 既に基本ゲート
circuit.cx(0,2)
# circuit.t(1)の置換
circuit.rz(pi/4, 1)
# circuit.t(2)の置換
circuit.rz(pi/4, 2)
# circuit.h(2)の置換
circuit.rz(pi/2, 2)
circuit.sx(2)
circuit.rz(pi/2, 2)
# 既に基本ゲート
circuit.cx(0,1)
# circuit.t(0)の置換
circuit.rz(pi/4, 0)
# circuit.tdg(1)の置換
circuit.rz(pi/4, 1)
# 既に基本ゲート
circuit.cx(0,1)

circuit.measure_all() # すべての量子ビットを測定します

qc = circuit

qc.draw(output='mpl')
_images/4ffface7c21613f6d0f5595d625360cc5b2e8c376a7ce18e4ead5f6c7f4131c8.png

解説

トフォリゲートの構築はいくつかの説明が可能ですが、 Qiskit Textbookにそのうちの一つである最も簡易なアダマールゲートとTゲートから構成する方法が掲載されています。

トフォリゲートのイラスト

この回路の最後に2つのSWAP操作があることがわかります。これはq0がq1と接続されていないことに起因し、一般的には実機のカップリング・マップに依存するので、SWAPを省く(CXとTダガーはq2ではなくq1にある）ことができます。

SWAPを除くと、回路は4つのゲート（アダマール、CX、TとTの共役転置）で構成されていることがわかります。Tゲートを調べると、実際にはZ軸周りに\(\pi/4\)の回転であり、\(\phi=\pi/4\)のRz(\(\phi\))ゲートとみなせます。Tダガーを構成するには、逆の回転を持つRzゲートを用います。

また、アダマールゲートは2つのRzゲートと1つのSXゲートを使って構成できます。以上、すべてのアダマールゲートとTゲート、Tダガーゲートを4つのゲートに分解すると回答のような回路が得られます。

今回はQiskit Textbookの説明を元にトフォリゲートを構成しましたが、コストを最小化するより効率的な回路は多数存在します。ぜひご自身でチャレンジしてみてください。

Additional information#

Created by: Marcel Pfaffhauser, Brahmani Thota, Junye Huang

翻訳者: 沼田祈史

Version: 1.0.1

Contents