An approach for calculating astrophysical opacities on quantum computers

技術要約：量子コンピュータによる天体物理学的不透明度の計算手法

問題提起
天体物理学における不透明度（オパシティ）は、周波数の関数として媒体が光をどの程度吸収・散乱するかを定量化するものであり、恒星構造、慣性閉じ込め核融合爆縮、および高エネルギー密度プラズマのモデリングにおいて極めて重要である。古典的なモデルは通常、半古典的な光結合を用いた構成間相互作用展開に依存している。これらの手法は、原子番号が増加するにつれて関連する電子配置の数が増大するため、特に高温・高密度物質の領域において深刻なスケーリングの問題に直面する。さらに、古典的な取り扱いは多光子効果を無視することが多く、平衡電子構造や遷移の記述を劣化させる打ち切り（截断）の問題を抱えている。太陽の鉄の不透明度測定などの高精度な実験によって、理論予測との乖離が明らかになっているが、現在の古典的な計算手法は、これらの差異を解決するために必要な多体電子間相互作用および電子・光子結合を完全に捉えることに苦慮している。

手法
著者らは、原子的な系を通過する光子の波束を直接モデリングすることにより、天体物理学的な不透明度をシミュレートするための量子アルゴリズム・プロトコルを提案している。このアプローチは、ハイブリッドな第1量子化および第2量子化表現を利用する：

• 電子サブシステム： 電子は実空間格子内での第1量子化を用いて扱われ、これにより、極限的な熱力学的条件に関連する電子間相互作用および高度に励起された配置を正確に表現することが可能となる。
• 光子サブシステム： 光子は第2量子化を用いて扱われ、電磁場を波束として符号化する。これにより、多光子効果を自然に含めることができ、測定結果を不透明度へと直接マッピングすることが可能となる。
• ハミルトニアン： 系は、単一のプラズマイオン近傍における半相対論的な電子・光子結合を記述するパーリー・フィアッツ（Pauli-Fierz）ハミルトニアンの下で進化する。このモデルは、基礎的な出発点として、相対論的な電子効果、イオン間相互作用、およびイオンの運動を無視しているが、これらを後ほど組み込めるよう設計されている。
• シミュレーション・アルゴリズム： プロトコルは、相互作用描像におけるハミルトニアン・シミュレーションを採用している。この選択は、光子の自由場ダイナミクスが高速にフォワード可能なこと、および相互作用項（電子・光子結合）のノルムが電子および光子のハミルトニアンよりも著しく小さいことに基づいている。シミュレーションは以下の手順で行われる：
  1. 状態準備： 電子レジスタを（ギブンス回転を介して）熱的アンサンブルからサンプリングされたスレーター行列式へと初期化し、光子レジスタを最大占有のガウス型波束へと初期化する。
  2. 時間発展： 相互作用描像アルゴリズムを用いて、電子シミュレーション、光子シミュレーション、および相互作用項のブロック符号化を交互に繰り返しながら、複合系を時間発展させる。
  3. 測定： 光子レジスタに対して運動量分解測定を行い、透過した光子の数を決定する。
  4. 統計的平均化： 統計的な不確実性を解消し、熱的アンサンブルをサンプリングするために、回路を $N_s$ 回繰り返す。

主要な貢献

• プロトコル設計： 本論文は、量子ビット数と不透明度計算のスペクトル範囲および解像度の間の直接的な関係を確立している。不透明度の推定を、原子的な光子透過測定の直接的なシミュレーションとして定式化している。
• リソース見積もり： 著者らは、太陽の鉄の不透明度（放射対流境界付近の条件：$T \approx 2.26 \times 10^6$ K, $\rho_e \approx 4 \times 10^{22}$ cm$^{-3}$）を計算するための詳細な論理リソース見積もりを提供している。
  • これらの見積もりは、論理量子ビット数 ($Q$)、トフォリ（Toffoli）ゲート数、およびショット数 ($N_s$) を考慮している。
  • K殻電子については、これらが関連する遷移に参加しないため、グリッドサイズと計算コストを削減するための擬似化スキームが採用されている。
• スケーリング分析： 研究では、リソース要件が最大光子占有数 ($\Lambda$) に対してどのようにスケールするかを分析している。著者らは、ショット数を減らす効果と、1ショットあたりのトフォリ・カウントの増加とのバランスをとる最適な $\Lambda$ の範囲（$10^3$ から $10^4$ の間）を特定している。
• 古典的手法との比較： 太陽の鉄の問題に関するリソース見積もりは、最近の文献で分析された別の高エネルギー密度物理学の問題（慣性核融合ターゲット設計のための阻止能）と比較して、同等の規模であることが示されている。

結果

• リソース要件： 太陽の鉄の不透明度シミュレーションにおいて、総トフォリ・カウントは、$\Lambda \ge 10^1$ の場合（悲観的な境界の下）にはハミルトニアン・シミュレーションのコストによって支配され、低い $\Lambda$ の場合（楽観的な境界の下）には状態準備のコストによって支配される。総コストは、同じ時間における純粋な電子構造ハミルトニアンのシミュレーションコストよりも約2桁大きく、これは明示的な光子場によって導入されるオーバーヘッドを浮き彫りにしている。
• スペクトル解像度： プロトコルは、7～13 Åの波長範囲において、$\sim 0.1$ Åの幅を持つスペクトル特徴を分解することができる。
• ショット数： 目標とする精度を達成するために必要なショット数は、吸収確率と1ショットあたりの光子数に依存する。著者らは、スペクトル線形状に関する事前情報を用いることで、特徴を意味のあるレベルで解明するために必要な1つの特徴あたりのショット数 ($N_{pf}$) を、約10まで削減できることを示している。

意義および主張
本論文は、将来の量子コンピュータを用いて天体物理学的な不透明度を予測するための基礎的なプロトコルを提供すると主張している。これは、既存の理論と実験の間の乖離を即座に解決することを目的としたものではなく、古典的なコンピュータでは不可能な多体効果および多光子効果を捉えるための経路を提示するものである。著者らは、この研究を、電子間相互作用を正確に保持しつつ他の物理的側面を簡略化した「基礎的な出発点」と位置付けている。彼らは、現在の見積もりが他の困難な高エネルギー密度物理学の問題と同等であることを指摘しており、本アプローチは高エネルギー物理学におけるシミュレーション手法と概念的に類似しており、他の光・物質相互作用（例：分子スペクトル）への転用が可能であるとしている。本論文は、提案されたプロトコルが、古典的な手法よりも効率的に、電磁場と結合したプラズマのヒルベルト空間を表現する量子コンピュータ独自の能力を活用しており、高温・高密度物質の領域において量子優位性を提供する可能性があることを強調している。