Efficient computation of real-time correlators using Pauli Propagation

技術要約：パウリ伝播を用いた実時間相関関数の効率的な計算

問題提起
実時間二点時間順序相関関数 $C(t) = \text{Tr}[\rho A(t)B]$ の計算は、多体系物理学における根本的な課題であり、微視的なハミルトニアンの詳細と、中性子散乱や電気伝導率などのマクロな実験的観測量との間の架け橋となるものである。従来の計算手法は大きなボトルネックに直面している。すなわち、厳密対角化はメモリ制限により限定され、実時間量子モンテカルロ法は符号問題に苦しみ、テンソルネットワーク法（TEBDやTDVPなど）は、時間発展に伴うボンド次元の急速な増大により、1次元を超えたスケールへの拡張が困難である。パウリ伝播（PP）は、ハイゼンベルク描像で直接演算子を時間発展させる有望な代替案であるが、本質的にパウリ文字列の指数関数的な増大に制約されており、これが切り捨て（truncation）によるバイアスを生じさせ、アクセス可能な時間窓を制限している。

手法
著者らは、短時間のパウリ伝播と高度な信号処理および基底状態近似技術を組み合わせることで、直接的なシミュレーション窓を超えて相関関数を拡張するハイブリッド計算戦略を提案している。この手法は、主に以下の3つのコンポーネントで構成される。

1. 切り捨てられたパウリ伝播（Truncated Pauli Propagation）: 時間発展した演算子 $A(t)$ は、パウリ文字列の基底を用いて $A(t) = \sum_i a_i(t) P_i$ と展開される。パウリ文字列の指数関数的な増大を管理するため、著者らは2つの切り捨て戦略を採用している。

   • 係数カットオフ ($\epsilon$): 係数が $|a_i| < \epsilon$ である文字列を破棄する。
   • 重みカットオフ ($w$): パウリの重み（非単位サイト数）が $w$ を超える文字列を破棄する。
     発展はトロッター化（Trotterization）を通じて行われ、各ステップにおいて、ハミルトニアンの項との交換関係に基づいてパウリ文字列が分岐する。

2. Clifford拡張行列積状態（Clifford-Augmented Matrix Product States: CAMPS）: もつれ状態にある基底状態（特に2次元の場合）に対する期待値 $\langle \Psi | A(t) | \Psi \rangle$ の評価の困難に対処するため、著者らはCAMPSを利用する。この手法では、行列積状態（MPS）にクリフォード変換を適用することで、パウリ和の係数構造を保持しつつ、必要なボンド次元を大幅に削減できる。これにより、標準的なMPS法では計算的に極めて困難となる複雑な基底状態に対する期待値の効率的な評価が可能になる。

3. 正定値（PSD）拡張: 相関関数 $C_{AA}(t) = \text{Tr}[\rho A^\dagger(t)A(0)]$ が正の関数であることを認識し、その離散時間トーピッツ行列が正定値（PSD）であるという性質を利用する。Carathéodory、Toeplitz、およびPisarenkoの研究に基づき、ランク $r$ のPSDトーピッツ行列は、$2r$ 個のノイズレスなデータ点によって一意に決定される。著者らは、短時間のPPデータを「超解像」するために低ランクPSDフィットを用い、さらなる演算子の時間発展を行うことなく、相関関数を遅い時間へと拡張する。これは、多くの物理系が少数の特徴的な周波数によって支配されているという観察に基づいている。

主な結果
著者らは、1次元および2次元ハイゼンバーグモデルを用いて、この複合的なアプローチの有効性を実証している。

• 1次元強磁性ハイゼンベルク鎖 ($N=64$): 攻撃的な切り捨てを行っても（時間トレースが厳密な結果と定性的に一致する程度であっても）、フーリエ変換されたスペクトルが正しい物理的特徴（マグノン分散など）を保持していることを示している。アドホックな減衰関数またはPSD拡張を用いることで、生の時間トレースにおける長時間の動的精度が失われても、スペクトル特性を回復できる。
• 1次元反強磁性ハイゼンベルク鎖 ($N=64$): DMRGで準備された基底状態を用い、PPをTDVPと比較した。攻撃的なPPの切り捨ては、MPS-TDVPと比較して大幅な周波数広がりと解像度の低下を招くが、本手法はスピノン連続体の一般的な特徴を捉えることに成功している。この研究は、もつれ状態にある基底状態の場合、ボトルネックが演算子の発展から期待値の評価へと移行することを浮き彫りにしている。
• 2次元反強磁性ハイゼンベルク格子 ($8 \times 8$): これは標準的なテンソルネットワーク（TNS）手法が非常に高価になる領域である。PPをCAMPS（ボンド次元 $D=240$）とランク1のPSD拡張と組み合わせることで、著者らはわずか7ステップの時間発展のみで動的スピン構造因子 $S(k, \omega)$ を計算した。得られたスペクトルは線形スピン波理論とよく一致しており、$\Gamma$ 点におけるゴールドストーンモード、 $X$ 点への分散、および $M$ 点での軟化を正しく特定している。
• スケーラビリティ: 2次元の場合、限界は（短時間の拡張によって回避される）演算子の無制限な増大ではなく、基底状態の準備と期待値の評価にあると著者らは指摘している。さらに、付録において、本手法が $100 \times 100$ の強磁性格子におけるマグノンスペクトルを回復できることを示している。

意義と主張
本論文は、この複合的なアプローチが、古典的な手法が信頼性を持って相互作用する量子多体系のリアルタイム力学を探索できる領域を広げると主張している。

• 指数関数的な増大の克服: 演算子のサポートが管理可能な範囲に収まる短時間のシミュレーションに限定し、物理的制約（正値性と低ランクの周波数支配）を用いてデータを拡張することで、パウリ演算子の指数関数的な増大を回避する。
• 高次元における効率性: CAMPSを統合することにより、著者らは、平衡相関関数が標準的なMPSの時間発展に伴う急激なボンド次元の増大を回避しつつ、他の手法が扱える限界領域の2次元系に対して計算可能であることを示した。
• スペクトルの忠実度: 高解像度のスペクトル情報が、極めて限られた時間データ（2次元の場合、わずか7ステップ）から抽出できることを結果は示しており、これはシステムのダイナミクスが限られた数の周波数によって支配されている場合に限られる。

著者らは、本手法の限界を謙虚に認めている。すなわち、低ランクの仮定は、連続的な周波数を持つ系（1次元スピノン連続体など）では成立しないこと、および、本手法の成功は、基礎となる周波数を特定できるほど切り捨てが十分に正確であるかどうかに依存することである。彼らは、この研究を1次元におけるテンソルネットワークの代替としてではなく、現在の最先端の手法が苦慮している大規模かつ高次元の格子における、動的相関関数を計算するための有望で実行可能な道筋として位置づけている。