Thermodynamic sampling of materials using neutral-atom quantum computers

技術要約：中性原子量子コンピュータを用いた材料の熱力学的サンプリング

問題提起
技術的に重要な材料は、しばしば構成的無秩序（例：固溶体、空孔、置換合金）を示します。それらの特性を予測するには、正確なエネルギーモデルだけでなく、特定の温度および化学ポテンシャルにおける平衡状態を決定するために、膨大な原子配置の空間をサンプリングできる手法が必要です。モンテカルロ・サンプリングやシミュレーテッド・アニーリングングといった古典的な手法は、システムサイズが増大するにつれて計算上のボトルネックに直面し、網羅的な列挙が困難になります。量子アニーリングはこの目的のために探索されてきましたが、現実的な材料モデル（特に密度汎関数理論（DFT）から導出されたもの）を現在の量子ハードウェア上にマッピングすることには、大きな課題が存在します。これらには、相互作用範囲、エネルギー尺度に関するハードウェアの制約、および中性原子アレイにおける周期境界条件のネイティブなサポートの欠如が含まれます。さらに、量子オプティマイザの出力分布に対する厳密な熱力学的解釈を確立することは、多くのデバイスが平衡熱アンサンブルよりもコヒーレントなユニタリダイナミクスとして記述されることが多いため、未解決の課題として残っています。

手法
著者らは、窒素ドープされたグラフェンの熱力学的特性を、中性原子量子ハードウェア（具体的にはQuEra Aquilaデバイス）に適したリドバーグ原子ハミルトニアンにマッピングするフレームワークを提案しています。その手法は以下のステップで進行します。

1. DFTデータの生成: 78原子のスーパーセルを用いた周期境界DFT（CRYSTAL23コード）を使用して、窒素ドープされたグラフェンの形成エネルギーを計算します。対称独立構成（SIC）を生成して緩和し、トレーニングおよびテストセットを拡張するために対称等価構成（SEC）を導出します。
2. ハミルトニアンのマッピング: DFTの形成エネルギーをリドバーグ・ハミルトニアンに適合させます。これには、オンサイト項（$V$）と距離依存のペア相互作用（$C_6/R^6$）からなるモデルと、DFTエネルギーとの間の誤差を最小化するプロセスが含まれます。このフィッティングにより、オンサイトエネルギー $V^{DFT}$ と有効な最近接距離 $R^{DFT}_{NN}$ が得られます。
3. ハードウェア・リスケーリング戦略: DFTから導出されたパラメータとハードウェアの能力との間には、決定的な不一致が存在します。ハードウェアは最小原子間隔（$R^{min}_{NN} = 4\,\mu\text{m}$）を必要としますが、フィッティングされた $R^{DFT}_{NN}$ はそれよりも大幅に小さい（$\approx 1.6\,\mu\text{m}$）値となります。その結果、ハードウェア上のネイティブな二体相互作用の強さは、対象となる材料の相互作用よりも数桁弱くなります。
   • これを解決するために、著者らはリスケーリング係数 $\alpha_v = (R^{min}_{NN} / R^{DFT}_{NN})^6 \approx 236.69$ を導入します。
   • ハードウェアのアクセス可能なエネルギーウィンドウ内に収めるため、オンサイトエネルギーとペア相互作用の両方を $\alpha_v$ でスケーリングします。
   • このリスケーリングにより、ハードウェアのサンプリング温度 $T$ と、有効な材料温度 $T' = \alpha_v T$ の間の対応関係が確立されます。
4. 化学ポテンシャルの制御: 量子アニーラにおけるグローバルなレーザーデチューニング（$\Delta_g$）は、材料のグランドカノニカル化学ポテンシャルの差（$\Delta\mu$）にマッピングされます。このリスケーリング係数により、ハードウェアの限られたデチューニング範囲が、材料モデルにおける化学的に関連のある $\Delta\mu$ の範囲をカバーできるようになります。
5. 検証:
   • 28サイト・システム: 全ての可能な構成の網羅的な古典的列挙に対して、マッピングの妥当性を検証します。
   • 78サイト・システム: 網羅的な列挙が不可能なより大きなシステムについては、ランダムなモンテカルロ（MC）サンプリングと比較してベンチマークを行います。
   • 温度チューニング: ハードウェア内の最近接間隔（$R_{NN}$）を変化させることで、有効温度を実験的に調整し、サンプリング分布の制御を実証します。

主な貢献

• リスケーリング・フレームワーク: 本論文は、DFT由来のエネルギー尺度と現在の中性原子ハードウェアの限定的なダイナミックレンジとの間のギャップを埋める、体系的なリスケーリング手順（$\alpha_v$）を導入しています。これにより、ハードウェアは有効温度 $T'$ におけるターゲット材料に対応する分布をサンプリングすることが可能になります。
• 熱力学的解釈: ハードウェアの制御パラメータ（レーザーデチューニングおよび原子間隔）を熱力学変数（化学ポテンシャルおよび温度）に結びつける厳密な導出を提供し、量子アニーリングの出力をボルツマン的なアンサンブルとして解釈することを可能にします。
• ハードウェア制約への対処: 幾何学的制約（最小原子間隔、2D配置）および中性原子デバイスにおける周期境界条件の欠如に明示的に対処し、有限のナノフレークをバルク挙動の代表として扱います。
• プルーフ・オブ・コンセプト: プログラマブルな量子アニーラを用いて、ドープされた材料システムの熱力学的特性（平均組成およびエネルギー分布）の抽出に成功したことを実証しました。

結果

• 28サイト・ベンチマーク: 28サイトのグラフェン・ナノフレークについて、量子アニーラの出力（平均窒素濃度対化学ポテンシャル）は、網羅的な列挙による正確な結果と密接に一致しました。最良の適合は、有効サンプリング温度 $T = 41\,\mu\text{K}$（これにより材料の有効温度 $T' \approx 9.7\,\text{mK}$ が得られる）に対応します。量子分布と古典的分布の間の全変動距離（TVD）は0.294であり、中程度の一致を示しています。この偏差は、低エネルギー状態の温度に対する敏感さに起因すると考えられます。
• 78サイト・ベンチマーク: 78サイトのシステムでは、量子アニーラは低エネルギーの構成を優先的にサンプリングしますが、ランダムなモンテカルロ・サンプリングはエネルギーランドスケープ全体を一様に探索します。量子デバイスは、組み合わせ論的に稀であり、ランダムサンプリングでは捉えるのが困難な、エネルギーがゼロの構成（相互作用のない窒素配置）を正常に特定しました。
• 温度制御: 原子間隔を $4\,\mu\text{m}$ から $5\,\mu\text{m}$ に増加させることで、有効サンプリング温度を $\approx 9.7\,\text{mK}$ から $\approx 36\,\text{mK}$ へと調整できることを実証し、リスケーリング戦略の熱力学的整合性を確認しました。

意義および主張
著者らは、本研究を、現実的な材料モデル（DFTから導出されたもの）をリドバーグ・ハミルトニアンにマッピングするための**プルーフ・オブ・コンセプト（概念実証）**として位置付けています。

• 主な意義は、深刻なエネルギー尺度の不一致にもかかわらず、利用可能な量子ハードウェア上で現実的な材料モデルの実装を可能にするワークフローを確立したことにあります。
• 本論文は、提案されたリスケーリング戦略により、ハードウェアが熱力学的な意味（具体的には有効温度におけるボルツマン的なアンサンブル）で解釈可能な分布をサンプリングできることを主張しています。
• 著者らは現在の限界についても認めており、本手法は二体相互作用および2Dシステムに限定されていることを述べています。彼らは、強い局在相互作用や多体効果（例：リン・窒素共ドープ）を持つシステムは、現在の二体ハミルトニアンでは正確に捉えることができず、高次の結合またはより柔軟な相互作用モデルが必要であることを明示しています。
• 本研究は、研究された特定のシステムにおいて古典的手法を凌駕することを主張するものではなく、むしろ将来の材料探索のためのハイブリッド・ワークフローの構成要素として、中性原子量子コンピュータを使用することの実現可能性を示すものです。