Eigenvalue analysis of SARS-CoV-2 viral load data: illustration for eight COVID-19 patients

Eigenvalue analysis is an important tool in economics and nonlinear physics to analyze industrial processes and instability phenomena, respectively. A model-based eigenvalue analysis of viral load data from eight symptomatic COVID-19 patients was conducted. The eigenvalues and eigenvectors of the instabilities were determined that give rise to COVID-19. For all eight patients, it was found that the virus dynamics followed the unstable eigenvectors until the viral load reached the respective peak values. At the peak virus values, the virus dynamics branched off from the directions specified by the eigenvectors. The temporal course of the unstable eigenvalues was determined as well. For all patients, it was found that the eigenvalues switched from positive to negative values just when the virus load reached peak values. These findings suggest that the fixed, instability-related eigenvalues and eigenvectors determine initial stages of SARS-CoV-2 infections during which virus load increases. In contrast, the time-dependent eigenvalues show a sign-switching phenomenon that indicates when the virus dynamics switches from the growth stage (increasing virus load) to the decay stage (decreasing virus load). The virus dynamics model was a standard three-variable virus dynamics model frequently used in the literature.

Characterizing stages of COVID-19 epidemics: a nonlinear physics perspective based on amplitude equations

The relevant dynamics underlying COVID-19 waves is described from an amplitude space perspective. To this end, the amplitude dynamics of infected populations is considered in different stages of epidemic waves. Eigenvectors and their corresponding amplitudes are derived analytically for low-dimensional models and by means of computational methods for high-dimensional models. It is shown that the amplitudes of all eigenvectors as functions of time can be tracked through the diverse stages of COVID-19 waves featuring jumps at the stage boundaries. In particular, it is shown that under certain circumstances the initial, outbreak stage and the final, subsiding stage of an epidemic wave are primarily determined by the unstable eigenvector of the initial stage and its corresponding remnant vector of the final stage. The corresponding amplitude captures most of the dynamics of the emerging and subsiding epidemics such that the problem at hand effectively becomes one dimensional leading to a dramatic reduction of the complexity of the problem at hand. Explicitly demonstrated for the first-wave COVID-19 epidemics of the year 2020 in the state of New York and Pakistan are given.

SARS-coronavirus-2 infections: biological instabilities characterized by order parameters

A four-variable virus dynamics TIIV model was considered that involves infected cells in an eclipse phase. The state space description of the model was transferred into an amplitude space description which is the appropriate general, nonlinear physics framework to describe instabilities. In this context, the unstable eigenvector or order parameter of the model was determined. Subsequently, a model-based analysis of viral load data from eight symptomatic COVID-19 patients was conducted. For all patients, it was found that the initial SARS-CoV-2 infection evolved along the respective patient-specific order parameter, as expected by theoretical considerations. The order parameter amplitude that described the initial virus multiplication showed doubling times between 30 minutes and 3 hours. Peak viral loads of patients were linearly related to the amplitudes of the patient order parameters. Finally, it was found that the patient order parameters determined qualitatively and quantitatively the relationships between the increases in virus-producing infected cells and infected cells in the eclipse phase. Overall, the study echoes the 40 years old suggestion by Mackey and Glass to consider diseases as instabilities.

Rise and Decay of the COVID-19 Epidemics in the USA and the State of New York in the First Half of 2020: A Nonlinear Physics Perspective Yielding Novel Insights

As of December 2020, since the beginning of the year 2020, the COVID-19 pandemic has claimed worldwide more than 1 million lives and has changed human life in unprecedented ways. Despite the fact that the pandemic is far from over, several countries managed at least temporarily to make their first-wave COVID-19 epidemics to subside to relatively low levels. Combining an epidemiological compartment model and a stability analysis as used in nonlinear physics and synergetics, it is shown how the first-wave epidemics in the state of New York and nationwide in the USA developed through three stages during the first half of the year 2020. These three stages are the outbreak stage, the linear stage, and the subsiding stage. Evidence is given that the COVID-19 outbreaks in these two regions were due to instabilities of the COVID-19 free states of the corresponding infection dynamical systems. It is shown that from stage 1 to stage 3, these instabilities were removed, presumably due to intervention measures, in the sense that the COVID-19 free states were stabilized in the months of May and June in both regions. In this context, stability parameters and key directions are identified that characterize the infection dynamics in the outbreak and subsiding stages. Importantly, it is shown that the directions in combination with the sign-switching of the stability parameters can explain the observed rise and decay of the epidemics in the state of New York and the USA. The nonlinear physics perspective provides a framework to obtain insights into the nature of the COVID-19 dynamics during outbreak and subsiding stages and allows to discuss possible impacts of intervention measures. For example, the directions can be used to determine how different populations (e.g., exposed versus symptomatic individuals) vary in size relative to each other during the course of an epidemic. Moreover, the timeline of the computationally obtained stages can be compared with the history of the implementation of intervention measures to discuss the effectivity of such measures.

SEIR order parameters and eigenvectors of the three stages of completed COVID-19 epidemics: with an illustration for Thailand January to May 2020

By end of October 2020, the COVID-19 pandemic has taken a tragic toll of 1150 000 lives and this number is expected to increase. Despite the pandemic is raging in most parts of the world, in a few countries COVID-19 epidemics subsided due to successful implementations of intervention measures. A unifying perspective of the beginnings, middle stages, and endings of such completed COVID-19 epidemics is developed based on the order parameter and eigenvalue concepts of nonlinear physics, in general, and synergetics, in particular. To this end, a standard susceptible-exposed-infected-recovered (SEIR) epidemiological model is used. It is shown that COVID-19 epidemic outbreaks follow a suitably defined SEIR order parameter. Intervention measures switch the eigenvalue of the order parameter from a positive to a negative value, and in doing so, stabilize the COVID-19 disease-free state. The subsiding of COVID-19 epidemics eventually follows the remnant of the order parameter of the infection dynamical system. These considerations are illustrated for the COVID-19 epidemic in Thailand from January to May 2020. The decay of effective contact rates throughout the three epidemic stages is demonstrated. Evidence for the sign-switching of the dominant eigenvalue is given and the order parameter and its stage-3 remnant are identified. The presumed impacts of interventions measures implemented in Thailand are discussed in this context.

Modeling the effects of prosocial awareness on COVID-19 dynamics: Case studies on Colombia and India

The ongoing COVID-19 pandemic has affected most of the countries on Earth. It has become a pandemic outbreak with more than 50 million confirmed infections and above 1 million deaths worldwide. In this study, we consider a mathematical model on COVID-19 transmission with the prosocial awareness effect. The proposed model can have four equilibrium states based on different parametric conditions. The local and global stability conditions for awareness-free, disease-free equilibrium are studied. Using Lyapunov function theory and LaSalle invariance principle, the disease-free equilibrium is shown globally asymptotically stable under some parametric constraints. The existence of unique awareness-free, endemic equilibrium and unique endemic equilibrium is presented. We calibrate our proposed model parameters to fit daily cases and deaths from Colombia and India. Sensitivity analysis indicates that the transmission rate and the learning factor related to awareness of susceptibles are very crucial for reduction in disease-related deaths. Finally, we assess the impact of prosocial awareness during the outbreak and compare this strategy with popular control measures. Results indicate that prosocial awareness has competitive potential to flatten the COVID-19 prevalence curve.

COVID-19: Analytic results for a modified SEIR model and comparison of different intervention strategies

The Susceptible-Exposed-Infected-Recovered (SEIR) epidemiological model is one of the standard models of disease spreading. Here we analyse an extended SEIR model that accounts for asymptomatic carriers, believed to play an important role in COVID-19 transmission. For this model we derive a number of analytic results for important quantities such as the peak number of infections, the time taken to reach the peak and the size of the final affected population. We also propose an accurate way of specifying initial conditions for the numerics (from insufficient data) using the fact that the early time exponential growth is well-described by the dominant eigenvector of the linearized equations. Secondly we explore the effect of different intervention strategies such as social distancing (SD) and testing-quarantining (TQ). The two intervention strategies (SD and TQ) try to reduce the disease reproductive number, R 0 , to a target value R 0 target < 1 , but in distinct ways, which we implement in our model equations. We find that for the same R 0 target < 1 , TQ is more efficient in controlling the pandemic than SD. However, for TQ to be effective, it has to be based on contact tracing and our study quantifies the required ratio of tests-per-day to the number of new cases-per-day. Our analysis shows that the largest eigenvalue of the linearised dynamics provides a simple understanding of the disease progression, both pre- and post- intervention, and explains observed data for many countries. We apply our results to the COVID data for India to obtain heuristic projections for the course of the pandemic, and note that the predictions strongly depend on the assumed fraction of asymptomatic carriers.

COVID-19 interventions in some European countries induced bifurcations stabilizing low death states against high death states: An eigenvalue analysis based on the order parameter concept of synergetics

Taking a dynamical systems perspective, COVID-19 infections are assumed to spread out in a human population via an instability. Conversely, government interventions to reduce the spread of the disease and the number of fatalities may induce a bifurcation that stabilizes a desirable state with low numbers of COVID-19 cases and associated deaths. The key characteristic feature of an infection dynamical system in this context is the eigenvalue that determines the stability of the states under consideration and is known in synergetics as the order parameter eigenvalue. Using a SEIR-like infection disease model, the relevant order parameter and its eigenvalue are determined. A three stage methodology is proposed to track and estimate the eigenvalue through time. The method is applied to COVID-19 infection data reported from 20 European countries during the period of January 1, 2020 to June 15. It is shown that in 15 out of the 20 countries the eigenvalue switched its sign suggesting that during the reporting period an intervention bifurcation took place that stabilized the desirable low death state. It is shown that the eigenvalue analysis also allows for a ranking of countries by the degree of the stability of the infection-free state. For the investigated countries, Ireland was found to exhibit the most stable infection-free state. Finally, a six point classification scheme is suggested with groups 5 and 6 including countries that failed to stabilize the desirable infection-free low death state. In doing so, tools for assessing the effectiveness of government interventions are provided that are at the heart of bifurcation theory, in general, and synergetics, in particular.