#### For citation:

Shabunin A. V. SIRS-model with dynamic regulation of the population: Probabilistic cellular automata approach. *Izvestiya VUZ. Applied Nonlinear Dynamics*, 2019, vol. 27, iss. 2, pp. 5-20. DOI: 10.18500/0869-6632-2019-27-2-5-20

# SIRS-model with dynamic regulation of the population: Probabilistic cellular automata approach

Aim. Construction a model of infection spread in the form of a lattice of stochastic cellular automata which can demonstrate nontrivial oscillating regimes; investigation of its dynamics and comparison with the mean-field model. Method. Numerical simulation of the square lattice of cellular automata by the Monte Carlo approach, theoretical and numerical study of the structure of the phase space of its mean-field model. Results. A modified SIRS-model of epidemic propagation has been proposed in the form of a lattice of stochastic cellular automata. It uses dynamic regulation of the population with limit on the maximum number of individuals and takes into account the influence of the disease on reproduction processes. The study has discovered that, depending on the control parameters, the model demonstrates four different steady-state regimes: (a) total dying of the population, (b) stationary course of the disease, (c) complete recovery of the deases and (d) self-sustaining fluctuations of the infected, which accompanied by fluctuations of the population. The last regime is manifested near the border of the zone of complete recovery and is characterized by irregular bursts in the number of ill with a clearly distinguishable periodic component. At periodic modulation of the parameters the model demonstrates noisy periodic or quasi-periodic oscillations. Discussion. In general, the cellular automata model behaves similarly to the mean field equations, however there are quantitative and qualitative differences.The first ones consist in small shifts in the average concentrations compared to theoretically predicted values. These shifts are significantly reduced when individuals in the population can migrate. This effect can be explained by the greater homogeneity of populations in this case. The qualitative difference is that, the model of cellular automata demonstrates the effects, which are absent in the mean-field equations: a threshold of ill individuals necessary for the epidemic spread, and the oscillatory regime. A suggested cause of the last effect is the probabilistic nature of cellular automata, which lead to random fluctuations of the concentrations. The last ones can be amplified just as it happens in excitable systems under external noise.

- Bailey N. The Mathematical Approach to Biology and Medicine. London: John Wiley and Sons, 1967.
- Hethcote H. W. The mathematics of infectious diseases. SIAM Review, 2000, vol. 42, no. 4, pp. 599–653.
- Anderson R., May R. Infectios Deseases of Humans – Dynamics and Control. Oxford: Oxford University Press, 1991.
- Kermack W., McKendrick A. A contribution to the mathematical theory of epidemics. Proc. R. Soc., 1927, vol. A115, pp. 700–721.
- Kobrinsky N.E., Trahtenberg B.A. Introduction to the Theory of Finite Automata. Moscow: Fizmatgiz, 1962 (in Russian).
- Tzetlin M.L. Research on the Theory of Automata and Modeling of Biological Systems. Moscow: Nauka, 1969 (in Russian).
- Toffoli T., Margolus N. Cellular Automata Machines: A New Environment for Modeling. Cambridge: MIT Press Cambridge, 1987.
- Boccara N., Cheong K. Automata network SIR models for the spread of infectious diseases in populations of moving individuals. Journal of Physics A: Mathematical and General, 1992, vol. 25, no. 9, p. 2447.
- Sirakoulis G.C., Karafyllidis I., Thanailakis A. A cellular automaton model for the effects of population movement and vaccination on epidemic propagation. Ecological Modelling, 2000, vol. 133, no. 3, pp. 209–223.
- Vanag V.K. Study of spatially extended dynamical systems using probabilistic cellular automata. Physics-Uspekhi, 1999, vol. 42, no. 5, p. 413.
- Provata A., Nicolis G., Baras F. Oscillatory dynamics in low-dimensional supports: A lattice Lotka-Volterra model. J of Chem. Phys. 1999, vol. 110, pp. 8361–8368.
- Shabunin A., Baras F., Provata A. Oscillatory reactive dynamics on surfaces: a lattice limit cycle model. Physical Review E, 2002, vol. 66, no. 3. 036219.
- Tsekouras G., Provata A., Baras F. Waves and their interactions in the lattice Lotka–Volterra model. Izvestiya VUZ. Applied Nonlinear Dynamics, 2003, vol. 11, no. 2, pp. 63–71.
- Wood K., Van den Broeck C., Kawai R., Lindenberg K. Universality of synchrony: Critical behavior in a discrete of stochastic phase-coupled oscillators model. Physical Review Letters, 2006, vol. 96, no. 14, 145701.
- Efimov A., Shabunin A., Provata A. Synchronization of stochastic oscillations due to long-range diffusion. Physical Review E, 2008, vol. 78, no. 5, 056201.
- Verhulst P. Notice sur la loi que la population suit dans son accroissement. Corr. Math. et Phys. 1838, vol. 10, pp. 113–121.
- Allee W. C. Cooperation Among Animals with Human Implications. New York: Henry Schuman, Publisher, 1951.
- Neverova G.P., Khlebopros R.G., Frisman E.Y. The allee effect in population dynamics with a seasonal reproduction pattern. Biophysics, 2017, vol. 62, no. 6, pp. 967–976.
- Glendinning P., Perry L. P. Melnikov analysis of chaos in a simple epidemiological model. J. Math. Biol. 1997, vol. 35, pp. 359–373.

- 1815 reads