Pauli-Sparse regularised Counterdiabatic Shortcuts for Linear-Ramp QAOA

技術要約：線形ランプ型QAOAのためのパウリ・スパース正則化反転対角（Counterdiabatic）ショートカット

問題提起
量子近似最適化アルゴリズム（QAOA）による組合せ最適化は、パラメータ効率と動的精度の間の決定的なトレードオフに直面している。線形ランプ型QAOA（LR-QAOA）は、層の角度を断熱的なスケジュールに従うよう制約することで、少数のパラメータで古典的な最適化の負担を軽減する。しかし、ハミルトニアンのスペクトルギャップが小さい、あるいは指数関数的に小さい場合、有限時間の断熱進化は強い非断熱遷移を引き起こし、システムが望ましい低エネルギー解の多様体から漏れ出す原因となる。これに対し、反転対角（CD）駆動は、断熱ゲージポテンシャル（AGP）から導出された補正項を加えることで、これらの遷移を抑制するための原理的な手法を提供する。しかし、標準的なアプローチは、有限のゲート予算下では実行不可能な、高密度な演算子や事前指定されたアンザッツを必要とすることが多い。さらに、困難な最適化問題における厳密なAGPは、低エネルギー多様体内の微視的かつ指数関数的に小さなエネルギー分裂を解明する必要がある場合が多く、これは計算コストが高く、潜在的に不要なこともある。

手法
著者らは、パウリ・スパースかつ正則化された反転対角拡張を備えたLR-QAOAのフレームワークを提案している。核となる手法は、完全な高密度リウヴィル行列を形成することなく、パウリ基底において正則化されたAGP方程式を直接解くことである。

1. 正則化されたAGP方程式： 本手法は、$(\mathcal{L}_H^2 + \eta I) A_\lambda^{(\eta)} = -i \mathcal{L}_H(G)$ という線形系を解く。ここで $\mathcal{L}_H(X) = [H, X]$ は交換子超演算子であり、$G = \partial_\lambda H$ である。$\eta > 0$ は正則化パラメータである。この方程式は、遷移のキャンセルとポテンシャルのノルムのバランスをとる、正則化された汎関数を最小化することから導出される。
2. エネルギー分解能メカニズム： パラメータ $\eta$ はエネルギー分解能スケールとして機能する。これは、エネルギー差 $|\omega| \ll \sqrt{\eta}$ に関連する遷移を抑制する一方で、$|\omega| \gg \sqrt{\eta}$ の遷移は保持する。これにより、特定の低エネルギー多様体内の指数関数的に小さな分裂を無視（厳密な断熱追従の必要性を回避）しつつ、高エネルギー励起状態への漏れを防ぐことが可能になる。
3. 不正確な共役勾配（CG）ソルバー： 著者らは、高密度なAGPを計算して圧縮する代わりに、パウリ座標における不正確な共役勾配法を採用している。
   • 行列フリー適用： 演算子 $(\mathcal{L}_H^2 + \eta I)$ は、パウリ交換子を用いて記号的に適用され、$4^n$ 次元の行列を構築することを回避する。
   • 切り捨て（Truncation）： CGイテレーション中、パウリ展開は、最大の係数を持つ $m$ 個の項を保持するように切り捨てられる（大きさによる切り捨て）。これにより、線形ソルバーは、実装可能なパウリ回転のサポート発見手順へと変換される。
   • 残差の再計算： 切り捨てによる誤差の蓄積を防ぐため、残差は再帰的に更新されるのではなく、各ステップで厳密に再計算される。
4. ガラーキン再適合（Galerkin Refit）と認証： 疎なサポートが発見された後、その特定のサポート上の係数を最適化するために、ガラーキン再適合が行われる。近似の質は、残差境界 $\|A - A_\lambda^{(\eta)}\|_{HS} \leq \|R\|_{HS} / \eta$ を用いて事後的に認証される。ここで $R$ は線形系の残差である。これは、厳密な解に関する知識を必要とせずに、計算可能な誤差範囲を提供する。
5. 回路への統合： 得られた疎なパウリ生成子は、LR-QAOA回路内の反転対角層として挿入される。変分最適化は、ランプ振幅 $(\Delta\beta, \Delta\gamma)$ と反転対角振幅 $\Delta_{CD}$ という削減されたパラメータ集合に対して行われる。

主な貢献

• 直接的なパウリ座標解法： 既存のCD-QAOA研究が、アンザッツ（局所項や入れ子状の交換子など）を規定したり、高密度なAGPを計算したりするのに対し、本研究は正則化されたAGP方程式をパウリ座標で直接解き、ポテンシャルの近似と実装可能な回転の疎なセットの選択を同時に行う。
• 物理的フィルタとしての正則化： 論文では、正則化パラメータ $\eta$ を単なる数値的な安定化手法としてではなく、特定のエネルギースケル以下の遷移を抑制する物理的なフィルタとして再定義している。これは、目標が単一の固有状態を追跡することではなく、低エネルギー多様体内に留まることである最適化問題において特に効果的である。
• 疎線形代数パイプライン： 著者らは、行列フリーの記号的交換子、サポート発見のための切断型CG、係数最適化のためのガラーキン再適合、および残差に基づく誤差認証を組み合わせた完全なパイプラインを開発した。このアプローチは、不正確なクリロフ部分空間法（Inexact Krylov methods）の概念を量子制御の文脈に導入したものである。
• 理論的保証： 論文では、正則化されたポテンシャルと未正則化のポテンシャルの間の距離、線形系の条件数、および厳密な正則化解に対する疎近似の誤差に関する厳密な境界を提供している。

結果
数値実験は、強磁性鎖（Ferromagnetic Chain: FC）のインスタンスおよび摂動を加えたFC–MaxCut/MarketSplitインスタンスに対して行われた。

• FCによる検証： FCインスタンスにおいて、標準的なLR-QAOAはシステムサイズが増大するにつれて性能が悪化し、秩序化された多様体に到達できなかった。提案されたLR-CD-QAOAは、厳密なAGPを用いることで、低深度でも高い秩序化多様体確率（$P_{ord} \approx 0.9$）を達成した。極めて重要なことに、本手法は様々な $\eta$ の値に対して堅牢であり、励起状態への漏れを防ぐために、指数関数的に小さな基底状態の分裂を解明する必要がないことを裏付けた。
• 疎CGの収束性： 実験により、少数のCGイテレーション（$K \approx 4$）と適度なパウリサポート予算（$m=64$）があれば、高密度なAGPの恩恵の大部分を回収できることが示された。近似比は、イテレーション回数とサポートサイズの両方に対して単調に改善した。
• 混合インスタンス： ハミルトニアンがFC成分によって支配されている摂動を加えたMaxCutおよびMarketSplitインスタンスにおいて、LR-CD-QAOAは補正のないLR-QAOAを大幅に上回った。LR-QAOAは深度 $p=8$ でも最適解から遠いままであったが、CD補正版は中程度の深度で厳密な最適解に近い近似比を達成した。

意義
本論文は、このフレームワークがQUBO最適化に対するQAOAの実用的な適用可能性を大幅に広げるものであると主張している。近接する低エネルギー構造や小さなスペクトルギャップを持つヘテロジニアスな問題のランドスケープにおける堅牢性を向上させることで、本手法は主要なボトルネックに対処している。これは、パラメータ効率的ではあるが動的に非効率な線形ランプの限界を超え、限られたゲート予算を持つ量子ハードウェア上で効果的な反転対角駆動を実装するための経路を提供する。本研究は、古典的な疎線形代数技術を用いて、問題のハミルトニアンの特定のスペクトル特性に動的に適応する量子回路を構築するという、ハイブリッド最適化パラダイムを確立している。