Efficient classical simulation of two-dimensional long-range systems: Rydberg arrays and beyond

技術要約：二次元長距離系の効率的な古典シミュレーション

問題提起
変分モンテカルロ法（VMC）は、量子シミュレータをベンチマークするための主要なツールであり、特に長距離相互作用（例：$1/r^3$ または $1/r^6$ 結合を持つリドベリ原子アレイ）を伴う二次元（2D）系において重要です。しかし、任意の全対全（all-to-all）二体相互作用を持つ系の「局所エネルギー」（$H_{\text{loc}}$）を評価する際に、重大な計算上のボトルネックが存在します。

• 標準的なVMCでは、テンソルネットワーク状態（TNS）の場合はサンプリングコストが $O(N)$、ニューラル量子状態（NQS）の場合は $O(N^2)$ となりますが、長距離の非対角結合を持つ $H_{\text{loc}}$ の評価には、$O(N^2)$ 個の非局所項の和が必要です。
• 各項の評価に $O(N)$ 時間を要するフルネットワーク評価が必要であるため、モンテカルロ（MC）スィープあたりの総コストは $O(N^3)$ にスケールします。この立方体のスケーリングにより、大規模な2D長距離系（例：$10 \times 10$ 格子）の正確な古典シミュレーションは困難となっています。
• さらに、製品状態（積状態）からTNSを進化させる際（これはリドベリ実験で一般的です）、量子幾何テンソル（QGT）のボンド次元が1より大きい場合、結合次元にヌル方向が生じ、時間依存変分原理（TDVP）方程式が特異になるという技術的な病理が発生します。

手法
著者らは、局所エネルギー評価の計算複雑性を $O(N^3)$ から $O(N)$ へと低減しつつ、一般的な不均一結合や現実的な幾何学的構造を扱えるように設計された、適応基底テンソルネットワーク tVMC フレームワークを導入しています。

1. 適応基底サンプリング:

   • 単一の固定基底（例：$z$ 基底）でサンプリングする代わりに、ハミルトニアンを特定の基底において対角となる成分に分解します。
   • 双極子XYモデル（$H = H_{XY} + H_Z$）の場合、進化を $H_X$（$x$ 基底で対角）と $H_Y + H_Z$（$y$ 基底で対角）に分割します。
   • 独立した2つのマルコフ連鎖を維持します。一つは $|\langle x|\psi_\theta\rangle|^2$ からサンプリングし、もう一つは $|\langle y|U|\psi_\theta\rangle|^2$ からサンプリングします（ここで $U$ は基底変換演算子です）。
   • 特定の基底内では長距離の非対角結合が消失するため、各成分の局所エネルギー評価は $O(N)$ となります。
   • パラメータ更新は各基底（$\dot{\theta}_X$ および $\dot{\theta}_Y$）について個別に計算され、$\dot{\theta} = \dot{\theta}_X + U^\dagger \dot{\theta}_Y$ によって結合されます。

2. 製品状態の初期条件への対処:

   • 製品状態（ここでボンド次元 $D > 1$）から始まる際の QGT の特異性を解決するために、著者らは初期時間ステップに対して 投影 tVMC (p-tVMC) アプローチを採用しています。
   • 特異な TDVP 方程式を解く代わりに、時間発展演算子の二次の展開である $V$ を用いて、フィデリティ $|\langle \psi_\theta | V | \psi(0) \rangle|^2$ を最大化することで最初のステップ（$\delta t$）を決定します。
   • 初期状態は製品状態であるため、この最適化は効率的であり、サンプリング分布 $|\langle s|V|\psi(0)\rangle|^2$ は計算可能です。この初期ステップの後、標準的な tVMC が再開されます。

3. 対称性の活用:

   • 本フレームワークは、TNS アンザッツにおいて格子の対角線に沿った反射対称性を明示的に課しています。これにより、独立した変分パラメータの数を4分の1に削減し、GPU メモリ要件を大幅に緩和し、パラメータ更新の分散を低減させています。

主な結果
本手法は、最近のリドベリ実験（Chen et al., Nature 2023）で実現された $10 \times 10$ 双極子XYモデル の断熱実時間プロトコルのシミュレーションに適用されました。これは、従来の古典シミュレーションの手が届かなかったシステムサイズです。

• ベンチマーク: $4 \times 4$ システムにおいて、TNS-tVMC 結果（ボンド次元 $D=4, 5$）は、厳密対角化の結果に系統的に収束し、駆動による緩和ダイナミクスと高速な振動を正確に再現しました。
• $10 \times 10$ シミュレーション:
  • シミュレーションはエネルギー進化を追跡し、基底状態エネルギーのわずかに上にある近傍定常プラトーへと緩和する様子を示しました。
  • 反強磁性的（AFM）相関: 本研究は、実験における不一致に対処しました。実験で準備された動的状態は、基底状態の行列積状態（MPS）計算と比較して、長距離AFM相関が抑制されていることが示されていました。
  • 著者らの理想的なダイナミクス・シミュレーション（より大きなデチューニング $\delta_0 = 19.48$ を持つ完全なネール状態から開始）は、規定された断熱ランプが、変分基底状態と一致する空間的に広がったAFM相関を実際に生成することを示しました。
  • これは、実験で観察された抑制が、理想的な断熱プロトコルの本質的な失敗ではなく、実験的な不完全さ（例：$10 \times 10$ アレイにおいて最適な $\delta_0$ の使用を妨げるレーザー出力の制限など）に起因している可能性が高いことを示唆しています。
• 収束性: 結果はサンプルサイズ（$N_s = 81,920$ 対 $163,840$）に対して頑健な収束性を示し、適応基底アプローチの信頼性を裏付けました。

意義および主張
本論文は、現在切実にツールを必要としているリドベリやイオントラップなどの長距離量子プラットフォームのための、スケーラブルな古典ベンチマークツールを提供すると主張しています。

• 技術的成果: $O(N^3)$ のボトルネックを解消し、$O(N)$ のコストで2D長距離系の実時間ダイナミクスおよび基底状態の正確なシミュレーションを可能にしました。
• 病理の解決: 製品状態から TNS を進化させるための堅牢なソリューションを提供しています。これは、実験プロトコルにおける一般的なシナリオです。
• 科学的影響: $10 \times 10$ リドベリ実験に対する直接的な古典的参照を提供することで、実験結果と理論的な基底状態予測との間の不一致の性質を明らかにしました。これは、理想的な断熱プロトコルが期待される長距離秩序を生成できることを証明しており、実験の最適化の焦点を、プロトコルの失敗から実現上の不完全さへとシフトさせました。
• 今後の展望: 著者らは、本手法を古典と量子のコラボレーションの架け橋として位置付けています。すなわち、エンタングルメントの成長がテンソルネットワークにとって過度に深刻になる前に、古典シミュレーションが興味深い物理的セットアップの特定をガイドし、その後、ダイナミクスを量子シミュレータに引き継ぐという流れです。また、本手法は、他の数値的手法が苦戦する、符号問題を持つ2D長距離系の研究における新たな道としても注目されています。