Variationally Optimized Imaginary-time Polynomial Filters for Ground State Projection

技術要約：量子基底状態準備のためのアンシラを用いた虚時間発展の変分最適化

問題提起
量子多体系の基底状態（$|GS\rangle$）を決定することは、平衡特性を予測する上で極めて重要である。厳密対角化（ED）のような古典的手法はヒルベルト空間の指数関数的なスケーリングによって制限され、他の古典的アプローチ（例：量子モンテカルロ法、テンソルネットワーク）は符号問題やアンザッツの表現力といった課題に直面するが、量子コンピュータはこれらのボトルネックを回避するためのネイティブな経路を提供する。量子基底状態準備における主要な課題は、非ユニタリ演算子 $e^{-\tau H}$ で定義される虚時間発展（ITE）を実装することである。この演算子は決定論的なユニタリゲートとして直接実現できないため、既存の確率的実装は、切断されたテイラー展開とトロッター・鈴木（TS）分解を組み合わせた固定係数近似に依存することが多い。これらの標準的な手法は、近似誤差を減らすには時間ステップを小さくするか展開次数を高くする必要があるという、回路の深さが増し、決定的なことに事後選択（post-selection）の成功確率が低下するという重大なトレードオフに直面する。このトレードオフは、近未来の量子デバイスにおける ITE の実用的な効率を制限している。

手法
著者らは、アンシラを用いた確率的 ITE を最適化するための「作用に基づく変分アプローチ（ABVA）」を提案する。このフレームワークは、事後選択を介して線形結合（LCU）を実現するために単一のアンシラクビットを使用するという標準的な回路アーキテクチャを維持しつつ、固定された解析的係数を、演算子レベルの作用原理から導出された変分最適化パラメータに置き換えるものである。

1. 変分フレームワーク： 本手法は、作用 $S = \int d\tau \text{Tr}[U^\dagger(\dot{U} + HU)]$ に基づいている。この作用を変分アンザッツに関するパラメータで極値化することにより、係数の運動方程式を導出する。
2. 変分多項式アンザッツ（ABVPA）： 虚時間更新において、標準的な一次テイラー展開（$1 - \tau H$）は、一般化された線形多項式 $c_0(\tau) + c_1(\tau)H$ に置き換えられる。係数 $c_0$ および $c_1$ は、ハミルトニアンの静的モーメント（$\mu_n = \text{Tr}[H^n]$）に依存する、量子幾何テンソルおよび一般化された力によって決定される。
3. 変分トロッター・鈴木アンザッツ（ABVTSA）： 同様に、LCU回路内の実時間プロパゲータ・ブロック（通常、二次トロッター・鈴木分解 $e^{-iA\delta t/2}e^{-iB\delta t}e^{-iA\delta t/2}$ を通じて実装されるもの）は、変分形式 $e^{id_0(t)A}e^{id_1(t)B}e^{id_2(t)A}$ に置き換えられる。パラメータ $d_i(t)$ は、正確なプロパゲータに対する誤差を最小化するように最適化される。
4. 実装： 最適化された係数は、アンシラの準備角度および制御ユニタリ・ブロック内の回転パラメータを修正する。極めて重要な点は、この最適化は回路の深さやゲート数を変更せず、パラメータの値のみを変更することである。
5. 誤差境界： 著者らは、線形フィルタリング（ABVPA）と実時間分解（ABVTSA）の両方について、厳密で状態に依存しない誤差境界を導出した。これらの境界は、ハミルトニアンの静的特性（トレースモーメント）のみに依存し、精度と成功確率のバランスをとるための最適な時間ステップ $\delta\tau$ を選択するための古典的な事前計算ガイドラインとなる。

主な貢献

• 作用に基づく変分最適化： ITE フィルタおよび実時間プロパゲータの係数を、回路の複雑さを増すことなく最適化する変分フレームワークの導入。
• 安定性と精度の向上： 変分パラメータを用いることで、標準的なテイラー展開が失敗する領域においても、大きな時間ステップを維持しながら安定性を保てることを実証した。
• 成功確率の改善： フィルタを最適化することで、ステップごとおよび累積の事後選択成功確率を高め、目標とする統計的精度を達成するために必要な回路ショット数を減少させた。
• リソース効率： 本手法は単一アンシラの要件を維持し、マルチアンシラ方式のオーバーヘッドを回避しつつ、固定係数法よりも優れた性能を提供する。

結果
本手法は、$N=5$ および $N=8$ 量子ビットの一次元横磁場イジングモデル（QIM）を用いてベンチマークを行った。

• 収束性： シミュレーションにおいて、ABVA 実装は標準的な TS-Taylor ITE よりも速く、かつ滑らかに基底状態エネルギーへと収束した。標準的な手法は、大きな時間ステップ（例：$\tau=1.8$）において偏差やオーバーシュートを示したが、ABVA は安定し、有効性を維持した。
• 忠実度（Fidelity）： ABVA は、特に後期時間領域において、正確な基底状態に対する高い忠実度を達成した。忠実度の分布は、標準的な手法と比較して、より狭く 1 に近い値を示しており、これは近似誤差に対する感度が低いことを示している。
• 成功確率： ステップごとの成功確率は、一貫して ABVA の方が高かった。その結果、全進化にわたる累積成功確率が大幅に向上し（一部の領域では最終的な成功確率において最大 1 桁の向上）、測定のオーバーヘッドを削減した。
• 誤差解析： 数値評価により、選択された時間ステップ（$\delta\tau = 0.4$）が変分近似の妥当な領域内にあり、グローバル誤差が標準的な TS-Taylor アプローチよりも大幅に低いことが確認された。

意義と主張
本論文は、ABVA フレームワークが、近未来の量子デバイスにおける高忠実度な基底状態準備への実用的な経路を提供すると主張している。静的な解析的近似からパラメータ最適化へと負担を移すことで、本手法は近似精度と事後選択成功確率の間に存在する固有のトレードオフを緩和する。著者らは、これが回路の深さやアンシラのオーバーヘッドを増やすことなく達成されることを強調しており、これはリソース効率の高い、既存のゲートベースの確率的 ITE スキームに対する改良であることを示している。本研究は、変分多項式フィルタリングが、厳密なエネルギー・フィルタリング手法と現在の量子プロセッサのハードウェア制約との間の溝を埋めることができる可能性を示唆している。著者らは、スケーラビリティに関して慎重な姿勢を維持しており、本手法は精度と確率のバランスを改善するものの、小さなスペクトルギャップを持つ系や複雑なハミルトニアンにおいては、依然として時間ステップの注意深いチューニングや高次の近似が必要となる可能性があり、それらは追加のリソース・トレードオフを伴う可能性があると述べている。