Spectral form factor of quadratic $R$-para-particle SYK model with Random Matrix Coupling

技術要約：ランダム行列結合を用いた二次型R-パラ粒子SYKモデルのスペクトル形式因子

問題提起
サックデフ・イェ・キタエフ（SYK）モデルは、量子カオス、ブラックホール、およびホログラフィーを研究するための典型的な系であり、通常、ランダムな$q$体相互作用を持つフェルミオンによって定義される。$q=2$（二次型）のSYKモデルは厳密に解けるが、そのスペクトル特性は、無秩序がスペクトル形式因子（SFF）におけるカオス的挙動を誘起しない限り、通常は自明であると考えられている。本論文では、R-パラ粒子を用いて構成された一般化された二次型SYKモデル（$R$-PSYK2）のSFFを調査する。これらの粒子は、ヤン・バクスター方程式を満たすR行列によって定義される一般化された交換統計に従い、標準的なボソンとフェルミオンの間を補間する。

主な課題は、ガウス型ユニタリアンアンサンブル（GUE）結合を用いた$R$-PSYK2モデルにおける指数関数的なランプ挙動の直接的な数値検証が欠如していることである。先行研究では、中間時間領域において指数関数的な成長 $K(t) \sim e^{C_0 t}$ が示唆されていたが、振動項や極めて大きなシステムサイズ（$N$）の必要性により、数値シミュレーションが困難であった。また、他のランダム行列アンサンブル、具体的にはガウス型直交アンサンブル（GOE）、ガウス型シンプレクティックアンサンブル（GSE）、およびそれらの円型アンサンブル（COE、CUE、CSE）におけるこのモデルの挙動は未探索のままであった。

手法
著者らは、解析的なランダム行列理論の手法と数値的検証を組み合わせて用いている：

1. モデルの定式化： ハミルトニアンは $H = \sum_{a,i,j} (h_{ij} - \mu \delta_{ij}) \psi^+_{i,a} \psi^-_{j,a}$ と定義される。ここで $h_{ij}$ は特定のアンサンブルから抽出されたランダム行列である。系はユニタリ変換によって対角化され、結合行列の固有値によって決定される非相互作用粒子の集まりへと還元される。
2. 分配関数とSFF： グランドカノニカル分配関数およびSFFは、$n$粒子ヒルベルト空間の縮退度 $d_n$ を用いて表現される。これは単一モード分配関数 $z_R(x)$ にエンコードされている。著者らは、2つの特定のR-パラフェルミオンのケースに焦点を当てる：例A（$z_R(x) = 1 + mx$）および例B（$z_R(x) = 1 + mx + x^2$）。
3. ガウス型アンサンブル（GUE, GSE, GOE）：
   • コヒーレント状態アプローチ： 初期時間挙動とプラトーの高さに使用されるが、中間ランプを捉えるには不十分であることが指摘されている。
   • クラスター関数アプローチ： ランプ挙動を導出するために利用される。GUEの場合、これは固有値の結合確率分布を展開し、転送行列のトレースを評価することを含む。
   • 行列式/パラフィアン定式化： SFFは、単一粒子エネルギー分布の積分から構成される行列 $\alpha$ の行列式（GUEの場合）、または行列式の平方根（GSEの場合）として表現される。
   • GOEの課題： 著者らは、GOEのケースが二重積分と非解析的なカーネル関数を伴うため、GUEやGSEよりも解析的な取り扱いが著しく困難であることを指摘している。
4. 円型アンサンブル（CUE, CSE, COE）：
   • 著者らは、円型アンサンブルの厳密な可解性を利用しており、カーネルは離散フーリエ級数に展開できる。
   • SFFは、置換の共役類に関する和として表現される転送行列形式を用いて導出される。
   • ガウス型アンサンブルと円型アンサンブルの比較は重要な手法である。CUEとGUEは、大きな$N$の極限において同じクラスター関数を共有するため、著者らはGUEの予測を検証するためのベンチマークとしてCUEの結果を使用している。

主要な結果

• GUEとCUEの対応： 例A（$z_R(x) = 1 + mx$）について、著者らはGUEとCUEの結果の間の正確な解析的一致を導出した。彼らは以下の時間領域において精密な対応関係を確立している：
  $$K_{\text{GUE}}(2t) \approx K_{\text{CUE}}(t)$$
  この時間スケーリング因子2は、両方のモデルがプラトーに達する時間を一致させる。GUEの成長率 $C_0$ は、大きな$N$の極限において、CUEランプの成長率のちょうど半分であることが示されている（$2\lambda_{\text{GUE}} = \lambda_{\text{CUE}}$）。
• 成長率（$C_0$）：
  • 例Aについては、$C_0$ が解析的に導出される。$m > 1$ の場合、$C_0 = \frac{1}{2} \ln(\frac{m^2}{m^2-1})$ である。$m=1$ の場合、$C_0 \sim \frac{1}{2} \ln N$ となる。
  • 例Bについては、解析解を得ることが困難である。数値シミュレーションは、$m \ge 3$ では $C_0$ が定数に収束することを示唆しており、$m=1, 2$ のケースでは異なるスケーリング挙動を示す。
• シンプレクティック・アンサンブル（GSE/CSE）： 著者らは、適切な時間再スケーリング（GSEでは $t_S \approx 1.36N$）の後、GSEとCSEの間に同様の挙動を観察している。CSEの数値データは、勾配の遷移、およびGUE/CUEのアナロジーと一致する成長率を裏付けている。ただし、GSEのカーネルが行列であるため、ランプの解析的導出はより複雑である。
• 直交アンサンブル（GOE/COE）：
  • GOE： 著者らは計算方法の概略を述べているが、二重積分カーネルの複雑さとランプの解析解の欠如により、数値結果は提供していない。
  • COE： CSEと同様の行列カーネルを用いた計算方法が概説されているが、数値結果は省略されている。
• 非エルミート系： 著者らは固有値が複素数であるケース（$\epsilon_j = e^{i\theta_j}$）を簡潔に検討し、SFFが時間依存しない（自明である）ことを発見した。これは、$\epsilon_j = \theta_j$ であるエルミート系のケースとは対照的である。

意義と主張
本論文は、一般化された統計（R-パラ粒子）を持つ二次型モデルにおける量子カオスの理解のための統一的な枠組みを提供することを主張している。その主な貢献は以下の通りである：

1. ベンチマーク： $K_{\text{GUE}}(2t) \approx K_{\text{CUE}}(t)$ という正確な対応関係を確立することにより、著者らは、振動や有限サイズ効果によってランプの検出が困難であった直接的な数値シミュレーションの限界を克服し、GUEの結果に対する厳密なベンチマークを提供した。
2. 一般化： 本研究は、標準的なフェルミオンからR-パラ粒子へ、またGUE結合から3つのガウス型および3つの円型アンサンブルへと、二次型SYKモデルの研究を拡張した。
3. 方法論的洞察： 円型アンサンブルのアプローチが、特にカオスのランプが発生する中間時間領域において、より複雑なガウス型アンサンブルのプロキシ（代理）として正確な結果を導き出すための強力なツールであることを示している。
4. 限定的な範囲： 著者らは、相互作用のある（$q>2$）ケースを解いたり実験的な提案を行ったりするのではなく、無秩序と一般化された統計によって誘起されるカオスの本質的な特徴を分離するために、二次型ケースに焦点を当てることを明示している。これは、相互作用のあるカオス系への将来の研究に向けた、モデルの単純化（例：交換可能な演算子や特定のスペクトル分布の使用）が実り多い道であることを示唆している。

結論として、直交アンサンブルは技術的な困難を伴うものの、円型アンサンブルのアプローチがR-PSYK2モデルのカオス的挙動（ランプとプラトー）を正常に検証することに成功し、指数関数的なランプがこれらの一般化された二次型系における堅牢な特徴であることを確認した。