Sublinear-depth Quantum Simulation of Electrons with Atomic Orbitals

技術要約：原子軌道を用いた電子の劣線形深度量子シミュレーション

問題提起
分子や材料における電子構造のシミュレーションは、デジタル量子コンピューティングの主要な応用分野である。トロッター積公式（Trotter product formulas）と量子化（qubitization）フレームワークの両方において、量子シミュレーションアルゴリズムの進化時間（$t$）および精度（$\epsilon$）に関する性能は十分に理解されているが、構造化が進んでいない系におけるシステムサイズ（$N$）に対するスケーリングは依然として課題となっている。

分子中の電子をモデル化するための標準的な手法は、原子が意味のあるサブユニットとして残るコンパクトで柔軟な表現を提供する原子軌道を利用する。しかし、電子の自由度を原子軌道へと離散化すると、最大で$O(N^4)$個の相互作用項を持つハミルトニアンが生じ、固有の構造が乏しくなる。これらすべての項を直接実装して量子コンピュータ上でシミュレートしようとする過去の試みは、非現実的なほど大きな多項式スケーリング（例：$O(N^{10})$）をもたらした。線形結合によるユニタリ演算（LCU）やハミルトニアン分解を用いても、従来の成果では、原子軌道は大規模シミュレーションにおいて$O(N^{5/3})$程度のゲート複雑度を実現する平面波法や実空間格子法に対して競争力がないことが示唆されていた。さらに、既存の効率的な手法の多くは量子化に依存しており、これはしばしば$\lambda \ge \|H\| = \Omega(N)$となるクエリ深度$\Omega(\lambda t)$を伴うため、劣線形深度のシミュレーションを阻んでいた。

手法
本研究では、劣線形深度を実現し、クリーンな補助量子ビット（ancilla qubits）を必要としない、原子軌道のためのトロッターベースの量子シミュレーションアルゴリズムを提案する。この手法は、以下の3つの核心的な技術的柱に基づいている。

1. 厳密な軌道の局在化と基底構築：
   著者らは、現実的な多原子系が満たすべき「軌道の局在化」の概念を正式に定義している。標準的な原子軌道は直交していないという課題があり、それらを直交基底（例：ハートリー・フォック基底）に変換すると空間的な局所性が失われることが多い。本論文では、ロドウィン直交化（Löwdin orthogonalization, $S^{-1/2}$）が、システムが熱力学的極限（原子あたりの最小軌道解像度が一定、かつ原子間距離の最小値が一定）にある条件下において、指数関数的な減衰特性を保持することを証明している。これにより、結果として得られる直交基底も、明確な減衰パラメータを持つ「原子軌道セット」であり続けることが確立される。

2. 階層的ハミルトニアン分解（高速多重極展開法に着想を得た手法）：
   軌道のオーバーラップの指数関数的な減衰を利用して、著者らはまずハミルトニアンを、重要な相互作用を行う軌道ペアが$O(N)$個となる疎なグラフへと切り詰める。その後、残りの長距離クーロン相互作用（$h_{text{long}}$）を、古典的な**高速多重極展開法（FMM）**に着想を得た戦略を用いて分解する。

   • システムは、ボックス（レベル $j$）を用いて階層的に分割される。
   • 相互作用は、短距離成分（直接処理）と長距離成分に分離される。
   • 空間的に離れたボックス間の長距離相互作用は、多重極展開（球面調和関数）を介して近似され、その結果、ハミルトニアン項の低ランク分解が行われる。
   • この分解により、長距離相互作用は、互いに素な量子ビット集合上で作用するランク1の演算子の和へと削減される。

3. 浅いトロッターステップのための量子回路：
   アルゴリズムは、これらの分解された項を処理することで各トロッターステップを実装する。

   • フェルミオン置換： ジョルダン・ウィグナー符号化において相互作用する軌道を整列させるため、フェルミオン置換を適用する。アルゴリズムは、補助量子ビットなしで$O(\log^2 N)$の深度でこれらの置換を実行するバージョンを利用する（$\Theta(N)$個の補助量子ビットを用いれば$O(\log N)$の深度が可能だが、提案手法はこのオーバーヘッドを回避する）。
   • ランク1のZZ相互作用： コアとなる革新は、ランク1のZZ相互作用（$e^{-ivwt}$）の進化を実装するための新しい量子回路である。著者らは、このような相互作用が、浅い量子フーリエ算術回路と少数の**「ダーティ（dirty）」な補助量子ビット**を用いることで、多項式対数深度（$O(\text{polylog}(Nt/\epsilon))$）で実装可能であることを示している（これらの補助量子ビットは、アイドル状態のシステム量子ビットから再利用可能であり、正味のクリーンな量子ビット数の増加を必要としない）。
   • 並列化： レベル内の相互作用を定数個のレイヤーに並列化できるように階層的分解を構成することで、1つのトロッターステージあたりの総深度を低く抑える。

主な貢献と結果
本論文は、以下の具体的な結果を提示している。

• 理論的基礎： ロドウィン直交化が原子軌道の指数関数的な減消を保持することを証明し、熱力学的極限における局在化された直交基座の使用を正当化した。
• アルゴリズムの効率性： 1ステップのトロッターステップを深度$O(\text{polylog}(Nt/\epsilon))$で実装するトロッターベースのシミュレーションアルゴリズム。
• 複雑度の境界：
  • ゲート数： 3次元系では全ゲート数が$N^{5/3 + o(1)}$（2次元系では$N^{3/2 + o(1)}$、1次元系では$N^{1 + o(1)}$）のスケーリングを示す。これは、あらゆる基底において達成可能な最良の知られているスケーリング（以前は平面波法でのみ達成されていた）と一致する。
  • 深度： シミュレーションの深度は、3次元系において$N^{2/3 + o(1)}$という劣線形スケーリングを示す。
  • 補助量子ビット： アルゴリズムはゼロのクリーンな補助量子ビットを要求し、システムから再利用される「ダーティ」な補助量子ビットのみを使用する。
• 誤差解析： 交換子構造を用いてトロラー誤差を抑え込み、最適なトロッター次数が上述のゲート複雑度および深度の複雑度を与えることを示している。

意義
本論文は、量子化学シミュレーションにおける長年の非効率性を解決したと主張している。厳密な軌道局在化、FMMに着想を得た分解、および浅い算術回路を組み合わせることで、著者らは、原子軌道が大規模な量子シミュレーションにおいて競争力を持つことを、従来の仮定に反して示した。

その意義は以下の点にある：

1. 深度の壁を打破： 以前の手法が線形またはそれ以上の深度スケーリングに制限されていたのに対し、原子軌道シミュレーションにおいて劣線形深度（$N^{2/3}$）を達成した。
2. リソース効率： 近傍または初期の誤り耐性デバイスにおいて希少なリソースであるクリーンな補助量子ビットを、再利用可能な「ダーティ」な補助量子ビットを用いることで排除した。
3. パラダイムシフト： 特定の指標において量子化よりも効率的とされることが多いトロッターベースの手法が、高度なハミルトニアン分解技術と組み合わせることで、深度および補助量子ビットの要件において量子化を凌駕し得ることを示した。

著者らは、漸近的なスケーリングは改善されているものの、定数因子や対数オーバーヘッドを考慮すると、本手法が他と比較して優位になるクロスオーバーポイントは特定のシステムサイズやハードウェアの制約に依存するため、さらなる研究が必要であると注記している。また、彼らのアプローチは低エネルギー部分空間への制限とは異なり、完全な電子構造問題に対処するものであることを明確にしている。