Certifying fermionic Gaussian states (and a little more) with optimal precision dependence

技術要約：最適な精度依存性を持つフェルミオン・ガウス状態の認証

問題提起
スレーター行列式やバーディン・クーパー・シュリーファー（BCS）対状態を含むフェルミオン・ガウス状態（fGS）は、フェルミオン系の量子シミュレーションにおける基本的な参照状態として機能する。これらの状態は量子ハードウェア上で効率的に準備できるが、デバイスが特定のターゲットとなるfGSを正しく準備したことを検証することは依然として困難である。既存の認証プロトコルは、サンプル複雑性、測定の適応性、および認証可能な状態のクラスの間でトレードオフに直面している：

• 忠実度推定（例：忠実度ウィットネス）は堅牢であるが、$O(\epsilon^{-2})$ のサンプルを必要とし、認証としては最適ではない。
• 量子状態検証（QSV） プロトコルは $O(\epsilon^{-1})$ の最適なスケーリングを実現できるが、多くの場合、もつれ測定を必要とするか、あるいは「ハール典型的な」状態に限定されており、fGSのような構造化された族をカバーできていない。
• オラクルベースのプロトコルは、任意の純粋状態を $O(\epsilon^{-1})$ のスケーリングで認証できるが、多数の量子ビット（$d-1$ 個）に対する適応的な測定と、任意の積基底における確率オラクルへのアクセスを必要とする。

中心となる問いは、「すべての純粋フェルミオン・ガウス状態に対して、単一量子ビット測定のみを使用し、効率的な古典的処理を維持しながら、同時に最適な精度スケーリング（$O(\epsilon^{-1})$）を達成できるか？」である。

手法
著者らは、完全な完備性（perfect completeness）を伴う量子状態検証（QSV）の枠組み内で動作する、適応型認証プロトコル（アルゴリズム1）を提案している。このプロトコルは単一量子ビット測定を利用し、フェルミオン・ガウス状態の特定の代数的構造に依拠している。

1. プロトコルの構造： 各ラウンドにおいて、プロトコルは以下の2つのテストのいずれかをランダムに選択する：

   • パリティ・テスト： すべての $d$ 個の量子ビットを計算基底（$Z$ 基底）で測定し、状態が正しいフェルミオン・パリティセクターにあることを検証する。
   • 相関テスト：
     • ランダムなモードのペア $\{k_1, k_2\}$ を選択する。
     • 残りの $d-2$ 個の量子ビットを計算基底で測定する。
     • 量子ビット $k_1$ を $X$ 基底で測定する。
     • 前の測定結果に基づく、最後の量子ビット $k_2$ の条件付き状態を計算する。
     • この条件付き状態を含む基底で量子ビット $k_2$ を測定する。
   • プロトコルは、すべてのテストに合格した場合のみ受理される。

2. 理論的枠組み：

   • プロトコルの効率は、検証演算子 $\Omega$ のスペクトルギャップ $\gamma(\Omega)$ によって支配される。
   • 著者らは、検証演算子を、$d \leftrightarrow d-2$ のダウンアップ・ウォーク（2サイト・グレイバー動力学に相当）である認証ウォークと呼ばれるマルコフ連鎖へと写像している。
   • スペクトルギャップは、ターゲット状態の計算基底分布 $\pi(x) = |\langle x|\psi\rangle|^2$ を分析することによって下限が抑えられる。
   • 決定的なのは、純粋なfGSに対して、$\pi(x)$ の生成多項式が**開いた右半平面において零を持たない（zero-free）**という性質を利用している点である。この性質は、状態の共分散行列（$\Gamma^T\Gamma = I$）のパラフィアン構造から導かれる。これは、高次元エキスパンダーの理論におけるエントロピー収縮引数を通じて、認証ウォークの急速な混合を意味している。

主な貢献と結果

1. 最適なサンプル複雑性： 本プロトコルは、不完全度（infidelity）$\epsilon$ および信頼度 $1-\delta$ で、$d$ モードの純粋フェルミオン・ガウス状態を以下の数のコピーを用いて認証する：
   $$N(\epsilon, \delta, d) \le \left\lceil \frac{d(d-1)}{\epsilon} \ln \frac{1}{\delta} \right\rceil$$
   $\epsilon^{-1}$ への依存性は最適であり、もつれ測定を用いたプロトコルのスケーリングと一致している。システムサイズに対する $O(d^2)$ の依存性は、従来の非適応型または推定ベースの手法による $O(d^4)$ や $O(d^3)$ のスケーリングに対する大幅な改善である。

2. 測定および処理の効率性：

   • 測定： 単一量子ビット測定のみが必要である。極めて重要な点として、適応的な測定が必要なのはたった1つの量子ビット（相関テストにおける最後の量子ビット）のみである。
   • 古典的処理： 1コピーあたりの古典的処理時間は、条件付き共分散行列とパラフィアンの計算によって支配される $O(d^3)$ である。

3. 境界のタイトさ：

   • 著者らは、$O(d^2)$ の境界が、特定の物理的に関連のある状態、すなわち完全二量体化されたSu–Schrieffer–Heeger (SSH) チェーンの基底状態およびBCS対状態に対してタイトであることを証明している。これらの状態は、スペクトルギャップの下限 $\gamma(L) = 2/[d(d-1)]$ を満たす。
   • $d=14$ までの数値シミュレーションは、「ハール典型的（ランダム）」なガウス状態に対して、スペクトルギャップが $O(d^{-1})$ でスケールすることを示唆しており、これはサンプル複雑性が $O(d\epsilon^{-1})$ であることを意味する。これは、ワーストケースの境界よりも因子 $d$ だけ優れた結果である。

4. フェーズ・ドレスト状態への拡張：

   • 本プロトコルは、ガウス状態に効率的に計算可能な対角フェーズ $U_\theta = \sum_x e^{i\theta(x)}|x\rangle\langle x|$ を適用して得られる**効率的にフェーズ・ドレストされた状態（efficiently phase-dressed states）**へと拡張される。
   • 認証ウォークは計算基底分布にのみ依存するため（これは $U_\theta$ に対して不変である）、フェーズ・ドレストされた状態を用いてもサンプル複雑性は変わらない。
   • このクラスには、4モードの $(|0000\rangle + |0011\rangle + |1100\rangle + e^{i\theta}|1111\rangle)/2$ のような、マッチゲート計算のための非ガウス「マジック状態」が含まれる。

意義と主張
本論文は、最適な精度スケーリングと最小限の測定適応性を用いて、すべての 純粋フェルミオン・ガウス状態を認証するという未解決問題の解決を主張している。高次元エキスパンダーの理論とマルコフ連鎖の混合時間の理論を、認証問題と結びつけることで、著者らはガウス状態の代数的性質（具体的には生成多項式の零フリー特性）と検証効率の間の厳密な関係を確立した。

本研究は、適応性とサンプル複雑性のトレードオフにおける中間領域に本プロトコルを位置づけている。すなわち、 $O(d^2\epsilon^{-1})$ のスケーリングを実現するために単一の適応量子ビットを使用し、$O(d^4\epsilon^{-2})$ のコストがかかる非適応型推定や、$O(d\epsilon^{-1})$ のコストがかかる $d-1$ 個の適応的測定を必要とするプロトコルを回避している。著者らは、ワーストケースの境界は特定の物理的状態に対してタイトであるが、ランダムな状態に対する平均的な挙動はより良好であることから、典型的なシナリオにおけるさらなる最適化の可能性を示唆している。フェーズ・ドレスト状態への拡張は、フェーズ関数が効率的に計算可能である限り、本プロトコルの効率性が広範な非ガウス的修正に対して頑健であることを示している。