Low-gate-count block encodings for second-quantized fermionic Hamiltonians

技術要約：第2量子化フェルミオン・ハミルトニアンのための低ゲート数ブロック符号化

問題提起
多体ハミルトニアンの効率的なブロック符号化は、科学計算における量子アルゴリズム、特に初期の誤り耐性量子シミュレーションにおいて不可欠な前提条件である。既存の手法（線形結合としてのユニタリ演算子（LCU）アプローチなど）における主要なボトルネックは、古典的なハミルトニアン係数を量子回路にロードするためのリソースコストである。一般的な第2量子化電子構造ハミルトニアンの場合、相互作用項の数 $L$ はスピン軌道数 $n$ に対して $\Theta(n^4)$ でスケールする。標準的な構成では、Tゲートの複雑さが $L$ に対して線形（すなわち $O(n^4)$）となり、さらに劣正規化因子 $\alpha$ も $O(n^4)$ でスケールする。さらに、多くのアプリケーションでは粒子数 $\eta$ が固定されたシステムを必要とするが、標準的なブロック符号化は全ヒルベルト空間に拡散してしまうため、リソースのオーバーヘッドを削減するために粒子数保存を利用することができない。

手法
著者らは、フェルミオン演算子をパウリ文字列へマッピングするのではなく、占有数表現上で直接作用する、第2量子化フェルミオン・ハミルトニアンの新しい明示的な構成法を提案している。この手法は、主に2つのオラクル、すなわち疎性オラクル（$O_C$）と振幅オラクル（$O_A$）に基づいている。

1. 疎性オラクル ($O_C$): 行列要素の位置を決定するためにマルチ制御ゲートを使用する代わりに、著者らはSWAPベースのデータルックアップ戦略を採用している。彼らは SWAP-UP ネットワークを利用して、占有されている特定の軌道を代表する量子ビットを最小次位ビットへと移動させ、効率的な制御ビット反転と位相符号化を可能にしている。このアプローチにより、SWAP-UP、ラダー回路（スタッガードCNOT）、およびフェルミオンの符号を扱うためのターゲット・パウリZゲートを用いて、列疎性オラクルを構築する。
2. 振幅オラクル ($O_A$): 係数の数値値をエンコードするために、著者らは SELECT-SWAP データルックアップ回路と直接サンプリング法を組み合わせている。SELECT-SWAP技術は、古典データのロードに伴うTゲート数を $O(L)$ から $\tilde{O}(\sqrt{L})$ へと削減するためにデータ項目をグループ化する。直接サンプリング法は、係数が $O(1)$ であるか、あるいはインデックスとともに減衰する場合において、複雑な算術演算やエイリアスサンプリングを必要とせずに、係数のバイナリ表現を確率振幅へと変換する。
3. 固定-$\eta$ 粒子構成: 粒子数 $\eta$ が固定されたシステムに対して、著者らは間接拡散戦略を導入している。すべての形式的な軌道インデックスに対して拡散するのではなく、拡散は入力状態 $|j\rangle$ に条件付けられ、占有されている軌道に対応するインデックスのみに対して重ね合わせを生成する。これには、第 $i$ 番目の占有量子ビットのアドレスを補助レジスタ内で一時的に計算するための占有検出オラクル（$O_{occ}$）が必要となる。この状態依存型のアプローチは、拡散を関連する部分空間に限定し、劣正規化因子を大幅に減少させる。
4. 構造化されたハミルトニアン: 本フレームワークは、並進不変性、近接相互作用、または係数の減衰を持つハミルトニアンへと拡張される。これらの場合、ロードすべき異なる係数の数は減少する。間接拡散は、すべての可能なペアの代わりに、有効なインデックスのペア（例：$|p\rangle|p \pm m\rangle$）の重ね合わせを生成するための算術回路と組み合わされる。

主な貢献と結果

• ゲート複雑性の二次的改善: 一般的な第2量子化ハミルトニアン（相互作用項 $L$）に対し、提案された構成は $\tilde{O}(\sqrt{L})$（具体的には $\tilde{O}(\sqrt{L \log(L/\epsilon)})$）のTゲート数を達成し、標準的なLCU/QROM手法の線形 $O(L)$ スケーリングに対して二次的な改善を実現している。電子構造ハミルトニアン（$L \sim n^4$）の場合、これによりTカウントは $O(n^4)$ から $\tilde{O}(n^2)$ へと減少する。
• 固定-$\eta$ セクターにおける劣正規化因子の削減: 全ヒルベルト空間の設定では劣正規化因子は依然として $O(L) \sim O(n^4)$ であるが、固定粒子数 $\eta$ の部分空間を標的とした構成では、一般的な二体相互作用に対して、この因子を $O(L)$ から $O(n^2 \eta^2)$ へと減少させる。$\eta \ll n$ の場合、これは量子化（qubitization）のコストにおいて大幅な削減となる。
• 最適化されたリソース見積もり:
  • ハバードモデル: 2つのユニークな係数のみを持つ標準的なハバードモデルでは、データルックアップのコストは消失し、複雑さは疎性オラクルによって支配される。これにより、全空間に対して $O(n + \log(n/\epsilon))$、$\eta$ 粒子セクターに対して $O(n \log \eta + \log(\eta/\epsilon))$ のTゲートが得られる。
  • 並進不変および局在化モデル: 並進不変性または局在化した軌道（係数が減衰するケース）を持つハミルトニアンでは、異なる係数の数 $L_{distinct}$ が減少する。本構成はこれらの構造に適応し、非構造化の場合と比較してTカウントとTデプスをさらに低下させる。
• 既存モデルとの比較: 本論文は、その結果をBabbushら[28]のLCUモデルおよびKivlichanら[49]のトロッターベースのアプローチと比較している。新しい手法は、量子化ベースのシミュレーションにおいて、より低いTカウントとTデプスを提供する。注目すべき点として、関連する固定粒子入力モデル（例：Tong, Anら[38]）が存在するものの、著者らの構成は、ハミルトニアンのシフトやインスタンス依存の古典的最適化に頼ることなく、回路構成自体を通じて最適な劣正規化スケーリングを達成し、一般的な第2量子化ハミルトニアン（非対角の一次および二体相互作用を含む）の固定-$\eta$ セクターを直接標的とする最初の明示的なブロック符号化である。

意義
著者らは、本研究が多体系の初期誤り耐性量子シミュレーションへの実用的な道筋を提供すると主張している。その意義は以下の点にある：

1. 大幅なリソース削減: 項の数（$L$）に対するTゲート複雑性の線形から平方根への減少、および固定粒子システムにおける劣正規化因子の削減は、時間発展および固有値問題のシミュレーションにおける誤り耐性オーバーヘッドを直接的に低下させる。
2. 直接的な部分空間標的化: 本構成は、一般的な第2量子化ハミルトニアンの固定-$\eta$ セクターを直接標的とする最初の明示的なブロック符号化であり、ハミルトニアンのシフトや古典的最適化ではなく、回路構成を通じて非対角相互作用に対する最適な劣正規化スケーリングを達成している。
3. スケーラビリティ: 本手法はシステムサイズ $n$ および粒子数 $\eta$ に対して効率的にスケールするため、標準的な入力モデルよりも大きな分子および凝縮系システムをシミュレートするのに適している。
4. モジュール性: 本フレームワークは汎用的であり、近接相互作用、並進不変、因子分解などの様々なハミルトニアン構造に適応可能である。これは、物理的な対称性や局所性を活用することで、コストを系統的に削減できることを示しており、さらなるコスト削減のためにブロック符号化の枠組みに統合できることを示している。

本論文は、改善が最も顕著なのは構造化されたハミルトニアンや固定粒子セクターであるが、一般的な構成においても、SWAP-UPネットワークとSELECT-SWAPデータルックアップを活用することで、Clifford+Tゲートのオーバーヘッドを最小限に抑え、従来のメソッドに対して既に大きな優位性を示していると結論付けている。