Kibble-Zurek Dynamics in Two-dimensional Frustrated Systems with a Neural Foundation-state Subspace Method

技術要約：ニューラル・ファンデーション状態部分空間法を用いた2次元フラストレート系におけるキブル・ズurek力学

問題提起
強相関な2次元多体系におけるリアルタイムの量子力学へのアクセスは、依然として重大な数値的課題である。1次元系はテンソルネットワーク法（例：DMRG）によって良好に扱われるが、2次元系ではエンタングルメントの急速な増大、量子モンテカルロにおける動的な符号問題、および時間依存変分モンテカルロ（tVMC）における統計的ノイズやコンディショニングの問題に直面する。この制限は、ネール反強磁性秩序と潜在的な量子スピン液体（QSL）相との間の遷移の性質が議論されている$J_1$-$J_2$ハイゼンベルク模型のような、フラストレートした量子磁性体において特に顕著である。具体的には、静的な基底状態の研究は普遍的な臨界指数を持つ連続的な遷移を示唆している一方で、この遷移を横断する非平衡力学（キブル・ズurek（KZ）機構によって支配されるもの）は、2次元における準断熱クエンチのシミュレーションの困難さゆえに、これまでほとんど探索されてこなかった。

手法：ニューラル・ファンデーション状態部分空間（NFS）法
著者らは、準断熱力学をシミュレートするためのニューラル・ファンデーション状態部分空間（NFS）法を導入している。このアプローチは、低エネルギー変分部分空間を構築することで多体系のシュレディンガー方程式を射影し、問題を少数の線形係数の発展へと縮小させる。

1. ファンデーション・ニューラルネットワーク量子状態（FNQS）： トランスフォーマー・アーキテクチャに基づく単一のニューラルネットワークが、クエンチがカバーする全結合範囲（初期$\gamma_0$から最終$\gamma_f$まで）にわたる基底状態多様体を表現するために一度だけ学習される。このFNQS、すなわち$|\Phi_\theta(\gamma)\rangle$は、「ファンデーション（基礎）」状態として機能する。
2. 部分空間の構築： FNQSから、離散的な結合点に沿って選択された忠実度の高い$M$個の基底状態のサブセット（$|\phi_k\rangle = |\Phi_\theta(\gamma_k)\rangle$）が抽出される。これらの状態は、時間発展のための固定された基底を形成する。
3. 線形変分アンザッツ： 時間発展する状態は、線形重ね合わせとして次のように表される：$|\Psi(t)\rangle = \sum_{k=1}^M c_k(t) |\phi_k\rangle$。
4. 射影力学： 時間依存シュレディンガー方程式は、この固定された部分空間に射影され、係数$c(t)$に関する連立線形常微分方程式系を与える：$S \dot{c}(t) = -i H(\gamma(t)) c(t)$。
   • $S$は（時間独立な）オーバーラップ行列である。
   • $H(\gamma(t))$は、事前計算された射影演算子$\langle \phi_i | \hat{h}_\mu | \phi_j \rangle$の線形結合として構成されるハミルトニアン行列である。
5. 効率性： 基底が固定されており、演算子が事前計算されているため、本手法は時間発展中のニューラルネットワークの最適化やサンプリングを回避できる。異なるランプ時間（$\tau_Q$）やプロトコルも、部分空間が確立された後は、追加のコストをほとんどかけることなくシミュレート可能である。

主要な結果
本手法は以下の2つのモデルで検証された：

1. 2次元横磁場イジング模型：

   • $6 \times 6$格子に対する厳密対角化と比較検証を行い、臨界点まで最大相対エネルギー誤差0.05%以内の合致を示した。
   • 有限サイズスケーリング解析（$6 \le L \le 12$）により、期待されるKZスケーリングを回収した。
   • 抽出された臨界指数：$\alpha = 0.362(7)$、$\nu = 0.567(18)$、および$\eta = 0.043(29)$（固定された動的指数$z=1$において）。これらの値は、先行する動的研究および静的推定値と一致している。

2. フラストレート正方格子$J_1$-$J_2$ハイゼンベルク模型：

   • $J_2/J_1 = 0.49$における遷移を研究するために、最大$16 \times 16$（$M=28$の基底状態）のクラスターに適用された。
   • 本手法は、ネール相からスピン液体相への遷移におけるKZスケーリングの強力な数値的証拠を提供した。
   • 抽出された指数：$\nu = 1.23(15)$および$\eta = 0.409(19)$（固定された$z=1$において）。これらの動的な推定値は、静的な有限サイズスケーリングおよびレベル分光の結果（$\nu \approx 1.21, \eta \approx 0.410$）と定量的に一致しており、連続的な臨界遷移の仮説を支持している。

意義および主張
本論文は、フラストレートした$J_1$-$J_2$ハイゼンベルク模型におけるキブル・ズurek機構の最初の証拠を提供するものであると主張している。普遍的な臨界指数を非平衡力学から抽出することに成功したことで、NFS法は遷移の性質に対する独立した動的なプローブを提供し、静的な基底状態の知見を裏付けている。

著者らは、NFS法が、直接的なリアルタイム手法が現在制限されている2次元フラストレート系における動力学的臨界性のテストのための実用的な経路を確立したことを強調している。また、本手法は（系の状態が瞬時の基底状態の近傍に留まる）遅いランプおよび臨界近傍の力学には自然に適しているものの、高速なクエンチや、状態が瞬時の基底状態が捉えていない方向へと探索するようなポスト臨界の長い発展においては、その精度が低下することに注意を促している。それにもかかわらず、本研究は、直接的なリアルタイム手法が現在困難な2次元フラストレート系における動力学的テストのための実用的なルートを確立した。