Diameter truncated operator evolution

技術要約：直径限定演算子進化 (Diameter-Truncated Operator Evolution, DTOE)

問題提起
非平衡多体系量子系のダイナミクスをシミュレートすることは、複雑性の急速な増大により根本的に困難である。系が進化するにつれ、もつれ（エンタングルメント）、非スタビライザー性、および非ガウス性が増大し、当初は局所的であった演算子が高度に非局所的で「スクランブル（かき混ぜられた）」状態となる。この成長により、従来のテンソルネットワークによるシミュレーションは、波動関数やあらゆる観測量の完全な動的進化を効率的に捉えることができず、極めて短時間のタイムスケールを超えて計算不能となる。後期的な熱化によって流体力学的記述が可能になる一方で、短時間および中間時間における平衡への接近を正確に計算することは依然として困難である。

手法
著者らは、有界な空間的広がりを持つ演算子の部分空間にシミュレーションを制限することで、演算子のダイナミクスをシミュレートするために設計された数値的手法、直径限定演算子進化 (DTOE) を導入する。

• 核となる概念: 演算子の重み（パウリ展開における非自明なパウリ項の数）に基づいて演算子を切り詰める既存の手法（DAOEやPauli Propagationなど）とは異なり、DTOEは直径に基づいて演算子を切り詰める。直径は、格子上の演算子が非自明にサポートされている実空間領域のサイズとして定義される。
• 実装: 本手法は、1+1次元ブリックワーク量子回路の枠組み内で定式化されている。演算子を倍化したヒルベルト空間へとベクトル化する「折り畳み（folded）」表記を利用する。時間発展は、直径 $d$ の演算子間の区間を写像する転送行列を適用することによって計算される。
• 切り詰めスキーム: ハードな切り詰めが採用されている。直径が固定のカットオフ $d$ を超えるすべての演算子は、シミュレーションから厳密に除外される。相関関数は、これらの転送行列を連鎖させることで、局所演算子間のパスを生成することにより再構成される。
• 理論的動機: このアプローチは、演算子のノルムの大部分は高度に非局所的な軌道へと流れ込む一方で、局所的な二点相関関数は、しばしば有界で一定の直径領域内に留まる「非典型的な」軌道によって支配されるという観察に基づいている。一般的な回路において、時空における「容易な（easy）」フェーズ（通常は臨界速度 $v_c$ 以上）が存在し、そこでは相関は $O(1)$ の一定直径を持つ演算子によって支配される。このレジームにおいて、DTOEは多項式時間計算量 $O(t^2)$ と、定数メモリスケーリング $O(t)$ を提供する。

主要な貢献と結果
本論文は、高性能計算リソースを使用することなく、可積分系および非可積分系の両方に対する広範な数値シミュレーションを通じて、DTOEの有効性を実証している。

1. 摂動を受けたデュアル・ユニタリ回路（キック・イジング模型）:

   • 著者らは、パラメータ $J=B=\pi/4+\eta$ および $h=1.2$ を持つキック・イジング模型を、摂動の強さ $\eta \in [0, 0.5]$ を変化させてシミュレートした。
   • 無限温度における二点相関関数 ($C_{XX}$) を、最大 $t=100$ 回の回路レイヤーまで計算した。
   • 結果: $\eta = 0.24588$ において、低直径近似 ($d=1, 2$) は厳密なデータから乖離したが、より高い直径 ($d \ge 9$) では $10^{-2}$ から $10^{-3}$ オーダーの相対誤差を達成した。本手法は、非常に非摂動的なレジームにおいても効果的であり、通常 $\sim 25$ フロッケサイクル程度しか扱えない現在のデジタル量子コンピュータのシミュレーションを凌駕している。

2. ハイゼンベルクXXZ模型:

   • 相関幅のスケール $\sigma(t) \sim t^\alpha$ に対する流体力学的輸送指数 ($\alpha$) を抽出するために、本手法をXXZ回路に適用した。
   • 弾道的レジーム ($|\Delta| < 1$): $d=10$ のDTOEは、$\alpha = 1.00 \pm 0.01$ を正確に抽出した。著者らは、手法が拡散へとドリフトする前の正確なタイムスケール $t^*$ が、直径 $d$ に対して指数関数的にスケールすること ($t^* \sim O(\exp(d))$) を観察した。
   • 超拡散レジーム ($|\Delta| = 1$): 等方的点において、本手法は理論的に予測された指数 $\alpha = 2/3$（具体的には $d=10$ で $\alpha \in [0.6663, 0.6684]$）を捉えた。これは、このレジームにおける長い準粒子ストリングの捕捉が困難であるという既知の事実にもかかわらず達成された。
   • 拡散レジーム ($|\Delta| > 1$): 本手法は、低い直径においても拡散指数 $\alpha = 1/2$ を正確に再現した。これは、このレジームにおける輸送が短い準粒子ストリングによって支配されているためである。

意義と主張
本論文は、DTOEが重みベースの手法と比較して、より「軽量な（leaner）」切り詰めプロトコルを提供すると主張している。理論的解析によれば、格子サイズ $L$ に対して、直径 $d$ を持つ演算子の数は $O(L)$ でスケールするのに対し、重み $d$ を持つ演算子の数は $O(L^d)$ でスケールする。したがって、大きな $L$ に対して、DTOEと重みベースの切り詰めにおける保持される演算子の比率は消失し、DTOEが本質的に効率的であることを示唆している。

著者らは、DTOEが局所的な相関が有界な直径の軌道によって支配される系に対して、物理的に妥当な動機付けを持つと主張している。彼らは、DTOEが現在の量子ハードウェアの能力を超えるタイムスケールにおいて、輸送特性や相関関数を抽出することに成功した点を強調している。さらに、彼らの結果は1+1次元で提示されているが、既存の文献は2次元の一般的な回路において時空格子全体が「容易な」フェーズに入ることを示唆しており、これはDTOEがより高次元において特に強力である可能性を示している（ただし、この一般化は本研究では明示的にテストされていない）。また、DAOEは特定の散逸超演算子を持つMPO手法として構成されているが、DTOEを同様に構成することは今後の課題であるとも述べている。