Communications in Nonlinear Science and Numerical Simulation

Dynamics and profiles of a degenerated reaction-diffusion host-pathogen model with apparent and inapparent infection period
Wang J and Lu H
Inapparent infection plays an important role in the disease spread, which is an infection by a pathogen that causes few or no signs or symptoms of infection in the host. Many pathogens, including HIV, typhoid fever, and coronaviruses such as COVID-19 spread in their host populations through inapparent infection. In this paper, we formulated a degenerated reaction-diffusion host-pathogen model with multiple infection period. We split the infectious individuals into two distinct classes: apparent infectious individuals and inapparent infectious individuals, coming from exposed individuals with a ratio of and , respectively. Some preliminary results and threshold-type results are achieved by detailed mathematical analysis. We also investigate the asymptotic profiles of the positive steady state (PSS) when the diffusion rate of susceptible individuals approaches zero or infinity. When all parameters are all constants, the global attractivity of the constant endemic equilibrium is established. It is verified by numerical simulations that spatial heterogeneity of the transmission rates can enhance the intensity of an epidemic. Especially, the transmission rate of inapparent infectious individuals significantly increases the risk of disease transmission, compared to that of apparent infectious individuals and pathogens in the environment, and we should pay special attentions to how to regulate the inapparent infectious individuals for disease control and prevention, which is consistent with the result on the sensitive analysis to the transmission rates through the normalized forward sensitivity index. We also find that disinfection of the infected environment is an important way to prevent and eliminate the risk of environmental transmission.
Assessing potential insights of an imperfect testing strategy: Parameter estimation and practical identifiability using early COVID-19 data in India
Bugalia S and Tripathi JP
A deterministic model with testing of infected individuals has been proposed to investigate the potential consequences of the impact of testing strategy. The model exhibits global dynamics concerning the disease-free and a unique endemic equilibrium depending on the basic reproduction number when the recruitment of infected individuals is zero; otherwise, the model does not have a disease-free equilibrium, and disease never dies out in the community. Model parameters have been estimated using the maximum likelihood method with respect to the data of early COVID-19 outbreak in India. The practical identifiability analysis shows that the model parameters are estimated uniquely. The consequences of the testing rate for the weekly new cases of early COVID-19 data in India tell that if the testing rate is increased by 20% and 30% from its baseline value, the weekly new cases at the peak are decreased by 37.63% and 52.90%; and it also delayed the peak time by four and fourteen weeks, respectively. Similar findings are obtained for the testing efficacy that if it is increased by 12.67% from its baseline value, the weekly new cases at the peak are decreased by 59.05% and delayed the peak by 15 weeks. Therefore, a higher testing rate and efficacy reduce the disease burden by tumbling the new cases, representing a real scenario. It is also obtained that the testing rate and efficacy reduce the epidemic's severity by increasing the final size of the susceptible population. The testing rate is found more significant if testing efficacy is high. Global sensitivity analysis using partial rank correlation coefficients (PRCCs) and Latin hypercube sampling (LHS) determine the key parameters that must be targeted to worsen/contain the epidemic.
Traveling wave speed and profile of a "go or grow" glioblastoma multiforme model
Tursynkozha A, Kashkynbayev A, Shupeyeva B, Rutter EM and Kuang Y
Glioblastoma multiforme (GBM) is a fast-growing and deadly brain tumor due to its ability to aggressively invade the nearby brain tissue. A host of mathematical models in the form of reaction-diffusion equations have been formulated and studied in order to assist clinical assessment of GBM growth and its treatment prediction. To better understand the speed of GBM growth and form, we propose a two population reaction-diffusion GBM model based on the 'go or grow' hypothesis. Our model is validated by data and assumes that tumor cells are more likely to leave and search for better locations when resources are more limited at their current positions. Our findings indicate that the tumor progresses slower than the simpler Fisher model, which is known to overestimate GBM progression. Moreover, we obtain accurate estimations of the traveling wave solution profiles under several plausible GBM cell switching scenarios by applying the approximation method introduced by Canosa.
Homogenisation for the monodomain model in the presence of microscopic fibrotic structures
Lawson BAJ, Dos Santos RW, Turner IW, Bueno-Orovio A, Burrage P and Burrage K
Computational models in cardiac electrophysiology are notorious for long runtimes, restricting the numbers of nodes and mesh elements in the numerical discretisations used for their solution. This makes it particularly challenging to incorporate structural heterogeneities on small spatial scales, preventing a full understanding of the critical arrhythmogenic effects of conditions such as cardiac fibrosis. In this work, we explore the technique of homogenisation by volume averaging for the inclusion of non-conductive micro-structures into larger-scale cardiac meshes with minor computational overhead. Importantly, our approach is not restricted to periodic patterns, enabling homogenised models to represent, for example, the intricate patterns of collagen deposition present in different types of fibrosis. We first highlight the importance of appropriate boundary condition choice for the closure problems that define the parameters of homogenised models. Then, we demonstrate the technique's ability to correctly upscale the effects of fibrotic patterns with a spatial resolution of into much larger numerical mesh sizes of 100- . The homogenised models using these coarser meshes correctly predict critical pro-arrhythmic effects of fibrosis, including slowed conduction, source/sink mismatch, and stabilisation of re-entrant activation patterns. As such, this approach to homogenisation represents a significant step towards whole organ simulations that unravel the effects of microscopic cardiac tissue heterogeneities.
The stochastic -SEIHRD model: Adding randomness to the COVID-19 spread
Leitao Á and Vázquez C
In this article we mainly extend a newly introduced deterministic model for the COVID-19 disease to a stochastic setting. More precisely, we incorporated randomness in some coefficients by assuming that they follow a prescribed stochastic dynamics. In this way, the model variables are now represented by stochastic process, that can be simulated by appropriately solving the system of stochastic differential equations. Thus, the model becomes more complete and flexible than the deterministic analogous, as it incorporates additional uncertainties which are present in more realistic situations. In particular, confidence intervals for the main variables and worst case scenarios can be computed.
On the economic growth equilibria during the Covid-19 pandemic
Bischi GI, Grassetti F and Sanchez Carrera EJ
The aim of this paper is to study the effects of the Covid-19 pandemic suppression policies (i.e. containment measures or lockdowns) on labor supply, capital accumulation, and so the economic growth. We merge an epidemic SIS population model and a Solow's type growth model, i.e. we propose a fusion between economics and epidemiology. We show the creation and the destruction of economic growth equilibria driven by the suppression policies and by the severity of the disease. The dynamic stability properties of the equilibria are mainly determined by (i) the stringency of the suppression policies, (ii) the proportion of infected workers, (iii) the recovery rate of workers, and (iv) the economy's saving rate. Thus, economies can fall into the stable equilibrium of the poverty trap if the propensity to save is low and the economic policies that reduce the spread of infection are severe enough with high levels of infection and low rates of illness recovery. Otherwise, with high savings rates and if the suppression policies perform in such a way that infection levels are low and recovery rates are high, then the economies converge towards the equilibrium of high economic growth with capital accumulation. The scenario is rather complex since there is a multiplicity of equilibria such that economies can be in one scenario or another, characterized by stability or (structural) instability, i.e. bifurcation paths. Numerical simulations corroborate our results.
Robust optimal control of compartmental models in epidemiology: Application to the COVID-19 pandemic
Olivares A and Staffetti E
In this paper, a spectral approach is used to formulate and solve robust optimal control problems for compartmental epidemic models, allowing the uncertainty propagation through the optimal control model to be represented by a polynomial expansion of its stochastic state variables. More specifically, a statistical moment-based polynomial chaos expansion is employed. The spectral expansion of the stochastic state variables allows the computation of their main statistics to be carried out, resulting in a compact and efficient representation of the variability of the optimal control model with respect to its random parameters. The proposed robust formulation provides the designers of the optimal control strategy of the epidemic model the capability to increase the predictability of the results by simply adding upper bounds on the variability of the state variables. Moreover, this approach yields a way to efficiently estimate the probability distributions of the stochastic state variables and conduct a global sensitivity analysis. To show the practical implementation of the proposed approach, a mathematical model of COVID-19 transmission is considered. The numerical results show that the spectral approach proposed to formulate and solve robust optimal control problems for compartmental epidemic models provides healthcare systems with a valuable tool to mitigate and control the impact of infectious diseases.
Epidemic spreading with migration in networked metapopulation
Wang NN, Wang YJ, Qiu SH and Di ZR
Migration plays a crucial role in epidemic spreading, and its dynamic can be studied by metapopulation model. Instead of the uniform mixing hypothesis, we adopt networked metapopulation to build the model of the epidemic spreading and the individuals' migration. In these populations, individuals are connected by contact network and populations are coupled by individuals migration. With the network mean-field and the gravity law of migration, we establish the N-seat intertwined SIR model and obtain its basic reproduction number . Meanwhile, we devise a non-markov Node-Search algorithm for model statistical simulations. Through the static network migration ansatz and formula, we discover that migration will not directly increase the epidemic replication capacity. But when , the migration will make the susceptive population evolve from metastable state (disease-free equilibrium) to stable state (endemic equilibrium), and then increase the influence area of epidemic. Re-evoluting the epidemic outbreak in Wuhan, top 94 cities empirical data validate the above mechanism. In addition, we estimate that the positive anti-epidemic measures taken by the Chinese government may have reduced 4 million cases at least during the first wave of COVID-19, which means those measures, such as the epidemiological investigation, nucleic acid detection in medium-high risk areas and isolation of confirmed cases, also play a significant role in preventing epidemic spreading after travel restriction between cities.
Modeling calcium dynamics in neurons with endoplasmic reticulum: existence, uniqueness and an implicit-explicit finite element scheme
Guan Q and Queisser G
Like many other biological processes, calcium dynamics in neurons containing an endoplasmic reticulum is governed by diffusion-reaction equations on interface-separated domains. Interface conditions are typically described by systems of ordinary differential equations that provide fluxes across the interfaces. Using the calcium model as an example of this class of ODE-flux boundary interface problems, we prove the existence, uniqueness and boundedness of the solution by applying comparison theorem, fundamental solution of the parabolic operator and a strategy used in Picard's existence theorem. Then we propose and analyze an efficient implicit-explicit finite element scheme which is implicit for the parabolic operator and explicit for the nonlinear terms. We show that the stability does not depend on the spatial mesh size. Also the optimal convergence rate in norm is obtained. Numerical experiments illustrate the theoretical results.
Nonlinear growth and mathematical modelling of COVID-19 in some African countries with the Atangana-Baleanu fractional derivative
Kolebaje OT, Vincent OR, Vincent UE and McClintock PVE
We analyse the time-series evolution of the cumulative number of confirmed cases of COVID-19, the novel coronavirus disease, for some African countries. We propose a mathematical model, incorporating non-pharmaceutical interventions to unravel the disease transmission dynamics. Analysis of the stability of the model's steady states was carried out, and the reproduction number , a vital key for flattening the time-evolution of COVID-19 cases, was obtained by means of the next generation matrix technique. By dividing the time evolution of the pandemic for the cumulative number of confirmed infected cases into different regimes or intervals, hereafter referred to as , numerical simulations were performed to fit the proposed model to the cumulative number of confirmed infections for different phases of COVID-19 during its first wave. The estimated declined from 2.452-9.179 during the first phase of the infection to 1.374-2.417 in the last phase. Using the Atangana-Baleanu fractional derivative, a fractional COVID-19 model is proposed and numerical simulations performed to establish the dependence of the disease dynamics on the order of the fractional derivatives. An elasticity and sensitivity analysis of was carried out to determine the most significant parameters for combating the disease outbreak. These were found to be the effective disease transmission rate, the disease diagnosis or case detection rate, the proportion of susceptible individuals taking precautions, and the disease infection rate. Our results show that if the disease infection rate is less than 0.082/day, then is always less than 1; and if at least 55.29% of the susceptible population take precautions such as regular hand washing with soap, use of sanitizers, and the wearing of face masks, then the reproduction number remains below unity irrespective of the disease infection rate. Keeping values below unity leads to a decrease in COVID-19 prevalence.
Network synchronization, stability and rhythmic processes in a diffusive mean-field coupled SEIR model
Verma T and Gupta AK
Connectivity and rates of movement have profound effect on the persistence and extinction of infectious diseases. The emerging disease spread rapidly, due to the movement of infectious persons to some other regions, which has been witnessed in case of novel coronavirus disease 2019 (COVID-19). So, the networks and the epidemiology of directly transmitted infectious diseases are fundamentally linked. Motivated by the recent empirical evidence on the dispersal of infected individuals among the patches, we present the epidemic model SEIR (Susceptible-Exposed-Infected-Recovered) in which the population is divided into patches which form a network and the patches are connected through mean-field diffusive coupling. The corresponding unstable epidemiology classes will be synchronized and achieve stable state when the patches are coupled. Apart from synchronization and stability, the coupled model enables a range of rhythmic processes such as birhythmicity and rhythmogenesis which have not been investigated in epidemiology. The stability of Disease Free Equilibrium (or Endemic Equilibrium) is attained through cessation of oscillation mechanism namely Oscillation Death (OD) and Amplitude Death (AD). Corresponding to identical and non-identical epidemiology classes of patches, the different steady states are obtained and its transition is taking place through Hopf and transcritical bifurcation.
Modeling the impact of SARS-CoV-2 variants and vaccines on the spread of COVID-19
Ramos AM, Vela-Pérez M, Ferrández MR, Kubik AB and Ivorra B
The continuous mutation of SARS-CoV-2 opens the possibility of the appearance of new variants of the virus with important differences in its spreading characteristics, mortality rates, etc. On 14 December 2020, the United Kingdom reported a potentially more contagious coronavirus variant, present in that country, which is referred to as VOC 202012/01. On 18 December 2020, the South African government also announced the emergence of a new variant in a scenario similar to that of the UK, which is referred to as variant 501.V2. Another important milestone regarding this pandemic was the beginning, in December 2020, of vaccination campaigns in several countries. There are several vaccines, with different characteristics, developed by various laboratories and research centers. A natural question arises: what could be the impact of these variants and vaccines on the spread of COVID-19? Many models have been proposed to simulate the spread of COVID-19 but, to the best of our knowledge, none of them incorporates the effects of potential SARS-CoV-2 variants together with the vaccines in the spread of COVID-19. We develop here a -SVEIHQRD mathematical model able to simulate the possible impact of this type of variants and of the vaccines, together with the main mechanisms influencing the disease spread. The model may be of interest for policy makers, as a tool to evaluate different possible future scenarios. We apply the model to the particular case of Italy (as an example of study case), showing different outcomes. We observe that the vaccines may reduce the infections, but they might not be enough for avoiding a new wave, with the current expected vaccination rates in that country, if the control measures are relaxed. Furthermore, a more contagious variant could increase significantly the cases, becoming the most common way of infection. We show how, even with the pandemic cases slowing down (with an effective reproduction number less than 1) and the disease seeming to be under control, the effective reproduction number of just the new variant may be greater than 1 and, eventually, the number of infections would increase towards a new disease wave. Therefore, a rigorous follow-up of the evolution of the number of infections with any potentially more dangerous new variant is of paramount importance at any stage of the pandemic.
A fractional-order compartmental model for the spread of the COVID-19 pandemic
Biala TA and Khaliq AQM
We propose a time-fractional compartmental model (SEI I HRD) comprising of the susceptible, exposed, infected (asymptomatic and symptomatic), hospitalized, recovered and dead population for the COVID-19 pandemic. We study the properties and dynamics of the proposed model. The conditions under which the disease-free and endemic equilibrium points are asymptotically stable are discussed. Furthermore, we study the sensitivity of the parameters and use the data from Tennessee state (as a case study) to discuss identifiability of the parameters of the model. The non-negative parameters in the model are obtained by solving inverse problems with empirical data from California, Florida, Georgia, Maryland, Tennessee, Texas, Washington and Wisconsin. The basic reproduction number is seen to be slightly above the critical value of one suggesting that stricter measures such as the use of face-masks, social distancing, contact tracing, and even longer stay-at-home orders need to be enforced in order to mitigate the spread of the virus. As stay-at-home orders are rescinded in some of these states, we see that the number of cases began to increase almost immediately and may continue to rise until the end of the year 2020 unless stricter measures are taken.
Stability analysis in COVID-19 within-host model with immune response
Almocera AES, Quiroz G and Hernandez-Vargas EA
The 2019 coronavirus disease (COVID-19) is now a global pandemic. Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is the causative pathogen of COVID-19. Here, we study an in-host model that highlights the effector T cell response to SARS-CoV-2. The stability of a unique positive equilibrium point, with viral load suggests that the virus may replicate fast enough to overcome T cell response and cause infection. This overcoming is the bifurcation point, near which the orders of magnitude for can be sensitive to numerical changes in the parameter values. Our work offers a mathematical insight into how SARS-CoV-2 causes the disease.
Nonlinear self-organized population dynamics induced by external selective nonlocal processes
Tumbarell Aranda O, Penna ALA and Oliveira FA
Self-organization evolution of a population is studied considering generalized reaction-diffusion equations. We proposed a model based on non-local operators that has several of the equations traditionally used in research on population dynamics as particular cases. Then, employing a relatively simple functional form of the non-local kernel, we determined the conditions under which the analyzed population develops spatial patterns, as well as their main characteristics. Finally, we established a relationship between the developed model and real systems by making simulations of bacterial populations subjected to non-homogeneous lighting conditions. Our proposal reproduces some of the experimental results that other approaches considered previously had not been able to obtain.
On the uncertainty of real-time predictions of epidemic growths: A COVID-19 case study for China and Italy
Alberti T and Faranda D
While COVID-19 is rapidly propagating around the globe, the need for providing real-time forecasts of the epidemics pushes fits of dynamical and statistical models to available data beyond their capabilities. Here we focus on statistical predictions of COVID-19 infections performed by fitting asymptotic distributions to actual data. By taking as a case-study the epidemic evolution of total COVID-19 infections in Chinese provinces and Italian regions, we find that predictions are characterized by large uncertainties at the early stages of the epidemic growth. Those uncertainties significantly reduce after the epidemics peak is reached. Differences in the uncertainty of the forecasts at a regional level can be used to highlight the delay in the spread of the virus. Our results warn that long term extrapolation of epidemics counts must be handled with extreme care as they crucially depend not only on the quality of data, but also on the stage of the epidemics, due to the intrinsically non-linear nature of the underlying dynamics. These results suggest that real-time epidemiological projections should include wide uncertainty ranges and urge for the needs of compiling high-quality datasets of infections counts, including asymptomatic patients.
Extreme events and emergency scales
Smirnov V, Ma Z and Volchenkov D
An event is extreme if its magnitude exceeds the threshold. A choice of a threshold is subject to uncertainty caused by a method, the size of available data, a hypothesis on statistics, etc. We assess the degree of uncertainty by the Shannon's entropy calculated on the probability that the threshold changes at any given time. If the amount of data is not sufficient, an observer is in the state of Lewis Carroll's who said "When you say hill, I could show you hills, in comparison with which you'd call that a valley". If we have enough data, the uncertainty curve peaks at two values clearly separating the magnitudes of events into three emergency scales: . Our approach to defining the emergency scale is validated by 39 years of Standard and Poor's 500 (S&P500) historical data.
Optimal control of a fractional order model for granular SEIR epidemic with uncertainty
Dong NP, Long HV and Khastan A
In this study, we present a general formulation for the optimal control problem to a class of fuzzy fractional differential systems relating to SIR and SEIR epidemic models. In particular, we investigate these epidemic models in the uncertain environment of fuzzy numbers with the rate of change expressed by granular Caputo fuzzy fractional derivatives of order  ∈ (0, 1]. Firstly, the existence and uniqueness of solution to the abstract fractional differential systems with fuzzy parameters and initial data are proved. Next, the optimal control problem for this fractional system is proposed and a necessary condition for the optimality is obtained. Finally, some examples of the fractional SIR and SEIR models are presented and tested with real data extracted from COVID-19 pandemic in Italy and South Korea.
Mathematical modeling of the spread of the coronavirus disease 2019 (COVID-19) taking into account the undetected infections. The case of China
Ivorra B, Ferrández MR, Vela-Pérez M and Ramos AM
In this paper we develop a mathematical model for the spread of the coronavirus disease 2019 (COVID-19). It is a new -SEIHRD model (not a SIR, SEIR or other general purpose model), which takes into account the known special characteristics of this disease, as the existence of infectious undetected cases and the different sanitary and infectiousness conditions of hospitalized people. In particular, it includes a novel approach that considers the fraction of detected cases over the real total infected cases, which allows to study the importance of this ratio on the impact of COVID-19. The model is also able to estimate the needs of beds in hospitals. It is complex enough to capture the most important effects, but also simple enough to allow an affordable identification of its parameters, using the data that authorities report on this pandemic. We study the particular case of China (including Chinese Mainland, Macao, Hong-Kong and Taiwan, as done by the World Health Organization in its reports on COVID-19), the country spreading the disease, and use its reported data to identify the model parameters, which can be of interest for estimating the spread of COVID-19 in other countries. We show a good agreement between the reported data and the estimations given by our model. We also study the behavior of the outputs returned by our model when considering incomplete reported data (by truncating them at some dates before and after the peak of daily reported cases). By comparing those results, we can estimate the error produced by the model when identifying the parameters at early stages of the pandemic. Finally, taking into account the advantages of the novelties introduced by our model, we study different scenarios to show how different values of the percentage of detected cases would have changed the global magnitude of COVID-19 in China, which can be of interest for policy makers.
Excitable dynamics in neural and cardiac systems
Barrio R, Coombes S, Desroches M, Fenton F, Luther S and Pueyo E
The application of mathematics, physics and engineering to medical research is continuously growing; interactions among these disciplines have become increasingly important and have contributed to an improved understanding of clinical and biological phenomena, with implications for disease prevention, diagnosis and treatment. This special issue presents examples of this synergy, with a particular focus on the investigation of cardiac and neural excitability. This issue includes 24 original research papers and covers a broad range of topics related to the physiological and pathophysiological function of the brain and the heart. Studies span scales from isolated neurons and small networks of neurons to whole-organ dynamics for the brain and from cardiac subcellular domains and cardiomyocytes to one-dimensional tissues for the heart. This preface is part of the Special Issue on "Excitable Dynamics in Neural and Cardiac Systems".
Key aspects for effective mathematical modelling of fractional-diffusion in cardiac electrophysiology: a quantitative study
Cusimano N, Gizzi A, Fenton FH, Filippi S and Gerardo-Giorda L
Microscopic structural features of cardiac tissue play a fundamental role in determining complex spatio-temporal excitation dynamics at the macroscopic level. Recent efforts have been devoted to the development of mathematical models accounting for non-local spatio-temporal coupling able to capture these complex dynamics without the need of resolving tissue heterogeneities down to the micro-scale. In this work, we analyse in detail several important aspects affecting the overall predictive power of these modelling tools and provide some guidelines for an effective use of space-fractional models of cardiac electrophysiology in practical applications. Through an extensive computational study in simplified computational domains, we highlight the robustness of models belonging to different categories, i.e., physiological and phenomenological descriptions, against the introduction of non-locality, and lay down the foundations for future research and model validation against experimental data. A modern genetic algorithm framework is used to investigate proper parameterisations of the considered models, and the crucial role played by the boundary assumptions in the considered settings is discussed. Several numerical results are provided to support our claims.