Epidemic models on lattices

Last updated
Spatial SIR model simulation. Each cell can infect its eight immediate neighbors. SIR model simulated using python.gif
Spatial SIR model simulation. Each cell can infect its eight immediate neighbors.

Classic epidemic models of disease transmission are described in Compartmental models in epidemiology. Here we discuss the behavior when such models are simulated on a lattice. Lattice models, which were first explored in the context of cellular automata, act as good first approximations of more complex spatial configurations, although they do not reflect the heterogeneity of space (e.g. population density differences, urban geography and topographical differentiations). [1] Lattice-based epidemic models can also be implemented as fixed agent-based models. [2]

Contents

Introduction

The mathematical modelling of epidemics was originally implemented in terms of differential equations, which effectively assumed that the various states of individuals were uniformly distributed throughout space. To take into account correlations and clustering, lattice-based models have been introduced. Grassberger [3] considered synchronous (cellular automaton) versions of models, and showed how the epidemic growth goes through a critical behavior such that transmission remains local when infection rates are below critical values, and spread throughout the system when they are above a critical value. Cardy and Grassberger [4] argued that this growth is similar to the growth of percolation clusters, which are governed by the "dynamical percolation" universality class (finished clusters are in the same class as static percolation, while growing clusters have additional dynamic exponents). In asynchronous models, the individuals are considered one at a time, as in kinetic Monte Carlo or as a "Stochastic Lattice Gas."[ citation needed ]

SIR model

In the "SIR" model, there are three states:

  • Susceptible (S) -- has not yet been infected, and has no immunity
  • Infected (I)-- currently "sick" and contagious to Susceptible neighbors
  • Removed (R), where the removal from further participation in the process is assumed to be permanent, due to immunization or death

It is to be distinguished from the "SIS" model, where sites recover without immunization, and are thus not "removed".

The asynchronous simulation of the model on a lattice is carried out as follows:

  • Pick a site. If it is I, then generate a random number x in (0,1).
  • If x < c then let I go to R.
  • Otherwise, pick one nearest neighbor randomly. If the neighboring site is S, then let it become I.
  • Repeat as long as there are S sites available.

Making a list of I sites makes this run quickly.

The net rate of infecting one neighbor over the rate of removal is λ = (1-c)/c.

For the synchronous model, all sites are updated simultaneously (using two copies of the lattice) as in a cellular automaton.

Latticezccλc = (1 - cc)/cc
2-d asynchronous SIR model triangular lattice60.199727(6) [5] 4.0068(2)
2-d asynchronous SIR model square lattice40.1765(5), [6] 0.1765005(10) [7] 4.66571(3)
2-d asynchronous SIR model honeycomb lattice30.1393(1) [5] 6.179(5)
2-d synchronous SIR model square lattice40.22 [8] 3.55
2-d asynchronous SIR model on Penrose lattice0.1713(2) [9]
2-d asynchronous SIR model on Ammann-Beenker lattice0.1732(5) [9]
2-d asynchronous SIR model on random Delaunay triangulations0.1963(3) [10]

Contact process (asynchronous SIS model)

I → S with unit rate; S → I with rate λnI/z where nI is the number of nearest neighbor I sites, and z is the total number of nearest neighbors (equivalently, each I attempts to infect one neighboring site with rate λ)

(Note: S → I with rate λn in some definitions, implying that lambda has one-fourth the values given here).

The simulation of the asynchronous model on a lattice is carried out as follows, with c = 1 / (1 + λ):

  • Pick a site. If it is I, then generate a random number x in (0,1).
  • If x < c then let I go to S.
  • Otherwise, pick one nearest neighbor randomly. If the neighboring site is S, then let it become I.
  • Repeat

Note that the synchronous version is related to the directed percolation model.

Latticezλc
1-d23.2978(2), [11] 3.29785(2) [12]
2-d square lattice41.6488(1), [13] 1.64874(2), [14] 1.64872(3), [11] 1.64877(3) [15]
2-d triangular lattice61.54780(5) [16]
2-d Delaunay triangulation of Voronoi Diagram6 (av)1.54266(4) [16]
3-d cubic lattice61.31685(10), [17] 1.31683(2), [11] 1.31686(1) [15]
4-d hypercubic lattice81.19511(1) [11]
5-d hypercubic lattice101.13847(1) [11]

See also

References

  1. von Csefalvay, Chris (2023), "Spatial dynamics of epidemics" , Computational Modeling of Infectious Disease, Elsevier, pp. 257–303, doi:10.1016/b978-0-32-395389-4.00017-7, ISBN   978-0-323-95389-4 , retrieved 2023-03-02
  2. von Csefalvay, Chris (2023), "Agent-based modeling" , Computational Modeling of Infectious Disease, Elsevier, pp. 305–375, doi:10.1016/b978-0-32-395389-4.00018-9, ISBN   978-0-323-95389-4 , retrieved 2023-03-02
  3. Grassberger, Peter (1983). "On the critical behavior of the general epidemic process and dynamical percolation". Mathematical Biosciences. 63 (2): 157–172. doi:10.1016/0025-5564(82)90036-0.
  4. Cardy, John; Grassberger, Peter (1985). "Epidemic models and percolation". J. Phys. A. 18 (6): L267. Bibcode:1985JPhA...18L.267C. doi:10.1088/0305-4470/18/6/001.
  5. 1 2 Tomé, Tânia; David de Souza; Robert M. Ziff (2013). "Correlations and thresholds for the SIR process on lattices". Preprint.
  6. de Souza, David; Tânia Tomé (2010). "Stochastic lattice gas model describing the dynamics of the SIRS epidemic process". Physica A. 389 (5): 1142–1150. arXiv: 0908.1296 . Bibcode:2010PhyA..389.1142D. doi:10.1016/j.physa.2009.10.039. S2CID   17145631.
  7. Tomé, Tânia; Robert Ziff (2010). "On the critical point of the Susceptible-Infected-Recovered model". Physical Review E. 82 (5): 051921. arXiv: 1006.2129 . Bibcode:2010PhRvE..82e1921T. doi:10.1103/PhysRevE.82.051921. PMID   21230514. S2CID   28861135.
  8. Arashiro, Everaldo; Tânia Tomé (2007). "The threshold of coexistence and critical behaviour of a predator–prey cellular automaton". J. Phys. A. 40 (5): 887–900. arXiv: cond-mat/0607360 . Bibcode:2007JPhA...40..887A. doi:10.1088/1751-8113/40/5/002. S2CID   54021762.
  9. 1 2 Santos, G. B. M.; Alves, T. F. A.; Alves, G. A.; Macedo-Filho, A. (2020). "Epidemic outbreaks on two-dimensional quasiperiodic lattices". Physics Letters A. 384 (2): 126063. arXiv: 1901.01403 . Bibcode:2020PhLA..38426063S. doi:10.1016/j.physleta.2019.126063. S2CID   119399157.
  10. Alves, T. F. A.; Alves, G. A.; Macedo-Filho, A. (2019-01-10). "Asynchronous SIR model on Two-Dimensional Random Delaunay Lattices". arXiv: 1901.03029 [cond-mat.stat-mech].
  11. 1 2 3 4 5 Sabag, Munir M. S.; Mário J. de Oliveira (2002). "Conserved contact process in one to five dimensions". Phys. Rev. E. 66 (3): 036115. Bibcode:2002PhRvE..66c6115S. doi:10.1103/PhysRevE.66.036115. PMID   12366192.
  12. Dickman, Ronald; I. Jensen (1993). "Time-dependent perturbation theory for non-equilibrium lattice models". J. Stat. Phys. 71 (1/2): 89–127. Bibcode:1993JSP....71...89J. CiteSeerX   10.1.1.540.2166 . doi:10.1007/BF01048090. S2CID   46519524.
  13. Moreira, Adriana; Ronald Dickman (1996). "Critical dynamics of the contact process with quenched disorder". Phys. Rev. E. 54 (4): R3090 –R3093. arXiv: cond-mat/9604148 . Bibcode:1996PhRvE..54.3090M. doi:10.1103/PhysRevE.54.R3090. PMID   9965620. S2CID   15905118.
  14. Vojta, Thomas; Adam Fraquhar; Jason Mast (2009). "Infinite-randomness critical point in the two-dimensional disordered contact process". Phys. Rev. E. 79 (1): 011111. arXiv: 0810.1569 . Bibcode:2009PhRvE..79a1111V. doi:10.1103/PhysRevE.79.011111. PMID   19257005.
  15. 1 2 Dickman, Ronald (1999). "Reweighting in nonequilibrium simulations". Phys. Rev. E. 60 (3): R2441 –R2444. arXiv: cond-mat/9902304 . Bibcode:1999PhRvE..60.2441D. doi:10.1103/PhysRevE.60.R2441. PMID   11970171. S2CID   17225358.
  16. 1 2 de Oliveira, Marcelo M.; S. G. Alves; S. C. Ferreira; Ronald Dickman (2008). "Contact process on a Voronoi triangulation". Phys. Rev. E. 78 (3): 031133. arXiv: 0810.0240 . Bibcode:2008PhRvE..78c1133D. doi:10.1103/PhysRevE.78.031133. PMID   18851019. S2CID   34027358.
  17. Moreira, Adriana G.; Ronald Dickman (1992). "Critical behavior of the three-dimensional contact process". Phys. Rev. E. 45 (2): R563 –R566. Bibcode:1992PhRvA..45..563J. doi:10.1103/PhysRevA.45.R563. PMID   9907104.

Further reading