[논문 리뷰] 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)|2022. 06. 08.

Theoretical and Computational Physics참고 문헌 52인용 수 3

한 줄 요약

이 연구는 2D 및 3D 격자에서 스토하스틱 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 값에 도달하지 못한 상태에서 발생함).

제안 방법

• 2D 및 3D 입방 격자에서 고정된 장애물 희석을 적용한 스토하스틱 SIR 모델의 몬테카를로 시뮬레이션.
• 면역화에 의해 유도된 격자 희석을 모의하기 위해 포아송 분포를 이용한 고정된 장애물 도입.
• 교차 행동을 탐지하기 위해 시간에 따라 변화하는 효과적 임계 지수 분석.
• 복구 비율을 체계적으로 변화시키고 근접 이웃으로의 이동을 도입하여 이동성 영향 평가.
• 동적 등방성 퍼콜레이션(DyIP) 보편성 클래스의 이론적 예측과 시뮬레이션 결과 비교.
• 소규모 분리된 감염 가능 집단의 전파 및 소멸 동역학 추적을 통해 회복 지연 효과 분리 분석.

실험 결과

연구 질문

• RQ1 해리스 기준이 예측하는 것처럼 고정된 장애물이 있는 상황에서, 현상이 SIR 모델의 임계 행동을 변화시키는가?
• RQ2 왜 고정된 장애물이 있는 SIR 모델 시뮬레이션에서 비보편적인 효과적 임계 지수를 보이는 장기적 교차 영역이 나타나는가?
• RQ3 왜 희박한 격자에서 소규모 분리된 감염 가능 개인 클러스터의 회복이 지연되는가?
• RQ4 복구 비율을 높이면 SIR 모델에서 교차 행동의 지속성이 어떻게 변화하는가?
• RQ5 근접 이웃으로의 이동을 통해 에이전트의 이동성을 도입하면 정적 SIR 모델에서 관찰된 교차 영역을 제거할 수 있는가?

주요 결과

• 해리스 기준은 점차적 DyIP 행동 수렴을 예측하지만, 희박한 SIR 체계에서는 DyIP 보편성 클래스와 다를 수 있는 효과적 임계 지수를 보이는 장기적 교차 영역이 관찰된다.
• 이 교차 영역은 포아송 분포를 통한 고정된 장애물에 의해 생성된 소규모 분리된 감염 가능 클러스터의 완전한 회복 지연으로 인해 발생한다.
• 격자 크기를 크게 증가시킬 때 비로소 효과적 지수는 기대되는 DyIP 값에 점차 수렴함을 확인하여, 이는 크기가 유한한 체계가 진정한 임계 영역에 도달하지 못할 수 있음을 시사한다.
• 복구 비율을 크게 높이면 교차 영역이 극적으로 단축되며, 이는 소규모 클러스터의 회복 지연이 근본 원인임을 뒷받침한다.
• 에이전트의 근접 이웃 이동을 허용함으로써도 교차 영역이 사라지며, 이는 클러스터의 회복 역학이 일시적 행동을 지배함을 추가로 확인한다.
• 결과적으로, 시뮬레이션 또는 관측 데이터에서 관찰되는 비보편적인 효과적 지수는 본질적인 보편성 클래스 변화 때문이 아니라, 유한 체적 효과와 완전한 회복 동역학 미흡으로 인한 것임을 시사한다.