はじめに
現在, ゲート型量子コンピューターはいくつかの方式で実際に開発されてはいる(IBM, Google, Rigetti, IonQ, Honeywell etc.)ものの, ノイズによる影響が大きく正しい計算が行えない. そこで, そのようなノイズの影響を取り除くことができる誤り耐性型量子コンピューターが開発されるまではノイズのある量子コンピューター(Noisy Intermediate Scaled Quantum, NISQ)と付き合っていく必要がある. NISQ用のアルゴリズムとして近年変分法による, 古典コンピューターと量子コンピューターによるハイブリッドアルゴリズムが注目されており, その1つとしてVQE (Variational Quantum Eigensolver; 変分量子固有値ソルバー)が知られている. 本記事ではVQEアルゴリズムの数学とその応用例を紹介する.
なお, 現在量子コンピューターとそのアルゴリズム開発・(産業)応用に対して世界的に関心が高まっており, 米欧中を中心に盛んに研究されている. 日本でも内閣府の政策:
ムーンショット目標6
において, "2050年までに、経済・産業・安全保障を飛躍的に発展させる誤り耐性型汎用量子コンピュータを実現" が掲げられており, 来たる量子情報社会に向けた取り組みが行われている.
準備
を有限次元複素ヒルベルト空間, の内積をで表し, この内積が誘導するノルムをで表す. ただし内積は, 量子力学あるいは量子情報科学の慣例に倣い, 右線型・左反線型とする. 本記事はVQEの数学についての解説を目的にしているためブラケット記法は採用せず, 通常の数学において取り扱われる記法を用いる.
は有限次元であるため, に対しである. 以後, を, 作用素を対応する行列, その他適宜読み換えてよい.
次の問題を考えよう.
上のエルミート作用素に対し, 最小固有値とその固有ベクトルを求めよ.
この問題は, 筆者は詳しくないが, 物理や量子化学において系の最小エネルギーとその基底状態を求める重要な問題であるという. これを数学的に表現したものが上記である.
このような固有値に関する問題に対しては, 誤り耐性型量子コンピューターでは量子位相推定アルゴリズム(_quantum phase estimation_)が知られているが, 必要となるゲート(ユニタリー作用素)の数が非常に多く, 精度に応じて必要な量子ビット数も増加するため, NISQデバイス用としては難しいといわれている.
量子計算のいくつかの基礎
ここに量子計算のための基礎的な事項を整理しておく. より詳しくはNielsen-Chuang [2], 石坂[4]などを参照されたい.
状態
でなるものを状態とよぶ. また, 状態があるに対しであるとき, これらを状態として区別しない.
をの双対空間, すなわち上の有界線型汎関数全体とする. に対し, リースの表現定理によるとの対応をで表す. すなわち, である.
エルミート作用素のスペクトル分解・固有値分解
を上のエルミート作用素とする. このとき, は実数の固有値を持ち, 次が成り立つ:
- の固有値, と対応する固有空間への射影をとすると, は
に分解される. この分解をスペクトル分解という. - の固有値と対応する固有ベクトルによって
と分解される. ただし, とみなす. この分解を固有値分解という.
を上のエルミート作用素, スペクトル分解とする. また, を状態とする. このとき,
とおくと, は右連続な単調非減少関数である. したがって, のルベーグ・スチルチェス測度はボレル可測空間上の確率測度である. 特にでに集中するディラック測度とすると,
である.
, , とおく.
と定めると, は上で定義された可測関数である. または確率空間上の確率変数である.
ボルンの確率規則
状態のエルミート作用素を測定するとは, 上で定めた確率変数をサンプルすることである. エルミート作用素は物理量と呼ばれる.
定義より直ちに
が分かる.
次は線型代数でもよく知られている命題である.
変分原理
を上のエルミート作用素とし, の最小固有値をとする. このとき,
が成り立つ.
固有値分解とする. このとき, 任意のに対し,
であることより主張を得る.
VQE
を有限次元複素ヒルベルト空間とする.
Ansatz
上のパラメータ付けされたユニタリー作用素の族をansatzという.
パラメータ付きのユニタリー作用素
, とすると,行列
はパラメータ付きユニタリー作用素である. これらの組合せはansatzに用いられる.
さて, VQEについて説明する.
初期状態, 上のansatzとし, とおく. ただし, とする. またを上のエルミート作用素とする.
とおく. このとき, がに関して強連続ならば, はに関して連続である. よって, がコンパクトであればはにおいて最小値を持つ. の計算をNISQマシン(量子コンピューター)が担当し, 最小化・最適化アルゴリズムを古典コンピューターが担当するハイブリッドアルゴリズムがVQEである. 以上の手順によって, の近似値を得られることが期待できる. まとめると次の通りである.
VQE:
- Ansatzを決定する.
- とする.
- とする.
- 初期状態にansatz を作用させる.
- 状態についてを測定し, を得る.
- が十分小さいならば終了, そうでなければ 7. へ
- 古典コンピューター上において最適化アルゴリズムを用いを求める.
- に1を加え, 3.に戻る.
7.で用いられる最適化アルゴリズムとしてBFGS法をはじめ, 深層学習でよく用いられるAdam等が挙げられる. 変分量子計算に特化した最適化アルゴリズムとしてはNakanishi-Fujii-Todo法(NFT法)が近年提案されている. NFT法はIBM QiskitにもAPIが実装されている.
変分原理により, ansatzが適切であればは最小固有値を, はその固有状態を近似することが期待される. 実際, 次が知られている.
の固有値をとする. の重複度を1とし, 対応する固有状態(単位固有ベクトル)をとする. さらに, とする. このとき次の不等式が成り立つ:
命題2と同じ仮定とする. このとき, を満たす任意のに対し,
である. したがって, に対し, 十分大きなでを満たし, かつならばである.
応用例:連立1次方程式への適用
, 正則行列, とする. 未知変数とする, 次の線型方程式を考える:
とおく. このとき次が成り立つ:
は半正定値エルミート行列でその最小固有値はである. さらに, 方程式(1)の解はの固有値に対応する固有ベクトルである.
がエルミートであることは明らか. 残りの主張は, がによって張られる部分空間の直交補空間への直交射影であることから分かる.
以上より, に対してVQEを適用し,
を求めることで方程式(1)の解の近似値が得られる.
その他への適用
他にも, 特に分子動力学や第一原理計算において基底状態を求めるのに有効とのことである. 筆者は当該分野に全く詳しくないので, 興味のある読者は当該分野の専門家や専門書籍を参照されたい.
その他, 金融・ファイナンス分野への適用については本記事に追記, または別記事にて紹介するかもしれない.
おわりに
量子コンピューターといえば大人気アニメ「神様になった日」が思い起こされる. エンディングテーマ「Goodbye Seven Seas」は神曲なので是非聞いてみて欲しい.
参考
- Nakanishi, K., Fujii, K. and Todo, S.
Sequential minimal optimization for quantum-classical hybrid algorithms
, 2020, Phys.Rev.Research, vol.2, 043158.
- Nielsen, M.A. and Chuang, I.L. Quantum Computation and Quantum Information, 2010, Cambridge University Press.
- Xu, X., Sun, J., Endo, S., Li, Y., Benjamin, S.C. and Yuan, X.
Variational algorithms for linear algebra
, 2021, Sci.Bull., vol.66, 21, pp.2181-2188.
- 石坂 智・小川 朋宏・河内 亮周・木村 元・林 正人, 量子情報科学入門, 2012, 共立出版.
- 内閣府, ムーンショット目標6,
https://www8.cao.go.jp/cstp/moonshot/sub6.html