Two-parameter Family-Vicsek scaling in a dissipative XXZ spin chain

以下は、提供された論文「Two-parameter Family–Vicsek scaling in a dissipative XXZ spin chain（散逸性 XXZ スピン鎖における二パラメータ型ファミリー・ヴィセク・スケーリング）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

非平衡状態における多体物理の普遍的な特徴を理解することは、現代物理学の中心的なテーマの一つです。特に、孤立量子系における「量子クエンチ（急激なパラメータ変化）」は、相関や揺らぎの非自明な時空間広がりを引き起こします。古典的な表面成長現象などで確立されたファミリー・ヴィセク（Family-Vicsek: FV）スケーリングは、揺らぎの粗さ（roughness）が時間とともに成長し、有限サイズ効果により飽和する様子を記述するスケーリング則として知られています。

しかし、量子系において散逸（環境との相互作用、ここではゲインとロスの過程）が存在する場合、この FV スケーリングがどのように変化するか、あるいは保存則が破れた場合にどのような普遍的な振る舞いが現れるかは十分に解明されていませんでした。本研究は、散逸性 XXZ スピン鎖において、非平衡定常状態（NESS）からのクエンチ後の転送されたセグメント磁化の揺らぎ（粗さ）が、どのように FV スケーリングに従うか、あるいは散逸によってどのように修正されるかを解明することを目的としています。

2. 手法 (Methodology)

• モデル: 一階の XXZ ハミルトニアン（一軸異方性 $\Delta$ を持つ）に、局所的な一様ゲイン（$\gamma_p$）とロス（$\gamma_l$）を加えたマルコフ過程を記述する Lindblad 方程式を用いています。
• 初期状態: ゲインとロスのバランスによって決定される解析的に扱いやすい「非平衡定常状態（NESS）」を用います。この状態は $\sigma^z$ 基底対角化された積状態であり、スピンの化学ポテンシャル（fugacity）$\zeta = \gamma_p/\gamma_l$ によって特徴づけられます。
• 観測量: セグメント（長さ $\ell$）内の磁化 $\Sigma_\ell$ の時間変化 $\Delta\Sigma_\ell(t) = \Sigma_\ell(t) - \Sigma_\ell(0)$ の揺らぎを定義し、その 2 次累積量から粗さ $W(\ell, t) = \sqrt{\kappa_2}$ を計算します。
• 数値・解析手法:
  • 量子生成関数（QGF）: 転送された磁化の統計を捉えるために、2 時間量子生成関数 $G_\ell(\lambda, t)$ を導入しました。これは仕事統計や全カウント統計の枠組みに類似しています。
  • ベクトル化 TEBD: 密度行列を Liouville 空間の状態ベクトルとして表現し、テンソルネットワーク（TEBD: Time-Evolving Block Decimation）アルゴリズムを用いて Lindblad 時間発展をシミュレーションしました。これにより、QGF の累積量を高精度に抽出しました。
  • 非相互作用極限の解析解: 相互作用がない場合（$\Delta=0$、XX 鎖）、Jordan-Wigner 変換により自由フェルミオン系に写像し、第三量子化法を用いて厳密解を導出しました。

3. 主要な貢献と結果 (Key Contributions & Results)

A. 非相互作用極限（$\Delta=0$）における厳密解と二パラメータ・スケーリング

非相互作用極限において、散逸がある場合の粗さ $W(\ell, t)$ に対する閉じた形式の式を導出しました。

• 結果式: 粗さは、コヒーレントなバルスチック（bessel 関数）構造と、散逸による指数関数的減衰の積として記述されます。
  $$\kappa_2(\ell, t) = 2\bar{n}(1-\bar{n}) \left[ \ell - e^{-\Gamma t} \sum_{r=-(\ell-1)}^{\ell-1} (\ell-|r|) J_r^2(Jt) \right]$$
  ここで $\Gamma = \gamma_l + \gamma_p$ は全散逸率です。
• 二パラメータ FV スケーリング: 散逸がない場合（$\Gamma \to 0$）は、従来の一パラメータ FV スケーリング（変数 $t/\ell^z$）に従いますが、散逸がある場合は二パラメータ・スケーリングが現れます。
  $$W(\ell, t) \sim \ell^\alpha f\left( \frac{t}{\ell^z}, \frac{t}{t_\Gamma} \right)$$
  ここで、$t_\Gamma \sim \Gamma^{-1}$ は散逸の緩和時間尺度です。
• 振る舞いの遷移:
  • 弱い散逸 ($t_\Gamma \gg t^*$): 系はまずコヒーレントなバルスチック成長（$z=1, \alpha=\beta=1/2$）を示し、その後有限サイズ効果で飽和します。
  • 強い散逸 ($t_\Gamma \ll t^*$): 散逸が支配的となり、コヒーレントな有限サイズ時間 $t^* \sim \ell/J$ に到達する前に、散逸時間 $t_\Gamma$ で粗さが飽和します。この領域では、従来の変数 $t/\ell^z$ によるスケーリング・カプセル化（collapse）は破綻し、散逸変数 $\Gamma t$ による支配が観測されます。

B. 相互作用系（$\Delta \neq 0$）における数値的発見

テンソルネットワークシミュレーションを用いて、相互作用がある場合の挙動を調べました。

• ユニタリ進化（散逸なし）:
  • 有限磁化密度（$\zeta \neq 1$）: 異方性 $\Delta$ の値に関わらず、粗さの成長は**バルスチック（$z \approx 1$）**に保たれます。相互作用は非普遍振幅や微小時間領域の振る舞いを変えますが、スケーリングの崩壊は起こらず、準粒子（マグノン）の伝播による有限速度輸送が支配的です。
  • 無限温度（$\zeta = 1$）: 既知の輸送相図に従い、$\Delta < 1$ でバルスチック、$\Delta = 1$ で KPZ 型超拡散（$z \approx 3/2$）、$\Delta > 1$ で拡散的（$z \approx 2$）な振る舞いを示します。
• Lindblad 進化（散逸あり）:
  • 相互作用の有無や初期状態の磁化に関わらず、散逸が支配的となります。
  • 粗さはコヒーレントな輸送変数 $t/\ell^z$ による一パラメータ・スケーリングに従わず、散逸時間尺度 $t_\Gamma$ によって支配される飽和を示します。
  • 散逸が導入されることで、初期状態のセグメント配置に関する「記憶」が指数関数的に失われるため、コヒーレントな FV 成長が抑制され、散逸による緩和がスケーリングを決定づけます。

C. 可積分性の破れ（次近接相互作用）の影響

可積分性を破る次近接結合（$J_2$）を追加した場合の検討も行いました。

• ユニタリ進化: 有限磁化密度では、可積分性が破れてもバルスチックなスケーリング（$z=1$）は頑健に維持されます。可積分性の破れは、実効的な伝播速度を再規格化し、遷移時刻 $t^*$ をシフトさせますが、スケーリングの崩壊を引き起こしません。
• Lindblad 進化: 可積分性の破れがあっても、散逸による支配的な振る舞い（$t_\Gamma$ による飽和）は変化しません。散逸がコヒーレントな記憶喪失を支配するため、可積分性の有無は Lindblad 進化のスケーリング構造には本質的な影響を与えません。

4. 意義 (Significance)

• 開いた量子系における FV スケーリングの拡張: 本研究は、FV スケーリングが散逸を含む開いた量子系においても有効であることを示し、従来の一パラメータ・スケーリングから、散逸時間尺度を含む二パラメータ・スケーリングへと一般化しました。
• 輸送普遍性クラスと散逸の競合: 散逸が輸送の普遍性クラス（バルスチック、拡散、KPZ など）をどのように変調するかを明らかにしました。特に、強い散逸下では、ハミルトニアンの対称性や可積分性に関わらず、散逸時間尺度がダイナミクスを支配し、コヒーレントな輸送特性が「隠蔽」されることを示しました。
• 実験への示唆: 冷原子ガスや量子シミュレータなど、散逸を制御可能な量子系において、磁化揺らぎの測定を通じて、散逸とコヒーレントなダイナミクスの競合を直接観測する道筋を提供します。

要約すると、この論文は、散逸性 XXZ 鎖における磁化揺らぎのダイナミクスを解析的・数値的に解明し、**「散逸が支配的な場合、コヒーレントな輸送スケーリングは破れ、散逸時間尺度による二パラメータ・スケーリングが現れる」**という重要な結論を導き出しました。