[論文レビュー] Self-organized criticality in a discrete model for Smoluchowski's equation

Mathieu Merle, Raoul Normand|arXiv (Cornell University)|Oct 30, 2014

Stochastic processes and statistical mechanics参考文献 24被引用数 7

ひとこと要約

本稿では、クラスターサイズの動的閾値 α(N) を持つ離散的凝集モデルを研究し、N²/³ log⁴(N) ≤ α(N) の条件下で、システムが自己組織的臨界性を示すことを示している。その結果、典型的な有限クラスターの分布は臨界ガルトン＝ワトソン木に収束する。モデルの小クラスターは時刻1以前に亜臨界的なエドーシュ＝レニー確率的グラフを形成し、以降は臨界的な性質を保ち、アルドウスの予想であるクラスター成分の局所的弱収束についての検証がなされる。

ABSTRACT

We study a discrete model of coagulation, involving a large number $N$ of particles. Pairs of particles are given i.i.d exponential clocks with parameter $1/N$. When a clock rings, a link between the corresponding pair of particles is created only if its two ends belong to small clusters, i.e. of size less than $α(N)$, with $1 \ll α(N) \ll N$. The concentrations of clusters of size $m$ in this model are known to converge as $N o \infty$ to the solution to Smoluchowski's equation with a multiplicative kernel. Under the additional assumption $N^{2/3} \log^γ(N) \le α(N)$, for some $γ> 1/3$, we study finer asymptotic properties of this model, namely the combinatorial structure of the graph consisting of small clusters. We prove that this graph is essentially an Erdos-Renyi random graph, which is subcritical before time 1, and remains critical after time 1. In particular, we show that our model exhibits self-organized criticality at a microscopic level: the limiting distribution of a typical finite cluster is that of a critical Galton-Watson tree. Our approach allows in particular to verify, under our additional assumption, a conjecture of Aldous.

研究の動機と目的

• 有限サイズの閾値 α(N) を持つ離散的凝集モデルにおける自己組織的臨界性の出現を調査すること。
• 小クラスターの組合せ的構造を確率的グラフとして分析し、その臨界性領域を同定すること。
• アルドウスが提起した、典型的な粒子を含む成分の局所的弱収束が臨界ガルトン＝ワトソン木に一致するかを検証すること。
• 限界においてゲル化時間 Z と成分サイズ C(Z) の依存関係を検討すること。

提案手法

• 各粒子を i.i.d. の指数分布時計（レート 1/N）でモデル化し、サイズが α(N) 未満のクラスター間でのみ凝集を許容する。
• 流体極限を用いて、クラスターカンセントレーションの収束をスモルゥッホフ方程式の解に示し、乗法的核を用いる。
• 確率的グラフ技法を用いて小クラスターグラフをエドーシュ＝レニー・モデルとして分析し、時刻1以前は亜臨的で、以降は臨的であることを示す。
• ER グラフの探索プロセスを用いて、典型的な粒子を含む成分内の過剰エッジの期待数を計算する。
• ブラウン運動近似における大きな逸脱の不在を条件づけることで、閾値制約をモデル化する。
• 確率的グラフ理論および分岐過程の収束結果を活用し、クラスター成分の局所的弱収束を確立する。

実験結果

リサーチクエスチョン

• RQ1 動的閾値 α(N) を持つ離散的凝集モデルは、微視的スケールで自己組織的臨界性を示すか？
• RQ2 小クラスターグラフの構造は時間経過とともにどのように変化するのか？ それは亜臨的か、臨的か？
• RQ3 アルドウスの予想である、典型的な粒子を含む成分が局所的弱収束して臨界ガルトン＝ワトソン木に一致するか？
• RQ4 限界において、ゲル化時間 Z と成分サイズ C(Z) は漸近的に独立か？
• RQ5 α(N) が N²/³ より小さい場合、特に臨界窓内では、モデルはどのように振る舞うか？

主な発見

• N²/³ log⁴(N) ≤ α(N) の仮定の下で、小クラスターのグラフは漸近的にエドーシュ＝レニー確率的グラフとなり、時刻1以前は亜臨的で、以降は臨的である。
• 典型的な有限クラスターの極限分布は臨界ガルトン＝ワトソン木の分布と一致し、微視的スケールでの自己組織的臨界性が確認される。
• 典型的な粒子を含む成分 C(Z) は局所的弱収束で臨界ガルトン＝ワトソン木に収束し、アルドウスの予想が検証される。
• C(Z) 内の過剰エッジの期待数は N → ∞ の極限で Z²/12 に収束し、Z と C(Z) が漸近的に独立でないことが示唆される。
• モデルは α(N) ≍ N²/³ まで閾値を低下させてもロバストであるが、臨界窓内では条件づけが非自明になる。