A fidelity metric for quantum annealing benchmarked by extreme scaling quantum Monte-Carlo simulations

技術要約：極限スケール量子モンテカルロ・シミュレーションによってベンチマークされた量子アニーリングのフィデリティ指標

問題提起
現在の量子アニーリング（QA）のベンチマークは、主に実用的な観点に偏っており、特定の最適化問題（通常、二次無制約バイナリ最適化（QUBO）として定式化される）を解く際の速度や成功率のみを測定している。著者らは、このアプローチは「本質的に脆弱である」と主張している。なぜなら、結果は特定の問題インスタンスや参照として使用される古典アルゴリズムに大きく依存するためである。さらに、この手法は、アニーリング過程の固有の限界（例：ランダウ・ゼネー遷移を招く指数関数的に消失するエネルギーギャップ）と、外的なハードウェアの限界（例：ノイズ、デコヒーレンス、有限温度）を区別できていない。その結果、最終的な解の質とは独立して、量子アニーリング過程自体のフィデリティを定量化するための堅牢な指標が欠如している。

手法
本論文は、新しいフィデリティ指標として、アニーリング過程における状態方程式 $E(h_x) \equiv \langle \hat{H} \rangle$ の精度（$\epsilon$）を提案している。この指標は、ゲート型量子コンピューティングにおける「ゲートあたりのフィデリティ」に類似している。実験的に状態方程式を得るには、アニーリング・スケジュールの様々な時点において、縦方向のエネルギー（$\langle \hat{H}_z \rangle$）と横方向の磁場エネルギー（$\langle \hat{H}_x \rangle$）の両方を測定する必要がある。

$\epsilon$ のベンチマーク値を確立し、古典シミュレーションの限界をテストするために、著者らはリドバーグ原子系を用いた量子モンテカルロ（QMC）技術による、2つの変分量子モンテカルロ（VMC）のバリアントを用いている：

1. 変分量子モンテカルロ (VMC):

   • アンザッツ: 著者らは、シミュレーテッド・サーマルアニーリングに着想を得た「熱的アンザッツ（thermal ansatz）」を使用している。波動関数は $\psi(n) \propto e^{-E_{\text{eff}}(n)/2}$ としてモデル化され、ここで有効エネルギー $E_{\text{eff}}(n)$ は修正された古典的エネルギー汎関数である。
   • バリアント: 彼らは、「単純な熱的アンザッツ」（$E_{\text{eff}} = \beta_{\text{eff}} F(n)$）と、対称性の破れを捉えるための「スタッガード（交互）熱的アンザッツ」（$E_{\text{eff}} = \beta_{\text{eff}} F(n) + \sum \beta_k N_k(n)$）をテストしている。
   • スケール: このアプローチは計算負荷が低いため、単一のCPUで最大 100,000,000 原子 ($N=10^8$) までのシミュレーションが可能である。
   • 精度: 推定精度 $\epsilon \sim 10^{-2} - 10^{-3}$ を達成している。

2. グリーン関数モンテカルロ (GFMC):

   • 手法: VMC の後処理として、虚時間発展を用いて変分アンザッツを真の基底状態へと投影させる手法である。これは偏りがなく、横磁場イジング模型に対して符号問題が存在しない。
   • スケール: 最大 100,000 原子 ($N=10^5$) までのシミュレーションを実施した。
   • 精度: $\epsilon \sim 10^{-4}$ というより高い精度を達成している。

著者らはこれらの手法を用いて、正方格子および三角格子上のリドバーグ原子、ならびに現実世界のQUBOインスタンス（金融リスク管理に関連する「フォールン・エンジェル問題」）に関する状態方程式を算出している。

主な結果

• スケーラビリティと精度: 熱的アンザッツを用いたVMCアプローチは、相対誤差 $\delta E/E_0 < 10^{-2}$ で最大 $10^8$ 個の原子のシステムをシミュレートすることに成功した。GFMCアプローチは、$N=10^5$ において $\epsilon \sim 10^{-4}$ に達した。
• 誤差の局在化: 状態方程式における相対誤差は、エネルギーギャップが閉じる位置（相転移点）で正確にピークに達する。これはランダウ・ゼネー遷移が発生する領域に対応している。これらの臨界点から離れた場所では、誤差は著しく低い。
• 実験との比較: 著者らは、現在のリドバーグ原子実験プラットフォームは、これらの古典的ベンチマークよりも大幅に低い精度で動作していると推定している。古典シミュレーションは、システムサイズと精度の両面において、現在の実験プラットフォームを数桁上回っている。
• QUBOのベンチマーク: 「フォールン・エンジェル」QUBOインスタンスに対し、単純な熱的アンザッツは $\approx 10^{-3}$ の相対誤差を示した。著者らは、最終的な解の誤差（$\delta F/F$）はショット数を増やすことで任意に減少させることができるが、状態方程式の誤差（$\delta E/E_0$）こそが、ハードウェアまたはアルゴリズムの能力を示すより本質的な尺度を提供すると指摘している。
• 熱的 vs 量子的: 本研究は、精度 $\epsilon \sim 10^{-2} - 10^{-3}$ の範囲内では、量子アニーラーは実質的にその古典的な熱的対応物と区別がつかないことを示唆しており、これは現在の実験ノイズレベルが真の量子優位性を覆い隠している可能性があることを意味している。

意義と主張
本論文は、状態方程式の精度を測定することが、QAのベンチマークにとって重要な前進であると主張している。その理由は以下の通りである：

1. 理想的な断熱進化からの逸脱を引き起こす支配的なメカニズム（固有のギャップの問題と外的なノイズの区別）を直接診断できるため。
2. 特定のQUBOインスタンスの難易度に依存しない、ハードウェアに依存しない指標を提供するため。
3. 将来のハードウェアに対して厳しい制約を課すため。すなわち、量子優位性を実証するためには、実験プラットフォームがこれらの古典的ベンチマーク（具体的には大規模システムにおける $\epsilon < 10^{-4}$）を超える状態方程式の精度を達成しなければならない。

著者らは、「一見単純に見える」熱的変分アンザッツがQAをシミュレートするための強力なツールであることを結論づけている。これは、量子優位性を理解するための道筋が、量子ゆらぎが最も洗練された古典的熱シミュレーションに対して明確な優位性を提供する領域を特定することにあることを示唆している。彼らは、量子アニーリングが現在、古典的なヒューリスティックよりも優れていると主張しているのではなく、将来的にその仮説を厳密にテストするためのツールを提供しているのである。