Digital Quantum Simulation of Nonequilibrium Dynamics in the Schwinger Model under a Strong External Electric Field

技術要約：強外部電場下におけるシュウィンガー模型の非平衡ダイナミクスのデジタル量子シミュレーション

問題提起
本論文は、強電場量子電磁力学（QED）、特に強電場が真空崩壊と粒子・反粒子対生成を誘起するシュウィンガー効果といった、強電場における非平衡ダイナミクスのシミュレーションにおける課題に取り組んでいる。従来の摂動論的手法では、この非摂動的な実時間現象を信頼性高く記述することは困難である。さらに、標準的なユークリッド経路積分モンテカルロ法は、実時間発展や複素重みを扱う際に「符号問題」に直面し、統計誤差が急速に増大するという問題がある。ハミルトニアン形式は符号問題を回避できるものの、指数関数的に増大するヒルベルト空間の取り扱いを必要とする。著者らは、量子多体系の状態を直接的に、かつユニタリな時間発展を量子ビット上で表現することで、これらの制限を克服するためにデジタル量子シミュレーションを用いることを提案している。具体的な焦点は、（1+1）次元の強外部電場下における（1+1）次元格子シュウィンガー模型であり、電場クエンチ（急激な変化）に伴う有限サイズでの非平衡ダイナミクスを調査することである。

手法
本研究では、開境界条件を持つ（1+1）次元格子シュウィンガー模型に基づくデジタル量子シミュレーションプロトコルを採用している。その手法は、主に以下の3つの段階で構成される。

1. モデルの定式化と量子ビットへのマッピング：

   • 連続体理論は、コグ・サスタインド（Kogut-Susskind）のスタガード・フェルミオン形式を用いて離散化される。
   • ガウスの法則を用いてゲージ自由度を排除し、各リンク上の電場を、物質電荷分布と境界電場の関数として表す。これにより、モデルは長距離相互作用を持つ純粋なフェルミオン・ハミルトニアンへと簡約される。
   • フェルミオン演算子はジョルダン・ウィグナー変換を通じてパウリ演算子へとマッピングされ、その結果、運動エネルギー、質量項、および電場エネルギーからなる量子ビット・ハミルトニアンが得られる。

2. 状態準備（VQE）：

   • ダイナミクスのための初期状態は、ゼロ電場真空状態（$\epsilon=0$）である。
   • この状態は、ハードウェア効率的なRealAmplitudesアンザッツとSLSQP古典最適化アルゴリズムを用いた変分量子固有値ソルバー（VQE）を用いて準備される。
   • 様々な外部電場下での静的性質を研究するために、著者らはまた、VQEを用いて基底状態を求め、変分量子デフレーション（VQD）を用いて第一励起状態のエネルギーを推定している。

3. 実時間発展（トロッター・鈴木分解）：

   • システムには、外部電場が$t=0$でスイッチオンされる「クエンチ」が課される。
   • 実時間発展は、時間発展演算子$e^{-iHt}$を近似するために、二次のトロッター・鈴木分解を用いてシミュレートされる。
   • シミュレーションは、量子ノイズを含めない状態で、アルゴリズムを厳密対角化（ED）と比較検証するために、Qiskitフレームワークを用いて古典コンピュータ上で実行される。

主な貢献と結果
本論文は、格子サイズ $N=8$ および $N=12$（格子単位 $a=m=g=1$）における数値結果を提示し、量子シミュレーション・プロトコルを厳密対角化（ED）ベンチマークと比較している。

• 静的基底状態の性質：

  • VQEによって準備されたゼロ電場基底状態は、EDの結果に対して高い忠実度（$>99.99\%$）を達成している。
  • 外部電場強度 $\epsilon$ に関するスキャンにより、カイラル凝縮体におけるステップ状の挙動、および特定の臨界電場強度における基底状態と第一励起状態の間のエネルギーギャップの消失が明らかになった。これは、電場による真空構造の再構成を示しており、有限体積における一次相転移に対応している。
  • 臨界電場強度 $\epsilon_c$ は、$1/N$ に対する明確な有限サイズ依存性を示す。線形補外により、無限体積における臨界電場は $\epsilon_c(\infty) \approx 0.469$ となる。著者らはこれを、開境界条件がエネルギーギャップを減少させている影響によるものとしている。

• 非平衡ダイナミクス：

  • 保存則： 二次のトロッター発展は、全$U(1)$電荷および全エネルギーを（数値精度内で）正常に保存しており、アルゴリズムの安定性を検証している。
  • 電荷分離： 電場がオンになると、空間的な電荷分布において、開境界付近で電荷・反粒子対の周期的な生成と消滅が見られる。電荷の振動の振幅と周期は、外部電場強度とともに増加する。
  • エネルギー再分配： 全電場エネルギー $H_E(t)$ は準周期的な挙動を示し、これは電場項とフェルミオン部門との間のエネルギー交換を反映している。これは電荷分布の単純な周期性とは異なり、高エネルギーモードからの複雑な寄与を示唆している。
  • 真空の忠実度： 真空状態の忠実度 $P_{vac}(t)$ は、システムが初期真空から離れていくにつれて減少する。初期時間窓（$0 \le t/a \le 1.0$）において、この減少は有効指数関数則 $P_{vac}(t) \sim e^{-\gamma_{eff} t}$ に従う。有効減衰率 $\gamma_{eff}$ は、外部電場強度とともに単調に増加する。
  • 有限サイズ効果： $N=8$ と $N=12$ のシステムの比較により、観察されたダイナミクスの特徴（境界での電荷分離、忠実度の減衰）が小さなシステムサイズの数値的アーティファクトではなく、より大きな格子でも持続することが確認された。より大きなシステムでは、初期真空からの偏差がより速くなる。

意義と主張
著者らは、VQEベースの状態準備と（トロッター・鈴木分解による）デジタル実時間発展を組み合わせることが、強電場格子ゲージ理論の非平衡ダイナミクスを研究する上で有効であることを、本研究が示していると主張している。

• 検証： 本プロトコルは、EDと比較して、電場誘起現象（具体的には、境界での電荷分離、真空状態の忠実度減衰、および準周期的なエネルギー再分配）の主要な特徴を正常に再現している。
• 手法の有用性： 本研究は、高エネルギー物理学におけるデジタル量子シミュレーションのベンチマークとして機能しており、明示的なデジタル・トロッター発展が、有限格子システムにおける複雑な強電場ダイナミクスを効果的に捉えられることを示している。
• 今後の展望： 著者らは、このアプローチが、より大きなシステムサイズへの拡張、改良された変分アンザッツ、ノイズのある量子ハードウェアへの実装、および時間依存する外部電場の研究のための基礎を提供するものであると控えめに示唆している。彼らは、現在のノイズのあるハードウェアへの即時の実験的実現を主張しているのではなく、古典シミュレータ上でのアルゴリズムの検証が、必要なステップであることを強調している。