Dissipative quantum North-East-Center model: steady-state phase diagram, universality and nonergodic dynamics

技術要約：散逸的量子北東中心（NEC）モデル

問題提起
本研究は、二次元スピン1/2格子系である散逸的量子北東中心（NEC）モデルの非平衡物理を調査するものである。このモデルは、カイラルで動力学的に制約された散逸と、コヒーレントな量子相互作用を組み合わせたものである。主な目的は、様々なハミルトニアンの下でのこの系の定常状態相図を決定し、その双安定相の性質を特徴付けることである。具体的には、古典的なNECモデルで見られる堅牢な双安定性が、量子ゆらぎが存在する場合でも持続するかどうか、また、異なる微視的な量子ノイズの源（制約あり vs 制約なし）が双安定領域の広がりにどのように影響するか、そして系の動的特性、特に少数スピンアイランドの非エルゴード的挙動について論じている。

手法
著者らは、系の開いた量子ダイナミクスを支配するリンドブラッド・マスター方程式を解くために、**クラスター平均場（CMF）**法を採用している。

• クラスター平均場（CMF）： 短距離相関を系統的に組み込むため、格子をサイズ $\ell \times \ell$ の正方形クラスターに分割する。特定のクラスターのダイナミクスは、隣接するクラスターから導出された境界項を平均場レベルで扱う、CMFハミルトニアンおよび散逸子によって支配される。著者らは、クラスターサイズ $\ell=2$ と $\ell=3$ の結果を比較することで収束性を検証している。
• 不均一CMF（iCMF）： 非並進対称な初期条件を持つ大規模な系の実時間ダイナミクスを研究するために、著者らはiCMFアンザッツを用いた不均一な一般化CMFを利用している。これにより、反対方向の向きを持つ背景に囲まれたスピンの「アイランド」のシミュレーションが可能となる。
• 安定性解析： 均質な定常状態の周囲に小さな空間ゆらぎ（平面波）を導入することにより、線形安定性解析が行われる。得られるスーパーオペレータの固有値は、相の安定性と、臨界点付近における不安定性の性質を決定する。
• 調査対象となったモデル： 本研究では、異なる量子ゆらぎの源を表す3つの異なるクラスのハミルトニアンを比較している：
  1. 制約なし： 均質な横磁場 ($H_X = \Omega \sum \sigma^x_j$)。
  2. 制約あり（2D PXP）： 全4つの隣接スピンに対してリドバーグ・ブロッケード制約を持つPXPモデル。
  3. NEC対称（カイラルPXP）： 制約がカイラルなNEC幾何学に従う（北および東の隣接スピンにのみ作用する）PXPモデル。
     さらに、制約のない散逸プロセス（非コヒーレントなスピン反転）に対する安定性についてもテストされている。

主要な結果

1. 定常状態相図：

   • 系は2つの異なる相を示す：一意の定常状態を持つノーマル相と、反対の巨視的磁化を持つ2つの定常状態によって特徴付けられる双安定相である。
   • 相図は、古典的なゆらぎの振幅（$T$）、そのバイアス（$h$）、および量子的なゆらぎの振幅（$\Omega$）によって定義される。
   • 普遍性： 双安定相の存在はNECモデルの普遍的な特徴であり、検討されたすべてのハミルトニアンにおいて持続する。しかし、双安定領域の広がりは、量子ゆらぎの微視的な詳細に依存する。
   • 制約の影響： 動学的な制約を持つハミルトニアン（PXPモデル）は、制約のないハミルトニアン（$H_X$）よりも双安定相をより効果的に保持する。制約なしの横磁場は、制約ありのモデルと比較して、より低いゆらぎ振幅で双安定性を破壊する。
   • 相転移： 双安定相からノーマル相への転移は、バイアスがゼロ（$h=0$）のときは連続的であるが、$h \neq 0$ のときは一次転移である。臨界線はバイアスに対して二次的な依存性、$T_c(h) \propto (1-|h|)^2$ を従う。

2. 安定性と臨界性：

   • 線形安定性解析により、臨界点付近で有限の波長ベクトルにおける不安定性が生じることが明らかになった。
   • これらの不安定性の構造は、駆動メカニズムによって異なる。量子ゆらぎは有限の波長ベクトル（$|k_x| = \pi/4$）で不安定性のピークを誘起するが、熱的ゆらぎは $k_x=0$ でピークを示す。

3. 非エルゴード的ダイナミクスとアイランドの再吸収：

   • ノーマル相では、少数派のアイランドは背景と混合して消失し、一意の定常状態へと至る。
   • 双安定相では、少数派のアイランドはアイランドのサイズに関わらず、常に周囲の多数派相によって再吸収される。これは、双安定領域における非エルゴード的な性質を裏付けており、系が初期条件（磁化の向き）を保持していることを示している。
   • 再吸収速度： 再吸収プロセスは、アイランドのサイズに依存しない一定の速度で行われる（緩和時間 $\tau$ はアイランドサイズ $\ell_\downarrow$ に対して線形にスケールする）。
   • 運動方程式： 著者らは、再吸収速度 $v$ に関する現象論的な方程式を提案している：
     $$\frac{dr}{dt} = -C + Bh + DT + E\Omega$$
     ここで、$C$ は古典的な吸収率、$B$ はバイアスに関連し、$D, E$ は熱的（$T$）および量子的（$\Omega$）なゆらぎによる速度の線形減少を定量化する。解析の結果、量子的なゆらぎは熱的なゆらぎ（$D \approx 2.4$）と比較して、より弱い影響（$E \approx 0.25$）を与えることが示唆されている。

意義と主張
本論文は、散逸的量子NECモデルが、量子多体系における真の双安定性を実現するための堅牢なメカニズムを提供することを主張しており、量子ゆらぎが必然的にそのような相を破壊するという概念に異を唱えている。知見は以下の点を強調している：

• 堅牢性： 双安定性は、様々な微視的なコヒーレント過程、さらには追加の制約なしの散逸チャネルに対しても安定しており、完全な隔離が不可能な実験的実装（例：リドバーグ原子）への潜在的な関連性を示唆している。
• 非エルゴード性： 本モデルは、エラーアイランドの吸収によって系の運命が決まるという非エルゴード的ダイナミクスを示しており、これは理論的に量子誤り訂正の戦略として機能し得るメカニズムである。
• 方法論的貢献： 本研究は、大規模な二次元散逸格子系、特にドメイン壁の運動のような非並進対称な現象を解析するための、不均一クラスター平均場（iCMF）法の有効性を実証している。

著者らは、カイラルな動力学的制約と散逸の相互作用が、熱平衡の振る舞いとは異なる、非平衡相のユニークなプレイグラウンドを作り出すと結論付けており、高次元への一般化やヒルベルト空間の断片化との関連性を含む将来の方向性を提示している。