[論文レビュー] Review of a Quantum Algorithm for Betti Numbers

Sam Gunn, Niels Kornerup|arXiv (Cornell University)|Jun 18, 2019

Topological and Geometric Data Analysis参考文献 10被引用数 9

ひとこと要約

本稿では、組合せラプラシアンのブロック対角行列表現に対する位相推定を用いて、単体複体のベッチ数を計算する簡素化された量子アルゴリズムを提示する。この手法は先行研究より効率が向上しているが、一般には多項式時間でない。特定の条件下（大きなベッチ数と有界なスペクトルギャップ）でのみ多項式実行時間が達成可能である。

ABSTRACT

We looked into the algorithm for calculating Betti numbers presented by Lloyd, Garnerone, and Zanardi (LGZ). We present a new algorithm in the same spirit as LGZ with the intent of clarifying quantum algorithms for computing Betti numbers. Our algorithm is simpler and slightly more efficient than that presented by LGZ. We present a thorough analysis of our algorithm, pointing out reasons that both our algorithm and that presented by LGZ do not run in polynomial time for most inputs. However, the algorithms do run in polynomial time for calculating an approximation of the Betti number to polynomial multiplicative error, when applied to some class of graphs for which the Betti number is exponentially large.

研究の動機と目的

• Lloyd ら（2014）が提唱したベッチ数計算の量子アルゴリズムを簡素化・改善すること。
• ヴェイトリス＝リプスフィルトレーションから得られる単体複体の文脈において、位相的不変量の量子計算を明確にすること。
• アルゴリズムの計算複雑度を解析し、多項式実行時間に到達可能な条件を同定すること。
• 実用的応用により関連性の高い恒常ホモロジーへのこの手法の拡張可能性を評価すること。

提案手法

• マーク済み要素の数が未知または推定値である場合に適応したグローバー法の変種を用いて、グラフ内のすべてのk-単体に均一な重ね合わせ状態を準備する。
• k-単体をハミング重みkを持つn-キュービット状態として符号化することで、キュービット数を削減し、近似カウントと振幅増強を用いた効率的な状態準備を可能にする。
• 各k-単体次元におけるラプラシアンΔkの二乗が得られる、ブロック対角ヘルミート行列Bに対する位相推定を用いる。
• 状態準備ユニタリP_kとその随伴を用いて、k-単体のヒルベルト空間に制限された修正されたグローバー拡散演算子を適用する。
• 位相推定の結果をオракルとして用い、量子カウントを用いてΔkの核の次元（ゼロ固有値の数）を推定する。
• 主要な要素として、全体のゲート複雑度と実行時間に影響を与える最小非ゼロ固有値λminの境界を求める。

実験結果

リサーチクエスチョン

• RQ1 LZG提案のオリジナル版よりも単純かつ効率的なベッチ数計算のための量子アルゴリズムを設計できるか？
• RQ2 提案された量子アルゴリズムが、どのような条件下でベッチ数計算において多項式実行時間に到達できるか？
• RQ3 実用的応用により関連性の高い恒常ホモロジーに向けたアルゴリズムの拡張は可能か？
• RQ4 ラプラシアンのスペクトルギャップλminがアルゴリズムの実行時間と実現可能性に果たす役割は何か？
• RQ5 なぜオリジナルのLZGアルゴリズムと改善版の両方とも、大多数の入力に対して多項式時間の複雑度を達成できないのか？

主な発見

• 提案されたアルゴリズムはオリジナルのLZGアルゴリズムよりも単純かつわずかに効率が良く、ゲート複雑度はÕ(√(βk|Sk|)(nk√(C(n,k)/|Sk|) + n²k/λmin))である。
• βk|Sk|の平方根に依存し、最小非ゼロ固有値λminの逆数に依存するため、大多数の入力に対して多項式実行時間にはならない。
• βkが指数的に大きくかつλminがゼロから離れている場合にのみ多項式実行時間が達成可能であり、これは実際にはまれな条件である。
• オリジナルのLZGアルゴリズムの複雑度の主張は、U_kにおける射影演算子の使用に関する誤った仮定に基づいており、著者らはこれらが不要であると主張する。
• βkの正確な計算にはδ ≈ 1/βkが必要であり、これにより提案されたアルゴリズムの複雑度に近いものになる。これはβkが小さくかつkが定数である場合を除き、依然として非多項式のままである。
• 恒常ベッチ数の計算へのこの手法の拡張は直感的ではない。恒常ホモロジーには類似のラプラシアン構造が存在せず、標準的な古典的手法も量子設定に容易に拡張できない。