Quantum Folding in Chebyshev-Based Quantum Subspace Diagonalization

技術要約：チェビシェフに基づく量子部分空間対角化における量子フォールディング（Quantum Folding）

1. 問題設定

量子部分空間対角化（QSD）、別名量子クリロフ法は、量子コンピュータ上で生成された低次元の部分空間にハミルトニアン $H$ を投影し、その結果得られる一般化固有値問題を古典的に解くことで、基底状態のエネルギーを推定する手法である。チェビシェフ多項式を用いるバリアントは、トロッター誤差のない厳密な実現が可能であるが、その理論的保証（具体的には収束率、ノイズフロア、および状態準備のコスト）は、参照状態 $|\psi_0\rangle$ と真の基底状態 $|g\rangle$ との初期オーバーラップ $|\gamma_0|^2 = |\langle g | \psi_0 \rangle|^2$ に決定的に依存する。

任意の固定された積状態（例：反強磁性体におけるネール状態）に対して、このオーバーラップは非自明なハミルトニアンに対してシステムサイズとともに指数関数的に減少する。この「直交性の破綻（orthogonality catastrophe）」により、従来のQSDは大規模な系に対して実用的ではなくなる。なぜなら、ノイズフロアの上に基底状態の信号を識別するために必要な測定予算と回路深さが指数関数的にスケールするためである。本論文は、参照状態の選択が、QSDがスケーラブルであるか否かを決定する主要なボトルネックであることを特定している。

2. 手法：量子フォールディング（Quantum Folding）

著者らは、局所的なハミルトニアンに対して指数関数的なオーバーラップの減少を回避する参照状態を構築するための戦略として、量子フォールディングを提案している。

• フォールディングされた状態（The Folded State）: システムを部分系 $A$ と $B$ に分割し、局所ハミルトニアンを $H_\Lambda = H_A \otimes \mathbb{1}_B + \mathbb{1}_A \otimes H_B + H_\partial$ （ここで $H_\partial$ は界面項を表す）とする場合、フォールディングされた参照状態は $|\psi_0\rangle = |g_A\rangle \otimes |g_B\rangle$、すなわち孤立した部分系の厳密な基底状態の積として定義される。
• 物理的正当性: この状態は、界面項を除去したフルハミルトニアンの厳密な基底状態である。したがって、真の基底状態との差異は界面に限定される。1本の鎖を半分に切断した場合、界面は単一の結合（ボンド）となる。ゆえに、オーバーラップはバルクの体積ではなく、この単一の結合の物理に支配される。これにより、ギャップのある系ではサイズに依存せず、ギャップレスな系では緩やかなべき乗則に従ってのみ減少する。
• 再帰的ラダー（Recursive Ladder）: 本手法は再帰的な構成（$L \to 2L \to 4L \dots$）を用いる。サイズ $L$ におけるQSD計算の出力が物理的な状態として準備され、サイズ $2L$ におけるQSD計算のシードとして使用される。

3. 主要な貢献と理論的展開

A. チェビシェフQSDの誤差境界
著者らは、有効ハミルトニアン解析（以前に実時間発展のために開発されたもの）をチェビシェフ・クリロフ基底へと移植している。

• 線形ノイズスケーリング: 著者らは、基底状態エネルギーに関する非漸近的な誤差境界を導出し、それが測定ノイズ（$\eta$）に対して線形であることを示した。これは、従来の劣線形スケーリングを改善するものである。
• フィルタノルムの制御: 重要な要素は、係数ベクトルのユークリッドノルムが $O(1)$ であるチェビシェフ低エネルギーフィルタの構築であり、これによりノイズや閾値処理による誤差が指数関数的に増幅されないことを保証している。
• 再スケーリングへの依存性: 解析により、ブロックエンコーディングによるハミルトニアンの再スケーリングが保証のパラメータであることが明らかになった。基底状態をチェビシェフ区間の端に配置すると、収束指数は $\sqrt{\Delta}$（ギャップ）に比例するが、内部に配置する場合（固定されたブロックエンコーディング正規化において一般的）は $\Delta$ への線形依存性を示す。

B. 状態準備とポストセレクション
本論文は、古典的なリッツベクトル（非直交基スの係数ベクトル）をいかにして物理的な量子状態に戻すかという課題に対処している。

• 実装: この状態は、ウォーク演算子の累乗を用いた、ユニタリの線形結合（LCU）または量子特異値変換（QSVT）を介して準備される。
• 閾値処理としてのフロア: 著者らは、エネルギー固有値問題を安定させるために必要な閾値処理が、ポストセレクション確率に下限（フロア）を課すことも証明している。閾値処理を行わない場合、オーバーラップ行列の近零部分空間におけるノイズが、成功確率が極めて低い偽のリッツベクトルを引き起こす。
• 認定されたフィルタ（Certified Filters）: フィルタの効率が悪いために基底状態生成確率が抑制されるのを防ぐため、著者らは「認定されたフィルタ」（リッツ値から古典的に構築された多項式）または $\ell_1$ ゲージ固定の使用を提案している。これにより、生成確率が $\Theta(|\gamma_0|^2)$ となり、再帰的なフォールディング・ラダーがスケーラブルになることが保証される。

C. スケーリング解析

• ギャップのある系: ギャップのある鎖の場合、フォールド・フィデリティ $F(L) = |\langle g_{2L} | g_L \otimes g_L \rangle|^2$ は非ゼロの定数 $F_\infty$ に収束する。再帰的ラダーの総コストは、システムサイズ $L$ に対して多項式スケールとなる。
• ギャップレスなスピン1/2鎖: ギャップレスなスピン1/2鎖では、オーバーラップは緩やかなべき乗 $L^{-0.107}$ で減少する。その結果得られるコストは準多項式（quasi-polynomial）であり、現実的なシステムサイズにおいては実質的に多項式スケールとなる。

4. 結果

著者らは、反強磁性ハイゼンベルク鎖（$S=1/2$ および $S=1$）およびAKLT鎖の厳密な数値シミュレーションを通じて、本手法を検証している。

• オーバーラップの優位性: フォールディングされた参照状態は $O(1)$（または緩やかなべき乗則による減少）のオーバーラップを維持するのに対し、ネール状態のオーバーラップは指数関数的に減少する。24サイトまでの $S=1/2$ 鎖において、フォールディングされたオーバーラップは $>0.5$ を維持しているが、ネール状態のオーバーラップは $\sim 10^{-3}$ である。
• 収束速度: ノイズのないシミュレーションにおいて、フォールディングされた状態は、与えられた精度に対して必要なクリロフ次元 $D$ をネール状態と比較して4〜5倍減少させる。
• ノイズ耐性: 測定ノイズの下で、フォールディングされた状態は、ネール状態よりも2〜5桁低い誤差を達成する。ネール状態の誤差は、その極めて小さなオーバーラップによって決定される高いフロアで飽和するが、フォールディングされた状態の誤差は、$D$ の増加とともに減少し続ける。
• フォールド・フィデリティ:
  • ギャップレス スピン1/2: $F(L) \approx 1.09 L^{-0.107}$ （開境界条件）。
  • ギャップのあるスピン1（ハイゼンベルク）: $F(L) \to F_\infty \approx 0.22$。
  • AKLT鎖: 行列積構造から導出された閉形式に従い、正確に $F(L) \to 1/4$ となる。

5. 意義と主張

本論文は、従来の積状態による初期化に伴う指数関数的に小さいオーバーラップを回避する、スケーラブルな量子基底状態準備の枠組みを提供すると主張している。

• スケーラビリティ: オーバーラップをバルク依存から界面依存へと置き換えることで、量子フォールディングはQSDのスケーリングをシステムサイズに対して指数関数的から多項式的へと変貌させる。
• 普遍性: このアプローチは1次元の鎖やスピンに限定されず、界面が部分系に対して小さいあらゆる局所ハミルトニアンに適用可能である。
• 実用性: 本手法は標準的なブロックエンコーディングおよび量子化（qubitization）技術に基づいているが、収束、ノイズ、および状態準備コストのトレードオフを管理するために、特定のプロトコル（再帰的フォールディング、認定フィルタリング、および閾値処理）を導入している。
• 理論的洞察: 本研究は、チェビシェフ法におけるブロックエンコーディングの再スケーリングの役割を明確にし、基底状態を内部に配置することが、後処理では回復できない重大な分解能およびフィルタ次数のペナルティを課すことを示している。

著者らは、量子フォールディングが、導出された誤差境界および状態準備技術と組み合さることで、誤り訂正可能な量子コンピュータ上での相関多体系のスケーラブルな量子シミュレーションへの実行可能な道筋を提供すると結論づけている。