On an SEIR Epidemic Model with Vaccination of Newborns and Periodic Impulsive Vaccination with Eventual On-Line Adapted Vaccination Strategies to the Varying Levels of the Susceptible Subpopulation

This paper investigates a susceptible-exposed-infectious-recovered (SEIR) epidemic model with demography under two vaccination effort strategies. Firstly, the model is investigated under vaccination of newborns, which is fact in a direct action on the recruitment level of the model. Secondly, it is investigated under a periodic impulsive vaccination on the susceptible in the sense that the vaccination impulses are concentrated in practice in very short time intervals around a set of impulsive time instants subject to constant inter-vaccination periods. Both strategies can be adapted, if desired, to the time-varying levels of susceptible in the sense that the control efforts be increased as those susceptible levels increase. The model is discussed in terms of suitable properties like the positivity of the solutions, the existence and allocation of equilibrium points, and stability concerns related to the values of the basic reproduction number. It is proven that the basic reproduction number lies below unity, so that the disease-free equilibrium point is asymptotically stable for larger values of the disease transmission rates under vaccination controls compared to the case of absence of vaccination. It is also proven that the endemic equilibrium point is not reachable if the disease-free one is stable and that the disease-free equilibrium point is unstable if the reproduction number exceeds unity while the endemic equilibrium point is stable. Several numerical results are investigated for both vaccination rules with the option of adapting through ime the corresponding efforts to the levels of susceptibility. Such simulation examples are performed under parameterizations related to the current SARS-COVID 19 pandemic.


Introduction
Epidemiology-from the Greek epi (about), demos (people), and logos (study)-tries to study the distribution and determinant factors of diseases and applies this study to their control and prevention. It integrates procedures and techniques from disciplines such as the biomedical sciences and the social sciences [1][2][3][4]. One of the first works on the subject was proposed by Kermack and McKendrick in [2] where an SIR (susceptible-infectious-recovered) epidemic model with the host population being split into three categories depending on the status of the individuals with respect to the disease was proposed. Multiple studies and investigations have been carried out, thus developing new variants of this model, such as SEIR (susceptible-exposed-infectious-recovered) or SIS (susceptible-infectious-susceptible) [3,4]. The use of mathematical models, epidemiological models, constitutes an essential part of epidemiology, mortality rate. Furthermore, Floquet´s theory on stability of periodic solutions is involved in the proposed analysis without an explicit calculation of the Floquet´s eigenvalues which conform some of the given conditions for the stability of the periodic oscillations. It is also proven mathematically in the paper that, in some cases, the local asymptotic stability ensures the global one, and that the periodic impulsive vaccination is slightly more efficient than the continuous traditional strategy.
In this paper, the application of both the newborns and impulsive vaccination strategies in the SEIR model are studied. We reconsider the problem of an SEIR model under an additional specific disease mortality rate under vaccination of newborns, which can be interpreted as an ad hoc modification of the recruitment in the susceptible subpopulation. Such a vaccination effort can be optionally accommodated to the current amounts of susceptible subpopulation. The basic objective is to study the stability of the equilibrium solutions of the system-both in the case of newborn vaccination and in the case of impulsive vaccination-for the proposed model, which helps to understand their long-term behavior. Concerning the newborn vaccination, it is demonstrated, both analytically and numerically, that the proposed vaccination strategy can lead to the eradication of an endemic disease. The disease-free and the endemic equilibrium points are both characterized explicitly, their local asymptotic stability properties are formally discussed based on the calculation of the next generation matrix and its eigenvalues, and the reproduction number is explicitly linked to the vaccination effort. Furthermore, the instability (respectively, the stability) of the disease-free equilibrium point is seen to coincide with the attainability and stability (respectively, the non-attainability) of the endemic one. The local asymptotic stability is seen to also be global by combining properties of non-negativity of the solution with the values of Poincaré index for individual or groups of singular related to the discussion on the existence of limit cycles. Concretely, the local asymptotic stability is proven to also be global on the basis of the fact that the stable attainable attractor is the unique disease-free equilibrium point if reproduction number is less than unity (in which case the endemic equilibrium point is not attainable). Contrarily, if the disease-free equilibrium point is unstable, in the case of a reproduction number exceeding unity, then it is shown that it cannot be surrounded by a limit cycle in any of the phase planes; the only stable and attainable attractor being then the endemic equilibrium point. Another alternative vaccination strategy, which is also discussed in the paper under the same basic model, consists of a periodic impulsive vaccination. Again, the vaccination effort can be accommodated to the current stocks of susceptible subpopulation. The stability properties of the periodic steady-state are performed via Floquet´s theory, with an explicit calculation of the eigenvalues of the monodromy matrix which are the Floquet multipliers. This method allows us to establish that, under a sufficiently small transmission rate, the local asymptotic stability is guaranteed according to a formula which makes explicit a trade-off between the vaccination interval between consecutive vaccination pulses and the susceptible population fraction to be vaccinated by a defining a reproduction-like number. To perform the above summarized study, this work is split into five sections which are organized as follows. Section 2 aims at establishing and discussing the SEIR model under eventual vaccination of newborns. The general formulation of the proposed SEIR model is given. Later on, the mathematical analysis of the model properties is discussed in detail in terms of local and global asymptotic stability concerns of both the disease-free and endemic equilibrium points, which are explicitly calculated as dependent of the vaccination of newborns levels. It is proven that the endemic equilibrium point is not reachable when the disease-free one is locally asymptotically stable, which happens as the basic reproduction number is less than unity. It is also proven that when the basic reproduction number exceeds unity then the endemic equilibrium point is reachable and locally asymptotically stable while the disease-free one is unstable. The reproduction number is improved with the levels of vaccination, which means that the reproduction number is kept below unity for wider ranges of potential values of the transmission rate as the vaccination effort increases. It is also proven that the local asymptotic stability of the equilibrium points also imply both respective global asymptotic stability properties. In the formulation of Section 2, the absence of vaccination of newborns might be understood as a particular case of the corresponding vaccination effort. On the other hand, in Section 3, an impulsive periodic vaccination strategy is proposed and studied, and the resulting model properties are formulated for the delay-free case. In this case, it is necessary to work with a resulting dynamic system of periodic coefficients, and, for this purpose, the Floquet theory is used to formulate the stability properties. Some numerical results are displayed and discussed in Section 4 for the COVID-19 pandemic, which are then performed to visualize the solutions for the two studied vaccination strategies. First, an endemic disease is simulated using the SEIR model, and then the mentioned vaccination strategies which are applied with the aim of eradicating the disease or, at least, softening its impact. To complete the section, two new control methods are designed by modifying the strategies presented by adapting them to the recruitment or susceptible levels in the sense that the control efforts are reinforced through time as such levels become larger. In this way, an attempt is made to propose new strategies that achieve eradication of the disease with less vaccination coverage, which is usually interesting in the context of the logistics of vaccination campaigns. The code for the programs developed for the simulations in this chapter has not been included in this report but can be found in the malenetxeberria/TFG-IngElec repository on GitHub. Finally, conclusions end the paper.
Compared to some related previous work, note that, in [10], the pulse vaccination strategy is programmed on an SIR epidemic model while in the current paper, the impulsive periodic vaccination is proposed for an SEIR epidemic model. There are recent studies for COVID-19 concerns, see, for instance [20][21][22][23][24][25]36]. However, those studies do not consider the potential application of vaccination which really does not yet exist in good condition for its mass application. Furthermore, most of the studies refer to short-term evaluations and prediction of the disease behavior. In the current work, two vaccination strategies are proposed for parameterizations associated to the COVID-19 pandemic. In particular, the second proposed strategy of periodic impulsive vaccination may be appropriate to deal with successive waves of serious incidence and later decrease of the infection. Furthermore, the vaccine administration strategies can be adapted to the levels of susceptible individuals which can facilitate the administration calendar.

The Epidemic Model
An SEIR epidemic model is proposed where the host population is split into four categories, namely: -Susceptible: individuals who have not contracted the disease but are at risk for it, since they do not have the necessary antibodies to cope with it. -Exposed: this is the group of individuals who are already infected, but who are not yet infectious, so they cannot transmit the infection. This state lasts for a limited period of time that varies according to the disease, and is known as the latent period [1] (Chap. Two). -Infectious: those individuals that are infected and can transmit the disease. -Recovered: individuals who have managed to defeat the disease. The body generates immunity as an adaptive response to the infectious agent [15], so that a recovered individual is no longer susceptible.
It is considered that the population is homogeneously mixed, that its spatial distribution is uniform, and that all susceptible individuals have the same probability of being infected, which simplifies the problem considerably. The mathematical representation of the SEIR model considered in this case [3] is given by the system of differential equations: . subject to an initial condition given by: S(0) ≥ 0, E(0) ≥ 0, I(0) ≥ 0 and R(0) ≥ 0. The variables S(t), E(t), I(t), and R(t) indicate the proportion of susceptible, exposed, infectious, and recovered individuals of the subpopulations, with respect to the total population, found in each compartment at a specific time.
The sum of the four subpopulations at the initial time instant is assumed to be unity, in general, so that the subpopulations of the model are fractions of the total initial population. However, this is not a real constraint since the initial conditions can be also fixed to the real values of the various subpopulations. Note the following facts on the above model: • New individuals enter the susceptible compartment per unit of time, either by birth or by immigration. βSI indicates the proportion of infections [4], that is, the proportion of individuals who are no longer susceptible and who become exposed; and finally, µS is the proportion of individuals who die from causes not related to the infection, also per unit of time.

•
The infected susceptible βSI pass into the exposed compartment and are removed from it once the latent period has passed, that is, σE individuals are removed per unit of time. Naturally, the fraction µE of those exposed individuals who die from causes unrelated to infection must also be considered.

•
At the end of the latent period, the exposed σE become infectious. Moreover, γI is the proportion of individuals who overcome the disease and who, therefore, pass to the recovered compartment; instead, αI represents the individuals who die from the infection, and who are removed from the population. • γI individuals enter the recovered compartment per unit of time, and µR are withdrawn due to mortality from causes not related to the infection.

•
The fraction of vaccinated newborns is directly added to the recovered compartment and withdrawn from the susceptible one.
In this way, the scheme of the SEIR model is represented in Figure 1.  Graphic representation of the susceptible-exposed-infectious-recovered (SEIR) model.

Non-Negativity of the Solution
The solutions of Equation (1) are non-negative for all time for any given non-negative initial condition. This is easily proven by calculating, in a closed analytic form, the unique solution of Equation (1) for any set of non-negative initial conditions. This property is needed for well-posedness of the epidemic model Equation (1). In parallel, the fact that the total population is proven to be bounded for all time guarantees also the boundedness for all time of all the subpopulations under any given arbitrary finite non-negative initial conditions. Theorem 1. The following properties hold for any q ∈ [0 , 1]: (i) The solution trajectory of the epidemic model (1) is non-negative for all time for any given non-negative initial conditions. (ii) The solution trajectory of Equation (1) is bounded for all time for any given finite non-negative initial conditions and lim sup t→∞ I(t) ≤ Amin(1/α , 1/µ). As a result, the epidemic model is globally stable for any given arbitrary finite non-negative initial conditions.
Proof. The solution to (1) is It is obvious that S(t) ≥ 0 for all time if S(0) ≥ 0 from the above first equation. Then, E(t) ≥ 0 for all time from the second equality of the above second equation if S(0) ≥ 0, E(0) ≥ 0, and I(0) ≥ 0 since S(t) ≥ 0 for all time. Subsequently, I(t) ≥ 0 for all time from the above third equation if I(0) ≥ 0, since E(t) ≥ 0 for all time. Finally, R(t) ≥ 0 for all time from the above fourth equation if R(0) ≥ 0, since I(t) ≥ 0 for all time. Property (i) has been proved. By summing-up the four subpopulations in Equation (1), one gets that the total population N(t) = S(t) + E(t) + I(t) + R(t) satisfies the differential Equation (3) .
whose solution is

Equilibrium Points
There are two equilibrium points, that is, the disease-free one, which is obtained from (1) by equalizing to zero all the time-derivatives of Equation (1) and equalizing the infectious to zero. This leads to On the other hand, the endemic equilibrium point is obtained from Equation (1) by equalizing again to zero all the time-derivatives while keeping a non-zero infectious subpopulation which yields: where: Recombining the above equations yields: From Equation (5) into Equation (3) equalized to zero yields the total endemic equilibrium population: Note that all the subpopulations of Equation (5), see also Equation (6), have to be non-negative for the endemic equilibrium point to exist in the first orthant of the state-space so that it be reachable. Otherwise, Theorem 1 would be violated. Thus, we have the subsequent result: Proposition 1. The endemic equilibrium point is reachable if and only if the transmission rate is large enough to be non less than a critical value according to: for any q ∈ [0 , 1). If the transmission rate equalizes its critical value, then the endemic equilibrium point is coincident with the disease-free one.

Basic Reproduction Number and Stability of the Equilibrium Points
One of the most important parameters in epidemiology is the basic reproduction number R 0 , which is defined as the average number of secondary cases caused by the entry of an infectious individual into a totally susceptible population [1] (Chap. 2). There are several techniques for calculating the basic reproduction number, and in this work, the next generation method will be shown [16], where R 0 is the spectral radius of the next generation matrix. So, to build this matrix, first, one defines the transmission and transition function matrices F and -V defined from the second and third equations of (1) by: which are evaluated at the disease-free equilibrium point to yield the matrix of dynamics of the linearized system around the disease-free equilibrium point, that is, the Jacobian matrix of the infective subsystem (E, I) at such a point, Ref. [38], where I 2 is the identity matrix of second order. Since −V =Ĵ EQ 0 β=0 is a stability matrix since the linearized system around the disease-free equilibrium point is globally asymptotically stable in the absence of disease, that is, if the transmission rate is zero, one follows from Equation (9) thatĴ EQ 0 is a stability matrix if and only if the spectral radius of V −1 F is less than one or, equivalently, that of that is referred to as the basic reproduction number, is less than one, that is, if and only if β < β c (q) and this Jacobian matrix is critically stable for β = β c (q) and unstable for β > β c (q). Comparing this result with Proposition 1, and since the eigenvalues of the Jacobian matrix Equation (9) of the infective pair (E, I) are continuous functions of the transmission rate, yields the result that the disease-free equilibrium point is asymptotically stable if and only if the endemic one is unreachable, which happens for a sufficiently small transmission rate as follows: The disease-free equilibrium point is the only reachable equilibrium point and it is globally asymptotically stable if and only if the transmission rate is small enough to satisfy β ∈ [0 , β c (q)). It is coincident with the endemic equilibrium point and critically stable if β = β c (q) and unstable if β > β c (q).
Proof. One has from the above discussion that for sufficiently close to the disease-free equilibrium point, E(t) → 0 and I(t) → 0 as t → ∞ if β < β c (q). On the other hand, note from Proposition 1 that, if β < β c (q), then the exposed subpopulation is negative which contradicts Theorem 1(i). Therefore, the endemic equilibrium is not attainable if the reproduction number is less than one so that the only global stable attractor is the disease-free equilibrium point. On the other hand, all the subpopulations are bounded for all time and any β ≥ 0 from Theorem 1. So that: (a) either S(t) → (1 − q)A/µ , R(t) → qA/µ as t → ∞ and the result is already proven; (b) at least one of S(t) and R(t) is unbounded as time increases, which is impossible since this unboundedness would contradict Theorem 1; (c) at least one of S(t) and R(t) exhibits asymptotic oscillations. However, then from the solutions in the proof of Theorem 1 and since I(t) ≥ 0 for all and t ≥ 0 as t → ∞ : therefore, S(t) and R(t) cannot oscillate asymptotically so that they converge to finite non-negative limits. Since the disease-free equilibrium point is unique, then S(t) → S 0 (q) = (1 − q)A/µ and R(t) → R 0 (q) = qA/µ as t → ∞ , and then the disease-free equilibrium point is (not only locally but also globally) asymptotically stable if and only if β ∈ [0 , β c (q)). If β > β c (q), then the spectral radius of the subsystem of the infective variables (E, I) exceeds unity so that the corresponding sub-matrix of the Jacobian matrix is unstable and then the whole Jacobian matrix at the disease-free equilibrium point is unstable. If β = β c (q), then some eigenvalue of the Jacobian matrix of the disease-free equilibrium point is at the stability boundary and the linearized system is critically stable.
Alternative Proof. Consider the whole Jacobian matrix at the disease-free equilibrium point: Note that λ 3 and λ 4 (< λ 3 ) are real eigenvalues of the above Jacobian matrix. Then, the Jacobian matrix at the disease-free equilibrium point is a stability matrix if and only if λ 3 < 0, that is, if and only if which holds if and only if the transmission rate is small enough to satisfy β < β c (q) or, equivalently, if the basic reproduction number Equation (11) is less than one, equivalently R 0 (q) < 1. If the inequality Equation (14), equivalently R 0 (q) < 1, becomes an equality, still equivalent to R 0 (q) = 1 in (11), then λ 3 = 0, that is, one eigenvalue of Equation (12) is critically stable and the disease-free equilibrium point is then critically stable. If the inequality in Equation (14) is reversed (R 0 (q) > 1) then the disease-free equilibrium point is unstable. Finally, note that the local asymptotic stability of the disease-free equilibrium point for R 0 (q) < 1 is also global since: (a) it is the unique attainable equilibrium point since the endemic one is not reachable from Proposition 1.
(b) the first right-hand-side term of Equation (2) tends asymptotically to zero as time tends to infinity so that, irrespective of the initial conditions, either the following limit exists: which implies that the susceptible subpopulation tends to a finite and non-negative limit (from Theorem 1) as time tends to infinity, which is the disease-free equilibrium susceptible subpopulation. In this case, the (non-negative) integrand of the second term of the above limit tends asymptotically to zero; or (c) it converges to a periodic solution of some period T > 0, that is, If the image of ρ : [0 , T) → R 0+ is a real constant, then the convergence is to a point for all ρ ∈ [0 , T) which is the susceptible component of the disease-free equilibrium point which is already discussed. If ρ : [0 , T) → R 0+ is not constant, then the susceptible subpopulation converges to a periodic solution. However, note that the second right-hand-side term of the first equation of (2) has a non-negative integrand so that either it converges to zero in order for S(t) to be bounded for all time or S(t) is asymptotically unbounded which violates Theorem 1. As a result, S(t) cannot exhibit an asymptotic oscillation while it converges to the disease-free susceptible equilibrium point as time tends to infinity. Now, consider the two-dimensional phase plane .
at the disease-free susceptible equilibrium point S * (q). Since the reproduction number is less than one, the equilibrium point cannot be a saddle point, then its Poincaré index is +1 and no stable limit cycle can surround it. Thus, the disease-free equilibrium point, which is a local asymptotic attractor and the unique stable global attractor, is also the unique global asymptotic attractor if the reproduction number is less than one. Remark 1. Note that, for any q ∈ [0 , 1), the basic reproduction number in the case of vaccination of newborns becomes reduced to R 0 (q) = (1 − q)R 0 (0) compared to the case of absence of vaccination (whose corresponding reproduction number is R 0 (0)), and then the disease-free equilibrium point is globally asymptotically stable if and only if R 0 (q) < 1. This implies that, under a vaccination of newborns effort, the global asymptotic stability of the disease-free equilibrium point together with the non-reachability of the endemic one are achieved if and only if R 0 (0) ∈ [0, 1/(1 − q)).

Remark 2.
An equivalent interpretation is that the stability of the disease-free equilibrium point and the parallel property of non-reachability of the endemic one are achievable for a larger range of transmission rate related to the absence of vaccination of newborns since the critical transmission rate satisfies β c (q) > β c (0) if q ∈ [0 , 1) (Proposition 1).

Remark 3.
As a result of the above considerations, one can argue that vaccination of newborns effort is useful in achieving the asymptotic eradication of the disease since the range of transmission rate values being compatible with the stability of the disease-free equilibrium increases related to its parallel range of values in the absence of vaccination.
The subsequent result, whose proofs follow by direct calculation, relies on alternative equivalent expressions to Equation (6) on the components of the endemic equilibrium point depending on the reproduction number.

Proposition 2.
The components of the endemic equilibrium point may be rewritten, equivalently to (6), as follows: It is reachable if and only if R 0 (q) ≥ 1 and it is coincident with the disease-free equilibrium point if and only if R 0 (q) = 1. Proof. The Jacobian matrix at the endemic equilibrium point is: One of the eigenvalues of (16) is λ 1 = −µ. Instead of calculating the other three eigenvalues of Equation (16), which is a difficult task, the Routh-Hurwitz stability criterion p(λ) = 3 i=0 a 3−i λ i = 0 is applied to the characteristic equation of the third order remaining submatrix [31], where The Routh-Hurwitz criterion is applied by first calculating the following determinants: The number of characteristic roots of p(λ) = 0 with positive real part is the number of sign changes in the set {∆ 0 , ∆ 1 , ∆ 2 /∆ 1 , ∆ 3 /∆ 2 }. Since ∆ 0 > 0 and ∆ 1 > 0, then there is no characteristic root with positive real part if and only if ∆ 2 /∆ 1 > 0 , ∆ 3 /∆ 2 > 0 [28]. Note that ∆ 3 /∆ 2 > 0 if and only if R 0 (q) > 1, and ∆ 2 /∆ 1 > 0 if and only if which always holds. Then, (16) is a stability matrix if and only if R 0 (q) > 1. Thus, the endemic equilibrium point is locally asymptotically stable if and only if R 0 (q) > 1. It is now proven that such a stability property is global by using contradiction arguments. Assume that the solution is oscillatory periodic with some arbitrary oscillation period T > 0. The period is supposed arbitrary since it is being rebutted by contradiction arguments in the following that some such oscillation period exits so that no limit oscillation exists as a result. Defineŷ(δ) = lim  (2) in the proof of Theorem 1, one gets for t = nT + δ, n ∈ Z 0+ = {z ∈ Z : z ≥ 0}, and δ ∈ [0 , T) that the limits below exist for any such a δ ∈ [0 , T): after using the mean value theorem for integrals since all the involved integrals in Equation (2) are continuous functions of time. Note thatŷ : [0, T) → R 0+ will define the periodic limit oscillation of the subpopulation z in the event that such an oscillation exists. By recombining the three first above equations, one gets for any given q ∈ (0 , 1], any δ ∈ (0 , T], and any supposed arbitrary oscillation period T > 0 that, sinceÊ(τ) =Ê M (δ) for some τ ∈ (0 , δ] and any δ ∈ (0 , T]: Hence, a contradiction follows such that the exposed subpopulation does not have a limit oscillation of any potential oscillation period, since a chaotic behavior cannot exist either since the model has no small time-varying parameter. Similarly, one gets a similar conclusion on no existence of limit oscillations of any period of the infectious subpopulation since the oscillation condition is violated for a range of sufficiently small inter-period time instants for any arbitrary claimed oscillation period according to:Î and similar conclusions can be deduced for the susceptible and the recovered individuals. Therefore, no limit cyclic exists surrounding the reachable endemic equilibrium point when the disease-free equilibrium point is unstable, that is, if the basic reproduction number exceeds unity. Furthermore, all solutions are bounded from Theorem 1. As a result, the endemic equilibrium point is globally asymptotically stable if the basic reproduction number exceeds unity, equivalently, if and only if β ∈ (β c (q) , ∞).

Remark 4.
It can be of interest to give some considerations which are based on known general results of the qualitative theory of differential equations which indicate that there is a unique global asymptotically stable attractor which is always one of the equilibrium points [43][44][45]. The global attractor is the disease-free equilibrium if the reproduction number is less than or equal to one and the endemic equilibrium point if such a parameter exceeds unity. In this way, there are no stable oscillations surrounding neither one of the equilibrium points nor jointly both of them. If the reproduction number is unity, then the asymptotic behavior is similar to that arising if it is less than unity so that the disease-free equilibrium point is a global asymptotically stable attractor. Thus, the constraints implying local asymptotic stability of any of the equilibrium points lead to their global asymptotic stability. Those considerations are deduced from the subsequent linked arguments: (a) If R 0 (q) < 1, then the disease-free equilibrium point is locally asymptotically stable and the endemic one is not reachable. A limit cycle surrounding it in any phase plane of two subpopulations, if any, would be unstable from the combined topology of equilibrium points and limit cycles. Therefore, it would be irrelevant for stability considerations. As a result, the disease-free equilibrium point is globally asymptotically stable (Theorem 2). (b) If R 0 (q) > 1, then the endemic equilibrium point is locally asymptotically stable and the disease-free equilibrium point is unstable. A limit cycle surrounding it leaving away the disease-free one, if any, would be unstable for stability considerations. No limit cycles can exist surrounding both equilibrium points since the combined Poincaré index of both of them would be 0 (if the disease-free equilibrium point were a saddle point) or +2 (if the disease-free equilibrium point were not a saddle point). So, in any case, such a combined Poincaré index is not +1 compatible with the existence of some limit cycle surrounding both equilibrium points. As a result, the endemic equilibrium point is globally asymptotically stable if R 0 (q) > 1 (Theorem 3). (c) If R 0 (q) = 1, then the disease-free and the endemic equilibrium points are coincident. From the non-negativity property of the solution for any given non-negative initial conditions (Theorem 1 (i)), it follows that a limit cycle surrounding the equilibrium point cannot exist in the phase plane (E, I). Also, by inspection of (1), it follows that S(t) and R(t) can only have periodic asymptotic solutions if I(t) and E(t) have asymptotic periodic solutions. As a result, the stability is also asymptotic and global in this case.

The Epidemic Model Under Periodic Impulsive Vaccination
Consider a modification of Equation (1) as follows with q = 0 (no vaccination of newborns) and under a periodic impulsive vaccination of gain p ∈ [0 , 1] and period T at the time instants t = nT for n ∈ Z 0+ : .
for all t ∈ [nT , (n + 1)T); ∀n ∈ Z 0+ with initial conditions S(0) ≥ 0, E(0) ≥ 0, I(0) ≥ 0, and R(0) ≥ 0 for n = 0, where δ(t) is the Dirac distribution which reflects in the model dynamics the presence of a very intensive vaccination (ideally of impulsive infinity gain) concentrated in very short time intervals (ideally instantaneously administered at the time instants t n = nT). The notation S(nT − ) stands for the left limit of S(t) at t n = nT while the right limit (i.e., the value as the impulse at t n is accounted for) is simply denoted as S(nT) = S(nT + ) = (1 − p)S(nT − ). The removed susceptible individuals at the impulsive time instants become recovered, which is reflected by the solution of the above fourth equation at the impulsive time instants, that is, R(nT) = R(nT + ) = R(nT − ) + pS(nT − ). Note that p = 0 describes the absence of impulsive vaccination while p = 1 describes the total vaccination of the susceptible at the time instants t n = nT, since no susceptible individual could be contacted by an infectious one because of the total removal of the first ones from the environment under study. It can be pointed out that impulsive vaccination is equivalent to intervention measures like isolation or confinement since a number of susceptible are transferred to the recovered subpopulation. The above-mentioned equivalence is at the level of transitions between subpopulations although no immunity is expected from such a state transition from susceptible to recovered. In this sense, the total recovered subpopulation is understood in a very wide sense as those individuals who are not susceptible to being infected, rather than just as those which are in fact immune, while those being immune are also included in the recovered subpopulation. Note that Equation (3) is still valid for the dynamics of the total population. That means that there are no finite jumps of the total population at the impulsive time instants (in particular, N * = N − * = A/µ at the disease-free equilibrium point) and there are no jumps either of the exposed or the infectious subpopulations at the disease-free equilibrium point.
By equalizing to zero the time-derivatives of the four subpopulations in Equation (21) while zeroing the exposed and infectious values, one gets the subsequent disease-free steady-state: with E * = E − * = I * = I − * = 0.

Remark 5.
Note that if q = p = 0 (i.e., in the absence of both newborns and impulsive vaccination efforts) then the disease-free equilibrium point discussed in the former section coincides with the above disease-free steady state. Note also that Equations (22) and (23) imply that the disease-free steady state is a periodic oscillation rather than an equilibrium point for nonzero p.
Then, any solution of Equation (21) for any given non-negative initial conditions is non-negative for all time. Property (i) has been proven. To prove Property (ii), first equalize to zero all the time-derivatives in Equation (21) to get E(t) = I(t) = 0, and that the solution of the first differential equation is: In particular, the above equation leads to at the impulsive time instants. Denoting S n = S(nT) so that the stroboscopic mapping F : S n−1 → S n ; ∀n ∈ Z + defined by with c = A/µ, has a unique fixed point satisfying S * = F(S * ) defined by the first relationship of Equation (22). In this way, the delay-free periodic solution of the susceptible, whose right and left limit values at the impulsive time instants are the two equations of Equation (22) is: The periodic disease-free solution of the recovered is given byR(t) =N(t) −Ŝ(t), whereN(t) is the periodic solution of the total population which satisfies under zero infection the differential equation: .
Property (ii) has been proven. Property (iii) is a direct consequence of Equations (22) and (23) for the disease-free solution. To prove Property (iv), the periodic oscillation is firstly perturbed with a small perturbation (s(t) , e(t) , i(t) , r(t)), that is, which results in the following modification of the differential system Equation (21) after neglecting the contribution to perturbation terms of orders higher than one: .
The periodic solution is stable if the perturbation vanishes through time. The functionŜ(t) is given by Equation (31) and it is a periodic coefficient of period T in Equation (35). The last equation is not needed for analysis if one considers the relationshipsR(t) =N(t) −Ŝ(t) = A/µ −Ŝ(t) for the disease-free solution. Therefore, a third-order differential system being built with the first three relations of Equation (35) might be compacted into a third-order differential system of the form .
x(t) = A(t)x(t) with A(t) = A(t + T) being a third-order periodic matrix function of dynamics for all t ≥ 0. Its entries follow directly from inspection of Equation (35). Such a system has to be solved on a time interval [0 , T]. A fundamental matrix X(t) = X(t + T) of the third-order differential system, that is being nonsingular for all time since its columns are particular linearly independent solutions to initial conditions satisfying X(0) = I 3 , is defined for all t = nT as X(T), the monodromy matrix, being of the form: where The eigenvalues of the monodromy matrix, that is, its diagonal entries in this case are the Floquet multipliers [29,30], and the matrix A(t) is a stability matrix if the moduli of their Floquet multipliers are less than one. One gets by inspection of Equations (36)-(38) that the two first eigenvalues have a modulus less than one since µ > 0, 0 ≤ p ≤ 1 and P + (t) > 0; ∀t ∈ [0 , T] and the third eigenvalue has a modulus less than one if and only if The condition Equation (39) implies that the periodic solution is locally asymptotically stable in the sense that the solution of the incremental system Equation (35) converges asymptotically to zero for sufficiently small initial conditions. Property (iv) has been proven. To prove Property (v), note that a sufficient condition for Equation (39) to hold is that which holds if, equivalently, Equation (26) holds. Property (v) has been proven. Remark 6. Theorem 4 (v) gives a testable trade-off between the sizes of the critical transmission rate β c and those of the impulsive period (p and T) to guarantee the local stability of the periodic solution. In particular, for small periods for the impulsive vaccination actions, the critical transmission rate value compatible with the local asymptotic stability of the disease-free solution is high. That means that, in the main range of values of interest, the tolerance to the transmission rate guarantees that such a stability tolerance increases as the impulsive vaccination period decreases as intuitively expected. Note that it is not difficult to extend the performed study to the combination of the two combined vaccination strategies of vaccination of newborns and periodic impulsive vaccination. On the other hand, the generalization to impulsive vaccination controls under time-varying inter-impulse periods are an open design tool for related future work based on the adaptive sampling rate designs [43,44].

Numerical Simulations
The SEIR model is now evaluated under numerical simulations. The system of differential equations has been implemented in Python based on Odeint of the SciPy package which solves ordinary differential equations. The epidemiological parameters were taken from COVID-19 pandemic registered data in the background literature which result to be β = 0.57, σ = 0.13, γ = 0.067 and α = 0.01 in day −1 [6][7][8].

Newborn Vaccination
An initial total population of N(0) = 4 × 10 7 individuals with the following demographic parameters: A = 0.005 (indivual × day) −1 and µ = 0.005 day −1 was assumed. The model was evaluated under normalization to unity. The initial conditions of the subpopulations were S(0) = 1 − 1/N(0), E(0) = R(0) = 0, I(0) = 1/N(0). The evolution in the absence of vaccination is displayed in Figure 2. The basic reproduction number was R 0 = 6.7, and the solution converges to the endemic equilibrium point which is asymptotically stable with an infectious proportion of I * = 0.05. Figure 3 exhibits the evolution of the infectious for different levels of vaccination effort. Note that for q > q c = 0.85, the disease-free equilibrium point again becomes a stable attractor.

Newborn Vaccination
An initial total population of N 0 ( ) = 4 ×10 7 individuals with the following demographic parameters: A = 0.005 (indivual × day) −1 and μ = 0.005 day -1 was assumed. The model was evaluated under normalization to unity. The initial conditions of the subpopulations were S 0 The evolution in the absence of vaccination is displayed in Figure 2.
The basic reproduction number was R 0 = 6.7 , and the solution converges to the endemic equilibrium point which is asymptotically stable with an infectious proportion of I * = 0.05 . Figure 3 exhibits the evolution of the infectious for different levels of vaccination effort. Note that for q > q c = 0.85, the disease-free equilibrium point again becomes a stable attractor.  It is now evaluated how to update the vaccination effort to the proportion of susceptible individuals through time. The basic idea is to increase the effort as the susceptible individuals increase. For that purpose, the following updating formula is used for the vaccination effort which adapts the vaccination effort to the proportion of susceptible through time:  It is now evaluated how to update the vaccination effort to the proportion of susceptible individuals through time. The basic idea is to increase the effort as the susceptible individuals increase. For that purpose, the following updating formula is used for the vaccination effort which adapts the vaccination effort to the proportion of susceptible through time: q(S(t)) = (q max − q min )S(t) + q min (40) Figures 4 and 5 show the model evolution with q max = 0.90, and q min = 0.55 and q min = 0.775, respectively, under the newborns vaccination effort designed with the above formula. The endemic proportion of infectious is I * = 0.015 in the first case and I * = 0.0036 in the second one. Note that Figure 5 exhibits an improvement in the maximum infectious level at the expense of larger vaccination efforts according to the above formula. It is now evaluated how to update the vaccination effort to the proportion of susceptible individuals through time. The basic idea is to increase the effort as the susceptible individuals increase. For that purpose, the following updating formula is used for the vaccination effort which adapts the vaccination effort to the proportion of susceptible through time:

Impulsive Vaccination
Firstly, one considers an impulsive vaccination with parameters of vaccination proportion and impulse period p = 0.45 and T = 120days . The evolution through time is displayed in Figure 6. It is seen that the evolution evolves to a periodic solution as expected.

Impulsive Vaccination
Firstly, one considers an impulsive vaccination with parameters of vaccination proportion and impulse period p = 0.45 and T = 120 days. The evolution through time is displayed in Figure 6. It is seen that the evolution evolves to a periodic solution as expected. Figure 5. Evolution with q max = 0.90 and q min = 0.775.

Impulsive Vaccination
Firstly, one considers an impulsive vaccination with parameters of vaccination proportion and impulse period p = 0.45 and T = 120days . The evolution through time is displayed in Figure 6. It is seen that the evolution evolves to a periodic solution as expected.       The performance of increasing the vaccination effort as the susceptible individuals increase is now evaluated. For that purpose, the following updating formula is used for the vaccination effort which adapts the vaccination effort to the proportion of susceptible over time to improve the efficiency of the vaccination strategy:  The performance of increasing the vaccination effort as the susceptible individuals increase is now evaluated. For that purpose, the following updating formula is used for the vaccination effort which adapts the vaccination effort to the proportion of susceptible over time to improve the efficiency of the vaccination strategy: p(S(t)) = (p max − p min )S(t) + p min (41) Some tested results are shown in Figures 9 and 10 with A = 0.0045, p max = 0.60, p min = 0.30, and T = 30 days. The test is repeated in Figures 11 and 12 by increasing the minimum vaccination effort to p min = 0.55 which improves the infection control and so improves the infection evolution, as expected.

Conclusions
In this work, the implementation of newborns and impulsive vaccination in the epidemiological model SEIR was studied. For this purpose, it was necessary to begin with an extensive theoretical framework in which the generalities of the model and the two vaccination strategies were explained, delving into some important concepts such as the existence and location

Conclusions
In this work, the implementation of newborns and impulsive vaccination in the epidemiological model SEIR was studied. For this purpose, it was necessary to begin with an extensive theoretical framework in which the generalities of the model and the two vaccination strategies were explained, delving into some important concepts such as the existence and location

Conclusions
In this work, the implementation of newborns and impulsive vaccination in the epidemiological model SEIR was studied. For this purpose, it was necessary to begin with an extensive theoretical framework in which the generalities of the model and the two vaccination strategies were explained, delving into some important concepts such as the existence and location of the equilibrium points or the state-steady periodic solutions and their stability properties, as well as their links to the ranges of values of the basic reproductive number. Although the central topic of the work is to study the influence of vaccination, it was considered important to devote a part of the paper body to these topics since feedback vaccination influences the allocation of the equilibrium states and the rates of convergence related to the absence of vaccination. In order to further study the vaccination model and control strategies, a part of the article has been devoted to stability analysis. In the case of newborn vaccination, the disease-free and endemic equilibrium points were found and Lyapunov's first method was applied to verify their local stability properties. It was proven that the endemic equilibrium point is not reachable as the disease-free one is locally asymptotically stable while the endemic one is reachable as the disease-free one is unstable. Those facts also lead to the conclusion that the respective local stability properties are also global since any solution is bounded for finite non-negative initial conditions.
In the case of impulsive vaccination, the calculations were more involved due to the impulsive nature of the system. It can be argued that the impulsive vaccination strategy may be of particular relevant interest in endemic situations when the disease exhibits successive decreasing and increasing levels of infection through time. Instead of equilibrium points, periodic solutions are found and characterized and, due to the complexity of the problem, only the disease-free periodic solution was obtained. To demonstrate the local stability of the calculated solution, Floquet's theory was applied, and in this way a stability condition containing an integral was obtained. Through some performed numerical simulations, it was possible to conclude, firstly, that both newborn and impulsive vaccination are effective strategies to eradicate the disease under sufficiently large levels of vaccination efforts. The required proportion of vaccination effort was seen to be relevant in both cases if the proportion of vaccines administered is proportional to the proportion of susceptible individuals. Secondly, it has been observed that the total coverage necessary to eradicate the disease was lower in the case of impulsive vaccination than in newborn vaccination, as mentioned in the theoretical framework. Regarding the strategy of impulsive vaccination, the local stability of a disease-free periodic solution was also demonstrated. Numerical simulations were developed by using background provided data for the disease parameterization of the SARS-COVID-19 pandemic. Although at this time there are no fully tested vaccines available for such a disease, the numerical examples adjusted to the formal setting indicate that the vaccination strategy improves the basic reproduction number and the disease is controlled more efficiently. It was also proven and tested through examples that the global vaccination strategy is better administrated if the vaccination effort is proportional to the susceptible subpopulation through time.