Zigzag ordering, defects, and anomalous relaxation in antiferromagnetic Kuramoto lattices

技術要約：反強磁性的蔵本格子におけるジグザグ秩序、欠陥、および異常緩和

問題提起
本論文は、負の（反強磁性的）最近接結合を持つ結合蔵部振動子の非平衡秩序ダイナミクスを調査している。正の結合を持つ蔵部モデルは同期転移に関してよく研究されているが、反位相の秩序（位相差 $\pi$）を好む斥力的結合のダイナミクスについては、あまり探索されていない。具体的には、著者らは、決定論的な非線形相互作用と幾何学的フラストレーションのみによって、どのようにシステムが秩序状態へと駆動されるのか、またそのプロセスが標準的な普遍性クラス（Isingモデルなど）に従うのか、あるいは異常なスケーリングを示すのかという点に関する理解の欠如に取り組んでいる。

手法
本研究では、数値シミュレーションと解析的な連続体近似の両方を用いている。

• 数値シミュレーション: 著者らは、$K < 0$ の蔵部方程式に従う $N$ 個の振動子の1次元鎖および2次元正方格子をシミュレートしている。運動方程式は、4次ルンゲ・クッタ法を用いて積分されている。
  • 欠陥の定義: 局所的な勾配の符号（$s_i = \text{sign}(\phi_i - \phi_{i-1})$）が交互に変化しない場合（すなわち $|s_i + s_{i+1}| = 2$）、それを欠陥と定義する。秩序パラメータは、正規化された欠陥密度 $D(t)$ である。
  • パーシスタンス（持続性）: 2種類のパーシスタンスを測定している。一つは勾配パーシスタンス（$P_s$）：局所的な勾配が初期の符号を保持する確率。もう一つは場（フィールド）パーシスタンス（$P_\phi$）：位相変数がゼロを横切らない確率である。
  • パラメータ: スケーリングの独立性を検証するために、システムサイズ ($N$)、結合強度 ($K$)、固有振動数 ($\omega$) を変化させてシミュレーションを行っている。
• 解析的手法: 離散方程式のテイラー展開を通じて、連続体極限が導出され、粗視化された偏微分方程式（PDE）が得られる。このPDEに対して線形安定性解析を行い、不安定モードと緩和を支配する主要な項を特定する。

主な結果

1. 一次元ダイナミクス（異常粗大化）

• 欠陥の減衰: 欠陥密度は、初期段階において代数的に $D(t) \sim t^{-1/4}$ と減衰する。これは、標準的な非保存型相秩序形成（例：Model A または Ising 系）で見られる典型的な $t^{-1/2}$ の減衰や、動力学指数 $z=2$ よりも著しく遅い。
• 動力学的スケーリング: 有限サイズスケーリング解析により、動力学指数 $z=2$ が明らかになった。これは、特性緩和時間が $t_c \sim N^2$ とスケールすることを示している。これは拡散的なドメイン成長（$L(t) \sim t^{1/2}$）を示唆している。
• 普遍性クラス: $z=2$ と欠陥減衰指数 $\delta = 1/4$ の共存は、独特の普遍性クラスを定義する。著者らは、$\delta = 1/z$（消滅するランダムウォーカーの場合）という標準的なスケーリング関係を、この組み合わせが満たしていないことを指摘している。
• パーシスタンス: 勾配および場のパーシスタンスは共に、引き伸ばされた指数関数形式 $P(t) \sim \exp(-c t^{1/4})$ に従う。引き伸ばしの指数 $\alpha = 1/4$ は、欠陥減衰指数 $\delta$ と正確に一致している。
• 飽和: 有限系においては、欠陥密度とパーシスタンスは $t_c \sim N^2$ の時間で飽和し、残留したほぼ等間隔な欠陥構成を残す。熱力学的極限では、長距離のジグザグ秩序が回復される。
• 振動数の役割: スケーリング指数は、固有振動数 $\omega$ や結合の大きさ $K$（$K < 0$ である限り）に依存せず、これらは時間の軸を再スケールさせるだけである。ただし、非ゼロの $\omega$ はパーシスタンスの飽和を防ぐ（記憶が最終的に失われる）一方で、$\omega = 0$ の場合は凍結した準安定状態を招く。

2. 二次元ダイナミクス（幾何学的フラストレーション）

• 2次元では、幾何学的フラストレーションが広域的な反強磁性的秩序の形成を阻害する。
• 局所的なチェッカーボードパターンは出現するものの、システムは点欠陥ではなく、拡張されたドメイン壁（欠陥線）を特徴とする準安定状態にトラップされる。
• 欠陥密度は初期の過渡的な減衰を示すが、最終的には非ゼロの値に飽和し、シミュレーションされた時間スケール内での完全な秩序形成を妨げる。これらの拡張された界面を排除するには、動力学的に抑制されている集団的な再配置が必要となる。

理論的説明
著者らは、位相場 $\phi(x,t)$ に対する連続体方程式を導出している：
$$\partial_t \phi = -|K| \partial_x^2 \phi - \frac{|K|}{12} \partial_x^4 \phi$$

• 不安定性: 負の拡散項（$\partial_x^2$）が不安定性を引き起こし、一方で双調和項（$\partial_x^4$）が短波長の波長を安定化させる。この競合がジグザグの波長を選択する。
• スケーリングの起源: $t^{-1/4}$ の減衰は、中間スケールにおいて（表面拡散のような挙動を制御する）双調和項が支配的であることに起因する。しかし、長波長のダイナミクスは、最終的には拡散的な $\partial_x^2$ 項によって支配され、これが全球的な平衡化時間における $z=2$ のスケーリングをもたらす。著者らは、欠陥の運動がジグザグ背景のより遅い平滑化に「隷属」しているため、拡散的な $z=2$ の背景にもかかわらず、異常な $\delta = 1/4$ 指数が現れるのだと主張している。
• パーシスタンスの導出: 運動論的な議論により、パーシスタンス指数と欠陥密度の関係が結び付けられる。欠陥の遭遇率が $\Gamma(t) \sim t^{-3/4}$ （$D(t) \sim t^{-1/4}$ およびドメインサイズ $\xi(t) \sim t^{1/2}$ から導出）とスケールすると仮定すると、累積ハザードは $t^{1/4}$ となり、引き伸ばされた指数関数的なパーシスタンスが得られる。

意義と主張
本論文は、決定論的な非線形ダイナミクスと幾何学的フラストレーションが、クエンチされた無秩序や確率的ノイズを必要とせずに、多体系における遅い緩和と異常なスケーリングを生み出すのに十分であることを主張している。

• 標準的なスケーリングの崩壊: 本研究は、反強磁性的蔵部格子における標準的な動力学的スケーリング関係の崩壊を示し、$\delta = 1/4$ および $z=2$ によって特徴付けられる独自の普遍性クラスを明らかにしている。
• 遅い緩和のメカニズム: 緩やかな緩和は、外部のランダム性からではなく、進化するジグザグ背景内における欠陥の制約された移動性に起因することを確立している。
• 次元依存性: 本研究は、次元がいかに根本的に結果を変えるかを強調している。すなわち、1次元系は異常なキネティクスを伴いつつ最終的に長距離秩序に到達するが、2次元系はフラストレーションを受けた、持続的なドメイン壁を持つ準安定状態に捕らえられる。
• 理論的洞察: 連続体近似は、四次の振幅減衰と二次的な時間スケーリングの共存をうまく説明しており、観察された指数に対する理論的枠組みを提供している。

著者らは、これらの結果が連続的な位相系における粗大化ダイナミクスに新しい視点を提供するものであり、フラストレーションそれ自体が、標準的なスピンモデルとは異なる複雑で遅い緩和現象を引き起こし得ることを示していると結論づけている。