Global stability of a SEIR discrete delay differential‐difference system with protection phase

Abstract

We consider an epidemiological model with the four classical compartments of susceptible, exposed, infected, and recovered population. We add a new compartment that is supposed to describe, for a limited time, individuals that are protected from the epidemic through vaccination or medication, for instance. We model the protection phase by an age‐structured partial differential equation. The age is the time since an individual entered the protection phase. The model is then reduced by integration on the characteristics to a differential‐difference system with delay. The discrete delay represents the limited duration of the protection phase. After establishing the basic properties of the model, we show that the disease‐free equilibrium (DFE) is globally asymptotically stable when the basic reproduction number is less than one and is unstable when this number is greater than one. Furthermore, we show that even if there is no mortality during the protection phase and the basic reproduction number is greater than one, the endemic equilibrium is globally asymptotically stable. The proofs of the global asymptotic stability of both equilibria are based on carefully constructed Lyapunov functions. To complete this study on the global dynamics, we discuss some results on weak and strong uniform persistence of the disease. Finally, numerical simulations are performed to illustrate and complete our main results.

Global Asymptotic Stability of a Hybrid Differential–Difference System Describing SIR and SIS Epidemic Models with a Protection Phase and a Nonlinear Force of Infection

“…(5). This work is a continuation of recently published articles [1][2][3]6]. These models are based on a coupling between differential equations and continuous-time difference equations.…”

“…In this section, we discuss the existence of steady-states and establish the basic reproduction number associated with the model (14). We assume that the function g is given by the expression (3) or that it satisfies the conditions (5). A steady-state Ē = ( S, Ī, ū) is a stationary solution of the system (14).…”

“…In this section, we are interested in the local asymptotic stability of the two steadystates of the hybrid system (14). We assume that the function g is given by the expression (3) or that it satisfies the conditions (5). We are now going to linearize the system ( 14) around the two steady-states.…”

“…The main result is the study of the global asymptotic stability of the two steady-states. In our previous articles, [1,2], we only considered incidence functions in the form of mass action (3).…”