Mathematical modeling and cellular automata simulation of infectious disease dynamics: Applications to the understanding of herd immunity.

The complexity associated with an epidemic defies any quantitatively reliable predictive theoretical scheme. Here, we pursue a generalized mathematical model and cellular automata simulations to study the dynamics of infectious diseases and apply it in the context of the COVID-19 spread. Our model is inspired by the theory of coupled chemical reactions to treat multiple parallel reaction pathways. We essentially ask the question: how hard could the time evolution toward the desired herd immunity (HI) be on the lives of people? We demonstrate that the answer to this question requires the study of two implicit functions, which are determined by several rate constants, which are time-dependent themselves. Implementation of different strategies to counter the spread of the disease requires a certain degree of a quantitative understanding of the time-dependence of the outcome. Here, we compartmentalize the susceptible population into two categories, (i) vulnerables and (ii) resilients (including asymptomatic carriers), and study the dynamical evolution of the disease progression. We obtain the relative fatality of these two sub-categories as a function of the percentages of the vulnerable and resilient population and the complex dependence on the rate of attainment of herd immunity. We attempt to study and quantify possible adverse effects of the progression rate of the epidemic on the recovery rates of vulnerables, in the course of attaining HI. We find the important result that slower attainment of the HI is relatively less fatal. However, slower progress toward HI could be complicated by many intervening factors.