전염병 확산 규모 SIS 모델 전파 횟수 정확 해와 임계 현상

초록

이 논문은 SIS(감수성‑감염‑감수성) 모델에서 전파 횟수의 평균을 정확히 구하고, 인구 규모가 큰 경우의 하위·상위·임계 영역별 점근적 거동을 분석한다. 특히 임계 영역에서 전파 횟수 분포를 구해 작은 n에서는 n³⁄² 특이성을, 큰 n에서는 지수적 감소를 보이며, 임계점으로부터의 거리와 관련된 지수의 변화를 제시한다.

상세 분석

논문은 먼저 SIS 모델을 마스터 방정식 형태로 기술하고, 전파 횟수 n을 상태 변수로 포함한 확률 분포 Pₙ(t)를 정의한다. 전파가 일어날 때마다 감염자가 하나 늘어나고 회복이 일어나면 감소하는 전이율을 이용해, 전파 횟수의 평균 ⟨n⟩에 대한 재귀 관계식을 도출한다. 이 재귀식을 생성함수 G(z)=∑ₙPₙzⁿ 형태로 변환하면, 비선형 미분 방정식이 얻어지고, 경계 조건을 이용해 정확한 해를 구한다. 해는 인구 규모 N과 전염성 파라미터 β/γ(감염률 대비 회복률)의 함수로 표현되며, 특히 전염성 λ=β/γ가 1을 기준으로 어떻게 변하는지에 초점을 맞춘다.

대규모 N에 대해 점근적 분석을 수행한다. λ<1(아래 임계)에서는 평균 전파 횟수가 O(1) 수준으로 제한되고, λ>1(위 임계)에서는 ⟨n⟩이 N에 비례해 선형적으로 성장한다. 임계점 λ≈1 근처에서는 스케일링 변수 x=(λ−1)N^{1/2}를 도입해 전파 횟수 분포가 보편적인 형태로 수렴함을 보인다. 이때 확률 분포 Pₙ는 n^{3/2} 형태의 특이성을 보이며, 작은 n에서 Pₙ∝n^{−3/2}가 된다. 이는 임계 현상에서의 급격한 플럭투에이션을 의미한다.

또한 논문은 큰 n 영역에서의 지수적 감소를 분석한다. 전파가 지속될 확률이 감소함에 따라 Pₙ∼exp(−α(λ) n) 형태가 되며, α(λ)라는 감쇠 계수는 λ와의 거리 |λ−1|에 따라 변한다. λ가 1보다 크게 떨어질수록 α는 작아져 긴 꼬리를 만들고, λ가 1보다 작게 멀어질수록 α는 급격히 커져 분포가 급격히 억제된다. 이러한 결과는 임계점 근처에서 전염병이 소멸하거나 폭발적으로 확산되는 두 가지 가능한 시나리오를 정량적으로 설명한다.

마지막으로 저자들은 정확 해와 점근식이 수치 시뮬레이션과 일치함을 Monte Carlo 실험을 통해 검증한다. 특히 임계 영역에서 관측된 n^{3/2} 특이성과 지수적 꼬리는 이론적 예측과 거의 일치한다는 점에서 모델의 신뢰성을 높인다. 전체 분석은 SIS 모델이 단순한 평균장 근사보다 더 복잡한 플럭투와 스케일링 행동을 내포하고 있음을 보여준다.