Analytical Modeling of the Temporal Evolution of Epidemics Outbreaks Accounting for Vaccinations

With the vaccination against Covid-19 now available, how vaccination campaigns influence the mathematical modeling of epidemics is quantitatively explored. In this paper, the standard susceptible-infectious-recovered/removed (SIR) epidemic model is extended to a fourth compartment, V, of vaccinated...

Full description

Bibliographic Details
Main Authors: Reinhard Schlickeiser, Martin Kröger
Format: Article
Language:English
Published: MDPI AG 2021-05-01
Series:Physics
Subjects:
Online Access:https://www.mdpi.com/2624-8174/3/2/28
Description
Summary:With the vaccination against Covid-19 now available, how vaccination campaigns influence the mathematical modeling of epidemics is quantitatively explored. In this paper, the standard susceptible-infectious-recovered/removed (SIR) epidemic model is extended to a fourth compartment, V, of vaccinated persons. This extension involves the time <i>t</i>-dependent effective vaccination rate, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>v</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math></inline-formula>, that regulates the relationship between susceptible and vaccinated persons. The rate <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>v</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math></inline-formula> competes with the usual infection, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>a</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math></inline-formula>, and recovery, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>μ</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math></inline-formula>, rates in determining the time evolution of epidemics. The occurrence of a pandemic outburst with rising rates of new infections requires <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>k</mi><mo>+</mo><mi>b</mi><mo><</mo><mn>1</mn><mo>−</mo><mn>2</mn><mi>η</mi></mrow></semantics></math></inline-formula>, where <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>k</mi><mo>=</mo><mi>μ</mi><mo>(</mo><mn>0</mn><mo>)</mo><mo>/</mo><mi>a</mi><mo>(</mo><mn>0</mn><mo>)</mo></mrow></semantics></math></inline-formula> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>b</mi><mo>=</mo><mi>v</mi><mo>(</mo><mn>0</mn><mo>)</mo><mo>/</mo><mi>a</mi><mo>(</mo><mn>0</mn><mo>)</mo></mrow></semantics></math></inline-formula> denote the initial values for the ratios of the three rates, respectively, and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>η</mi><mo>≪</mo><mn>1</mn></mrow></semantics></math></inline-formula> is the initial fraction of infected persons. Exact analytical inverse solutions <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>t</mi><mo>(</mo><mi>Q</mi><mo>)</mo></mrow></semantics></math></inline-formula> for all relevant quantities <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>Q</mi><mo>=</mo><mo>[</mo><mi>S</mi><mo>,</mo><mi>I</mi><mo>,</mo><mi>R</mi><mo>,</mo><mi>V</mi><mo>]</mo></mrow></semantics></math></inline-formula> of the resulting SIRV model in terms of Lambert functions are derived for the semi-time case with time-independent ratios <i>k</i> and <i>b</i> between the recovery and vaccination rates to the infection rate, respectively. These inverse solutions can be approximated with high accuracy, yielding the explicit time-dependences <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>Q</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math></inline-formula> by inverting the Lambert functions. The values of the three parameters <i>k</i>, <i>b</i> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mi>η</mi></semantics></math></inline-formula> completely determine the reduced time evolution of the SIRV-quantities <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>Q</mi><mo>(</mo><mi>τ</mi><mo>)</mo></mrow></semantics></math></inline-formula>. The influence of vaccinations on the total cumulative number and the maximum rate of new infections in different countries is calculated by comparing with monitored real time Covid-19 data. The reduction in the final cumulative fraction of infected persons and in the maximum daily rate of new infections is quantitatively determined by using the actual pandemic parameters in different countries. Moreover, a new criterion is developed that decides on the occurrence of future Covid-19 waves in these countries. Apart from in Israel, this can happen in all countries considered.
ISSN:2624-8174