[論文レビュー] Mean-field forest-fire models and pruning of random trees

Xavier Bressaud, Nicolas Fournier|arXiv (Cornell University)|Jan 31, 2012

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

ひとこと要約

本稿では、1次元のフォレストファイプロセスにインspiredされた平均場コアグラーション・フラグメンテーションモデルを研究している。粒子は2のレートでコアグラートし、質量iの粒子は (i−1)/n の割合で単位質量にフラグメンテーションする。n→∞のとき、典型的な粒子質量の平衡分布は、臨界的な2分岐ガルトン＝ワトソン木の葉のサイズ分布に収束する。一方、サイズバイアスをかけた典型的な粒子は、エアリー関数とその導関数の零点を含む極限プロファイルに収束する。

ABSTRACT

We consider a family of discrete coagulation-fragmentation equations closely related to the one-dimensional forest-fire model of statistical mechanics: each pair of particles with masses $i,j \in n$ merge together at rate 2 to produce a single particle with mass $i+j$, and each particle with mass $i$ breaks into $i$ particles with mass 1 at rate $(i-1)/n$. The (large) parameter $n$ controls the rate of ignition and there is also an acceleration factor (depending on the total number of particles) in front of the coagulation term. We prove that for each $n\in n$, such a model has a unique equilibrium state and study in details the asymptotics of this equilibrium as $n o \infty$: (I) the distribution of the mass of a typical particle goes to the law of the number of leaves of a critical binary Galton-Watson tree, (II) the distribution of the mass of a typical size-biased particle converges, after rescaling, to a limit profile, which we write explicitly in terms of the zeroes of the Airy function and its derivative. We also indicate how to simulate perfectly a typical particle and a size-biased typical particle, which allows us to give some probabilistic interpretations of the above results in terms of pruned Galton-Watson trees and pruned continuum random trees.

研究の動機と目的

• パrameter n によって制御される点火レートを持つ平均場コアグラーション・フラグメンテーションモデルの平衡状態の存在と一意性を確立すること。
• n→∞ のときの平衡分布の漸近的挙動を分析し、特に典型的およびサイズバイアスをかけた粒子質量分布に注目すること。
• ガルトン＝ワトソン木および連続的ランダム木における pruning 手続きを用いて、極限分布の確率的解釈を提供すること。
• エアリー関数に関連するドリフトを持つ拡散過程の局所時間の逆のラプラス指数の正確な表現を導出することにより、確率過程とスケーリング極限を結びつけること。

提案手法

• 質量 i, j の粒子対がレート2でコアグラートする平均場コアグラーション・フラグメンテーション系 (CF_n) を形式化する。質量 i の粒子は、レート (i−1)/n で i 個の単位質量粒子にフラグメンテーションする。
• 元のフォレストファイモデルにおける空間的相関を反映させるために、コアグラーション項に 1/∑c_l^n(t) の加速因子を導入する。
• 解析的手法を用いて、各有限 n に対して平衡測度の存在と一意性を証明する。
• スケーリング極限と漸近的解析を適用し、n→∞ における平衡の挙動を、特に典型的およびサイズバイアスをかけた粒子の質量分布に注目して研究する。
• ガルトン＝ワトソン木および連続的ランダム木における pruning 手続きを基にした完全なシミュレーションアルゴリズムを開発し、典型的およびサイズバイアスをかけた粒子の平衡分布からの正確なサンプリングを可能にする。
• 極限分布と臨界的な2分岐ガルトン＝ワトソン木の輪郭過程との関係を確立し、スケーリング極限を介して連続的ランダム木 (CRT) と関連付ける。

実験結果

リサーチクエスチョン

• RQ1 n→∞ のとき、平均場フォレストファイモデルの平衡状態における典型的な粒子の質量の漸近的分布は何か？
• RQ2 n→∞ の極限において、サイズバイアスをかけた典型的な粒子の質量分布はどのように振る舞い、その極限プロファイルは何か？
• RQ3 サイズバイアスをかけた粒子の極限分布は、エアリー関数やその導関数を用いて明示的に表現できるか？
• RQ4 コアグラーション・フラグメンテーション過程の平衡状態と、ガルトン＝ワトソン木または連続的ランダム木の pruning の間にどのような関係があるか？
• RQ5 極限プロファイルから導出されたラプラス指数と一致する最初の通過時間分布を持つ確率過程は存在するか？そのプロセスの性質は何か？

主な発見

• 各有限 n に対して、平均場コアグラーション・フラグメンテーション系 (CF_n) は一意な平衡状態を持つ。
• n→∞ のとき、典型的な粒子質量の分布は、臨界的な2分岐ガルトン＝ワトソン木の葉の数の分布に弱収束する。
• 適切なスケーリングを施したサイズバイアスをかけた典型的な粒子の質量分布は、エアリー関数 Ai とその導関数 Ai′ の零点を用いて明示的に表せる極限プロファイルに収束する。
• ガルトン＝ワトソン木における pruning 手続きを基にした完全なシミュレーションアルゴリズムが構築され、典型的およびサイズバイアスをかけた粒子の両方の平衡分布からの正確なサンプリングが可能になった。
• サイズバイアスをかけた粒子の極限過程は、ドリフト係数 β(x) = −sign(x)α Ai′(α|x|+a₁′)/Ai(α|x|+a₁′) を持つ拡散過程と関連しており、その逆局所時間のラプラス指数は ψ(λ) = β(2λ/α³) に等しい。
• 驚くべきことに、α = 2^{1/3} のとき、ドリフト係数 β と逆局所時間のラプラス指数 ψ(λ) が一致し、2つの異なる確率的対象の間の深いつながりが明らかになった。