Criticality without Temperature in an Ising Spin System

技術要約：イジング・スピン系における温度に依存しない臨界性

問題提起
従来のイジングモデルにおける臨界性は、ハミルトニアンの動力学によって生成され、温度によって制御される。そこでは熱ゆらぎが自発的対称性の破れを駆動する。非平衡システムにおいても、動的な規則（例：吸収状態転移）を通じて秩序が生じることがあるが、イジング系がそのハミルトニアンから完全にデカップリングされ、温度のようなパラメータも持たない状態で、いかにして臨界的な挙動を示すことができるのかという具体的なメカニズムは、未解決の問いである。本論文は、記憶依存的な連結性のみによって、基礎となるエネルギー動力学なしに、イジング・スピン系において集団的な臨界挙動を生成できるか否かについて検討する。

手法：静止時間クラスターモデル (Quiet-Time Cluster Model: QTCM)
著者らは、正方格子（サイズ $N \times N$）上のイジング・スピン ($s_i = \pm 1$) に対する、新しい適応型クラスター動力学を提案している。この核心的な革新は、ダイナミクスをイジング・ハミルトニアン ($H = -\sum_{\langle ij \rangle} s_i s_j$) から切り離している点にある。温度依存的な結合確率を用いる代わりに、システムは「静止時間」変数 $q(i, t)$ に基づいて進化する。ここで $q(i, t)$ は、スピン $i$ が最後に更新されてからの経過時間として定義される。

ダイナミクスは以下のように進行する：

1. クラスター構築: クラスターは標準的なクラスター・アルゴリズムと同様に構築されるが、連結性はスピンの静止時間によって支配される。クラスター内のスピン $s_i$ の隣接スピン $s_j$ がクラスターに加えられるのは、$s_j = s_i$ である場合のみである。このスピンが時刻 $t$ において追加される確率は以下の通りである：
   $$P_{con}(j, t) = 1 - \exp\left(-\alpha \frac{q(j, t)}{\langle q(t) \rangle}\right)$$
   ここで $\alpha$ は連結強度を制御する制御パラメータであり、$\langle q(t) \rangle$ はグローバルな平均静止時間である。
2. メモリ更新: クラスターが形成された後、その中のすべてのスピンは、エネルギー的な受理基準なしに反転する。これらのスピンの静止時間はゼロにリセットされ ($q(i, t+1) = 0$ for $i \in C(t)$)、それ以外のスピンの静止時間は1ずつ増加する ($q(i, t+1) = q(i, t) + 1$)。
3. 観測量: イジング・ハミルトニアンは、得られる構成を特徴付けるための観測量としてのみ扱われる。著者らは、磁化 ($m$)、平均エネルギー ($E$)、クラスターサイズ分布 ($n(s)$)、および平均連結確率 ($P$) を含むグローバルな観測量を監視する。

シミュレーションは、$32 \times 32$ から $256 \times 256$ の格子サイズを用い、$10^6$ から $4 \times 10^6$ のアルゴリズム・ステップにわたって、$\alpha$ を $[0.5, 1.2]$ の範囲で変化させて行われた。

主要な結果
連結パラメータ $\alpha$ の変化は、システムを明確に異なる動的レジームへと導き、最終的に臨界的な挙タビズムへと到達する：

• 自発的対称性の破れと磁化のダイナミクス:

  • 低 $\alpha$（例：0.5）では、磁化はゼロの周囲で極めて小さな振幅で変動する。
  • $\alpha$ が増加するにつれ、ゆらぎが増大する。$\alpha \approx 1.0$ 付近では、熱的な臨界点にある強磁性体と同様のダイナミクスを示し、磁化分布は平坦化する。
  • $\alpha \approx 1.12$ において、システムは磁化分布が $m = \pm 1$ 付近に2つの支配的な対称ピークを持つ転移を起こす。これは自発的対称性の破れを示唆している。しかし、平衡系とは異なり、システムは単一の状態に落ち着くのではなく、正負の値の間を連続的にジャンプし続ける。
  • $\alpha > 1.12$（例：1.2）では、磁化のタイムシリーズはほぼ周期的になり、$m \approx -1$ と $m \approx 1$ の間を急速にジャンプする。

• 特異なエネルギー応答:
  ハミルトニアンはダイナミクスを駆動していないにもかかわらず、スピンあたりの平均エネルギー $\langle E \rangle$ は $\alpha$ に対して特異な応答を示す。$-\partial \langle E \rangle / \partial \alpha$ は、$\alpha_s \approx 1.12$ 付近で $(\alpha - \alpha_s)^{-0.93}$ という形の片側べき乗則の発散を示す。

• クラスターサイズ分布:

  • 小さな $\alpha$ では、クラスターサイズに指数関数的な減衰が見られる。
  • $\alpha$ が臨界領域 $[1, 1.12]$ に近づくと、分布は拡張されたべき乗則の領域 ($n(s) \sim s^{-\tau}$) を発達させる。ここで指数は $\tau \approx -1.08$ から $-1.11$ である。
  • このべき乗則の挙動は、クラスターサイズの数桁にわたって広がり、システムサイズ ($N^2$) と共にスケールし、最大スケールにおける小さな有限サイズ・ハンプ（隆起）によって終端する。
  • $\alpha > 1.12$ では、べき乗則の領域は消失し、分布は指数関数的な減衰の後に非常に大きなサイズでの鋭いハンプを持つ、不連続なものとなる。

• 連結性とパーコレーション:
  平均連結確率 $P$ は $\alpha$ と共に進化する。低 $\alpha$ では $P$ は一定である。$\alpha$ が増加するにつれ、$P$ の分布は広がる。臨界点 ($\alpha \approx 1.12$) 付近では、分布は正方格子のサイト・パーコレーション閾値 ($P \approx 0.59$) 付近にピークを持つ。$\alpha > 1.12$ では、臨界パーコレーション領域内にある $P$ の統計的重みは消失し、これがべき乗則クラスター分布の消失と相関している。

意義と主張
本論文は、イジング・スピン系において、ハミルトニアンのダイナミクスや熱力学的な制御パラメータなしに、臨界的な集団挙動が生成され得ることを示すと主張している。特徴的なメカニズムは、記憶依存的な連結性ルールである。すなわち、個々のスピンの時間的持続性（静止時間）が、マクロなスケールでシステムを組織化するフィードバックループとして機能する。

著者らは、この構成が従来のイジング・ハミルトニアンの役割を逆転させていると断言している。すなわち、ハミルトニアンはダイナミクスを生成するのではなく、適応プロセスによって生成された構成を特徴付ける観測量として機能している。これらの結果は、適応的な時間的記憶のみによって集団的な臨界性が生成され得ることを示唆しており、このような複雑なシステムにおける臨界挙動は必ずしも有効なエネルギー景観に根ざしていなければならないという仮定に疑問を投げかけている。論文では、この転移はパーコレーションの性質（特に連結確率がパーコレーション閾値に近いこと）を共有しているものの、特定のパーコレーション・クラスとの普遍性を確立するには、さらなる有限サイズ・スケーリング解析が必要であり、それは今後の課題であると述べている。