Ground State Energy via Adiabatic Evolution and Phase Measurement for a Molecular Hamiltonian on an Ion-Trap Quantum Computer

技術要約：イオントラップ量子コンピュータにおける、断熱進化と位相測定を用いた分子ハミルトニアンの基底状態エネルギー推定

問題提起
分子の基底状態エネルギー（GSE）を正確に推定することは量子計算における主要な応用分野であるが、正確な量子状態の準備、高精度な観測量の推定、およびハードウェアノイズの緩和が困難であるため、実用化への道のりは険しい。古典的な手法であるハートリー・フォック（HF）法や結合クラスター法は、強相関系における本質的な電子相関効果を捉えることができず、より高度な古典的手法は指数関数的なスケールアップに直面する。量子アルゴリズム、特に断熱状態準備（ASP）と量子位相推定（QPE）の組み合わせは系統的な改善を提供するが、パラメータの最適化、バレン・プラトー（平坦な勾配）、および制御された$U^{2^k}$演算による高い回路深さといった課題に直面している。さらに、ソフトウェアによって緩和可能なエラーと、ハードウェアレベルでの削減を必要とするエラーを区別することは、この分野を進展させる上で極めて重要である。

手法
著者らは、Quantinuum社のH1-1捕捉イオン量子コンピュータを用いて、$H_3^+$分子（6個の量子ビットにエンコード）のGSEを推定するための、純粋に量子的なab initioプロトコルを実装した。ワークフローは、主に以下の3つのコンポーネントで構成される：

1. 断熱状態準備 (ASP): 離散化誤差を導入する標準的なトロッター・鈴木分解の代わりに、著者らはTETRIS（Time Evolution Through Random Independent Sampling）アルゴリズムを採用した。この手法は、平均的に時間発展演算子を離散化誤差なしに正確に再現するランダム化ユニタリを生成する。システムは、単純な初期ハミルトニアン（$H_Z$）の基底状態から、ターゲットとなる分子ハミルトニアン（$H$）へと進化する。効率を向上させるため、著者らは保存量項（粒子数演算子）を付加することでハミルトニアンの相互作用ノルム（$\mu_I$）を最適化し、回路の深さを短縮するためにパリティベースおよび占有数ベースの演算子削減を含む古典的な前処理技術を適用した。
2. エネルギー推定: 著者らは、ノイズに強い反復量子位相推定（iQPE）の変種を利用した。期待値を直接サンプリングするのではなく、アダマールテストを用いて、重なり $\rho(s) = \text{Im} \langle \psi(T) | e^{is(E_{\text{trial}} - H)} | \psi(T) \rangle$ の虚部を測定する。2つの試行エネルギー（$E_{\text{trial}}^\pm$）を選択し、固定された中心時刻 $s$ において $\rho(s)$ を測定することで、逆正接（arctangent）公式を介してGSEを抽出する。このアプローチは振幅情報を活用しており、デポラリゼーション・ノイズに対して堅牢である。
3. エラー緩和と後処理: 実験では、2つの古典的な後選択技術を用いた：
   • パリティ後選択: 物理量子ビットの総パリティが期待されるセクターから逸脱した場合、それはリークまたはエラーを示唆するため、その測定結果を破棄する。
   • ハートリー・フォック（HF）投影: 非アンシラ量子ビットが初期のHF状態に戻らなかったショットに対して、アンシラの期待値をゼロに設定する。これは、ノイズレスな設定ではそのような寄与が平均的にゼロになるという理論的期待に基づいている。

主な結果

• エネルギー推定: 実験により、厳密なGSEからの偏差が $-25.5 \pm 5.3$ mHaである最終的なエネルギー推定値が得られた。これは、古典的なハートリー・フォック推定の偏差である $52.8$ mHaと比較して改善されている。量子的な結果は、5標準偏差を超えてHFエネルギーから統計的に区別可能である。
• ノイズ分析: 著者らは、異なるノイズ源の影響を分離するために広範な古典シミュレーションを実施した：
  • コヒーレントおよびインコヒーレント・ノイズ: 現実的なコヒーレントエラー（例：$Z$回転）およびインコヒーレントエラー（デポラリゼーション）のレートを組み込んだシミュレーションでは、化学的精度（$\sim 2$ mHa以内）に近いエネルギー偏差が予測された。これは、アルゴリズムがこれらの標準的なノイズに対して耐性があることを示唆している。
  • リークエラー: シミュレーションにおいてリークエラー（計算部分空間からの遷移）を導入すると、顕著な系統的バイアスが生じた。シミュレーションにおけるバイアスの大きさは、リーク率を含めた場合にのみ、実験的な偏差（$\approx -25.5$ mHa）と一致した。
• リークの兆候: 測定統計の分析により、ハードウェアのデータにおいて、ノイズレスなシミュレーションと比較して、物理量子ビットおよびアンシラ量子ビットの両方で $|1\rangle$ 状態を測定する強いバイアスがあることが明らかになった。この不均衡は、リーク事象、特に不釣り合いなシャトリングやゲート操作を受けるアンシラ量子ビットがリークに対して非常に脆弱であることを示唆している。

意義と主張
本論文は、非スケーラブルなタスクを古典コンピュータにオフロードすることに依存しない、実用的かつスケーラブルな基底状態エネルギー推定のための量子アルゴリズムを実証している。主要な貢献は、現在のイオントラップ・ハードウェアにおいて、化学的精度を達成するための主要な障害はコヒーレントまたはインコヒーレント・ノイズではなく、リークエラーであることを特定した点にある。

著者らは、リークエラーがなければ、彼らの実験設定とノイズに強いiQPEプロトコルは、ショットノイズを含めても化学的精度（2倍の範囲内）に近い結果を達成できていたはずであると主張している。したがって、本研究は、将来のイオントラップ・デバイスにおける量子化学アプリケーションの進展は、標準的なゲート忠実度の低減だけでなく、リーク抑制（例：アンシラ量子ビットのためのリーク検出コード）に向けたハードウェア設計の改善や、アルゴリズムの適応をターゲットにすることに決定的に依存していることを強調している。