key: cord-1009981-tnuvjqra
authors: Majee, Suvankar; Adak, Sayani; Jana, Soovoojeet; Mandal, Manotosh; Kar, T. K.
title: Complex dynamics of a fractional-order SIR system in the context of COVID-19
date: 2022-01-14
journal: J Appl Math Comput
DOI: 10.1007/s12190-021-01681-z
sha: 9ffef328b9d66015c6e87f227e39c163d0d1f04c
doc_id: 1009981
cord_uid: tnuvjqra

This paper proposes and analyses a new fractional-order SIR type epidemic model with a saturated treatment function. The detailed dynamics of the corresponding system, including the equilibrium points and their existence and uniqueness, uniform-boundedness, and stability of the solutions are studied. The threshold parameter, basic reproduction number of the system which determines the disease dynamics is derived, and the condition of occurrence of backward bifurcation is also determined. Some numerical works are conducted to validate our analytical results for the commensurate fractional-order system. Hopf bifurcations for the fractional-order system are studied by taking the order of the fractional differential as a bifurcation parameter.

Disease is an integral part of modern civilization. As long as human lives, illnesses will last in our society. But, the extent of these diseases can be reduced, and even we can eradicate the disease from the population. Among various types of diseases, infectious diseases are the most vulnerable diseases. Any contagious disease, caused by the pathogen of bacteria, viruses, fungi, or some parasites can transmit from an infected human to a susceptible human either through direct contact (viz. influenza, tuberculosis, rubella, measles, HIV-AIDS, etc.) or via some medium (viz. malaria, dengue, chikungunya, etc.). Among all infectious diseases, the recent pandemic COVID-19 caused by the virus SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the first pandemic of the twenty-first century has affected the lifestyle of human population of every country throughout the globe. The disease COVID-19 was first noticed in the Wuhan city of the Republic of China in December of last year [38] . Till today except for very few countries, COVID-19 pandemic has worsened the situation all over the globe [37] ). It should be noted that till May 2020, more than six million world populations have been affected by the disease COVID-19, where almost 50% closed cases, including 12% death (see [37] ).

The rate of progression of infectious diseases including COVID-19 depends on the amount of pathogen living within a host, its rate of growth, and the interaction with the host's immunity. A mathematical model can perfectly capture these relations. Mathematical models help to represent all the scenarios in a compact form. It also captures how disease propagates, its long-term behavior, and ultimately helps to determine possible interventions to restrain the infectious disease. The mathematical model of an infectious disease was first introduced by Kermack and McKendrick [19] . They have used a simple SIR model to formulate the mathematical model. After this theoretical epidemiology has come a long way. Several works on control strategies like vaccination, treatment, isolation were taken into consideration by many researchers [17, 20, 24] . But, in all the works done so far, researchers have considered the integer order system while constructing the mathematical model. Recently some researchers have also developed some model-based works to examine the detailed dynamical behavior of the COVID-19. The researchers including Ahmed et al. [3] , Kucharski et al. [21] , Prem et al. [32] , etc. have also developed some mathematical models and presented rigorous analysis. In an article, Ribeiro et al. [33] have used some stochastic-based regression models to forecast the phenomena in ten most affected states of Brazil. Noting that isolation is a good measure to control the disease, Hellewell et al. [14] have used the operative techniques of COVID-19 diseases using separation. In the recent works, Mandal et al. [25, 27] have constructed and analyzed a mathematical model on the COVID-19 in the pandemic scenario of three states of India. In other recent work, a balance between lockdown and compliance to investigate the COVID-19 scenes has been proposed by Zegarra et al. [2] . Some more theoretical works on COVID-19 can be found in [25, 27, 34] .

The shortcoming of the integer order system is that it does not rely on the previous history of the system. But, when a disease transmits, the susceptible population uses their memory to prevent infection. A dynamical system involving fractional-order derivatives bears the information regarding its present and also its past states [12] .

Hence, fractional-order systems give more in-depth knowledge about the system than integer order systems. Also, as there are many developments in fractional calculus (See the articles [1, 4, 28, 29] ), so, it is logical to apply a fractional-order differential equation to represent a biological phenomenon through mathematical modeling. According to Du et al. [10] , if the parameter α(∈ (0, 1]), the order of the fractional-order system tends to 0, the model has full memory, and when α = 1, the model system is known as memoryless. Hartley et al. [13] used fractional differentiation for the first time for a physical system in demonstrating the behavior of the semi-infinite lossy (RC) line. Torvik et al. [35] have studied that fractional derivatives appear naturally to describe some properties of motions of a Newtonian fluid. At present, there are several works involving the application of fractional-order systems in various fields including epidemic models [15, 26] .

In human society, any epidemic can not be thought without the effect of memory. Whenever an infectious disease propagates in human society, the knowledge or experience of each people about that particular disease should affect their response [36] . People apply several precautions (e.g. vaccination, if it is available) for a disease if they have a past experience of the disease. In this regard, some important control measures can suppress the spreading of the disease. But, mere knowledge about a disease may not protect all the time. That is why people are willing to adopt some new paths to control these diseases. The memory of an infectious disease from the previous situation has more contribution on the recent situation as compared to no memory. It is anticipated that the long-time impact of memory decreases in future time more tardily than exponential decay, but generally, its behavior is a similar way of power-law damping function.

In this work, we have formulated a new SIR type epidemic model with the help of fractional differential equations. Although there are many ways to define a fractionalorder differential equation in the sense of Riemann-Liouville, Grünwald-Letnikov, and the Caputo. In this work, we use the Caputo definition due to its similarity in using the initial conditions with the integer-order differential equation. From this perspective, we develop a mathematical model considering a system of fractional-order differential equations and analyze the thorough dynamics of the proposed model. Moreover, we apply this model to predict the COVID-19 situation of India after the second wave.

The whole manuscript is structured in the following way: An SIR type model with the fractional differential equation is formulated in Sect. 2. The equilibrium points, existence and uniqueness of solution, their dynamical behavior, stability of the system, and the occurrence of Hopf bifurcation are discussed in Sect. 3. Several numerical simulations are performed in Sect. 4, and the application of the model on some realworld data is discussed in Sect. 5, and lastly, in Sect. 6, we present a concise discussion and conclusion.

In this section, a new SIR type model is considered. We divide the whole population into susceptible S(t), infected I (t), and recovered R(t), the three mutually exclusive timedependent classes. First, we consider that at any time t, A be the newly recruitment rate, β be the disease transmission rate, d be the natural death rate, δ be the rate of death due to the disease, be the natural recovery rate. We assume that some infected person will be recovered due to available treatment control. Let v denote the treatment control parameter for the infected population. The rate of recovery is considered as γ v 1+v I with γ representing the effectiveness of the treatment due to saturated type treatment control (see Jana et al. [16] ). Every country, including the most developed countries in the world, has limitations in providing medical resources like proper medicines, availability of beds in hospitals, etc. Particularly, the contagious disease like COVID-19, which has appeared all of a sudden, it is not possible to provide treatment to every patient when the number of patients becomes very large. From this point of view, we have used the saturated type recovered rate due to treatment. A schematic diagram of the above assumptions is described in Fig. 1 . Based on these assumptions, we define our model in the following manner:

with initial conditions

According to Du et al. [10] a memory process usually consists of two stages: the fresh stage and the working stage. The fresh memory is used to just remember things while the working memory helps to do cognitive works. In the fractional model, we capture the working memory. The critical point between the fresh stage and the working stage is usually not the origin. This observation is quite different from the traditional fractional models of one stage. For example, the fractional Maxwell's model is a one-stage model. As the combination of two simple models, it has a more complicated expression than Eq. (2.1). We also find that the order of fractional. Also, in reality, no individual has the same memories, it changes from person to person. Once infected an individual gains knowledge about the disease which helps to obstruct the transmission of infectious disease. We consider this fact into the fractional form of the system (2.1). Here we assume that α 1 , α 2 and α 3 are such properties corresponding to the susceptible, infected and recovered population respectively, where α i ∈ (0, 1], i = 1, 2, 3. Whereas, the recovered will have better working memory than infected one but less working memory than susceptible. So, the fractional order form of the system (2.1) takes the form

Here D α i t is the fractional-order derivative in Caputo sense [30, 31] for i = 1, 2, 3. In the system (2.2) we replace the parameters of the system (2.1) in such a way that the problem of dimension mismatch is solved.

For simplicity, we set

Using these transformations into the system (2.2), we have

is the fractional-order derivative in Caputo sense [30, 31] for i = 1, 2, 3. The system (2.3) is called to be an incommensurate fractional-order system. If α 1 = α 2 = α 3 = α ∈ (0, 1], then the system is said to be a commensurate fractional-order system which is given by

is the fractional-order derivative in Caputo sense [30, 31] .

At this stage, we analyze the dynamical behavior of the fractional-order system (2.4).

Now we discuss the existence of nonnegative equilibrium points of the system (2.4). We observe that the system (2.4) has three equilibrium points. One is disease free equilibrium (DFE) E 0 = A d , 0, 0 and the other two equilibrium points are endemic equilibrium (EE) if they exist E 1 (S * , I * , R * ) and E 2 (S * * , I * * , R * * ) where I * and I * * are the roots of the following quadratic equation

Here the parameter R 0 is recognized as the basic reproduction number. Hence we can derive S * = 1

Thus we can say that the system (2.1) has a positive unique EE (S * ,

We formulate the next theorem which is grounded on the above discussions.

From the above theorems we experience that the system (2.4) has two different equilibrium points even when R 0 < 1. Hence we can conclude that the proposed model experiences a backward bifurcation at R 0 < 1 which shown in Fig. 2 . Thus the condition R 0 < 1 is not sufficient to eradicate the disease from the population. undergoes a backward bifurcation at R 0 = 1. Hence it can be explained that the system possesses two EE point in some interval

. Now, this case occurs only when a 1 < 0 and R 0 < 1. Thus we get,

From the above conditions we have derived at the critical value of v as v * =

. Thus the new threshold value at which the diseases eradicates

. Therefore, depending on the above discussion we can state the next theorem. (ii) two EE points if R 0 < R 0 < 1 which implies that the system passes through a backward bifurcation. (iii) a unique DFE point if R 0 < R 0 < 1 which implies that the disease will die out.

We now state and prove the following lemma which helps us to conclude that our constructed model has a unique solution.

Lemma 3.4 [22] Consider the fractional-order system

If the Lipschitz condition with respect to x is satisfied by the function f (t, x), then the system (3.2) has a unique solution on the interval [t 0 , ∞) × R.

Theorem 3. 5 The system (2.4) has a unique solution in the region

Proof We prove the existence and uniqueness criterion for the fractional-order system (2.4) . For that purpose we consider the region

be arbitrary two points in the region R and define the mapping H :

For any X , Y ∈ R, we have

Therefore, the function H(X) fulfills the criteria of the Lipschitz's condition with respect to X = (S, I , R) ∈ R. Then by using the Lemma 3.4, we can conclude that our system (2.4) has a unique solution X ∈ R with initial conditions X t 0 = (S t 0 , I t 0 , R t 0 ) ∈ R.

In mathematical biology, our main concern is on the nonnegative and bounded solutions. Hence, we try to find the solutions that possess the (i) non-negativity and (ii) bounded properties. For this purpose, we consider the following region R + = {(S, I , R) ∈ R : S, I , R ∈ R + }.

Theorem 3.6 Every solutions of the fractional system (2.4) whose initial values start in R + 3 i.e. S t 0 ≥ 0, I t 0 ≥ 0, R t 0 ≥ 0 are all nonnegative and bounded.

Proof To prove the nonnegativity of the solution, we use Lemma 3.2 [15] . We have discussed the existence of solution of the system (2.1) in Theorem 3.1. As (2.4) is a homogeneous system of equations with initial conditions S(t) ≥ 0, I (t) ≥ 0, R(t) ≥ 0 we can conclude that the solutions of (2.4) is non-negative. Furthermore, the solutions are nondecreasing in t. Now to prove the uniform boundedness of the solutions, we consider a mapping F(t) = S(t) + I (t) + R(t). Then we have,

Using the Lemma 3.2 [15] , we have

Therefore, all the solution of the fractional system (2.4) whose initial conditions start in the region R + is bounded and it lies in the region {(S + I + R) ∈ R + : S + I + R ≤ A d }.

In this portion of the research article, we study the stability conditions of each equilibrium point. Here, we assume that the fractional system (2.4) possesses unique EE point. The local stability investigation of the system (2.4) is done with the help of making linearization around each equilibrium point. The Jacobian matrix of the fractional system (2.4) at the DFE point ( A d , 0, 0) is

The eigenvalues of the above Jacobian matrix are λ 1 = −d,

Again we have established that the system undergoes backward bifurcation iff R 0 < R 0 . Hence, by the Lemma 3.3 [15] , we can conclude that the system is asymptotically stable around the DFE if R 0 < min{1, R 0 }. Thus, we can state the next theorem.

Next, we discuss the local stability criteria of EE point (S * , I * , R * ). J , the Jacobian matrix at (S * , I * , R * ) of the system (2.4) is given by

One eigenvalue of the matrix is λ 1 = −d and the other two eigenvalues are obtained from the matrix

The characteristic equation of the Jacobian matrix at the EE, (S * , I * , R * ) is given by

Let λ 2 and λ 3 be the eigenvalues of the matrix J * . Now, a 1 will be negative if

Then the eigenvalues are

Previously, we have determined one eigenvalue λ 1 = −d whose argument i.e. |arg(−d)| = π > απ 2 , for α ∈ (0, 1). Now, we consider different cases considering different values of b 1 and a 2 .

(i) If a 2 1 ≥ a 2 and a 1 < 0, then both the eigenvalues are negative. So, |arg(λ 2,3 )| = π > απ 2 , α ∈ (0, 1). Hence from Lemma 3.3 and Lemma 3.4 of the article by Jana et al. [15] , the EE is asymptotically stable. (ii) If a 2 1 ≥ a 2 and a 1 ≥ 0, then one of the eigenvalues will be non negative. Then, |arg(λ i )| = 0 < απ 2 for i=2 or, 3, and α ∈ (0, 1). Hence, by Lemma 3.3 and 3.4 [15] , the EE is unstable. (iii) If the conditions a 2 1 < a 2 and a 1 > 0 hold, then both the eigenvalues will be complex conjugates. λ 2,3 = a 1 ± i a 2 − a 2 1 , where i = √ −1. Therefore, the EE point will be asymptotically stable if |arg(λ 2,

Hence, for the order of differentiation α ∈ 0, 2 π tan −1 | 

(v) If a 2 1 < a 2 and a 1 < 0, then |arg(λ 2,3 )| = π 2 > απ 2 . Then EE is locally asymptotically stable.

From the above discussions we may state the next theorem about the stability of the EE point.

(i) The EE is locally asymptotically stable if a 2 1 ≥ a 2 and a 1 < 0. (ii) The EE is locally unstable if a 2 1 ≥ a 2 and a 1 ≥ 0. (iii) The EE is locally asymptotically stable if a 2 1 < a 2 and a 1 > 0. (iv) The EE is locally asymptotically stable if a 2 1 < a 2 and a 1 < 0. (v) The EE is locally asymptotically stable if a 2 1 < a 2 and a 1 = 0.

We now demonstrate the asymptotic global stability of the DFE point. Then the solution of the fractional system (2.4) and the constructed function V together imply that

The Lemma 4.3 [15] suggests that any solution of (2.1) starting in the region R tends to largest invariant set S = {(S, I , R) ∈ R : t 0 D α t V = 0}. Hence, lim t→∞ I (t) = 0, the fractional system (2.4) changes to the equations The solution of (3.6) is

Then we have S(t) = A d and R(t) = 0 as t → ∞. Further, we have already established that the DFE is locally asymptotic stable if R 0 < min{1, R 0 }. Combining both results, we may conclude that the DFE point of the limit set (3.6) is globally stable asymptotically for R 0 < min{1, R 0 }. Hence the proof.

The existence and analysis of Hopf bifurcation for a system is modeled by ordinary differential equations is studied by many researchers [17, 18] . In this portion, we analyze Hopf bifurcation of the fractional-order system (2.4) for different values of α ∈ (0, 1). Here we assume the effectiveness of the saturated treatment control parameter γ and the fractional-order α as bifurcation parameters.

It is known that Hopf bifurcation may take place in spite of stability of a system at the critical value of bifurcation parameter for an integer order system. Therefore, all the roots of the Eq. (3.3) are all real. Hence, a 1 (γ ) = 0 delivers us the necessary critical value γ 0 of the bifurcation parameter

In this case, the Hopf bifurcation occurs at γ 0 if the succeeding conditions are fitted i.e. Here we take α as the bifurcation parameter. Using the Lemma 3.4 [15] , we can say that stability of the fractional system (2.4) depends on sign of ( απ 2 − |arg(λ i )|) i = 1, 2, 3. Based on this observation, we define the function F(α) by

The Lemma 3.4 [15] implies that the EE point will be stable asymptotically if f 2 (α) < 0 and it will be unstable if f 2 (α) > 0. (i) the Jacobian J at the EE point of the system (2.4) induces complex conjugate eigenvalues

Let us assume the incommensurate fractional system (2.4), α i ∈ (0, 1], i = 1, 2, 3 and ζ = 1 L where L = lcm(q 1 , q 2 , q 3 ) and α i = p i q i , and p i , q i are relatively prime for i = 1, 2, 3. Then by Lemma 3.5 [15] , an equilibrium point E(S, I , R) is locally stable asymptotically iff |arg(λ)| > ζ π 2 for all the roots λ's of the characteristic equation

where M = diag(λ Lα 1 , λ Lα 2 , λ Lα 3 ) and J , Jacobian matrix evaluated at the same equilibrium point E(S, I , R).

In the numerical experiments, we apply predictor-corrector P(EC) m E (Predict, multiterm(Evaluate, Correct), Evaluate) method which is the modified method of PECE (Predict, Evaluate, Correct, Evaluate) and Adams-Moulton algorithm [8, 9, 11] in the environment of MATLAB-16 software. For this purpose, we use the solver function as implicit fractional linear multi-step methods (FLMMs) for fractional-order systems.

To perform numerical works, we choose the set of parameters as

These parameters yield the EE E * of the fractional system (2.4) as E * = (5.2 + 10γ 1+I * , I * , (1 + 5 1+I * )I * ) where I * = 0.1(3 − 10γ + 169 − 80γ + 100γ 2 ). E * exists if I * > 0 which implies that the feasible region of γ is 0 < γ < 8. The EE E * of the ordinary system will be stable if a 1 (γ ) = 0 and a 1 (γ ) = 0 gives γ = 0.47 (say, γ 0 ). For γ < γ 0 , the eigenvalues contain negative real parts and hence |arg(λ 2,3 )| > απ 2 . Therefore the EE point is asymptotically stable and converges to E * = (7.343, 1.056, 22.79) for all α ∈ (0, 1). If γ > γ 0 , for α = 1, the EE is unstable as the eigenvalues contain positive real parts. Now, we depict some figures for the solutions of the fractional system (2.4) in the presence of the initial conditions S(0) = 0.045, I (0) = 0.04, R(0) = 0.037. If γ > γ 0 and α ∈ (0, 1], then it can be said that the EE point is stable by Theorem 3.8. If we choose γ = 0.49 then the eigenvalues are λ 2,3 = 0.004 ± 0.257i. When α < 2 π |arg(λ 2,3 )| = 0.987, EE point will be asymptotically stable and it will be unstable whenever α > 0.987. These behaviors are shown in Fig. 3 for γ = 0.44, Fig. 4 for γ = 0.47, and Fig. 5 for γ = 0.49.

Next we assume γ = 0.9, the eigenvalues are 0.068 ± 0.216i. When the value of α < 2 π |arg(λ 2,3 )| = 0.803, EE point will be asymptotically stable and for α > 0.803 the EE point will be unstable. Figure 6 depicts the behavior.

From the above discussion it can be said that with changing values of γ , we can get different region of α for which the EE points remain asymptotically stable. The fractional order system will be asymptotically stable if

is a decreasing function. The stability region in terms of α and γ are shown in Fig. 7 . Now we demonstrate Hopf bifurcation for the fractional system (2.4). To satisfy the singularity condition (i) of Theorem 3.10, we must have a 1 (γ ) > 0 and a 1 (γ ) 2 − a 2 (γ ) < 0 of the equation (3.5) . We see that Hopf bifurcation of the system (2.4) occurs for α ∈ (0, 1) if the parameter γ representing the effectiveness of treatment control lies in the interval (0, 1.5). We choose γ = 0.97 then the critical value of α as α * is determined from condition (ii) of Theorem 3.10 and we get the critical value α * = 0.802. The third condition holds good since, d dα f 2 (α) = π 2 at α = α * . Hence, we see that Hopf bifurcation undergoes at α * = 0.802 and with increasing value of α from α * the orbits of the limit cycles also increase. Figure 8 . depicts such behavior. The eigenvalues at EE of the incommensurate fractional system (2.3) can be determined by the help of the equation (3.7) which is, in general, a higher-order polynomial based on the values of α 1 , α 2 , α 3 . Now we choose γ = 0.95. Then we obtain the EE point (10.89, 0.703, 2.769). We set the fractional orders as α 1 = 0.9, α 2 = 0.7, α 3 = 0.8 and obtain the eigenvalues from equation (3.7). We see that the absolute value of the minimum of arguments of the eigenvalues is 0.133 < πζ 2 = 0.157. Hence, in this case, the eigenvalues are unstable. Again we choose α 1 = 0.7, α 2 = 0.5, α 3 = 0.6 and find the absolute value of the minimum of argument of the eigenvalues of the equation (3.7) as 0.180 > πζ 2 = 0.157. Hence, in this case, the EE point is asymptotically stable. These behavior are presented in Fig. 9 . Now we state the theorem for existence criteria of Hopf Bifurcation due to the parameter γ .

Here we define a function f 3 (γ ) by From Theorem 2 in [7] , we can say that the EE point will be stable asymptotically if f 3 (γ ) < 0 and it will be unstable if f 3 (γ ) > 0. Now we state the theorem for existence criterion of Hopf bifurcation about the parameter γ . 

We now choose γ as a bifurcation parameter. Here also we obtain the eigenvalues from the equation (3.5). The first condition of Theorem 4.1 implies that the eigenvalues obtained from equation (3.5) must contain positive real parts. Then we must have a 1 > 0 and a 2 2 − a 1 < 0. Hence, a Hopf bifurcation of the system (2.4) occurs if γ ∈ (0, 1.5). Thus, the existence of Hopf bifurcation confirms at γ * = 0.587 and with increasing value of γ from γ * the orbits of the limit cycles also increase. Figure 10 . depicts such behavior.

Here we consider the incommensurate system (2.3) where the fractional orders α 1 , α 2 , α 3 ∈ (0, 1) are taken as bifurcation parameter. We study the stability criteria of the fractional system (2.3) by using Theorem 2 in [7] . It shows that the stability depends on the sign of ( αζ 2 − |arg(λ i )|). Here we see that the number of eigenvalues, i depends on the term ζ = 1 L which depends on the different fractional orders (generally, i is greater than the number of equations of the given system). Here we consider the orders of fractional derivative α 1 , α 2 , α 3 as the bifurcation parameters. Now we define a function

From Theorem 2 in [7] , we can see that EE point will be stable asymptotically if w(ζ ) < 0 and it will be unstable if w(ζ ) > 0.

To examine the existence of the Hopf bifurcation for incommensurate fractional system, we state the pursuing theorem. 

).ae = 0, where ae is the directional derivative of the vector curve w(ζ(α 1 , α 2 , α 3 )) = 0 at the critical value of bifurcation parameter

, the critical value of the parameter for fixed γ = 0.97. Let us choose fixed α * 2 = 0.94 and α * 3 = 0.96. Then the equation w(ζ(α * 1 , α * 2 , α * 3 ) = 0 determine (α * 1 , α * 2 , α * 3 ) as (0.692, 0.94, 0.96). At the critical value α * 1 , Fig. 11 establishes that the Hopf bifurcation occurs. Also if α 1 increases beyond the critical value α * 1 , the limit cycles will be an attractor with larger radius.

In this section, we apply the above system (2.1) to predict the COVID-19 situation of India after the second wave. For this, we collect the data of daily (cumulative) active infected of India from 1st March 2021 to 15th August 2021 from the official websites of the Government of India [6] and taking the parameter set as A = 85620, β = 0.084, d = 0.001, = 0.00125, γ = 0.0825, v = 0.1, δ = 0.00157 and α = 0.90. At this environment, to fit the above data with the above model, we use the MATLAB-16 software to drow the prediction graph up to February 2022 which is shown in Fig. 12 and we see that the number of daily active infected will be decreased and at 28th February 2022, the number of active infected will be approximately 19100. From this graph, we can predict the daily active infected of any day between this period of time. disease treatment control get saturated after a certain time. Both commensurate and incommensurate fractional differential equations are used to investigate the complete nature and dynamics of infectious diseases. The uniqueness, positiveness, and uniform boundedness of the solutions of the system (2.4) have been established. The threshold parameter to discuss the infectious disease dynamics is commonly known as basic reproduction number (R 0 ) is derived, and it has also been confirmed that the DFE is locally and also globally stable when the threshold condition R 0 < min{1, R 0 } holds. On the other hand, the system also passes through a backward bifurcation at the EE point for the threshold range R 0 < R 0 ≤ 1. Moreover, several computer simulations works enrich our theoretical results. A stability region for the corresponding fractional system is identified. Along with the Hopf bifurcation phenomenon considering the fractional derivative i.e., α as bifurcation parameter is visualized graphically. We have further observed that with the increasing value of α, the length of the cycle increases, and this phenomenon established that by increasing the importance of α, the system repels away from the equilibrium point. Moreover, the Hopf bifurcation scenario for the incommensurate fractional system is also derived with bifurcation parameters α 1 , α 2 , and α 3 . We also experienced the significance of α with the help of the model fitting. Also, we have drow a prediction curve of daily active infected of India from 1st March 2021 to 28 February 2022, i.e., of one year with help of the existing data from 1st March 2021 to 15th August 2021. From this prediction curve, we can conclude that if the present situation (i.e., the value of the parameter that we have assumed here) does not change, there will be no chance of a third wave before February 2022.

Computational algorithm for solving singular Fredholm time-fractional partial integrodifferential equations with error estimates

Modeling behavioral change and COVID-19 containment in Mexico: a trade-off between lockdown and compliance

A mathematical model of Coronavirus Disease (COVID-19) containing asymptomatic and symptomatic classes

The RKHS method for numerical treatment for integrodifferential algebraic systems of temporal two-point BVPs

Backward bifurcations in simple vaccination models

Dynamic analysis of fractional-order predator-prey biological economic system with Holling type II functional response

A predictor-corrector approach for the numerical solution of fractional differential equations

Braunschweig: efficient solution of multi-term fractional differential equations using P(EC) m E methods

Measuring memory with the order of fractional derivative

On linear stability of predictor-corrector algorithms for fractional differential equations

Front dynamics in fractional-order epidemic models

Chaos in a fractional order Chua's system

Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts

Analysis of a fractional-order SIS epidemic model with saturated treatment

Complex dynamics of an SIR epidemic model with saturated incidence rate and treatment

Application of three controls optimally in vector borne disease-a mathematical study

Stability analysis of a prey-predator model incorporating a prey refuge

Contribution to the mathematical theory of epidemics

Mathematical modelling of infectious diseases

Early dynamics of transmission and control of COVID-19: a mathematical modelling study

Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and and generalized Mittag-Leffler stability

Hopf bifurcation analysis of a new commensurate fractional-order hyperchaotic system

A domain decomposition approach to a SIR epidemic model with constant vaccination strategy

A model based study on the dynamics of COVID-19: prediction and control

Modelling and control of a fractional-order epidemic model with fear effect

Modelling and control of COVID-19: a short term forecasting in the context of India

Piecewise optimal fractional reproducing Kernel solution and convergence analysis for the Atangana-Baleanu-Caputo model of the Lienard's equation

The reproducing kernel algorithm for numerical solution of Van der Pol damping model in view of the Atangana-Baleanu fractional approach

Fractional-Order Nonlinear Systems: Modelling Analysis and Simulation

Fractional Differential Equations

The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study

The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study

Mathematical modeling based study and prediction of COVID-19 epidemic dissemination under the impact of lockdown in India

On the appearance of the fractional derivative in the behavior of real materials

Knowledge epidemics and population dynamics models for describing idea diffusion

Estimating the unreported number of novel coronavirus (2019-nCoV) cases in China in the first half of January 2020: a data-driven modelling analysis of the early outbreak

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

 19-02-2019) . Moreover, the authors are very much thankful to the anonymous reviewers and Chin-Hong Park, the editor in chief of the journal, for their constructive comments and helpful suggestions to improve both the quality and presentation of the manuscript significantly.