Nonequilibrium Phase Transition To Temporal Oscillations In Mean-Field Spin Models

技術要約：平均場スピンモデルにおける時間振動への非平衡相転移

問題提起
相互作用する単位系の集合における自発的な振動の出現は、生化学的時計からアクティブマターに至るまで、非平衡系において観察される顕著な特徴である。熱力学的極限におけるこのような振動の開始は、通常、決定論的なホップ分岐によって記述されるが、メゾスコピックな系においては、ゆらぎが重要な意味を持つため、確率論的な取り扱いが必要となる。非平衡相転移に関する既存の理論的枠組みは、エントロピー生成を一般化された熱力学的ポテンシャルとして特定したり、磁化のみを用いてランダウ理論を拡張したりすることに依存することが多い。しかし、これらの手法は、時間並進対称性の破れを伴う時間振動の具体的な開始を特徴付けることが困難であり、磁化とその時間微分（の確率的な値）の両方のダイナミクスを捉える一般化されたランダウ自由エネルギーを提供できていない。

手法
著者らは、$N$個のスピン（$s_i = \pm 1$）および潜在的に補助変数を持つ、駆動されたキネティック・スピンモデルのための平均場理論を提案している。この手法の核心は以下の通りである：

1. 確率的微分の定義： 磁化 $m(t)$ の直接の時間微分に伴う発散するホワイトノイズのゆらぎを避けるため、著者らは、マイクロスコピックな配置 $C$ に対する平滑化された確率的時間微分 $\dot{m}(C)$ を定義する。これはマルコフ・ジャンプ過程の遷移速度 $W(C'|C)$ から導かれる：
   $$\dot{m}(C) = \sum_{C' \neq C} (m(C') - m(C)) W(C'|C)$$
   この定義により、$\langle \dot{m} \rangle = d\langle m \rangle / dt$ が保証され、ゆらぎは $m$ と同程度にスケールする。

2. 同時分布と大偏差理論： 理論は、磁化 $m$ とその確率的微分 $\dot{m}$ の同時確率分布 $P_N(m, \dot{m})$ に焦点を当てる。大きな $N$ の極限において、この分布は、大偏差形式 $P_N(m, \dot{m}) \sim \exp[-N\phi(m, \dot{m})]$ に従うと仮定される。ここで $\phi$ はレート関数（非平衡ランダウ自由エネルギーと解釈される）である。

3. マスター方程式の展開： $P_N(m, \dot{m})$ に関する粗視化されたマスター方程式から出発し、著者らはレート関数 $\phi(m, \dot{m})$ に関する定常状態方程式を導出する。時間並進対称性が破れる臨界点付近でこの方程式を展開することにより、制御パラメータ $u_0$（$\dot{m}$ に関するドリフト項の微分に関連するもの）を特定する。

4. ハミルトニアン構造： 転移の近傍において、著者らはレート関数 $\phi(m, \dot{m})$ が有効ハミルトニアン $H(m, \dot{m}) = \frac{1}{2}\dot{m}^2 + V(m)$ の関数 $f(H)$ として表現できることを示す。関数 $f(H)$ はランダウ自由エネルギーの役割を果たすであり、その最小値がマクロな状態を決定する。

主要な貢献と結果

• 一般化されたランダウ自由エネルギー： 著者らは、磁化と、その確率的時間微分に依存する非平衡一般化ランダウ自由エネルギー $\phi(m, \dot{m})$ を構築した。これは、磁化のみに依存していた従来の一般化とは一線を画すものである。
• 秩序パラメータ： 振動状態への転移の秩序パラメータは、有効ハミルトニアン $H$ であると特定される。
  • 時間に依存しない相（常磁性または強磁性）では、$f(H)$ の最小値は $H^* = 0$ で起こる。
  • 振動相では、最小値は $H^* > 0$ へとシフトし、$(m, \dot{m})$ 平面におけるリミットサイクルの開始を合図する。
• スケーリングと臨界挙動：
  • 連続転移： 一般的なパラメータにおいて、転移は連続的である。リミットサイクルの大きさは $H^* \sim \epsilon$ とスケールする（ここで $\epsilon$ は臨界点からの距離）。観測量は $\langle m^2 \rangle \sim \langle \dot{m}^2 \rangle \sim \epsilon$ とスケールする。
  • 三重点： 常磁性、強磁性、および振動相が合流する三重点付近では、ポテンシャル $V(m)$ が 4 次（$m^4$）となる。この領域では、自由エネルギー $f(H)$ は非解析的な形式（$f(H) \sim -\epsilon H + c H^{3/2}$）をとり、異なるスケーリング則（$H^* \sim \epsilon^2$、$\langle m^2 \rangle \sim \epsilon$、$\langle \dot{m}^2 \rangle \sim \epsilon^2$）をもたらす。リミットサイクルは平坦化し、周期は $\tau \sim \epsilon^{-1/2}$ として発散する。
  • 不連続転移： 理論は、常磁性状態と振動状態が $f(H)$ の局所的最小値として共存する不連続転移も考慮に入れている。この場合、観測される相は $f(H)$ の大域的最小値によって決定される。
• オーバーラップ分布とレプリカ対称性の破れ： 著者らは、2つの独立したスピン配置間のオーバーラップ分布 $P(q)$ を評価する。振動相において、$P(q)$ は $q=0$ で対数的な発散を伴う連続的な支持（サポート）を示す。この挙動は、無秩序系で見られる連続的なレプリカ対称性の破れに類似しているが、現在のモデルにはクエンチされた無秩序は存在しない。
• 明示的なモデル： 理論は、2つの変数群（スピンと場）を持ち、詳細釣合が破れたパラメータ $\mu$ によって制御されるキネティック平均場イジングモデルを用いて示されている。モンテカルロ・シミュレーションは、常磁性、強磁性、および振動相を含む相図、ならびに共存領域と三重点を裏付けている。

意義と主張
本論文は、平均場スピンモデルにおける自発的振動への非平衡相転移を記述するための一貫した理論的枠組みを提供すると主張している。秩序パラメータの確率的時間微分をランダウ理論に組み込むことで、著者らは振動の開始を特徴付けるハミルトニアンに基づく秩序パラメータを特定することに成功した。

本研究は、振動相が単なる動力学的不安定性ではなく、特定の観測量の結合分布によって特徴付けられる熱力学的相であることを強調している。さらに、無秩序が存在しないにもかかわらず、無秩序系に見られるものと同様の非自明なオーバーラップ分布を発見したことは、駆動された非平衡振動系と無秩序平衡系との間の深い構造的類似性を示唆している。著者らは、彼らの摂動論的枠組みが、臨界温度からの小さな偏差および制御パラメータ $\mu$ の特定の範囲において有効であることを述べており、有限次元系におけるこれらの転移を特徴付けるための将来の繰り込み群を用いた研究の必要性を指摘している。