Block Encoding Non-Abelian Lattice Gauge Theory

技術要約：非アーベル・格子ゲージ理論のブロック符号化

問題提起
量子コンピュータ上で、量子色力学（QCD）のような非アーベル格子ゲージ理論をシミュレーションすることは、ハミルトニアンの磁気セクターにおいて大きな課題となっている。電気的部分のハミルトニアンは既約表現（irrep）基底において対角的で単純であるが、磁気部分は、複数のリンクおよびサイトの自由度の間で複雑かつ相関した遷移を生成するプラケット演算子を伴う。ギブンス回転法に依存する従来のアプローチは、打ち切り（truncation）に伴って行列要素が爆発的に増加し、物理的に妥当な打ち切りレベルでは古典的な事前計算やストレージが不可能になるという「スケーリングの壁」に直面している。既存のリソース見積もりは、解析的なカットオフ・スケーリング公式、あるいは計算コストの高いインレジスタ・クレブシュ・ゴルダン（CG）算術を含む不完全な構成に依存していることが多い。

手法
本論文では、疎な符号化フレームワーク、特に「既約表現（irrep）」または「簡約電気（reduced electric）」基底を利用した、プラケット演算子（$\square_p + \square_p^\dagger$）の効率的なブロック符号化オラクルを提案する。この基底はゲージ不変であり、量子ビット効率が高く、内部のカラー指数をリンクの既約表現ラベルとサイトのシングレット多重度指数に置き換えることで、内部カラー指数を排除している。

手法の核心は、プラケット行列要素の因子分解特性を利用することにある。行列要素を単一のモノリシックな値として扱うのではなく、行列要素が4つの「サイト因子」（プラケットの各頂点に対応）と既約表現次元の事前係数に因子分解できることを利用する：
$$\langle t | \square_p | s \rangle = \left( \prod_{k=1}^4 \sqrt{\frac{\dim R_{\ell_k}}{\dim R'_{\ell_k}}} \right) \prod_{v=1}^4 S_{v}$$
ここで、$S_v$ は、頂点 $v$ におけるアクティブなリンク既約表現、制御リンク既約表現、およびシングレット多重度指数に依存する。

構成は以下のステップで行われる：

1. 古典的事前計算： 公開ツール（例：pyclebsch）を用いて、サイト因子を事前計算し、それらを回転角のシーケンスに変換する。決定的なのは、制御クラスのルックアップ（ADDR回路経由）を導入し、制御リンクのパターンを等価類にグループ化することである。この圧縮により、量子ルックアップテーブル（QROM）に必要な異なるエントリの数を大幅に削減できる。
2. 提案（Proposal）ステージ： 「提案」オラクルが、ソースとなるリンクコアに基づき、ターゲットとなる状態と方向（順方向または逆方向）を提案する。このステージでは、回転ツリーを使用して、ソースに対する最大可能な振幅で重み付けされた、潜在的な遷移の重ね合わせを作成する。
3. コーナー（Corner）ステージ： 4つの独立したコーナーステージが提案を精緻化する。ソースのデータ、制御クラス、および提案されたターゲットに基づき、これらのステージは振幅を回転させて正確なサイト因子に一致させる。提案された遷移が局所的なガウス則または制御リンクと矛盾する場合、振幅は「拒絶（reject）」フラグへと回転される。
4. ブロック符号化の組み立て： 完全なオラクル $U_p$ は以下によって構築される：
   • 重み付き準備（$PREP$）を適用して、隣接状態の重ね合わせを作成する。
   • ターゲットコアとソースコアを入れ替える（演算子の非対角性を実装するため）。
   • 方向ビットを反転させる。
   • アジュイント準備（$PREP'^\dagger$）を適用してアンシラをアンコンピュートし、正しい行列要素の符号を正確に一度だけ適用することを保証する。
5. ルーティング： 単一プラケット・オラクルは、プラケット・サブレジスタをオラクルの位置へとスワップするルーティング回路を通じて、格子全体のブロック符号化へと埋め込まれる。

主な貢献

• 明示的な構成： 本論文は、インレジラCG算術を必要とせずに、既約表現基底における磁気プラケット演算子のブロック符号化オラクルに関する完全かつ明示的な構成を提供している。
• 因子分解と圧縮： 行列要素の因子分解を利用し、制御クラスの等価性を導入することで、行列要素を直接テーブル化する場合と比較して、ルックアップテーブルのエントリ数を数桁単位で削減した。
• リソース見積もり： 2つの特定の打ち切り（$B=6$ および $B=9$。ここで $B$ は頂点あたりの全電気カシミール・カットオフ）における、SU(3) ゲージ理論の具体的なTゲート数と論理量子ビット数を見積もった。
  • $B=6$ の場合、プラケット・オラクルの1回の呼び出しに約 $10^5$ Tゲートを要する。
  • $B=9$ の場合、コストは約 $10^6$ Tゲートに上昇する。
  • 論理量子ビット数は $Q \approx 12L^3 + O(10^3)$ と推定され、格子状態レジスタによって支配される。
• 代替手法との比較： 本論文は、これまでのギブンス回転法を直接拡張した場合、膨大な数の行列要素により、これらの打ち切りレベルでは少なくとも $10^5$ 倍高いコストが発生することを指摘している。

結果と意義
主要な結果は、「プラケット行列要素のスケーリングの壁」が、因子分解と古典的事前計算を通じて緩和可能であることを示したことである。構築されたオラクルは、量子特異値変換（QSVT）のようなフォールトトレラント・アルゴリズムにとって管理可能なサブノーマライゼーション係数（$\alpha_p \approx 56$ at $B=9$）を持つ、非アーベル・ゲージ理論のシミュレーションを可能にする。

本論文は、ゲージ・マターセクターに関する最近の研究を補完する「磁気部分」を提供し、フォールトトレラントなQCDシミュレーションの完全なリソース見積もりを可能にした点に意義があると主張している。また、現在の構成は行列要素の指数関数的な増大を緩和しているものの、大規模な格子に対する総リソースコストは依然として体積スケーリング（純粋なQSVTルーティングの場合 $L^6$）に支配されており、これが将来のアルゴリズム開発における重要な未解決問題であることを明らかにしている。著者は、自身の研究がノイズのあるハードウェア上でのデモンストレーションではなく、フォールトトレラント領域におけるアルゴリズム設計とリソース見積もりの演習であることを強調している。

今後の方向性
本論文は、以下の将来の研究の方向性を提案している：

• 連続体極限に到達するために、ルックアップテーブルのさらなる最適化（例：ハッシングやインレジスタ算術）を用いた、打ち切り限界の拡張（例：$B=16$ への拡張）。
• スケーリングを $L^6$ から $L^{4.5}$ へ改善するための、ハイブリッド・トロッター-QSVTアルゴリズムの調査。
• 同様の因子分解アプローチを用いた、マターセクター（クォーク）の統合。
• 提案された既約表現ベースのブロック符号化と、ループ・ストリング・ハドロン（LSH）基底やオービフォールド・アプローチとの比較。