Solution of the Equation-of-Motion Phonon Method eigenvalue problems on the D-Wave quantum annealer

本論文は、D-Wave量子ハードウェア上で式に基づいたフォノン法（Equation-of-Motion Phonon Method）による大規模な固有値問題を反復的に解くために、量子アニーリングと古典的なデフレーションを組み合わせたハイブリッド量子・古典アルゴリズムを提示しており、核多体系理論における近未来の量子デバイスの可能性と現在の限界の両方を実証している。

要約

あなたは、非常に巨大で、信じられないほど複雑なパズルを解こうとしているところだと想像してください。核物理学の世界において、このパズルとは、原子核の中で陽子と中性子（核子）がどのように振る舞うかを解明することです。このパズルの「ピース」は、ハミルトニアン行列と呼ばれる巨大な数学的格子の中に配置されています。原子核が大きくなればなるほど、ピースの数は増え、その格子はあまりにも巨大になるため、世界最速のスーパーコンピュータでさえ、妥当な時間内にすべての解（固有値と固有ベクトル）を見つけ出すことが困難になります。

本論文は、これらのパズルを解くための新しい方法を提示しています。それは、古典的なコンピュータと、量子アニーリング（具体的にはD-Waveマシン）と呼ばれる特殊なタイプの量子コンピュータを提携させるという手法です。

以下に、彼らのアプローチを簡単な比喩を用いて解説します。

1. 問題点：解の連なる山脈

原子核のエネルギー状態を、広大で霧に包まれた山脈だと考えてください。

• 目標： あなたは、最も低い谷（基底状態）を見つけ、次にそれ以外のすべての谷や峰（励起状態）を順番に見つけ出したいと考えています。
• 従来の方法（古典的コンピュータ）： 従来のアルゴリズムは、一歩一歩、注意深く確認していくハイカー（登山者）のようなものです。彼らは優秀ですが、山脈が巨大になると、疲れ果てて時間が足りなくなったり、そこが底だと思い込んで局所的な窪みに捕まってしまったりします。
• 量子の方法（量子アニーリング）： 魔法の霧が、一瞬にして山脈全体の形を「感じ取る」ことができる様子を想像してください。量子アニーラーは、人間よりもはるかに速く、霧を通り抜けて（トンネル効果を利用して）最も低い地点を見つけ出すことができるハイカーなのです。

2. 戦略：パズルをバイナリ・ゲームに変える

量子アニーラーは、複雑な数式を直接理解することはできません。彼らは、0と1（スイッチのオンとオフのようなもの）という、より単純な言語を話します。

• 翻訳（QUBO）： 著者らは、複雑な核物理学の方程式を、「二次無制約バイナリ最適化（QUBO）」問題へと翻訳する必要がありました。これは、複雑なレシピを、「オン／オフ」のチェックリストへと変換する作業に似ています。量子マシンは、このスイッチの組み合わせをさまざまに試し、最も良い（エネルギーが最も低い）結果をもたらす組み合わせを見つけ出します。

3. 革新：玉ねぎの皮をむく（デフレーション）

最大の課題は、現在の量子アニーラーは、たった一つの解（絶対的な最低の谷）を見つけるのが最も得意であるという点です。しかし、科学者たちが本当に必要としているのは、最初のひとつだけでなく、解のリスト全体です。

• 解決策： 著者らは「ハイブリッド」な手法を作り上げました。
  1. ステップ1： 量子アニーラーを使用して、最初の解（最も低いエネルギー）を見つけます。
  2. ステップ2： 古典的コンピュータを使用して、「デフレーション（収縮）」を行います。これは、山脈の中で最も低い谷を見つけた後、ハイカーがそこに戻れないように、その谷を一時的にコンクリートで埋めてしまうようなものです。
  3. ステップ3： 「埋められた」地図を再び量子アニーラーに送り、次の低い地点を見つけさせます。
  4. 繰り返す： 全ての解のスペクトルが見つかるまで、玉ねぎの皮をむくように、一層ずつこの工程を繰り返します。

4. 結果：速度と精度

チームは、この手法を実際の量子コンピュータ（D-Wave Advantage）でテストし、標準的な古典的シミュレーション（シミュレーテッド・アニーリング）と比較しました。

• レース： 彼らは、異なるサイズのパズルを解くために、「量子ハイカー」と「古典的ハイカー」の間のレースを設定しました。
• 結果：
  • 小規模なパズルでは、どちらの手法も良好でした。
  • より大規模で複雑なパズルにおいては、量子ハイカーの方が大幅に高速でした。あるケースでは、古典的な手法が答えに近づくまでに数百ステップを要したのに対し、量子手法はわずか数十ステップで到達しました。
  • 量子手法は、より高いレベルの精度（正確さ）に、より迅速に到達しました。

5. 注意点：すべての道具がすべての仕事に適しているわけではない

論文では、谷を「埋める」ための3つの異なるデフレーション手法についてもテストを行いました。

• ホテリング（Hotelling）および直交射影（Orthogonal Projection）： これらはうまく機能しました。数学的な整合性を損なうことなく、量子マシンが次の解を見つけるのを効果的に助けました。
• ハウスホルダー（Householder）： この手法は単純なパズルには非常に有効でしたが、パズルが複雑になると（特に「一般化」固有値問題において）、破綻し始めました。それは、時計を修理するためにスレッジハンマー（大槌）を使おうとするようなものでした。全体像を捉えることには長けていますが、後のステップにおいて誤差を生じさせ、正確性を損なってしまうのです。

まとめ

この論文は、核物理学を永遠に解決したと主張しているわけではありません。その代わりに、近未来の量子コンピュータ（現在私たちが手にしている、ノイズが多く不完全なもの）が、有用なパートナーになり得ることを証明しています。探索を整理するための古典的コンピュータの信頼性と、最適な答えを見つけ出すための量子アニーリングのスピードを組み合わせることで、彼らは、これら特定の巨大な核のパズルに対して、古典的コンピュータ単独で使用する場合よりも高速で正確な手法を作り上げました。

これは、完璧でエラーのない量子コンピュータが登場する前であっても、実世界の物理学の問題に対して量子マシンを活用できる段階に近づいていることを示す、概念実証（プルーフ・オブ・コンセプト）なのです。

技術概要

技術要約：D-Wave 量子アニーラーにおける方程式運動フォノン法の固有値問題の解法

問題提起
原子核の多体問題には、核子の数とともに指数関数的に増大するハミルトニアン行列の固有値方程式を解くことが含まれる。これらはしばしば、厳密な対角化が計算上不可能となる次元に達する。古典的な反復法（例：ランチョス法、デビッドソン法）は疎性を利用するが、大規模な計算においてはストレージのボトルネックに直面する。量子コンピューティングは潜在的な代替案を提供し得るが、量子位相推定（QPE）アルゴリズムは、現在利用可能な誤り耐性を持つハードウェアを必要とする理論的に強力な手法である。したがって、近年のノイズあり中規模量子（NISQ）デバイスで大規模な固有値問題に対処できる、即時利用可能な戦略が必要とされている。具体的には、本研究では、核構造理論における方程式運動フォノン法（EMPM）から生じる標準的（SEVP）および一般化（GEVP）固有値問題の両方について、全固有スペクトル（基底状態だけでなく）を計算するという課題に取り組む。

手法
著者らは、量子アニーリングと古典的な行列デフレーション技術を統合したハイブリッド量子・古典アルゴリズムを提案している。コアとなる手法は、以下の3つのコンポーネントで構成される。

1. QUBO 定式化： レイリー商の連続最適化を離散化し、二次無制約バイナリ最適化（QUBO）問題へとマッピングする。固有ベクトルの成分を表す実変数は、特定の精度ベクトルを用いたバイナリベクトルによって近似される。これにより、結果として得られる二次形式のグローバルな最小値（または最大絶対値）を見つけることによって、固有値問題を量子アニーラー上で解くことが可能になる。
2. 反復降下と初期推測： アルゴリズムは2段階のアプローチを採用する。第一に、シフトされたQUBO問題を解くことで固有ベクトルの初期推測を生成する。第二に、目的関数の局所的な二次近似（勾配とヘッセ行列を含む）をQUBOとして符号化したものを最小化することで、この解を精緻化する反復降下フェーズを行う。離散化グリッドは、精度を高めるために段階的にスケールダウンされる。
3. 全スペクトルのための行列デフレーション： 主要な固有対だけでなく、全固有スペクトルを抽出するために、逐次的なデフレーション戦略を実装する。量子アニーラーが固有対を特定した後、古典的なデフレーション技術をハミルトニアン行列に適用して、見出された状態の影響を「剥ぎ取り」、次の反復が後続の固有状態をターゲットにできるようにする。評価された3つのデフレーション手法は以下の通りである：
   • ホテリング・デフレーション（Hotelling Deflation）： 減算的な手法（$A' = A - \lambda vv^T$）。
   • 直交射影（Orthogonal Projection）： 見出された固有ベクトルに対して直交する部分空間へ行列を投影する減算的手法。
   • ハウスホルダー・デフレーション（Householder Deflation）： ハウスホルダー反射を用いて行列を直交化する手法。

一般化固有値問題（GEVP）については、レイリー商はメトリック行列 $B$ を含むように修正され、デフレーションのステップは、可能な限り $B$-内積構造を保持するように適応される。

主な結果
本手法は、$^4$He 核の EMPM 計算から導出されたハミルトニアン行列を用いて、D-Wave Advantage System 6.4（5,614 量子ビット）および比較用の古典的シミュレーテッド・アニーリング（SA）を用いてベンチマークが行われた。

• 基底状態の性能： 量子アニーリング（QA）は、SA と比較して優れた収束性を示した。より大きな行列次元（例：$d=25$）および高いビット解像度（$b=4$）において、QA は約20〜30 回の反復でマシン精度（$10^{-8}$）を達成したが、SA は大幅に多くの反復（最大1,200回）を必要としたか、あるいは同等の精度に到達できなかった。行列サイズが増大するにつれて、性能の差は拡大した。
• 全スペクトルとデフレーション：
  • SEVP： 3 つのデフレーション手法（ホテリング、直交射影、ハウスホルダー）すべてが、高い精度（$b=2$ で最大13桁、$b=4$ で8桁）で全固有スペクトルを回収することに成功した。ハウスホルダー・デフレーションは、有効な部分行列のサイズを縮小する戦略により、計算実行時間の面で最も効率的であった。
  • GEVP： ホテリングおよび直交射影法は精度を維持したが、ハウスホルダー・デフレーションは GEVP に対して不安定であることが判明した。ハウスホルダー変換の繰り返しの適用は、メトリック行列 $B$ の対称性を歪め、丸め誤差を増幅させ、近接した固有値において数値的な劣化と収束の失敗を招いた。

意義と主張
本論文は、核構造計算において近未来の量子ハードウェアを利用するための実用的な戦略を提供すると主張している。主な意義は、ハイブリッド QUBO ベースのフレームワークを古典的なデフレーションと組み合わせることで、現在の変分量子アルゴリズムに共通する「基底状態のみの計算」という制限を超え、現実的な核ハミルトニアンの全固有スペクトルを系統的に再構築できることを示した点にある。

著者らは、固有値問題を既存のアニーラーで解可能な QUBO 定式化へと再構成することで、本アプローチが量子アルゴリズムとノイズのあるハードウェアの間の溝を埋めることに成功したと強調している。本研究は、現在の量子ビット数や結合性に制限されているため、扱う問題の規模については控えめなものであるが、スペクトル推定のための量子アニーリングの生存能力を確立している。本研究は、将来的な実装において、デフレーションのステップ自体を QUBO 問題として符号化することで、古典的なオーバーヘッドをさらに削減し、完全な量子アニーリングベースのスペクトルソルバーを目指すことを示唆している。