key: cord-1022172-ynr9hhc6 authors: Natiq, Hayder; Saha, Asit title: In search of COVID-19 transmission through an infected prey date: 2022-03-16 journal: Eur Phys J Spec Top DOI: 10.1140/epjs/s11734-022-00429-5 sha: 9a08b3246c6f47c6ce0d29def6ff171d638aa5dc doc_id: 1022172 cord_uid: ynr9hhc6 This paper considers a nonlinear dynamical model of an ecosystem, which has been established through combining the classical Lotka–Volterra model with the classic SIR model. This nonlinear system consists of a generalist predator that subsists on two prey species in which disease is becoming endemic in one of them. The dynamical analysis methods prove that the system has a chaotic attractor and extreme multistability behavior, where there are infinitely many attractors that coexist under certain conditions. The occurrence of extreme multistability demonstrates the high sensitivity of the system for the initial conditions, which means that tiny changes in the original prey species could enlarge and be widespread, and that could confirm through studying the complexity of the time series of the system’s variables. Simulation results of the sample entropy algorithm show that the changes in the system’s variables expand over time. It is reasonable now to consider the endemic in the prey species of the system could evolve to be pandemic such as COVID-19. Consequently, our results could provide a foresight about the unpredictability of the COVID-19 outbreak in its original host species as well as after the transmission to other species such as humans. In the second week of December 2019, unknown viral infection was identified in a small local wild animal market in Wuhan city, Hubei province in China [1] . The Chinese Center for Disease Control and Prevention (CCDC) identified this infection as a novel coronavirus infection, which is announced as COVID-19 [2, 3] and spreads with several deadly waves [4, 5] . COVID-19 is the new virus in the coronavirus family and mostly spreads through transmission among human. In this family, it has been reported that most of these types of viruses have an origin in animals. The infection spreading from animal to human can also be observed as a starting points of many outbreaks, from influenza to HIV and from SARS to COVID-19 [6] . The coronavirus causing COVID-19 is highly related to those that caused the SARS (severe acute respiratory syndrome) pandemic in 2003 [7] . Despite speculation, it is early to predict whether COVID-19 will disappear or stay with us permanently like the flu [8] . However, a predator-prey model was published recently by A. Eilersen et. al [9] . In this model, there is 100% fatal disease in both prey and predator species. In fact, the dynamics of predator-prey-pathogen interaca e-mail: asit sah123@rediffmail.com (corresponding author) tions in general have received significant attention in recent decades [10] . Most attention has been given to the interaction between predators and single-host epidemics [11] [12] [13] . Therefore, studying a dynamical model, which describes the interplay between predation and multi-host fatal disease, could be a reasonable option to understand the unpredictable behavior of COVID- 19 . Numerous predator-prey models in the literature show chaotic and quasi-periodic attractors for given sets of parameters and initial conditions. For discrete-time systems, the 1D logistic map [14] , which describes animal reproduction, is one of the most well-known examples of chaos. Fewer instances of chaos were found in ecological models in continuous time. However, over and above investigations of the dynamics of the nonlinear system revealed some unexpected results. It observed that a dynamical system can exhibit ambiguous behavior, which is more complex than existing chaos in the system [15] [16] [17] . The coexisting attractors or multistability is the name of this behavior [18, 19] . Multistability has reported in both continuous and discrete systems [20, 21] . In this article, we have considered the 4D dynamical model of an ecosystem, which constructs by combining the classical Lotka-Volterra model with the classic SIR model [9] . This model describes the interplay between a generalist predator subsists on two prey species in which disease is becoming endemic in one of them. Several numerical analysis methods, including phase space, Lyapunov exponents, and bifurcation diagram, are used in the research to investigate the system's dynamics. This research, through numerical simulations, demonstrates that the predator-prey system can exhibit various types of multistability behaviors, in which the coexistence of infinity many quasi-periodic attractors, and the coexistence of infinity many chaotic attractors can be observed for the same parameters values. Moreover, the complex multistability behavior of coexisting chaotic and quasi-periodic is identified for specific parameters values. The different types of multistability behaviors reflect the high sensitivity of the system. Thus, tiny changes in the original prey species could enlarge and be widespread, and that could confirm through studying the complexity of the time series of the system's variables. In this regard, a contour plot based on the sample entropy algorithm has been applied to estimate the complexity of each variable. The complexity results revealed that the coexisting attractors of the same set of parameters exhibit different complexity values; hence, the changes in the system's variables expand over time. Therefore, it is reasonable now to consider the endemic in the prey species of the system could evolve to be pandemic such as COVID-19. The rest of this paper is organized as follows. In Sect. 2, we briefly describe the dynamical model of the combination of Lotka-Volterra model with the classic SIR model, and we then investigate its stability and parameters-related chaotic regions. Section 3 provides an overview of the dynamics changes in the system when initials-boosted, and this can illustrate the apparent unpredictability of epidemic outbreaks in prey species and humans. In Sect. 4, we conduct numerical experiment based on sample entropy algorithm to illustrate the complexity for the time series of predator and prey species. This section studies a 4D dynamical model of an ecosystem [9] , which consists of a generalist predator subsists on two prey species in which a disease is becoming endemic in one of them. Generally, a disease can be spillover from an infected prey to human when the infection is widespread in the prey species [6] . If we consider that one of the coronaviruses, such as COVID- 19 , is widespread in the prey species of the 4D ecosystem. Consequently, the mode of transmission from the infected prey to human might be occurred in two sce-narios, the consumption of infected prey as a source of food, and direct contact or contact through an intermediate host. However, the 4D ecosystem model is established by combining classical Lotka-Volterra model with the classic SIR model, which can be defined as follows: The above dynamical system is a predator-prey model which consists of a generalist predator subsists on two prey species in which a disease is becoming endemic in one of them. In this model, a, b, c, d, R are the positive constant parameters, and x 1 , x 2 , x 3 , x 4 are the state variables defined as follows: x 1 is the healthy populations of the susceptible prey species; x 2 is infected populations of the prey species; x 3 is the populations of the immune prey species; and x 4 is the population of predators. By setting the left-hand side of differential equations in the system (1) equal to zero, one can obtain five different equilibria as follows: Linearizing the system (1) at an arbitrary equilib- , the following Jacobian matrix is obtained: where i = 1, 2, 3, 4, 5. To determine the stability of each equilibrium point, one has to compute the corresponding eigenvalues. The following statements are true. Proof It is easy to obtain the eigenvalues of the system at E i by solving the following equation: (1) Therefore, the eigenvalues of the system at E 1 are given by Hence, the equilibrium point E 1 is unstable saddle. (2) The eigenvalues of the system at E 2 are given by It is well known that an equilibrium point is unstable only when there is at least one eigenvalue with positive real part. Thus, E 2 is unstable point if R − a > 1 and/or b > ac. (3) The eigenvalues of the system at E 3 are given by Since all the parameters of the system are greater then zero, there is only one eigenvalue can have positive real part, which is λ 1 . As a result, E 2 is unstable point only when a > b c . (4) The eigenvalues of the system at E 4 are given by Since b > 0, then E 4 is unstable equilibrium point. Finally, the characteristic equation of the system at E 5 has mixed signs. By Routh-Hurwitz criterion, the real parts of the corresponding eigenvalues have also positive and negative signs, which means that E 5 is unstable equilibrium point. As can be seen in the previous subsection, the 4D infected prey model (1) has five different equilibria in which there are always at least three unstable equilibria. Consequently, two important matters can be concluded: (1) it is quite possible to have chaos in this system; (2) in case of presence chaotic attractor, this attractor is called self-excited. Maximum Lyapunov Exponents (MLE) is one of the fundamental algorithms used to identify chaos in nonlinear dynamical systems. Generally, the MLE of a dynamics system is calculated versus one parameter varying. In this paper, we calculate the MLE of the system (1) with two parameters varying at the same time, which can provide a wider vision about the chaotic regions of the system. Therefore, four sets of parameters are chosen for simulating four cases of MLE with step size of 0.01, as shown in Fig. 1. Based on this figure and with the support of time series of each case, we observe that the system shows chaotic behavior only when MLE > 0.001. Otherwise, the system shows no chaos. However, it crucial to clarify that the chaotic regions of the 4D epidemic outbreaks system (1) in Fig. 1 are valid only when the initial conditions are chosen to be (1, 1, 1, 1). In the previous section, we have demonstrated that the system (1) shows chaotic attractors as well as nonchaotic attractors. These two different behaviors are appeared with only one specific set of initial conditions. In fact, initiating from a one set of initial values may not give a full picture about the dynamics of a chaotic system. Therefore, this section investigates the dynamics of the system (1) through choosing randomly several sets of initial conditions in which the chaotic and nonchaotic regions of the system are examined with these sets of initial conditions. When the parameters of the system (1) are selected as a = 0.3, b = 0.4, c = 2, d = 0.3, and R to be a control parameter for the range 1.26 ≤ R ≤ 1.32, the bifurcation diagram of the state variable x 1 and MLE are depicted in Fig. 2a , b, respectively. In this figure, we have chosen nine different sets of initial conditions to examine the non-chaotic and chaotic regions of the system. As can be seen in Fig. 2 , the system exhibits extreme multistability behaviors in which the coexistence of nine different quasi-periodic attractors is primarily occurred within the range 1.26 ≤ R ≤ 1.3. Meanwhile, the coexistence of several chaotic and quasiperiodic attractors is appeared for R ∈ (1.3, 1.32]. To further illustrate the case of coexisting quasi-periodic attractors, we depict the phase portraits and the time series of the system for R = 1.3, as shown in Fig. 3 . The coexistence of nine quasi-periodic attractors is quite clear in the time series of the state variable x 2 as well as two different projections. However, it is important to mention that choosing many other sets of initial conditions can produce infinitely many other coexisting quasi-periodic attractors. Additionally, Fig. 4a . The evolution of these orbits begin with coexisting quasi-periodic attractors and then spread over a larger region of the equilibrium points to produce coexistence of chaotic attractors, as illustrated in Fig. 5 . As can be seen in Fig. 5 , five different chaotic attractors are generated by the system via selecting a certain set of parameters and different sets of initial conditions. In fact, except for these coexisting chaotic attractors, there are infinitely many others, which can be generated by selecting appropriate sets of initial conditions. The most complicated behavior in nonlinear dynamical systems is the transit from chaotic to non-chaotic behaviors (and vice versa) for a certain sets of sys-tem parameters. This is due to inefficient techniques for analyzing the dynamical systems with different sets of initial conditions, especially when the systems exhibit infinitely many coexisting attractors even in a small region of these sets. To check whether the system (1) shows this complicated transitions, we depict Fig. 6a , b, respectively. The orbits in Fig. 6 colored in blue, green, and purple begin with the initial conditions (2.5, 0.2, 0.2, 0.2), (2, 0.2, 0.2, 0.2), and (0.5, 0.2, 0.2, 0.2), respectively. It can be seen that there are several transitions from chaotic to quasi-periodic attractors for all three orbits. For instance, the orbit of purple color drops from chaos to quasi-periodic when R = 1.799, as illustrated in Fig. 7a -c. Therefore, it can be observed that the behavior of system (1) is unpredictable in which choosing only three sets of initial conditions can generate coexistence between chaotic and non-chaotic attractors. Different chaotic systems have different degrees of sensitivity. Based on our simulation results, it can be concluded that the system (1) is very sensitive to its parameters and initial conditions, in which a tiny alteration in the initial values can drive the system to exhibit completely different behavior. In other words, each set of system parameters can produce infinity many coexisting attractors by choosing appropriate sets of initial conditions. As a result, the high-sensitivity performance of the system (1) provides a clear explanation about the unpredictability of the epidemic among infected preys. However, the next section would describe the complexity of the attractors that generate by the system (1). This section investigates the complexity performance of each variable x i (i = 1, 2, 3, 4) of system (1) using one of the fundamental algorithms, which is Sample Entropy (SamEn) [22] . Generally, the SamEn measure is used to estimate how much extra information is required to predict the (t + 1)th output of a trajectory using its previous (t) outputs. Obviously, larger positive SamEn values indicate that the trajectories need a more additional information to be predicted. In our simulations, the parameters of SamEn measure are set as the embedding dimension m = 2 and tolerance parameter r = 0.2×Standard Deviation. As aforementioned in the previous section, the system (1) shows extreme multistability behaviors, which means that tiny changes in the original prey species could enlarge and be widespread, and that could confirm through studying the complexity of the time series of the system's variables. Therefore, let us display the complexity of the variables x i (i = 1, 2, 3, 4) of the 1. For all values of the parameter R, the system with varying a single initial condition exhibits different complexity values, especially for the variables x 1 , x 2 , and x 4 , in which the complexity values of the first case increase with the increasing of the initial condition k 2 ; meanwhile, the complexity values of the second case decrease with the increasing of the initial condition k 3 ; 2. The variable of healthy prey (x 1 ) is the most complicated variable in the system; meanwhile, it is obvi-ous that the variable of infected prey (x 2 ) is more complex than the variable of immune prey (x 3 ). In conclusion, the complexity results of the system's variables confirm the high sensitivity to initial conditions; hence, the changes in the system's variables expand over time. Therefore, it is reasonable now to consider the endemic in the prey species of the system could evolve to be pandemic such as COVID-19. This paper investigated the 4D dynamical model of an ecosystem, which constructs by combining the classical Lotka-Volterra model with the classic SIR model. This model describes the interplay between a generalist predator subsists on two prey species in which disease is becoming endemic in one of them. Several numerical analysis methods find that this model shows chaotic behavior and extreme multistability behaviors. The occurrence of extreme multistability comes from the high sensitivity of the system. That means a tiny change in the system's variable could enlarge over time. This might interpret the widespread disease among prey species till becoming endemic. Since disease spillover from animals to humans could occur when the infec-tion is widespread in the prey species, it is sensible to consider that the endemic in the infected prey of the 4D model can jump to humans. However, the complexity of the disease in infected prey species can be evaluated by the Sample entropy algorithm. Contour plot results have demonstrated the unpredictability of disease in the prey species. The results of this paper could provide a demonstration of the randomness of the COVID-19 outbreak among prey species and humans. Outbreak of pneumonia of unknown etiology in Wuhan, China: the mystery and the miracle Recent progress and challenges in drug development against COVID-19 coronavirus (SARS-CoV-2)-an update on the status Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges Can India develop herd immunity against COVID-19? The second and third waves in India: when will the pandemic be culminated? Origins of major human infectious diseases COVID-19 infection: origin, transmission, and characteristics of human coronaviruses A model based study on the dynamics of COVID-19: prediction and control Chaos in disease outbreaks among prey The uneasy coexistence of predators and pathogens A predator-prey model with disease in the prey A predatorprey model with infected prey. Theor Effects of a disease affecting a predator on the dynamics of a predator-prey system Simple mathematical models with very complicated dynamics. Theory Chaotic Attractors Degenerating the butterfly attractor in a plasma perturbation model using nonlinear controllers Dynamics of a new hyperchaotic system and multistability Multistability and chaotic scenario in a quantum pair-ion plasma Multistability and chaos in a noise-induced blood flow Multistability and dynamical properties of quantum ion-acoustic flow Self-excited and hidden attractors in a novel chaotic system with complicated multistability Can hyperchaotic maps with high complexity produce multistability? Physiological time-series analysis using approximate entropy and sample entropy Acknowledgements Asit Saha is thankful to SMIT (SMU) for providing research funding (Ref. Nos. 6100/ SMIT/ R&D/ Project/ 05/ 2018). All authors have contributed equally to the paper. The authors declare that they have no conflict of interest.