On the generalized eigenvalue problem in subspace-based excited state methods for quantum computers

量子コンピュータを用いた励起状態計算において、重なり行列の条件数が大きい場合、統計的サンプリング誤差により標準的な一般化固有値問題に基づく手法（QSE や qEOM）が不安定になるのに対し、固有値問題として定式化される q-sc-EOM 法は誤差に対してより安定であり、より適した候補であると示されています。

要約

🧐 背景：量子コンピュータと「分子の歌」

まず、量子コンピュータは、分子の動きやエネルギーをシミュレーションする「未来の計算機」として期待されています。
分子には、普段の安定した状態（基底状態）だけでなく、エネルギーを吸収して高ぶった状態（励起状態）があります。これを計算することは、新しい薬や太陽電池の開発に不可欠です。

しかし、現在の量子コンピュータは「ノイズ（雑音）」が多く、完全ではありません。この雑音の中で、いかに正確に分子の「歌（エネルギー状態）」を聞き取るかが課題でした。

⚔️ 対決：2 つの戦士

この論文では、励起状態を計算する 2 つの主要な方法（アルゴリズム）を比較しました。

1. QSE（量子部分空間拡張）: 従来の有力な方法。
2. q-sc-EOM（量子自己無撞着運動方程式）: 比較対象の新しい方法。

これらを「歪んだ鏡」と「完璧な鏡」に例えてみましょう。

1. QSE：歪んだ鏡（Generalized Eigenvalue Problem）

QSE という方法は、計算をする際に**「重み付けされた鏡（重なり行列）」**を使います。

• 仕組み: 鏡に映った像（分子の状態）を計算する際、鏡自体が少し歪んでいたり、重みがついていたりします。
• 問題点: 量子コンピュータには「ノイズ（雑音）」が必ず混入します。この歪んだ鏡にノイズが乗ると、**「歪みがノイズを何倍にも増幅させてしまう」**という現象が起きます。
  • 例え話: 小さな音（ノイズ）を、歪んだメガホンで増幅すると、耳を劈くような大音量の雑音になってしまいます。
  • 結果: 鏡の歪み（条件数）が大きいと、計算結果がバラバラになったり、計算そのものが破綻（特異点）して答えが出せなくなったりします。

2. q-sc-EOM：完璧な鏡（Standard Eigenvalue Problem）

一方、q-sc-EOM という方法は、**「歪みのない、完璧な鏡（単位行列）」**を使います。

• 仕組み: 鏡の歪みという要素を最初から排除しています。
• メリット: ノイズが混入しても、歪んだ鏡のように増幅されません。
  • 例え話: 完璧な鏡に小さな傷（ノイズ）がついても、映る像はそれほど乱れません。
  • 結果: 計算が非常に安定しており、少ない計算回数（ショット数）でも正確な答えが得られます。

---

📊 実験結果：3 つのシナリオ

研究者たちは、分子（H4 や NH3）を使って、鏡の歪みの度合い（条件数）を変えながら実験を行いました。

① 歪みが小さい場合（低条件数）

• 状況: 鏡がほとんど歪んでいない時。
• 結果: QSE も q-sc-EOM も、どちらもよく似ていました。どちらも正確に計算できました。

② 歪みが中くらいの場合（中条件数）

• 状況: 鏡が少し歪んできた時。
• 結果:
  • QSE: 急激に不安定になりました。ノイズが増幅され、計算結果が大きくブレます。正確な答えを出すために、10 倍もの計算回数が必要になりました。
  • q-sc-EOM: 全く影響を受けず、安定して正確な結果を出し続けました。

③ 歪みがひどい場合（高条件数）

• 状況: 鏡が極端に歪んで、割れそうになっている時。
• 結果:
  • QSE: 計算が完全に破綻しました。答えが出ません。
  • 無理やり直す方法（しきい値法）: 壊れかけの鏡の「歪んだ部分」を切り取って計算を強行しました。
    • 代償: 計算はできましたが、「一部の分子の状態（励起状態）が見えなくなりました」。まるで、鏡の割れた部分に映っていた人物が突然消えてしまったような状態です。化学研究では「見えない状態」は致命的な問題です。
  • q-sc-EOM: 歪みがないため、すべての状態を正確に、安定して計算し続けることができました。

---

💡 この研究が伝えるメッセージ

この論文の結論は非常にシンプルで重要です。

「量子コンピュータで分子の励起状態を計算するなら、QSE よりも q-sc-EOM の方が、ノイズに強く、信頼性が高い」

• QSE の弱点: 計算の過程で「歪んだ鏡（重なり行列の逆行列）」を使うため、量子コンピュータ特有のノイズが増幅されやすく、計算が不安定になります。
• q-sc-EOM の強み: 「歪みのない鏡」を使うため、ノイズの影響を最小限に抑えられ、必要な計算回数も少なく済みます。また、重要な「見えない状態（欠落した励起状態）」の問題も起きません。

🎯 未来への展望

この研究は、量子コンピュータが実用化される「ノイズの多い時代（NISQ 時代）」において、**「どのアルゴリズムを選ぶべきか」**という指針を示しました。

化学者や研究者にとって、**「正確なスペクトル（すべての状態）」**を知ることは、新しい材料や薬の開発に不可欠です。q-sc-EOM のような、ノイズに強く、すべての状態を捉えられる方法は、将来の量子化学計算において、より有望な候補であると言えます。

一言でまとめると：
「歪んだ鏡（QSE）で雑音の中で計算するのは危険ですが、完璧な鏡（q-sc-EOM）を使えば、雑音に負けないで正確な分子の姿が見える！」という発見です。

技術概要

論文の技術的概要：量子コンピュータにおける基底状態励起状態法のための一般化固有値問題の解析

この論文は、量子コンピュータを用いた量子化学、特に励起状態の計算において、サブスペースベースの手法（QSE, qEOM, q-sc-EOM）が統計的サンプリング誤差（ショットノイズ）に対してどのように耐性を持つかを、理論的および数値的に詳細に分析したものです。

1. 問題提起 (Problem)

量子コンピュータを用いた励起状態計算において、量子サブスペース展開（QSE）や量子運動方程式（qEOM）などの主要な手法は、古典コンピュータ上で一般化固有値問題（$H\mathbf{C} = S\mathbf{C}E$）を解く必要があります。ここで、$H$ はハミルトニアン行列、$S$ は重なり行列（オーバーラップ行列）です。

• 核心課題: 量子コンピュータでは測定に統計的誤差（ショットノイズ）が避けられません。重なり行列 $S$ の条件数（$\kappa(S)$）が大きくなると、$S$ の逆行列を計算する過程で誤差が劇的に増幅され、固有値の精度が著しく低下します。
• 深刻な状況: 条件数が非常に高い場合、行列が特異（singular）になり、数値的に解けなくなります。これを回避するために「しきい値処理（thresholding）」を行うと、一部の励起状態が欠落する（スペクトルが不完全になる）という重大な問題が発生します。

2. 手法と理論的枠組み (Methodology)

著者らは、以下の 2 つの手法を比較対象として分析しました。

1. QSE (Quantum Subspace Expansion):
   • 基底状態波動関数に励起演算子を作用させてサブスペースを構築。
   • 一般化固有値問題 $H\mathbf{C} = S\mathbf{C}E$ を解く。
   • $S$ の逆行列計算が必要であり、誤差増幅のリスクが高い。
2. q-sc-EOM (Quantum Self-Consistent Equation-of-Motion):
   • 自己無撞着な演算子を用いた運動方程式理論に基づく手法。
   • 基底関数の直交性により、重なり行列 $S$ が解析的に単位行列 ($I$) となる。
   • したがって、解くべき式は通常の固有値問題 $H\mathbf{C} = \mathbf{C}E$ となり、$S$ の逆行列計算が不要。

理論的解析:
摂動論を用いて、標準固有値問題と一般化固有値問題における固有値の誤差増幅を定式化しました。

• 標準固有値問題: 誤差はショットノイズ $\epsilon$ に比例 ($O(n\epsilon)$)。
• 一般化固有値問題: 誤差は $\frac{n(1+|\lambda|)}{\lambda_{\min}(S)} \epsilon$ に比例。
  • ここで $\lambda_{\min}(S)$ は重なり行列の最小固有値です。
  • 条件数 $\kappa(S)$ が大きい（$\lambda_{\min}(S)$ が小さい）場合、誤差が $1/\lambda_{\min}(S)$ 倍に増幅され、ショットノイズの影響が極めて大きくなります。

3. 主要な貢献 (Key Contributions)

1. 誤差増幅の定量的診断: ショットノイズ、条件数、および励起状態の信頼性を結びつける実用的な診断基準を提示しました。
2. 中程度の条件数領域での不安定性の示唆: 条件数が「中程度」（例：100〜1000 程度）であっても、QSE ではショットノイズにより大きな不安定性が生じることを数値的に実証しました。
3. しきい値処理の限界の明確化: 高条件数領域でしきい値処理を用いて解を安定化させることは可能ですが、その代償として励起状態の一部が欠落し、スペクトルの完全性が損なわれることを示しました。
4. q-sc-EOM の優位性の立証: 直交基底（$S=I$）を採用する q-sc-EOM は、重なり行列の条件数に依存せず、ショットノイズに対してロバストであることを示しました。

4. 結果 (Results)

• 低条件数 ($\kappa(S) < 100$):
  • H4 分子の例において、QSE と q-sc-EOM の両方でショットノイズに対する耐性は同程度でした。
• 中条件数 ($100 < \kappa(S) < 10^4$):
  • H4 分子（$\kappa(S) \approx 592.5$）および NH3 分子（$\kappa(S) \approx 124 \sim 427$）の計算において、QSE は誤差と分散が急激に増大しました。
  • 同等の精度を得るために、QSE は q-sc-EOM よりも約 10 倍（1 桁）多くのショット数を必要としました。
  • 一方、q-sc-EOM は条件数に関わらず安定した結果を提供しました。
• 高条件数 ($\kappa(S) > 10^4$):
  • H4 分子（$\kappa(S) \approx 2.05 \times 10^{12}$）において、QSE は通常の計算では解けなくなりました。
  • しきい値処理を適用して解いた場合、解は得られましたが、特定の励起状態（縮退状態を含む）が欠落していることが確認されました。
  • q-sc-EOM はしきい値処理を必要とせず、すべての励起状態を安定して取得しました。

5. 意義と結論 (Significance)

• 量子化学計算への示唆: 量子コンピュータを用いた励起状態計算において、重なり行列の条件数による誤差増幅は避けられない課題です。QSE や qEOM のような一般化固有値問題を解く手法は、分子幾何構造や系によって精度が不安定になるリスクがあります。
• 手法の選択: 励起状態のスペクトル完全性（すべての状態を捉えること）と精度が求められる化学研究（例：二重励起状態、コニカル交差、励起状態ダイナミクス）においては、q-sc-EOMのような、重なり行列が単位行列となる手法が、ショットノイズに対してより頑健で、より適した候補であることが示されました。
• 将来的な展望: この研究は、一般化固有値問題を用いる他の量子アルゴリズム（量子クリロフ法など）においても同様の課題が存在することを示唆しており、誤差耐性の高いアルゴリズム設計の重要性を強調しています。

要約すると、本論文は「一般化固有値問題における重なり行列の条件数が、量子計算の統計的誤差を増幅させる主要な要因であり、これを回避する q-sc-EOM のような手法が、実用的な量子化学計算においてより安定した選択肢である」ことを理論的・数値的に証明したものです。