Next Article in Journal
The Spectral Distribution of Random Mixed Graphs
Next Article in Special Issue
Weighted Generalized Fractional Integration by Parts and the Euler–Lagrange Equation
Previous Article in Journal
On the Semi-Group Property of the Perpendicular Bisector in a Normed Space
Previous Article in Special Issue
Determining COVID-19 Dynamics Using Physics Informed Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Impact of the Imperfect Vaccination of the COVID-19 with Optimal Containment Strategy

1
Fundamental and Applied Mathematics Laboratory, Department of Mathematics and Computer Science, Faculty of Sciences Ain Chock, Hassan II University of Casablanca, Maarif B.P 5366, Casablanca 20100, Morocco
2
Department of Mathematical Sciences, United Arab Emirates University, Al Ain P.O. Box 15551, United Arab Emirates
3
Laboratory of Analysis Modelling and Simulation, Department of Mathematics and Computer Science, Faculty of Sciences Ben M’sik, Hassan II University of Casablanca, Sidi Othman Casablanca B.P 7955, Casablanca 20100, Morocco
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(3), 124; https://doi.org/10.3390/axioms11030124
Submission received: 16 January 2022 / Revised: 21 February 2022 / Accepted: 25 February 2022 / Published: 10 March 2022

Abstract

:
Since the beginning of the COVID-19 pandemic, vaccination has been the main strategy to contain the spread of the coronavirus. However, with the administration of many types of vaccines and the constant mutation of viruses, the issue of how effective these vaccines are in protecting the population is raised. This work aimed to present a mathematical model that investigates the imperfect vaccine and finds the additional measures needed to help reduce the burden of disease. We determine the R 0 threshold of disease spread and use stability analysis to determine the condition that will result in disease eradication. We also fitted our model to COVID-19 data from Morocco to estimate the parameters of the model. The sensitivity analysis of the basic reproduction number, with respect to the parameters of the model, is simulated for the four possible scenarios of the disease progress. Finally, we investigate the optimal containment measures that could be implemented with vaccination. To illustrate our results, we perform the numerical simulations of optimal control.

1. Introduction

Since the beginning of the ongoing COVID-19 pandemic, the world has been racing to develop a vaccination that helps protect the populations around the world and bring human life to a normal status. This race to find a vaccine was not only challenged by the fast spread of the disease, but also the high rate of mutation of COVID-19. As a result, we witness many vaccination types with different biotechnological approaches and different efficacy [1]. These efficacies are based on clinical trials that might have some limitations as their samples do not necessarily cover a wide population from different parts of the world. These facts make the question of the efficacy of vaccines legitimate and need to be investigated.
This problem was investigated using mathematical modeling to study the possible measure that needed to be implemented to reduce the impact of an imperfect vaccine.
The mathematical model of the imperfect vaccination of infectious disease was started by the work of Arino et al. [2]. Many papers have followed up on this work, looking into the various spreads of imperfect vaccination for various diseases. For example, the work of Abu-Raddad [3] studied the mathematical model of a possible HIV vaccination where the authors investigated the impact of the defectiveness of vaccination on the progress of the disease. Liu et al. [4] studied a general SIR with an added vaccination compartment and imperfect vaccination. This study showed how vaccination reduced the infected population but could not eradicate the infection. In fact, the eradication required an additional necessary condition. If vaccination efficacy improves, this condition may be alleviated. A mathematical model with the imperfect vaccination of birds in the case of avian influenza was studied in [5]. This model considered age-since-vaccination structure and symptomatically infected birds. This study showed that the only way eradicate the disease was by the full coverage of the bird population or by full efficacy. A time delay model of imperfect vaccination was studied in [6] with a possible loss of immunity. The study showed the existence of the critical vaccination coverage needed to eliminate the infection. In the case of imperfect vaccination, the authors showed that a critical proportion of the population needed vaccination. Another delay model with distributed delay [7] and the delay model with a generalized incidence function were studied in [8].
Regarding the ongoing pandemic, there are some studies that have investigated imperfect vaccination in the USA ([9,10]) and the UK [11] but without finding optimal measures that could help contain the pandemic, as the use of an imperfect vaccine cannot achieve the low endemicity of COVID-19. On the other hand, many studies (see [12,13,14,15,16,17,18,19,20,21,22]) used optimal control to find the optimal way to allocate vaccination and the best strategy to vaccinate the population, depending on the age or comorbidity of the population. The goal of this paper was to investigate a mathematical model of the imperfect vaccination of COVID-19. The aim was to study the dynamics of this model and present the possible control measures that need to be implemented in order to reduce the impact of the vaccine’s imperfection.
To our knowledge, our work is the only one to date to have studied the potential dynamics of imperfect vaccination and the optimal use of other public health measures that help to reduce the effect of administering imperfect vaccination. The only work that combined these two problems was used in the case of possible malaria vaccination [23].
The structure of this paper is summarized as follows. In Section 2, the mathematical model is formed and the existence conditions of the system are verified. Section 3 takes into account the basic reproduction number. Section 4 and Section 5 analyze the local and global stability at the disease-free equilibrium point, respectively. The optimal control problem of vaccination and additional measures to reduce the disease spread are presented in Section 6. In Section 7, we fit our model to data from Morocco to estimate the parameters of the model. We also discussed four possible scenarios of the dynamic of the model via the elasticity of the basic reproduction number and we give the illustration of the optimal solution via numerical simulations. The conclusion and discussion of the results are given in Section 8.

2. The Mathematical Model

Nine epidemiological compartments were recognized in the population: susceptible (S); vaccinated (V); exposed (E) (asymptomatic); infected with mild symptoms (I); quarantined in a home with mild symptoms (Q); hospitalized (H); quarantined in a hospital with complications (C) (i.e., isolated in a hospital with breathing assistance); and mortality due to disease (D):
d S d t = β 1 S E N β 2 S I N λ 1 S d V d t = λ 1 S β 3 V E N β 4 V I N λ 2 V d E d t = β 1 S E N + β 2 S I N + β 3 V E N + β 4 V I N θ E d I d t = θ E ( γ 1 + γ 2 + γ 3 ) I d Q d t = γ 1 I ( σ 1 + δ 1 ) Q d H d t = γ 2 I + σ 1 Q ( σ 2 + δ 2 ) H d C d t = γ 3 I + σ 2 H ( μ + δ 3 ) C d R d t = δ 1 Q + δ 2 H + δ 3 C + λ 2 V d D d t = μ C
In this model, we assume that vaccination does not provide complete protection against COVID-19. As a result, the rate of being vaccinated is λ 1 , and λ 2 is the rate of vaccinated persons who have recovered and developed immunity. We suppose that each infectious sub-population (E and I) infected the healthy population at varying densities, where β i is the infection density per capita with i = 1 , 2 . In reality, our major assumption about vaccination’s imperfection is that some vaccinated persons may only receive partial protection and may become sick if they are exposed to multiple infections. The imperfection can be due to the mutation of the virus. In fact, when people are vaccinated, they tend to relax their guard and take fewer protection measures again the virus. θ is the proportion of infected individuals. Some people show mild symptoms with the per capita rate γ 1 and can stay at home with treatment; whilst others develop hard symptoms and must be monitored in the hospital with the per capita rate γ 2 and another critical condition that requires penetrating breathing with per capita rate γ 3 . The parameter δ i , with i = 1 , 2 , 3 , represents the recovery rates of quarantined, hospitalized and critical infected persons (respectively). Average quarantine and hospitalization times are 1 / σ 1 and 1 / σ 2 , respectively. Finally, μ represents the mortality rate due to disease. The flow chart of the model is given in the Figure 1. When the nine equations in (1) are combined together, the total population size N remains constant.

2.1. Positivity and Boundedness

This part is dedicated to establishing the positivity and boundedness of the model (1) solutions under non-negative conditions.
Theorem 1.
If S ( 0 ) 0 , V ( 0 ) 0 , E ( 0 ) 0 , I ( 0 ) 0 Q ( 0 ) 0 , H ( 0 ) 0 , C ( 0 ) 0 , R ( 0 ) 0 and D ( 0 ) 0 , then the solution S ( t ) , V ( t ) , E ( t ) , I ( t ) Q ( t ) , H ( t ) , C ( t ) , R ( t ) , D ( t ) of system (1) is non-negative and the solutions exist in Ω for all t 0 .
The proof of positivity follows the standard argument, as can be seen, for example, in [24].
The solution of the model (1) exists in the positively invariant region:
Ω = { ( S ( t ) , V ( t ) , E ( t ) , I ( t ) , Q ( t ) , H ( t ) , C ( t ) , R ( t ) , D ( t ) ) R + 9 : S ( t ) + V ( t ) + E ( t ) + I ( t ) + Q ( t ) + C ( t ) + R ( t ) + D ( t ) = N N ( 0 ) } .
Since all the parameters and sub populations of the system are non-negative and N(t) = constant t 0 .
Straightforwardly, we obtain S ( t ) N ( t ) N ( 0 ) for all t 0 . The other variables yielded the same result. As a result, the overall population is finitely upper bounded. For system (1), the region Ω is positively invariant.

2.2. Existence and Uniqueness of Solutions

Theorem 2.
The system (1) that fulfills a given initial condition ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) Q ( 0 ) , H ( 0 ) , C ( 0 ) , R ( 0 ) , D ( 0 ) ) has a unique solution.
Proof. 
The system (1) may be expressed as follows:
X ˙ = A X + B ( X )
where X ( t ) = [ S ( t ) , V ( t ) , E ( t ) , I ( t ) , Q ( t ) , H ( t ) , C ( t ) , R ( t ) , D ( t ) ]
A = λ 1 0 0 0 0 0 0 0 0 λ 1 λ 2 0 0 0 0 0 0 0 0 0 θ 0 0 0 0 0 0 0 0 θ γ 1 γ 2 γ 3 0 0 0 0 0 0 0 0 γ 1 σ 1 δ 1 0 0 0 0 0 0 0 γ 2 σ 1 σ 2 δ 2 0 0 0 0 0 0 γ 3 0 σ 2 μ δ 3 0 0 0 λ 2 0 0 δ 1 δ 2 δ 3 0 0 0 0 0 0 0 0 μ 0 0 ,
B ( X ) = β 1 S E N β 2 S I N β 3 V E N β 4 V I N β 1 S E N + β 2 S I N + β 3 V E N + β 4 V I N 0 0 0 0 0 0 .
Equation (2) is a non-linear system that can be written as
Σ ( X ) = A X + B ( X ) .
We have:
| B ( X 1 ) B ( X 2 ) | 2 β 1 N | S 2 E 2 S 1 E 1 | + 2 β 2 N | S 2 I 2 S 1 I 1 | + 2 β 3 N | V 2 E 2 V 1 E 1 | + 2 β 4 N | V 2 I 2 V 1 I 1 | 2 β 1 N | S 1 ( E 2 E 1 ) + E 2 ( S 2 S 1 ) | + 2 β 2 N | S 1 ( I 2 I 1 ) + I 2 ( S 2 S 1 ) | + 2 β 3 N | V 1 ( E 2 E 1 ) + E 2 ( V 2 V 1 ) | + 2 β 4 N | V 1 ( I 2 I 1 ) + I 2 ( V 2 V 1 ) | 2 β 1 ( | E 2 E 1 | + | S 2 S 1 | ) + 2 β 2 ( | I 2 I 1 | + | S 2 S 1 | ) + 2 β 3 ( | E 2 E 1 | + | V 2 V 1 | ) + 2 β 4 ( | I 2 I 1 | + | V 2 V 1 | ) 2 ( β 1 + β 3 ) | E 2 E 1 | + 2 ( β 1 + β 2 ) | S 2 S 1 | + 2 ( β 2 + β 4 ) | I 2 I 1 | + 2 ( β 3 + β 4 ) | V 2 V 1 | 4 M ( | S 2 S 1 | + | V 2 V 1 | + | E 2 E 1 | + | I 2 I 1 | ) , where M = max { β 1 , β 2 , β 3 , β 4 } 4 M X 2 X 1
then, we obtain | Σ ( X 1 ) Σ ( X 2 ) | M ¯ X 1 X 2 , where M ¯ = max { M , A } < . As a result, the function Σ ( t ) is uniformly Lipschitz continuous. The restriction on S ( t ) 0 , V ( t ) 0 , E ( t ) 0 , I ( t ) 0 , Q ( t ) 0 , H ( t ) 0 , C ( t ) 0 , R ( t ) 0 , and D ( t ) 0 . Thus, a solution to the system (2) exists [25].  □

3. The Basic Reproduction Number

The basic reproduction number R 0 is the average number of persons in a susceptible population that one person infected with COVID-19 is expected to infect, and it is calculated using the next generation matrix approach [26]. The disease compartments are thus E and I. Therefore, the all-time disease-free equilibrium point E 0 = ( N 0 , 0 , 0 , 0 , 0 , 0 , 0 ) .
F = β 1 S * E * N + β 2 S * I * N + β 3 V * E * N + β 4 V * I * N 0 , V = V V + = θ E ( γ 1 + γ 2 + γ 3 ) I * θ E *
The Jacobian matrices of F and V computed at E 0 are provided by F and V, respectively, such that:
F = β 1 β 2 0 0 , V = θ 0 θ ( γ 1 + γ 2 + γ 3 )
The inverse of V is given by
V 1 = 1 θ 0 1 ( γ 1 + γ 2 + γ 3 ) 1 ( γ 1 + γ 2 + γ 3 ) and F V 1 = β 1 θ + β 2 γ 1 + γ 2 + γ 3 β 2 γ 1 + γ 2 + γ 3 0 0
Therefore, the domainant eigenvalue of F V 1 :
R 0 = ρ ( F V 1 ) = β 1 θ + β 2 ( γ 1 + γ 2 + γ 3 )

4. Local Stability Analysis at Disease-Free Equilibrium (DFE) E 0

The DFE’s local stability is investigated as follows.
The Jacobian matrix of the system (1) at E 0 = ( N , 0 , 0 , 0 , 0 , 0 , 0 ) is:
J * ( E 0 ) = λ 1 0 β 1 β 2 0 0 0 λ 1 λ 2 0 0 0 0 0 0 0 β 1 θ β 2 0 0 0 0 0 θ ( γ 1 + γ 2 + γ 3 ) 0 0 0 0 0 0 γ 1 ( σ 1 + δ 1 ) 0 0 0 0 0 γ 2 σ 1 ( σ 2 + δ 2 ) 0 0 0 0 γ 3 0 σ 2 ( μ + δ 3 ) .
The eigenvalues of the Jacobian matrix J * ( E 0 ) are the roots of the following characteristic equation:
( λ 1 λ ) ( λ 2 λ ) ( σ 1 δ 1 λ ) ( σ 2 δ 2 λ ) ( μ δ 3 λ ) ( λ 2 + a 1 λ + a 0 ) = 0
where:
a 0 = ( θ δ 1 ) ( γ 1 + γ 2 + γ 3 ) θ β 2 a 1 = ( θ δ 1 ) + ( γ 1 + γ 2 + γ 3 )
The roots of ( λ 2 + a 1 λ + a 0 ) = 0 are given by
λ 1 = β 1 θ ( γ 1 + γ 2 + γ 3 ) ( β 1 θ + γ 1 + γ 2 + γ 3 ) 2 + 4 θ β 2 2 = α 1 γ ( α 1 + γ ) 2 + 4 θ β 2 2 λ 2 = β 1 θ ( γ 1 + γ 2 + γ 3 ) + ( β 1 θ + γ 1 + γ 2 + γ 3 ) 2 + 4 θ β 2 2 = α 1 γ + ( α 1 + γ ) 2 + 4 θ β 2 2
with α 1 = β 1 θ and γ = ( γ 1 + γ 2 + γ 3 ) .
It is straightforward that if α 1 γ < 0 R * = β 1 θ + γ < 1 and λ 1 < 0 .
In the case for λ 2 , we write the equation:
1 R 0 = α 1 γ θ β 2 γ θ
From the previous Equation (8), if R 0 < 1 , then:
α 1 γ + θ β 2 < 0 4 α 1 γ + 4 θ β 2 < 0 2 α 1 γ + 4 θ β 2 < 2 α 1 γ α 1 2 + 2 α 1 γ + 4 θ β 2 + γ 2 < α 1 2 2 α 1 γ + γ 2 ( α 1 + γ ) 2 + 4 θ β 2 < ( γ α 1 ) 2 ( α 1 + γ ) 2 + 4 θ β 2 < γ α 1 α 1 γ + ( α 1 + γ ) 2 + 4 θ β 2 < 0 λ 2 < 0
Using the same steps as in Equation (8), if R 0 > 1 , then:
α 1 γ + θ β 2 > 0 4 α 1 γ + 4 θ β 2 > 0 2 α 1 γ + 4 θ β 2 > 2 α 1 γ α 1 2 + 2 α 1 γ + 4 θ β 2 + γ 2 > α 1 2 2 α 1 γ + γ 2 ( α 1 + γ ) 2 + 4 θ β 2 > ( γ α 1 ) 2 ( α 1 + γ ) 2 + 4 θ β 2 > γ α 1 α 1 γ + ( α 1 + γ ) 2 + 4 θ β 2 > 0 λ 2 > 0
Notice that the fact that R 0 < 1 implies that β 1 θ < 1 and since R * < β 1 θ . Then, we conclude that it is enough to have R 0 < 1 to ensure that the roots of ( λ 2 + a 1 λ + a 0 ) = 0 are negative. On the other hand, if R * > 1 , then R 0 > 1 . As a result, we have just proven the following theorem:
Theorem 3.
If R 0 < 1 , the disease-free equilibrium E 0 of the system (1) is locally asymptotically stable, but unstable if R * > 1 , where: R * = β 1 θ + γ .

5. Global Stability Analysis at Disease-Free Equilibrium

The global stability of the disease-free equilibrium point E 0 was found in this part by creating the Lyapunov function as follows:
Theorem 4.
The disease-free equilibrium E 0 of the model (1) is globally asymptotically stable whenever R 0 T o t a l 1 , where R 0 T o t a l = R 0 + R 0 v .
Proof. 
Consider the Lyapunov function L in the trivial equilibrium point E 0 , which has non-negative coefficients g 1 and g 2 :
L = g 1 E + g 2 I
Differentiating Equation (9) with respect to time t, and substituting both d E d t and d I d t from Equation (1) yields the result:
L ˙ = g 1 E ˙ + g 2 I ˙ = g 1 β 1 S E N + g 1 β 2 S I N + g 1 β 3 V E N + g 1 β 4 V I N g 1 θ E + g 2 θ E g 2 ( γ 1 + γ 2 + γ 3 ) I
By simplifying Equation (10) by collecting similar terms of E and I, then by solving for coefficient g 1 and g 2 , this yields:
L ˙ g 1 β 1 E + g 1 β 3 E + g 2 θ E g 1 θ E + g 1 β 2 I + g 1 β 4 I g 2 ( γ 1 + γ 2 + γ 3 ) I g 1 θ β 1 + β 3 θ + g 2 g 1 1 E + g 1 β 2 + β 4 g 2 ( γ 1 + γ 2 + γ 3 ) I g 1 β 2 + β 4 g 1 ( 1 β 1 + β 3 θ ) ( γ 1 + γ 2 + γ 3 ) I g 1 2 β 2 + β 4 g 1 ( γ 1 + γ 2 + γ 3 ) + β 1 + β 3 θ ( γ 1 + γ 2 + γ 3 ) I g 1 2 ( γ 1 + γ 2 + γ 3 ) β 2 g 1 ( γ 1 + γ 2 + γ 3 ) + β 1 θ 1 + β 4 g 1 ( γ 1 + γ 2 + γ 3 ) + β 3 θ I ( γ 1 + γ 2 + γ 3 ) R 0 + R 0 v 1 I
where g 1 = 1 , g 2 g 1 = 1 β 1 + β 3 θ .
R 0 T o t a l = β 2 ( γ 1 + γ 2 + γ 3 ) + β 1 θ + β 4 ( γ 1 + γ 2 + γ 3 ) + β 3 θ = R 0 + R 0 v .
Since g 2 g 1 > 0 , then β 1 + β 3 θ < 1 . Therefore, L ˙ is negative if R 0 T o t a l < 1 . Furthermore, L ˙ = 0 if and only if I = 0 . It can thus be investigated whether singleton E 0 is the highest compact invariant set for the model (1). Thus, by LaSalle’s invariance principle [27], the DFE is globally asymptotically stable in a region Ω around E 0 .  □
Remark 1.
The above result shows that R 0 T o t a l and R 0 can be reduced to less than a unit so that the disease disappears. Obviously, R 0 < R 0 T o t a l means that if R 0 < 1 , then the complete eradication of disease is guaranteed.

6. The Optimal Imperfect Vaccination

When imperfect vaccination is administered to a population, there is a need to find the optimal approach to use it in order to reduce the burden of the disease in the population. The goal of this section is to implement the best control strategy possible in the situation of an imperfect vaccination. Three types of control are used for this purpose. First, the control u 1 represents the awareness of taking the vaccine via media, as well as creating knowledge of the positive effects of vaccination to gain herd immunity in the population. The second control u 2 is the movement restrictions for susceptible and vaccinated individuals by adhering to a preventative protocol, avoiding the exposure of the vaccinated people to the coronavirus via non-pharmaceutical measures. The third one u 3 seeks to improve the efficacy of the vaccine.
Therefore, the model with control strategies is given by
d S d t = S N ( 1 u 2 ) ( β 1 E + β 2 I ) ( λ 1 + u 1 ) S d V d t = ( λ 1 + u 1 ) S V N ( 1 u 2 ) ( β 3 E + β 4 I ) ( λ 2 + u 3 ) V d E d t = S N ( 1 u 2 ) ( β 1 E + β 2 I ) + V N ( 1 u 2 ) ( β 3 E + β 4 I ) θ E d I d t = θ E ( γ 1 + γ 2 + γ 3 ) I d Q d t = γ 1 I ( σ 1 + δ 1 ) Q d H d t = γ 2 I + σ 1 Q ( σ 2 + δ 2 ) H d C d t = γ 3 I + σ 2 H ( μ + δ 3 ) C d R d t = δ 1 Q + δ 2 H + δ 3 C + ( λ 2 + u 3 ) V d D d t = μ C
With:
( u 1 ( t ) , u 2 ( t ) , u 3 ( t ) ) U a d T
and U a d T is a set of admissible controls defined by
U a d T = u | ( u 1 ( t ) , u 2 ( t ) , u 3 ( t ) ) are measurable , 0 u 1 ( t ) 1 λ 1 , 0 u 2 ( t ) 1 , 0 u 3 ( t ) 1 λ 2 , t [ 0 , T ]
The objective function to minimize is:
J ( u 1 ( t ) , u 2 ( t ) , u 3 ( t ) ) = 0 T [ A 1 V ( t ) + A 2 I ( t ) A 3 R ( t ) + 1 2 τ 1 u 1 2 ( t ) + 1 2 τ 2 u 2 2 ( t ) + 1 2 τ 3 u 3 2 ( t ) ] d t
The positive weight constants A 1 , A 2 and A 3 , respectively, keep the sizes of V ( t ) , I ( t ) , and R ( t ) in balance. Positive weight parameters: τ 1 , τ 2 , and τ 3 are related with the controls u 1 ( t ) , u 2 ( t ) , and u 3 ( t ) in the objective functional.
To solve the problem, we first compute the Lagrangian and Hamiltonian. Equation (16) in order to identify an optimal solution. The optimal problem’s Lagrangian is:
L = A 1 V ( t ) + A 2 I ( t ) A 3 R ( t ) + 1 2 τ 1 u 1 2 ( t ) + 1 2 τ 2 u 2 2 ( t ) + 1 2 τ 3 u 3 2 ( t ) .
For the control problem, we may define the Hamiltonian H ¯ as follows:
H ¯ = L + ζ 1 ( t ) d S d t + ζ 2 ( t ) d V d t + ζ 3 ( t ) d E d t + ζ 4 ( t ) d I d t + ζ 5 ( t ) d Q d t + ζ 6 ( t ) d H d t + ζ 7 ( t ) d C d t + ζ 8 ( t ) d R d t + ζ 9 ( t ) d D d t
where ζ 1 , , ζ 9 are the adjoint functions to be found.
We have the existence result:
Theorem 5.
The optimal control problem, defined by Equations (13)–(16), has a solution ( u 1 * , u 2 * , u 3 * ) that satisfies
J ( u 1 * , u 2 * , u 3 * ) = min ( u 1 , u 2 , u 3 ) U a d T J ( u 1 , u 2 , u 3 )
Proof. 
We use the result [28] to show that an optimal control exists. The control and the state variables are both non-negative. This minimization problem satisfies the convexity requirement of the objective functional.
The control space previously defined as (15) is both convex and closed by definition. In order for the optimal control to exist, it is necessary for the optimal system to be compact. The boundedness of the optimal system determines the compactness needed. Additionally, an integrand throughout the functional (16) is convex on the control ( u 1 ( t ) , u 2 ( t ) , u 3 ( t ) ) . It can be concluded that the constant ρ > 1 exists, as do positive integers w 1 , w 2 and w 3 such that J ( u 1 , u 2 , u 3 ) w 2 + w 1 ( ( u 1 , u 2 , u 3 ) 2 ) ρ 2 . This leads us to conclude that optimal control exists.  □

Characterization of the Optimal Control

We then investigate the necessary optimal control conditions. For this purpose, the maximum principle of Pontryagin to Hamiltonian [29] can be applied.
Theorem 6.
Let S * ( t ) , V * ( t ) , E * ( t ) , I * ( t ) , Q * ( t ) , H * ( t ) , C * ( t ) , R * ( t ) and D * ( t ) represent optimal state solutions with optimal control variables ( u 1 * ( t ) , u 2 * ( t ) , u 3 * ( t ) ) for the optimal control problem (16). ζ 1 , , ζ 9 are thus adjoint variables that satisfy:
ζ ˙ 1 ( t ) = ( 1 u 2 ) ( β 1 E ( t ) N + β 2 I ( t ) N ) ( ζ 1 ( t ) ζ 3 ( t ) ) + ( λ 1 + u 1 ) ( ζ 1 ( t ) ζ 2 ( t ) ) ζ ˙ 2 ( t ) = A 1 + ( 1 u 2 ) ( β 3 E ( t ) N + β 4 I ( t ) N ) ( ζ 2 ( t ) ζ 3 ( t ) ) + ( λ 2 + u 3 ) ( ζ 2 ( t ) ζ 8 ( t ) ) ζ ˙ 3 ( t ) = ( 1 u 2 ) β 1 S ( t ) N ( ζ 1 ( t ) ζ 3 ( t ) ) + ( 1 u 2 ) β 3 V ( t ) N ( ζ 2 ( t ) ζ 3 ( t ) ) + θ ( ζ 3 ( t ) ζ 4 ( t ) ) ζ ˙ 4 ( t ) = A 2 + ( 1 u 2 ) β 2 S ( t ) N ( ζ 1 ( t ) ζ 3 ( t ) ) + ( 1 u 2 ) β 4 V ( t ) N ( ζ 2 ( t ) ζ 3 ( t ) ) + γ 1 ( ζ 4 ( t ) ζ 5 ( t ) ) + γ 2 ( ζ 4 ( t ) ζ 6 ( t ) ) + γ 3 ( ζ 4 ( t ) ζ 7 ( t ) ) ζ ˙ 5 ( t ) = σ 1 ( ζ 5 ( t ) ζ 6 ( t ) ) + δ 1 ( ζ 5 ( t ) ζ 8 ( t ) ) ζ ˙ 6 ( t ) = σ 1 ( ζ 6 ( t ) ζ 7 ( t ) ) + δ 2 ( ζ 6 ( t ) ζ 8 ( t ) ) ζ ˙ 7 ( t ) = μ ( ζ 7 ( t ) ζ 9 ( t ) ) + δ 3 ( ζ 7 ( t ) ζ 8 ( t ) ) ζ ˙ 8 ( t ) = A 3 ζ ˙ 9 ( t ) = 0
Conditions of transversality.
ζ 1 ( T ) = ζ 3 ( T ) = ζ 5 ( T ) = ζ 6 ( T ) = ζ 7 ( T ) = ζ 9 ( T ) = 0 , ζ 2 ( T ) = A 1 , ζ 4 ( T ) = A 2 , ζ 8 ( T ) = A 3 ,
Moreover, the optimal control ( u 1 * , u 2 * , u 3 * ) is provided by
u 1 * ( t ) = max { min { S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) , 1 λ 1 } , 0 } . u 2 * ( t ) = max { min { S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) , 1 } , 0 } . u 3 * ( t ) = max { min { V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) , 1 λ 2 } , 0 }
Proof. 
By using Hamiltonian (18), Pontryagin’s maximum principle and setting S ( t ) = S * ( t ) , V ( t ) = V * ( t ) , E ( t ) = E * ( t ) , I ( t ) = I * ( t ) , Q ( t ) = Q * ( t ) , H ( t ) = H * ( t ) , C ( t ) = C * ( t ) , R ( t ) = R * ( t ) and D ( t ) = D * ( t ) to obtain the following:
d ζ 1 d t = ( 1 u 2 ) ( β 1 E ( t ) N + β 2 I ( t ) N ) ( ζ 1 ( t ) ζ 3 ( t ) ) + ( λ 1 + u 1 ) ( ζ 1 ( t ) ζ 2 ( t ) ) d ζ 2 d t = A 1 + ( 1 u 2 ) ( β 3 E ( t ) N + β 4 I ( t ) N ) ( ζ 2 ( t ) ζ 3 ( t ) ) + ( λ 2 + u 3 ) ( ζ 2 ( t ) ζ 8 ( t ) ) d ζ 3 d t = ( 1 u 2 ) β 1 S ( t ) N ( ζ 1 ( t ) ζ 3 ( t ) ) + ( 1 u 2 ) β 3 V ( t ) N ( ζ 2 ( t ) ζ 3 ( t ) ) + θ ( ζ 3 ( t ) ζ 4 ( t ) ) d ζ 4 d t = A 2 + ( 1 u 2 ) β 2 S ( t ) N ( ζ 1 ( t ) ζ 3 ( t ) ) + ( 1 u 2 ) β 4 V ( t ) N ( ζ 2 ( t ) ζ 3 ( t ) ) + γ 1 ( ζ 4 ( t ) ζ 5 ( t ) ) + γ 2 ( ζ 4 ( t ) ζ 6 ( t ) ) + γ 3 ( ζ 4 ( t ) ζ 7 ( t ) )
d ζ 5 d t = σ 1 ( ζ 5 ( t ) ζ 6 ( t ) ) + δ 1 ( ζ 5 ( t ) ζ 8 ( t ) ) d ζ 6 d t = σ 1 ( ζ 6 ( t ) ζ 7 ( t ) ) + δ 2 ( ζ 6 ( t ) ζ 8 ( t ) ) d ζ 7 d t = μ ( ζ 7 ( t ) ζ 9 ( t ) ) + δ 3 ( ζ 7 ( t ) ζ 8 ( t ) ) d ζ 8 d t = A 3 d ζ 9 d t = 0
Using optimality conditions, we conclude that:
d H ¯ ( t ) d u 1 ( t ) = τ 1 u 1 * ( t ) + S ( t ) ( ζ 2 ( t ) ζ 1 ( t ) ) d H ¯ ( t ) d u 2 ( t ) = τ 2 u 2 * ( t ) + S ( t ) N β 1 E ( t ) + β 2 I ( t ) ( ζ 1 ( t ) ζ 3 ( t ) ) + V ( t ) N β 3 E ( t ) + β 4 I ( t ) ( ζ 2 ( t ) ζ 3 ( t ) ) d H ¯ ( t ) d u 3 ( t ) = τ 3 u 3 * ( t ) + V ( t ) ( ζ 8 ( t ) ζ 2 ( t ) ) .
Hence:
d H ¯ ( t ) d u 1 ( t ) = 0 u 1 * ( t ) = S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) , d H ¯ ( t ) d u 2 ( t ) = 0 u 2 * ( t ) = S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) d H ¯ ( t ) d u 3 ( t ) = 0 u 3 * ( t ) = V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) .
By applying the control space property, we obtain that:
u 1 * = 0 if S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) 0 u 1 * = S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) if 0 < S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) < 1 u 1 * = 1 if S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) 1 λ 1 .
which means that:
u 2 * = 0 if S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) 0 u 2 * = ω * if 0 < S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) < 1 u 2 * = 1 if S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) 1 .
Since:
ω * = S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) .
We have:
u 3 * = 0 if V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) 0 u 3 * = V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) if 0 < V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) < 1 u 3 * = 1 if V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) 1 λ 2 .
Thus, optimal control is defined as
u 1 * = max { min { S ( t ) τ 1 ( ζ 1 ( t ) ζ 2 ( t ) ) , 1 λ 1 } , 0 } , u 2 * = max { min { S ( t ) τ 2 N β 1 E ( t ) + β 2 I ( t ) ( ζ 3 ( t ) ζ 1 ( t ) ) + V ( t ) τ 2 N β 3 E ( t ) + β 4 I ( t ) ( ζ 3 ( t ) ζ 2 ( t ) ) , 1 } , 0 } . u 3 * = max { min { V ( t ) τ 3 ( ζ 2 ( t ) ζ 8 ( t ) ) , 1 λ 2 } , 0 }
 □

7. Numerical Simulation

The goal of this section is to show how the control strategies can be used to improve outcomes in vaccination campaigns in Morocco. After fitting the model to the data, the estimated parameters are taken to perform sensitivity analysis and determine the optimal control.

7.1. Parameter Estimation

We used data from a vaccination campaign in Morocco between 1 February 2021 and 25 March 2021 to validate our findings. We consider the data from the COVID-19 data [30] and the initial conditions within the values of available data on 1 February 2021 are ( S 0 , V 0 , E 0 , I 0 , Q 0 , H 0 , C 0 , R 0 , D 0 ) = (36,202,000, 200,081, 25,543, 13,099, 9824, 1572, 131, 450,052, 8287). These initial values were estimated from the data, apart from E 0 and Q 0 , which were assumed. To obtain the best fitting curve for actual data, we applied the least-squares fitting technique [31]. The parameter values of the model are estimated based on this fitting and are given as follows: λ 1 = 3.23× 10 3 , λ 2 = 1.5× 10 4 , θ = 7.4385× 10 2 , σ 1 = 3.9601× 10 2 , σ 2 = 5.1954× 10 2 , δ 1 = 2.04× 10 2 , δ 2 = 1.01× 10 2 , δ 3 = 0.00203 , β 1 = 9.813× 10 3 , β 2 = 2.4× 10 3 , β 3 = 0.21 , β 4 = 0.21 , γ 1 = 9.81× 10 2 , γ 2 = 1.2× 10 2 , γ 3 = 1.02× 10 4 and μ = 1.501× 10 3 . The value of the basic reproduction number in this case is R 0 = 0.1536999 . The fit of the number of individuals infected with COVID-19 in Morocco is described in Figure 2.

7.2. Sensitivity Analysis

In this part, we aimed to study the sensitivity of the model parameters with respect to R 0 T o t a l . The goal is to determine the impact of these parameters on the endemicity of the disease. The sensitivity analysis for this outbreak threshold demonstrates the importance of each parameter in the spread of COVID-19, allowing us to determine which parameter to make more sensitive on R 0 T o t a l with respect to a specific parameter, ρ , via the sensitivity index define by
ζ ρ R 0 T o t a l = R 0 T o t a l ρ ρ R 0 T o t a l .
Using the previous definition:
ζ β 1 R 0 T o t a l = 1 θ β 1 R 0 T o t a l , ζ β 2 R 0 T o t a l = 1 ( γ 1 + γ 2 + γ 3 ) β 2 R 0 T o t a l ζ β 3 R 0 T o t a l = 1 θ β 3 R 0 T o t a l , ζ β 4 R 0 T o t a l = 1 ( γ 1 + γ 2 + γ 3 ) β 4 R 0 T o t a l ζ γ 1 R 0 T o t a l = ( β 2 + β 4 ) ( γ 1 + γ 2 + γ 3 ) 2 γ 1 R 0 T o t a l , ζ γ 2 R 0 T o t a l = ( β 2 + β 4 ) ( γ 1 + γ 2 + γ 3 ) 2 γ 2 R 0 T o t a l ζ γ 3 R 0 T o t a l = ( β 2 + β 4 ) ( γ 1 + γ 2 + γ 3 ) 2 γ 3 R 0 T o t a l , ζ θ R 0 T o t a l = ( β 1 + β 3 ) θ 2 θ R 0 T o t a l .
Each parameter’s sensitivity index, corresponding to the basic reproductive numbers R 0 T o t a l , was computed and displayed in Table 1, and the graphical bar-graph findings were generated in Figure 3. The sensitivity indices indicate the importance of each parameter in disease transmission and prevalence. These are some examples: if ζ β 1 R 0 T o t a l = 0.931 , it means that if β 1 went up (or decreased) by 93.1%, R 0 T o t a l is also likely to increase or decrease by 93.1%. Similarly, for ζ θ R 0 T o t a l = 0.8202 , the decrease in the parameter θ by 82.02% will fall (or increase ) R 0 T o t a l by a similar proportion.
Our goal is to simulate the elasticity of R 0 T o t a l with respect to model parameters in four scenarios that represent the different scenarios of the dynamics of the model as follows:
  • Scenario 1: There is no disease persistence ( R 0 < 1 , R 0 v < 1 and R 0 T o t a l < 1 ).
  • Scenario 2: Persistence of diseases with low threshold values ( R 0 < 1 , R 0 v < 1 and R 0 T o t a l > 1 ).
  • Scenario 3: The disease persists, with a high endemicity among vaccinated people and a low endemicity among non-vaccinated people ( R 0 < 1 , R 0 v > 1 and R 0 T o t a l > 1 ).
  • Scenario 4: The disease persists, with low endemicity among the vaccinated and high endemicity among the unvaccinated ( R 0 > 1 , R 0 v < 1 and R 0 T o t a l > 1 ).
All of the sensitivity scenarios described above demonstrate that the basic reproductive number R 0 T o t a l is more sensitive to some parameters than others—particularly in the cases of β 1 , β 3 , γ 1 and θ . Scenario 1 has no persistent disease, but Scenario 3 has persistent disease with significant endemicity among the vaccinated, as can be seen from the sensitivity indices of γ 1 (rate of infected getting isolated) and θ (incubation period). The same observation applies to the scenarios of the persistence of disease with low endemicity in Scenario 2 and the persistence of disease with high endemicity among the vaccinated in Scenario 4. The level of sensitivity of these parameters is higher when endemicity is low among the non-vaccinated population. Our simulations revealed that the value of the sensitivity index changes depending on the Scenarios (1, 2, 3, 4), with β 3 being the highest value for Scenarios 1 and 2. However, in β 1 , we can see more domination for the index in Scenario 4.

7.3. Simulation of Optimal Control

The time series of variables in the model without and with optimal control are shown in the figures below. The goal is to compare the effect of control on the different variables of the model.
The three variables are shown in Figure 4: susceptible, infected and recovered without and with optimal control effect. These simulations show that optimal control increases the number of susceptible and recovered people while decreasing the number of infected individuals. The effect of control on the susceptible and recovered populations is obviously more significant since the susceptible and recovered populations rose four-fold in 50 days while the infected population declined by one fold during the same period. This means that the control is very effective for all of these three compartments.
The simulation of the time series of the variables representing the vaccinated, exposed and deceased individuals is shown in Figure 5. Clearly, the vaccinated population benefits more from the optimal control than the exposed population since the control aims to improve vaccination effectiveness. The control strategy, on the other hand, has no obvious effect on the number of deaths.
Similarly, as seen in Figure 6, the control strategy reduces the number of isolated populations (at home or in hospital). However, there is more benefit to controlling the population with mild symptoms compared to people in the hospital with severe or critical symptoms. When applying the control to all three categories of quarantined, hospitalized, and critical cases of COVID-19, it is clear that the control is obvious.
The three types of optimal controls are given in Figure 7, Figure 8 and Figure 9. Those figures show the intensity of each measure needed to be implemented in the case of imperfect vaccination. The awareness campaign should stay as long as 7 weeks with maximum intensity. At the same time, non-pharmaceutical measures must take place, including mobility restrictions or lockdown for at least 24 days. With regard to the third control measure, the improvement of the vaccine population should increase to reach possible efficacy during the first week of vaccination and stay above 30 % within the first 50 days of vaccination. The outcome of this control approach shows that these three measures must be simultaneously implemented to deal with imperfect vaccination.

8. Conclusions

COVID-19 is still taking a toll on people’s lives all over the world. Countries are rushing to implement the vaccination to gain herd immunity, to contain the spread of the disease, and to bring the fatality rate of the disease to the lowest possible level. However, the administered vaccines have different levels of efficacy among the same population, which means the idea of relying on vaccination alone to control the pandemic would not provide total protection for the population against further waves of COVID-19. In this work, we aimed to present a mathematical model of the imperfect vaccination of COVID-19 and to study the dynamics of this model. One further element of this model is that susceptible and vaccinated people have different infection rates and the fatality rate of the disease (rate of death due to the infection) is related to the percentage of the isolated population that is in critical condition. We derive a threshold R 0 which is the basic reproduction number of disease transmission among the population. We showed that this threshold does not give us sharp epidemiological properties of the model. In fact, we prove, via the Lyapunov method, that the disease is globally asymptotically stable if R 0 T o t a l 1 with R 0 T o t a l is the sum of R 0 and R 0 v , which is the threshold of transmission of the disease among the vaccinated population. This finding shows that the increase in the efficacy of the vaccination should lead to protecting the population from infection (low β 3 and β 4 ) which will help control the pandemic.
To make our analysis more realistic, we estimated the parameters of our model using data from Morocco between 1 February 2021 and 25 March 2021. Within the range of the estimated parameters, we performed a sensitivity analysis of R 0 T o t a l with respect to the parameters of the model to find the elasticity index with respect to each parameter. Depending on the disease status, our simulation (Figure 3) produced four outcomes. Scenario 1: There is no disease persistence. Scenario 2: Persistence of diseases with low threshold values. Scenario 3: The disease persists with a high endemicity among vaccinated people and a low endemicity among non-vaccinated people. Scenario 4: The disease persists, with low endemicity among the vaccinated and high endemicity among the unvaccinated.
To further investigate the possible additional measures that help with vaccination. We introduce an optimal control problem with the goal of increasing the awareness of vaccination, limiting the probability of infection by adhering to a preventative protocol and increasing the efficacy of the vaccine. The public health authorities can easily implement these measures. In fact, as the virus mutates, many governments are pushing their populations to get vaccinated and asking people to reduce their contact and wear masks. Moreover, there is a constant effort to increase the efficacy of the vaccine by producing new ones or by boosters.
Our solution of optimal control showed that, to reduce the impact of imperfect vaccination, we needed a longer awareness campaign to engage the population in vaccination. On the other hand, the restriction on population mobility should not be long, since our simulation showed a drop of u 2 from 1 (full restriction) after 24 days. To ensure the full protection of the health population, vaccination efficacy must increase by 30 % in the first 50 days.
In conclusion, our work showed that facing the imperfection of the vaccination of COVID-19, we mainly have to focus on two measures. The first one is to increase awareness of the importance of vaccination, which will increase the number of people vaccinated. The second is to work on developing vaccines with high efficacy that give more protection to the population instead of making the symptoms of viruses less severe. With these two measures, our study showed that population mobility restrictions have the lowest impact on controlling the virus spread.

Author Contributions

Conceptualization, L.B. and A.T.; methodology, M.L. and A.T.; software, L.B. and M.L.; validation, A.T., M.L. and O.Z.; formal analysis, L.B. and A.T.; investigation, L.B.; resources, M.L. and M.R.; data curation, L.B. and M.R.; writing—original draft preparation, L.B. and A.T.; writing—review and editing, L.B. and A.T.; visualization, L.B. and M.L.; supervision, M.L., A.T. and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the reviewers for their comments and inputs that help improve the quality of our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Corey, L.; Mascola, J.R.; Fauci, A.S.; Collins, F.S. A strategic approach to COVID-19 vaccine R&D. Science 2020, 368, 948–950. [Google Scholar] [PubMed]
  2. Arino, J.; Cooke, K.; Van Den Driessche, P.; Velasco-Hernández, J. An epidemiology model that includes a leaky vaccine with a general waning function. Discret. Contin. Dyn. Syst.-B 2004, 4, 479. [Google Scholar] [CrossRef]
  3. Abu-Raddad, L.J.; Boily, M.C.; Self, S.; Longini, I.M., Jr. Analytic insights into the population level impact of imperfect prophylactic HIV vaccines. JAIDS J. Acquir. Immune Defic. Syndr. 2007, 45, 454–467. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, X.; Takeuchi, Y.; Iwami, S. SVIR epidemic models with vaccination strategies. J. Theor. Biol. 2008, 253, 1–11. [Google Scholar] [CrossRef] [PubMed]
  5. Gulbudak, H.; Martcheva, M. A structured avian influenza model with imperfect vaccination and vaccine-induced asymptomatic infection. Bull. Math. Biol. 2014, 76, 2389–2425. [Google Scholar] [CrossRef] [PubMed]
  6. Al-Darabsah, I. Threshold dynamics of a time-delayed epidemic model for continuous imperfect-vaccine with a generalized nonmonotone incidence rate. Nonlinear Dyn. 2020, 101, 1281–1300. [Google Scholar] [CrossRef] [PubMed]
  7. Djilali, S.; Bentout, S. Global dynamics of SVIR epidemic model with distributed delay and imperfect vaccine. Results Phys. 2021, 25, 104245. [Google Scholar] [CrossRef]
  8. Al-Darabsah, I. A time-delayed SVEIR model for imperfect vaccine with a generalized nonmonotone incidence and application to measles. Appl. Math. Model. 2021, 91, 74–92. [Google Scholar] [CrossRef] [PubMed]
  9. Iboi, E.A.; Ngonghala, C.N.; Gumel, A.B. Will an imperfect vaccine curtail the COVID-19 pandemic in the US? Infect. Dis. Model. 2020, 5, 510–524. [Google Scholar] [PubMed]
  10. Mancuso, M.; Eikenberry, S.E.; Gumel, A.B. Will vaccine-derived protective immunity curtail COVID-19 variants in the US? Infect. Dis. Model. 2021, 6, 1110–1134. [Google Scholar] [CrossRef]
  11. Webb, G. A COVID-19 Epidemic Model Predicting the Effectiveness of Vaccination in the US. Infect. Dis. Rep. 2021, 13, 62. [Google Scholar] [CrossRef] [PubMed]
  12. Libotte, G.B.; Lobato, F.S.; Platt, G.M.; Neto, A.J.S. Determination of an optimal control strategy for vaccine administration in COVID-19 pandemic treatment. Comput. Methods Programs Biomed. 2020, 196, 105664. [Google Scholar] [CrossRef] [PubMed]
  13. Olivares, A.; Staffetti, E. Optimal control-based vaccination and testing strategies for COVID-19. Comput. Methods Programs Biomed. 2021, 211, 106411. [Google Scholar] [CrossRef] [PubMed]
  14. Shim, E. Optimal allocation of the limited COVID-19 vaccine supply in South Korea. J. Clin. Med. 2021, 10, 591. [Google Scholar] [CrossRef] [PubMed]
  15. Choi, W.; Shim, E. Vaccine Effects on Susceptibility and Symptomatology Can Change the Optimal Allocation of COVID-19 Vaccines: South Korea as an Example. J. Clin. Med. 2021, 10, 2813. [Google Scholar] [CrossRef] [PubMed]
  16. Jordan, E.; Shin, D.E.; Leekha, S.; Azarm, S. Optimization in the Context of COVID-19 Prediction and Control: A Literature Review. IEEE Access 2021, 9, 130072–130093. [Google Scholar] [CrossRef]
  17. Viana, J.; van Dorp, C.H.; Nunes, A.; Gomes, M.C.; van Boven, M.; Kretzschmar, M.E.; Veldhoen, M.; Rozhnova, G. Controlling the pandemic during the SARS-CoV-2 vaccination rollout. Nat. Commun. 2021, 12, 3674. [Google Scholar] [CrossRef] [PubMed]
  18. Matrajt, L.; Eaton, J.; Leung, T.; Brown, E.R. Vaccine optimization for COVID-19: Who to vaccinate first? Sci. Adv. 2021, 7, eabf1374. [Google Scholar] [CrossRef] [PubMed]
  19. Choi, Y.; Kim, J.S.; Kim, J.E.; Choi, H.; Lee, C.H. Vaccination Prioritization Strategies for COVID-19 in Korea: A Mathematical Modeling Approach. Int. J. Environ. Res. Public Health 2021, 18, 4240. [Google Scholar] [CrossRef] [PubMed]
  20. Makhoul, M.; Ayoub, H.H.; Chemaitelly, H.; Seedat, S.; Mumtaz, G.R.; Al-Omari, S.; Abu-Raddad, L.J. Epidemiological impact of SARS-CoV-2 vaccination: Mathematical modeling analyses. Vaccines 2020, 8, 668. [Google Scholar] [CrossRef] [PubMed]
  21. Wagner, C.E.; Saad-Roy, C.M.; Morris, S.E.; Baker, R.E.; Mina, M.J.; Farrar, J.; Holmes, E.C.; Pybus, O.G.; Graham, A.L.; Emanuel, E.J.; et al. Vaccine nationalism and the dynamics and control of SARS-CoV-2. medRxiv 2021, 373, eabj7364. [Google Scholar] [CrossRef] [PubMed]
  22. Piraveenan, M.; Sawleshwarkar, S.; Walsh, M.; Zablotska, I.; Bhattacharyya, S.; Farooqui, H.H.; Bhatnagar, T.; Karan, A.; Murhekar, M.; Zodpey, S.; et al. Optimal governance and implementation of vaccination programmes to contain the COVID-19 pandemic. R. Soc. Open Sci. 2021, 8, 210429. [Google Scholar] [CrossRef] [PubMed]
  23. Mohammed-Awel, J.; Numfor, E.; Zhao, R.; Lenhart, S. A new mathematical model studying imperfect vaccination: Optimal control analysis. J. Math. Anal. Appl. 2021, 500, 125132. [Google Scholar] [CrossRef]
  24. Thieme, H.R. Princeton series in theoretical and computational biology. In Mathematics in Population Biology; Princeton University Press: Princeton, NJ, USA, 2003. [Google Scholar]
  25. Birkhoff, G.; Rota, G. Ordinary Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 1989. [Google Scholar]
  26. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef]
  27. LaSalle, J. Stability theory for ordinary differential equations. J. Differ. Equ. 1968, 4, 57–65. [Google Scholar] [CrossRef] [Green Version]
  28. Fleming, W.; Rishel, R. Stochastic Differential Equations and Markov Diffusion Processes. In Deterministic and Stochastic Optimal Control; Springer: Berlin/Heidelberg, Germany, 1975; pp. 106–150. [Google Scholar]
  29. Lee, D.; Milroy, I.P.; Tyler, K. Application of Pontryagin’s maximum principle to the semi-automatic control of rail vehicles. In Proceedings of the Second Conference on Control Engineering, Newcastle, UK, 25–27 August 1982; pp. 233–236. [Google Scholar]
  30. Roser, M.; Ritchie, H.; Ortiz-Ospina, E.; Hasell, J. Coronavirus Pandemic (COVID-19). Our World in Data. 2022. Available online: https://ourworldindata.org/coronavirus (accessed on 15 January 2022).
  31. Weisstein, E.W. Least Squares Fitting. 2002. Available online: https://mathworld.wolfram.com/ (accessed on 15 January 2022).
Figure 1. Flow chart of the model.
Figure 1. Flow chart of the model.
Axioms 11 00124 g001
Figure 2. Fitting the infected population with real data from 1 February 2021 to 25 March 2021.
Figure 2. Fitting the infected population with real data from 1 February 2021 to 25 March 2021.
Axioms 11 00124 g002
Figure 3. Scenarios of R 0 T o t a l sensitivity respecting the parameters.
Figure 3. Scenarios of R 0 T o t a l sensitivity respecting the parameters.
Axioms 11 00124 g003
Figure 4. Time series of the states S, I and R without and with control.
Figure 4. Time series of the states S, I and R without and with control.
Axioms 11 00124 g004
Figure 5. The time series of states V, E and D without and with control.
Figure 5. The time series of states V, E and D without and with control.
Axioms 11 00124 g005
Figure 6. Time series of the states Q, H and C with and without control.
Figure 6. Time series of the states Q, H and C with and without control.
Axioms 11 00124 g006
Figure 7. Evolutionary dynamics of the control u 1 ( t ) .
Figure 7. Evolutionary dynamics of the control u 1 ( t ) .
Axioms 11 00124 g007
Figure 8. Evolutionary dynamics of the control u 2 ( t ) .
Figure 8. Evolutionary dynamics of the control u 2 ( t ) .
Axioms 11 00124 g008
Figure 9. Evolutionary dynamics of the control u 3 ( t ) .
Figure 9. Evolutionary dynamics of the control u 3 ( t ) .
Axioms 11 00124 g009
Table 1. Sensitivity index for each parameter that has a direct correlation to R 0 T o t a l .
Table 1. Sensitivity index for each parameter that has a direct correlation to R 0 T o t a l .
ScenariosScenario 1Scenario 2Scenario 3Secnario 4
ParametersValuesS. IndexValuesS. IndexValuesS. IndexValuesS. Index
β 1 6.71 × 10 3 2.55 × 10 1 9.81 × 10 3 1.28 × 10 1 9.81 × 10 3 2.70 × 10 2 9.81 × 10 2 6.56 × 10 1
β 2 2.70 × 10 3 6.93 × 10 2 2.40 × 10 3 2.11 × 10 2 2.40 × 10 3 4.46 × 10 3 2.40 × 10 2 1.08 × 10 1
β 3 1.10 × 10 2 4.18 × 10 1 5.10 × 10 2 6.66 × 10 1 2.10 × 10 1 5.78 × 10 1 2.10 × 10 2 1.40 × 10 1
β 4 1.00 × 10 2 2.57 × 10 1 2.10 × 10 2 1.85 × 10 1 2.10 × 10 1 3.90 × 10 1 2.10 × 10 2 9.48 × 10 2
γ 1 9.81 × 10 2 −2.90 × 10 1 9.81 × 10 2 −1.84 × 10 1 9.81 × 10 2 −3.51 × 10 1 9.81 × 10 2 −1.81 × 10 1
γ 2 1.20 × 10 2 3.08 × 10 1 1.20 × 10 2 1.06 × 10 1 1.20 × 10 2 2.23 × 10 2 1.20 × 10 2 5.42 × 10 2
γ 3 2.00 × 10 6 2.62 × 10 3 1.02 × 10 4 8.99 × 10 4 1.02 × 10 4 1.90 × 10 4 1.02 × 10 4 4.61 × 10 4
θ 7.44 × 10 2 −2.20 × 10 1 7.44 × 10 2 −1.39 × 10 1 7.44 × 10 2 −2.66 × 10 1 7.44 × 10 2 −1.37 × 10 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Benahmadi, L.; Lhous, M.; Tridane, A.; Zakary, O.; Rachik, M. Modeling the Impact of the Imperfect Vaccination of the COVID-19 with Optimal Containment Strategy. Axioms 2022, 11, 124. https://doi.org/10.3390/axioms11030124

AMA Style

Benahmadi L, Lhous M, Tridane A, Zakary O, Rachik M. Modeling the Impact of the Imperfect Vaccination of the COVID-19 with Optimal Containment Strategy. Axioms. 2022; 11(3):124. https://doi.org/10.3390/axioms11030124

Chicago/Turabian Style

Benahmadi, Lahbib, Mustapha Lhous, Abdessamad Tridane, Omar Zakary, and Mostafa Rachik. 2022. "Modeling the Impact of the Imperfect Vaccination of the COVID-19 with Optimal Containment Strategy" Axioms 11, no. 3: 124. https://doi.org/10.3390/axioms11030124

APA Style

Benahmadi, L., Lhous, M., Tridane, A., Zakary, O., & Rachik, M. (2022). Modeling the Impact of the Imperfect Vaccination of the COVID-19 with Optimal Containment Strategy. Axioms, 11(3), 124. https://doi.org/10.3390/axioms11030124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop