Spectral Gaps via Imaginary Time

技術要約：虚時間によるスペクトルギャップの推定

問題設定
スペクトルギャップ（$\Delta E = E_1 - E_0$、基底状態と第一励起状態のエネルギー差と定義される）は、蛍光波長や材料における励起子エネルギー移動から、断熱的な実行時間や状態転送の忠実度に至るまで、多様な物理的文脈において極めて重要なパラメータである。このギャップを決定することは計算上困難であり、最悪の場合、局所ハミルトニアンがギャップを持つかどうかの判定は決定不能である。また、低エネルギー部のスペクトルを推定することは一般にNP困難である。既存の数値的手法（ファームイオン量子モンテカルロ法は符号問題に悩まされ、テンソルネットワーク法はボンド次元に制限される）は、システムサイズ、基底サイズ、および測定オーバーヘッドに関して大きなトレードオフに直面している。虚時間発展（ITE）は基底状態を見つけるための堅牢な古典的ツールであるが、全状態密度の分解や量子位相推定を行わずに特定のスペクトルギャップを抽出することは、依然として特有の課題である。

手法
著者らは、虚時間軌道に沿って評価された入れ子状の交換子の期待値の比から、スペクトルギャップ $\Delta E$ を直接推定する手法を提案している。

1. 理論的枠組み:
   有限次元のハミルトニアン $H$ とそのスペクトル分解 $H = \sum E_n \Pi_n$、および基底（$\Pi_0$）と第一励起部分空間（$\Pi_1$）の両方にサポートを持つ初期状態 $|\phi_0\rangle$ について、系は虚時間 $\tau$ に沿って $|\psi(\tau)\rangle = e^{-\tau H}|\phi_0\rangle$ と進化する。
   本手法は、基底および第一励起部分空間を結合する「調整用」観測量 $O$ を利用する。著者らは、$M$ 次の入れ子状の交換子 $[H, O]_M$ を定義する。
   核となる結果（定理1）は、$\tau \to \infty$ のとき、期待値の比がエネルギーギャップの二乗に収束することを示している：
   $$\frac{\langle [H, O]_{M+2} \rangle_\tau}{\langle [H, O]_M \rangle_\tau} = (E_0 - E_1)^2 + \mathcal{O}(e^{-\tau(E_2 - E_1)})$$
   このエスティメータは、個々の固有ベクトルではなく射影形式に基づいているため、高次励起レベルの縮退に対して非依存である。

2. 対数エスティメータ:
   また、本論文では、$E_0 + E_1$（$\ln|\langle [H, O]_M \rangle_\tau|$ 対 $\tau$ の傾きから）および第二のギャップ $E_2 - E_1$（比のフィッティングの残差から）を抽出するための対数エスティメータを導出しており、数値的な小さな数の除算が不安定な場合の代替案を提供している。

3. 量子実装 (ITQDE):
   ユニタリ発展がネイティブである量子ハードウェア上でこれを実装するために、著者らは虚時間量子力学的エミュレーション（ITQDE）を利用している。彼らは、ユークリッド核 $e^{-\tau H^2}$ を実時間ユニタリのガウス平均として表現するために、ハバード・ストラトノビッチ変換を用いている。この核はガウス・エルミート求積法を用いて離散化され、発展をユニタリの加重和として近似する。
   期待値は、$\langle \phi_0 | U_j^\dagger P_\alpha U_k | \phi_0 \rangle$ という形のオーバーラップの実部の加重和へと簡約される（ここで $P_\alpha$ はパウリ文字列である）。これらのオーバーラップは、標準的なアダマールテストを用いて測定できる。著者らは、核の対称性により、これらのオーバーラップの実部のみが必要であることを指摘している。

主な結果
著者らは、QuSpinライブラリを用いた古典シミュレーションを通じて、2つの標準的なモデルで本手法のベンチマークを行った。

• 1次元横磁場イジングモデル (TFIM): システムサイズ $L=4$ において、推定されたギャップの相対誤差は、数値的なノイズ以下に第一励起状態が抑制される前に、$\tau$ とともに指数関数的に減少した（$\sim 10^{-8}$ に到達）。
• 1次元フェルミ・ハバードモデル: 半充填の鎖（$L=4$）において、同様の誤差の指数関数的減少が観察された（相対誤差は約 $\sim 10^{-4}$）。本手法は、非自明なスピンおよび電荷構造を持つ相互作用のあるフェルミオンの設定を正常に扱った。
• 量子シミュレーション: QiskitのAerバックエンドと $10^6$ ショットを用いて、TFIMの原理証明シミュレーションが行われた。本手法はITQDEフレームワークと互換性があったが、達成された最小誤差は約 $\sim 10^{-1}$ であった。著者らは、これを「帯域幅のトレードオフ」に起因すると考えている。すなわち、高次状態を抑制するために必要なガウスフィルタが、$\tau$ と求積次数 $K$ が注意深く調整されていない場合、関連するオーバーラップをショットノイズの床へと押し下げてしまうためである。

意義と主張
本論文は、スペクトル推定における「控えめながらも真の構造的利点」を提供すると主張している。

• 簡潔性と範囲: 本手法は、個々の固有状態を分解したり、量子位相推定を行ったり、全状態密度を再構成したりすることなく、低次のギャップを推定する方法を提供する。これは、スペクトルの問題を、一パラメータの流れに沿った緩和の問題へと変換するものである。
• 堅牢性: このエスティメータは、高エネルギーレベルの縮退に対して堅牢であり、古典的な疎行列演算と量子的な和の項形式のハミルトニアンの両方に適用可能である。
• 謙虚さ: 著者らは、このアプローチが局所ハミルトニアン問題の最悪ケースの複雑性を回避するものではないことを明示している。これは、詳細な状態密度が必要な場合の、完全なスペクトル再構成や位相推定の代わりとなるものではない。
• 今後の展望: 本研究は原理証明としての役割を担っている。著者らは、量子ハードウェア上で本手法を競争力のあるものにするには、パラメータ（$\tau, K$）を最適化するために先行研究からの完全な帯域幅分析を導入し、より優れた調整用観測量を設計し、より測定効率の高いオーバーラップエスティメータを採用する必要があると述べている。本論文は、このギャップに焦点を当てた視点を、緩和率やギャップが主要な量となる、微分方程式を解くための有効ハミルトニアンを用いる広範なプログラムと並行して活用できるツールとして位置づけている。