2025-11-21T20:49:16.205065

Slow dynamics from a nested hierarchy of frozen states

MariÄ, Paljk, Zadnik

We identify the mechanism of slow heterogeneous relaxation in quantum kinetically constrained models (KCMs) in which the potential energy strength is controlled by a coupling parameter. The regime of slow relaxation includes the large-coupling limit. By expanding around that limit, we reveal a \emph{nested hierarchy} of states that remain frozen on time scales determined by powers of the coupling. The classification of such states, together with the evolution of their Krylov complexity, reveal that these time scales are related to the distance between the sites where facilitated dynamics is allowed by the kinetic constraint. While correlations within frozen states relax slowly and exhibit metastable plateaus that persist on time scales set by powers of the coupling parameter, the correlations in the rest of the states decay rapidly. We compute the plateau heights of correlations across all frozen states up to second-order corrections in the inverse coupling. Our results explain slow relaxation in quantum KCMs and elucidate dynamical heterogeneity by relating the relaxation times to the spatial separations between the active regions.

academic

凍結状態のネストされた階層構造からの遅い動力学

基本情報

論文ID: 2510.03159
タイトル: Slow dynamics from a nested hierarchy of frozen states
著者: Vanja Marić, Luka Paljk, Lenart Zadnik (リュブリャナ大学)
分類: cond-mat.stat-mech, quant-ph
発表時期: 2025年10月 (arXiv プレプリント)
論文リンク: https://arxiv.org/abs/2510.03159

要約

本論文は、量子動力学制約モデル(KCMs)における遅い不均質緩和のメカニズムを明らかにしており、ポテンシャル強度は結合パラメータによって制御される。大結合極限の周りの展開を通じて、結合強度の冪乗で決定される時間スケール上で凍結したままである状態のネストされた階層構造を明らかにする。これらの状態の分類とそのKrylov複雑度の進化は、時間スケールが動力学制約によって動力学を促進することが許可されるサイト間距離に関連していることを示唆している。凍結状態内の相関は緩やかに緩和し、結合パラメータの冪乗によって設定される持続時間を持つメタ安定プラットフォームを示し、一方、他の状態の相関は急速に減衰する。

研究背景と動機

問題定義

中心的問題：量子動力学制約モデルにおける遅い緩和と動力学的不均質性の微視的メカニズムを理解すること
重要性：
- 量子技術の急速な発展は、量子多体系における遅い緩和メカニズムの理論研究に強い動機を与えている
- KCMsはRydberg原子プラットフォームで実験的に実現可能であり、重要な応用価値を持つ
- これらのシステムは遅い緩和を示すだけでなく、豊かな非平衡現象もサポートしている

既存手法の限界

先行研究は主に量子多体スカーと強結合極限における強化された動力学制約を通じて動力学の遅さを説明していた
全体的な緩和時間階層構造の微視的メカニズムに対する完全な理解が欠けていた
すべてのプラットフォームの時間スケールと高さの統一的なフレームワークを確立できていなかった

研究動機

本論文は強結合極限を超え、プラットフォーム時間スケールと高さを決定でき、全体的な緩和時間の階層構造を説明できる完全なフレームワークを確立することを目指している。

中核的貢献

大結合展開の完全なフレームワークの確立：Hamiltonianの完全な大結合展開に基づいて、ネストされた計算基底状態の階層構造を定義した
凍結状態の分類法則の解明：Ckが任意の2つの↓スピンが少なくともk個の↑スピンで分離されているスピン配置を含むことを証明した
時間スケール階層の決定：Krylov複雑度分析を通じて、level-k状態が時間スケールt∼Δᵏ上で凍結したままであることを証明した
相関プラットフォーム高さの計算：凍結状態における相関関数プラットフォーム高さの解析表現を提供し、Δ⁻²次の修正まで正確である
空間-時間相関の確立：緩和時間を動力学を促進する活性領域間の空間分離に関連付けた

方法の詳細

タスク定義

1次元XPXモデルの動力学を研究する： $H_{XPX} = \sum_{j=1}^L \sigma^x_{j-1}(1-\sigma^z_j)\sigma^x_{j+1} + \Delta\sigma^z_j$

ここでΔはポテンシャル強度を制御する結合パラメータであり、|Δ|>1の場合、システムは遅い動力学を示す。

モデルアーキテクチャ

1. 大結合展開

再帰的スキームを採用して大結合展開を計算し、反Hermite演算子Sとしてのアルゴリズムを得る： $H_{XPX} = e^{-S}He^S, \quad H = H_F + \Delta V, \quad [H_F,V] = 0$

ここで：

$S = \sum_{n≥1} \Delta^{-n}S_n$
$H_F = \sum_{n≥0} \Delta^{-n}H_{F,n}$
切断： $H^{(k)} := \Delta V + \sum_{n=0}^{k-1} \Delta^{-n}H_{F,n}$

2. 凍結状態の階層構造

凍結状態集合を定義する： $C_k := \{|s⟩ ≡ |s_1...s_L⟩ | H^{(k)}|s⟩ = ε(s)|s⟩\}$

ネストされた階層を形成する： $C_1 ⊇ C_2 ⊇ ... ⊇ C_k ⊇ ...$

3. 主要な発見

凍結状態の特性：真空状態を除き、Ckは任意の2つの↓スピンが少なくともk個の↑スピンで分離されている配置を含む
状態数公式： $|C_k| = 1 + \sum_{N=⌈kL/(k+1)⌉}^L \frac{L}{kN-(k-1)L}\binom{kN-(k-1)L}{L-N}$
漸近的振る舞い： $|C_k| \sim χ_k^L$ 、ここでχkはχk(χ-1) = 1を満たす

技術的革新点

完全な大結合展開：XPXモデルの任意の次数の大結合展開を初めて提供した
ネストされた階層構造：凍結状態のネストされた構造とその時間スケールとの関係を革新的に識別した
Krylov複雑度分析：Krylov複雑度を利用して状態の進化を定量化し、凍結時間スケールを検証した
空間-時間相関：促進サイト間距離と緩和時間の指数関係を確立した

実験設定

数値手法

QuSpinパッケージを使用して正確な対角化計算を実施
システムサイズ：主な結果はL=14のシステムに基づいている
結合強度：Δ = 4, 8, 16, 32などの異なる値

評価指標

時間平均相関関数： $c_t(s) = t^{-1}\int_0^t dτ c_τ(s), \quad c_t(s) := \frac{1}{L}\sum_{j=1}^L ⟨s|n_j(t)n_j(0)|s⟩$
Krylov複雑度： $K_t(ψ) = \sum_{n=0,1,2,...} n|⟨ψ(t)|B_n⟩|^2$
プラットフォーム高さ推定： $c_{pl}(s) = \frac{N_↑(s)}{L} + O(Δ^{-2})$

比較分析

異なるlevel-k状態の動力学的振る舞いを比較
理論予測と数値結果の一致を検証
異なるΔ値下でのスケーリング振る舞いを分析

実験結果

主要な結果

時間スケール階層：
- Level-k状態は時間t∼Δᵏ上で凍結したままである
- Krylov複雑度は重スケール時間t/Δᵏ下でスケーリング崩壊を示す
- 理論予測の指数時間スケール分離を検証した
相関関数プラットフォーム：
- 凍結状態の相関関数は持続的なプラットフォームを示す
- プラットフォーム高さは初期↑スピン密度によって決定される： $c_{pl}(s) ≈ N_↑(s)/L$
- 非凍結状態の相関は急速に減衰する
状態数スケーリング：
- 解析公式|Ck|の正確性を数値的に検証した
- 漸近スケーリング $|C_k| \sim χ_k^L$ を確認した
- Krylovセクタ数と凍結状態数の相関を観察した