Optimal surveillance mitigation of COVID'19 disease outbreak: Fractional order optimal control of compartment model.

In present time, the whole world is in the phase of war against the deadly pandemic COVID'19 and working on different interventions in this regard. Variety of strategies are taken into account from ground level to the state to reduce the transmission rate. For this purpose, the epidemiologists are also augmenting their contribution in structuring such models that could depict a scheme to diminish the basic reproduction number. These tactics also include the awareness campaigns initiated by the stakeholders through digital, print media and etc. Analyzing the cost and profit effectiveness of these tactics, we design an optimal control dynamical model to study the proficiency of each strategy in reducing the virulence of COVID'19. The aim is to illustrate the memory effect on the dynamics of COVID'19 with and without prevention measures through fractional calculus. Therefore, the structure of the model is in line with generalized proportional fractional derivative to assess the effects at each chronological change. Awareness about using medical mask, social distancing, frequent use of sanitizer or cleaning hand and supportive care during treatment are the strategies followed worldwide in this fight. Taking these into consideration, the optimal objective function proposed for the surveillance mitigation of COVID'19, is contemplated as the cost function. The effect analysis is supported through graphs and tabulated values. In addition, sensitivity inspection of basic reproduction number is also carried out with respect to different values of fractional index and cost function. Ultimately, social distancing and supportive care of infected are found to be significant in decreasing the basic reproduction number more rapidly.

Introduction

A deadly coronavirus that basically initiated from Wuhan city of China, all of a sudden incarcerated the people all around the world. This strain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected more than 210 countries and territories. It has brought devastating consequences on public health as well as on social and economic activities. Governments around the world prompted surveillance on mitigating the global spread of COVID'19. Among many of these dramatic measures, majority are substantiating to be effective in reducing the virus transmission. Imposing curfew and locking down the cities in addition public awareness campaigns such as, stay-at-home, encouraging social distancing, cleanliness that include frequent washing hand, using sanitizers through digital and print media are the key measures in restraining this virus. On enforcing these policies and engaging communities in these campaigns, undoubtedly enormous social and economic cost is expected. But until an effectual vaccine or treatment becomes available, these strategies may play important roles [1], [2], [3], [4], [5], [6].

Variety of research has been conducted at an extraordinary pace to analyze the COVID'19 in different perspectives [7], [8], [9]. Epidemiological dynamical systems to control the breakout of this pandemic through basic reproduction number has been obtained by various researchers [10], [11], [12]. Clinical studies to determine therapeutic solutions through the findings of the biological features of this virus [13]. Perceptions on impact of government's preventing strategies on other environmental, social and economic activities [14]. Decision making models to consider an effective managing prevention strategy of COVID'19 transition [15]. Machine learning models to predict the high risk and efficiently triage the patients with high accuracy [16]. Mathematical models fit out to be substantial contrivances in investigating the dynamical controls of the infectious diseases[17]. Research articles, based on optimal control models can be found in the literature to a great extent in this regard [18], [19], [20], [21], [22], [23], [24]. In the recent times of battle against COVID'19, numerous authors have added their valuable contributions in this connection. Grigorieva et al. formulate two SEIR-type model to investigate the cost-effective quarantine strategies and analyzed the optimal solutions numerically [25]. Analysis of interventions of COVID'19 through transmission model and observing the most effective non-pharmaceutical strategies to lessen the disease nuisance in Pakistan is found in the literature [26]. In particular, plenty of endeavors have been carried out in different context of cost-effective strategy, to control the transmission of this deadly pandemic [27], [28].

In this attempt, we design mathematical model that covers two major areas, epidemiology together with dynamical optimal control. Firstly, the compartmental model is taken into account with the control variables and stability analysis are carried out. Secondly, optimizing cost functions is subjected to the compartmental model to assess the cost-effectiveness of the prevention strategies. As aforementioned, there exist significant mathematical efforts in this connection, but the novelties that invigorate the proposed assessment can be classified as:

This study is not only susceptible, exposed, quarantine, infected and recovery compartments, but also the isolation and precautions. Thus, the model is named as SEQIMRP i.e. susceptible-expose-quarantined-infected-isolated-recovered-protected.

The non-pharmaceutical control variables, awareness campaigns about using mask, encouraging social distancing, signifying frequent use of sanitizer and washing hands, supportive care during treatment.

Regulatory of basic reproduction number through these campaigns.

Incorporating fractional order derivative for dynamical scrutiny of the model with.

This significant contribution will undoubtedly add great perspicacity of COVID'19 interventions. The proposed SEQIMRP model with proportional fractional [29] signifies the broader application of the fractional definition. Its expansion elegantly converts the fractional order derivative operator into integer order that the fractional order index reallocates linearly in the equations. By virtue of this, the dynamics of COVID’19, for instance the basic reproduction number and equilibrium points can be interpreted with memory effects. Subsequently, historical values of these parameters or the compartmental functions will enable to devise defensive precautionary steps, revealed from the past experiences. In addition, the effect of memory on the optimality of awareness strategies is also illustrated through the proportional fractional derivative. The designed system provides a novel contribution in epidemiological study of epidemic and pandemic diseases. It will instruct the healthcare researchers a new mode of generating results and might be capable to investigating prior information about the risk factors or transmission rate for preparatory measures. The remaining paper contains sections of formulating the dynamical system, stability analysis of equilibrium points and optimality assessment. Furthermore, numerical discussions are also carried out to evidently establish an effective conclusion.

Model formulation for COVID'19 optimal control

Susceptible-expose-quarantined-infected-isolated-recovered-protected (SEQIMRP)

Mathematical models based on disease dynamics are quite helpful in studying the functional behavior of any virus, which then helps to overcome or lessen its contaminating breakout. The destructive coronavirus converted into a pandemic within a few months and affected billions of peoples around a globe. Early laboratory research and scientific experiments to construct a drug or vaccine could not triumph. Many epidemiological models also expressed significant contributions in this connection to determine the basic reproduction number and predict the dispersion, recovery and mortality rates [11], [12], [30]. Here, to analyze the dynamical behavior and impact of COVID'19 pandemic, a system of differential equations is designed with respect to compartmental classes and prevention measures on the basis of following assumptions.

Regardless of different risk rate of COVID'19 for different age-group and pre-existing disease carriers, the model assumes a homogeneous mixing of individuals in the population.

Prevention strategies: Usage of medical mask , social distancing , frequently cleaning hands and supportive care during treatments are taken into account as control variables of the optimal system.

The individuals in any compartment, following the operational prevention strategies, are assumed as will not get infected and are defined by means of protected compartment.

Susceptible is outlined in the form of logistic growth that encompasses maximum sustainability to survive in the available resources in an environment.

Exposed are quarantined that might recover and use prevention measures later to insulate themselves from virus.

Treatment of infected COVID'19 patients is the isolation process, which is explained in the isolation compartment. These compartments with the supportive care from staff might recover and move to recovery compartment.

Treated individuals after recovery do not participate in transmitting the disease as they use the operational prevention strategies.

The assessments of basic reproduction number and stability analysis are carried out in fractional calculus environment.

Hence, the fractional order epidemiological model, susceptible-expose-quarantined-infected-isolated-recovered-protected (SEQIMRP) is mathematically signified as:

with initial conditions,

where, for. Table 1 further elaborates the dimensions of all the variables and parameters of system (1). Moreover, pictorial demonstration of the compartmental system, representing the flow of the diseases transmission is also given in Fig. 1 . Assume is the total population density of individuals that can be structured as:

Table 1

Variables and parameters of the SEQIMRP.

Compart-mental functions	Descriptions	Units	Initial values(population = in millions & time = days)	Source
Nt	Total population	Population / day	113	Estimated
St	Susceptible	Population / day	111	Estimated
Et	Exposed	Population / day	0	Estimated
Qt	Quarantined	Population / day	0	Estimated
It	Infected	Population / day	2	Estimated
Mt	Infected isolated	Population / day	0	Estimated
Rt	Recovered	Population / day	0	Estimated
Pt	Protected	Population / day	0	Estimated
α	Order of fractional derivative	Dimensionless	0<α⩽1	Fitted
Parameters	Descriptions	Units	Value	Source
ξ	R is taking part in social distancing	Individuals/(individuals × day)	14.771	Fitted
β	Contact rate of susceptible with infected	Individuals/(individuals × day)	14.781	Fitted
γ	Rate of exposed individuals quarantined	Individuals/(individuals × day)	1.887 × 10-7	Fitted
η	Rate of treated infected and quarantined	Individuals/(area × day)	0.13266	Fitted
σ	Rate of susceptible exposed to quarantine	Individuals/(individuals × day)	0.0714	Fitted
r_S	Intrinsic Growth rate of susceptible individuals	Individuals/(individuals × day)	30	Fitted
k_S	Carrying capacity of susceptible individuals	Individuals/(individuals × day)	100,000	Fitted
σ	Rate of susceptible exposed to infection	Individuals/(individuals × day)	0.1259	Fitted
ρ	Death due to COVID'19 disease	Individuals/(individuals × day)	1.782 × 10-5	[35]
ψ_Q	Recovery rate of the quarantine individuals	Individuals/(individuals × day)	0.11624	Fitted
ψ_M	Recovery rate of the isolated infected individuals	Individuals/(individuals × day)	0.33029	Fitted
d_S	Susceptible death rate	Individuals/(individuals × day)	0.15	Fitted
d_E	Exposed death rate	Individuals/(individuals × day)	0.84	Fitted
d_Q	Quarantined death rate	Individuals/(individuals × day)	0.84	Fitted
d_I	Infected death rate	Individuals/(individuals × day)	0.9	Fitted
d_R	Recovered death rate	Individuals/(individuals × day)	0.11	Fitted
d_P	Precautionary death rate	Individuals/(individuals × day)	0.84	Fitted
mm	Rate of individuals using medical mask	Individuals/(individuals × day)	0–1	Fitted
sd	Rate of individuals taking part in social distancing	Individuals/(individuals × day)	0–1	Fitted
ch	Rate of individuals frequently cleaning hand	Individuals/(individuals × day)	0–1	Fitted
sc	Rate of individuals who follow step of supportive care during treatment.	Individuals/(individuals × day)	0–1	Fitted

Fig. 1

Pictorial illustration of SEQIMRP model.

Moreover, articulates proportional fractional derivative of order [29], which can be expanded as for any continuous function ,

where, for , with and . Additionally, for , with and. Let, and, so Eq. (4) becomes

Assume that all control functions (prevention steps) are constant within time, therefore, by applying expansion (5) on system (1), we get the system as:

with the same initial conditions (2). System (6) evidently depicts the lucidity of the proportional fractional derivative, which greatly reduces the manipulation complexities of system (1).

(Boundedness)

Let is the set of all feasible solutions of the system (6), then there exists uniformly bounded subset of such that:

By applying proportional fractional derivative and its expansion, as defined in the Eqs. (4)-(5), on Eq. (3), we get the expression of the form:

On simplifying by using system (6) and suppose be total proportion of deaths in all compartments i.e.

In addition, since

where , so the above inequality reduces to

On integrating

Therefore as , we obtained the final statement of boundedness as

(Existence and Uniqueness)

Assume the matrix of right hand side of system (6) be the real-valued function , such that and are continuous and

Then, satisfying the initial conditions (2), there exists a unique, non-negative and bounded solution of the system (6).

Boundedness of system (6) can be followed from Theorem 1, now assume, the system (6) can be expressed as:

where,

and

Eq. (16) can be further expanded into:

such that

and . Then, Eq. (17) can be rewritten as,

Let , so the final statement is achieved as for ,

where . Next, we prove the non-negativity of the solutions by using the positivity of initial conditions (2) i.e., for . Considering first equation of system (6), it can be deduced to:

On manipulating, we get

Since for , therefore Eq. (18) reduces to,

Thus, proved the non-negativity of . Analogously, all the remaining equations of system (6) can be proved to have non-negative solutions with the assumption of positive initial conditions.

Optimal control problem

Furthermore, the dynamical model (6) of COVID'19 would be incomplete if the assumption of optimal control of infection and intervention cost is not incorporated. Therefore, we formulate optimal control problem by means of the cost function type of quadratic function as:

where, for are replace by respectively. Moreover, here , for , are the weights of human population cost, whereas , for , are the weights of undertaken intervention cost for COVID'19. At this juncture, intervention cost comes from government campaigns of using mask, social distancing and frequently washing hand. In addition, the hospitalization cost for drugs, ventilators and trained medical staffs for supportive care of the COVID'19 infected individuals also become higher with the increase in number of patients. Therefore, if greater cost is implemented of campaigns of enforcing the people on usage of mask, social distancing and frequently washing hand will reduce the COVID'19 transmission, which on the other hand it reduces the supportive care cost. Thus, we assume , for . Analogously, the objective of the present scenario is to control the spread out of COVID'19, which ultimately leads to minimize the infected individuals, therefore we consider and remaining equal to zero.

Basic reproduction number

In this sequel, we utilize the next generation method, to structure the for the governing model (6). For this purpose, a sub-model of the SEQIMRP is considered that includes the four infected classes i.e. exposed, quarantine, infected and isolated individuals. Therefore, the equation:

will have as a vector of the,, and , which is outlined as,

with, expressed as,

On the other hand,, can be further split down as,

Taking Jacobian matrix of Eq. (20) at disease free equilibrium point,, we get,

where,

, ,

From Eq. (21), we can extract and manipulate,

The spectral radius is the required basic reproduction number, so after some simplification we get

Consequently, the generated contains the fractional derivative index as well, which advantageously enables to inspect . The health care researchers will be capable to investigate the trajectory of basic reproduction number for the COVID'19 at small change.

Dynamical anatomization

In this section, on the strength of proportional fractional derivative, dynamical analysis of equilibrium points and optimality conditions are discussed in fractional environment as follows:

Systematic stability analysis

Trivial Equilibrium Point

The trivial equilibrium solution, , of system (6), is asymptotically unstable, for .

It can be easily proved by eigenvalues of at , for all ,

,,,,

, , .

Since, it is clear that, for. Thus, is unstable.

Disease Free Equilibrium Point

The disease-free equilibrium of the system (6)

For , is locally asymptotically stable if and unstable when , for .

On manipulating Jacobian at , the negative eigenvalues i.e. for , are attained as:

,,,,

with the equation,

where

where

On applying Routh-Hurwitz criteria [31], [32], [33], [34] i.e. if , and then polynomial (23) is greater than zero and thus all the real part of the eigenvalues must be negative. It can be evidently seen that for , now the thing which left to prove is . Hence, is locally asymptotically stable if and if , implies that is Eq. (23) must have a nonnegative real part, thus becomes unstable.

Endemic Equilibrium Point

The endemic equilibrium is locally asymptotically stable if and only if, , for .

The Jacobian at , generates the negative real eigenvalues,

with the polynomial equation,

where,

where,

and

The factor , if the magnitude of , which implies that becomes positive if and only if . Thus with reference to Lemma 5.1 of [20], the positive constant of the polynomial implies is locally asymptotically stable if .

Characterization of optimal control

It is evidently clear from Theorem 1 that there exist a unique solution of system (6). Now to optimize the solution, we define the Lagrangian by

In addition, describing the Hamiltonian as the inner product of the right hand side of the state system (6) and the adjoint variables , we get

where is to be determined. Now, utilizing the Pontryagin’s maximum principle for the Hamiltonian , following theorem is obtained to determine the adjoint variables.

(Existence of adjoint variable)

For the controlling functions , , and together with the solution of the corresponding system (6), there exists adjoint variables that satisfy,

Furthermore, the optimal control pairs are descripted as:

By using Pontryagin’s maximum principle in state, the adjoint equations with transversality conditions is stated as:

with transversality where . By using optimality condition, we deduce the optimal control pairs as:

Further, taking into account the property of the control space, we achieve,

, ,

where,

Ultimately, the control pair and state variables are found by using the following composed systems:

and

Numerical simulation and deliberation

In this segment, numerical investigations of the aforementioned system are carried out by considering some numerical values of the parameters, as shown in Table 1. The graphical predisposition analysis of with respect to the strategies are also added in the discussion. Moreover, the simulations of all compartmental class, with prevention and without prevention campaign cases are plotted and tabulated by using Mathematica 11.0.

Sensitivity analysis of parameters with optimality

The sensitivity analysis of by means of control variables are described in Table 2 and Fig. 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7 for the parameters mentioned in Table 1 and at different values of . These control variables define the strategic campaigns utilized to prevent the deadly transmission of the COVID'19. It can be clearly seen from the Fig. 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7 that at each value of , the influential strength of each campaign together minimizes the significance of . The generation of colorized output in these figures, ranging from light to dark, indicates the gradual decrease in from largest to lowest value. The obtained value of without any awareness campaign is greater than 1, which gradually reduces to less than 1 on increasing awareness campaigns that can be seen from the Table 2 and Fig. 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7. Furthermore, the lines of on the Fig. 2, which are attained by fixing and and varying and represent decrease in value starting from 1.4 to 0.4. In the same way, Fig. 3, Fig. 5 plotted for and , respectively which demonstrate the same pattern of decrease in . On the other hand, Fig. 4 exhibit the decrease in starting from 1.8 to 0.4, for . Similarly, same sketches are found in Fig. 6, Fig. 7, which are produced by fixing and , accordingly. Besides, Table 2 explains the sensitivity of with some different values of intervention strategies, which elucidates that for , the value of decreases more rapidly than the other combinations, at each value of . In addition, the last column of Table 2 also elaborates the minimum cost function against each mitigation strategy for weights and . Evidently, the optimal values of for the cost efforts of surveillance mitigations, , which greatly reduce are 9068.92, 9090.57 and 9098.23 for , respectively and . According to these values, increasing the awareness about social distancing and supportive care of the infected individuals will significantly affect the transmission of COVID'19 with optimal cost efforts, comparative to other combinations of mitigations.

Table 2

Sensitivity inspection of and optimal surveillance based on prevention strategies for weights , and at different values of .

α	Intervention Strategies	R_0	J
0.8	mm=0,sd=0,ch=0,sc=0	2.08323	6171.69
	mm=0.3,sd=0.5,ch=0.65,sc=0.9	0.236717	9702.47
	mm=0.3,sd=0.7,ch=0.5,sc=0.9	0.227723	9068.92
	mm=0.6,sd=0.5,ch=0.3,sc=0.9	0.246215	9017.57
0.95	mm=0,sd=0,ch=0,sc=0	2.83546	9622.94
	mm=0.3,sd=0.5,ch=0.65,sc=0.9	0.266884	9724.52
	mm=0.3,sd=0.7,ch=0.5,sc=0.9	0.252657	9090.57
	mm=0.6,sd=0.5,ch=0.3,sc=0.9	0.278141	9040.02
1	mm=0,sd=0,ch=0,sc=0	3.17529	10815.9
	mm=0.3,sd=0.5,ch=0.65,sc=0.9	0.278141	9732.32
	mm=0.3,sd=0.7,ch=0.5,sc=0.9	0.266884	9098.23
	mm=0.6,sd=0.5,ch=0.3,sc=0.9	0.290079	9047.98

Fig. 2

Sensitivity inspection of with respect to and for , at .

Fig. 3

Sensitivity inspection of with respect to and for , at .

Fig. 4

Sensitivity inspection of with respect to and for , at .

Fig. 5

Sensitivity inspection of with respect to and for , at .

Fig. 6

Sensitivity inspection of with respect to and for , at .

Fig. 7

Sensitivity inspection of with respect to and for , at .

Equilibrium states and optimality

Moreover, solving SEQIMRP system different plots are attained that define the stability of and . In the current scenario, evaluations of these equilibrium points are produced on the basis of the prevention campaigns. Commencing from Table 3 , the values are generated for and , at and . Manifestly, it can be seen when no prevention measures are taken increases gradually and endemic state of the pandemic becomes stable. Additionally, Fig. 8, Fig. 9, Fig. 10, Fig. 11, Fig. 12, Fig. 13, Fig. 14 also plotted for and , at and represent the stability of the deadly endemic state of COIVD'19 for the current rate of transmission, recovery and mortality. Since, no prevention measures are taken at initial spread stage of COVID'19, therefore the curve of protected population yields a constant straight line on zero. This further elaborates the circumstances where everyone is at high risk of being infected that the pandemic situation becomes worst.

Table 3

Basic reproduction number and endemic equilibrium points , for parameters describe in Table 1, for and , at different values of and .

α	R_0	St	Et	Qt	It	Mt	Rt	Pt
0.4	1.07658	90564.3	130,860	0.0166747	0.14077	0.0108021	0.00775503	0
0.5	1.2482	78379.6	341,367	0.0437936	0.394839	0.0339465	0.0267258	0
0.6	1.46208	67141.9	503,966	0.0650954	0.629689	0.0615493	0.0546976	0
0.7	1.73309	56834.8	623,165	0.0810462	0.84565	0.0957691	0.100128	0
0.8	2.08323	47,442	703,299	0.0921023	1.04305	0.140393	0.184117	0
0.9	2.54613	38947.9	748,529	0.0987099	1.22222	0.202724	0.373485	0
1.	3.17529	31335.7	762,847	0.101305	1.38348	0.298912	1.00458	0

Fig. 8

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 9

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 10

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 11

Dynamics of of SEQIMRP, for parameters described in Table 1 for and , at and .

Fig. 12

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 13

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 14

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Contrarily Table 4 depicts the values, which are generated for and , at and . Evidently from Table 4, when prevention measures are taken into account to some extent, we attain the disease-free state of the dynamics at each value of . In addition, it also shows the value of to be less than one, which is proved in Theorem 4. Fig. 15, Fig. 16, Fig. 17, Fig. 18, Fig. 19, Fig. 20, Fig. 21 graphically demonstrate the stability of for and , at and . Contrary to endemic, in disease free case the infected cells become zero whereas susceptible and protected individuals remain at a population level other than zero. Running awareness campaigns about using mask, social distancing, hand wash and also invigorating supportive care of the patients will decrease the basic reproduction number and eventually the deadly spread of COVID'19.

Table 4

Basic reproduction number and disease free equilibrium points , for parameters describe in Table 1, for and , at different values of and .

α	R_0	St	Et	Qt	It	Mt	Rt	Pt
0.4	0.297661	94566.7	0	0	0	0	0	97904.3
0.5	0.324551	94,900	0	0	0	0	0	111,349
0.6	0.354992	95233.3	0	0	0	0	0	128,931
0.7	0.389628	95566.7	0	0	0	0	0	152,907
0.8	0.429253	95,900	0	0	0	0	0	187,538
0.9	0.474857	96233.3	0	0	0	0	0	241,958
1.	0.527691	96566.7	0	0	0	0	0	339,915

Fig. 15

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 16

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 17

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 18

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 19

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 20

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Fig. 21

Dynamics of of SEQIMRP, for parameters described in Table 1 and and , at and .

Conclusion

The declaration of PHEIC by the WHO about the COVID'19 outbreak, agitate the scientific community and the healthcare professionals of the countries. After the failure of several experiments on the inoculations, the only operational plan of action to decelerate the spread of COVID'19 is to adopt non-pharmaceutical restrictions. For this purpose, different unprecedented measures are taken into account such as lockdown, closure of institutions and initiating different awareness campaigns. Here, we discussed the cost and public effectiveness of the awareness campaigns taken into consideration by the stakeholders. These maneuvers include the strict imposition of using medical mask in public places, social distancing of 6 feet, frequent use of hand wash and sanitizers, training medical staffs and officers for extraordinary supportive care of COVID'19 patients in hospitals. The optimal control function was designed with the epidemic dynamical system SEQIMRP to mutually study its dynamical stability and the feasibility of the prevention tactics. The system was formulated with the proportional fractional derivative, in order to analyze the basic reproduction number at each chronological change. Ultimately, through the aforementioned analytical and numerical illustrations, the following propitious facts can be extracted:

The strategies of using medical mask, social distancing, frequently sanitizing hands and supportive care of COVID'19 for speedy recovery are significant attempts to win this battle against this pandemic.

The awareness and necessitating of these lines of attacks may change the state of pandemic into a stable disease-free environment.

These can greatly lesser the basic reproduction number from to .

The optimal surveillance mitigation with respect to cost effectiveness, social distancing and supportive care may reduce the diffusion of COVID'19 more hastily.

Illustrations at different fractional derivative index show systematic reading in the susceptible, expose, quarantined, infected, isolated, recovered and protected population.

Without precautions, as the fractional derivative approaches the whole change, the readings represent step by step increase in susceptible, expose, quarantined, infected, isolated and recovered population.

Following precautions, as the fractional derivative approaches the whole change, the number individuals in protection increases gradually, while expose, quarantined, infected, isolated and recovered remain zero.

Competency in prior recognition of the track of COVID'19 transmission risk through the proportional fractional derivative model.

Proficiently trace the basic reproduction number and take preparatory measures before becoming a deadly pandemic.

In the current phase, understanding the epidemiological characteristics is a serious bone of contention question for researchers and health professionals. The successful investigations may significantly help out the stakeholders in making effective standard operational procedures of interventions. The designed model SEQIMRP will categorically aid a great contribution in dynamically scrutinizing and exhibiting the optimal strategy to control the deadly escalation of COVID'19.

Funding

This research received no external funding.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.