End-to-End Quantum Algorithms for the Jones Polynomial

技術要約：ジョーンズ多項式のためのエンドツーエンド量子アルゴリズム

問題定義
本論文は、結び目および絡み目の位相不変量であるジョーンズ多項式を、5の1の冪根（$t = e^{i 2\pi/5}$）において近似するという計算上の課題に取り組んでいる。この問題は、任意の結び目が「引き抜き（closure）」によって形成される「ブレイド（編み目）」の観点から定式化されている。具体的には、以下の2つの定式化が検討されている：

1. マルコフ・クロージャ（DQC1完全）： ブレイドを表すユニタリ行列の重み付きトレースを推定すること。これは「1個のクリーンな量子ビット」計算モデルに対応する。
2. プラット・クロージャ（BQP完全）： 特定の量子振幅（ユニタリ行列の要素）を推定すること。

DQC1はBQPの中に厳密に包含されていると広く信じられているため（つまりDQC1は「より量子性が低い」）、著者らは、DQC1版（トレースの推定）の古典的リソース要件は、BQP版（振幅の推定）と同等かそれ以上に高いことを指摘している。したがって、著者らは両方のパイプラインをサポートしつつも、より大きな量子優位性を示す可能性が高いとして、DQC1の定式化（マルコフ・クロージャ）に焦点を当てている。

手法
著者らは、ノイズのあるゲート型量子コンピュータ（具体的にはQuantinuum社のH2-2）向けに設計された、エンドツーエンドのアルゴリズム・パイプラインを提示している。コアとなる手法は以下の通りである：

• ユニタリ表現： アルゴリズムは、ブレイド生成元（$\sigma_i^{\pm 1}$）のフィボナッチ・ユニタリ表現を利用する。この表現は、$n$個の量子ビットのヒルベルト空間内の、「連続するゼロを持たないフィボナッチ文字列」によってスパンされる部分空間上で作用する。この部分空間の次元は、$n$番目のフィボナッチ数（$\phi^n$）のようにスケールする（ここで$\phi$は黄金比）。
• アルゴリズム・プリミティブ (cfev)： 標準的なアダマール・テスト（制御ユニタリ操作を必要とし、高いゲートオーバーヘッドを伴う）の代わりに、著者らは**制御フリー・エコー検証（control-free echo-verification: cfev）**プロトコルを採用している。このプロトコルは以下の特徴を持つ：
  • キャット状態 $|\psi_{cat}\rangle = \frac{1}{\sqrt{2}}(|0\rangle^{\otimes n} + |s\rangle)$ を使用する。ここで $|0\rangle^{\otimes n}$ はブレイド・ユニタリ $U_B$ の既知の固有状態である。
  • 制御-$U$ゲートを不要にし、2量子ビットゲートの数を大幅に削減する。
  • 最後に逆回路（$V_{cat}^\dagger$）を用いて、システムが初期状態に戻ったことを検証し、エラー検出を可能にする。
• エラー緩和：
  • 非フィボナッチ・エラー検出： ポストプロセッシング関数（$g_2$）により、フィボナッチ部分空間から外れた測定結果を破棄し、特定のハードウェアエラーをフィルタリングする。
  • 共役トリック（Conjugate Trick）： コヒーレントな位相エラー（特にトラップイオン量子ビットのメモリ・エラー）を緩和するために、ブレイド $B$ とその共役 $B^*$（生成元が反転したもの）の両方でアルゴリズムを実行する。これらの結果を組み合わせることで、追加のショットを必要とせずにグローバル位相エラーを除去する。
• 古典的ベースライン： 著者らは、同一の問題に対する最先端の古典的アルゴリズムを開発し、ベンチマークを行っている。これには以下が含まれる：
  • tn-proj: 直接的なテンソルネットワーク収縮。
  • sv-shor: フィボナッチ部分空間の疎性を利用した状態ベクトル・シミュレーション。
  • mpo-proj: ボンド次元圧縮を伴う行列積演算子（MPO）収縮。
  • bracket: ブラケット多項式（Kauffman）に基づく古典的手法。

主な貢献

1. エンドツーエンド・パイプライン： ランダムなブレイドを生成し、それらを特定のハードウェア（Quantinuum H2）向けに最適化された量子回路へとコンパイルし、実行し、古典的な前処理および後処理を行う、完全に構成可能なパイプライン。
2. 効率的に検証可能なベンチマーク： ジョーンズ多項式の位相不変性を活用し、単純なブレイド $B$ をランダムなブレイド $A$ で共役させることで生成された複雑なブレイド $B'$ を構築する。これにより、$B'$ は $B$ と位相的に等価となる。$VM(B)$は古典的に無視できる時間で計算可能であるため、$VM(B')$ は量子プロセッサのノイズとエラーのスケーリングを測定するための検証可能なグラウンドトゥルース（正解）として機能する。
3. 最適化されたコンパイル： フィボナッチ・ユニタリ生成元の最適化された回路分解を提供し、同様の問題に対する従来の推定と比較して、2量子ビットゲートの数を約15倍削減した。
4. リソース推定フレームワーク： ノイズモデルと古典的ハードウェアの限界（A100 GPUやFrontierスーパーコンピュータなど）を考慮した、量子と古典のアプローチ間の「解決までの時間（time-to-solution）」およびエネルギー消費のクロスオーバーポイントを推定するための手法。

結果

• 実験的実証： アルゴリズムは、15ストランド、104クロッシングを持つブレイド（16量子ビット、340個の2量子ビットゲート）に対して、QuantinuumのH2-2量子コンピュータ上で実行された。
  • 緩和なしの場合、相対誤差は75%、標準偏差は9%であった。
  • 非フィボナッチ・エラー検出と共役トリックを適用することで、相対誤差は43%に減少した（ただし標準偏差は14%に増加した）。
• スケーリング分析： 脱分極ノイズモデル（$\epsilon_{2q} = 10^{-4}$ および $5 \cdot 10^{-4}$）を用いた24,000個のブレイド（10–15ストランド、50–600クロッシング）のデータセットによるシミュレーションにより、以下が明らかになった：
  • 相対誤差は主にクロッシング数（ゲート数）とともにスケールする。
  • $\epsilon_{2q} = 5 \cdot 10^{-4}$ の場合、量子エラーが非常に大きくなり、すべてのテストされたサイズにおいて古典的なMPO近似（mpo-proj）の方が高速になる。
  • $\epsilon_{2q} = 10^{-4}$ の場合、量子アルゴリズム（cfev）は、約2,800クロッシングのブレイドサイズにおいて、最高の古典的手法（mpo-proj）をFrontierスーパーコンピュータ上で上回ると予測される。
• エネルギー効率： 分析によれば、量子プロセッサの消費電力を75kWと仮定した場合、$\epsilon_{2q} = 10^{-4}$ において、量子アプローチは2,400クロッシング付近で古典的なmpo-proj法よりもエネルギー効率が高くなる。

意義と主張
本論文は、抽象的な計算複雑性クラスや検証されていないランダム回路のベンチマークに頼るのではなく、計算トポロジーにおける意味のある問題に対して、量子優位性を特定するための実用的かつ定量的なアプローチを提供することを主張している。

• 優位性に関する控えめな主張： 著者らは、自らの結果が量子優位性が可能となる場所の**保守的な推定値（下限）**であることを明示している。彼らは、ベンチマークに使用したブレイドが特定の構造を持っており、これらが「平均的なケース」や「最も困難な」インスタンスを代表していない可能性があることを認めている。しかし、もし一般的なブレイドがこれらのベンチマークセットと同等に困難であるならば、特定されたクロスオーバーポイントは、一般的なインスタンスにおいて優位性が観察される可能性のある有効な上限値であると論じている。
• 議論の転換： 本研究は、量子優位性の議論を「問題のクラス」から、エラー、時間、エネルギーの厳密かつ多次元的な最適化を用いた「具体的な問題のインスタンス」へとシフトさせることを目的としている。
• 有用性： このパイプラインは、量子プロセッサのノイズを特性化し、現在の古典的スーパーコンピュータにとって手に負えない非自明なジョーンズ多項式の問題を解くために必要な最小限のハードウェア仕様（量子ビット数、ゲート忠実度、深さ）を決定するためのツールとして機能する。

著者らは、高精度なアプリケーションには誤り耐性のある領域とさらなる誤り訂正が必要であるが、現在のNISQ（Noisy Intermediate-Scale Quantum）アプローチは、結び目理論における実用的な量子有用性を特定するための実行可能な道筋を示していると結論付けている。