Fast convergence of Majorana Propagation for weakly interacting fermions

技術要約：弱相互作用フェルミオンにおけるマヨラナ伝播の高速収束

問題設定
局所的な観測量のハミルトニアン発展下での時間ダイナミクスをシミュレーションすることはBQP完全な問題であり、これは制限された設定を除いて、効率的な古典アルゴリズムは一般に期待できないことを意味する。本研究は、$N$個のマヨラナモード $\gamma_1, \dots, \gamma_N$ に作用する四次ハミルトニアン $H$ で記述される相互作用するフェルミオン系に焦点を当てている。目的は、観測量 $A$ のハイゼンベルク発展 $A(t) = e^{iHt}Ae^{-iHt}$ を近似することである。厳密な発展は、直ちに $t>0$ において次数 $N$（非ゼロの係数を持つ最長のマヨラナ文字列の長さとして定義される次数）の観測量を生成し得るが、本研究では、低次の近似が存在する場合の設定において、それを効率的に見つけられるかどうかを調査している。

手法：マヨラナ伝播 (Majorana Propagation, MP)
本論文は、トロッター化と次数ベースの切り捨てを組み合わせたマヨラナ伝播 (MP) アルゴリズムを分析している。アルゴリズムは以下のように進行する：

1. ハミルトニアンの分解: 四次ハミルトニアン $H$ を、$G \le 4\Delta$ 個のグループ $H = H_1 + \dots + H_G$ に分割する。ここで $\Delta$ はスパース性パラメータであり（各マヨラナモードは高々 $\Delta$ 個の項に含まれる）、各グループ $H_g$ 内の各マヨラナ文字列は互いに支持集合が離れており、したがって互いに可換である。
2. トロッター化: 時間発展をステップサイズ $\delta t$ で離散化する。発展演算子は、一次のトロッター積によって近似される：$\tau^{\text{trot}}_{\delta t}(A) = e^{iH_G \delta t} \dots e^{iH_1 \delta t} A e^{-iH_1 \delta t} \dots e^{-iH_G \delta t}$。
3. 伝播と切り捨て: 観測量は各トロッターステップを通じて伝播される。極めて重要な点として、各ステップの後、観測量の展開における次数が閾値 $\ell$ を超えるすべてのマヨラナ文字列を破棄する切り捨て写像 $\text{Trunc}_\ell$ が適用される。
4. 反復: このプロセスを $t/\delta t$ 回繰り返し、$A(t)$ を近似する。

主要な技術的貢献
本論文は、MPアルゴリズムを支える2つの主要な理論的柱を確立している：

1. 次数依存のノルム不等式: 著者らは、システムサイズ $N$ ではなく、観測量の次数に依存するフロベニウスノルムの交換子およびトロッター誤差に関する新しい境界を導出した。

   • 交換子境界: $\Delta$-スパースな四次ハミルトニアンおよび次数 $d$ の観測量 $A$ に対して、$\|[H, A]\|_F \le 2\Delta\sqrt{d(d+2)}\|A\|_F$ が成立する。これは、標準的な境界が $N$ に比例するのと対照的である。
   • トロッター誤差境界: 低次の観測量に対する厳密な発展とトロッター近似との誤差は、$N$ に依存せず、$d, \Delta, G$ に依存することが示されている。具体的には、$\|\tau_t(A) - \tau^{\text{trot}}_t(A)\|_F \le C \cdot (d\Delta G)^2 t^2$ である。

2. 収束解析: 主要な定理（定理 1.1/3.3）は、MPアルゴリズムの誤差が、トロッター離散化誤差と、真の時間発展した観測量の「最良の次数-$\ell$ 近似」との誤差の和によって抑えられることを証明している。
   $$\|A^{\ell, \delta t}_{\text{MP}}(t) - A(t)\|_F = O\left( \ell^2 t \delta t + (t/\delta t)\eta^\star_F(H, A, t, \ell) \right)$$
   ここで $\eta^\star_F$ は最良の次数-$\ell$ 近似器への距離を定量化する。この結果は、低次の近似が存在する場合に、MPがそれを効率的に見つけるという最初の証明可能な保証を提供する。

弱相互作用ハミルトニアンに対する結果
本論文は、これらの結果を $H = H_0 + uV$（$H_0$ は二次形式、$V$ は四次形式、$u \le 1$ は小さな相互作用強度）という形式のハミルトニアンに適用している。

• 低次近似器の存在: 摂動論を用いて、著者らは $t < t_{\max}(u) \propto \log(e/u)/\Delta$ の時間において、$A(t)$ が次数 $\ell$ に対して誤差が指数関数的に減少する低次の多項式によって近似可能であることを証明している。
• 効率的なシミュレーション: これらの近似器の存在とMPの収束保証を組み合わせることで、著者らは（系従 1.3）、正規化されたフロベニウスノルムにおいて誤差 $\varepsilon$ を達成するための実行時間が $N^{O(\log((1+t)/\varepsilon))}$ で、任意の $\Delta$-スパースな弱相互作用ハミルトニアンのダイナミクスを時間 $t_{\max}(u)$ までシミュレートできることを証明している。
• 極限挙動: 相互作用強度 $u \to 0$ となるにつれ、$t_{\max}(u) \to +\infty$ となり、二次形式（自由）フェルミオン系が一定の次数で正確にシミュレート可能であるという既知の事実を回収する。

数値検証
著者らは、1次元および2次元格子におけるフェルミ・ハバードモデルを用いて、これらの理論的知見を検証している。数値実験は以下を示している：

• 切り捨て次数 $\ell$ が増加するにつれて、MPの出力が理想的なトロッター発展へと指数関数的に収束すること。
• 固定された $\ell$ によるシミュレーションが、相互作用強度 $U$ が増加するにつれて、厳密な解から乖離し始める特性的な時間 $t_{\max}(u)$ の挙動。
• 様々な相互作用強度における2次元格子でのホール密度ダイナミクスの高いシミュレーション精度。

意義と主張
本論文は、ハミルトニアン発展のためのマヨラナ伝播の効率性に関する最初の証明可能な保証を提供すると主張している。これは、高次の項を切り捨てるアルゴリズムが、弱相互作用系において長い時間軸にわたって正確であるという直感を定式化したものである。本研究は、幾何学的に局所的な設定や短時間のダイナミクスを超えて、古典的なシミュレーション技術の適用範囲を、相互作用グラフのスパース性と相互作用の弱さに依存することで拡張している。著者らは、結果はフェルミオン系から導出されているものの、スピン系（パウリ伝播）にも容易に一般化できると述べている。その意義は、古典的シミュレーションが依然として効率的である特定の領域（弱相互作用、スパースな四次ハミルトニアン）を特定した点にある。これは、量子優位性の研究におけるベンチマークを提供することになる。