技術要約：変分量子固有値ソルバーのためのリー代数生成子プールの最適化された構築

問題定義

変分量子固有値ソルバー（VQE）、特にADAPT-VQEのような適応型アルゴリズムは、量子アンザッツを構成するために使用されるパウリ演算子の集合である「オペレーター・プール」に依存しています。これらのアルゴリズムが厳密な基底状態への収束を保証するためには、プールは最小完全プール（Minimal Complete Pool: MCP）、すなわち、システムのハミルトニアンに関連する完全な動的リー代数（Dynamical Lie Algebra: DLA）を最小の基数で生成する演算子の集合である必要があります。

現在の手法における主要なボトルネックは、完全性の検証における計算量的な困難さです。従来の手法は、DLAの基底を明示的に構築するか、積群（product group）を列挙することを必要とし、これらは量子ビット数 $N$ に対して指数関数的（$O(4^N)$）にスケールします。この指数関数的なスケーリングにより、従来の手法は小規模なシステム（通常 $N \le 14$）への適用に限定されており、より大規模で複雑な分子への固定アンザッツ手法（NI-DUCC-VQEなど）の適用を妨げています。さらに、一般的なMCPは理論的には存在しますが、量子化学においては、分子ハミルトニアンに固有の対称性（粒子数、スピン、時間反転対称性など）を尊重しないことが多く、その結果、勾配の消失やアルゴリズムの停滞を招きます。

手法

著者らは、パウリ演算子の代数的な性質を利用し、それらをバイナリベクトル空間へとマッピングすることで、MCPを構築および検証するための一般的な多項式時間フレームワークを提案しています。

1. バイナリ行列への代数的マッピング

この手法は、パウリ文字列の集合 $A$ を直接分析する代わりに、それをバイナリ隣接行列 $\Gamma_A \in M_{|A| \times |A|}(\mathbb{F}_2)$ へとマッピングします。

• 定義: エントリ $\Gamma_A(i, j)$ は、パウリ演算子 $\alpha_i$ と $\alpha_j$ が反交換する場合に 1、そうでない場合に 0 となります。
• 意義: この行列はすべての交換関係をエンコードします。著者らは、「縮約（contraction）」（反交換グラフを簡略化するために使用されるグラフ理論的操作）が、$\mathbb{F}_2$ 上の $\Gamma_A$ に対する基本行/列操作に対応することを証明しています。

2. MCP合同定理（MCP Congruence Theorem）

中心となる理論的貢献は、MCP合同定理です。これは、ブラケット独立なパウリ文字列の集合 $A$ が、ターゲットとなるリー代数 $\mathfrak{g}$ に対するMCPであるための必要十分条件は、その行列 $\Gamma_A$ が $\mathbb{F}_2$ 上で標準的な行列 $\Gamma_{\text{canonical-MCP}}$ と**合同（congruent）**であることだと定めています。

• ランク条件: 対角成分がゼロである対称バイナリ行列（$\Gamma_A$ は常にこれに該当する）において、合同性は行列のランクのみによって決定されます。
• 検証: 完全性の検証は、$\text{rank}(\Gamma_A)$ がターゲットとなる代数の理論的ランクと一致するかどうかを確認することに簡略化されます。この検証は $O(N^3)$ の多項式時間でスケールし、従来の手法の指数関数的なスケーリングから劇的に改善されます。

3. 対称性を保持する構築

量子化学の特定のニーズに対処するため、本フレームワークは物理的な制約を取り入れています。

• 対称性の保持: プールは、分子の対称性（例：実数ウェーブファンクションを必要とする時間反転対称性による奇パウリ文字列、および粒子数保存）を尊重するように構築されます。
• ターゲットランク: 化学に関連する特定の部分代数（$\mathfrak{so}(2N)$ の部分代数）に関して、著者らはターゲットとなるMCPのランクが、一般的な $2N - 2$ ではなく $2N - 4$ であることを数値的に特定しました。
• ワークフロー: アルゴリズムは、「物理的に動機付けられたスターター（physically motivated starters）」（対称性を保持するよう事前選別されたUCCSDからの演算子）から開始します。ランク $\Gamma_A$ がターゲット（$2N - 4$）に達するまで、ブラケット独立な演算子を反復的に追加していきます。ランクに達しない場合は、指数関数的な積群を明示的に列挙することなく、探索空間を拡張します。

主な貢献

理論的進展

• 多項式時間の検証: 本論文は、ターゲットとなるリー代数に対する演算子プールの最小性と完全性を検証するための、初の証明可能な多項式時間（$O(N^3)$）の手法を提供します。
• 合同分類: 演算子プールの完全性とバイナリ行列の合同類との関係を定式化し、網羅的な探索ではなくランク評価を通じてMCPを特定することを可能にしました。
• 部分代数の特性評価: 著者らは、化学に関連する対称性制限された部分代数に関する特定のランク条件（$2N-4$）を特定し、抽象的なリー代数理論と実用的な分子シミュレーションの制約との間の溝を埋めました。

アルゴリズムの応用

• MB-ADAPT-VQE: 構築されたMCPを利用したバッチ型ADAPT-VQEのバリアントを導入しました。標準的なQEB（Qubit Excitation-Based）プール（$O(N^4)$ スケール）と比較して、プールのサイズを $O(N)$ に抑えることで、勾配評価のための測定オーバーヘッドを大幅に削減しました。
• スケーラブルなNI-DUCC-VQE: 本フレームワークは、構築のボトルネックにより以前は不可能であった大規模なシステム（H$_2$Oまで）のシミュレーションを成功させ、NI-DUCC-VQE（非反復型離散化ユニタリ結合クラスターVQE）の適用を可能にしました。

結果

本フレームワークは、Hyperionステートベクターエミュレータを用いて、一連の分子システム（LiH, H$_6$, H$_8$, H$_2$O）に対して検証されました。

1. 次元削減: 構築されたプールは、標準的なQEBプールと比較して、演算子数を1桁以上削減しました。H$_2$O（26量子ビット）の場合、プールサイズは、約15,000個の演算子から、48個（MCP）または667個（MCP + スターター）のコアへと削減されました。
2. 収束効率:
   • MB-ADAPT-VQE: 削減されたプールを用いることで、大きなバッチサイズ（$k=20, 30$）を用いた場合、標準的なADAPT-VQEと比較して、特に強相関系（引き伸ばされたH$_6$）や大きな分子（H$_2$O）において、関数評価回数の面で収束が大幅に加速されました。
   • NI-DUCC-VQE: 本手法は約1,500回の関数評価でH$_2$Oの化学精度を達成し、大規模シミュレーションへの実行可能性を実証しました。
3. スターターの重要性: 数値実験により、理論的な完全性は必要条件であるものの、物理的に動機付けられたスターターが極めて重要であることが確認されました。ランダムに選択された完全なプールや、（スターターのない）厳密に最小限のプールは、しばしば停滞や化学精度への到達失敗を招きました。これは、集合の「完全性」と同じくらい、演算子の「選択」が重要であることを浮き彫りにしています。

意義と主張

本論文は、リー代数生成子プールの構築を制限していた古典的な計算上のボトルネックを取り除く、基礎的なフレームワークを提供すると主張しています。検証問題を指数関数的な群の列挙から多項式時間の行列ランク評価へとシフトさせることで、著者らは以下を可能にしました。

• 以前は固定アンザッツ手法ではアクセス不可能であった、複雑な強相関系のシミュレーション。
• 測定負荷を軽減することによる、適応型VQEのリソース割り当ての効率化。
• 対称性適応型プールの完全性を保証するための数学的に厳密な手法の提供、および「対称性の壁ロードブロック」による勾配消失の防止。

著者らは、理論的フレームワークは完全性を保証するものの、実用的なVQEの収束においては、最小限のコアにさらなる物理的に動機付けられた演算子を加えることが有益であると述べています。彼らは、本研究を量子化学シミュレーションの領域内での主要な検証として位置づけつつも、量子エラー訂正、機械学習、ハードウェア制御への応用可能性を含む、より広い分野に応用可能な汎用的なツールとして提示しています。