[論文レビュー] Aging dynamics in interacting many-body systems

Lloyd P. Sanders, Michael A. Lomholt+4 more|arXiv (Cornell University)|Nov 15, 2013

Theoretical and Computational Physics参考文献 76被引用数 27

ひとこと要約

本稿は、ハードコア的排除とべき則分布の待ち時間を持つ一次元で相互作用する多粒子系における超遅い老化ダイナミクスを調査する。スケーリング理論とシミュレーションを用いて、待ち時間がスケールフリーのべき則分布に従う場合（0 < α < 1）に、ラベル付き粒子の平均二乗変位 ⟨x²(t)⟩ が対数的増加 (log t)^{1/2} に比例することを示し、固定された近傍粒子による持続的ブロッキングに起因する極めて遅い拡散を明らかにする。

ABSTRACT

Low-dimensional, complex systems are often characterized by logarithmically slow dynamics. We study the generic motion of a labeled particle in an ensemble of identical diffusing particles with hardcore interactions in a strongly disordered, one-dimensional environment. Each particle in this single file is trapped for a random waiting time $τ$ with power law distribution $ψ(τ)\simeqτ^{-1- α}$, such that the $τ$ values are independent, local quantities for all particles. From scaling arguments and simulations, we find that for the scale-free waiting time case $02$ we recover Harris law $\simeq t^{1/2}$.

研究の動機と目的

• 排他体積相互作用を有する強い無秩序系における、老化と長尾待ち時間のトレーサー拡散に与える影響を理解すること。
• 単一ファイル運動とべき則待ち時間分布を有する連続時間ランダムウォーク（CTRW）ダイナミクスの相乗的相互作用を調査すること。
• 平均二乗変位 ⟨x²(t)⟩ ∝ t^γ の動的指数 γ が待ち時間指数 α にどのように依存するかを特定すること。
• 標準的ハリススケーリング（⟨x²(t)⟩ ∝ t^{1/2}）が回復されるか、あるいは超遅い対数的ダイナミクスに置き換えられる条件を明確にすること。

提案手法

• 一次元格子上にN個の同一粒子からなる単一ファイルをモデル化し、ハードコア的排除を課す。各粒子はべき則待ち時間分布 ψ(τ) ∝ τ^{−1−α} に従うCTRWを実行する。
• 確率的更新ルールを用いる：各ステップで、最も早く待ち時間が満了する粒子を選択し、0 < α < 1 の範囲でスケールフリーな分布から τ を抽出する。
• 待ち時間 τ を τ = τ* [r^{−1/α} − 1] の形で生成するCTRW更新機構を実装する。ここで r ∈ [0,1] は一様乱数であり、べき則統計を保証する。
• 長時間にわたるシミュレーションを実行し、ラベル付きトレーサー粒子の平均二乗変位（MSD）を計算してスケーリング指数を抽出する。
• スケーリング議論を適用し、異なる α の領域におけるMSD挙動の解析的予測を導出する。
• 既知のモデル（標準的単一ファイル：ハリス、バイアス付きCTRW、フェイク粒子のチェンジングモデル）と比較し、排他体積の役割を特定する。

実験結果

リサーチクエスチョン

• RQ1 0 < α < 1 の範囲でべき則分布の待ち時間が単一ファイル系におけるトレーサー粒子の平均二乗変位に与える影響は何か？
• RQ2 排他体積相互作用は、長尾待ち時間を持つ系における超遅いダイナミクスの出現にどのような影響を及ぼすか？
• RQ3 α のどの値の範囲で、標準的ハリススケーリング（⟨x²(t)⟩ ∝ t^{1/2}）が回復されるか？
• RQ4 中間領域 1 < α < 2 において、⟨x²(t)⟩ ∝ t^γ の動的指数 γ は α にどのように依存するか？
• RQ5 MSDにおける対数的スケーリング (log t)^{1/2} の起源は何か？ 他の対数的拡散モデルとはどのように異なるか？

主な発見

• 0 < α < 1 の場合、トレーサー粒子の平均二乗変位は ⟨x²(t)⟩ ∝ (log t)^{1/2} に比例し、長期間にわたる近傍粒子による固定的ブロッキングに起因する超遅い対数的拡散を示す。
• 1 < α < 2 の領域では、平均二乗変位は ⟨x²(t)⟩ ∝ t^γ に比例し、γ < 1/2 となる。これは有限な平均待ち時間を持つにもかかわらず、異常な準拡散を示す。
• α > 2 の場合、標準的ハリススケーリング ⟨x²(t)⟩ ∝ t^{1/2} が回復され、系は通常の単一ファイル系と同様に振る舞う。
• 超遅い対数的スケーリングは、単一粒子の捕獲（シナイ拡散のものとは異なり）ではなく、排他体積効果に起因する強い多体相関に起因する。
• 長時間におけるシステムサイズに依存しないスケーリング挙動が、N = 201 個の粒子と L = 600 の格子サイズを用いた広範なシミュレーションで確認された。
• フェイク粒子のチェンジングモデル（⟨x²(t)⟩ ∝ t^{α/2}） やクラスタリングCTRWモデルとは対照的に、排他体積の役割が超遅いダイナミクスを誘発することの重要性が浮き彫りになった。