Quantum speed-up for solving the one-dimensional Hubbard model using quantum annealing

本論文は、1 次元ハバードモデルに対するゲート型量子アニーリングシミュレーションが、最大 40 量子ビットの半充填系における基底状態の探索において、古典的なベテアンサッツアルゴリズムに対して実質的な量子加速を達成することを示す。

要約

広大で霧に包まれた山脈の中で、絶対的な最低点を見つけようとしていると想像してください。この「最低点」は、物質中を移動する電子（電気を運ぶ微小な粒子）の系が持つ、最も安定した静かな状態を表しています。物理学において、この特定の山脈はハバードモデルと呼ばれます。何十年もの間、科学者たちは複雑な数学を用いてこれらの山脈を地図化してきましたが、山脈が大きくなる（電子の数が増える）につれて、数学はあまりにも重くなり、世界の最速のスーパーコンピュータでさえ、膨大な時間をかけずに底を見つけるのに苦労しています。

この論文は、単純な問いを投げかけます：量子コンピュータは、従来の数学よりもこの最低点を速く見つけることができるでしょうか？

以下に、日常の比喩を用いて著者たちがどのようにこの問題に取り組んだかを説明します。

1. 問題点：「ベテ・アンサッツ」の山脈

この電子問題の一次元版（電子の単一の列）については、科学者たちはすでにベテ・アンサッツ方程式と呼ばれる地図を持っています。

• 従来の方法： これは、複雑な結び目で絡み合った巨大なジグソーパズルのピースを解こうとするようなものです。解くことは可能ですが、パズルが大きくなるにつれて、結び目を解くのに要する時間は非常に急速に増大します。論文は、エネルギーは比較的速く計算できるものの、実際にはすべての個々の電子の特定の配置（「基底状態」）を特定するには、指数関数的な数の詳細を計算する必要があると指摘しています。それは、潮が最も引く正確な場所を見つけるために、砂浜のすべての砂粒を数えようとするようなものです。

2. 解決策：量子アニーリング（「氷を溶かす」方法）

ピースごとにパズルを解く代わりに、著者たちは量子アニーリングと呼ばれる手法を用いました。

• 比喩： 中に隠れた物体が凍りついた氷のブロックを持っていると想像してください。その物体を壊さずに取り出したいとします。
  • ステップ 1： 物体を見つけやすい単純で平坦な氷のブロック（「初期ハミルトニアン」）から始めます。
  • ステップ 2： 氷をゆっくりと溶かし、その形を徐々に変化させて、関心のある複雑でギザギザした山脈（「ハバード・ハミルトニアン」）と全く同じになるようにします。
  • ルール： 氷を十分にゆっくり溶かせば、内部の物体は形が変わるにつれて自然に可能な限り低い点へと滑り落ちます。小さな障壁をすり抜けられる「量子」的な性質があるため、高いピークに留まることは決してありません。

3. 実験：溶かす過程のシミュレーション

彼らの研究所には巨大な量子コンピュータが存在しなかったため、量子コンピュータがどのように振る舞うかをシミュレートするために、強力な古典的なスーパーコンピュータを使用しました。

• 彼らは、溶かす過程を模倣するデジタル「回路」（一連の指示）を構築しました。
• 彼らは最大40 量子ビット（ビットに相当する量子の単位）を持つシステムでこれをテストしました。これを理解しやすく言えば、40 量子ビットをシミュレートすることは、小さな部屋にあるすべての粒子の位置を同時に追跡しようとするようなもので、通常のコンピュータにとっては信じられないほど困難なタスクです。
• 彼らは、底を見つけるのにどれくらい時間がかかるかを見るために、異なる「溶かす速度」（アニーリング時間）でシミュレーションを実行しました。

4. 結果：速度向上

この論文は驚くべき結果を見つけました。

• 従来の数学： システムが大きくなるにつれて、従来の方程式を用いて基底状態を見つけるのに必要な時間は爆発的に（指数関数的に）増大します。それは、電子を一つ追加するたびに山脈が突然二倍の高さになるようなものです。
• 量子的手法： 量子アニーリング手法で基底状態を見つけるのに必要な時間は、線形的（あるいはそれ以下）に増大しました。これは、システムのサイズを倍にすれば、答えを見つけるために必要な時間も倍（あるいはわずかに増加）するだけで済むことを意味します。
• 結論： 電子の半分の列という特定のケースにおいて、量子手法は大幅な速度向上を提供します。それは、一歩進むごとに高さが二倍になる山を登るのと、わずかに高くなるだけ丘を登るのとの違いです。

5. なぜこれが重要なのか（論文によると）

著者たちは、これは「玩具的な問題」（単純化されたモデル）であると強調していますが、それは重要な点を証明しています。

• すでに数学によって「解かれている」系（可積分系）であっても、量子コンピュータは解を見つける方法において巨大な利点を提供する可能性があります。
• このスケーリングが真実であれば、量子アニーリングは、電子の実際の状態を見つけるための最良の古典的手法と比較して、これらの問題を指数関数的な速度向上で解決できることが示唆されています。
• また、彼らはこれが機能する理由は、彼らが登っている「山脈」（1 次元ハバードモデル）が、系を閉じ込めるような急峻で危険な崖（相転移）を持っていないからだと指摘しています。

要約：
この論文は、シミュレートされたコンピュータ上で量子の「溶かす」手法（アニーリング）を使用することで、従来の数学が許すよりもはるかに速く電子の最も安定した状態を見つけられることを実証しています。この特定のモデルは単純化された電子の列ですが、これは概念実証であり、量子コンピュータが最終的に、現在の最速のスーパーコンピュータでは処理しすぎてしまう複雑な材料科学の問題を解決できることを示しています。

技術概要

技術的概要：量子アニーリングを用いた一次元ハバードモデルの解法における量子加速

問題定義
一次元ハバードモデルは、凝縮系物理学において相関電子系を理解するための基本的な枠組みである。このモデルは「可積分」であり、すなわちベテ・アンサッツ方程式を通じて熱力学的極限における基底状態の厳密な解析的解が存在するが、開放境界条件を持つ有限系对于这些方程式を解くことは、計算上の重大な課題を呈する。具体的には、基底状態のエネルギーは連立非線形ベテ・アンサッツ方程式を解くことで多項式時間で求められる一方、完全な基底状態波動関数（係数）を決定するには、指数関数的に大量の項の計算が必要となり、指数関数的な計算資源を要する。著者らは、有限の半充填一次元ハバード系の基底状態を求める際、量子アルゴリズム、特に量子アニーリングがこれらの古典的手法に対して加速を提供しうるかを調査する。

手法
著者らは、汎用量子コンピュータシミュレータ（JUQCS）上でゲートベースの量子アニーリングシミュレーションを実装し、ハバードハミルトニアンを解いた。手法は以下の 3 つの主要段階からなる：

1. ハミルトニアンのマッピング: フェルミオンハバードハミルトニアンを、ジョルダン・ウィグナー変換を用いて量子ビットハミルトニアンにマッピングする。これにより、フェルミオンの生成・消滅演算子が $2L$ 個の量子ビット（$L$ は格子サイトの数）に作用するパウリ行列に変換され、反交換関係が保存される。
2. 初期状態の準備: アニーリング過程は、非相互作用のホッピングハミルトニアン（$H_I$）の基底状態から開始される。著者らはギベンズ回転に基づくアルゴリズムを用いてこの初期状態を準備する。これは、与えられたスピンアップおよびスピンダウン電子の数に対して、最低エネルギーの運動量状態を充填することに相当する。
3. アニーリングダイナミクス: 系は時間依存ハミルトニアン $H(s) = (1-s)H_I + sH_H$ の下で進化し、ここで $s$ は 0 から 1 まで変化する。時間発展は、連続的な進化を離散時間ステップに分割するために、スズキ・トロター分解（具体的には 2 次近似）を用いてシミュレートされる。各ステップのユニタリ演算子は、汎用的な 1 量子ビットおよび 2 量子ビットゲート（アダマード、$R_Z$、$R_{ZZ}$、および $X$ ゲート）に分解される。
4. スケーリング解析: シミュレーションは $L=20$（40 量子ビット）までの系サイズで行われた。性能指標は、残留エネルギー $\Delta E = \langle \psi(1) | H_H | \psi(1) \rangle - E_0$ であり、ここで $E_0$ はベテ・アンサッツを通じて計算された厳密な基底状態エネルギーである。著者らは、特定の精度を達成するために必要な総アニーリング時間（$T_A$）が系サイズ（$L$）とどのようにスケーリングするかを分析する。

主要な貢献と結果

• ゲートベース実装: 本論文は、ゲートベース量子コンピュータ上でハバードモデルに対する量子アニーリングをシミュレートするための具体的な構成を提供し、回路の分解と資源要件（ゲート数）を詳述している。
• アニーリング時間のスケーリング: 半充填系（$2N_\downarrow = N = L$）において、著者らは線形アニーリングスケジュールの場合、残留エネルギーが $\Delta E \propto L^{1.234} / T_A^2$ としてスケーリングすることを発見した。特定の精度 $\Delta E^*$ を達成するために必要なアニーリング時間は、ある閾値 $L^*$ 以下の系サイズにおいて $T_A \propto L^{0.62}$ としてスケーリングする。
• 強結合領域: 強い相互作用強度（$U/t_H = 8, 16$）に対するシミュレーションは、スケーリング領域の開始が穏やかであることを確認し、任意に大きな $L$ に対してアニーリング時間が系サイズの線形関数（$T_A \propto L^{0.995}$）によって有界であることを示した。
• 古典的手法との比較: 著者らは、ベテ・アンサッツ方程式に対する古典的な根発見アルゴリズムとの比較を行った。エネルギーをベテ・アンサッツで求めることは多項式時間（$T_{BA} \propto L^{3.3}$）であるが、完全な基底状態係数を取得することは指数関数的に高価である。これに対し、量子アニーリングアプローチは、系サイズに対して時間および資源のサブリニアからリニアなスケーリングで基底状態を導き出す。

意義と主張
本論文は、一次元ハバードモデルの基底状態を求める際、量子アニーリングが古典アルゴリズムに対して実質的な量子加速を提供すると主張している。

• 状態準備における指数関数的加速: 著者らは、エネルギーが多項式時間で発見されたとしても、古典的手法が基底状態係数を決定するために指数関数的な資源を必要とするため、量子アニーリング時間の多項式（具体的にはサブリニアからリニア）スケーリングは、基底状態の準備という文脈において指数関数的な量子加速を表すと論じている。
• 断熱計算の妥当性検証: 結果は、可積分な量子多体問題に対する断熱量子計算の有効性を支持する。観測された多項式スケーリングは、非ゼロの相互作用強度において一次元ハバードモデルに相転移が存在しないことに起因しており、これによりエネルギーギャップが指数関数的に閉じることを防いでいる。
• 将来展望: 本研究は 1 次元可積分ケースに焦点を当てているが、著者らはこれらの結果が非可積分系（例えば 2 次元ハバードモデル）へのさらなる調査を促すと示唆している。高次元における相転移がギャップの挙動を変化させる可能性があることに留意しつつ、1 次元ケースで実証された優位性は、複雑な多体問題を解くための量子アニーリングの可能性を強化するものである。

著者らは、即時の応用については慎重であり、現在の数千ゲート規模の計算は近未来のハードウェア能力を超えていることを認めつつも、資源のスケーリングは将来のフォールトトレラントアーキテクチャにとって有望であると述べている。