Magnetic domains stabilized by symmetry-protected zero modes

この論文「Magnetic domains stabilized by symmetry-protected zero modes（対称性保護されたゼロモードによって安定化された磁気ドメイン）」の技術的サマリーを以下に日本語で提供します。

1. 研究の背景と問題設定

閉じた量子多体系における熱化の破れ（非エルゴード性）のメカニズムの理解は、量子多体物理学の中心的な課題です。従来の熱化の破れを説明する主要なメカニズムには、積分可能性、多体局在（MBL）、ヒルベルト空間の断片化（HSF）、量子多体傷（QMBS）の 4 つが知られていますが、これらとは質的に異なる新しいメカニズムの存在が示唆されてきました。

本研究は、XX モデル（スピン 1/2 粒子の二次元正方格子）において、対称性によって保護された「ゼロモード（エネルギー E=0 の固有状態）」の存在が、磁気ドメインの無限寿命化（非エルゴードな振る舞い）を引き起こすことを発見しました。特に、2 本の鎖が結合したスピンラダー系において、ドメインウォール（磁化の境界）初期状態が、長時間にわたり均一でない磁化プロファイルを維持する現象を解明しました。

2. 手法とモデル

• モデル: $N_x \times N_y$ の正方格子上のスピン 1/2 粒子系を記述する XX モデル（式 1）。
  $$\hat{H} = \sum_{\lambda=\parallel,\perp} \sum_{\langle i,j\rangle_\lambda} J_\lambda \left[ \hat{S}^x_i \hat{S}^x_j + \hat{S}^y_i \hat{S}^y_j \right]$$
  ここで、$J_\parallel$ と $J_\perp$ はそれぞれ鎖方向（$x$ 方向）と鎖間結合（$y$ 方向）の結合定数です。
• 初期状態: 完全なドメインウォール状態 $|\Psi(0)\rangle = |\Uparrow \cdots \Uparrow \Downarrow \cdots \Downarrow\rangle$（左半分がスピンアップ、右半分がスピンダウン）。
• 解析手法:
  1. ランチョス法（Lanczos algorithm）: 初期状態からクリロフ基底を生成し、有効的な 1 次元緊密結合モデルを導出。これにより、ゼロモードと初期状態の重なり（オーバーラップ）を熱力学極限で評価。
  2. 対角化（Exact Diagonalization）: 小規模系（最大 $8 \times 2$）での完全対角化を行い、エンタングルメントエントロピーやスペクトル構造を解析。
  3. 対称性の解析: キラル対称性（Chiral symmetry）の存在とその破れが熱化に与える影響を調査。

3. 主要な発見と結果

A. 非エルゴードな磁気ドメインの安定化

• 鎖が結合していない場合（$J_\perp = 0$）、Jordan-Wigner 変換により自由フェルミオン系となり、スピンはバリスティックに移動し、時間平均磁化はゼロになります（熱化）。
• 鎖間結合 $J_\perp$ を導入すると、臨界結合 $J_c^\perp \approx 0.5 J_\parallel$ を超えた領域で、磁化プロファイルが均一化せず、ドメインウォールが長時間安定して残存します（図 1）。
• この現象は、初期状態が各ラング（段）で強磁性相関（$C_{ZZ}=1$）を持つ場合にのみ観測され、反強磁性欠陥を含む初期状態では熱化が回復します。

B. 局在転移とランチョス係数の振動

• ランチョス法を用いた解析により、クリロフ基底展開における係数 $\beta_j$ の振動が局在の鍵であることが判明しました。
• $J_\perp > J_c^\perp$ の場合、ゼロエネルギー固有状態の展開係数 $|c_j|^2$ が $j^{-\gamma}$ ($\gamma > 1$) より速く減衰し、ゼロモードが初期状態と有限の重なりを持ちます（局在相）。
• 逆に $J_\perp \le J_c^\perp$ では減衰が遅く（$\gamma \le 1$）、オーバーラップは熱力学極限でゼロになります（非局在相）。
• この振動は、ラング強磁性状態（スピンが同じ）ではハミルトニアンの作用が抑制され、移動度が低下する「構造化された移動度（structured mobility）」に起因します。

C. 対称性保護されたゼロモードの多様性

• XX モデルはキラル対称性 $\hat{C} = \hat{X}\hat{I}\hat{S}$ を持ちます（$\hat{X}$: スピン反転、$\hat{I}$: 空間反転、$\hat{S}$: サブラattice 演算子）。
• この対称性により、ゼロ磁化セクターにおけるゼロモード（$E=0$）の数は、格子サイズ $N$ に対して指数関数的に多くなります（下限：$2^{N/2}$）。
• この指数関数的に多いゼロモードの縮退が、ラング強磁性状態の熱化を妨げる主要因です。初期状態はこれらのゼロモード部分空間に有限の重なりを持っており、その結果として磁気ドメインが安定化されます。
• 格子の縦方向の長さ $N_y$ が偶数の場合にのみこの効果が顕著に現れ、$N_y$ が奇数では消失します。

D. 摂動に対する頑健性と熱化の回復

• 対称性を保存する摂動: キラル対称性を保つ摂動（例：次近接結合など）を加えても、非エルゴードなドメインは頑健に維持されます。
• 対称性を破る摂動:
  • ZZ 相互作用（$\hat{S}^z_i \hat{S}^z_j$）や、キラル対称性を破る結合を加えると、ゼロモードの縮退は解けます。
  • しかし、完全な XXZ モデル（ZZ 結合あり）では、ゼロモードの代わりに「量子多体傷（Quantum Many-Body Scars）」に似た低エンタングルメント状態が残り、磁化は非常に長い時間スケールで振動（リバイバル）を示します。
  • 最終的に、ZZ 結合が十分に強い場合、ドメインは熱化し、磁化はゼロに緩和します。

4. 意義と結論

• 理論的意義: 積分可能性や MBL などの既存の枠組みとは異なり、「対称性保護されたゼロモードの指数関数的な縮退」が、熱力学的に安定な非エルゴードダイナミクスを生み出す新しいメカニズムを確立しました。
• 実験的実現性: このモデルは、光格子中の超低温原子、超伝導量子プロセッサ、リドバーグ原子アレイ、双極性分子など、現在の量子シミュレーターで容易に実装可能です。
• 将来的展望: ランチョス法に基づくこのアプローチは、熱力学極限における局在現象の安定性を評価するための強力な診断ツールとなり得ます。また、キラル対称性による非エルゴード性の安定化は、他の多くの量子多体系にも普遍的に適用可能な原理である可能性があります。

要約すると、この論文は、XX モデルにおける特定の対称性とゼロモードの多様性が、どのようにして磁気ドメインの「凍結」を引き起こし、熱化を回避するかを、ランチョス法と対称性の観点から詳細に解明した画期的な研究です。