Loading...
7 results
Search Results
Now showing 1 - 7 of 7
- Scaling law in Saddle-node bifurcations for one-dimensional maps: a complex variable approachPublication . Duarte, Jorge; Januário, Cristina; Martins, Nuno; Sardanyés, JosepThe study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, tau, scales according to an inverse square-root power law, tau similar to (mu-mu (c) )(-1/2), as the bifurcation parameter mu, is driven further away from its critical value, mu (c) . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.
- Complex Dynamics of defective interfering baculoviruses during serial passage in insect cellsPublication . Zwart, Mark P.; Pijlman, Gorben P.; Sardanyés, Josep; Duarte, Jorge; Januário, Cristina; Elena, Santiago F.Defective interfering (DI) viruses are thought to cause oscillations in virus levels, known as the ‘Von Magnus effect’. Interference by DI viruses has been proposed to underlie these dynamics, although experimental tests of this idea have not been forthcoming. For the baculoviruses, insect viruses commonly used for the expression of heterologous proteins in insect cells, the molecular mechanisms underlying DI generation have been investigated. However, the dynamics of baculovirus populations harboring DIs have not been studied in detail. In order to address this issue, we used quantitative real-time PCR to determine the levels of helper and DI viruses during 50 serial passages of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) in Sf21 cells. Unexpectedly, the helper and DI viruses changed levels largely in phase, and oscillations were highly irregular, suggesting the presence of chaos. We therefore developed a simple mathematical model of baculovirus-DI dynamics. This theoretical model reproduced patterns qualitatively similar to the experimental data. Although we cannot exclude that experimental variation (noise) plays an important role in generating the observed patterns, the presence of chaos in the model dynamics was confirmed with the computation of the maximal Lyapunov exponent, and a Ruelle-Takens-Newhouse route to chaos was identified at decreasing production of DI viruses, using mutation as a control parameter. Our results contribute to a better understanding of the dynamics of DI baculoviruses, and suggest that changes in virus levels over passages may exhibit chaos.
- On chaos, transient chaos and ghosts in single populations models with allee effectsPublication . Duarte, Jorge; Januário, Cristina; Martins, Nuno; Sardanyés, JosepDensity-dependent effects, both positive or negative, can have an important impact on the population dynamics of species by modifying their population per-capita growth rates. An important type of such density-dependent factors is given by the so-called Allee effects, widely studied in theoretical and field population biology. In this study, we analyze two discrete single population models with overcompensating density-dependence and Allee effects due to predator saturation and mating limitation using symbolic dynamics theory. We focus on the scenarios of persistence and bistability, in which the species dynamics can be chaotic. For the chaotic regimes, we compute the topological entropy as well as the Lyapunov exponent under ecological key parameters and different initial conditions. We also provide co-dimension two bifurcation diagrams for both systems computing the periods of the orbits, also characterizing the period-ordering routes toward the boundary crisis responsible for species extinction via transient chaos. Our results show that the topological entropy increases as we approach to the parametric regions involving transient chaos, being maximum when the full shift R(L)(infinity) occurs, and the system enters into the essential extinction regime. Finally, we characterize analytically, using a complex variable approach, and numerically the inverse square-root scaling law arising in the vicinity of a saddle-node bifurcation responsible for the extinction scenario in the two studied models. The results are discussed in the context of species fragility under differential Allee effects. (C) 2011 Elsevier Ltd. All rights reserved.
- Activation of effector immune cells promotes tumor stochastic extinction: A homotopy analysis approachPublication . Sardanyés, Josep; Rodrigues, Carla; Januário, Cristina; Martins, Nuno; Gil-Gómez, Gabriel; Duarte, JorgeIn this article we provide homotopy solutions of a cancer nonlinear model describing the dynamics of tumor cells in interaction with healthy and effector immune cells. We apply a semi-analytic technique for solving strongly nonlinear systems – the Step Homotopy Analysis Method (SHAM). This algorithm, based on a modification of the standard homotopy analysis method (HAM), allows to obtain a one-parameter family of explicit series solutions. By using the homotopy solutions, we first investigate the dynamical effect of the activation of the effector immune cells in the deterministic dynamics, showing that an increased activation makes the system to enter into chaotic dynamics via a period-doubling bifurcation scenario. Then, by adding demographic stochasticity into the homotopy solutions, we show, as a difference from the deterministic dynamics, that an increased activation of the immune cells facilitates cancer clearance involving tumor cells extinction and healthy cells persistence. Our results highlight the importance of therapies activating the effector immune cells at early stages of cancer progression.
- How complex, probable, and predictable is genetically driven red queen chaos?Publication . Duarte, Jorge; Rodrigues, Carla; Januário, Cristina; Martins, Nuno; Sardanyés, JosepCoevolution between two antagonistic species has been widely studied theoretically for both ecologically- and genetically-driven Red Queen dynamics. A typical outcome of these systems is an oscillatory behavior causing an endless series of one species adaptation and others counter-adaptation. More recently, a mathematical model combining a three-species food chain system with an adaptive dynamics approach revealed genetically driven chaotic Red Queen coevolution. In the present article, we analyze this mathematical model mainly focusing on the impact of species rates of evolution (mutation rates) in the dynamics. Firstly, we analytically proof the boundedness of the trajectories of the chaotic attractor. The complexity of the coupling between the dynamical variables is quantified using observability indices. By using symbolic dynamics theory, we quantify the complexity of genetically driven Red Queen chaos computing the topological entropy of existing one-dimensional iterated maps using Markov partitions. Co-dimensional two bifurcation diagrams are also built from the period ordering of the orbits of the maps. Then, we study the predictability of the Red Queen chaos, found in narrow regions of mutation rates. To extend the previous analyses, we also computed the likeliness of finding chaos in a given region of the parameter space varying other model parameters simultaneously. Such analyses allowed us to compute a mean predictability measure for the system in the explored region of the parameter space. We found that genetically driven Red Queen chaos, although being restricted to small regions of the analyzed parameter space, might be highly unpredictable.
- Controling delayed transitions with applications to prevent single species extinctionsPublication . Duarte, Jorge; Januário, Cristina; Martins, Nuno; Sardanyés, JosepA new method is proposed to control delayed transitions towards extinction in single population theoretical models with discrete time undergoing saddle-node bifurcations. The control method takes advantage of the delaying properties of the saddle remnant arising after the bifurcation, and allows to sustain populations indefinitely. Our method, which is shown to work for deterministic and stochastic systems, could generally be applied to avoid transitions tied to one-dimensional maps after saddle-node bifurcations.
- Optimal homotopy analysis of a chaotic HIV-1 model incorporating AIDS-related cancer cellsPublication . Duarte, Jorge; Januário, Cristina; Martins, Nuno; Ramos, Carlos; Rodrigues, Carla; Sardanyés, JosepThe studies of nonlinear models in epidemiology have generated a deep interest in gaining insight into the mechanisms that underlie AIDS-related cancers, providing us with a better understanding of cancer immunity and viral oncogenesis. In this article, we analyze an HIV-1 model incorporating the relations between three dynamical variables: cancer cells, healthy CD4 + T lymphocytes, and infected CD4 + T lymphocytes. Recent theoretical investigations indicate that these cells interactions lead to different dynamical outcomes, for instance to periodic or chaotic behavior. Firstly, we analytically prove the boundedness of the trajectories in the system’s attractor. The complexity of the coupling between the dynamical variables is quantified using observability indices. Our calculations reveal that the highest observable variable is the population of cancer cells, thus indicating that these cells could be monitored in future experiments in order to obtain time series for attractor’s reconstruction. We identify different dynamical behaviors of the system varying two biologically meaningful parameters: r 1, representing the uncontrolled proliferation rate of cancer cells, and k 1, denoting the immune system’s killing rate of cancer cells. The maximum Lyapunov exponent is computed to identify the chaotic regimes. Considering very recent developments in the literature related to the homotopy analysis method (HAM), we calculate the explicit series solutions of the cancer model and focus our analysis on the dynamical variable with the highest observability index. An optimal homotopy analysis approach is used to improve the computational efficiency of HAM by means of appropriate values for the convergence control parameter, which greatly accelerate the convergence of the series solution. The approximated analytical solutions are used to compute density plots, which allow us to discuss additional dynamical features of the model.