Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

技術要約：周期境界条件およびディリクレ境界条件下における相対論的量子シミュレーション

問題提起
量子系のシミュレーションは、システムサイズとともに複雑性が指数関数的に増大するため、計算負荷が高い。これは、物理現象を正確に記述するために相対論的効果を組み込む必要がある場合に、さらに深刻化する課題となる。摂動論は、非相対論的な枠組みを超えて相対論的効果を研究するための標準的な手法であるが、近未来の量子デバイス上でこれらのシミュレーションを実装することは依然として困難である。既存の手法は、量子化学において標準的な第二量子化に依存することが多いが、非線形性を伴うシステムの直感的な取り扱いには欠ける場合がある。したがって、異なる境界条件下にあるシステムの相対論的な基底状態エネルギーを、近未来のハードウェア上で推定できる、実用的な第一量子化フレームワークが必要とされている。

手法
著者らは、システム量子ビットによって表される有限の格子上に波動関数を離散化する、相対論的量子シミュレーションのための第一量子化フレームワークを提案している。コアとなる手法は、以下のコンポーネントで構成される。

1. 離散化と並進演算子:
   連続的な1次元空間の波動関数は、$L$ 個の量子ビットを用いて $2^L$ 個の点からなる離散格子にマッピングされ、位置 $x$ は領域 $0 \le x < 1$ にスケーリングされる。二乗運動量演算子 $\hat{p}^2$ は、量子並進演算子に基づく有限差分法を用いて近似される。並進演算子 $\hat{A}_x$（量子加算器）は、位置状態 $|j\rangle$ を $|j+\delta x\rangle$ へシフトさせ、その共役である $\hat{A}_x^\dagger$（量子減算器）は逆方向にシフトさせる。二乗運動量は、$\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$ と表現される。ここで、$L_m = \lambda_m / \delta x$ は、縮退したコンプトン波長を含む無次元パラメータである。

2. 運動エネルギーの摂動展開:
   相対論的な運動エネルギーハミルトニアン $\hat{H}_{re}^K$ は、相対論的なエネルギー・運動量関係の平方根の形式から導出される。これは、$\langle \hat{p}^2 \rangle \ll m^2 c^2$ の領域で有効な、高次の運動量項（$\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$）の摂動級数へと展開される。全相対論的運動エネルギーは、クライン-ゴルドン方程式から導かれた係数を取り入れつつ、この級数を打ち切ることで近似される。

3. 境界条件の取り扱い:
   本フレームワークは、2つの異なる境界条件に対処する。

• 周期境界条件 (PBC): 並進演算子は、格子をラップさせることでユニタリ性を保持する（$\hat{A}_x |2^L-1\rangle = |0\rangle$）。運動量の偶数次の期待値は、$\langle \hat{A}^{(l)} \rangle$ の組み合わせとして表現される。
• ディリクレ境界条件 (DBC): ここでは、波動関数はラップしない（$\psi(0) \neq \psi(1)$）。著者らは、境界を考慮するために補正演算子 $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$ を導入する。DBCにおける期待値は、PBCの項に、格子の境界における確率振幅に依存する $\hat{E}_0$ および関連する演算子（$\hat{E}_1, \hat{E}_2$）を含む追加項を加えることで定式化される。

4. 量子シミュレーションアルゴリズム:
   著者らは、必要な期待値を推定するための量子回路を設計している。

• 運動エネルギー: 制御並進ゲートが制御量子ビットに対して反復的に適用される。制御量子ビットをZ基底で測定することで、$\langle (\hat{A}_x^\dagger)^l \rangle$ の実部が得られ、これが必要な運動量の項に対応する。
• ポテンシャルエネルギー: ポテンシャルを表す密度行列 $\hat{\rho}_V$ を用い、系の状態の確率を通じて $\langle \hat{V} \rangle$ を推定する回路を用いる。
• 変分最適化: 全エネルギーは、アンザッツ状態 $|\psi(\{\bar{\kappa}\}, x)\rangle$ における一連の変分パラメータ $\{\bar{\kappa}\}$ を最小化することによって推定される。

主な貢献

• 第一量子化のレシピ: 本論文は、相対論的シミュレーションのための具体的なレシピを第一量子化を用いて提供しており、特に非線形性を伴うシステムに有益な、広く用いられている第二量子化法に代わる選択肢を提示している。
• 統一された境界フレームワーク: 周期境界条件とディリクレ境界条件の両方を扱うための統一的なアプローチを提示している。決定的な点として、PBCで得られたデータが、特定の境界補正項を加えることでDBCのエネルギーを推定するために再利用できることを示しており、計算オーバーヘッドを削減している。
• 近未来デバイスへの適合性: 本手法は、標準的な量子回路プリミティブ（制御加算器、アダマールゲート、および測定）を利用しており、近未来の量子デバイスへの実装に適している。
• 明示的な回路設計: 並進演算子の期待値および摂動展開に必要な境界補正項を計算するための詳細な回路設計を提供している。

結果と主張
本論文は、近未来の量子デバイス上で相対論的効果をシミュレートするための「実用的な経路」を提供すると主張している。理論的な結果は以下の通りである。

• 特定の次数の摂動までの相対論的運動エネルギー項は、並進演算子（$\langle \hat{A}^{(l)} \rangle$）と境界演算子（$\langle \hat{E}_d \rangle$）の期待値の線形結合に分解できる。
• これらの展開係数は、物理パラメータ（$m, c, \hbar, \delta x$）に基づいて事前計算が可能である。
• 本手法は、アンザッツ状態の変分最適化を通じて、非相対論的および相対論的な両方の基底状態エネルギーを推定することを可能にする。

意義
著者らは、本研究を将来の量子物理学および量子化学の発展を支援するための基礎的な一歩として位置づけている。相対論的効果を近未来のハードウェア上でシミュレート可能にすることで、現在の量子シミュレーション能力におけるギャップを埋めるものである。本論文は、格子解像度を高めるために量子ビット数を増やすことは可能であるが、摂動展開が有効であり続けるように、無次元パラメータ $L_m$ が十分に小さくなるよう注意する必要があることを強調している。本研究は、ディラック方程式の完全な解決や極限的な相対論的領域を扱うことを目的としたものではなく、非相対論的な枠組みに対する相対論的な補正が重要となる領域に特化して焦点を当てている。今後の課題として、高次のラプラシアン演算子の探索や、より正確なアンザッツ構築手法の検討が示唆されている。