Fractional epidemic model of coronavirus disease with vaccination and crowding effects

Suhail Saleem,Muhammad Rafiq,Nauman Ahmed,Muhammad Shoaib Arif,Ali Raza,Zafar Iqbal,Shafiullah Niazai,Ilyas Khan
DOI: https://doi.org/10.1038/s41598-024-58192-7
IF: 4.6
2024-04-10
Scientific Reports
Abstract:Most of the countries in the world are affected by the coronavirus epidemic that put people in danger, with many infected cases and deaths. The crowding factor plays a significant role in the transmission of coronavirus disease. On the other hand, the vaccines of the covid-19 played a decisive role in the control of coronavirus infection. In this paper, a fractional order epidemic model (SIVR) of coronavirus disease is proposed by considering the effects of crowding and vaccination because the transmission of this infection is highly influenced by these two factors. The nonlinear incidence rate with the inclusion of these effects is a better approach to understand and analyse the dynamics of the model. The positivity and boundedness of the fractional order model is ensured by applying some standard results of Mittag Leffler function and Laplace transformation. The equilibrium points are described analytically. The existence and uniqueness of the non-integer order model is also confirmed by using results of the fixed-point theory. Stability analysis is carried out for the system at both the steady states by using Jacobian matrix theory, Routh–Hurwitz criterion and Volterra-type Lyapunov functions. Basic reproductive number is calculated by using next generation matrix. It is verified that disease-free equilibrium is locally asymptotically stable if and endemic equilibrium is locally asymptotically stable if . Moreover, the disease-free equilibrium is globally asymptotically stable if and endemic equilibrium is globally asymptotically stable if . The non-standard finite difference (NSFD) scheme is developed to approximate the solutions of the system. The simulated graphs are presented to show the key features of the NSFD approach. It is proved that non-standard finite difference approach preserves the positivity and boundedness properties of model. The simulated graphs show that the implementation of control strategies reduced the infected population and increase the recovered population. The impact of fractional order parameter is described by the graphical templates. The future trends of the virus transmission are predicted under some control measures. The current work will be a value addition in the literature. The article is closed by some useful concluding remarks.
multidisciplinary sciences
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to evaluate the impact of vaccination and crowd - gathering effects on the transmission dynamics of the COVID - 19 virus, especially within the framework of fractional - order calculus. Specifically, the author studies the following aspects by constructing a fractional - order SIVR (Susceptible - Infected - Vaccinated - Recovered) epidemic model: 1. **Vaccination and crowd - gathering effects**: Explore how these two factors affect the spread of the COVID - 19 virus. 2. **Non - linear incidence rate**: Introduce a non - linear incidence rate to more accurately describe the dynamic behavior of virus transmission. 3. **Mathematical properties of the model**: Ensure the positivity and boundedness of the model and verify its existence and uniqueness. 4. **Equilibrium point analysis**: Determine the disease - free equilibrium point and the endemic equilibrium point of the system and conduct stability analysis. 5. **Basic reproduction number \(R_0\)**: Calculate and analyze the value of \(R_0\) to judge the local and global asymptotic stability of the disease - free equilibrium point and the endemic equilibrium point. 6. **Numerical simulation**: Use the non - standard finite - difference method (NSFD) to conduct numerical simulation on the model to verify the effectiveness of control strategies. ### Model background Since the COVID - 19 virus was first discovered in Wuhan, China at the end of 2019, it has spread rapidly around the world, leading to a serious public health crisis. Vaccination has played an important role in controlling the epidemic, but crowd - gathering is still an important factor in virus transmission. In order to more accurately describe the impact of these factors on the epidemic, the author chooses to use fractional - order calculus to construct the SIVR model. ### Main contributions 1. **Application of fractional - order calculus**: By introducing fractional - order derivatives, the model can better capture the memory and hereditary characteristics of the system, thus more realistically reflecting the dynamic behavior of infectious diseases. 2. **Theoretical analysis**: Through strict mathematical analysis, the positivity, boundedness, existence and uniqueness of the model and the stability of the equilibrium points are proved. 3. **Numerical simulation**: Through numerical simulation, the effects of different control measures are shown, especially the role of vaccination and reducing crowd - gathering in reducing the number of infected people and increasing the number of recovered people. ### Conclusion This research provides valuable references for governments and public health departments to help them formulate more effective epidemic prevention and control strategies. At the same time, this research also provides new ideas for future mathematical modeling and the application of fractional - order calculus in epidemiology. ### Formula summary 1. **Definition of fractional - order derivative**: \[ {}_0^C D_t^\alpha f(t)=\frac{1}{\Gamma(n - \alpha)}\int_a^t(t - x)^{n-\alpha - 1}f^{(n)}(x)dx,\quad n - 1<\alpha <n\in\mathbb{N} \] 2. **Laplace transform**: \[ \mathcal{L}\{{}_0^C D_t^\alpha f(t)\}=s^\alpha F(s)-\sum_{k = 0}^{n - 1}s^{k-\alpha}f^{(k)}(0),\quad n - 1<\alpha <n\in\mathbb{N} \] 3. **Mittag - Leffler function**: \[ E_{\alpha,\beta}(z)=\sum_{k = 0}^\infty\frac{z^k}{\Gamma(\alpha k+\beta)},\quad\alpha,\beta > 0,\alpha,\beta\in\mathbb{R},z\in\mathbb{C} \] 4. **Basic reproduction number \(R_0\)**: \[ R_0=\frac{\beta\Lambda}{(\delta_2+\mu)(\gamma+\delta_1+\mu)} \] Through these formulas and theoretical analysis, the author comprehensively explores the dynamic characteristics of the COVID - 19 virus transmission and provides for the future.