A case study of 2019-nCOV cases in Argentina with the real data based on daily cases from March 03, 2020 to March 29, 2021 using classical and fractional derivatives

In this study, our aim is to explore the dynamics of COVID-19 or 2019-nCOV in Argentina considering the parameter values based on the real data of this virus from March 03, 2020 to March 29, 2021 which is a data range of more than one complete year. We propose a Atangana–Baleanu type fractional-order model and simulate it by using predictor–corrector (P-C) method. First we introduce the biological nature of this virus in theoretical way and then formulate a mathematical model to define its dynamics. We use a well-known effective optimization scheme based on the renowned trust-region-reflective (TRR) method to perform the model calibration. We have plotted the real cases of COVID-19 and compared our integer-order model with the simulated data along with the calculation of basic reproductive number. Concerning fractional-order simulations, first we prove the existence and uniqueness of solution and then write the solution along with the stability of the given P-C method. A number of graphs at various fractional-order values are simulated to predict the future dynamics of the virus in Argentina which is the main contribution of this paper.


Introduction
The coronavirus disease 2019 (COVID-19) is considered as the most dangerous epidemic disease of this decade, which appeared for the first time in Wuhan (China) in the last month of 2019. It has been observed that coronavirus disease 2019 (COVID-19) is a vital health concern as it can be fatal specially in old-aged people. SARS-CoV-2 virus causes COVID-19 disease. Researchers and medical practitioners have much information about the death due to the clinical disease but they do not know much about its pathobiology. It has been seen that the characteristics of cellular responses to COVID-19 are not clearly known and understood, but based on the previous studies on SARS-CoV, a predictable sequence of events can be hypothesized. According to the cells that SARS-CoV infects, COVID-19 can be separated into three periods that coincide with the clinical parts of the disease [1].
The structure of SARS-CoV-2 needs to be understood. It has been found that SARS-CoV-2 is part of the family of beta-coronavirus, some of which caused two other epidemics named as MERS-CoV and SARS-CoV as can be found in Ref. [2]. In [3] Hui et al. analyzed the structure of such coronaviruses. According to the World Health Organisation (WHO), the rate of death due to SARS-CoV-2 is lower than the one due to SARS-CoV though the rate of transmission is higher in the case of SARS-CoV-2 than SARS-CoV. The RNA structure of different coronaviruses has been stated in [4]. It has been observed that SARS-CoV-2 has distinctive spike proteins and presents a specific peptide, namely PrrA, which allows one to divide the spike protein using cellular protease enzymes in order to disseminate the virus from the cell of the host easily [5,6]. So, the rate of transmission of SARS CoV-2 is much higher and this is the prime cause of its spreading throughout the globe in a very short period of time.
It has been observed that as soon as SARS-CoV-2 makes its entrance into the epithelial cell of alveoli of the respiratory tract of human being, the virus activates the immune response of human being due to the fast multiplication rate of SARS-CoV-2 cells. After that, the pulmonary tissue from the respiratory tract of human being is damaged causing the stimulation of more and more white blood cells. This process is known as the cytokine release syndrome [7], which can cause a multiple organ failure [8].
There are a number of stages included in the infection of COVID-19. The initial stage is asymptomatic and covers 1 to 2 days after getting infected. It has been noticed that after inhaling the virus it attaches to epithelial cells from the nasal cavity of human being and the virus starts the process of replication [9]. In this initial stage, infected people are able to spread the infection. During the second stage the virus is able to spread and move quickly down to the respiratory tract along with the conducting passage of air [10]. It has been observed that for 80% of the infected individuals the COVID-19 illness will be lenient and it is restricted to upper as well as conducting passages of air. The third stage of the disease involves hypoxia and ground glass infiltrates along with the growth of acute respiratory distress syndrome (ARDS). It has been seen that 20% of the infected individuals will move to third stage of the disease and they will develop pulmonary invasions which will lead to severe disease. From the initial estimation the rate of severity is about 2% but it varies with age [1]. In this stage SARS-Cov-2 virus is able to reach the units of the lung which are responsible for exchange of gas and infect the type II cells of alveoli. SARS-CoV is able to spread within type II cells and huge viral particles are liberated [11]. Old-aged people are at high risk since they have a low immune response, which can lead to severe consequences [12].
Mathematical models are the best way to describe the dynamics of the disease. Many mathematical models have been provided to illustrate the behavior of COVID-19 by using integer-and fractional-order derivatives [13][14][15][16][17][18]. A fractional-order COVID-19 model with delay has been proposed in [19]. Atangana et al. in [20] have discussed a mathematical model of COVID-19 with deterministic and stochastic approaches. The authors in [21][22][23] have also simulated the dynamics of 2019-nCOV by using effective mathematical models. A study of a fractional-order model of HIV is given in [24]. The au-thors in [25] proposed a research on the dynamics of a population model. In [26], a COVID-19 model along with necessary awareness programs has been explored. The roles of isolation and quarantine measures for COVID-19 outbreaks are given in [27]. The authors in [28] have solved a model of huanglongbing transmission with an effective numerical algorithm. A psychological model in fractional-order sense has been introduced in [29]. Kumar et al. in [30] have described a time-delay fractional-order mathematical model of oncolytic virotherapy. Many other epidemics have been described by using fractional mathematical models [31,32]. Recently, the authors in [33] studied the structure of rabies, and canine distemper virus epidemics by using generalised Caputo nonclassical derivative. Two fractional-order mathematical models for describing mosaic disease have been given in [34]. A environmental study in the fractional-order sense is given in [35]. This paper is formulated to trace the dynamics of COVID-19 cases in Argentina by using Atangana-Baleanu type fractional mathematical model. The study is organized as follows: In Sect. 2, we describe the integer-order model dynamics. In Sect. 3, we plot the real-data cases of COVID-19 in Argentina and calculate the parameter values for our simulations.

Model description
In this paper, we consider a compartmental mathematical model [36]. Natural births and death rates have not been considered in the model as those have no impact on the shortterm outbreaks of COVID-19. The mentioned model focuses on two distinct groups of susceptible individuals: susceptible individuals given by S(t) and confined individuals who follow lockdown or confinement intervention partially as the confinement is not perfect.
We symbolize the class of confined individuals as C(t). The rest of the population is compartmentalized as follows: exposed individuals E(t) at time t, asymptomatic (having no clinical symptoms or mild symptoms) individuals A(t), quarantined symptomatic infectious individuals Q(t), hospitalized or isolated individuals H(t) and recovered individuals In the given model, p is the transmission rate from confined susceptible humans to rejoined unconfined susceptible humans, qγ is the exposure rate of the asymptomatic humans, r 4 σ 2 is the progression rate from severely infected to hospitalization and quarantine rate is given by (1q)γ . r 1 σ 1 denotes the transmission rate of the asymptomatic humans to become severely infectious. r 2 σ 1 shows the transmission rate of confinement or hospitalization to isolation of the asymptomatic group. The natural recovery of soft symptomatic humans is represented by the rate (1r 1r 2 )σ 1 , whereas (1r 3r 4 )σ 2 is the recovery rate for Q classes. In the model, δ h , δ q , and δ a are stood for the COVID-19 death rates in the given model classes, respectively. The structure of the non-linear model based on the given parameters explanation is provided as follows: The parameters are given in Table 1 and the model structure is provided in Fig. 1. A brief discussion on the model characteristics like boundedness and positivity of solutions is given in reference [36]. The authors in [36] calculated the greatest disease-free equilibrium point, which is defined by Rate of transition from hospitalized to recovered group r 1 Fraction of A(t) becoming quarantined humans r 2 Rate of unquarantined infected humans going to be hospitalized r 3 Rate of quarantined infected humans moving to unquarantined infected humans r 4 Rate of quarantined infected humans going to be hospitalized Using the procedure given in [37], the matrices F and V for calculating the basic reproductive number are given by As a consequence, they obtained the basic reproductive number as the spectral-radius of the generation matrix, FV -1 : Here ρ(·) denotes the spectral-radius operator.
Proof First, we observe that the last two equations of (1) are not coupled to the rest of equations of the model. As the total population size, N 0 , is constant, we get S + C = N 0 -(E + A + Q + H + R). Then the local stability of the given system (1) can be analyzed via remaining model of variables (E, A, Q). As a consequence, we establish that the Jacobian matrix related to these variables is written by The roots of the next characteristic polynomial correspond to the eigenvalues of J : where Note that a 2 is always positive. a 0 and a 1 are positive as long as R c < 1. Therefore, all eigenvalues of J have negative real parts. Consequently, the epidemic-free equilibrium x 0 , is locally asymptotically stable if R c < 1.
Excluding confinement measures, i.e. if = 0, R c is convergent to the basic reproductive number, R 0 , written by Using (3), it follows that

Model calibration and forecasting
A recently proposed optimization scheme, which is an evolution of Levenberg-Marquardt one [36,38] and depends on trust-region-reflective (TRR) scheme, is employed to perform the model (1) calibration. This robust optimization method can be applied for simulating non-linear least-squares problems. The implementation of the scheme is done with the help of lsqcurvefit function, which is included in Optimization Toolbox from MATLAB.
Possible parameter values are calculated by this function. The real data of daily collected cases in Argentina are collected from the trusted data website which can be verified in [39]. We used 7-day running average of the given reported COVID-19 cases analyzed in the calibration of model because of changeable structure of real data as can be seen in where ψ = {β, σ 1 , σ 2 , σ 3 , r 1 , r 2 , r 3 , r 4 , δ a , δ q , δ h , η, p, c, , γ , q} is shown in Table 2.
In Figs. 2 and 3 it is shown that our model fits really well to the given real data of Argentina. The estimated value of the basic reproductive number R 0 is about ∼1.41 (95% CI: 1.2-1.6) as of March 10, 2021. This value could change (increase or decrease) in future cause of the new wave of Covid-19 which will totally depend on our health care measures. The disease could fade out at the end of December 2021. All necessary parameter values aspects applied to justify this scenario are mentioned in Table 2. In Fig. 4 the new daily and reported cases projected in Argentina from March 2020 to late May 2021 are given. In Fig. 5 cumulative infected cases and projected in the same dates are provided.

Lemma 2 ([40])
The solution of the given system for 0 < Φ < 1 is stated as

Unique solution existence for fractional-order model
Let consider the Banach space

R(t),
Applying (12) and (13), the two operators F, G are defined as We now state some basic axioms and the Lipschitzian hypothesis for showing the existence and uniqueness of a solution: (H1) There are C , D > 0 such that (H2) There is L > 0 such that ∀A,Ā ∈ Z and it follows that Proof The proof is given in two parts.
Step I: ConsiderĀ ∈ B, where B = {A ∈ Z : A ≤ ρ, ρ > 0} is closed and convex. The operator F provided in (14) gives Thus, F is a contraction.
Step-II: We want G to be relatively compact. Clearly, it is sufficient if G is equicontinuous and bounded. Indeed, G is continuous as is continuous and for all A ∈ B, one has Hence (16) shows the boundedness of G. For equicontinuity we assume t 1 > t 2 ∈ [0, τ ], so that The right-hand side in (17) goes to zero as t 1 → t 2 and then G(A(t 1 ) -G(A(t 1 ) → 0, as t 1 → t 2 since G is continuous. Because of the boundedness and continuity of G, then G is uniformly continuous and bounded. According to the Arzelá-Ascoli theorem, G is thus relatively compact and therefore entirely continuous. Consequently, the integral equation (12) and also the method have at least one solution.
We now proceed to show the uniqueness of the solution.

Theorem 3 Assuming (H 2 ), Eq. (12) has a unique solution and this implies that model (9) has a unique solution only if
Proof Consider T : Z → Z defined as Given A,Ā ∈ Z, then one can take Thus, T is a contraction. Therefore, Eq. (12) possesses a unique solution and so does model (9).

Derivation of solution
To date, number of numerical approximation algorithms have been proposed to solve the various kind of models describing the real world problems. When we apply such algorithms the analysis of convergence and stability are two important aspects of the effectiveness of the method. For solving our model, we are going to use predictor-corrector method, which has been defined probably for all fractional-order operators. Here we implement this algorithm in the sense of Atangana-Baleanu operator for writing the solution of the proposed COVID-19 model. For more information about this technique see Ref. [43]. We first recall Eq. (11) and consider The fractional Volterra integral equation is stated as According to the method proposed in [43] for Φ ∈ [0, 1], 0 ≤ t ≤ T and considering h = T/N and t n = nh, for n = 0, 1, 2, . . . , N ∈ Z + , the corrector formula for the given system is where and The predictor formula is attained as where Hence the corrector formulas for the given model (9) are stated as

Graphical justifications
Now we plot the fractional-order solution by considering the parameter values simulated in Table 2. This time we used Mathematica software. In Fig. 6, we plotted the susceptible population S(t) and in Fig. 7, the plot of confined individuals C(t) is given. We can observe that when time increases then the population of susceptible individuals is decreasing which is a good indication for Argentina. In Fig. 8, we plotted the nature of exposed individuals E(t). We notice that when the fractional-order values decreases then the peak of E(t) shifts to the next month. At fractional order Φ = 1, a peak occurs at near around 1 July 2020 and when we shifted to lower values of fractional order then we receive the peak near around 31 July 2020, at Φ = 0.95, near around 29 Sep 2020, at Φ = 0.85, and near around   28 Dec 2020 when Φ = 0.75. This shows how the fractional-order derivatives give us more varieties to predict the real-data structure for the future. Similarly, in Fig. 9, we plotted the nature of asymptomatic individuals A(t) where we can see that at fractional-order Φ = 1, peak occurs at near around 31 July 2020 and when we shift to lower values of fractional order then we receive the peak near around 15 Aug 2020 at Φ = 0.95, near around 15 Oct 2020 at Φ = 0.85, and near around 27 Jan 2021 when Φ = 0.75. In this series, Fig. 10 plotted the nature of quarantined individuals Q(t) where different peaks can be seen. Here at fractional order Φ = 1, peak occurs at near around 15 July 2020 and when we shifted to lower values of fractional order then we obtain the peak around 31 July 2020 at Φ = 0.95, around 29 Sep 2020 at Φ = 0.85, and around 28 Dec 2020 at Φ = 0.75. Figure 11 described the behavior of hospitalized individuals H(t) for the given time range. Here at fractional order Φ = 1, a peak occurs around 15 July 2020 and when we shifted to lower values of fractional order then we receive the peak around 31 July 2020 at Φ = 0.95, around 29 Sep

Conclusion
In our observations, we have explored the dynamics of the most deathly disease of this decade, COVID-19 in Argentina by considering parameter values based on real data. Firstly, we considered the reported cases of this virus from March 03, 2020 to March 29, 2021 and projected it for the future time period by using our models. We proposed a Atangana-Baleanu type fractional-order model and solved it employing predictorcorrector (P-C) algorithm. After analyzing the biological nature of this virus, we formulated a mathematical structure to define its dynamics. We used a well-known effective optimization method based on the renowned trust-region-reflective (TRR) scheme to perform the model calibration. We plotted the real cases of COVID-19 and fitted our integerorder model with the data along with the calculation of basic reproductive number. In the simulations of fractional-order analysis, first we proved the existence of unique solution and then wrote the solution of the model along with the stability of the given P-C method. We performed separate graphs for every model class at different fractional-order values to predict the future dynamics of the virus in Argentina. Finally, we conclude that this virus is under control in Argentina for future under the given data range and all health care measures. In the future, the given estimated parameter values and models can be used to provide further predictions on the transmission of this virus. The proposed mathematical model is effective and trustable for future use. Not only from a mathematical point of view but also from a medical perspective, this research study may become very useful for the medical authorities to simulate the outbreaks of COVID-19 in Argentina for a future time period. Still, all the peak predictions which have been simulated in this paper, do not claim with 100% certainty that the dynamics of COVID-19 will be as it is. It may differ because of many health care measures (like isolation strategies, hospital facilities, and vaccine availability) which have been mentioned in this study.