Impact of periodic vaccination in SEIRS seasonal model.

We study three different strategies of vaccination in an SEIRS (Susceptible-Exposed-Infected-Recovered-Susceptible) seasonal forced model, which are (i) continuous vaccination; (ii) periodic short-time localized vaccination, and (iii) periodic pulsed width campaign. Considering the first strategy, we obtain an expression for the basic reproduction number and infer a minimum vaccination rate necessary to ensure the stability of the disease-free equilibrium (DFE) solution. In the second strategy, short duration pulses are added to a constant baseline vaccination rate. The pulse is applied according to the seasonal forcing phases. The best outcome is obtained by locating intensive immunization at inflection of the transmissivity curve. Therefore, a vaccination rate of 44.4% of susceptible individuals is enough to ensure DFE. For the third vaccination proposal, additionally to the amplitude, the pulses have a prolonged time width. We obtain a non-linear relationship between vaccination rates and the duration of the campaign. Our simulations show that the baseline rates, as well as the pulse duration, can substantially improve the vaccination campaign effectiveness. These findings are in agreement with our analytical expression. We show a relationship between the vaccination parameters and the accumulated number of infected individuals, over the years, and show the relevance of the immunization campaign annual reaching for controlling the infection spreading. Regarding the dynamical behavior of the model, our simulations show that chaotic and periodic solutions as well as bi-stable regions depend on the vaccination parameters range.

Multistability and chaos in SEIRS epidemic model with a periodic time-dependent transmission rate.

In this work, we study the dynamics of a susceptible-exposed-infectious-recovered-susceptible epidemic model with a periodic time-dependent transmission rate. Emphasizing the influence of the seasonality frequency on the system dynamics, we analyze the largest Lyapunov exponent along parameter planes finding large chaotic regions. Furthermore, in some ranges, there are shrimp-like periodic structures. We highlight the system multistability, identifying the coexistence of periodic orbits for the same parameter values, with the infections maximum distinguishing by up one order of magnitude, depending only on the initial conditions. In this case, the basins of attraction have self-similarity. Parametric configurations, for which both periodic and non-periodic orbits occur, cover 13.20% of the evaluated range. We also identified the coexistence of periodic and chaotic attractors with different maxima of infectious cases, where the periodic scenario peak reaches approximately 50% higher than the chaotic one.

Nontwist field line mapping in a tokamak with ergodic magnetic limiter.

For tokamaks with uniform magnetic shear, Martin and Taylor have proposed a symplectic map which has been used to describe the magnetic field lines at the plasma edge perturbed by an ergodic magnetic limiter. We propose an analytical magnetic field line map, based on the Martin-Taylor map, for a tokamak with arbitrary safety factor profile. With the inclusion of a nonmonotonic profile, we obtain a nontwist map which presents the characteristic properties of degenerate systems, such as the twin islands scenario, shearless curve, and separatrix reconnection. We estimate the width of the islands and describe their changes of shape for large values of the limiter current. From our numerical simulations about the shearless curve, we show that its position and aspect depend on the control parameters.

Basin Entropy and Shearless Barrier Breakup in Open Non-Twist Hamiltonian Systems.

We consider open non-twist Hamiltonian systems represented by an area-preserving two-dimensional map describing incompressible planar flows in the reference frame of a propagating wave, and possessing exits through which map orbits can escape. The corresponding escape basins have a fractal nature that can be revealed by the so-called basin entropy, a novel concept developed to quantify final-state uncertainty in dynamical systems. Since the map considered violates locally the twist condition, there is a shearless barrier that prevents global chaotic transport. In this paper, we show that it is possible to determine the shearless barrier breakup by considering the variation in the escape basin entropy with a tunable parameter.

Stickiness and recurrence plots: An entropy-based approach.

The stickiness effect is a fundamental feature of quasi-integrable Hamiltonian systems. We propose the use of an entropy-based measure of the recurrence plots (RPs), namely, the entropy of the distribution of the recurrence times (estimated from the RP), to characterize the dynamics of a typical quasi-integrable Hamiltonian system with coexisting regular and chaotic regions. We show that the recurrence time entropy (RTE) is positively correlated to the largest Lyapunov exponent, with a high correlation coefficient. We obtain a multi-modal distribution of the finite-time RTE and find that each mode corresponds to the motion around islands of different hierarchical levels.

Suppression of chaotic bursting synchronization in clustered scale-free networks by an external feedback signal.

Oscillatory activities in the brain, detected by electroencephalograms, have identified synchronization patterns. These synchronized activities in neurons are related to cognitive processes. Additionally, experimental research studies on neuronal rhythms have shown synchronous oscillations in brain disorders. Mathematical modeling of networks has been used to mimic these neuronal synchronizations. Actually, networks with scale-free properties were identified in some regions of the cortex. In this work, to investigate these brain synchronizations, we focus on neuronal synchronization in a network with coupled scale-free networks. The networks are connected according to a topological organization in the structural cortical regions of the human brain. The neuronal dynamic is given by the Rulkov model, which is a two-dimensional iterated map. The Rulkov neuron can generate quiescence, tonic spiking, and bursting. Depending on the parameters, we identify synchronous behavior among the neurons in the clustered networks. In this work, we aim to suppress the neuronal burst synchronization by the application of an external perturbation as a function of the mean-field of membrane potential. We found that the method we used to suppress synchronization presents better results when compared to the time-delayed feedback method when applied to the same model of the neuronal network.

Curry-Yorke route to shearless attractors and coexistence of attractors in dissipative nontwist systems.

The routes to chaos play an important role in predictions about the transitions from regular to irregular behavior in nonlinear dynamical systems, such as electrical oscillators, chemical reactions, biomedical rhythms, and nonlinear wave coupling. Of special interest are dissipative systems obtained by adding a dissipation term in a given Hamiltonian system. If the latter satisfies the so-called twist property, the corresponding dissipative version can be called a "dissipative twist system." Transitions to chaos in these systems are well established; for instance, the Curry-Yorke route describes the transition from a quasiperiodic attractor on torus to chaos passing by a chaotic banded attractor. In this paper, we study the transitions from an attractor on torus to chaotic motion in dissipative nontwist systems. We choose the dissipative standard nontwist map, which is a non-conservative version of the standard nontwist map. In our simulations, we observe the same transition to chaos that happens in twist systems, known as a soft one, where the quasiperiodic attractor becomes wrinkled and then chaotic through the Curry-Yorke route. By the Lyapunov exponent, we study the nature of the orbits for a different set of parameters, and we observe that quasiperiodic motion and periodic and chaotic behavior are possible in the system. We observe that they can coexist in the phase space, implying in multistability. The different coexistence scenarios were studied by the basin entropy and by the boundary basin entropy.

Ratchet current in nontwist Hamiltonian systems.

Non-monotonic area-preserving maps violate the twist condition locally in phase space, giving rise to shearless invariant barriers surrounded by twin island chains in these regions of phase space. For the extended standard nontwist map, with two resonant perturbations with distinct wave numbers, we investigate the presence of such barriers and their associated island chains and compare our results with those that have been reported for the standard nontwist map with only one perturbation. Furthermore, we determine in the control parameter space the existence of the shearless barrier and the influence of the additional wave number on this condition. We show that only for odd second wave numbers are the twin island chains symmetrical. Moreover, for even wave numbers, the lack of symmetry between the chains of twin islands generates a ratchet effect that implies a directed transport in the phase space.

Basin of attraction for chimera states in a network of Rössler oscillators.

Chimera states are spatiotemporal patterns in which coherent and incoherent dynamics coexist simultaneously. These patterns were observed in both locally and nonlocally coupled oscillators. We study the existence of chimera states in networks of coupled Rössler oscillators. The Rössler oscillator can exhibit periodic or chaotic behavior depending on the control parameters. In this work, we show that the existence of coherent, incoherent, and chimera states depends not only on the coupling strength, but also on the initial state of the network. The initial states can belong to complex basins of attraction that are not homogeneously distributed. Due to this fact, we characterize the basins by means of the uncertainty exponent and basin stability. In our simulations, we find basin boundaries with smooth, fractal, and riddled structures.

Spike-burst chimera states in an adaptive exponential integrate-and-fire neuronal network.

Chimera states are spatiotemporal patterns in which coherence and incoherence coexist. We observe the coexistence of synchronous (coherent) and desynchronous (incoherent) domains in a neuronal network. The network is composed of coupled adaptive exponential integrate-and-fire neurons that are connected by means of chemical synapses. In our neuronal network, the chimera states exhibit spatial structures both with spike and burst activities. Furthermore, those desynchronized domains not only have either spike or burst activity, but we show that the structures switch between spikes and bursts as the time evolves. Moreover, we verify the existence of multicluster chimera states.

Using rotation number to detect sticky orbits in Hamiltonian systems.

In Hamiltonian systems, depending on the control parameter, orbits can stay for very long times around islands, the so-called stickiness effect caused by a temporary trapping mechanism. Different methods have been used to identify sticky orbits, such as recurrence analysis, recurrence time statistics, and finite-time Lyapunov exponent. However, these methods require a large number of map iterations and to know the island positions in the phase space. Here, we show how to use the small divergence of bursts in the rotation number calculation as a tool to identify stickiness without knowing the island positions. This new procedure is applied to the standard map, a map that has been used to describe the dynamic behavior of several nonlinear systems. Moreover, our procedure uses a small number of map iterations and is proper to identify the presence of stickiness phenomenon for different values of the control parameter.

Efficient manifolds tracing for planar maps.

In this work, we introduce an exact calculation method and an approximation technique for tracing the invariant manifolds of unstable periodic orbits of planar maps. The exact method relies on an adaptive refinement procedure that prevents redundant calculations occurring in other approaches, and the approximated method relies on a novel interpolation approach based on normal displacement functions. The resulting approximated manifold is precise when compared to the exact one, and its relative computational cost falls like the inverse of the manifold length. To present the tracing method, we obtain the invariant manifolds of the Chirikov-Taylor map, and as an application we illustrate the transition from homoclinic to heteroclinic chaos in the Duffing oscillator that leads from localized chaos to global chaotic motion.

Recurrence-based analysis of barrier breakup in the standard nontwist map.

We study the standard nontwist map that describes the dynamic behaviour of magnetic field lines near a local minimum or maximum of frequency. The standard nontwist map has a shearless invariant curve that acts like a barrier in phase space. Critical parameters for the breakup of the shearless curve have been determined by procedures based on the indicator points and bifurcations of periodical orbits, a methodology that demands high computational cost. To determine the breakup critical parameters, we propose a new simpler and general procedure based on the determinism analysis performed on the recurrence plot of orbits near the critical transition. We also show that the coexistence of islands and chaotic sea in phase space can be analysed by using the recurrence plot. In particular, the measurement of determinism from the recurrence plot provides us with a simple procedure to distinguish periodic from chaotic structures in the parameter space. We identify an invariant shearless breakup scenario, and we also show that recurrence plots are useful tools to determine the presence of periodic orbit collisions and bifurcation curves.

The dose-dense principle in chemotherapy.

Chemotherapy is a cancer treatment modality that uses drugs to kill tumor cells. A typical chemotherapeutic protocol consists of several drugs delivered in cycles of three weeks. We present mathematical analyses demonstrating the existence of a maximum time between cycles of chemotherapy for a protocol to be effective. A mathematical equation is derived, which relates such a maximum time with the variables that govern the kinetics of the tumor and those characterizing the chemotherapeutic treatment. Our results suggest that there are compelling arguments supporting the use of dose-dense protocols. Finally, we discuss the limitations of these protocols and suggest an alternative.

Dynamical properties of the soft-wall elliptical billiard.

Physical systems such as optical traps and microwave cavities are realistically modeled by billiards with soft walls. In order to investigate the influence of the wall softness on the billiard dynamics, we study numerically a smooth two-dimensional potential well that has the elliptical (hard-wall) billiard as a limiting case. Considering two parameters, the eccentricity of the elliptical equipotential curves and the wall hardness, which defines the steepness of the well, we show that (1) whereas the hard-wall limit is integrable and thus completely regular, the soft wall elliptical billiard exhibits chaos, (2) the chaotic fraction of the phase space depends nonmonotonically on the hardness of the wall, and (3) the effect of the hardness on the dynamics depends strongly on the eccentricity of the billiard. We further show that the limaçon billiard can exhibit enhanced chaos induced by wall softness, which suggests that our findings generalize to quasi-integrable systems.

Suppression of phase synchronisation in network based on cat's brain.

We have studied the effects of perturbations on the cat's cerebral cortex. According to the literature, this cortex structure can be described by a clustered network. This way, we construct a clustered network with the same number of areas as in the cat matrix, where each area is described as a sub-network with a small-world property. We focus on the suppression of neuronal phase synchronisation considering different kinds of perturbations. Among the various controlling interventions, we choose three methods: delayed feedback control, external time-periodic driving, and activation of selected neurons. We simulate these interventions to provide a procedure to suppress undesired and pathological abnormal rhythms that can be associated with many forms of synchronisation. In our simulations, we have verified that the efficiency of synchronisation suppression by delayed feedback control is higher than external time-periodic driving and activation of selected neurons of the cat's cerebral cortex with the same coupling strengths.

Mathematical model of brain tumour with glia-neuron interactions and chemotherapy treatment.

In recent years, it became clear that a better understanding of the interactions among the main elements involved in the cancer network is necessary for the treatment of cancer and the suppression of cancer growth. In this work we propose a system of coupled differential equations that model brain tumour under treatment by chemotherapy, which considers interactions among the glial cells, the glioma, the neurons, and the chemotherapeutic agents. We study the conditions for the glioma growth to be eliminated, and identify values of the parameters for which the inhibition of the glioma growth is obtained with a minimal loss of healthy cells.

Intermingled basins in coupled Lorenz systems.

We consider a system of two identical linearly coupled Lorenz oscillators presenting synchronization of chaotic motion for a specified range of the coupling strength. We verify the existence of global synchronization and antisynchronization attractors with intermingled basins of attraction such that the basin of one attractor is riddled with holes belonging to the basin of the other attractor and vice versa. We investigated this phenomenon by verifying the fulfillment of the mathematical requirements for intermingled basins and obtained scaling laws that characterize quantitatively the riddling of both basins in this system.

Local predictability and nonhyperbolicity through finite Lyapunov exponent distributions in two-degrees-of-freedom Hamiltonian systems.

By using finite Lyapunov exponent distributions, we get insight into both the local and global properties of a dynamical flow, including its nonhyperbolic behavior. Several distributions of finite Lyapunov exponents have been computed in two prototypical four-dimensional phase-space Hamiltonian systems. They have been computed calculating the growth rates of a set of orthogonal axes arbitrarily pointed at given intervals. We analyze how such distributions serve or not for tracing the orbit nature and local flow properties such as the unstable dimension variability, as the axes are allowed or not to tend to the largest stretching direction. The relationship between the largest and closest to zero exponent distribution is analyzed. It shows a linear dependency at short intervals, related to the number of degrees of freedom of the system. Finally, the hyperbolicity indexes, associated to the shadowing times, are calculated. They provide interesting information at very local scales, even when there are no Gaussian distributions and the values cannot be regarded as random variables.

Using recurrences to characterize the hyperchaos-chaos transition.

Hyperchaos occurs in a dynamical system with more than one positive Lyapunov exponent. When the equations governing the time evolution of the dynamical system are known, the transition from chaos to hyperchaos can be readily obtained when the second largest Lyapunov exponent crosses zero. If the only information available on the system is a time series, however, such method is difficult to apply. We propose the use of recurrence quantification analysis of a time series to characterize the chaos-hyperchaos transition. We present results obtained from recurrence plots of coupled chaotic piecewise-linear maps and Chua-Matsumoto circuits, but the method can be applied as well to other systems, even when one does not know their dynamical equations.