[論文レビュー] An example of Brunet-Derrida behavior for a branching-selection particle system on $\Z$

Jean Bérard|ArXiv.org|Oct 30, 2008

Stochastic processes and statistical mechanics参考文献 10被引用数 6

ひとこと要約

この論文は、ℤ 上の離散的ブランチャイング・セレクション粒子系において、Brunet-Derrida の挙動を厳密に確立している。$0 < p < 1/2$ のとき、無限大人口速度 $v_\infty(p)$ と有限 $N$ における速度 $v_N(p)$ の差が $\log(N)^{-2}$ の割合で減少することを証明し、長年の数値的・ヒューリスティックな研究から予想されていた結果を裏付けた。この結果は、カップリングの議論、マルティンゲール不等式、およびブランチャイング過程における大偏差推定に基づいて得られた。

ABSTRACT

We consider a branching-selection particle system on $\Z$ with $N \geq 1$ particles. During a branching step, each particle is replaced by two new particles, whose positions are shifted from that of the original particle by independently performing two random walk steps according to the distribution $p δ_{1} + (1-p) δ_{0}$, from the location of the original particle. During the selection step that follows, only the N rightmost particles are kept among the 2N particles obtained at the branching step, to form a new population of $N$ particles. After a large number of iterated branching-selection steps, the displacement of the whole population of $N$ particles is ballistic, with deterministic asymptotic speed $v_{N}(p)$. As $N$ goes to infinity, $v_{N}(p)$ converges to a finite limit $v_{\infty}(p)$. The main result is that, for every $0

研究の動機と目的

• ℤ 上の離散的粒子系における Brunet-Derrida の挙動を、数学的に厳密に証明すること。特に、$v_N(p)$ が無限大人口における極限 $v_\infty(p)$ に収束する際の速度が、$\log(N)^{-2}$ の遅いレートで行われることを示すこと。
• 有限の集団サイズ $N$ がブランチャイング・セレクション過程におけるフロント速度に与える影響を、特に $0 < p < 1/2$ の領域で理解を深めること。
• 速度差 $v_\infty(p) - v_N(p)$ に対して、鋭い上界および下界を確立し、Brunet と Derrida が非厳密な議論によって予測した $\log(N)^{-2}$ スケーリングを確認すること。
• $p = 1/2$ および $p > 1/2$ における挙動の遷移を分析し、収束レートがそれぞれ $N^{-1}$ および $N$ に指数的に依存するものに変わることを明らかにすること。

提案手法

• 各粒子が $p\delta_1 + (1-p)\delta_0$ に従って独立に移動し、2つに分岐した後、右端から $N$ 個の粒子が選別される、ℤ 上の離散時間マルコフ過程として粒子系をモデル化すること。
• カップリング技術を用いて、粒子系をブランチャイング過程 $V_k$ および二項分布マルティンゲール $Z_k$ と比較し、確率的支配関係とモーメント推定を可能にすること。
• ドーブの最大不等式を適用し、$m \propto \log(N)^3$ の時間区間内に粒子数が平均から著しく逸脱する確率を抑え込むこと。
• 測度の変更とブラウン運動スケーリングを用いて、右端の粒子の位置に関連するレアイベントの確率を推定し、$\log(N)^{-2}$ スケーリングを導出すること。
• 大偏差推定と、漂移 $p$ と有限 $N$ によるカットオフ効果の間の相互作用を用いて、フロント速度 $v_N(p)$ の漸近的挙動を分析すること。
• 系の緩和に重要な役割を果たす $\log(N)^3$ の時間スケールを特定し、関連するモデルにおける系統的共通祖先の時間と一致させること。

実験結果

リサーチクエスチョン

• RQ1 $0 < p < 1/2$ のとき、有限集団におけるフロント速度の $\log(N)^{-2}$ 収束レートが、厳密に成り立つかどうか。
• RQ2 $p = 1/2$ または $p > 1/2$ のとき、$v_N(p)$ から $v_\infty(p)$ への収束レートはどのように変化するか。
• RQ3 粒子位置のレアフラクチュエーションやブランチャイング過程といった確率的メカニズムによって、$\log(N)^{-2}$ スケーリングが説明可能か。
• RQ4 時間スケール $m \propto \log(N)^3$ は、系の緩和ダイナミクスにおいてどのような役割を果たしているか。

主な発見

• $0 < p < 1/2$ に対して、すべての $N \to \infty$ において正の定数 $C_*(p)$ と $C^*(p)$ が存在し、$v_\infty(p) - v_N(p) \asymp \log(N)^{-2}$ が成り立つ。
• 収束レート $v_\infty(p) - v_N(p)$ は、下界 $C_*(p)\log(N)^{-2}$ と上界 $C^*(p)\log(N)^{-2}$ で抑えられ、$\log(N)^{-2}$ スケーリングが確認された。
• $p = 1/2$ のとき、速度差は $N^{-1}$ の割合で減少し、$p < 1/2$ の領域よりも速い収束を示している。
• $p > 1/2$ のとき、速度差は $N$ に対して指数的に減少し、$-\log(v_\infty(p) - v_N(p)) \asymp N$ となる。これは極めて速い収束を示している。
• 時間スケール $m \propto \log(N)^3$ は、系の緩和にとって臨界的であることが特定され、系統的共通祖先の時間と一致した。
• $\log(N)^{-2}$ スケーリングは、漂移 $\gamma^{\pm\Delta}$ と大偏差コスト $\exp(-m/\Delta^2)$ のバランスによって生じる。これにより $\Delta \propto m^{1/3}$ となり、結果として $v \propto \Delta/m \propto \log(N)^{-2}$ となる。