Dynamical spectral functions from bitstring-sampled quantum subspaces: entanglement, not one-body magic, tracks the sampling cost

技術要約：ビットストリング・サンプリングによる量子部分空間からの動的スペクトル関数

問題提起
サンプリングに基づく量子対角化（SQD）および量子選択的構成間相互作用（QSCI）は、主に基底状態エネルギーを標的とした、近未来の電子構造計算における主要な量子アルゴリズムとして台頭している。しかし、基底状態エネルギーに関する古典的なフロンティアはすでに閉ざされている。ヒューリスティックな古典的手法（例：マッチゲート・シミュレータ、自己回帰型ニューラル・サンプラー、熱浴構成間相互作用など）は、この特定の観測量において、量子サンプリングによる部分空間と同等またはそれ以上の性能を示す。したがって、基底状態の化学における一般的な指数関数的な量子優位性の主張は、ますます維持が困難になっている。

本論文は、決定的なギャップを特定している。すなわち、グリーン関数に必要とされる荷電セクター（$N \pm 1$）において、ビットストリング・サンプリングされた量子部分空間から、周波数分解された一粒子スペクトル関数（$A(\omega)$）や動的構造因子（$S(q, \omega)$）を再構成する手法が公開されていないことである。さらに、どのリソースが実際にサンプリングコストを決定しているのかについても明確さが欠けている。自然な仮説は、「一粒子フェルミオンの魔法（one-body fermionic magic）」（非ガウス性の尺度）が困難さを駆動しているというものだが、これは未証明であり、誤解を招く可能性がある。

手法
著者は、計算基底の測定のみを用いて、SQD/QSCIのツールキットを静的なエネルギーから連続的な周波数分解された動力学へと昇華させる、統一されたプリミティブを提案している。

1. サンプリング・プリミティブ： 本手法は、アダマール・テスト、制御ユニタリ、またはオーバーラップ推定を用いない、浅い実時間発展回路を利用する。
   • シード（種）： 一粒子スペクトルの場合、システムは $(N \pm 1)$ セクターにおける応答シード $|\phi^\pm\rangle = \hat{c}^\dagger_p |0\rangle$ および $\hat{c}_p |0\rangle$ を準備する。中性の動的構造因子の場合は、数保存的な密度演算子またはスピン演算子 $\hat{n}_q |0\rangle$ または $\hat{S}^z_q |0\rangle$ をシードとする。
   • 発展： これらのシードは、時間発展（$e^{-i\hat{H}t_k}$）させ、計算基底で測定されることでビットストリングを生成する。
   • 部分空間の構築： 発展時間を通じてサンプリングされた構成の和集合が、行列式部分空間 $\mathcal{S}$ を形成する。
2. 古典的再構成： 量子プロセッサはビットストリングのみを提供する。ハミルトニアンは $\mathcal{S}$ 上に投影され、古典的に対角化される。レマン表現を用いて、後退グリーン関数とスペクトル関数が組み立てられる：
   $$G_{pq}(\omega) = \sum_m \frac{\langle 0|\hat{c}_p|\tilde{m}^+\rangle \langle \tilde{m}^+|\hat{c}^\dagger_q|0\rangle}{\omega - (\tilde{E}^+_m - E_0) + i\eta} + \dots$$
   これにより、$A(\omega)$、$k$ 分解された $A(k, \omega)$、および動的構造因子 $S(q, \omega)$、$S^{zz}(q, \omega)$ が得られる。
3. ノイズ緩和： 本手法は、ノイズの多いビットストリングを正しい粒子数セクターへと投影し、デバイスノイズ下での部分空間の質を向上させるために、自己整合的な構成回復（S-CoRe）を採用している。

主な貢献

1. 動的スペクトルの再構成： 本論文は、ビットストリング・サンプリングされた部分空間から、周波数分解された一粒子スペクトル関数および動的構造因子を再構成することを初めて示した。これは、SQD/QSCIのパラダイムを離散的なエネルギー準位から連続的な動的観測量へと拡張するものであり、ハバード鎖および19種類の分子を用いた厳密な対角化によって検証されている。
2. 魔法とコストのデカップリング： 著者は、一粒子フェルミオンの魔法（フェルミオン的反平坦性 $F_1 = 4 \text{tr}[\gamma(1-\gamma)]$ によって定量化される）が、サンプリングコストからデカップリングされていることを厳密に証明した。
   • $F_1$ は軌道回転不変（ガウス不変）であるが、サンプリングコスト（行列式サポート $|S|$）は基底に依存する。
   • $F_1$ は、静的な多参照相関（結合解離中に上昇するなど）を忠実に診断するが、サンプリングコストを予測するには失敗する。例えば、ハバード鎖においてオンサイト斥力を増大させると、$F_1$ は上昇するが、局在化により $|S|$ は減少する。
   • 真のコスト要因は、基底依存的なもつれ（entanglement）、具体的には最小ボンド次元 $\chi$（テンソルネットワークのボンド次元）であると特定された。これは経験的に $|S|$ と高い相関（$\rho \approx 0.90$）を示す。
3. リソース診断： 本研究は、真の量子優位性が存在するならば、それは固定されたサンプリング基底における高ランクかつ非ガウス的な高次累積量の中に存在することを確立した。

結果

• 検証： 本手法は、ハバード鎖および19種類の分子に対する厳密なスペクトル関数を再現しており（全構成間相互作用との比較において相対-$L_1$誤差が $10^{-3}$ 未満）、精度は $10^{-5}$ Ha を下回っている。
• ハードウェア実証： アルゴリズムは、IBM Heron プロセッサ（$L=6$ ハバードモデルには $ibm\_fez$、$伸長した$N_2$には$ibm_marrakesh$）上で実行された。
  • $L=6$ ハバードのスペクトル関数は、50,000 ショットから再構成され、広がり分解能の範囲内で厳密な結果と一致した。
  • 伸長した $N_2$ については、ノイズを補完する構成回復により、基底状態エネルギー誤差 $0.59 \pm 0.14$ mHa（化学的精度以下）を達成し、同じ浅い回路のノイズレスなシミュレーションを上回った。
• スケーリングと比較：
  • サンプリングコスト（必要なセクターの割合）は、固定された精度に対してシステムサイズとともに指数関数的に減少するが、絶対的な行列式サポート $|S|$ は指数関数的に増大する。
  • 時間発展セレクターは、強結合領域（$U/t \ge 12$）において、決定論的な古典的 Krylov および CIPSI セレクターよりも優れた性能を示すが、アクセス可能なサイズにおいては熱浴 CI を凌駕することはない。
  • ハバード鎖と（一致した一粒子魔法を持つ）二脚ラダーとの制御された比較では、ラダーの方が高いボンド次元 $\chi$ のために、大幅に大きなサンプリング分画を必要とすることが示された。これは、コストが $F_1$ ではなく $\chi$ に追従することを裏付けている。
• AI/ML の限界： 「テイル・コンプリーター（裾の補完）」として学習された生成モデル（GFlowNet）を用いてデバイスのサンプルを補完したが、安価な古典的 CIPSI ベースラインを統計的に上回ることはできず、現在の生成モデルは、この文脈において量子サンプリングを代替する段階にはないことが示唆された。

意義と主張

本論文は、電子構造における量子優位性の議論を再定義している：

• ターゲットの転換： 著者は、正当な量子優位性は、基底状態エネルギー（古典的手法が追いついている領域）ではなく、材料探索や相関電子物理において中心的な役割を果たすが、1次元を超える次元では古典的テンソルネットワークにとって困難な、動的観測量（スペクトル関数、構造因子）にあると主張している。
• リソースの明確化： 一粒子魔法がサンプリングの困難さを保証するという誤解を正している。代わりに、正確なリソースに関する記述を提供している。すなわち、一粒子魔法は「相関の種類」（多参照性）の診断指標であるが、コストは基底依存的なもつれ（ボンド次元 $\chi$）によって決定される。
• 実用的な有用性： 本手法は、アンシラや複雑なゲートを必要としない、ハードウェアネイティブなアプローチを提供する。これは、現在のエネルギーのみを扱うハードウェアでは提供できない観測量を計算するためのものである。これは「トリアージ（選別）」ツールとして機能する。効率的に測定可能なフェルミオンの魔法は、平均場理論が失敗する強相関ターゲットをフラグ立てできるが、量子プロセッサを使用するかどうかの決定は、一粒子魔法では予測できない基底依存的なコスト（ボンド次元）を推定することに基づかなければならない。

著者は、本手法がテストされた特定の事例において（DMRGが1次元や狭いラダーに対して依然として有効であるため）「古典超越的」な優位性を主張するものではないとしつつも、量子サンプリングが必要となる境界線（自由フェルミオン推定や低ボンド次元のテンソルネットワークの手が届かない、固定サンプリング基底における高ランク・非ガウス領域）を明確に定義することに成功したと結論づけている。