Fractional SIAQR model with time dependent infection rate

Abstract In this study we have proposed the SIAQR model with time-dependent infection rate. The model’s well-posedness is demonstrated, and the proposed model is then extended to the fractional mathematical model by employing the Caputo-Fabrizio fractional derivative operator. We have also used the Lipschitz condition and the linear growth condition to determine the conditions under which the model has a unique solution. The numerical solutions are presented. Furthermore, using graphical representations we will see how the number of infected, recovered individuals vary as the order of fractional derivative varies.


Introduction
As we know that mathematical models are capable of decision-making, saving lives, assisting in policy, and many more.These are helpful in understanding the conditions needed to sustain lives and provide us with ways to study and predict the behavior of the spread.The concept of derivatives and integrals plays a lot in the formulation of these mathematical models.In this work, we will be concerned with the SIAQR mathematical model.A SIAQR model is an epidemiological model that determines the number of people infected with a transmissible infection in a closed population over a period of time.This class of model derives from the fact that they involve equations relating the number of susceptible people S(t), the number of people infected symptomatically I(t), the number of people who are asymptomatically infected A(t), the number of people who have been quarantined Q(t), and who have recovered R(t).Recently the concept of fractional derivatives (Kilbas et al., 2006;Miller and Ross, 1993;Oldham and Spanier, 1974) and integrals are proved to be very useful in making a better understanding of realworld problems which exhibit non-local behaviors.The first and most widely accepted definition of the singular kernel was proposed by Riemann-Liouville, later on, many new developments were made in the field of fractional calculus which involves the definition of the fractional derivative with the non-singular kernel, one of which was given by Caputo.This definition is based on the concept of the power law.The literature is full-fledged with many definitions of fractional derivatives and integrals.Applications can be seen in the literature (Ahmad et al., 2021;Ain & He, 2019;Ain et al., 2021Ain et al., , 2022aAin et al., , 2022bAin et al., , 2022c;;Alkahtani & Alzaid, 2020, 2021;Alshomrani et al., 2021;Al-Smadi et al., 2021;Anjum et al., 2021aAnjum et al., , 2021b;;Atangana & Aguilar, 2018;Atangana & Araz, 2021;Atangana & Atangana, 2020;Djeddi et al., 2020;Hasan et al., 2022;Jarad et al., 2018;Karaji & Nyamoradi, 2020;Khan & Atangana, 2020;Modi et al., 2021;Momani et al., 2021;Naik et al., 2020Naik et al., , 2021Naik et al., , 2022;;Ucar et al., 2019).In this article, we consider the SIAQR model with time-dependent infection rate.An analysis of this SIAQR model is presented in detail.Using the Caputo definition of fractional derivative, we will study the fractional SIAQR model.In Section 1, an introduction along with some definitions of the fractional derivatives are given.In Section 2, we have presented the SIAQR model with a time-dependent infection rate.In Section 3, we have presented the existence and the uniqueness of the system of solution of the SIAQR model.In Section 4, the SIAQR model with time-dependent infection rate in sense of the Caputo fractional derivative operator is given and numerical simulations are performed for a better understanding of the model.
Following are some definitions of fractional derivatives and integrals.
Definition 1.For an integrable function f on R, c > 0 the Riemann-Liouville fractional integral of order c is given as where RL I c t denotes the Riemann-Liouville fractional integral of order c (Miller and Ross, 1993).
Definition 2. For an integrable function f on R, 0 < c < 1 the Riemann-Liouville fractional derivative of order c is given as where RL v c t denotes the Riemann-Liouville fractional derivative of order c (Miller and Ross, 1993).
Definition 3.For an integrable function f on R, 0 < c < 1 the Caputo fractional derivative of order c is given as where C v c t denotes the Caputo fractional derivative of order c (Miller and Ross, 1993).

SIAQR model with time dependent infection rate
In the SIAQR model, we have considered the five classes of individuals as S: where, d defines the birth and natural death rate, r 1 , r 2 , r 3 , represents the recovery rates for infected symptomatically, asymptomatically infected, and quarantined groups of individuals respectively, q is the rate by which infected people are quarantined.b 1 and b 2 are the infection rate for the individuals infected symptomatic and asymptomatically.Further, a time dependence on the infection rate b 1 and b 2 is introduced with an exponential decay behaviour e Àzt , where, z is the rate of decrease.

SIAQR model with time-dependent infection rate in sense of Caputo fractional derivative operator
We now extend the above formulated SIAQR model to the fractional SIAQR model in sense of Caputo fractional derivative, which is given as follows: Here C v c t denotes the Caputo fractional derivative of order c: Theorem 2. Assuming that all the parameters defined in the above formulated model are positive, the solution of the model SðtÞ, IðtÞ, AðtÞ, QðtÞ, RðtÞ with nonnegative initial conditions are positive 8t > 0: Proof.We consider the first equation On solving, we get SðtÞ ¼ C 1 e Àdt , where C 1 is constant of integration.
Using the fact that Sð0Þ ¼ S 0 > 0, we get In a similar way, we consider the second equation On solving, we get IðtÞ ¼ C 2 e Àðr 1 þqþdÞt , where C 2 is constant of integration.
Using the fact that Ið0Þ ¼ I 0 > 0, we get On a similar account, we get

Existence and uniqueness of the fractional SIAQR model
In this section we will prove the existence and the uniqueness for the SIAQR model.For convenience, we write the SIAQR model as (4.2) To prove the existence and the uniqueness of the solution of the SIAQR model, we use the concept which was recently proposed by Atangana (2020) and hence prove the following theorem.
Theorem 3. Assuming that there exists positive constants L 1 , L 2 , L 3 , L 4 , L 5 and L 1 , L 2 , L 3 , L 4 , L 5 such that the following holds i.
jg 1 ðt, SðtÞÞ À g 1 ðt, S 1 ðtÞÞj L 1 jS À S 1 j, jg 2 ðt, IðtÞÞ À g 2 ðt, I 1 ðtÞÞj L 2 jI À I 1 j, Define the norm as jjwjj 1 ¼ Sup t2½0, T jwj, we get Again defining the norm as jjwjj 1 ¼ Sup t2½0, T jwj, we get Next consider, On a similar account, we get where Similarly, we can show and jg 5 ðt, RÞ À g 5 ðt, R 1 Þj ¼ L 5 jR À R 1 j: (4.14) We now prove the second part of the above stated theorem.We first show jg We now show On taking L 2 ¼ jjNjj 2 1 , and under the condition On taking L 3 ¼ jjNjj 2 1 , and under the condition On taking L 4 ¼ 3q 2 jjIjj 2 1 þ 3q 2 jjAjj 2 1 and under the condition and under the condition We get Hence, by using Theorem 1, the existence and uniqueness of the SIAQR model is proved.

Derivation of numerical scheme for
Caputo fractional SIAQR model using Atangana-Seda method
where g is a non-linear function.
In Figure 2, we have presented the numerical simulation at c ¼ 0:7: We see that as the order of fractional derivative increases from c ¼ 0:5 to c ¼ 0:7, the decrease in the number of individuals who are infected symptomatically decreases well and the number of quarantined individuals increases.
In Figure 3, we have presented the numerical simulation at c ¼ 0:899 and z ¼ 0.001.We see that as the order of fractional derivative increases from c ¼ 0:7 to c ¼ 0:899, the decrease in the decrease in the number of individuals who are symptomatically infected is much greater compared to c ¼ 0:5 and c ¼ 0:7:: Furthermore, we can see that the number of asymptomatically infected individuals is also decreased.
In Figure 4, we have presented the numerical simulation at c ¼ 0:5 and z ¼ 0.01.We see that as the value of z increases from 0.001 to 0.01, the decrease in the decrease in the number of  In Figure 5, we have presented the numerical simulation at c ¼ 0:7 and z ¼ 0.01.We see that as the value of z increases from 0.001 to 0.01.The number of recovered individuals increases well and the number of asymptomatically infected individuals is quite small.In Figure 6, we have presented the numerical simulation at c ¼ 0:899 and z ¼ 0.01.We see that as the value of z increases from 0.001 to 0.01.The number of recovered individuals is much more than the number of individuals who are symptomatically and asymptomatically infected.In Figure 7, we have presented the numerical simulation at c ¼ 0:899 and z ¼ 0.1.We see that due to the effect of the parameter z (change in infection time), there is a sharp increase in the number of quarantined individuals, and hence the recovery rate is also much higher.

Conclusion
In this work, we have formulated the SIAQR model for Covid pestilence with a time-dependent infection rate.The SIAQR model is then extended to the fractional mathematical model within the framework of the Caputo derivative operator.The existence and uniqueness of the solution of the SIAQR system is carried out using the Lipschitz condition and the linear growth condition.The numerical scheme for the Caputo fractional SIAQR model is derived using the newly proposed Atangana-Seda method.Numerical simulations are performed for the SIAQR model with a time-dependent infection rate for different fractional order derivatives.Furthermore, we see that as the value of z increases, the number of infected individuals decreases.This is because the rate of infection is lower as the value of z increases.We also see that as the order of fractional derivatives increases, the number of infected populations decreases, and the number of recovered populations increases.Hence, the fractional derivatives provide us with a better understanding of the model.