[論文レビュー] Phase transition for a contact process with random slowdowns

Kuoch, Kevin|arXiv (Cornell University)|Jan 5, 2015

Stochastic processes and statistical mechanics参考文献 57被引用数 5

ひとこと要約

本稿では、動的確率的環境下での接触過程に、不妊個体が存在する際に繁殖速度が低下するランダムな遅延を導入したモデルを検討する。グラフィカル表現とブロック構成を用いて、野生個体群の生存に関する相転移を証明した。不妊個体の放出率が臨界閾値を超えると絶滅が発生するが、それ未満であれば生存が可能である。これは、害虫駆除における不妊昆虫法の厳密な確率的モデルを提供する。

ABSTRACT

Motivated by a model of an area-wide integrated pest management, we develop an interacting particle system evolving in a random environment. It is a generalised contact process in which the birth rate takes two possible values, determined either by a dynamic or a static random environment. Our goal is to understand the phase diagram of both models by identifying the mechanisms that permit coexistence or extinction of the process.

研究の動機と目的

• 不妊昆虫法を、動的確率的環境下の格子上における確率的空間過程としてモデル化すること。
• 野生個体群が不妊個体との競争によって生存または絶滅する条件を分析すること。
• 野生個体群の生存と絶滅を分かつ臨界的放出率を同定すること。
• 確率的漸近的理論と再帰的スケーリング技術を用いて、ランダムな遅延を伴う接触過程における相転移行動を確立すること。

提案手法

• 状態 {0,1,2,3}（空、野生個体のみ、不妊個体のみ、共存）を用いて、$\mathbb{Z}^d$ 上の多型接触粒子系としてプロセスを定式化する。
• 遷移率の定義：状態1では出生率 $\lambda_1$、状態3では $\lambda_2 < \lambda_1$ の出生率、死亡率は1、不妊個体の放出率は $r$。
• 確率的漸近的理論を用いるために、確率的グラフ構造を用いてプロセスを定義するグラフィカル表現を構築する。
• 長時間挙動と生存確率の分析に、ブロック構成と動的再帰的スケーリング技術を適用する。
• カップリング論法と単調性性質を用いて、臨界閾値 $r_c$ の存在と一意性を確立する。
• 環境のモーメント条件とエルゴディック理論を用いて、相転移の明示的境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 不妊昆虫の放出に起因するランダムな遅延を伴う接触過程において、野生個体群が生存または絶滅する条件は何か？
• RQ2 臨界的放出率 $r_c$ は、出生率 $\lambda_1$ と $\lambda_2$ にどのように依存するか？
• RQ3 CPRSモデルの対称的および非対称的バージョンにおいて、相転移が存在するか？
• RQ4 環境のモーメント条件とエルゴディック性質を用いて、臨界閾値を境界づけることができるか？
• RQ5 不妊個体の競争能力が、野生個体群の生存に与える影響は何か？

主な発見

• 固定された $\lambda_1 > \lambda_c(d) > \lambda_2$ に対して、一意な臨界的放出率 $r_c \in (0, \infty)$ が存在し、$r < r_c$ のとき野生個体群は生存し、$r \geq r_c$ のとき絶滅する。
• 非対称的状況では、不妊個体が新たな不妊個体を生むことも許容されるが、一意な臨界閾値 $s_c$ が存在し、$r < s_c$ のとき生存が成立する。
• $\lambda_1 > 1 + \sqrt{2} \approx 2.41$ かつ $\lambda_2 < 1 + \sqrt{2}$ のとき、$r < \frac{\lambda_2(\lambda_1 - \sqrt{2} - 1)}{\lambda_1(\lambda_2 - \sqrt{2} - 1)}$ であれば生存が保証され、明示的な上界が得られる。
• 1次元の基本接触過程の臨界値 $\lambda_c$ は、モデルの推定に基づき $\lambda_c \leq 1 + \sqrt{2}$ を満たす。
• 数値的境界として、相転移区間が提示される：$\lambda_1 = 1000, \lambda_2 = 0.8$ の場合、遷移は $[0.49, 4995)$ に発生する；$\lambda_1 = 2, \lambda_2 = 1.4$ の場合、すべての $\mathbb{R}_+$ で遷移が発生する。
• 野生個体群の適応度が上昇するか、不妊個体の競争力が低下するほど、臨界閾値は上昇する。これは生物学的直観を裏付ける。