Emergence of continuously varying critical exponents in coupled map lattice as an effect of quenched disorder

技術要約：結合マップ格子における連続的に変化する臨界指数の出現：クエンチされた無秩序の影響として

問題提起
本論文は、非平衡吸収相転移（特に指向性パーコレーション（DP）クラス）における普遍性クラスの堅牢性を調査している。普遍性クラスは通常、次元と対称性によって決定される固定された臨界指数によって定義されるが、著者らはクエンチされた空間的無秩序がこれらのクラスにどのように影響するかを探求している。先行研究では、無秩序が活性化スケーリング（対数減衰）や新しい普遍性クラスをもたらすことが示されているが、1次元（1D）非平衡系において、臨界点での「連続的に変化する」臨界指数が生じる具体的なシナリオは依然として関心の対象となっている。著者らは、既知のDPモデルに非対称なクエンチ無秩序を導入することで、臨界指数が無秩序の強さに応じて連続的に変化するレジームを誘起できるかどうかを明らかにすることを目的としている。

手法
著者らは、JanakiとSinhaによって研究された、層流（吸収）状態から乱流（活性）状態への転移を示す1次元結合マップ格子（CML）モデルを利用している。

• 基本モデル: システムは、$N$個のサイト（$N = 2 \times 10^5$）と変数 $x_i(t)$ で構成される。進化は、固定パラメータ $k=3.2, \omega=0.29$ を持つ円写像 $f(x) = x + \omega - \frac{k}{2\pi}\sin(2\pi x)$ によって制御される。秩序パラメータは、乱流サイト（$x_i > 0.5$）の割合 $m(t)$ である。
• 無秩序の実装: 非対称な結合を通じて空間的なクエンチ無秩序が導入される。結合パラメータ $\epsilon$ は、$\epsilon_1 = \epsilon/2 + 0.1$ と $\epsilon_2 = \epsilon/2 - 0.1$ の2つの成分に分割される。サイトは、「タイプA」（確率 $p$）または「タイプB」（確率 $1-p$）としてランダムに割り当てられる。
  • タイプAのサイトは左隣のサイトとより強く結合（$\epsilon_1$）し、右のサイトとは弱く結合（$\epsilon_2$）する。
  • タイプBのサイトは右のサイトとより強く結合（$\epsilon_1$）し、左のサイトとは弱く結合（$\epsilon_2$）する。
  • パラメータ $p$ は0から0.5まで変化させる。$p=0$ は完全な非対称結合（すべての活動が左へ駆動される）を表し、$p=0.5$ は局所的な非対称性にもかかわらず、対称な平均を表す。
• シミュレーション: システムは様々な $p$ 値に対してシミュレートされる。臨界指数は、有限サイズスケーリング、オフクリティカル・スケーリング、および臨界結合 $\epsilon_c$ における秩序パラメータの冪乗則減衰 $m(t) \sim t^{-\delta}$ の解析を通じて決定される。接続行列の固有値スペクトルも、スペクトルの変化と臨界挙動を相関させるために解析される。

主な結果

1. 標準的な普遍性の破壊: $p=0$ の場合、システムは標準的なDP普遍性クラスに属し、臨界指数は $\delta \approx 0.159, \beta \approx 0.253, \nu_{\parallel} \approx 1.732, z \approx 1.58$ となる。しかし、$p > 0$ の場合、システムは連続的に変化する臨界指数を示す。
2. 指数の連続的な変化:
   • 臨界点における減衰指数 $\delta$ は、$p$ とともに連続的に変化する。$p=0.1$ で $\delta = 0.034$ から、$p=0.5$ で $\delta = 0.158$ まで増加する。
   • 定常状態密度 $m(\infty) \sim |\epsilon - \epsilon_c|^\beta$ を支配する指数 $\beta$ も $p$ とともに連続的に増加し、0.125から0.607の範囲をとる。
   • 相関長指数 $\nu_{\parallel}$ および動的指数 $z$ も顕著な変化を示す。特に、$z$ は非常に大きく（$p=0.5$ で $z \approx 3.90$ であり、$p < 0.5$ ではさらに高くなる可能性が高い）、標準的なDPよりも擾乱の移動がはるかに遅いサブディフューシブ（劣拡散）的な挙動を示している。
3. 臨界における冪乗則減衰: 無秩序のあるシステムの中には臨界点で対数減衰（活性化スケーリング）を示すものもあるが、このモデルはすべての $p$ において冪乗則減衰 $m(t) \sim t^{-\delta}$ を維持する。ただし、その指数は無秩序の割合に依存する。
4. 固有値スペクトルとの相関: 著者らは、変化する臨界挙動を接続行列の固有値スペクトルに関連付けている。
   • $p=0$ のとき、スペクトルは複素平面上に楕円を形成する。
   • $0 < p < 0.5$ のとき、スペクトルは「穴」を発達させ、固有値は曲線上に厳密に位置するのではなく、境界の周囲に広がる。
   • $p=0.5$ のとき、局所的な非対称性にもかかわらず、スペクトルは実数（線分上）となる。これは、接続が「平均的に対称」であることを示唆している。
5. 吸収相における非普遍的挙動: 吸収相（$\epsilon_c$ 未満）において、一連の結合パラメータ $\epsilon$ に対して、一連の冪乗則減衰指数が観察される。この非普遍的な挙動の範囲は、$p$ が減少するにつれて縮小する。

意義と主張
本論文は、臨界指数が無秩序の強さに応じて連続的に変化するという、非平衡系における非普遍的な臨界挙動のユニークな事例を実証したと主張している。

• 新しい普遍性クラス: 著者らは、指数が固定されず $p$ の関数であるため、このシステムは既知のどの普遍性クラス（DPや活性化スケーリングクラスを含む）にも適合しないことを示唆している。
• 弱普遍性との区別: 本挙動は、スケーリング関係は保持されるが指数が変化する「弱普遍性」のシナリオとは異なる。ここでは、すべての臨界指数（$\delta, \beta, \nu_{\parallel}, z$）が連続的に変化し、システムは従来の意味での標準的なハイパースケーリング関係に従わない（ただし、著者らは $\nu_{\parallel} \approx \beta/\delta$ が近似的に成立することには言及している）。
• メカニズム: 連続的な変化は、非対称な無秩序によって反対方向へと駆動される活動の間の「葛藤（tussle）」または競争に起因する。この競争が、低速で移動する欠陥と極めて遅い緩和時間（大きな $z$）をもたらし、システムが標準的な固定点に落ち着くのを妨げている。
• 理論的含意: 本研究は、特定の種類のクエンチ無秩序が存在する場合、固定された普遍性クラスという概念が崩壊し、固定点の列や連続的に変化する指数が生じ得ることを浮き彫りにしている。これは以前、平衡統計力学（例：8頂点モデル）においてより一般的であった現象であるが、今や駆動された非平衡CMLにおいて観察された。

著者らは、このモデルが、クエンチされた空間的無秩序がいかにしてシステムの臨界力学を根本的に変え、臨界指数が普遍的な定数ではなく無秩序分布の具体的な詳細に依存するレジームをもたらすかを示す明確な例を提供すると結論づけている。