An efficient algorithm for approximate shadow Hamiltonian simulation

問題の記述
多体系量子系の実時間ダイナミクスのシミュレーションは、物性物理学から量子化学に至るまで、幅広い分野における中心的な課題である。テンソルネットワークやパウリ伝播（Pauli propagation）といった古典的手法は、絡み合い（エンタングルメント）の急速な増大によるボトルネックに直面している。量子アルゴリズムも解決策を提示しているが、既存のアプローチの多く（例：積公式、変分法）は、$N$個の量子ビットの物理レジスタに対して$N$個の量子ビットによる直接的な符号化に依存している。

文献[41]で導入されたシャドウ・ハミルトニアン・シミュレーションは、代替案を提供する。これは、ハミルトニアンとの交換子によって生成される代数内の演算子の期待値の正規化されたベクトルとして、量子状態を符号化するものである。非相互作用系においては、この代数が多項式的に成長するため効率的であるが、相互作用系においてはこの手法は実用的ではなくなる。このような場合、不変性を正確に満たすために必要な代数はシステムサイズ（$N$）に対して指数関数的に増大し、「シャドウ状態」を符号化するために最大で$2^N$個の量子ビットを必要とする。この指数関数的なスケーリングは、物理的に重要なモデルに対する量子優位性の可能性を打ち消してしまう。

手法
著者らは、近似的なシャドウ・ハミルトニアン・シミュレーションのための効率的なアルゴリズムを提案している。核心となるアイデアは、全代数を符号化するのではなく、特定のターゲットとなる観測量に最も関連性の高い演算子の要素のみを系統的に特定し、保持することである。これは、「演算子プルーニング（operator pruning）」を通じて、ルート観測量 $\bar{O}$ から出発して、切り詰められた正規直交演算子基底 $S_{\bar{O}}$ を構築することによって達成される。

論文では、この切り詰められた基底を構築するための3つの異なるプルーニング・スキームを紹介している：

1. 既定の基底（グラフベース）: この手法は、頂点が既定の正規直交基底（例：パウリ文字列）を表す重み付き層状グラフを構築する。エッジは交換子超演算子 $A_H(\cdot) = [H, \cdot]$ の作用を表す。頂点は、ルート観測量に対する関連性を定量化する累積重み閾値 $\epsilon$ に基づいて保持される。このアプローチは、摂動的なハミルトニアン（$H = H_0 + \xi H_1$）において特に効果的であり、その際、閾値はディゾン展開の打ち切り次数に対応する。
2. クライロフ基底（演算子ランチョス）: このアプローチは、演算子ランチョス・アルゴリズムを用いてクライロフ部分空間を反復的に構築する。これは、前の正規直交演算子に交換子超演算子を適用し、その結果を直交化することで新しい演算子を生成する。乗法的打ち切り基準が用いられ、直交成分のノルムの積が閾値 $\epsilon$ を下回った時点で反復が終了する。この手法は適応的に最も関連性の高い演算子を見つけ出すが、相互作用系においてクライロフ・ベクトルを疎な基底（パウリ文字列など）で表現しようとすると、ベクトルが次第に高密度になり、指数関数的な古典メモリコストに苦しむことになる。
3. ハイブリッド・アプローチ: クライロフ法のメモリ・スケーリングの問題を克服するために、著者らはハイブリッド・スキームを提案している。まず、既定の基ース法を用いて、物理的なサイズに対して多項式スケールする、プルーニングされた演算子集合 $S_{\bar{O}}$ を特定する。その後、この制限された部分空間内でクライロフ基底を構築する。これにより、クライロフ・ベクトルが（$S_{\bar{O}}$ の要素のみを含むため）疎な状態に保たれる一方で、ランチョス手順の適応的な効率性の恩恵を受けることができる。

主要な貢献

• アルゴリズムの枠組み: 本論文は、制御された近似フレームワークを提示し、指数関数的な量子ビット・オーバーヘッドのボトルネックを克服することで、相互作用する多体系へのシャドウ・ハミルトニアン・シミュレーションの拡張を実現した。
• プルーニング・スキーム: 古典的な前処理コストと要求される量子レジスタのサイズとのバランスを取る、3つの実行可能なプルーニング手法（既定、クライロフ、およびハイブリッド）の開発と比較。
• スケーラビリティ: 特定の格子モデルにおいて、シャドウ・レジスタのサイズが物理システムのサイズが増加しても飽和することを示し、固定された量子ビット数で任意の大きさのシステムをシミュレートできることを実証した。
• 汎用性: 本フレームワークは、局所的な観測量（例：磁化）だけでなく、スピン電流自己相関関数やアウト・オブ・タイム・オーダー相関関数（OTOC）を含む高次の相関関数も扱うことができる。

結果
アルゴリズムは、1次元および2次元の格子スピンモデルでベンチマークが行われた：

• 1D 混合磁場イジングモデル (MFIM): 適度な横磁場を持つ系において、著者らは100量子ビットの物理系における磁化ダイナミクスを、わずか10量子ビットのシャドウ・レジスタで追跡できることを示した。ハイブリッド法を用いることで、12スピンの系に対して、既定の基底法と同等の精度を維持しながら、必要量を7量子ビットにまで削減できた。
• 2D MFIM: 本手法は、$4 \times 4$（16スピン）の系における磁化ダイナミクスを、14量子ビットのレジスタ（既定の基底）および7量子ビット（ハイブリッド法）を用いてシミュレートすることに成功し、物理的なシステムサイズと比較して大幅なリソース節約を示した。
• 縦磁場XXZモデル: 本アプローチは、スピン輸送（スピン電流自己相関関数の計算）および量子カオス（OTOCの計算）に適用された。ハイブリッド法により、クライロフ・ベクトルの効率的な格納が可能となり、メモリ・オーバーヘッドが削減された。
• 精度: 正確な対角化法との数値比較により、本手法が、弾道的領域における拡散係数のスケーリングや、カオス系におけるOTOCの特徴的な指数関数的成長を含む、ダイナミクスを正確に捉えていることが示された。

意義と主張
本論文は、シャドウ・ハミルトニアン・シミュレーションの主要なボトルネックを取り除き、相互作用系に対する実用的な量子アルゴリズムにすることを目的としている。状態を近似的に閉じられた代数上に符号化することで、本手法は、時間発展に必要な量子ビット・レジスタのサイズを劇的に削減する「スピン・ベクトル表現」を提供する。

著者らは、本アプローチが（ライトコーン伝播法とは異なり）動的なキャリブレーションを必要とせずにシャドウ・ハミルトニアンを構築する体系的な方法を提供すること、また（一部のテンソルネットワーク手法とは異なり）2次元および3次元モデルへ自然に一般化可能であることを強調している。主な目的は量子シミュレーションであるが、著者らは、古典的な前処理アルゴリズムも、結果として得られるシャドウ・ハミルトニアンの次元と構造に応じて、ダイナミクスのための古典的シミュレーション手法としても機能し得る旨を述べている。本研究は、特徴的な小さなパラメータを持つモデルに対しては、プルーニングが非常に効率的であり得ることを示唆しているが、その有効性は特定のモデルや観測量の選択に依存する。