量子状態と状態ベクトル
量子コンピュータで利用されるデータの扱いは今のコンピュータと違って、1量子ビットでは球を使って表現されます。量子ビットは|0>状態と|1>状態を持っていてベクトルで表現できます。
∣ 0 ⟩ = [ 1 0 ] , ∣ 1 ⟩ = [ 0 1 ] \mid 0 \rangle =
\begin{bmatrix}
1 \\
0
\end{bmatrix},
\mid 1 \rangle =
\begin{bmatrix}
0 \\
1
\end{bmatrix} ∣ 0 ⟩ = [ 1 0 ] , ∣ 1 ⟩ = [ 0 1 ] また、任意の量子ビットは重ね合わせて量子状態を表現できます。
∣ ψ ⟩ = α ∣ 0 ⟩ + β ∣ 1 ⟩ = α [ 1 0 ] + β [ 0 1 ] = [ α β ] \mid \psi \rangle = \alpha \mid 0 \rangle + \beta \mid 1 \rangle =
\alpha
\begin{bmatrix}
1 \\
0
\end{bmatrix}
+ \beta
\begin{bmatrix}
0 \\
1
\end{bmatrix}
=
\begin{bmatrix}
\alpha \\
\beta
\end{bmatrix} ∣ ψ ⟩ = α ∣ 0 ⟩ + β ∣ 1 ⟩ = α [ 1 0 ] + β [ 0 1 ] = [ α β ] また、下記条件を満たします。
∣ α ∣ 2 + ∣ β ∣ 2 = 1 |\alpha|^2+|\beta|^2 = 1 ∣ α ∣ 2 + ∣ β ∣ 2 = 1 ∣ ψ ⟩ \mid \psi \rangle ∣ ψ ⟩ α \alpha α β \beta β ∣ α ∣ 2 |\alpha|^2 ∣ α ∣ 2 ∣ β ∣ 2 |\beta|^2 ∣ β ∣ 2
量子状態の極座標表示
ブロッホ球は二つの直交する純粋状態の重ね合わせで表現できる量子状態を単位球面上に表す表記法。https://ja.wikipedia.org/wiki/%E3%83%96%E3%83%AD%E3%83%83%E3%83%9B%E7%90%83
∣ ψ ⟩ = [ c o s θ − s i n θ e − i ϕ s i n θ e i ϕ c o s θ ] [ 1 0 ] = [ c o s θ s i n θ e i ϕ ] \mid \psi \rangle =
\begin{bmatrix}
cos\theta & -sin\theta e^{-i\phi}\\
sin\theta e^{i\phi} & cos\theta
\end{bmatrix}
\begin{bmatrix}
1\\0
\end{bmatrix}
=
\begin{bmatrix}
cos\theta \\ sin\theta e^{i\phi}
\end{bmatrix} ∣ ψ ⟩ = [ cos θ s in θ e i ϕ − s in θ e − i ϕ cos θ ] [ 1 0 ] = [ cos θ s in θ e i ϕ ] のように角度を利用したパラメータ表記もできます。量子データは角度で表現されるということを覚えてください。
状態ベクトルの操作と量子ゲート
状態ベクトルを操作するには、量子ゲートを利用します。変分アルゴリズムでは基本のパウリゲート、任意回転ゲート、アダマールゲートを主に利用します。
パウリゲートはXYZの３種類があり、それぞれブロッホ球でのXYZ軸周りでの180度の回転を意味し、対応する行列は、
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 = [ 0 1 1 0 ] , Y = [ 0 i − i 0 ] , Z = [ 1 0 0 − 1 ] となります。また、XYZ軸周りに任意回転にしたのが、RX,RY,RZゲートで、
R x ( θ ) = [ cos ( θ 2 ) − i sin ( θ 2 ) − i sin ( θ 2 ) cos ( θ 2 ) ] , R y ( θ ) = [ cos ( θ 2 ) − sin ( θ 2 ) sin ( θ 2 ) cos ( θ 2 ) ] , R z ( θ ) = [ e − i θ 2 0 0 e i θ 2 ] Rx(\theta) = \begin{bmatrix} \cos\left(\frac{\theta}{2}\right) &
-i\sin\left(\frac{\theta}{2}\right)\\
-i\sin\left(\frac{\theta}{2}\right) &
\cos\left(\frac{\theta}{2}\right) \end{bmatrix},
Ry(\theta) = \begin{bmatrix} \cos\left(\frac{\theta}{2}\right) &
-\sin\left(\frac{\theta}{2}\right)\\
\sin\left(\frac{\theta}{2}\right) &
\cos\left(\frac{\theta}{2}\right) \end{bmatrix},
Rz(\theta) = \begin{bmatrix} e^{-i\frac{\theta}{2}} & 0\\ 0 & e^{i\frac{\theta}{2}} \end{bmatrix} R x ( θ ) = [ cos ( 2 θ ) − i sin ( 2 θ ) − i sin ( 2 θ ) cos ( 2 θ ) ] , R y ( θ ) = [ cos ( 2 θ ) sin ( 2 θ ) − sin ( 2 θ ) cos ( 2 θ ) ] , R z ( θ ) = [ e − i 2 θ 0 0 e i 2 θ ] となります。状態ベクトルにこれらの量子ゲートを作用させることで量子状態を変化させることができます。
X [ 1 0 ] = [ 0 1 1 0 ] [ 1 0 ] = [ 0 1 ] X
\begin{bmatrix}
1 \\
0
\end{bmatrix}
=
\begin{bmatrix}
0&1 \\
1&0
\end{bmatrix}
\begin{bmatrix}
1 \\
0
\end{bmatrix}
=
\begin{bmatrix}
0 \\
1
\end{bmatrix} X [ 1 0 ] = [ 0 1 1 0 ] [ 1 0 ] = [ 0 1 ] アダマールゲートはXZ軸45度周りに回転させる便利なゲートです。
H = 1 2 [ 1 1 1 − 1 ] H=
\frac{1}{\sqrt{2}}
\begin{bmatrix}
1&1 \\
1&-1
\end{bmatrix} H = 2 1 [ 1 1 1 − 1 ]
固有値・固有ベクトル問題
ある行列に対する固有値と固有ベクトルを見つけることができれば様々な問題を解くことができます。与えられた行列H H H E 0 E_0 E 0 ∣ ψ ⟩ \mid \psi \rangle ∣ ψ ⟩
H ∣ ψ ⟩ = E 0 ∣ ψ ⟩ H \mid \psi \rangle = E_0 \mid \psi \rangle H ∣ ψ ⟩ = E 0 ∣ ψ ⟩ を満たすような、E 0 E_0 E 0 ∣ ψ ⟩ \mid \psi \rangle ∣ ψ ⟩ ∣ ψ ⟩ \mid \psi \rangle ∣ ψ ⟩
量子コンピュータの問題のほとんどがこの固有値問題に対応しています。行列Hに問題を設定し、固有値Eを探します。
量子変分原理
固有値の探し方はNISQでは主に最適化を使います。インターンでは二回目に学びました。
⟨ ψ ( θ ) ∣ H ∣ ψ ( θ ) ⟩ ≥ E 0 \langle \psi (\theta) \mid H \mid \psi (\theta) \rangle \geq E_0 ⟨ ψ ( θ ) ∣ H ∣ ψ ( θ )⟩ ≥ E 0 固有値を直接求めるのは厳しいので、間接的に固有値の期待値にできるだけ近づける最適化計算に落とし込みます。ハミルトニアンH H H ∣ ψ ( θ ) ⟩ \mid \psi (\theta) \rangle ∣ ψ ( θ )⟩
ハミルトニアンとその期待値
ハミルトニアンH H H ∣ ψ ⟩ \mid \psi\rangle ∣ ψ ⟩ H H H E E E ⟨ ψ ∣ \langle \psi \mid ⟨ ψ ∣
⟨ ψ ∣ H ∣ ψ ⟩ = ⟨ ψ ∣ E ∣ ψ ⟩ ⟨ ψ ∣ H ∣ ψ ⟩ = E ⟨ ψ ∣ ψ ⟩ ⟨ ψ ∣ H ∣ ψ ⟩ = E \langle \psi \mid H \mid \psi \rangle = \langle \psi \mid E \mid \psi \rangle\\
\langle \psi \mid H \mid \psi \rangle = E \langle \psi \mid \psi \rangle\\
\langle \psi \mid H \mid \psi \rangle = E ⟨ ψ ∣ H ∣ ψ ⟩ = ⟨ ψ ∣ E ∣ ψ ⟩ ⟨ ψ ∣ H ∣ ψ ⟩ = E ⟨ ψ ∣ ψ ⟩ ⟨ ψ ∣ H ∣ ψ ⟩ = E のように計算できます。1量子ビットの例として、任意の量子状態∣ ψ ⟩ \mid \psi \rangle ∣ ψ ⟩ H H H
⟨ ψ ∣ H ∣ ψ ⟩ = [ α ∗ β ∗ ] [ a b c d ] [ α β ] \langle \psi \mid H \mid \psi \rangle =
\begin{bmatrix}
\alpha^* & \beta^*
\end{bmatrix}
\begin{bmatrix}
a&b\\
c&d
\end{bmatrix}
\begin{bmatrix}
\alpha\\
\beta
\end{bmatrix} ⟨ ψ ∣ H ∣ ψ ⟩ = [ α ∗ β ∗ ] [ a c b d ] [ α β ] のように計算できます。
ハミルトニアンZの期待値を求める
任意の状態ベクトルにおけるハミルトニアンがH = Z H = Z H = Z
⟨ ψ ∣ H ∣ ψ ⟩ = [ α ∗ β ∗ ] [ 1 0 0 − 1 ] [ α β ] = ∣ α ∣ 2 − ∣ β ∣ 2 \langle \psi \mid H \mid \psi \rangle =
\begin{bmatrix}
\alpha^* & \beta^*
\end{bmatrix}
\begin{bmatrix}
1&0\\
0&-1
\end{bmatrix}
\begin{bmatrix}
\alpha\\
\beta
\end{bmatrix}
= |\alpha|^2 - |\beta|^2 ⟨ ψ ∣ H ∣ ψ ⟩ = [ α ∗ β ∗ ] [ 1 0 0 − 1 ] [ α β ] = ∣ α ∣ 2 − ∣ β ∣ 2 となります。∣ α ∣ 2 |\alpha|^2 ∣ α ∣ 2 ∣ β ∣ 2 |\beta|^2 ∣ β ∣ 2
ハミルトニアンがH = X H = X H = X H H H
⟨ ψ ∣ H ∣ ψ ⟩ = [ α ∗ β ∗ ] [ 0 1 1 0 ] [ α β ] = α ∗ β + α β ∗ \langle \psi \mid H \mid \psi \rangle =
\begin{bmatrix}
\alpha^* & \beta^*
\end{bmatrix}
\begin{bmatrix}
0&1\\
1&0
\end{bmatrix}
\begin{bmatrix}
\alpha\\
\beta
\end{bmatrix}
= \alpha^* \beta + \alpha \beta^* ⟨ ψ ∣ H ∣ ψ ⟩ = [ α ∗ β ∗ ] [ 0 1 1 0 ] [ α β ] = α ∗ β + α β ∗ となりますが、通常の測定ではこの値を求めることができませんので、ハミルトニアンの変形を通じて解を得ます。
ハミルトニアンの変形
ハミルトニアンH H H Z Z Z
より、
⟨ ψ ∣ X ∣ ψ ⟩ = ⟨ ψ ∣ H Z H ∣ ψ ⟩ = ⟨ ψ ′ ∣ Z ∣ ψ ′ ⟩ \langle \psi \mid X \mid \psi \rangle \\
= \langle \psi \mid HZH \mid \psi \rangle\\
= \langle \psi' \mid Z \mid \psi' \rangle ⟨ ψ ∣ X ∣ ψ ⟩ = ⟨ ψ ∣ H Z H ∣ ψ ⟩ = ⟨ ψ ′ ∣ Z ∣ ψ ′ ⟩ のように状態ベクトルを変形することで対応できます。Yの場合は、
Y = R X ( − π / 2 ) Z R X ( π / 2 ) Y = RX(-\pi/2) Z RX(\pi/2) Y = RX ( − π /2 ) ZRX ( π /2 ) と、なり同様に、
⟨ ψ ∣ Y ∣ ψ ⟩ = ⟨ ψ ∣ R X ( − π / 2 ) Z R X ( π / 2 ) ∣ ψ ⟩ = ⟨ ψ ′ ′ ∣ Z ∣ ψ ′ ′ ⟩ \langle \psi \mid Y \mid \psi \rangle \\
= \langle \psi \mid RX(-\pi/2) Z RX(\pi/2) \mid \psi \rangle\\
= \langle \psi'' \mid Z \mid \psi'' \rangle ⟨ ψ ∣ Y ∣ ψ ⟩ = ⟨ ψ ∣ RX ( − π /2 ) ZRX ( π /2 ) ∣ ψ ⟩ = ⟨ ψ ′′ ∣ Z ∣ ψ ′′ ⟩ とすることができます。これらは、量子回路の測定の直前に上記のHやRXゲートを入れて測定を行うことで期待値を計算することができます。
線形結合
実際にはハミルトニアンH H H
⟨ ψ ∣ a H 1 + b H 2 ∣ ψ ⟩ = ⟨ ψ ∣ a H 1 ∣ ψ ⟩ + ⟨ ψ ∣ b H 2 ∣ ψ ⟩ = a ⟨ ψ ∣ H 1 ∣ ψ ⟩ + b ⟨ ψ ∣ H 2 ∣ ψ ⟩ \langle \psi \mid aH_1 + bH_2 \mid \psi \rangle \\ = \langle \psi \mid aH_1 \mid \psi \rangle + \langle \psi \mid bH_2 \mid \psi \rangle \\ = a\langle \psi \mid H_1 \mid \psi \rangle + b\langle \psi \mid H_2 \mid \psi \rangle ⟨ ψ ∣ a H 1 + b H 2 ∣ ψ ⟩ = ⟨ ψ ∣ a H 1 ∣ ψ ⟩ + ⟨ ψ ∣ b H 2 ∣ ψ ⟩ = a ⟨ ψ ∣ H 1 ∣ ψ ⟩ + b ⟨ ψ ∣ H 2 ∣ ψ ⟩ これを使えば、複雑なハミルトニアンの場合でも、
H = 1.2 X 0 Z 2 + 2.5 Z 0 X 1 Y 2 − 3.4 Z 2 X 1 H = 1.2 X_0 Z_2 + 2.5 Z_0 X_1 Y_2 - 3.4 Z_2 X_1 H = 1.2 X 0 Z 2 + 2.5 Z 0 X 1 Y 2 − 3.4 Z 2 X 1 期待値は、
⟨ ψ ∣ 1.2 X 0 Z 2 + 2.5 Z 0 X 1 Y 2 − 3.4 Z 2 X 1 ∣ ψ ⟩ = 1.2 ⟨ ψ ∣ X 0 Z 2 ∣ ψ ⟩ + 2.5 ⟨ ψ ∣ Z 0 X 1 Y 2 ∣ ψ ⟩ − 3.4 ⟨ ψ ∣ Z 2 X 1 ∣ ψ ⟩ \langle \psi \mid 1.2 X_0 Z_2 + 2.5 Z_0 X_1 Y_2 - 3.4 Z_2 X_1 \mid \psi \rangle\\
= 1.2\langle \psi \mid X_0 Z_2 \mid \psi \rangle + 2.5 \langle \psi \mid Z_0 X_1 Y_2\mid \psi \rangle - 3.4 \langle \psi \mid Z_2 X_1 \mid \psi \rangle ⟨ ψ ∣ 1.2 X 0 Z 2 + 2.5 Z 0 X 1 Y 2 − 3.4 Z 2 X 1 ∣ ψ ⟩ = 1.2 ⟨ ψ ∣ X 0 Z 2 ∣ ψ ⟩ + 2.5 ⟨ ψ ∣ Z 0 X 1 Y 2 ∣ ψ ⟩ − 3.4 ⟨ ψ ∣ Z 2 X 1 ∣ ψ ⟩ のように個別に計算をして和を求めることで計算できます。
量子断熱計算
量子断熱計算は、状態ベクトルを断続的に変化させることで基底状態をキープしたまま時間発展を行うことができる理論です。
初期状態のハミルトニアンをH s t a r t H_{start} H s t a r t H f i n a l H_{final} H f ina l t t t T T T
H t e m p = ( 1 − t T ) H s t a r t + t T H f i n a l H_{temp} = (1-\frac{t}{T})H_{start} + \frac{t}{T}H_{final} H t e m p = ( 1 − T t ) H s t a r t + T t H f ina l
としたときにT → ∞ T\rightarrow\infty T → ∞ λ \lambda λ
H t e m p ∣ ψ ⟩ = E 0 t e m p ∣ ψ ⟩ H_{temp}\mid \psi \rangle = E_{0temp}\mid \psi \rangle H t e m p ∣ ψ ⟩ = E 0 t e m p ∣ ψ ⟩ 時間発展計算は、
∣ ψ t + 1 ⟩ = U ∣ ψ t ⟩ = e − i H t ∣ ψ t ⟩ \mid \psi_{t+1} \rangle = U \mid \psi_t \rangle = e^{-iHt} \mid \psi_t \rangle ∣ ψ t + 1 ⟩ = U ∣ ψ t ⟩ = e − i H t ∣ ψ t ⟩ となります。課題は基底状態と第一励起状態が最接近する部分ですが、E 1 ( t ) − E 0 ( t ) E_1(t)-E_0(t) E 1 ( t ) − E 0 ( t )
QAOA回路
QAOAは上記の量子断熱計算の時間発展計算をansatzとして変分アルゴリズムに適用したものです。
|0> --H-- --RZZ--RZ-- --RX-- --RZZ--RZ-- --RX--
| |
|0> --H-- --RZZ--RZ-- --RX-- --RZZ--RZ-- --RX--
全体構成を大きく分けて、
1.mixerを選ぶと、初期状態が決まる
です。2番目の定式化は自動的に時間発展されます。
一番左は初期固有状態|+>を作っています。これに対応するハミルトニアンXは時間発展の形でRxで出てきます。通常はmixerはデフォルトではXを選び、初期状態はそれに対応する|+>を作るために、H ⊗ N H^ {\otimes N} H ⊗ N
また、CX-Rz-CXは問題設定のハミルトニアンのウェイトに対応し、次のRzはハミルトニアンのバイアスに対応しています。
blueqatのQaoaAnsatzはハミルトニアンとステップ分割数を決めれば、自動的にこの時間発展のansatzを構成してくれます。ライブラリ側で適切な計算をしてくれているので、私たちがするのはハミルトニアンの定式化とmixerの選択です。簡単な定式化を解いてみましょう。
## イジング定式化
イジングモデルと呼ばれる物理モデルを利用して組合せ最適化問題を解きます。イジングモデルでは通常スピンの上下を利用して問題を定式化します。定式化では演算子としてZを利用します。Zは+1と-1を期待値としてもとめられます。
H = − Z 0 − Z 0 ∗ Z 1 H = -Z_0 -Z_0*Z_1 H = − Z 0 − Z 0 ∗ Z 1 ## QUBO
つぎに、上記ハミルトニアンZの定式化から派生して、0と1で定式化をするQUBOを確認します。QUBOはZの定式化を書き換えた形になりますが、blueqatでは自動変換してくれるツールがあります。
定式化は、
H = 2 q 1 − 4 q 0 q 1 H = 2q_1-4q_0q_1 H = 2 q 1 − 4 q 0 q 1 としてみたいと思います。
Copy
from blueqat import Circuit
from blueqat.utils import qaoa
from blueqat.pauli import qubo_bit as q
from blueqat.pauli import X,Y,Z,I
hamiltonian = 2*q(1)-4*q(0)*q(1)
step = 1
result = qaoa(hamiltonian, step)
result.circuit.run(shots=100)Counter({'11': 47, '00': 36, '10': 17})
