COVID-19 pandemic in India: a mathematical model study

Modeling the effects of cross immunity and control measures on competitive dynamics of SARS-CoV-2 variants in the USA, UK, and Brazil

Mutation in the SARS-CoV-2 virus may lead to the evolution of new variants. The dynamics of these variants varied among countries. Identification of the governing factors responsible for distinctions in their dynamics is important for preparedness against future severe variants. This study investigates the impact of cross immunity and control measures on the competition dynamics of the Alpha, Gamma, Delta, and Omicron variants. The following questions are addressed using an n-strain deterministic model: (i) Why do a few variants fail to cause a wave even after winning the competition? (ii) In what scenarios a new variant cannot replace the previous one? The model is fitted and cross-validated with the data of COVID-19 and its variants for the USA, UK, and Brazil. The model analysis highlights implementations of the following measures against any deadlier future variant: (i) an effective population-wide cross-immunity from less lethal strains and (ii) strain-specific vaccines targeting the novel variant. The system exhibits a fascinating dynamical behavior known as an endemic bubble due to Hopf bifurcation. It is observed that the actual situation in which Omicron won the competition from Delta followed by no wave due to Delta may turn into a competitive periodic coexistence of two strains due to substantial disparity in fading rates of cross-immunity. Global sensitivity analysis is conducted to quantify uncertainties of model parameters. It is found that examining the impact of cross-immunity is as crucial as vaccination.

Nonlocal models in biology and life sciences: Sources, developments, and applications

Mathematical modeling is one of the fundamental techniques for understanding biophysical mechanisms in developmental biology. It helps researchers to analyze complex physiological processes and connect like a bridge between theoretical and experimental observations. Various groups of mathematical models have been studied to analyze these processes, and the nonlocal models are one of them. Nonlocality is important in realistic mathematical models of physical and biological systems when local models fail to capture the essential dynamics and interactions that occur over a range of distances (e.g., cell-cell, cell-tissue adhesions, neural networks, the spread of diseases, intra-specific competition, nanobeams, etc.). This review illustrates different nonlocal mathematical models applied to biology and life sciences. The major focus has been given to sources, developments, and applications of such models. Among other things, a systematic discussion has been provided for the conditions of pattern formations in biological systems of population dynamics. Special attention has also been given to nonlocal interactions on networks, network coupling and integration, including brain dynamics models that provide an important tool to understand neurodegenerative diseases better. In addition, we have discussed nonlocal modeling approaches for cancer stem cells and tumor cells that are widely applied in the cell migration processes, growth, and avascular tumors in any organ. Furthermore, the discussed nonlocal continuum models can go sufficiently smaller scales, including nanotechnology, where classical local models often fail to capture the complexities of nanoscale interactions, applied to build biosensors to sense biomaterial and its concentration. Piezoelectric and other smart materials are among them, and these devices are becoming increasingly important in the digital and physical world that is intrinsically interconnected with biological systems. Additionally, we have reviewed a nonlocal theory of peridynamics, which deals with continuous and discrete media and applies to model the relationship between fracture and healing in cortical bone, tissue growth and shrinkage, and other areas increasingly important in biomedical and bioengineering applications. Finally, we provided a comprehensive summary of emerging trends and highlighted future directions in this rapidly expanding field.

A six-compartment model for COVID-19 with transmission dynamics and public health strategies

The global crisis of the COVID-19 pandemic has highlighted the need for mathematical models to inform public health strategies. The present study introduces a novel six-compartment epidemiological model that uniquely incorporates a higher isolation rate for unreported symptomatic cases of COVID-19 compared to reported cases, aiming to enhance prediction accuracy and address the challenge of initial underreporting. Additionally, we employ optimal control theory to assess the cost-effectiveness of interventions and adapt these strategies to specific epidemiological scenarios, such as varying transmission rates and the presence of asymptomatic carriers. By applying this model to COVID-19 data from India (30 January 2020 to 24 November 2020), chosen to capture the initial outbreak and subsequent waves, we calculate a basic reproduction number of 2.147, indicating the high transmissibility of the virus during this period in India. A sensitivity analysis reveals the critical impact of detection rates and isolation measures on disease progression, showing the robustness of our model in estimating the basic reproduction number. Through optimal control simulations, we demonstrate that increasing isolation rates for unreported cases and enhancing detection reduces the spread of COVID-19. Furthermore, our cost-effectiveness analysis establishes that a combined strategy of isolation and treatment is both more effective and economically viable. This research offers novel insights into the efficacy of non-pharmaceutical interventions, providing a tool for strategizing public health interventions and advancing our understanding of infectious disease dynamics.

A mathematical model for the transmission of co-infection with COVID-19 and kidney disease

The world suffers from the acute respiratory syndrome COVID-19 pandemic, which will be scary if other co-existing illnesses exacerbate it. The co-occurrence of the COVID-19 virus with kidney disease has not been available in the literature. So, further research needs to be conducted to reveal the transmission dynamics of COVID-19 and kidney disease. This study aims to create mathematical models to understand how COVID-19 interacts with kidney diseases in specific populations. Therefore, the initial step was to formulate a deterministic Susceptible-Infected-Recovered (SIR) mathematical model to depict the co-infection dynamics of COVID-19 and kidney disease. A mathematical model with seven compartments has been developed using nonlinear ordinary differential equations. This model incorporates the invariant region, disease-free and endemic equilibrium, along with the positivity solution. The basic reproduction number, calculated via the next-generation matrix, allows us to assess the stability of the equilibrium. Sensitivity analysis is also utilised to understand the influence of each parameter on disease spread or containment. The results show that a surge in COVID-19 infection rates and the existence of kidney disease significantly enhances the co-infection risks. Numerical simulations further clarify the potential outcomes of treating COVID-19 alone, kidney disease alone, and co-infected cases. The study of the potential model can be utilised to maximise the benefits of simulation to minimise the global health complexity of COVID-19 and kidney disease.

Within-host delay differential model for SARS-CoV-2 kinetics with saturated antiviral responses

The present study discussed a model to describe the SARS-CoV-2 viral kinetics in the presence of saturated antiviral responses. A discrete-time delay was introduced due to the time required for uninfected epithelial cells to activate a suitable antiviral response by generating immune cytokines and chemokines. We examined the system's stability at each equilibrium point. A threshold value was obtained for which the system switched from stability to instability via a Hopf bifurcation. The length of the time delay has been computed, for which the system has preserved its stability. Numerical results show that the system was stable for the faster antiviral responses of epithelial cells to the virus concentration, i.e., quick antiviral responses stabilized patients' bodies by neutralizing the virus. However, if the antiviral response of epithelial cells to the virus increased, the system became unstable, and the virus occupied the whole body, which caused patients' deaths.

COVID-19 transmission dynamics and the impact of vaccination: modelling, analysis and simulations

Despite the lifting of COVID-19 restrictions, the COVID-19 pandemic and its effects remain a global challenge including the sub-Saharan Africa (SSA) region. Knowledge of the COVID-19 dynamics and its potential trends amidst variations in COVID-19 vaccine coverage is therefore crucial for policy makers in the SSA region where vaccine uptake is generally lower than in high-income countries. Using a compartmental epidemiological model, this study aims to forecast the potential COVID-19 trends and determine how long a wave could be, taking into consideration the current vaccination rates. The model is calibrated using South African reported data for the first four waves of COVID-19, and the data for the fifth wave are used to test the validity of the model forecast. The model is qualitatively analysed by determining equilibria and their stability, calculating the basic reproduction number R 0 and investigating the local and global sensitivity analysis with respect to R 0 . The impact of vaccination and control interventions are investigated via a series of numerical simulations. Based on the fitted data and simulations, we observed that massive vaccination would only be beneficial (deaths averting) if a highly effective vaccine is used, particularly in combination with non-pharmaceutical interventions. Furthermore, our forecasts demonstrate that increased vaccination coverage in SSA increases population immunity leading to low daily infection numbers in potential future waves. Our findings could be helpful in guiding policy makers and governments in designing vaccination strategies and the implementation of other COVID-19 mitigation strategies.

Stability analysis and optimal control of a fractional-order generalized SEIR model for the COVID-19 pandemic

In view of the spread of corona virus disease 2019 (COVID-19), this paper proposes a fractional-order generalized SEIR model. The non-negativity of the solution of the model is discussed. Based on the established threshold R 0 , the existence of the disease-free equilibrium and endemic equilibrium is analyzed. Then, sufficient conditions are established to ensure the local asymptotic stability of the equilibria. The parameters of the model are identified based on the statistical data of COVID-19 cases. Furthermore, the validity of the model for describing the COVID-19 outbreak is verified. Meanwhile, the accuracy of the relevant theoretical results are also verified. Considering the relevant strategies of COVID-19 prevention and control, the fractional optimal control problem (FOCP) is proposed. Numerical schemes for Riemann-Liouville (R-L) fractional-order adjoint system with transversal conditions is presented. Based on the relevant statistical data, the corresponding FOCP is numerically solved, and the control effect of the COVID-19 outbreak under the optimal control strategy is discussed.

Transition dynamics between a novel coinfection model of fractional-order for COVID-19 and tuberculosis via a treatment mechanism

In this paper, a fractional-order coinfection model for the transmission dynamics of COVID-19 and tuberculosis is presented. The positivity and boundedness of the proposed coinfection model are derived. The equilibria and basic reproduction number of the COVID-19 sub-model, Tuberculosis sub-model, and COVID-19 and Tuberculosis coinfection model are derived. The local and global stability of both the COVID-19 and Tuberculosis sub-models are discussed. The equilibria of the coinfection model are locally asymptotically stable under certain conditions. Later, the impact of COVID-19 on TB and TB on COVID-19 is analyzed. Finally, the numerical simulation is carried out to assess the effect of various biological parameters in the transmission dynamics of COVID-19 and Tuberculosis coinfection.

Fractional order mathematical model for B.1.1.529 SARS-Cov-2 Omicron variant with quarantine and vaccination

In this paper, a fractional order nonlinear model for Omicron, known as B.1.1.529 SARS-Cov-2 variant, is proposed. The COVID-19 vaccine and quarantine are inserted to ensure the safety of host population in the model. The fundamentals of positivity and boundedness of the model solution are simulated. The reproduction number is estimated to determine whether or not the epidemic will spread further in Tamilnadu, India. Real Omicron variant pandemic data from Tamilnadu, India, are validated. The fractional-order generalization of the proposed model, along with real data-based numerical simulations, is the novelty of this study.

A new fractional mathematical model to study the impact of vaccination on COVID-19 outbreaks

This study proposes a new fractional mathematical model to study the impact of vaccination on COVID-19 outbreaks by categorizing infected people into non-vaccinated, first dose-vaccinated, and second dose-vaccinated groups and exploring the transmission dynamics of the disease outbreaks. We present a non-linear integer order mathematical model of COVID-19 dynamics and modify it by introducing Caputo fractional derivative operator. We start by proving the good state of the model and then calculating its reproduction number. The Caputo fractional-order model is discretized by applying a reliable numerical technique. The model is proven to be stable. The classical model is fitted to the corresponding cumulative number of daily reported cases during the vaccination regime in India between 01 August 2021 and 21 July 2022. We explore the sensitivities of the reproduction number with respect to the model parameters. It is shown that the effective transmission rate and the recovery rate of unvaccinated infected individuals are the most sensitive parameters that drive the transmission dynamics of the pandemic in the population. Numerical simulations are used to demonstrate the applicability of the proposed fractional mathematical model via the memory index at different values of 0 . 7 , 0 . 8 , 0 . 9 and 1. We discuss the epidemiological significance of the findings and provide perspectives on future health policy tendencies. For instance, efforts targeting a decrease in the transmission rate and an increase in the recovery rate of non-vaccinated infected individuals are required to ensure virus-free population. This can be achieved if the population strictly adhere to precautionary measures, and prompt and adequate treatment is provided for non-vaccinated infectious individuals. Also, given the ongoing community spread of COVID-19 in India and almost the pandemic-affected countries worldwide, the need to scale up the effort of mass vaccination policy cannot be overemphasized in order to reduce the number of unvaccinated infections with a view to halting the transmission dynamics of the disease in the population.

An SEQAIHR model to study COVID-19 transmission and optimal control strategies in Hong Kong, 2022

During the COVID-19 pandemic, one of the major concerns was a medical emergency in human society. Therefore it was necessary to control or restrict the disease spreading among populations in any fruitful way at that time. To frame out a proper policy for controlling COVID-19 spreading with limited medical facilities, here we propose an SEQAIHR model having saturated treatment. We check biological feasibility of model solutions and compute the basic reproduction number ( R 0 ). Moreover, the model exhibits transcritical, backward bifurcation and forward bifurcation with hysteresis with respect to different parameters under some restrictions. Further to validate the model, we fit it with real COVID-19 infected data of Hong Kong from 19th December, 2021 to 3rd April, 2022 and estimate model parameters. Applying sensitivity analysis, we find out the most sensitive parameters that have an effect on R 0 . We estimate R 0 using actual initial growth data of COVID-19 and calculate effective reproduction number for same period. Finally, an optimal control problem has been proposed considering effective vaccination and saturated treatment for hospitalized class to decrease density of the infected class and to minimize implemented cost.

Complex dynamics and control analysis of an epidemic model with non-monotone incidence and saturated treatment

In this manuscript, we consider an epidemic model having constant recruitment of susceptible individuals with non-monotone disease transmission rate and saturated-type treatment rate. Two types of disease control strategies are taken here, namely vaccination for susceptible individuals and treatment for infected individuals to minimize the impact of the disease. We study local as well as global stability analysis of the disease-free equilibrium point and also endemic equilibrium point based on the values of basic reproduction number R 0 . Therefore, disease eradicates from the population if basic reproduction number less than unity and disease persists in the population if basic reproduction number greater than unity. We use center manifold theorem to study the dynamical behavior of the disease-free equilibrium point for R 0 = 1 . We investigate different bifurcations such as transcritical bifurcation, backward bifurcation, saddle-node bifurcation, Hopf bifurcation and Bogdanov-Takens bifurcation of co-dimension 2. The biological significance of all types of bifurcations are described. Some numerical simulations are performed to check the reliability of our theoretical approach. Sensitivity analysis is performed to identify the influential model parameters which have most impact on the basic reproduction number of the proposed model. To control or eradicate the influence of the emerging disease, we need to control the most sensitive model parameters using necessary preventive measures. We study optimal control problem using Pontryagin's maximum principle. Finally using efficiency analysis, we determine most effective control strategy among applied controls.

Optimal control strategies for the reliable and competitive mathematical analysis of Covid-19 pandemic model

To understand dynamics of the COVID-19 disease realistically, a new SEIAPHR model has been proposed in this article where the infectious individuals have been categorized as symptomatic, asymptomatic, and super-spreaders. The model has been investigated for existence of a unique solution. To measure the contagiousness of COVID-19, reproduction numberR 0 is also computed using next generation matrix method. It is shown that the model is locally stable at disease-free equilibrium point whenR 0 < 1 and unstable forR 0 > 1 . The model has been analyzed for global stability at both of the disease-free and endemic equilibrium points. Sensitivity analysis is also included to examine the effect of parameters of the model on reproduction numberR 0 . A couple of optimal control problems have been designed to study the effect of control strategies for disease control and eradication from the society. Numerical results show that the adopted control approaches are much effective in reducing new infections.

A multicompartment mathematical model to study the dynamic behaviour of COVID-19 using vaccination as control parameter

To analyse novel coronavirus disease (COVID-19) transmission in India, this article provides an extended SEIR multicompartment model using vaccination as a control parameter. The model considers eight classes of infection: susceptible ( S ), vaccinated ( V ), exposed ( E ), asymptomatic infected ( A ), symptomatic infected ( I ), isolated ( J ), hospitalised ( H ), recovered ( R ). To begin, a mathematical study is performed to demonstrate the suggested model's uniform boundedness, epidemic equilibrium, and basic reproduction number. The findings indicate that if,R 0 < 1 , the disease-free equilibrium is locally asymptotically stable; but, if,R 0 > 1 the equilibrium is unstable. Secondly, we examine the effect on those who have received vaccinations with what are deemed optimal values. The suggested model is numerically simulated using MATLAB 14.0, and the results confirm the capacity of the proposed model to provide an accurate forecast of the progress of the epidemic in India. Finally, we examine the impact of immunisation on COVID-19 dissemination.

Dynamics of SEIR model: A case study of COVID-19 in Italy

COVID-19 takes a gigantic form worldwide in a short time from December, 2019. For this reason, World Health Organization (WHO) declared COVID-19 as a pandemic outbreak. In the early days when this outbreak began, the coronavirus spread rapidly in the community due to a lack of knowledge about the virus and the unavailability of medical facilities. Therefore it becomes a significant challenge to control the influence of the disease outbreak. In this situation, mathematical models are an important tool to employ an effective strategy in order to fight against this pandemic. To study the disease dynamics and their influence among the people, we propose a deterministic mathematical model for the COVID-19 outbreak and validate the model with real data of Italy from 15th Feb 2020 to 14th July 2020. We establish the positivity and boundedness of solutions, local stability of equilibria to examine its epidemiological relevance. Sensitivity analysis has been performed to identify the highly influential parameters which have the most impact on basic reproduction number (R0). We estimate the basic reproduction number (R0) from available data in Italy and also study effective reproduction numbers based on reported data per day from 15th Feb 2020 to 14th July 2020 in Italy. Finally, the disease control policy has been summarized in the conclusion section.

Optimal control design incorporating vaccination and treatment on six compartment pandemic dynamical system

In this paper, a mathematical model of the COVID-19 pandemic with lockdown that provides a more accurate representation of the infection rate has been analyzed. In this model, the total population is divided into six compartments: the susceptible class, lockdown class, exposed class, asymptomatic infected class, symptomatic infected class, and recovered class. The basic reproduction number (R0) is calculated using the next-generation matrix method and presented graphically based on different progression rates and effective contact rates of infective individuals. The COVID-19 epidemic model exhibits the disease-free equilibrium and endemic equilibrium. The local and global stability analysis has been done at the disease-free and endemic equilibrium based on R0. The stability analysis of the model shows that the disease-free equilibrium is both locally and globally stable when R0<1, and the endemic equilibrium is locally and globally stable when R0>1 under some conditions. A control strategy including vaccination and treatment has been studied on this pandemic model with an objective functional to minimize. Finally, numerical simulation of the COVID-19 outbreak in India is carried out using MATLAB, highlighting the usefulness of the COVID-19 pandemic model and its mathematical analysis.

Optimal Drug Regimen and Combined Drug Therapy and Its Efficacy in the Treatment of COVID-19: A Within-Host Modeling Study

The COVID-19 pandemic has resulted in more than 524 million cases and 6 million deaths worldwide. Various drug interventions targeting multiple stages of COVID-19 pathogenesis can significantly reduce infection-related mortality. The current within-host mathematical modeling study addresses the optimal drug regimen and efficacy of combination therapies in the treatment of COVID-19. The drugs/interventions considered include Arbidol, Remdesivir, Interferon (INF) and Lopinavir/Ritonavir. It is concluded that these drugs, when administered singly or in combination, reduce the number of infected cells and viral load. Four scenarios dealing with the administration of a single drug, two drugs, three drugs and all four are discussed. In all these scenarios, the optimal drug regimen is proposed based on two methods. In the first method, these medical interventions are modeled as control interventions and a corresponding objective function and optimal control problem are formulated. In this framework, the optimal drug regimen is derived. Later, using the comparative effectiveness method, the optimal drug regimen is derived based on the basic reproduction number and viral load. The average number of infected cells and viral load decreased the most when all four drugs were used together. On the other hand, the average number of susceptible cells decreased the most when Arbidol was administered alone. The basic reproduction number and viral load decreased the most when all four interventions were used together, confirming the previously obtained finding of the optimal control problem. The results of this study can help physicians make decisions about the treatment of the life-threatening COVID-19 infection.

Data analysis and prediction of the COVID-19 outbreak in the first and second waves for top 5 affected countries in the world

In this paper, we introduce a SEIATR compartmental model to analyze and predict the COVID-19 outbreak in the Top 5 affected countries in the world, namely the USA, India, Brazil, France, and Russia. The officially confirmed cases and death due to COVID-19 from the day of the official confirmation to June 30, 2021 are considered for each country. Primarily, we use the data to make a comparison between the cumulative cases and deaths due to COVID-19 among these five different countries. This analysis allows us to infer the key parameters associated with the dynamics of the disease for these five different countries. For example, the analysis reveals that the infection rate is much higher in the USA, Brazil, and France compared to that of India and Russia, while the recovery rate is found almost the same for these countries. Further, the death rate is measured higher in Brazil as opposed to India, where it is found much lower among the remaining countries. We then use the SEIART compartmental model to characterize the first and second waves of these countries, as well as to investigate and identify the influential model parameters and nature of the virus transmissibility in respective countries. Besides estimating the time-dependent reproduction number (Rt) for these countries, we also use the model to predict the peak size and the time occurring peak in respective countries. The analysis demonstrates that COVID-19 was observed to be much more infectious in the second wave than the first wave in all countries except France. The results also demonstrate that the epidemic took off very quickly in the USA, India, and Brazil compared to two other countries considered in this study. Furthermore, the prediction of the epidemic peak size and time produced by our model provides a very good agreement with the officially confirmed cases data for all countries expect Brazil.

Modeling of COVID-19 spread with self-isolation at home and hospitalized classes

This work examines the impacts of self-isolation and hospitalization on the population dynamics of the Corona-Virus Disease. We developed a new nonlinear deterministic model eight classes compartment, with self-isolation and hospitalized being the most effective tools. There are (Susceptible S C ( t ) , Exposed E ( t ) , Asymptomatic infected I A ( t ) , Symptomatic infected A S ( t ) , Self-isolation T M ( t ) , Hospitalized T H ( t ) , Healed H ( t ) , and Susceptible individuals previously infected H C ( t ) ). The expression of basic reproduction number R 0 comes from the next-generation matrix method. With suitably constructed Lyapunov functions, the global asymptotic stability of the non-endemic equilibria Σ 0 for R 0 < 1 and that of endemic equilibria Σ ∗ for R 0 > 1 are established. The computed value of R 0 = 3 . 120277403 proves the endemic level of the epidemic. The outbreak will lessen if a policy is enforced like self-isolation and hospitalization. This is related to those policies that can reduce the number of direct contacts between infected and susceptible individuals or waning immunity individuals. Various simulations are presented to appreciate self-isolation at home and hospitalized strategies if applied sensibly. By performing a global sensitivity analysis, we can obtain parameter values that affect the model through a combination of Latin Hypercube Sampling and Partial Rating Correlation Coefficients methods to determine the parameters that affect the number of reproductions and the increase in the number of COVID cases. The results obtained show that the rate of self-isolation at home and the rate of hospitalism have a negative relationship. On the other hand, infections will decrease when the two parameters increase. From the sensitivity of the results, we formulate a control model using optimal control theory by considering two control variables. The result shows that the control strategies minimize the spread of the COVID infection in the population.

Study of COVID-19 epidemiological evolution in India with a multi-wave SIR model

The global pandemic due to the outbreak of COVID-19 ravages the whole world for more than two years in which all the countries are suffering a lot since December 2019. In this article characteristics of a multi-wave SIR model have been studied which successfully explains the features of this pandemic waves in India. Origin of the multi-wave pattern in the solution of this model is explained. Stability of this model has been studied by identifying the equilibrium points as well as by finding the eigenvalues of the corresponding Jacobian matrices. In this model, a finite probability of the recovered people for becoming susceptible again is introduced which is found crucial for obtaining the oscillatory solution in other words. Which on the other hand incorporates the effect of new variants, like delta, omicron, etc in addition to the SARS-CoV-2 virus. The set of differential equations has been solved numerically in order to obtain the variation of susceptible, infected and removed populations with time. In this phenomenological study, some specific sets of parameters are chosen in order to explain the nonperiodic variation of infected population which is found necessary to capture the feature of epidemiological wave prevailing in India. The numerical estimations are compared with the actual cases along with the analytic results.

Epidemiological modeling for COVID-19 spread in India with the effect of testing

A novel coronavirus has resulted in an outbreak of viral pneumonia in China. Person-to-person transmission has been demonstrated, but, to our knowledge, the spreading of novel coronavirus takes place due to an asymptomatic carrier. Most models are not considering testing and underlying network topology that shows the spreading pattern. By failing to integrate testing into the epidemiological model, models missed a vital opportunity to better understand the role of asymptomatic infection in transmission. In this work, we propose a model considering testing as well as asymptomatic infection considering underlying network topology. We extract the transmission parameters from the data set of COVID-19 of India and apply those parameters in our proposed model. The simulation results support our theoretical derivations, which show the impact of testing and asymptomatic carrier in infection spreading.

Transmission in home environment associated with the second wave of COVID-19 pandemic in India

India has suffered from the second wave of COVID-19 pandemic since March 2021. This wave of the outbreak has been more serious than the first wave pandemic in 2020, which suggests that some new transmission characteristics may exist. COVID-19 is transmitted through droplets, aerosols, and contact with infected surfaces. Air pollutants are also considered to be associated with COVID-19 transmission. However, the roles of indoor transmission in the COVID-19 pandemic and the effects of these factors in indoor environments are still poorly understood. Our study focused on reveal the role of indoor transmission in the second wave of COVID-19 pandemic in India. Our results indicated that human mobility in the home environment had the highest relative influence on COVID-19 daily growth rate in the country. The COVID-19 daily growth rate was significantly positively correlated with the residential percent rate in most state-level areas in India. A significant positive nonlinear relationship was found when the residential percent ratio ranged from 100 to 120%. Further, epidemic dynamics modelling indicated that a higher proportion of indoor transmission in the home environment was able to intensify the severity of the second wave of COVID-19 pandemic in India. Our findings suggested that more attention should be paid to the indoor transmission in home environment. The public health strategies to reduce indoor transmission such as ventilation and centralized isolation will be beneficial to the prevention and control of COVID-19.

Solution of a dynamical memory effect COVID-19 infection system with leaky vaccination efficacy by non-singular kernel fractional derivatives

In this paper, the recent trends of COVID-19 infection spread have been studied to explore the advantages of leaky vaccination dynamics in SEVR (Susceptible Effected Vaccinated Recovered) compartmental model with the help of Caputo-Fabrizio (CF) and Atangana-Baleanu derivative in the Caputo sense (ABC) non-singular kernel fractional derivative operators with memory effect within the model to show possible long-term approaches of the infection along with limited defensive vaccine efficacy that can be designed numerically over the closed interval ranging from 0 to 1. One of the main goals is to provide a stepping information about the usefulness of the aforementioned non-singular kernel fractional approaches for a lenient case as well as a critical case in COVID-19 infection spread. Another is to investigate the effect of death rate on state variables. The estimation of death rate for state variables with suitable vaccine efficacy has a significant role in the stability of state variables in terms of basic reproduction number that is derived using next generation matrix method, and order of the fractional derivative. For non-integral orders the pandemic modeling sense viz, CF and ABC, has been compared thoroughly. Graphical presentations together with numerical results have proposed that the methodology is powerful and accurate which can provide new speculations for COVID-19 dynamical systems.

Transmission dynamics model and the coronavirus disease 2019 epidemic: applications and challenges

Since late 2019, the beginning of coronavirus disease 2019 (COVID-19) pandemic, transmission dynamics models have achieved great development and were widely used in predicting and policy making. Here, we provided an introduction to the history of disease transmission, summarized transmission dynamics models into three main types: compartment extension, parameter extension and population-stratified extension models, highlight the key contribution of transmission dynamics models in COVID-19 pandemic: estimating epidemiological parameters, predicting the future trend, evaluating the effectiveness of control measures and exploring different possibilities/scenarios. Finally, we pointed out the limitations and challenges lie ahead of transmission dynamics models.

Dynamical analysis of coronavirus disease with crowding effect, and vaccination: a study of third strain

Countries affected by the coronavirus epidemic have reported many infected cases and deaths based on world health statistics. The crowding factor, which we named "crowding effects," plays a significant role in spreading the diseases. However, the introduction of vaccines marks a turning point in the rate of spread of coronavirus infections. Modeling both effects is vastly essential as it directly impacts the overall population of the studied region. To determine the peak of the infection curve by considering the third strain, we develop a mathematical model (susceptible-infected-vaccinated-recovered) with reported cases from August 01, 2021, till August 29, 2021. The nonlinear incidence rate with the inclusion of both effects is the best approach to analyze the dynamics. The model's positivity, boundedness, existence, uniqueness, and stability (local and global) are addressed with the help of a reproduction number. In addition, the strength number and second derivative Lyapunov analysis are examined, and the model was found to be asymptotically stable. The suggested parameters efficiently control the active cases of the third strain in Pakistan. It was shown that a systematic vaccination program regulates the infection rate. However, the crowding effect reduces the impact of vaccination. The present results show that the model can be applied to other countries' data to predict the infection rate.

Age Structured Mathematical Modeling Studies on COVID-19 with respect to Combined Vaccination and Medical Treatment Strategies

In this study, we develop a mathematical model incorporating age-specific transmission dynamics of COVID-19 to evaluate the role of vaccination and treatment strategies in reducing the size of COVID-19 burden. Initially, we establish the positivity and boundedness of the solutions of the non controlled model and calculate the basic reproduction number and do the stability analysis. We then formulate an optimal control problem with vaccination and treatment as control variables and study the same. Pontryagin’s Minimum Principle is used to obtain the optimal vaccination and treatment rates. Optimal vaccination and treatment policies are analysed for different values of the weight constant associated with the cost of vaccination and different efficacy levels of vaccine. Findings from these suggested that the combined strategies (vaccination and treatment) worked best in minimizing the infection and disease induced mortality. In order to reduce COVID-19 infection and COVID-19 induced deaths to maximum, it was observed that optimal control strategy should be prioritized to the population with age greater than 40 years. Varying the cost of vaccination it was found that sufficient implementation of vaccines (more than 77 %) reduces the size of COVID-19 infections and number of deaths. The infection curves varying the efficacies of the vaccines against infection were also analysed and it was found that higher efficacy of the vaccine resulted in lesser number of infections and COVID induced deaths. The findings would help policymakers to plan effective strategies to contain the size of the COVID-19 pandemic.

A spread model of COVID-19 with some strict anti-epidemic measures

In the absence of specific drugs and vaccines, the best way to control the spread of COVID-19 is to adopt and diligently implement effective and strict anti-epidemic measures. In this paper, a mathematical spread model is proposed based on strict epidemic prevention measures and the known spreading characteristics of COVID-19. The equilibria (disease-free equilibrium and endemic equilibrium) and the basic regenerative number of the model are analyzed. In particular, we prove the asymptotic stability of the equilibria, including locally and globally asymptotic stability. In order to validate the effectiveness of this model, it is used to simulate the spread of COVID-19 in Hubei Province of China for a period of time. The model parameters are estimated by the real data related to COVID-19 in Hubei. To further verify the model effectiveness, it is employed to simulate the spread of COVID-19 in Hunan Province of China. The mean relative error serves to measure the effect of fitting and simulations. Simulation results show that the model can accurately describe the spread dynamics of COVID-19. Sensitivity analysis of the parameters is also done to provide the basis for formulating prevention and control measures. According to the sensitivity analysis and corresponding simulations, it is found that the most effective non-pharmaceutical intervention measures for controlling COVID-19 are to reduce the contact rate of the population and increase the quarantine rate of infected individuals.

Time Optimal Control Studies on COVID-19 Incorporating Adverse Events of the Antiviral Drugs

COVID -19 pandemic has resulted in more than 257 million infections and 5.15 million deaths worldwide. Several drug interventions targeting multiple stages of the pathogenesis of COVID -19 can significantly reduce induced infection and thus mortality. In this study, we first develop SIV model at within-host level by incorporating the intercellular time delay and analyzing the stability of equilibrium points. The model dynamics admits a disease-free equilibrium and an infected equilibrium with their stability based on the value of the basic reproduction number R 0. We then formulate an optimal control problem with antiviral drugs and second-line drugs as control measures and study their roles in reducing the number of infected cells and viral load. The comparative study conducted in the optimal control problem suggests that if the first-line antiviral drugs show adverse effects, considering these drugs in reduced amounts along with the second-line drugs would be very effective in reducing the number of infected cells and viral load in a COVID-19 infected patient. Later, we formulate a time-optimal control problem with the goal of driving the system from any initial state to the desired infection-free equilibrium state in finite minimal time. Using Pontryagin’s Minimum Principle, it is shown that the optimal control strategy is of the bang-bang type, with the possibility of switching between two extreme values of the optimal controls. Numerically, it is shown that the desired infection-free state is achieved in a shorter time when the higher values of the optimal controls. The results of this study may be very helpful to researchers, epidemiologists, clinicians and physicians working in this field.

Ranking non-pharmaceutical interventions against Covid-19 global pandemic using global sensitivity analysis-Effect on number of deaths

In this study, we use Global Sensitivity Analysis (GSA) to rank four non-pharmaceutical interventions (NPIs) in a deterministic compartmental model that might control Covid-19 related deaths in the United States. The NPIs are social distancing, isolation of infected individuals, identifying asymptomatically infected individuals through testing, and the use of face masks. The model uses a fear-based behavioral model that leads unmasked susceptible individuals to wear masks. The model parameters are estimated from the reported deaths for the United States of America from March 1, 2020 to November 26, 2020. Two GSA tools, the Sobol' sesntivity indices and Partial Rank Correlation Coefficient are used to obtain the rankings of the input parameters at different stages of the disease propagation. We found that social distancing and outward mask efficiency alone decreases the output uncertainty by 25-45%. Sobol' second order indices show that the combined effect of social distancing with increased mask usage and identifying and isolating asymptomatically infected individuals decreases uncertainty an additional 10%.

Numerical simulation and stability analysis of a novel reaction-diffusion COVID-19 model

In this study, a novel reaction-diffusion model for the spread of the new coronavirus (COVID-19) is investigated. The model is a spatial extension of the recent COVID-19 SEIR model with nonlinear incidence rates by taking into account the effects of random movements of individuals from different compartments in their environments. The equilibrium points of the new system are found for both diffusive and non-diffusive models, where a detailed stability analysis is conducted for them. Moreover, the stability regions in the space of parameters are attained for each equilibrium point for both cases of the model and the effects of parameters are explored. A numerical verification for the proposed model using a finite difference-based method is illustrated along with their consistency, stability and proving the positivity of the acquired solutions. The obtained results reveal that the random motion of individuals has significant impact on the observed dynamics and steady-state stability of the spread of the virus which helps in presenting some strategies for the better control of it.

The prediction of the lifetime of the new coronavirus in the USA using mathematical models

The World Health Organization (WHO) on December 31, 2019, was informed of several cases of respiratory diseases of unknown origin in the city of Wuhan in the Chinese Province of Hubei, the clinical manifestations of which were similar to those of viral pneumonia and manifested as fever, cough, and shortness of breath. And, the disease caused by the virus is named the new coronavirus disease 2019 and it will be abbreviated as 2019-nCoV and COVID-19. As of January 30, 2020, the WHO classified this epidemic as a global health emergency (Chung et al. in Radiology 295(1):202-207, 2020). It is an international real-life problem. Due to deaths, globally everyone is under fear. Now, it is the responsibility of researchers to give hope to the people. In this article, we aim to better protect people and general pandemic preparedness by predicting the lifetime of the disease-causing virus using three mathematical models. This article deals with a complex real-life problem people face all over the world, an international real-life problem. The main focus is on the USA due to large infection and death due to coronavirus and thereby the life of every individual is uncertain. The death counts of the USA from February 29 to April 22, 2020, are used in this article as a data set. The death counts of the USA are fitted by the solutions of three mathematical models and a solution to an international problem is achieved. Based on the death rate, the lifetime of the coronavirus COVID-19 is predicted as 1464.76 days from February 29, 2020. That is, after March 2024 there will be no death in the USA due to COVID-19 if everyone follows the guidelines of WHO and the advice of healthcare workers. People and government can get prepared for this situation and many lives can be saved. It is the contribution of soft computing. Finally, this article suggests several steps to control the spread and severity of the disease. The research work, the lifetime prediction presented in this article is entirely new and differs from all other articles in the literature.

COVID-19 outbreak and Urban dynamics: regional variations in India

India was the second highest COVID-19 affected country in the world with 2.1 million cases by 11th August. This study focused on the spatial transmission of the pandemic among the 640 districts in India over time, and aimed to understand the urban-centric nature of the infection. The connectivity context was emphasized that possibly had inflicted the outbreak. Using the modes of transmission data for the available cases, the diffusion of this disease was explained. Metropolitans contributed three-fourths of total cases from the beginning. The transport networks attributed significantly in transmitting the virus from the urban containment zones. Later, there was a gradual shift of infections from urban to rural areas; however, the numbers kept increasing in the former. The massive reverse migration after lockdown spiked the infected cases further. Districts with airports reported more with influx of international passengers. A profound east-west division in April with higher infections in the southern and western districts existed. By mid-May eastern India saw a steep rise in active cases. Moran's I analysis showed a low autocorrelation initially which increased over time. Hotspot clustering was observed in western Maharashtra, eastern Tamil Nadu, Gujarat and around Kolkata by the second week of August. The diffusion was due to travel, exposure to infected individuals and among the frontline workers. Spatial regression models confirmed that urbanization was positively correlated with higher incidences of infections. Transit mediums, especially rail and aviation were positively associated. These models validated the crucial role of spatial proximity in diffusion of the pandemic.

Exactly solvable SIR models, their extensions and their application to sensitive pandemic forecasting

The classic SIR model of epidemic dynamics is solved completely by quadratures, including a time integral transform expanded in a series of incomplete gamma functions. The model is also generalized to arbitrary time-dependent infection rates and solved explicitly when the control parameter depends on the accumulated infections at time t. Numerical results are presented by way of comparison. Autonomous and non-autonomous generalizations of SIR for interacting regions are also considered, including non-separability for two or more interacting regions. A reduction of simple SIR models to one variable leads us to a generalized logistic model, Richards model, which we use to fit Mexico's COVID-19 data up to day number 134. Forecasting scenarios resulting from various fittings are discussed. A critique to the applicability of these models to current pandemic outbreaks in terms of robustness is provided. Finally, we obtain the bifurcation diagram for a discretized version of Richards model, displaying period doubling bifurcation to chaos.

Global dynamics and control strategies of an epidemic model having logistic growth, non-monotone incidence with the impact of limited hospital beds

In this paper, we have considered a deterministic epidemic model with logistic growth rate of the susceptible population, non-monotone incidence rate, nonlinear treatment function with impact of limited hospital beds and performed control strategies. The existence and stability of equilibria as well as persistence and extinction of the infection have been studied here. We have investigated different types of bifurcations, namely Transcritical bifurcation, Backward bifurcation, Saddle-node bifurcation and Hopf bifurcation, at different equilibrium points under some parametric restrictions. Numerical simulation for each of the above-defined bifurcations shows the complex dynamical phenomenon of the infectious disease. Furthermore, optimal control strategies are performed using Pontryagin's maximum principle and strategies of controls are studied for two infectious diseases. Lastly using efficiency analysis we have found the effective control strategies for both cases.