Covid-19: From Limit Cycle to Stable Focus

We investigate a new fundamental property of infectious diseases with natural adaptive immunity that weakens over time. Numerical experiments with a model of the Covid-19 epidemic in Moscow have demonstrated that when the reproduction number R0 is about 4, a qualitative change (bifurcation) occurs in the behavior of the virus-human system. Below this value, the long-term forecast tends to undamped oscillations; above it, the forecast shows damped oscillations: the amplitudes of epidemic waves decrease gradually, with a constant, very high background level of morbidity that keeps the natural immunity near 100%. To confirm this result analytically, we use an original modification of the Euler-Lotka renewal equation, which describes the dynamics of infected patients distributed by disease duration (time since infection) and accounts for immunity. To construct a bifurcation diagram, which illustrates the dependence of the equilibrium stability on the parameter R0, we linearize the equation in the vicinity of the equilibrium point and examine its numerical approximation (discrete form). This approximation can be interpreted as a Leslie model, with the matrix elements dependent on the parameter R0. By examining the roots of the corresponding Lotka polynomial, we can assess the stability of the equilibrium point. For the bifurcation diagram, we use the functions obtained from the simulation of the Covid-19 epidemic in Moscow. However, observations of the epidemic in other cities and countries support the primary finding of our study regarding the attenuation of epidemic waves.