Quantum Simulation of Generalized Parton Distributions in the Schwinger Model

技術要約：シュウィンガーモデルにおける一般化パルトン分布の量子シミュレーション

問題の定義
一般化パルトン分布（GPD）は、横方向のインパクトパラメータと縦方向の運動量の相関を符号化しており、ハドロンの内部構造を理解する上で極めて重要である。しかし、量子色力学（QCD）におけるGPDの第一原理からの計算は、ユークリッド格子QCDにおける符号問題によって阻まれており、これが実時間のリコン・ライトコーン相関関数の直接的な計算を妨げている。LMET（Large Momentum Effective Theory）のような古典的手法は回避策として準PDF（quasi-PDFs）を用いるが、これらは極端な運動量ブーストを必要とし、べき補正の問題を抱えている。量子コンピューティングは、ハミルトニアン力学の直接的なシミュレーションと実時間発展を可能にすることでパラダイムシフトをもたらし、符号問題を回避できる可能性がある。それにもかかわらず、GPDのシミュレーションは、離散対称性（電荷共役 $C$、パリティ $P$、時間反転 $T$）の保持や、非ゼロの運動量を持つハドロン状態の効率的な準備に関して、依然として困難な課題となっている。

手法
著者らは、フルQCDのプロキシ（代理）として(1+1)次元のシュウィンガーモデル（量子電磁力学）におけるGPDをシミュレートするための包括的な量子アルゴリズムを提案している。本手法は、以下の主要な技術的選択によって特徴付けられる。

1. ウィルソン・フェルミオン： 現在の量子シミュレーションで主流となっているスタガード・フェルミオンとは異なり、本研究では格子離散化にウィルソン・フェルミオンを採用している。この選択は、格子上での厳密な電荷共役（$C$）対称性の保持において極めて重要である。著者らは、$C$対称性が中性メソンの研究において不可欠であり、相関関数における非物理的な実部による汚染を排除し、GPDが運動量分率 $x$ の奇関数であることを保証すると主張している。
2. ハミルトニアン定式化： シュウィンガーモデルはアキシャルゲージで定式化されている。連続体のハミルトニアンは、ウィルソン・フェルミオンを用いて離散化されており、$C, P, T$ 対称性を維持しつつフェルミオンのダブリング問題を回避するためにウィルソン項が付加されている。時間方向のゲージ場はガウスの法則によって消去され、その結果、ハミルトニアンには非局所的な相互作用項が生じる。
3. 状態準備： アルゴリズムには、量子数分解型変分量子固有値ソルバー（VQE）を利用している。特定の非ゼロ運動量を持つハドロン状態（$|h(k)\rangle$）および真空（$|\Omega\rangle$）を準備するために、著者らは一連の直交参照状態を構築している。運動量注入回路 $D(k)$ は、参照状態を確定した運動量を持つ状態へと写像する。VQEのアンスァッツ $U(\theta)$ は対称性を保持するように設計されており、無電荷セクター（$Q=0$）内での真空と励起ハドロン状態を同時に抽出するために、マルチ状態損失関数を最適化する。
4. 相関関数の測定： アルゴリズムの中核は、ウィルソン線を含むゲージ不変なライトコーン相関関数を評価することである。著者らは、補助量子ビットを用いた制御された運動量並進演算子および制御されたゲージリンク（$W_0$）を実装している。実時間発展後の補助量子ビットの非対角期待値（$\sigma_x$ および $\sigma_y$）を測定することで、空間相関関数 $\tilde{H}_{q/h}(N, \xi)$ を抽出する。
5. 複雑度解析： リソース要件は、時間発展のための2次スズキ・トロッター積公式を用いて分析されている。総複雑度は、システムサイズ $N$ と逆精度 $\epsilon$ に対して多項式オーダーでスケールし、具体的には $O(N^p + N^{5/2}/\epsilon^{1/2})$ となる。

主な貢献

• 統一されたフレームワーク： 本論文は、PDFとGPDを別個の評価対象として扱うのではなく、両方を計算可能な統一された量子アルゴリズムのフレームワークを確立している。
• 対称性の保持： ウィルソン・フェルミオンを量子コンピューティングの枠組みに実装しても計算複雑度が増大しないこと、そして優れた対称特性を提供することを実証している。具体的には、厳密な$C$対称性を保持しており、これが中性メソンの物理的に意味のある結果を得るために不可欠であることを示している。
• アルゴリズムの構築： 運動量注入（$D(k)$）、ゲージリンクの実装、およびシュウィンガーモデルに必要な対称性保持型VQEアンスァッツのための具体的な量子回路の構築を詳述している。

結果
著者らは、空間格子サイズ $N=18$（特定のチェックでは最大 $N=22$）における厳密対角化（ED）を用いて、提案手法のベンチマークを行った。

• 質量スペクトル： 計算された最も軽い $q\bar{q}$ 束縛状態の質量スペクトルは、線形外挿を通じて既知の解析的連続体極限（$m_h/g \approx 0.564$）を再現しており、ハミルトニアン定式化の妥当性を検証している。
• 相関関数： 空間ライトコーン相関関数 $\tilde{H}_{q/h}(N, \xi)$ が抽出された。結果は、格子体積が増大するにつれて相関関数の実部が漸近的に消失することを確認しており、これは$C$対称性と一致している。
• スタガード・フェルミオンとの比較： 直接比較により、ウィルソン・フェルミオンはスタガード・フェルミオンと比較して、相関関数の実部と虚部の比が著しく抑制されることが示された。これは、ウィルソン定式化が有限体積効果や対称性の破れに起因する非物理的な汚染を効果的に排除していることを裏付けている。
• GPDおよびPDF： 相関関数のフーリエ変換により、運動量空間のGPDが得られた。得られたPDF（$\xi=0$ の極限）は $x=0.5$ でピークを持ち、これは1+1次元における束縛状態の2つの価ルク（valence quark）の性質を正しく反映している。GPDは $x$ の奇関数であることが判明し、電荷中性状態の弾性的電荷フォームファクターがゼロにならなければならないという理論的要件を満たしている。

意義と主張
本論文は、本研究がフル(3+1)次元QCDシミュレーションへの重要な前駆体であることを主張している。ウィルソン・フェルミオンを計算複雑度を増大させることなく量子ハードウェア上に効率的に実装できることを示すことで、著者らは、基本対称性を厳格に遵守しながら非摂動現象をシミュレートするための堅牢な道筋を提示している。システムサイズと精度に対するアルゴリズムの多項式スケーリングは、ライトコーン相関を直接シミュレートする上での量子コンピュータの理論的優位性を強調している。著者らは、本研究を、非アーベルゲージ理論（例：SU(2) または SU(3) ヤン＝ミルズ理論）へのアルゴリズムの一般化、そして最終的には将来のフォールトトレラント量子ハードウェア上でのフルQCDへの展開に向けた基礎的なステップとして位置づけている。