[論文レビュー] The Contact Process on Random Graphs and Galton-Watson Trees

Xiangying Huang, Rick Durrett|arXiv (Cornell University)|Oct 14, 2018

Complex Network Analysis Techniques参考文献 31被引用数 8

ひとこと要約

本稿は、ランダムグラフおよびギャルトン＝ウォーソン木における接触過程を分析し、スターチャートにおける生存時間の精密な推定を導入することで、臨界感染率の tighter な境界を導出する。部分指数的（subexponential）な子孫分布では局所的生存が確実（$ olambda_2 = 0$）であり、幾何分布の子孫では $ olambda_2$ が従来の推定より顕著に小さいことが示され、配置モデルグラフにおける $ olambda_c(n)$ の収束速度も向上している。

ABSTRACT

The key to our investigation is an improved (and in a sense sharp) understanding of the survival time of the contact process on star graphs. Using these results, we show that for the contact process on Galton-Watson trees, when the offspring distribution (i) is subexponential the critical value for local survival $λ_2=0$ and (ii) when it is geometric($p$) we have $λ_2 \le C_p$, where the $C_p$ are much smaller than previous estimates. We also study the critical value $λ_c(n)$ for "prolonged persistence" on graphs with $n$ vertices generated by the configuration model. In the case of power law and stretched exponential distributions where it is known $λ_c(n) o 0$ we give estimates on the rate of convergence. Physicists tell us that $λ_c(n) \sim 1/Λ(n)$ where $Λ(n)$ is the maximum eigenvalue of the adjacency matrix. Our results show that this is not correct.

研究の動機と目的

• スターチャートにおける接触過程の生存時間の理解を深めること。これは、より複雑なランダムグラフや木構造を分析するための基盤的ツールである。
• 特に部分指数的および幾何分布の子孫分布に対して、ギャルトン＝ウォーソン木における局所的生存の臨界値 $ olambda_2$ を特定すること。
• 配置モデルで生成される有限ランダムグラフにおける、長期間にわたる持続のための臨界値 $ olambda_c(n)$ を分析すること。特に、パワー則および伸びた指数的（stretched exponential）度分布に対して。
• 物理学的ヒューリスティック $ olambda_c(n) \sim 1/\Lambda(n)$（$ Lambda(n)$ は隣接行列の最大固有値）が、さまざまな度分布に対してどれほど正確であるかを評価すること。
• $ olambda_c(n)$ に対する改善された上界を確立し、パワー則および伸びた指数的ケースにおける $0$ への収束速度を定量化すること。

提案手法

• 連続時間の接触過程を、中心ノードが感染している瞬間にのみ注目する離散時間の増減過程（birth-death chain）に還元することで、解析を単純化する。
• 感染率を $pk$（$p = \lambda/(1+2\lambda)$）で上限化する変更版の過程を導入し、感染したリーブの数が有界に保たれ、マルティンググールの技法が適用可能になるようにする。
• 埋め込み連鎖から導かれる上向きマルティンググール（supermartingale）に「停止定理」を適用し、絶滅確率の上界を導出し、生存時間の推定を行う。
• 個々のスターチャートの生存を偏ったランダムウォークと比較するカップリングの議論を用い、感染したスターチャートの数が正のドリフトを持つ超臨界ランダムウォークを支配することを示す。
• 大偏差推定および指数的バウンドを用いて、スターチャート間の感染移動サイクル中の絶滅確率を制御する。
• スターチャートの生存時間と、感染移動に成功するのに必要な時間とを等置することで、$ olambda_c(n)$ の漸近的境界を導出し、$ lambda^2 k$ スケーリングに基づく臨界値推定を行う。

実験結果

リサーチクエスチョン

• RQ1 $k$ 個のリーブを持つスターチャートにおける接触過程の生存時間の正確な漸近的挙動は何か？
• RQ2 特に部分指数的および幾何分布の場合に、ギャルトン＝ウォーソン木における局所的生存の臨界値 $ olambda_2$ は子孫分布にどのように依存するか？
• RQ3 パワー則および伸びた指数的度分布を持つ配置モデルランダムグラフにおいて、$ olambda_c(n)$ が $0$ に収束する速度は何か？
• RQ4 物理学的ヒューリスティック $ olambda_c(n) \sim 1/\Lambda(n)$ は、すべての度分布に対して正確か、それとも特定の状況で失敗するか？
• RQ5 スターベースの構成要素とランダムウォークカップリングを用いて、ランダムグラフにおける接触過程の生存時間の下界を確立できるか？

主な発見

• ギャルトン＝ウォーソン木における部分指数的子孫分布では、局所的生存の臨界値は $ olambda_2 = 0$ であり、任意の $ lambda > 0$ に対して正の確率で生存することが保証される。
• パラメータ $p$ の幾何分布子孫では、臨界値 $ olambda_2$ は $C_p$ で上界が与えられ、$C_p$ は従来の推定より顕著に小さい。
• パワー則度分布（$P(d \geq m) \sim m^{-a}$）を持つ配置モデルグラフでは、$ olambda_c(n) \sim n^{-1/2a}$ であり、本稿は $n$ が十分に大きい場合に高い確率でこのレートを確認している。
• 伸びた指数的度分布（$P(d \geq m) \sim \exp(-m^{1/b})$）では、$ olambda_c(n) \sim \log^{-b/2}n$ であり、本稿は物理学的ヒューリスティック $ olambda_c(n) \sim 1/\Lambda(n)$ がこの場合に失敗することを示している。
• このようなグラフにおける接触過程の生存時間は、$ exp(O(n^{\epsilon}))$ 以上である（ある $ epsilon > 0$ に対して）、これは感染スターチャートの数が超臨界ランダムウォークを支配することを示すことで達成される。
• 本稿は、伸びた指数的グラフにおいて $ olambda_c(n)$ が $1/\Lambda(n)$ よりも遅く減少することを確立しており、ヒューリスティックがこの領域で不正確であることを示している。