★想定読者は,理系学部の授業を受けた学生以上です.
こちらの本をもとに説明していきます.
この本は量子アニーリングの提唱者である西森秀稔さんによって書かれており,過不足なく,かなり分かりやすくまとめてあります.
大学の基礎物理を学んだ人には理解できるようになっていると思います.
概要:『量子アニーリング』は最適化手法のひとつ
量子アニーリングとは何かを説明する前に,最適化という言葉を説明しておく必要があります.
最適化問題とは,制約の下で何かを最大もしくは最小にしたいとき,その解を見つけることです.
最大,最小にしたい関数を目的関数や評価関数と呼びます.
つまり,下の図のようなグラフが最小になる矢印のところを見つける問題です.

縦軸が目的関数,横軸が解の候補です.
特に,変数が離散的である場合を組合せ最適化問題と言います.すなわち,上の図において解はとびとびの値のどれかであり,それを探す問題が組合せ最適化問題です.
組合せ最適化問題には,いくつかの点を訪れて戻ってくる巡回セールスマン問題,箱にできるだけ良いものを詰めるナップザック問題などがあります.
組合せ最適化は,多くの場合,問題の規模が大きくなると解の候補が莫大に増え,一つ一つ解を見ていくというのはコンピュータでも実質的に不可能となります.
そのため,一つ一つの解を見るのではなく,実質的な時間で「それなりに」良い解を求める手法がたくさん提案されており,
量子アニーリングも,この組合せ最適化問題を解くための手法のひとつです.
後で述べるように,組合せ最適化問題の解はスピンという物理学の概念で置き換えて考えることができます.
そして,物理学の自然な法則をまねてスピンを従わせる(アルゴリズムを組む)ことで,スピンの状態(解)を計算しようという考えに基づいた最適化手法として,アニーリング法が提案されています.
低エネルギー状態に落ち着いた時のスピンの状態を最適解とするのがアニーリング法です.
最適解の候補 ⇔ スピン
目的関数 ⇔ エネルギー(ハミルトニアン)
計算アルゴリズム ⇔ 物理法則
すでに良く知られているシミュレテッドアニーリングは温度を下げていって,熱ゆらぎを持つエネルギーが低くなったときの状態を解とする,古典的な物理学を利用した方法です.

赤の線は解の動きを表しています.
一方,量子アニーリングは磁場を下げていって,量子ゆらぎを持つエネルギーが小さくなったときの状態を解とする,量子力学の効果を利用した方法です.

青の部分は解の存在確率を表しています.
ハードとソフトの両面からの研究がある
量子アニーリングに関して述べるには,ハードウェアとソフトウェアの両方から議論する必要があります.
上で述べた話はただ単に物理現象を模した計算方法ではなく,実際にスピンと磁場を利用した量子アニーリングマシン(量子コンピュータの一種)を用いて計算されるからです.
量子アニーリングマシンを用いて計算することで,古典的なコンピュータ(いま我々が使っている,0と1ですべてが記述されるコンピュータ)よりも高速に解くことができるとされています(これについての議論は,他の記事でします).
※量子アニーリングマシンは組合せ最適化のみ解けることに注意.連続最適化や他の計算については量子ゲート方式の量子コンピュータを使用する.
量子アニーリング式のマシンは,D-Wave社によって開発されており,利用することもできます.
一方で,量子アニーリングのアルゴリズム自体にも議論の余地があるとされています.
古典コンピュータで量子アニーリング計算を行うと,シミュレテッドアニーリング(古典物理的なアルゴリズム)よりも高速に解ける場合があるからです.
つまり,ハードウェアとして古典コンピュータを用いてシミュレーションするとしても,アルゴリズム(量子効果を模した計算方法)自体に研究の余地があるとされています.
量子アニーリングのアルゴリズム
いよいよ,量子アニーリングのアルゴリズムについて説明していきます.
まずざっくりと
以降で説明する内容を,ざっくりと説明します.
量子アニーリングは,以下の図のように,はじめは横磁場(外部磁場)によって平坦なポテンシャル(エネルギー)をつくり,そこにどの解も等しい確率で存在するとします.

そして,外部磁場をだんだんと弱めていき,もとのポテンシャルに近づけていき,自然と最適解の確率が最大となるという仕組みです.
この解の動きはどのように計算するのかというと,シュレディンガー方程式に従って発展させます.つまり,各計算ステップにおけるポテンシャルに対して,固有ベクトルを求めます(詳しくはあとで説明します).
最終的に外部磁場が0になったときのポテンシャルに対する固有ベクトルが最適解となります.
(準備)量子力学を簡単に
量子アニーリングの理解に必要な量子力学にの知識について少し解説しておきます.
重ね合わせ
量子力学では,ミクロな物の状態(例えば,電子の位置や電子スピンの向き)ははっきりと観測できないために,状態を確率的に扱います.つまり,「電子は右のあたりに0.4の確率で存在して,左のあたりに0.6の確率で存在する」といったように扱います.
そして,その状態は次のように重ね合わせで表します.
|電子の状態>=
この|電子の状態>のことを波動関数と呼び,電子の存在確率密度の片割れ(片割れなので係数は平方根になっている)を意味しますが,量子力学を学んだことがないひとは,とりあえず電子の状態を表すものなんだなと思っておけばOKです.
電子や原子にはスピンと呼ばれるものがあり,それもミクロな状態のひとつなので,同様に状態の重ね合わせで表現します.
上では電子の位置について述べましたが,位置は固定されているとして,スピンの向きだけが分からないと仮定します.
例えば,一つの電子を考えましょう.電子スピンの候補として,上向きスピンと下向きスピンの2つの状態があります.電子がどちらのスピンをもっているか分からず,上向きスピンと下向きスピンが等しい確率で存在する場合は,
|電子> =
とします.
スピンを行列で表現する
上で述べたや
はケットベクトルと呼ばれ,各々は
です.一方,や
のような書き方はブラベクトルと呼ばれ,ケットベクトルの複素共役をとって転置したものです.
これらスピンは,パウリ行列と呼ばれるスピン演算子の固有ベクトルとなっています.
z方向(上下方向)のパウリ行列は
であり,これに上下スピンを作用させると固有値は
のように,1とー1になります.
後で述べるように,ほとんどの組合せ最適化は1とー1をとる変数で記述できるので,これら上下スピンを組合せ最適化を定式化することができます.
とにかくパウリ行列にスピンを作用させたら,1もくしは-1の固有値が出るとだけ理解しておけばOKです.
シュレディンガー方程式
量子力学の基礎方程式として,シュレディンガー方程式があります.
シュレディンガー方程式は,以下のように書かれます.
ただし,はハミルトニアン,
はエネルギー演算子です.どちらも演算子で,波動関数
を作用させて初めて意味を持ちます.
このシュレディンガー方程式は,自然が従う方程式なので,実際に世に存在する電子や原子などのミクロなものはこれに従って存在しています.
この方程式により,波動関数(量子状態)がハミルトニアンに対してどのように時間発展するのかが分かります.
上のシュレディンガー方程式は時間に依存しますが,時間に依存しない場合は,以下のように書けます.
ほとんど変わらないように見えますが,エネルギー演算子が定数
になっています.この式はハミルトニアン
に対して,
を固有状態,
を固有値とした固有方程式となっていることがわかります.
この時間依存しないシュレディンガー方程式は単にシュレディンガー方程式と呼ばれることが多いです.
量子アニーリングでは,時間依存しないシュレディンガー方程式が常に成り立つ程度にゆっくりと外部磁場を弱め,ハミルトニアンを変化させます.
つまり,量子アニーリングの過程では上の固有方程式(時間依存しないシュレディンガー方程式)を解き,固有値,固有状態を求めればよいことになります.
それでは,いよいよ量子アニーリングの流れについて説明していきます.
とはいっても,おおまかな流れはすでに上で述べました.
解きたい組合せ最適化を量子アニーリングで解ける形に直し(Step1),外部磁場が加わったハミルトニアンに対するシュレディンガー方程式を,外部磁場を弱めながら解いていき,外部磁場が0になったときの固有状態を最適解とする(Step2)流れです.
量子アニーリング Step1.問題をイジング形にする
量子アニーリングを実行するには,解きたい最適化問題の目的関数を,イジング模型と呼ばれる,スピンが並んだもののハミルトニアンの形にします.
以下は,イジング模型の説明です.
イジング模型の説明
一番簡単なイジング模型のハミルトニアンは,以下の形で表します.
はスピン
とスピン
の相互作用の強さを表します.
はスピン
のパウリ行列であり,先ほど述べたように,波動関数を作用させると固有値1もしくは-1を吐き出します.
は具体的には,
です.
は2×2の単位で,
は直積です.直積の計算方法については,こちらをご覧ください.https://atsblog.org/tyokuseki-matrix/
とはいえ,直積のような,難しい行列計算をわざわざする必要はありません.は単にスピン
にだけ作用する演算子です.
例えば,2つのスピンがあったとして,それを(1つ目のスピンが上向き,2つ目のスピンが下向き)のように表すと,
(スピン1だけに作用して,固有値1を吐き出した)
(スピン2だけに作用して,固有値ー1を吐き出した)
となるだけです.
ハミルトニアンを見ればわかるように,2個のスピン
に対して
のときは,2個のスピンは同じ向きを向いたほうがエネルギーが小さくなります.反対に,
のときは,反対向きを向いたほうがエネルギーが小さくなります.
具体例
ここで,具体例で見てみましょう.スピンが3つである場合を考えます.パターンとしては,2×2×2=8通りあり,この8通りからハミルトニアンが最小になるものを見つけます.

各スピンの相互作用は,
とします.
すると,この場合,ハミルトニアンは
となります.では,各スピンを作用せます.例えば,を作用させると,
となり,固有値(イジング模型のエネルギーに相当するもの)は8になります.
同様に,すべてのスピン状態について調べると,
となり,最小の固有値-6をとるの状態が最適解になります.
このように,3つの変数(スピン)だけでも,最適解を探すのは大変ですね(対称性を利用すればもう少し計算量が減りますが).
このような全てのスピン状態を書き出す大変な作業をしなくても,シュレディンガー方程式という自然法則を応用して最適解を見つけるというのが,量子アニーリングです.
ここでは述べたのはスピン同士の2体相互作用のみを考慮したハミルトニアンですが,ハミルトニアンをパウリ行列で表すことさえできれば,他の形のハミルトニアンでも量子アニーリングを行うことができます.
さて,解きたい組合せ最適化問題の目的関数を,上で説明したイジング模型のハミルトニアンの形にします.
ここでは例として,代表的な組合せ最適化問題である巡回セールスマン問題(TSP)を扱いましょう.
巡回セールスマン問題をイジング模型に書き直す
巡回セールスマン問題は,複数の都市を訪れ,同じ場所に戻ってくる問題です.

この組合せの数は都市の数に対して
通りあるので,しらみつぶしに調べていくのは困難なため,アニーリング法などの最適化手法を用いて計算します.
目的関数は,その総移動距離であり,都市と都市
の間の距離を
とすると,
です.ただし,は都市
と都市
に行く選択をしたなら1,行かない選択をしたなら0をとる変数です.つまり,この
のパターンの中から,目的関数が最小になるものを探すことが目標となります.
上の目的関数は,以下のようにも書くことができます.
は都市を周る順番を表しており,
は
番目に都市
を訪れたなら1,そうでないなら0です.このように書くことで,
と
がどちらも1をとるときにのみ
が目的関数に加算されます.
例えば,上の図について書くと,は以下のように値を取ります.
0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | |
1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 |
3 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
4 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 |
5 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 |
6 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 |
7 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 |
8 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
列番号は都市,行番号は周る順番です.
この表を見てわかるように,この問題では,1周する間に各都市は一度だけ訪れ(表において各列に1は1つだけ),各時点で訪れる都市は1つだけ(表において各行において1は1つだけ)という制約があります.
すなわち,制約条件は,全ての都市について,全ての時刻において
が成り立つことです.
目的関数を最小化したときにこれらの制約が自然に満たされるように,目的関数を
とします.は制約が満たされるくらい十分大きな係数です.この目的関数を最小化したとき,第二項,第三項は0に近づくので,制約が満たされることになります.
さて,この目的関数は0と1の値をとる変数のみで記述されており,イジング模型のハミルトニアンと似ていることがわかります.
ただし,スピンはー1,1をとるため,0,1に変換しなければなりません.-1,1を0,1に変換するには
とすれば,のときに
,
のときに
となります.
この変換は最適化を行った後に変換してもよいですし,を目的関数に代入して
の関数で表しても,
と,変数によらない定数
がくっつく形になるだけあり,最小化する過程ではこの定数は意味をなさないため,目的関数は変数をー1,1だと思っても本質的には変わらず,結局,巡回セールスマン問題の目的関数と同じ形の
をハミルトニアンとして用いればよいことになります.
ただし,量子アニーリングでは目的関数ではなく,ハミルトニアン(演算子)に対する固有値が最小になるように計算するので,は演算子であるパウリ行列に書き直しました.
さて,これで巡回セールスマン問題の目的関数をイジング模型のハミルトニアンに書き直すことができました(とはいえ,大きな変更は加えていませんが).
このように,あらゆる組合せ問題はイジング模型のハミルトニアンとして書き直すことができ,書き直すことで量子アニーリングを実行することができます.
外部磁場を加える
量子アニーリングでは,このイジング模型のハミルトニアンに,外部磁場によるハミルトニアンを加えます.
もとの解きたい問題のハミルトニアンをとすると(先ほどの巡回セールスマン問題の例で言えば,
は
),量子アニーリングのハミルトニアンは
のようにします.は横磁場のハミルトニアンです.
は磁場の強さ,
はx方向横方向のパウリ行列で,
です.z成分のパウリ行列と同様,はスピン
のみに作用します.一応定義を書いておくと,
です.
ここで,一つのスピンに対するx成分のパウリ行列の意味について考えてみます.上向きスピン,下向きスピンを作用させてみると,
となり,スピンを上下反転させることが分かります.
また,上下スピンを重ね合わせたスピン状態
をに作用させると,
となり,は
の固有状態になっていることが分かります.
すなわち,横磁場を受けたスピンは,上下スピンが等確率で重なっている状態であるということが分かります.

量子アニーリングでは,この横磁場を弱めていくことで,本来求めたいハミルトニアンの形に徐々に近づけていきます.
量子アニーリングの Step2.外部磁場を下げながらシュレディンガー方程式を解く
量子アニーリングのアルゴリズムの手順は次の通りです.
0.ステップで
とし,初期解を
とする.
1.とし,ハミルトニアン
に対してシュレディンガー方程式を用いて固有値,固有状態を求める(正確には,状態を時間発展させる).このとき,
は小さくなる.
2.終了条件(が十分小さくなる,
となる)を満たせば終了.そうでないなら,1に戻る.
これを詳しく説明していきます.
初期解
量子アニーリング開始時のハミルトニアンは,
,
は十分大きな値にしておきます.つまり,初期ステップ(時刻)でのハミルトニアンは
です.例えば,十分大きな定数
を用いて
としておけば,,
となります.この場合,ステップ
が終了条件となります.
初期解は,
のみに対する解であり,これはすべてのスピン状態を等確率で重ね合わせたものになるので,
です(2個下の図の左の青い部分).
これは,の±1の固有状態
を書くサイト
について積を取ったものであり,これを展開すると,すべてのスピン状態を等確率
で足し合わせたものになっています.
磁場を弱める
を小さくしていき,その都度,シュレディンガー方程式を解いて固有状態(スピン状態,解)を求めいていきます(ただし,後で述べるように,古典コンピュータではシュレディンガー方程式を毎回解く(状態を時間発展させる)のは現実的ではないので,他の方法を用います).
そして,,つまり横磁場が0となったときの各スピンの状態(z方向に上向きか下向きか)が,(不思議と?)最適解となっているというわけです.

ハミルトニアンと各スピンのz成分の確率に注目して言えば,ハミルトニアンは始めはであり,初期解はすべてのスピン状態は等確率をとっているので,平坦なハミルトニアンに等確率で状態が存在していることになり,最終状態では
は本来の形になり,特定のスピン状態が高い確率で存在しているというイメージです.

実際に古典コンピュータでシミュレーションするには
さて,これまでで,量子アニーリングのアルゴリズムについて説明しました.
しかし,以上で説明したのは,量子アニーリング用マシンで実現できるものです.なぜなら,シュレディンガー方程式は現実の自然な法則なので,放っておけば固有値,固有状態が自然と実現されるからです.
古典コンピュータで量子アニーリングをシミュレーションするには,シュレディンガー方程式を毎回計算する(状態を時間発展させる)必要がありますが,シュレディンガー方程式をそのまま毎回計算するのは膨大な計算量がとなるため,現実的ではありません.
したがって,古典コンピュータで量子アニーリングをシミュレーションするには,量子モンテカルロ法などの方法を用いて計算します.
つまり,量子アニーリングをお手元のコンピュータで動かしたいなら,この記事に書いた量子アニーリングの流れだけ知っていても,シミュレーションできず,量子モンテカルロ法について学ぶ必要があるということです.
量子モンテカルロ法については,次の記事で詳しく説明しています.
コメント