Green's Functions from Sample-based Krylov Quantum Diagonalization: An Impurity Solver for Dynamical Mean-Field Theory

超伝導体のように電気抵抗なしで電気を運ぶものや、ハードドライブの動力源となる磁性材料のように、奇妙で強情な振る舞いを見せる物質は、通常の固体を説明するために用いる単純な法則をしばしば無視します。これらの材料では、電子は独立して動いているわけではありません。電子同士が互いに押し引きし合う力が非常に強いため、一つの電子の振る舞いは隣接する電子の振る舞いと不可分に結びついています。これらの材料がどのように機能するかを理解するために、科学者たちは「グリーン関数」という、混雑した混沌とした環境の中を電子がどのように移動するかを示す地図のような役割を果たす数学的ツールを計算しなければなりません。この地図は、材料のエネルギー準位や、光や電気に対する反応を明らかにしますが、大規模なシステムに対してこれを計算することは極めて困難な作業です。電子の可能な配置の数は、システムのサイズとともに爆発的に増加するため、最も強力なスーパーコンピュータであってもすぐにメモリ不足に陥り、研究者は実際に注目している材料の、非常に小さく簡略化されたバージョンしか研究できない状況に置かれます。

研究チームは、システム全体を一度にシミュレートしようとするのではなく、巧妙で小さな「スナップショット」を連続的に撮ることで、量子コンピュータを用いてこの問題に取り組む新しい方法を提案しました。金属中の電子の振る舞いを模したモデルシステムに焦点を当てた研究において、彼らは「サンプルベースのクリロフ量子対角化（sample-based Krylov quantum diagonalization）」と呼ばれる手法を応用しました。量子コンピュータに対して、システムの巨大で膨大な状態をすべてメモリに保持するように求めるのではなく、システムを極めて短い時間だけ進化させ、その結果を測定するように求めたのです。このプロセスを何度も繰り返し、最も一般的な結果を収集することで、彼らは、従来の方法で必要とされる不可能な量のデータを保存することなく、電子の動きの完全な地図を再構成することができました。このアプローチは、古い手法で必要とされるような長く複雑な回路を実行することがまだできない、現在利用可能なノイズが多く不完全な量子コンピュータでも動作するように設計されています。

研究者たちは、この手法を、単一の原子が電子の海と相互作用する様子を描写する「単一不純物アンダーソン・モデル」として知られる特定のモデルを用いてテストしました。彼らは、中央の原子が周囲の環境を表す12サイトの鎖に接続されているシミュレーションを設定しました。この設定では、中央の原子が強い反発力を持ち、二つの電子が同じ場所に存在することを防いでいます。これは、金属（電気を通すもの）と絶縁体（電気を通さないもの）の間の転移を駆動する重要な特徴です。チームは、量子デバイスをシミュレートするために古典的なコンピュータを使用し、これにより、新しい手法をこのシステムの正確で完璧な解と比較することができました。その結果、量子進化をサンプリングすることで、全状態のほんの一部を利用するだけで、電子のエネルギーマップを高精度に再構成できることがわかりました。12サイトのシステムに対して、この手法は必要な計算空間を3桁から4桁減少させました。つまり、完全な計算に必要とされる数億もの状態の代わりに、わずか数千の状態を扱うだけで済んだのです。

この研究により、この手法は、システムが完全に単純でも完全に混沌としてもいない状態にあるときに最も効果的に機能することが明らかになりました。電子間の相互作用が非常に弱いとき、システムは独立した粒子の穏やかな海のように振る舞い、相互作用が非常に強いとき、電子は硬直したパターンにロックされます。これら両方の極端なケースにおいて、手法は容易に答えを見つけ出しました。しかし、中間領域、つまり材料が金属から絶縁体へと変化する境界にあるとき、電子は複雑で変動する状態にあります。ここで研究者たちは、手法が完全な全体像を捉えるためにより多くのサンプルを必要とすることを発見しましたが、それでも従来のメソッドが失敗するところで成功を収めました。結果は、材料の特性を定義するエネルギーマップにおける独特なピークや谷である「スペクトル特徴」が、幅広い相互作用の強さにわたって忠実に再現されていることを示しました。

最も重要な発見の一つは、この手法が、現在のデバイスにおけるエラーの要因となりやすい「アンシラ（ancilla）」と呼ばれる補助ビットを用いた複雑な操作を量子コンピュータに要求しないことでした。代わりに、量子デバイスは状態を準備し、短時間進化させ、そして電子の位置を測定するだけで済みます。これらの測定値を最終的な答えへと結びつけるための重労働はすべて、古典的なコンピュータによって行われました。この役割分担により、量子回路は浅く単純なものに保たれ、近い将来利用可能となるハードウェアにとって現実的な候補となります。研究者たちは、この手法は現在はシミュレーションであるものの、システムが大きくなるにつれて必要なサンプル数の増え方が有望であることを指摘しました。鎖のサイト数が8から12へと増加するにつれ、正確な答えを得るために必要な全状態の割合は実際に減少しており、これは手法が現在手の届かないさらに大きなシステムにもスケールできる可能性を示唆しています。

この研究はまた、これらの計算においてどの出発点が最適かという一般的な問いにも対処しました。チームは、システムの初期状態を設定するための2つの異なる方法をテストしました。一つは単純な非相互作用に基づくものであり、もう一つはより複雑な自己整合的な図に基づくものです。彼らは、この特定のタスクにおいては、より複雑な出発点と同じくらい単純な出発点が効果的であることを発見しました。これは、将来の実験におけるセットアップを簡素化します。さらに、彼らは、電子がシステムに追加される、あるいは取り除かれる様子を計算するために必要な数学的な「分岐（branching）」を扱うことができることを実証しました。サンプリングされた状態のうち、どの状態を保持し、どの状態を破棄するかを慎重に選択することで、最終的な結果が確率の保存といった物理学の基本原則を遵守するようにしました。

このアプローチは、強相関材料を研究するために用いられる強力な枠組みである「動力学的平均場理論（dynamical mean-field theory）」への有望な道筋を提供します。この枠組みでは、複雑な原子格子が単一の不純物問題へとマッピングされ、材料全体の性質を見出すために、この問題を繰り返し解かなければなりません。ボトルネックは常に「不純物ソルバー（impurity solver）」、つまり単一のサイトにおける電子の振る舞いを計算するステップでした。もしこのソルバーが、現在よりも多くのバスサイト（bath sites）を持つ量子コンピュータ上で実行できるようになれば、科学者はより高い精度で材料をモデル化できるようになります。研究者たちは、彼らの手法が新しい不純物ソルバーとして機能し、これまでよりもはるかに大きく、より現実的なシステムの研究を可能にすると示唆しています。この研究は特定のモデルに限定されており、まだ理論の完全な自己整合ループを含んでいませんが、その結果は、このような計算が実現可能であるという強力な概念実証を提供しています。

また、研究では、電子の反発力の強さを変化させるなど、異なる条件下での手法の性能についても調査しました。特定の精度レベルに達するために必要なサンプル数は、単純な直線関係にはならず、中間的な相互作用の強さでピークに達した後、再び低下することがわかりました。この挙動は物理的に理にかなっています。なぜなら、中間領域はシステムが最も複雑であり、記述するために最も多くの情報を必要とする場所だからです。しかし、このピークがあるにもかかわらず、必要な総サンプル数は、最も困難なケースにおいても全状態のわずかな割合にとどまりました。研究者たちは、この手法が堅牢であり、シミュレーションで見られた良好なスケーリングは、近い将来の量子ハードウェアにとって実用的なツールになり得ることを結論付けました。

結局のところ、この研究は、最も困難な問題の一部を解決するために、完璧でエラーのない量子コンピュータを待つ必要はないということを示しています。確率的な量子測定を受け入れ、巧妙なサンプリング技術を用いることで、研究者は複雑な材料の中で電子がどのように動くかについての不可欠な情報を抽出することができます。この手法は、以前は手に負えないと考えられていた問題を、管理可能なものへと変貌させ、高度な材料の設計や理解を助ける次世代のシミュレーションへの扉を開きます。今回の知見は、適切なアルゴリズムがあれば、現在のノイズの多い量子コンピュータであっても、強相関物質の最も深い謎に光を当て始めることができることを示唆しています。