Comparative analysis of classical and Caputo models for COVID-19 spread: vaccination and stability assessment

Several epidemiological models use the Caputo fractional-order diﬀerential operator without establishing its signiﬁcance. This study veriﬁes a Caputo operator-based fractional-order epidemiological model of the SAIVR type. COVID-19 kills. Infection weakens the immune system. The fractional Caputo operator describes COVID-19 immunization. Fundamental system characteristics are determined using fractional calculus. Our analysis included the fractional system’s Hyers–Ulam–Rassias stability and stable states. The uniqueness and existence of fractional Caputo system solutions are explored. The least-squares approach determines system parameters. The Caputo fractional-order α value is optimized to 6.757e–01, indicating that the system best ﬁts real-life medical data for infection. Caputo and classical systems were compared for absolute mean errors. The Box-Whisker chart case summaries show the Caputo operator superiority. When α → 1, the memory traces and hereditary traits are also observed. Finally, the Caputo fractional framework simulates COVID-19 using strong numerical methods.


Introduction
The World Health Organization (WHO) called COVID-19 a pandemic after it spread to many countries on all continents in the first quarter of 2020.Different mutations of the new coronavirus have been observed in different regions, resulting in a wide range of symptoms.Fever, dry cough, and fatigue are the most common signs of illness, but other signs, such as a sore throat, diarrhea, loss of taste or smell, or a rash, are less prevalent.COVID-19 can spread through direct contact, indirect contact, droplet spray (like when someone sneezes), short-range transmission, airborne transmission (through aerosol), and long-range transmission.
Integer-order differential equations have been extensively explored for their application in epidemiological models of infectious diseases.Evidence from mathematical modeling of epidemics in the literature suggests that nonlinear dynamical equations can shed light on disease transmission dynamics.Constructing realistic nonlinear compartmental math-ematical models that are data-driven to better explain the transmission dynamics of epidemics has garnered much attention in the wake of the recent global COVID-19 outbreaks [1][2][3][4][5][6][7][8], including several other important research studies on mathematical modeling with the Caputo and other memory operators, as discussed in [9][10][11][12].
Musa et al. [13] propose a novel eight-compartmental deterministic model for the COVID-19 epidemic in Nigeria, which accounts for awareness campaign initiatives and hospitalization control tactics for severe and moderate instances of infections.Infected cumulative instances in Nigeria were used to fit the suggested nonlinear dynamical model from 29 March 2020 to 12 June 2020.Their research shows that an increase in infections could happen in the population if programs to raise awareness are not done well.Some authors in [14] developed and tested a SEQIJR epidemic model that looks at how well quarantine and isolation work to stop COVID-19 outbreaks in Pakistan.
Caputo fractional differential operators are used in mathematical epidemiology to model the dynamics of infectious diseases in populations where memory effects play an important role.Unlike classical differential equations that assume the instantaneous response to perturbations, fractional calculus considers the possibility of memory effects, meaning that the system's response to a perturbation depends on the entire history of the system.In epidemiology, Caputo fractional differential operators have been used to model the spread of infectious diseases that have a long-lasting impact on individuals, such as measles or rubella.The Caputo fractional derivative of a function can be interpreted as the fractional order derivative of the function's memory, which can be used to capture the long-term effects of infection on an individual's immune system.Caputo fractional differential operators have also been used to model the effect of vaccination campaigns, where individuals who have been vaccinated have a lower probability of becoming infected in the future.Using fractional calculus, it is possible to model the long-term effects of vaccination on the dynamics of the disease, which can help inform public health policy decisions.Overall, the use of Caputo fractional differential operators in mathematical epidemiology allows for a more accurate representation of the dynamics of infectious diseases, taking into account the complex interplay between memory effects, disease transmission, and control measures.When modeling biological and engineering processes mathematically, fractional order differential equations are extremely valuable and powerful tools.This is due to the fact that the differential operators in such equations or models are associated with systems with memory dynamics, a property shared by the vast majority of biological and technical systems [15,16].To solve the problem of fractional-order models with control functions that change over time, authors in [17] develop a new and effective numerical method based on hybrid Chelyshkov functions.
The fractional calculus Caputo-Fabrizio derivative operator was used by Rezapour and coauthors in [18] to expand the nonlinear integer-order anthrax illness model established and assessed by Githire et al. [19].The existence criterion of solutions was presented for the suggested fractional-order anthrax disease epidemic model using the Picard-Lindelof technique.Several infectious disease transmission dynamics, including those of HIV/AIDS [20], tuberculosis [21], malaria [22], dengue fever [23], Zika [24], Ebola [25], and hepatitis B [26], have been studied and analyzed using differential equations characterized by the Caputo fractional-order derivative.
In [27], the dynamics of reverse bifurcation are investigated in a simple but realistic model of a vaccine pandemic.To back up the theoretical results, the author provides qual-itative and numerical simulations of the formulated mathematical model.The mathematical model of the cholera outbreak that the author developed in [28] was extended and investigated by Javidi et al. [29] to capture the Caputo fractional-order derivative.In [30], the author examined the uniform asymptotic stability of a few simple epidemic models (SIS, SIR, SIRS), as well as the well-known Ross vector-borne diseases in the Caputo sense using Lyapunov functions of the Volterra type.
In [31], the Caputo fractional derivative is used to analyze the fundamental mathematical model of SEIR, predicated on stochastic population dynamics.The authors provided an in-depth qualitative stability investigation of their novel and plausible deterministic model.The authors in [32] used a nonlinear system of differential equations in the sense of the Caputo fractional order derivative operator to form an epidemiological model for Zika virus infection, dividing the total human and mosquito populations into two compartmental classes (susceptible people and infected people; susceptible mosquitoes and infected mosquitoes).The authors of [33] came up with and studied an SEIR-type epidemic model based on classical and Caputo fractional-order differential operators to describe how the Rubella epidemic in Pakistan changed over time.In [34], a fractional system models liquid surface-stripe wave evolution.The Riemann-Liouville fractional derivative solves the fractional system analytically.
Recent research has attempted to mathematically model this devastating COVID-19 pandemic using some of these valuable differentiable operators, which are, in turn, derived from fractional-order differential equations [35].Parameter estimation and numerical simulations for a nonlinear COVID-19 epidemiological model built with Caputo and Atangana-Baleanu fractional derivative operators were supplied by Naik et al. [36], along with a comprehensive qualitative analysis.Another study [37] used a nonlinear Atangana-Baleanu fractional-order differential equation model to analyze the COVID-19 pandemic in Nigeria.The latest research by Baleanu et al. [38] explored a Caputo-Fabrizio derivative version of the integer-order epidemic model introduced and studied by Chen et al. [39].The uniqueness of the solution to the nonlinear Caputo-Fabrizio fractional order COVID-19 model was demonstrated using fixed point theory.They used the transform method of homotopy analysis and developed a convergent series approximation to the model problem.In [40], a Caputo fractional order deterministic epidemic model for COVID-19 infection was constructed and analyzed.To prove that there is a unique solution to the mathematical model, they turned to the well-known Banach contraction mapping concept.To add the importance of fractional theory, authors in [41] solve the (2 + 1)-dimensional elliptic nonlinear Schrodinger problem with three fractional operators and derive fractional analytical solutions.
It is difficult to predict how the present pandemic may influence people's choices to enroll in a COVID-19 vaccine clinical trial or to get immunized against the COVID-19 vaccine.France, which has the highest number of people who refuse to get vaccinated, is particularly concerned about this.The Centers for Disease Control and Prevention (CDC) classifies four of the top five counties in the United States with the highest percentage of totally vaccinated people (84.3%) as "high" transmission jurisdictions.Many countries consider vaccinated people unlikely to be a source of disease transmission.When thinking about ways to manage public health, it seems like it would be irresponsible not to consider the vaccinated population as a possible and major source of transmission.
In terms of the reasons for the lower COVID-19, the majority of research (25/35 = 71.4%)found that the male gender was an enabling factor in vaccine apprehension and increased desire for COVID-19 vaccinations.When subjects were older or had a PhD or higher education, the acceptance rate of the COVID-19 vaccination was higher in the 23 trials (23/35 = 65.7%).Vaccination was approved in 51.4 percent of the studies (18/35 = 51.4%).In addition, white or Asian ethnicity, higher income/education, history of chronic illness, familiarity with COVID-19 infection and sickness, employment in non-rural regions, and confidence that vaccinations may protect family and community members were all associated with greater acceptance of the COVID-19 immunization.
The present literature's use of fractional-order operators provided the impetus for developing an epidemiological model of the COVID-19 infection, which includes a class for the vaccinated population.This paper proposes an alternate SAIVR model to the model originally proposed in [42] under classical calculus, in which the spread of COVID-19 occurs directly between one class of vulnerable individuals and the infected patients being vaccinated against.The model uses the Caputo fractional-order derivative and is based on the coupling of nonlinear ordinary differential equations with real-world properties.
The following outline constitutes the paper's structure: The mathematical construction of a SAIVR nonlinear dynamical system is carried out in Sect. 2 with some explanations for the usage of the Caputo fractional derivative operator in Sect.3. In Sect.4, the biological parameters of the model are estimated.The qualitative analysis of the nonlinear model is discussed in detail in Sect.5, which includes an examination of the existence, uniqueness, and steady states of the Caputo model.The simulations of the model are included in Sect.6. Section 7 gives the conclusion and future research directions.

Mathematical formulation
Recently, in [42], a standard epidemic model for explaining COVID-19 transmission has been developed using the first-order time derivatives from classical calculus.The vaccinated SAIVR nonlinear dynamical model is a deterministic type of system based on five compartments whose flow chart is in Fig. 1, wherein the dynamics for the movement of individuals can easily be comprehended.The model is presented as follows: The Caputo operator, a typical fractional-order operator, is used to analyze the integerorder deterministic model.The COVID-19 epidemic SAIVR model is fractionally evaluated using the Caputo-type operator under the premise that solutions exist and that the operator utilized for nonlinear differential equations is of this kind.As a result of its widespread application in epidemiological modeling, the Caputo fractional operator has gained popularity.The SAIVR model in a deterministic system has been separated into five categories: susceptible, asymptomatic/undetected, infected, vaccinated, and recovered individuals at any point in time t.Nonlinear differential equations are shown below under the Caputo operator, with α being the fractional-order operator in the Caputo sense.
In this study, we examine how fractional calculus in mathematical epidemiology can be used to create a model of COVID-19, among other areas.In fact, accurate modeling of those infectious disease systems makes use of fractional derivative formulations.In any case, fractional order differential equation systems provide an interesting modeling technique in the context of epidemiology, as one could expect, given that they permit higher degrees of freedom and incorporate memory effects in the model.In fact, the fractionalorder derivative is a good tool for expressing the memory and heredity features of many materials and processes, making fractional differential equations more sufficient for de-scribing phenomena related to non-locality than integer-order derivatives.Therefore, epidemic-based fractional derivatives have also been employed to address some epidemic tendencies.In general, this basic/classical model does not yield adequate or satisfactory results, as demonstrated by the failure of the classical first-order differential equations to accurately reproduce the statistical data collected during an actual epidemic of the disease.In our work, we have explored a more precise and intricate set of fractional differential equations in an effort to achieve better findings that are more in line with reality.

Why COVID-19 model with the Caputo operator?
Traditional models employ local differential and integral operators, which ignore the details of the epidemic being studied.Therefore, the memory features of the underlying system are ignored by conventional calculus.Some recent studies have shown that nonlocal operators are superior to classical ones, making them the only choice for incorporating memory effects into the deterministic model of the epidemic [43].Recent scientific literature has seen the proposal of multiple epidemic models for the spread of the COVID-19 virus.Several have been created with the help of fractional-order differential operators, the most popular of which is called Caputo.However, most research articles practically never explain why the Caputo operator was used.In this paper, we have attempted to describe some of the more salient reasons why the Caputo operator should be considered when modeling an infectious disease with nonlinear differential equations.The following are some of the explanations: • Because it is common knowledge that infectious diseases contain hereditary components, the most appropriate alternatives for modeling infectious diseases are fractional-order operators that can potentially preserve the memory of the systems being considered [44].In this way, the Caputo operator is the best choice to replace the integer-order time derivative in the COVID-19 model that is currently being talked about.• Substituting α ∈ C for n ∈ N in the Cauchy formula for repeated integration yields the well-known Riemann-Liouville integral formula, which has served as the basis for developing numerous numerical methods for solving fractional ordinary and partial differential equations.• Recent works have successfully used Caputo's revision of traditional epidemiological models, which are paired with particulars regarding the existence of a single solution and a study of the system's stability [45].In that part, simulations were used to show why the Caputo variant is better than the standard (classical: integer-order derivatives) one.• Epidemiological research published recently [46] demonstrates that older models failed to account for the complex and unpredictable dynamics of an infection's spread.Instead, we used actual data regarding the outbreak, largely from reputable sources like the World Health Organization and scientifically published articles, to confirm and corroborate the Caputo accounts.To further clarify the basic reproductive number, which describes the typical number of secondary infected cases produced when an infectious individual enters a completely susceptible class, we can use Caputo's differential operator to analyze the disease's behavior under different values for biological parameters.
• Finally, using a fractional-order operator for any physical or biological model, one can verify that there is no physical or geometrical meaning of the fractional parameter α while keeping in mind that there is no physical or geometrical meaning of even integer-order derivatives, such as the fourth-order ODE for the deflection of a beam [47] wherein the fourth derivative itself does not have any meaning.

Best fitting of biological parameters
An epidemic's dynamics can be explained by a system of nonlinear ordinary differential equations, but there are challenges associated with validating the suggested differential equations and achieving appropriate values for the working (or biological) parameters of the proposed equations.Since most parameters are fully indeterminate from the available information about the nonlinear system under study, the models with parameters that need to be optimized and fitted are spread out across several domains.Infectious diseases in mathematical epidemiology are typically represented using a set of nonlinear differential equations with a number of continuous parameters whose fitted values are crucial to get using various methods currently accessible in the literature.However, with the aid of demography and the studied disease's past history, several of the biological parameters connected with the suggested epidemic can be simply computed.It is typical practice to utilize the biological parameters by employing the values discovered through study or making educated assumptions about them.However, this method can sometimes make the disease act in ways that are hard to predict, which means that the proposed model cannot be tested.

Abstract analysis
In this section, a detailed qualitative analysis of the fractional-order Caputo model ( 2) has been carried out wherein major focus is placed upon the existence and uniqueness of solutions of the system, including its stability analysis via Hyers-Ulam-Rassias stability criterion.

Existence and uniqueness results
In the context of this study, the Banach contraction principle supports the uniqueness of the solution by providing a mathematical foundation to prove that the proposed Caputo SAIVR model converges to a single, unique solution, indicating the reliability and stability of the model in predicting the dynamics of the COVID-19 pandemic.This section will discuss the Lipchitz condition, existence, uniqueness, and stability of the model (2).As a result, we first assume the following five kernel values for simplicity and clarity: ( Applying the fractional-order Caputo operator on the above five kernels given in (3), where (α) and (α) are positive real constants.Now, we will prove the Lipchitz condition for the Caputo system (2).
Proof First, the Lipchitz condition is justified for 1 kernel.Take S and S * as two functions, then the corresponding norm is as follows: Simplifying and applying the norm property, we get taking κ 1 = where N and δ α N are bounded functions.Similarly, the norms can be created for rest of the model equations.

Existence of the solution
In this subsection, we will prove that the Caputo model (2) under consideration has at least one solution.Thus, the recursive formula for (4) becomes the one shown below: The positive initial conditions are first iterative values.The difference between two consecutive terms is as follows It is worth noticing that Taking first equation of system (8), we assess the following Applying the triangular inequality reduces the above equation to The kernel 1 (t, S) satisfies the Lipchitz condition on the evidence of Eq. ( 6).Therefore, we can write it as follows On the evidence of (9), we can reduce the above inequality in the following manner Similarly, we can get the following results: Theorem 2 The analytical solution exists for the fractional-order Caputo model (2) under the given condition at 0 : Proof Since the functions S(t), A(t), I(t), V (t), R(t) are bounded and the condition for the Lipschitz is satisfied, using Eq. ( 13) and the recursive relation yields As a result, the solutions above will continue to exist.On the other hand, to show that the aforementioned functions represent the suggested model's solution, we take into account Therefore, we have Using the Lipchitz condition, This gives, Then, at t 0 , we have As n tends to ∞, we attain Similarly, we can derive This justifies the existence of the solution.

Uniqueness of the solution
Now, we must demonstrate whether or not the answer is unique.So, suppose that another solution to the suggested model exists and it is S(t), Equation ( 22) with norm, Theorem 3 The analytical solution is unique for the Caputo fractional model under the following condition that is Proof Take into account that (24) holds so that from ( 23) Hence, we can say that S(t) -S * (t) = 0.It implies that S(t) = S * (t) and the solution is unique.We now apply the same approach to another function yielding the following results Thus, this proof shows that the proposed version of model ( 2) in the sense of the Caputo operator has a unique solution.

Hyers-Ulam-Rassias stability
The Hyers-Ulam-Rassias stability is a concept in functional analysis, particularly in the theory of functional equations.It refers to the stability of a functional equation under small perturbations of its arguments.In other words, it measures how close a solution of a functional equation is to be a solution of a perturbed version of that equation.The stability of a functional equation is a useful property for proving the existence and uniqueness of solutions, and it has important applications in a wide range of fields, including mathematics, physics, engineering, and economics.This subsection deals with stability analysis of fractional model ( 2) under the concepts of Hyers-Ulam-Rassias stability analysis.Let us rewrite model (3) as follows: where, the vector ϕ = {S, A, I, V , R} and = ( 1 , 2 , 3 , 4 , 5 ) is a continuous vector function.
Proof Let ϕ ∈ P([0, T], R 5 ) be a solution of model (27).Then, the unique solution of model (27) from Theorem (3) is On the evidence of (30), we can say that The Gronwall's inequality yields On setting Inequality (35) authenticated that model ( 27) is generalized Hyers-Ulam-Rassias stable with respect to .

Numerical results from simulations
This section analyzes the numerical dynamics of the COVID-19 epidemic as modeled with the Caputo fractional operator.When doing so, it is important to keep in mind that the fractional order parameter that is used for simulations is the one that was optimized in the Sect.4, while the rest of the parameters are taken from the third column of the Table 1.An explicit numerical method of the predictor-corrector type that was developed [48] specifically for use in simulations of fractional Caputo types of ordinary differential equations is used in this particular instance.One may find the method itself as well as a comprehensive study of it based upon convergence and error analysis in [49].The numerical method developed specifically for this objective has garnered much praise in the most recent body of research because of its ease of use and adaptability.Furthermore, the author of [50] has thoroughly discussed the numerical approach in MATLAB implementation, including predictor-corrector capabilities.The simulations for the present research work for the Caputo COVID-19 model (2) were significantly simplified by utilizing such readily available routines on MathWorks.It should be noted that we run the MATLAB software on Windows with an Intel(R) Core(TM) i7-1065G7 CPU running at 1.30 GHz, 1.50 GHz, and 24 GB of RAM.The version number of this software is "9.8.0.1323502 (R2020a)." For the COVID-19 model selected to be operated on by the Caputo operator, there are many crucial parameters whose values need careful attention to determine whether they are growing or shrinking.In this regard, we have selected a few factors, such as the individuals' exposure rate β 1 , the asymptomatic infection rate β 2 , and the Caputo fractional order optimized parameter α.We have also shown the numerical dynamics of each of the five state variables using the Caputo operator, taking into account the fitted parameters and the value of α that has been optimized with the least-squares technique.
Several observations have been produced as a consequence of running simulations using the Caputo version of the COVID-19 model.Figure 5 shows that the dynamical behavior of the four classes (susceptible, asymptotic, vaccinated, and recovered) reflects the theoretical analysis that is being observed in the community after it was simulated using the four state variables (susceptible, asymptotic, vaccinated, and recovered).The simulations have been carried out with the parameters that best match the Caputo model.If the parameters shown in the table are used, the number of individuals who are susceptible to the disease and those who have recovered from it will both increase with the passage of time, while the number of individuals in the asymptotic and vaccinated groups will begin to decrease.In addition, it can be seen from Fig. 6 that if we slightly increase the exposure rate, then the infection substantially rises, demonstrating the sensitivity of this fundamental parameter.This is shown to be the case.This demonstrates that if one can manage the number of people exposed to the infection, one can also control the infection.The behav- ior of a comparable nature is seen to occur with regard to the infection rate, as shown in Fig. 7.
Last, the fractional-order parameter is important to consider while comprehending the disease transmission pattern.As can be seen in Fig. 8, it is not advised to use values of this parameter that are lower than those shown, and determining the parameter's optimal value is of the utmost significance.In the current research investigation, this particular procedure was carried out.It should also be mentioned that quite a few studies have been conducted that have carried out this kind of analysis, which finds the optimal value of the fractional-order parameter in the Caputo sense.This is something to consider.

Memory trace and hereditary traits
Using fractional-order differentiation to represent memory effects in the system, the Caputo fractional-order epidemic model is a mathematical model used to predict the transmission of infectious illnesses [51].The term "memory trace" is used here to describe how previous infections have affected the current pandemic.To account for genetic factors that may influence the dynamics of illness transmission, the model can be expanded to include hereditary features.Fractional-order differential equations (FODE), in which the derivative order is not an integer but a fractional value, are used in epidemic models to add memory effects.This helps the model more accurately reflect the course of the pandemic.Fractional derivatives capture the effect of prior infections in the setting of the Caputo fractional-order epidemic model, where fractional-order differential equations characterize the dynamics of the epidemic.Therefore, the current status of the disease is dependent not only on the current affected population but also on its past levels.This model's memory trace can be used to capture phenomena like the population's slow but steady development of immunity, the disease's tenacity, and the lasting effects of therapies.
Traits that can be passed down from one generation to the next are said to be hereditary.These characteristics are relevant to epidemiology because they can influence an individual's vulnerability to infection, transmission rate, and prognosis.Incorporating genetic elements into the differential equations that regulate the dynamics of the epidemic is what is meant by including hereditary qualities in the Caputo fractional-order epidemic model.Transmission rates, recovery rates, and contact rates are just some of the model parameters that might be affected by genetic factors.Disease dynamics may be drastically altered by inherited characteristics.Whether an individual contracts a disease, is able to transfer it to others, or benefits from a vaccination program may all depend on their genetic makeup.Integrating genetic data and population genetics principles into an epidemiological model to account for heritable features is a hard endeavor.In conclusion, the Caputo fractionalorder epidemic model is an effective tool for modeling the role of memory in the spread of disease.It takes into account the impact of previous infections on the current epidemic by employing fractional-order derivatives.Genetic factors influencing illness transmission and susceptibility can also be accounted for by including hereditary features in the model.These two features can be used to develop more precise and all-encompassing models of infectious illness dynamics in communities.
To delve into the behavior of model ( 2), we utilize the Caputo operator defined in [52] for our analysis.For α, 0 < α ≤ 1 derivative, let the fractional derivative of variable χ(t) be Using one of the most widely used numerical methods, namely, the L1 scheme [52], the numerical approximation of the fractional-order derivative (FOD) of χ(t) is as follows: where dt = h stands for the stepsize.One of the most effective numerical methods for discretizing the Caputo fractional-order derivative (CFOD) in the time domain is the L1 scheme.Although the memory component is present in other numerical approaches, its integration is more clearly represented in the L1 scheme.Taking ( 36) and ( 37), we get the following numerical solution to (36): - Therefore, the difference between the Markov term and the memory trace [53,54] can be thought of as the FODE solution.Here is how the Gamma function affects the Markov term: The memory trace (MT) (χ -memory trace since it is related to variable χ(t)) is The memory is adept at combining all prior acts, and this includes the system's tremendous historical evolution.When α = 1, the memory trace is zero at all times t.The behavior of memory traces changes significantly over time.Nonlinearly increasing from zero, the memory trace increases when α is decreased.Therefore, fractional-order systems behave considerably differently from integer systems.We now present numerical simulations and elaborate biological interpretations of memory traces.To achieve this, the numerical approximation of the fractional-order derivative of S(t) is as follows: Using (41) and the first compartment of the Caputo fractional system (2), the numerical solution of Susceptible individuals S(t) is given by: By following the same steps, the numerical approximations of the fractional-order derivative of A(t), I(t), V (t), and R(t) can be achieved.Numerical simulations were conducted to visually illustrate the influence of memory trace on specific sub-populations inside the Caputo fractional system, as defined by Equation ( 2), employing the aforementioned methodology.Figure 9 illustrates the impact of memory trace on population dynamics across a range of fractional-order α values.Based on the observations depicted in these plots, it can be inferred that the absence of a memory effect is evident when the value of α is set to 1.As the value of α declines from 1 to 0.7, it becomes evident that the fractional order and the presence of a memory effect have discernible effects.The phenomenon of the memory effect has the potential to yield accurate outcomes and predictions pertaining to the COVID-19 pandemic.Therefore, the influence of memory plays a crucial role in epidemiological models.After doing an analysis in Fig. 9, it becomes evident that the memory effect tends to approach zero after a certain period of time.The observed results are consistent with the expected outcomes of real biological mechanisms occurring within the human organism.The results obtained from the graphical representations suggest that FODEs successfully capture the memory impact of the system, eliminating the need for additional components.It is widely recognized that fractional-order derivatives are favored due to their inherent memory effect.The activation of the memory effect plays a crucial role in the system effectiveness.

Conclusion
We have fractionalized the SAIVR epidemic model of order α using the well-known nonlocal character of the Caputo differential operator, which is ideal for studying the dynamics of disease transmission.The proposed Caputo SAIVR model is used for the study of the COVID-19 pandemic.Based on the Banach contraction principle, it is proven that the model has a unique solution.The stability of the Caputo model is established using Ulam-Hyers and its generalized form.One of the main contributions of this work is the leastsquares method used to acquire the model's fitted parameters from clinical samples of the virus (March-April, 2022); this method also optimizes the fractional order α (6.757e-01).
It has been shown that the Caputo model performs better than its classical counterpart.According to the results of a number of numerical simulations, the exposure rate must be reduced to control the pandemic effectively, and this is attainable when people practice social distancing and use protective masks.However, if such regulations are not carefully adhered to, a semi-developed country like Turkey will face serious difficulties.To effectively mitigate the spread of the pandemic, it is imperative to decrease the rate of exposure.This can be achieved by implementing social distancing measures and utilizing protective masks by individuals.Memory traces and hereditary traits are used to show the vanishing behavior of the model when α → 1.In future work, we plan to look into how non-singular differential operators change the classic SAIVR framework.

Figure 1 A
Figure 1 A flow chart for the dynamics of the SAIVR model given in (1) Unknown parameters are one of the most challenging aspects of infectious disease mathematics study.In epidemiology, deterministic models may employ several approaches to fit values for unknown parameters, including principle component iterated sorting and nonlinear least-squares curve fitting.In the COVID-19 SAIVR model, scientists need to know eleven epidemiological parameters: the inflow susceptible rate, asymptomatic/undetected, vaccinated rate, and recovery rate.Data from the WHO's official website was used to predict the outbreak of COVID-19 in Turkey (March-April 2022).The Caputo derivative has an initial recovery of R(0) = 50,269, people are initially asymptomatic/undetected A(0) = 23,954, people initially infected I(0) = 59,885, V (0) = 40,658 is the number of people who are vaccinated, and S(0) = 8,506,527 people who are susceptible to getting infected.Both the standard SAIVR model (1) and the proposed Caputo fractional SAIVR model (2) are used to obtain the COVID-19 pandemic's behavior.The Caputo fractional derivative operator clearly outperforms the standard approach as shown by the absolute mean errors: |E C | = 1.9861e-01 for the classical model and |E Cap | = 1.0474e-01 for the Caputo model.The biological parameters obtained with the technique of least-squares are stored in Table 1 wherein the optimized parameters, including the optimal value of α (the fractional order), are also shown.In addition, Figs. 2 and 3 compare the simulations obtained from the models and the real medical data.It is clear that the Caputo simulations agree better with the real medical data.It is also evident from the Box-Whisker chart in Fig. 4, wherein the Caputo version of the model describes summary statistics that match well enough with the real ones as shown in the first and third boxes in Fig. 4.

ηFigure 2 Figure 3 Figure 4
Figure 2 The behavior of the Caputo model in comparison with the statistical cases with the residual plot

Figure 5
Figure 5The dynamical behavior of the susceptible, asymptomatic, vaccinated, and removed individuals included in model (2) over a period[0, 100]

Figure 6 1 Figure 7 2 Figure 8
Figure 6The dynamical behavior of the infectious individuals included in model (2) over a period[0, 100]    under the effects of increasing exposure rate β 1

Figure 9
Figure 9The effects of memory trace on each population of model(2)

Table 1
The biological parameters used in the model