[論文レビュー] Effects of lattice dilution on the non-equilibrium phase transition in the stochastic Susceptible-Infectious-Recovered model

Ruslan I. Mukhamadiarov, Uwe C. Täuber|arXiv (Cornell University)|Jun 8, 2022

Theoretical and Computational Physics参考文献 52被引用数 3

ひとこと要約

本研究は、2次元および3次元格子上の確率的SIRモデルにおける、非平衡な活性状態から吸収状態への転移に、ワクチン接種を模倣するサイト希釈が与える影響を調査する。ハリス基準は普遍的なDyIP臨界行動を予測しているが、シミュレーションでは、小規模な分離した感受性集団の回復が遅れるため、漸近的なDyIP値とは異なる有効臨界指数を示す長期間のクロスオーバー領域が観測される。このクロスオーバー領域は、回復率を上昇させたり、エージェントの移動性を導入したりすることで消失し、有限系における顕在的な非普遍的スケーリングが解消される。

ABSTRACT

We investigate how site dilution, as would be introduced by immunization, affects the properties of the active-to-absorbing non-equilibrium phase transition in the paradigmatic Susceptible-Infectious-Recovered (SIR) model on regular cubic lattices. According to the Harris criterion, the critical behavior of the SIR model, which is governed by the universal scaling exponents of the dynamic isotropic percolation (DyIP) universality class, should remain unaltered after introducing impurities. However, when the SIR reactions are simulated for immobile agents on two- and three-dimensional lattices subject to quenched disorder, we observe a wide crossover region characterized by varying effective exponents. Only after a sufficient increase of the lattice sizes does it becomes clear that the SIR system must transition from that crossover regime before the effective critical exponents asymptotically assume the expected DyIP values. We attribute the appearance of this exceedingly long crossover to a time lag in a complete recovery of small disconnected clusters of susceptible sites which are apt to be generated when the system is prepared with Poisson-distributed quenched disorder. Finally, we demonstrate that this transient region becomes drastically diminished when we significantly increase the value of the recovery rate or enable diffusive agent mobility through short-range hopping.

研究の動機と目的

• クエンチド不純度としてのサイト希釈が、SIRモデルにおける非平衡相転移に与える影響を理解すること。
• クエンチド不純度が存在するSIRモデルにおいて、ハリス基準が成立するかどうかを検証すること。
• シミュレーションで観測された臨界指数の長期間クロスオーバー領域の原因を特定すること。
• 回復率とエージェントの移動性がクロスオーバー行動の持続性に与える影響を特定すること。
• 有限系における有効指数が、漸近的なDyIP値に達していないために非普遍的であると誤解される理由を明確にすること。

提案手法

• クエンチドサイト希釈を用いた2次元および3次元立方格子上の確率的SIRモデルのモンテカルロシミュレーション。
• ワクチン接種に起因する格子希釈を模倣するため、ポアソン分布を用いたクエンチド不純度の使用。
• クロスオーバー行動を検出するため、時間依存の有効臨界指数の分析。
• 回復率の系統的変動と近隣格子への移動（ホッピング）の導入により、移動性効果を評価。
• ダイナミックアイソトロープペルコレーション（DyIP）普遍クラスからの理論的予測とシミュレーション結果を比較。
• 小規模で分離した感受性集団における感染拡散および消滅ダイナミクスの追跡により、回復遅延効果を分離。

実験結果

リサーチクエスチョン

• RQ1 ハリス基準が予測するように、クエンチド不純度が存在するSIRモデルにおいて、サイト希釈が臨界行動を変えるか？
• RQ2 クエンチド不純度を伴うSIRモデルのシミュレーションが、非普遍的有効臨界指数を示す長期間クロスオーバー領域を示す理由は何か？
• RQ3 希釈格子における小規模で分離した感受性個体のクラスタの回復が遅れる原因は何か？
• RQ4 回復率を上昇させると、SIRモデルにおけるクロスオーバー行動の持続性にどのような影響を与えるか？
• RQ5 エージェントが近隣格子へのホッピングを可能にすることで、静的SIRモデルで観測されたクロスオーバー領域を排除できるか？

主な発見

• ハリス基準が漸近的DyIP行動への収束を予測しているにもかかわらず、希釈SIR系ではDyIP普遍クラスとは異なる有効臨界指数を示す長期間のクロスオーバー領域が観測される。
• このクロスオーバーは、ポアソン分布を用いたクエンチド不純度によって生成される小規模で分離した感受性個体クラスタの完全回復に時間遅れが生じるため生じる。
• 格子サイズを大幅に増加させることで、有効指数が漸近的に期待されるDyIP値に近づくことが示され、有限系では真の臨界領域に到達しない可能性がある。
• 回復率を上昇させると、クロスオーバー領域が著しく短縮され、小クラスタの回復遅延が根本的原因であるという仮説が裏付けられる。
• エージェントに近隣格子へのホッピングを許可することで、クロスオーバー領域が消失し、クラスタ回復ダイナミクスが一時的挙動を支配することをさらに確認できる。
• これらの結果は、シミュレーションや観測データにおける顕在的な非普遍的臨界指数が、普遍クラスの本質的変化ではなく、有限サイズ効果および完全な回復ダイナミクスの不完全さに起因する可能性があることを示唆する。