Matrix Product Operators In The Age of Block Encoding

技術要約：ブロックエンコーディング時代における行列積演算子（MPO）

問題提起
量子シミュレーションアルゴリズム、特に線形結合によるユニタリ（LCU）およびブロックエンコーディングの枠組みに基づく手法は、ハミルトニアンの時間発展をシミュレートする際に、しばしば重大なリソースオーバーヘッドに直面する。従来のアプローチでは、時間発展演算子 $e^{-itH}$ をパウリ文字列の明示的な和（テイラー展開）へと展開する。局所的なスピン鎖において、サイト数 $N$ および相互作用の次数 $T=O(N)$ とすると、この明示的な展開は、次数 $K$ の多項式に対して $O(N^K)$ という組合せ爆発的なパウリ項の増加をもたらす。さらに、まずハミルトニアン $H$ を行列積演算子（MPO）としてブロックエンコードしてから量子固有値変換（QET）を適用するという代替手法は、しばしば指数関数的な正規化ペナルティを招く。このペナルティは、非ユニタリなハミルトニアンMPOをブロックエンコードするために必要な局所的なユニタリ拡張において、正規化因子 $\alpha_{dil}(H)$ がシステムサイズに対して指数関数的に増大する（例：$\sim 1.43^N$）ことに起因しており、成功確率を著しく低下させ、過剰な振幅増幅を必要とする。

手法
著者らは、「コンパイル済み多項式（compiled-polynomial）」アプローチを提案している。これは、MPOを単にブロックエンコードされるべき圧縮されたオブジェクトとしてではなく、量子線形代数のための構造化された中間表現（IR）として扱うものである。コアとなる手法は、以下の3つの明確なステップで構成される：

1. 圧縮されたLCUとしてのMPO: 著者らは、MPOを「仮想パス（virtual paths）」上の圧縮されたLCUプログラムとして解釈する。各局所テンソルを個別にユニタリ化する代わりに、各局所テンソルを局所ユニタリ基底（例：パウリ基底）で展開する。MPOの仮想ボンドは、局所演算子のシーケンスを生成する有限状態オートマトンとして解釈される。これにより、MPOを、各パスが特定のパウリ文字列と係数 $c_\gamma$ に対応する、仮想パス $\gamma$ に関する和として捉えることが可能になる。
2. 条件付きコンパイル（Conditional Compilation）: 著者らは、MPO構造から直接LCU回路を構築するブロックエンコーディング・コンパイラを開発した。標準的なLCU回路（$PREP^\dagger SELECT PREP$）を、一連の条件付きステージへと一般化している。
   • 条件付き準備（CPREP）: $N$ 個の条件付き状態準備が行われる。各ステージでは、入力される仮想ボンドのラベル $a_{n-1}$ に依存して、局所的な演算子のラベル $\mu_n$ と出力される仮想ボンドのラベル $a_n$ を準備する。確率は、MPOテンソルから導出された「後方メッセージ（backward messages）」$\beta_n$ を用いてエンコードされ、これによりグローバルな分布がMPOの係数と一致するように制御される。
   • 条件的選択（SELECT）: 準備されたパスラベルに基づいて、局所ユニタリが適用される。
   • 正規化: 正規化因子 $\alpha_{MPO}$ は、絶対値転送行列のテンソル縮約を通じて効率的に計算される。これにより、すべてのパスを明示的に合計する必要がなくなる。
3. 多項式合成: ターゲットとなるユニタリ $e^{-itH}$ は、まずMPO表現内において多項式 $P_\chi(t)$ として近似される（テイラー展開を使用）。この多項式MPOが、その後直接ブロックエンコードされる。これは、先に $H$ をブロックエンコードしてからQETを適用するという標準的な経路とは対照的である。

主な貢献

• 仮想パスLCU: 著者らは、MPOを仮想パス上の圧縮されたLCUプログラムとして解釈する手法を導入し、テンソルネットワークのデータ構造からブロックエンコーディング・オラクルを組み立てる体系的な方法を提供した。
• コンパイラ・アーキテクチャ: 条件付き準備および選択ステージを用いて、MPOデータ構造を量子回路へと翻訳する新しいコンパイラを提示した。これは、テンソルネットワークを量子回路のための有効な中間表現として扱うものである。
• ボトルネックの回避: 著者らは、多項式変換を（ブロックエンコーディングの前段階で）MPOレベルで コンパイルすることが、以下の2つの主要なボトルネックを回避できることを示した：
  1. 明示的なパウリ文字列の組合せ爆発 $O(N^K)$。
  2. 非ユニタリなハミルトニアンMPOを直接ブロックエンコードすることに伴う指数関数的な正規化ペナルティ $\alpha_{dil}(H)$。

結果
著者らは、ハイゼンベルク模型およびその摂動版（サイト依存場および二量化を含む）を用いて、システムサイズ最大 $N=64$ まで、提案手法の数値的検証を行った。

• ボンド次元のスケーリング: ハミルトニアンの冪乗 $H^K$ に対して、必要なMPOボンド次元 $\chi$ は緩やかに増大する（例：誤差 $10^{-6}$ において $K=4$ で $\chi=24$）。これはシステムサイズ $N$ に依存しない。これは、圧縮されていないスケーリングである $5^K$ よりも大幅に低い。
• 正規化: コンパイルされたユニタリMPOの正規化因子 $\alpha_{MPO}(P_\chi)$ は、大きな $N$ や $K$ に対しても、小さな定数倍の範囲内に留まる。これは、$H$ を先にブロックエンコードしてMPO拡張を行うルートで見られる指数関数的な増大とは極めて対照的である。
• 複雑性: 全体の回路コストのプロキシは $C_3 \sim \alpha_{MPO}(P_\chi) N \chi^2$ とスケーリングする。テストされた領域において、これは明示的なパウリ文字列LCU展開の $O(N^K)$ スケーリングと比較して、多項式的な高速化を実現している。具体的には、場が摂動を与えたモデルにおいて $\chi=4$ のとき、コストは $N^{1.17}$ とスケーリングし、二量化モデルでは $N^{1.35}$ となり、線形スケーリングに近づいている。

意義および主張
本論文は、テンソルネットワークを単なるコンパクトな表現としてではなく、コンパイラ解析のための構造化されたデータとして扱うことが、量子アルゴリズムを加速させる新たな道を切り開くことを主張している。主な意義は、ハミルトニアンを多項式MPOへと古典的に前処理してからコンパイルするという「操作の順序」が、リソース要件を劇的に変化させ得ることを示した点にある。

著者らは、本手法が明示的なパウリ文字列の増大と、従来のMPOベースのブロックエンコーディング戦略における指数関数的な正規化ペナルティの両方を回避することで、多項式的な高速化を提供すると断言している。結論として、このアプローチは、MPOボンド次元とパスの正規化が緩やかな、近似的にユニタリな実時間伝搬関数に対して特に効果的であるとしている。また、オラクルコスト、反復回数、および成功確率のスケーリングを包括的に分析する必要があるため、量子化（qubitization）手法との完全な比較は依然として未解決の課題であるとも述べている。さらに、将来の研究として、このコンパイラをマルチプロダクト公式（multiproduct formulas）と組み合わせることや、非ユニタリ関数へ適用することの可能性についても示唆している。