Low-Depth Initial-State Preparation for Ground-State Energy Estimation of Two-Dimensional Strongly Correlated Systems

技術要約：二次元強相関系の基底状態エネルギー推定のための低深度初期状態準備

問題提起
量子位相推定（QPE）などの基底状態エネルギー推定のための量子アルゴリズムには、真の基底状態に対して非無視できない忠実度を持つ初期状態が必要である。成功確率は、この忠実度（$F_0$）に依存する。もし $F_0$ がシステムサイズに対して指数関数的に減少する場合、アルゴリズムは必要な反復回数（$1/F_0$）の増大により、計算不能となる。変分量子固有値ソルバー（VQE）はハイブリッドなアプローチを提供するが、勾配評価のコストや浅い回路の古典的シミュラビリティといった課題に直面している。さらに、正確な古典近似（例：行列積状態）を量子回路にコンパイルすることは、過大な回路深度をもたらすことがあり、これは論理量子ビットやマジックステートのリソースが限られている初期の耐故障性量子コンピュータ（FTQC）にとって致命的である。中心となる課題は、深い最適化ループや古典的な計算困難性に依存することなく、強相関系（特に二次元ハバードモデル）に対して高忠実度の初期状態を準備するための、浅くリソース効率の高い量子回路を設計することである。

手法
著者らは、シュリーファー・ヴォルフル（SW）変換を介した、ハイゼンベルクモデルと半充填ハバードモデルの間の強結合関係を利用することで、浅い量子回路を構築する戦略を提案している。この手法は主に以下の3つの段階で進行する：

1. ハイゼンベルク回路の構築:

   • アンザッツ: 回路は、反強磁性ハイゼンベルクモデルの $S=0$（シングレット）部分空間を対象とする。初期状態は、格子上の完全マッチング上に定義されたシングレット・ダイマー（$| \Phi \rangle$）の積である。
   • ダイナミクス: 状態は、指数関数的なeSWAPゲートの層を用いて進化する。これらのゲートはSU(2)対称性を保持するため、状態はシングレット部分空間内に留まる。
   • 構造: 格子の結合は、4つの互いに素な部分集合（$B_{x0}, B_{x1}, B_{y0}, B_{y1}$）に分割される。回路は、複数のブロック（$N_b$）にわたって、特定の順序（$B_{y1} \to B_{y0} \to B_{x1} \to B_{x0}$）でeSWAPゲートを適用する。この構造は、ハイゼンベルク基底状態の共鳴価結合（RVB）物理を模倣している。
   • フェルミオン写像: スピン状態は、補助量子ビットを導入する局所的なユニタリ埋め込み（$V$）を通じて、単一占有フェルミオン部分空間へと写像される。

2. SWゲートによるハバードへの拡張:

   • 有限の $U/t$ におけるハバードモデルに存在する電荷揺らぎ（ダブロンとホロン）を捉えるために、ハイゼンベルク回路を電荷揺らぎゲートで「ドレス」する。
   • これらのゲートは、一次のSW生成子（$\Omega$）から直接導出され、二サイトユニタリ（$W_{\langle i,j \rangle}$）への鈴木・トロッター分解によって近似されている。
   • 最終的な変分アンザッツは、ハイゼンベルク発展（$U_B$）と電荷揺らぎブロック（$W_B$）を組み合わせ、回転角 $\theta$ と $\phi$ によってパラメータ化される。

3. 忠実度の推定とベンチマーク:

   • パラメータ転送: 小規模な $4 \times 4$ 格子で最適化されたパラメータは、再最適化なしに、より大きな格子（$6 \times 6$ から $10 \times 10$）へと転送される。
   • シミュレーション: 回路状態は、信念伝播（BP）ガイド付き単純アップデートを用いた投影絡み合いペア状態（PEPS）を用いてシミュレートされる。
   • 忠実度の下限: 大規模な系における厳密な基底状態は未知であるため、著者らは基底状態への忠実度の下限を推定する。これは以下によって達成される：
     • 回路状態と変分モンテカルロ（VMC）参照状態との間の忠実度（$F_{ref}$）の計算。
     • 参照エネルギー、厳密な基底状態エネルギー（ストーク・シリーズ展開による）、およびシングレット励起ギャップ（$\Delta_0$）を用いた、VMC状態の基底状態への忠実度（$F_{ref,0}$）の境界設定。
     • フビニ・スティューディ角を用いた三角不等式による、$F_0$ の厳密な下限の導出。

主要な貢献と結果

• 効果的なパラメータ転送: ハイゼンベルクモデルにおいて、$4 \times 4$ 格子で最適化されたパラメータは、再最適化なしに最大 $10 \times 10$ の格子に対して堅牢な忠実度下限を与える。具体的には、$N_b=3$ の場合、推定された基底状態への忠実度下限は $10 \times 10$ 格子においても正の値（$\approx 0.028$）を維持しており、この初期状態が振幅増幅に適していることを示唆している。
• ハバード忠実度の向上: $4 \times 4$ 格子における厳密な計算により、SW由来の電荷揺らぎゲートを組み込むことで、ハバード基底状態への忠実度が大幅に改善されることが示された。$N_b=4$ の場合、忠実度は $\approx 0.44$（ハイゼンベルクのみ）から $\approx 0.89$（最適化された電荷角を含む場合）へと増加する。
• リソース効率: 著者らは、$10 \times 10$ 格子の状態準備コストを推定している。回転あたりのエラーを $10^{-5}$ と仮定したClifford+$T$合成を想定すると、総コストは $10^5$ 個の $T$ ゲートのオーダーとなる。これは、状態準備コストを初期のFTQCにおける「メガクオップ（megaquop）」領域内に収めており、その後の基底状態エネルギー推定アルゴリズムに伴うコストよりも大幅に低い。
• スケーラビリティ: 回路の深さは格子結合数（$|B|$）に対して線形にスケールし、パラメータ共有戦略により、小規模系の最適化に基づいた大規模系の回路構築が可能となる。

意義と主張
本論文は、強相関系における基底状態エネルギー推定のための低コストな初期状態準備への実行可能な経路を提示していると主張している。ハイゼンベルクモデルとハバードモデルの間の物理的な接続を利用することで、著者らは、浅い回路が厳密な古典的対角化の手が届かない格子に対しても、非無視できない基底状態への忠実度を達成できることを示している。これらの結果は、回路の深さと $T$ ゲート数を最小化することが極めて重要となる初期FTQCレジームにおいて、これらの回路が魅力的な候補であることを示唆している。本研究は、パラメータが転送可能であり、かつ回路の深さが低く抑えられている限り、小規模格子のシミュレーションが、古典的に困難な大規模系のための量子回路設計を効果的に導くことができるという仮説を支持している。著者らは、提示された境界は厳密なものであるが、実際の忠実度はこれよりも高い可能性があることを明記しており、今後の研究では、より深いハミルトニアン変分アンザッツや、より正確なシミュレーション手法（例：フルアップデートPEPS）を用いて、これらの推定をさらに精緻化できる可能性があるとしている。