The world is still addressing the Coronavirus illness 2019 (COVID-19), which is caused by the new Coronavirus SARSCoV-2, a highly aggressive virus that attacks the individual respiratory system. The hospitalized individuals' ailments were linked to the marine and moist animal industries in Wuhan, Hubei Province, China. COVID-19 spreads from person to person by touching contaminated surfaces and inhalation of infected persons' respiratory droplets. Those who have been infected with COVID-19 have reported high fever, persistent cough, and exhaustion. Nonetheless, depending on the immune system, COVID-19 symptoms and consequences differ from person to person. People with a strong immune response seem to be more likely to get mild:to:moderate illnesses as well as recover avoid going to the hospital. Various investigations, however, have identified other symptoms such as neurological illnesses and gastroenteritis of different severity. With so many waves of infection, the illness caused numerous deaths in many countries. COVID-19 outbreaks have occurred in Italy, with the population suffering the effects of the consequences. The number of confirmed incidence and mortality in every phase has been published, and there appears to be an increasing incidence. On February 21, 2020, the first Italian victim of COVID-19, a 38-year-old male hospitalized at Codogno Hospital in Lodi, was diagnosed. On the 12th of February, 2022, it has infected over 424,636,034 people over the world, resulting in 5903,485 deaths and 349,857,774 recoveries. According to reports, the mortality rate in waves 1 and 2 was 1%. Many social programmers and events have been discontinued or extended as a result of the epidemic. The T-20 cricket world cup will be hosted in Australia in 2020, while the Summer Olympics, which were scheduled to be held in Tokyo, have been postponed. The Indian Premier League, one of the most popular cricket events, has been relocated from India to the United Arab Emirates.
Its importance has been demonstrated by the construction of mathematical models in the fields of epidemiology and physics. The Coronavirus infection has been examined by several researchers from various perspectives. While biologists and mathematicians working on the systems of the COVID-19 disease analyzed and constructed mathematical systems based on real-world cases from various countries, and offered information on the infection's peak and clearance. In this context, are some mathematical models that have been developed for this disease. The information from Italy is taken into account, and a mathematical model for the COVID-19 disease is developed, with its study reported in. Examines the number of genuine instances from the Mexican population using a mathematical model. Proposes a fractional SEIR model utilizing the wavelet approach. The authors investigated the influence of social distance and other factors that might be regarded important for the reduction of COVID-19 infection in. The authors used a mathematical modeling technique to evaluate genuine infected patients from Saudi Arabia and generated results on disease eradication in the nation, Describes a comparative study of Coronavirus infection dynamics. In suggests some additional relevant work on COVID-19 modeling and associated illness outcomes. In, Paul et al. analyzed the scenario analysis of COVID-19 pandemic using SEIR epidemic model. A deeper understanding of the pandemic dynamics, including the characteristics of Covid-19 transmission, was made possible by the modeling technique. We have previously published research on fractional order phenomena. The novel fractional operator has shown to be quite effective in solving a variety of mathematical modeling problems as well as some recent work on COVID-19. From the study of data obtained from Wuhan, Li et al. calculated the epidemiology and discovered the mean incubation time was 5.2 days. The authors of highlighted some interesting outcomes of Corona virus disease. In introduced deeper investigation of modified epidemiological computer virus model containing the Caputo operator. Furthermore, the researchers in analyzed fractional derivatives of the COVID-19 infection models and presented some recommendations for infection minimization in the form of lockdown and control measures.
Fractional order modeling is a useful tool that has been used to explore the nature of diseases since the fractional derivative is an extension of the integer-order derivative. In order to replicate real-world issues, several innovative fractional operators with various properties have been designed. In addition, the integer derivative has a local identity, whereas the fractional derivative has a global character. Numerous varieties of fractional derivatives, both with and without singular kernels, are available today. Leibniz's query from 1695 marks the beginning of the fractional derivative. The fractional derivative also improves in the improvement of the system's consistency domain. We have the derivatives of Caputo, Riemann-Liouville, and Katugampola for singular kernels. There are two varieties of fractional derivatives without singular kernels: the Caputo-Fabrizio fractional derivative, which has an exponential kernel, and the Atangana-Baleanu fractional derivative, which has a Mittag-Leffler kernel. While memory and genetic properties are involved, working with fractional-order derivatives is crucial because it provides a more accurate technique to describe COVID-19 outbreaks. Numerous academic articles, monographs, and novels have provided evidence to support this claim; for instance,. Motivated by the current research, we present and analyze the SEIQRD model in Caputo sense. The Caputo derivative is particularly useful for discussing real-world situations since it permits traditional beginning and boundary conditions to be used in the derivation, and the derivative of a constant is zero, whereas the Riemann-Liouville fractional derivative does not. It is quite challenging to genuinely create an appropriate mathematical model using classical differentiation in the situation of COVID-19 because to the large number of uncertainties, unknowns, and disinformation. Generally, non-local operators are better suited for such circumstances because, depending on whether power law, fading memory, or overlap effects are taken into account, they can represent non-localities and certain memory effects.
We present the reader with some important definitions and characteristics of fractional derivatives in Section 2. In Section 3, we have established the SEIQRD epidemic model of Covid-19 in Caputo sense. We have investigated the existence, uniqueness, non-negative, boundedness criterion and stability analysis of the solution of model in Section 4. In Section 5, the fractional-order Taylor's approach in Caputo derivative is utilized to approximate the solution to the proposed model. The numerical study is given using MATLAB (2018a) in Section 6. Finally, the paper's conclusion is found in Section 7.
We provide the reader with some useful definitions and characteristics of fractional derivatives.
Definition 1 "The Caputo fractional derivative of order 0<phi<=1 for the function u: Cn[0,  ] R is defined as  where Cn [0,  ] is a n tines continuously differentiable function and the Gamma function is defined by Gamma() such that n-1<phi<n".
Theorem 1 "If  is piecewise continuous, then , where the Laplace transform is denoted by L(g(t))".
Theorem 2 "One-parametric and two-parametric Mittag-Leffler functions are described as follows:  and , where a1, a2 R+".
Lemma 1 "Let 0<phi<=1, u(t) C[p, q] and if  is continuous in [p, q], then ,
where 0<=z<=x,  x (p, q]".
Note 1 "If  then u(t) is a non-decreasing (non-increasing) function for t [p, q]".
Lemma 2 "Let us consider the fractional order system as  with  and . For calculate the equilibrium points, we have Psi(Y) = 0. These equilibrium points are locally asymptotically stable iff each eigen value lambdaj of the Jacobian matrix  calculated at the equilibrium points satisfies ".
Lemma 3 "Assume that u(t) R+ is a differentiable function. Then, for any t>0, 
The mathematical model of COVID-19 transmission formulated in this study was motivated by the study of. In the present study, the model will be divided into six compartments [see Fig 1]. The total human population to be considered is denoted as N(t), and at any time, it comprises of the susceptible (S), exposed (E), infected (I), quarantined (Q), recovered (R), and death (D) compartments, respectively.
Now we formulate the SEIQRD model with fractional order derivatives with Caputo operator of order 0<phi<=1.
Now for the sake of convenience of calculation, we redefine the parameters [see Table 1] as   Thus, the modified model system (3.2) can be finally written in the following form with initial conditions,      
The initial conditions are 
The following are the necessary and sufficient conditions for a fractional order system's solution to exist and be unique:
Theorem 4.1.1. For each initial condition, there exists a unique solution of fractional order system (3.3).
Proof We are looking for a sufficient condition for the presence and uniqueness of system (3.3) solutions in the region Pix(0, T] where
. The method employed in is used. Consider a mapping  where Y = (S, E, I, Q, R, D) and :      
For any      
Where G = max{G1, G2, G3, G4, G5} and 
As a result, F(Y) fulfils the Lipschitz requirement. As a consequence, fractional order system (3.3) exists and is unique.
Proposition The region  is non-negative invariant for the model (3.3)   t>=0.
Proof We have   
Using Laplace transform and Theorem 7.2 in, we have
 where z is the Laplace transform parameter.
Appling inverse Laplace transform, we have 
According to Mittag-Leffler function, 
Hence, .
As a result, the functions S, E, I, Q, R, and D are all non-negative.
The system's equilibrium may be found by solving the model (3.3) i.e., 
The model (3.3) has two equilibrium points namely, the infection free equilibrium  and the epidemic equilibrium point E1 = (S*, E*, I*, Q*, R*, D*), where , , , 
The next-generation matrix technique is used to calculate the model's basic reproduction number RCovid 19, which may be obtained from the maximum eigen value of the matrix  where,  
The Jacobian matrix of the model (3.3) at E0 is given by  where  with A11 = -a, A22 = -(e+c), A = e,  A43 = f, A53 = g, A63 = h, A44 = -(k+l+c), A54 = k, A64 = l, A55 = -c.
Theorem 4.5.1. When RCovid 19< 1, the system (3.3) is globally asymptotically stable, and unstable when RCovid 19>1 at E0.
Proof Using the appropriate Lyapunov function 
The aforementioned function's time derivative is 
From (3.3) we get, 
Now,  
Since , it follows that  
Hence if RCovid 19<1, then .
As a result of LaSalle's use of Lyapunov's concept, the point E0 is globally asymptotically stable and unstable if RCovid 19>1.
Theorem 4.6.1. If RCovid 19>1, the system (3.3) is globally asymptotically stable at E1.
Proof The Lyapunov function of the Goh-Volterra form's is as follows: 
Using Lemma 3 and taking Caputo derivative, we get 
Using (3.3) we get, 
Eq (3.3) gives us the steady state, 
Substituting Eq (4.8) into (4.7) we have 
Further simplification gives, 
Taking all infected classes that do not have a single star (*) from (4.9) and equal to zero: 
The steady state was slightly perturbed between (3.3) and (4.10), resulting in: 
Using (4.11) into (4.9) gives: 
Using A. M>=G. M., we have , .
Thus, .
The point E1 is globally asymptotically stable if RCovid 19>1.
As discussed in Theorem 4.1.1, the solution of the system (3.3) is unique. To obtain the numerical solution of the system (3.3), Taylor's theorem will be used.
As a result, we proceed with the model's 1st equation as follows: 
Consider the set of points [0, A] as the points on which we are prepared to approximate the system's solution. Actually, we are unable to calculate S(t), which will be the system's necessary solution. We divide [0, A], into P subintervals [tr, tr+1] of length, i.e., , by using the nodes tr = rm, for r = 0, 1, 2,...,P. We extend the Taylor's theorem at about t = t0, we have a constant k [0,A], such that 
Now substitute , and t = t1 in (5.2), which provides 
If m is small, we ignore the higher terms, then (5.3), implies 
A general formula of expanding about tr = tr+m, is 
In similar way, we get     
Numerical simulations employing Taylor's theorem are carried out with the help of MATLAB software to support the mathematical study of the system (3.3). This section is divided into four parts. The stability of our proposed model is discussed at E0 and E1 in Part 1. Part 2 delves into the dynamical behavior of all individuals of various fractional orders. Part 3 is to explore the varying effects of face masks. Part 4 is to determine whether the model (3.3) fits the data. One of the key components in the verification of an epidemiological model is the fitting of the parameters. We ran numerical simulations to contrast the output of our model with actual data from a number of reports released by the WHO and worldometer. Italy has a population of about 60,278,248 people. In Italy, there are 7.2 births per 1000 people. The computed recruiting rate is  per day.
The stability of our suggested model is discussed in this section. The parameter values used for the numerical simulations in Part 1 is provided in Table 2. Fig 2(A)-2(E) depict the nature of all cases corresponding to phi = 0.98. From the following figures, we have observed that the system is locally asymptotically stable at E0.
To analyze the dynamical behavior of all people, the values of the parameters in Table 2 are employed. Fig 3(A)-3(F) depict all individuals' behavior over time for various fractional orders phi. Fig 3(A) depicts that the number of susceptible individuals increases when phi changes from 0.8 to 0.95. An increase value of phi leads to decrease in the exposed rate in the exposed population in Fig 3(B). We see in Fig 3(C) that number of infected individuals increases when phi changes from 0.8 to 0.95. Fig 3(D) depicts that the number of quarantined individuals increases with time when phi decreases. The number of recovered individuals increases when phi changes from 0.8 to 0.95 in Fig 3(E). Fig 3(F) depicts that the number of death individuals increase with time when phi increases.
Part 3 of the numerical simulation investigates at how changing value impacts the fundamental reproduction number calculated in this work. Table 2 shows the parameter values utilized in the numerical simulations for Part 3. The acquired findings are shown in Table 3 after computing the fundamental reproduction numbers and utilizing the model parameters from Table 2. Table 3 shows that if a higher number of individuals in a community constantly utilize face masks, the COVID-19 epidemic can be decreased.
Fig 4 depicts that the values of RCovid 19 decrease when alpha increase. The various consequences of wearing face masks were also investigated in this study, and it was discovered that wearing face masks on a consistent and suitable basis can inhibit the spreading of the COVID-19 pandemic.
The impact of alpha and phi on the Infected individuals (I(t)) is depicted in Fig 5(A) and 5(b). Based on the following figures, it can be noted that the implementation of maximum portion of population who use a face masks in order to effectively reduce COVID-19 transmission.
This section describes the data matching and model validation of the system (3.3) for Infected instances. Table 2 depicts the parametric values.
Fig 6 depicts the graphical representation of the infected cases respectively of the model (3.3) and the real infected cases in Italy [see Table 4] from 1st January 2022 to 31st January 2022.
The present study's possible goal is to develop a mathematical model for studying COVID-19 transmission patterns using actual pandemic cases in Italy, assisted by epidemiological modeling. The fractional order SEIQRD model was constructed and explored in this article in order to better explain the dynamics of the COVID-19 epidemic in Italy. We employed nonlinear analysis to demonstrate the model's existence and uniqueness. The model's fundamental reproduction number was also calculated using the next generation matrix technique. In order to stop the virus from spreading throughout the nation, our main goal is to establish the fundamental reproductive number and equilibrium. Furthermore, the global stability at the points E0 and E1 has been demonstrated. The results reveal that if RCovid 19<1, the point E0 is globally asymptotically stable. Also if RCovid 19>1, the point E1 is global asymptotic stable. Furthermore, using the fractional Taylor's approach, numerical analysis was done to establish an approximate solution for the suggested model. From the 1st of January 2022 to the 31st of January 2022, we have compared model values with real-world scenarios in Italy. Real data was also used to fit the model in order to forecast infected population instances in real life. In real-world dynamical processes, such as epidemic propagation, fractional calculus plays a vital role. The strength of memory effects, which is regulated by the order of fractional derivatives, is discovered to be dependent on the system dynamics. If the order of derivatives for the same set of parametric values is changed, the results will be changed (Fig 3). This study looked at the many consequences of wearing face masks, and it was discovered that wearing face masks on a consistent and suitable basis can help reduce the propagate of the COVID-19 disease (Figs 4 and 5). Currently, research on a vaccine to avert the COVID-19 pandemic is showing promising results, with Pfizer claiming that their vaccine has a 95% effectiveness rate. However, it will be some time before the vaccinations are widely distributed around the world. As a result, wearing a face mask should be made mandatory until everyone has access to vaccinations. Ministries and public health professionals may be able to develop strategic strategies to close vaccination gaps and stop outbreaks in the future with the use of the research findings from the current study. Future studies should use the methodology provided in this work to the third wave of infected patients in Italy to assess the efficacy of existing COVID-19 prevention strategies.