Showing 1–50 of 111 results for author: Caldas, I L

Identifying ballistic modes via Poincaré sections

Authors: André Farinha Bósio, Iberê Luiz Caldas, Ricardo Luiz Viana, Yves Elskens

Abstract: Exploring chaotic systems via Poincaré sections has proven essential in dynamical systems, yet measuring their characteristics poses challenges to identify the various dynamical regimes considered. In this paper, we propose a new approach that uses image processing to distinguish chaotic and regular regions of area-preserving dynamics, and then classify the transport regime. We characterize differ… ▽ More Exploring chaotic systems via Poincaré sections has proven essential in dynamical systems, yet measuring their characteristics poses challenges to identify the various dynamical regimes considered. In this paper, we propose a new approach that uses image processing to distinguish chaotic and regular regions of area-preserving dynamics, and then classify the transport regime. We characterize different transport regimes in the standard map with the proposed method based on image reconstruction techniques, identifying the superdiffusion much faster than the usual mean square displacement method. The procedure is also applied to a two-wave, time-dependent Hamiltonian to investigate superdiffusion in function of two parameters. △ Less

Submitted 19 September, 2024; originally announced September 2024.

Control, bi-stability and preference for chaos in time-dependent vaccination campaign

Authors: Enrique C. Gabrick, Eduardo L. Brugnago, Ana L. R. de Moraes, Paulo R. Protachevicz, Sidney T. da Silva, Fernando S. Borges, Iberê L. Caldas, Antonio M. Batista, Jürgen Kurths

Abstract: In this work, effects of constant and time-dependent vaccination rates on the Susceptible-Exposed-Infected-Recovered-Susceptible (SEIRS) seasonal model are studied. Computing the Lyapunov exponent, we show that typical complex structures, such as shrimps, emerge for given combinations of constant vaccination rate and another model parameter. In some specific cases, the constant vaccination does no… ▽ More In this work, effects of constant and time-dependent vaccination rates on the Susceptible-Exposed-Infected-Recovered-Susceptible (SEIRS) seasonal model are studied. Computing the Lyapunov exponent, we show that typical complex structures, such as shrimps, emerge for given combinations of constant vaccination rate and another model parameter. In some specific cases, the constant vaccination does not act as a chaotic suppressor and chaotic bands can exist for high levels of vaccination (e.g., $> 0.95$). Moreover, we obtain linear and non-linear relationships between one control parameter and constant vaccination to establish a disease-free solution. We also verify that the total infected number does not change whether the dynamics is chaotic or periodic. The introduction of a time-dependent vaccine is made by the inclusion of a periodic function with a defined amplitude and frequency. For this case, we investigate the effects of different amplitudes and frequencies on chaotic attractors, yielding low, medium, and high seasonality degrees of contacts. Depending on the parameters of the time-dependent vaccination function, chaotic structures can be controlled and become periodic structures. For a given set of parameters, these structures are accessed mostly via crisis and in some cases via period-doubling. After that, we investigate how the time-dependent vaccine acts in bi-stable dynamics when chaotic and periodic attractors coexist. We identify that this kind of vaccination acts as a control by destroying almost all the periodic basins. We explain this by the fact that chaotic attractors exhibit more desirable characteristics for epidemics than periodic ones in a bi-stable state. △ Less

Submitted 29 August, 2024; originally announced September 2024.

Effective transport barriers in the biquadratic nontwist map

Authors: Gabriel C. Grime, Ricardo L. Viana, Yves Elskens Iberê L. Caldas

Abstract: Nontwist area-preserving maps violate the twist condition at specific orbits, resulting in shearless invariant curves that prevent chaotic transport. Plasmas and fluids with nonmonotonic equilibrium profiles may be described using nontwist systems, where even after these shearless curves breakdown, effective transport barriers persist, partially reducing transport coefficients. Some nontwist syste… ▽ More Nontwist area-preserving maps violate the twist condition at specific orbits, resulting in shearless invariant curves that prevent chaotic transport. Plasmas and fluids with nonmonotonic equilibrium profiles may be described using nontwist systems, where even after these shearless curves breakdown, effective transport barriers persist, partially reducing transport coefficients. Some nontwist systems present multiple shearless curves in phase space, increasing the complexity of transport phenomena, which have not been thoroughly investigated until now. In this work, we examine the formation of effective transport barriers in a nontwist area-preserving mapping with multiple shearless transport barriers. By quantifying the effectiveness of each transport barrier in phase space, we identified two scenarios where particular barriers dominate over others. Our results also reveal configurations where the interplay of two transport barriers creates regions in phase space with significant orbit trapping, thereby influencing the overall transport dynamics. △ Less

Submitted 1 September, 2024; originally announced September 2024.

Shearless bifurcations for two isochronous resonant perturbations

Authors: Bruno B. Leal, Matheus J. Lazarotto, Michele Mugnaine, Alfredo M. Ozorio de Almeida, Ricardo L. Viana, Iberê L. Caldas

Abstract: In nontwist systems, primary shearless curves act as barriers to chaotic transport. Surprisingly, the onset of secondary shearless curves has been reported in a few twist systems. Meanwhile, we found that, in twist systems, the onset of these secondary shearless curves is a standard process that may appear as control parameters are varied in situations where there is resonant mode coupling. Namely… ▽ More In nontwist systems, primary shearless curves act as barriers to chaotic transport. Surprisingly, the onset of secondary shearless curves has been reported in a few twist systems. Meanwhile, we found that, in twist systems, the onset of these secondary shearless curves is a standard process that may appear as control parameters are varied in situations where there is resonant mode coupling. Namely, we analyze these shearless bifurcations in two-harmonic systems for the standard map, the Ullmann map, and for the Walker-Ford Hamiltonian flow. The onset of shearless curves is related to bifurcations of periodic points. Furthermore, depending on the bifurcation, these shearless curves can emerge alone or in pairs, and in some cases, deform into separatrices. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Shrinking shrimp-shaped domains and multistability in the dissipative asymmetric kicked rotor map

Authors: Matheus Rolim Sales, Michele Mugnaine, Edson Denis Leonel, Iberê L. Caldas, José Danilo Szezech Jr

Abstract: An interesting feature in dissipative nonlinear systems is the emergence of characteristic domains in parameter space that exhibit periodic temporal evolution, known as shrimp-shaped domains. We investigate the parameter space of the dissipative asymmetric kicked rotor map and show that, in the regime of strong dissipation, the shrimp-shaped domains repeat themselves as the nonlinearity parameter… ▽ More An interesting feature in dissipative nonlinear systems is the emergence of characteristic domains in parameter space that exhibit periodic temporal evolution, known as shrimp-shaped domains. We investigate the parameter space of the dissipative asymmetric kicked rotor map and show that, in the regime of strong dissipation, the shrimp-shaped domains repeat themselves as the nonlinearity parameter increases while maintaining the same period. We analyze the dependence of the length of each periodic domain with the nonlinearity parameter, revealing that it follows a power law with the same exponent regardless of the dissipation parameter. Additionally, we find that the distance between adjacent shrimp-shaped domains is scaling invariant with respect to the dissipation parameter. Furthermore, we show that for weaker dissipation, a multistable scenario emerges within the periodic domains. We find that as the dissipation gets weaker, the ratio of multistable parameters for each periodic domain increases, and the area of the periodic basin decreases as the nonlinearity parameter increases. △ Less

Submitted 13 August, 2024; originally announced August 2024.

Larmor radius effect on the control of chaotic transport in tokamaks

Authors: L. A. Osorio-Quiroga, M. Roberto, R. L. Viana, Y. Elskens, I. L. Caldas

Abstract: We introduce a second-order gyro-averaged extension to a guiding-center transport model for test particles subject to an $\mathbf{E}\times\mathbf{B}$ drift in a tokamak. Using this approximation, we numerically investigate the influence of finite Larmor radius effects on the chaotic transport and the formation of transport barriers. For non-monotonic plasma profiles, we show that the twist conditi… ▽ More We introduce a second-order gyro-averaged extension to a guiding-center transport model for test particles subject to an $\mathbf{E}\times\mathbf{B}$ drift in a tokamak. Using this approximation, we numerically investigate the influence of finite Larmor radius effects on the chaotic transport and the formation of transport barriers. For non-monotonic plasma profiles, we show that the twist condition of the system is violated along a curve, usually called shearless transport barrier, which, under non-equilibrium conditions, exhibits significant resilience to destruction, thereby inhibiting chaotic transport. We analyze bifurcation diagrams of the shearless barriers and escape rates to explore the impact of the Larmor radius on controlling chaotic transport while varying the amplitude of electrostatic perturbations. Our findings show that increasing the Larmor radius enhances the robustness of transport barriers, as a larger electrostatic perturbation amplitude is required to disrupt them. Additionally, even in the absence of transport barriers, we observe a decrease in chaotic transport with an increase in the Larmor radius. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: Submitted to Physics of Plasmas

Shearless effective barriers to chaotic transport induced by even twin islands in nontwist systems

Authors: M. Mugnaine, J. D. Szezech Jr., R. L. Viana, I. L. Caldas, P. J. Morrison

Abstract: For several decades now it has been known that systems with shearless invariant tori, nontwist Hamiltonian systems, possess barriers to chaotic transport. These barriers are resilient to breakage under perturbation and therefore regions where they occur are natural places to look for barriers to transport. We describe a novel kind of effective barrier that persists after the shearless torus is bro… ▽ More For several decades now it has been known that systems with shearless invariant tori, nontwist Hamiltonian systems, possess barriers to chaotic transport. These barriers are resilient to breakage under perturbation and therefore regions where they occur are natural places to look for barriers to transport. We describe a novel kind of effective barrier that persists after the shearless torus is broken. Because phenomena are generic, for convenience we study the Standard Nontwist Map (SNM), an area-preserving map that violates the twist condition locally in the phase space. The novel barrier occurs in nontwist systems when twin even period islands are present, which happens for a broad range of parameter values in the SNM. With a phase space composed of regular and irregular orbits, the movement of chaotic trajectories is hampered by the existence of both shearless curves, total barriers, and a network of partial barriers formed by the stable and unstable manifolds of the hyperbolic points. Being a degenerate system, the SNM has twin islands and, consequently, twin hyperbolic points. We show that the structures formed by the manifolds intrinsically depend on period parity of the twin islands. For this even scenario the novel structure, named a torus free barrier, occurs because the manifolds of different hyperbolic points form an intricate chain atop a dipole configuration and the transport of chaotic trajectories through the chain becomes a rare event. This structure impacts the emergence of transport, the escape basin for chaotic trajectories, the transport mechanism and the chaotic saddle. The case of odd periodic orbits is different: we find for this case the emergence of transport immediately after the breakup of the last invariant curve, and this leads to a scenario of higher transport, with intricate escape basin boundary and a chaotic saddle with non-uniformly distributed points. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Ratchet current and scaling properties in a nontwist mapping

Authors: Matheus Rolim Sales, Daniel Borin, Leonardo Costa de Souza, José Danilo Szezech Jr., Ricardo Luiz Viana, Iberê Luiz Caldas, Edson Denis Leonel

Abstract: We investigate the transport of particles in the chaotic component of phase space for a two-dimensional, area-preserving nontwist map. The survival probability for particles within the chaotic sea is described by an exponential decay for regions in phase space predominantly chaotic and it is scaling invariant in this case. Alternatively, when considering mixed chaotic and regular regions, there is… ▽ More We investigate the transport of particles in the chaotic component of phase space for a two-dimensional, area-preserving nontwist map. The survival probability for particles within the chaotic sea is described by an exponential decay for regions in phase space predominantly chaotic and it is scaling invariant in this case. Alternatively, when considering mixed chaotic and regular regions, there is a deviation from the exponential decay, characterized by a power law tail for long times, a signature of the stickiness effect. Furthermore, due to the asymmetry of the chaotic component of phase space with respect to the line $I = 0$, there is an unbalanced stickiness which generates a ratchet current in phase space. Finally, we perform a phenomenological description of the diffusion of chaotic particles by identifying three scaling hypotheses, and obtaining the critical exponents via extensive numerical simulations. △ Less

Submitted 13 August, 2024; v1 submitted 10 June, 2024; originally announced June 2024.

When climate variables improve the dengue forecasting: a machine learning approach

Authors: Sidney T. da Silva, Enrique C. Gabrick, Paulo R. Protachevicz, Kelly C. Iarosz, Iberê L. Caldas, Antonio M. Batista, Jürgen Kurths

Abstract: Dengue is a viral vector-borne infectious disease that affects many countries worldwide, infecting around 390 million people per year. The main outbreaks occur in subtropical and tropical countries. We study here the influence of climate on dengue in Natal (2016-2019), Brazil, Iquitos (2001-2012), Peru, and Barranquilla (2011-2016), Colombia. For the analysis and simulations, we apply Machine Lear… ▽ More Dengue is a viral vector-borne infectious disease that affects many countries worldwide, infecting around 390 million people per year. The main outbreaks occur in subtropical and tropical countries. We study here the influence of climate on dengue in Natal (2016-2019), Brazil, Iquitos (2001-2012), Peru, and Barranquilla (2011-2016), Colombia. For the analysis and simulations, we apply Machine Learning (ML) techniques, especially the Random Forest (RF) algorithm. In addition, regarding a feature in the ML technique, we analyze three possibilities: only dengue cases (D); climate and dengue cases (CD); humidity and dengue cases (HD). Depending on the city, our results show that the climate data can improve or not the forecast. For instance, for Natal, D induces a better forecast. For Iquitos, it is better to use CD. Nonetheless, for Barranquilla, the forecast is better, when we include cases and humidity data. For Natal, when we use more than 64\% and less than 80\% of the time series for training, we obtain results with correlation coefficients ($r$) among 0.917 and 0.949 and mean absolute errors (MAE) among 57.783 and 71.768 for the D case in forecasting. The optimal range for Iquitos is obtained when 79\% up to 88\% of the time series is considered for training. For this case, the best case is CD, having a minimum $r$ equal to 0.850 and maximum 0.887, while values of MAE oscillate among 2.780 and 4.156. For Barranquilla, the optimal range occurs between 72\% until 82\% of length training. In this case, the better approach is HD, where the measures exhibit a minimum $r$ equal to 0.942 and a maximum 0.953, while the minimum and maximum MAE vary between 6.085 and 6.669. We show that the forecast of dengue cases is a challenging problem and climate variables do not always help. However, when we include the mentioned climate variables, the most important one is humidity. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Adaptive Exponential Integrate-and-Fire Model with Fractal Extension

Authors: Diogo L. M. Souza, Enrique C. Gabrick, Paulo R. Protachevicz, Fernando S. Borges, José Trobia, Kelly C. Iarosz, Antonio M. Batista, Iberê L. Caldas, Ervin K. Lenzi

Abstract: The description of neuronal activity has been of great importance in neuroscience. In this field, mathematical models are useful to describe the electrophysical behaviour of neurons. One successful model used for this purpose is the Adaptive Exponential Integrate-and-Fire (Adex), which is composed of two ordinary differential equations. Usually, this model is considered in his standard formulation… ▽ More The description of neuronal activity has been of great importance in neuroscience. In this field, mathematical models are useful to describe the electrophysical behaviour of neurons. One successful model used for this purpose is the Adaptive Exponential Integrate-and-Fire (Adex), which is composed of two ordinary differential equations. Usually, this model is considered in his standard formulation, i.e., with integer order derivatives. In this work, we propose and study the fractal extension of Adex model, which in simple terms corresponds to replacing the integer derivative by non-integer. As non-integer operators, we choose the fractal derivatives. We explore the effects of equal and different orders of fractal derivatives in the firing patterns and mean frequency of the neuron described by the Adex model. Previous results suggest that fractal derivatives can provide a more realistic representation due to the fact that the standard operators are generalized. Our findings show that the fractal order influences the inter-spike intervals and changes the mean firing frequency. In addition, the firing patterns depend not only on the neuronal parameters but also on the order of respective fractal operators. As our main conclusion, the fractal order below the unit value increases the influence of the adaptation mechanism in the spike firing patterns. △ Less

Submitted 9 February, 2024; originally announced March 2024.

Journal ref: Chaos 34 2024

Spiral wave dynamics in a neuronal network model

Authors: Diogo L. M. Souza, Fernando S. Borges, Enrique C. Gabrick, Lucas E. Bentivoglio, Paulo R. Protachevicz, Vagner dos Santos, Ricardo L. Viana, Ibere L. Caldas, Kelly C. Iarosz, Antonio M. Batista, Jürgen Kurths

Abstract: Spiral waves are spatial-temporal patterns that can emerge in different systems as heart tissues, chemical oscillators, ecological networks and the brain. These waves have been identified in the neocortex of turtles, rats, and humans, particularly during sleep-like states. Although their functions in cognitive activities remain until now poorly understood, these patterns are related to cortical ac… ▽ More Spiral waves are spatial-temporal patterns that can emerge in different systems as heart tissues, chemical oscillators, ecological networks and the brain. These waves have been identified in the neocortex of turtles, rats, and humans, particularly during sleep-like states. Although their functions in cognitive activities remain until now poorly understood, these patterns are related to cortical activity modulation and contribute to cortical processing. In this work, we construct a neuronal network layer based on the spatial distribution of pyramidal neurons. Our main goal is to investigate how local connectivity and coupling strength are associated with the emergence of spiral waves. Therefore, we propose a trustworthy method capable of detecting different wave patterns, based on local and global phase order parameters. As a result, we find that the range of connection radius (R) plays a crucial role in the appearance of spiral waves. For R < 20 μm, only asynchronous activity is observed due to small number of connections. The coupling strength (gsyn ) greatly influences the pattern transitions for higher R, where spikes and bursts firing patterns can be observed in spiral and non-spiral waves. Finally, we show that for some values of R and gsyn bistable states of wave patterns are obtained. △ Less

Submitted 28 February, 2024; originally announced March 2024.

Impact of periodic vaccination in SEIRS seasonal model

Authors: Enrique C. Gabrick, Eduardo L. Brugnago, Silvio L. T. de Souza, Kelly C. Iarosz, José D. Szezech Jr., Ricardo L. Viana, Iberê L. Caldas, Antonio M. Batista, Jürgen Kurths

Abstract: We study three different strategies of vaccination in a 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 ra… ▽ More We study three different strategies of vaccination in a 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, the 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 the intensive immunization at inflection of the transmissivity curve. There, 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 immunisation campaign annual reaching for controlling the infection spreading. Regarding the dynamical behaviour of the model, our simulations shows that chaotic and periodic solutions, as well as bi-stable regions, depend on the vaccination parameters range. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Isochronous bifurcations in a two-parameter twist map

Authors: Michele Mugnaine, Bruno B. Leal, Alfredo M. Ozorio de Almeida, Ricardo L. Viana, Iberê L. Caldas

Abstract: Isochronous islands in phase space emerge in twist Hamiltonian systems as a response to multiple resonant perturbations. According to the Poincaré-Birkhoff theorem, the number of islands depends on the system characteristics and the perturbation. We analyze, for the two-parameter standard map, also called two-harmonic standard map, how the island chains are modified as the perturbation amplitude i… ▽ More Isochronous islands in phase space emerge in twist Hamiltonian systems as a response to multiple resonant perturbations. According to the Poincaré-Birkhoff theorem, the number of islands depends on the system characteristics and the perturbation. We analyze, for the two-parameter standard map, also called two-harmonic standard map, how the island chains are modified as the perturbation amplitude increases. We identified three routes for the transition from one chain, associated with one harmonic, to the chain associated with the other harmonic, based on a combination of pitchfork and saddle-node bifurcations. These routes can present intermediate island chains configurations. Otherwise, the destruction of the islands always occurs through the pitchfork bifurcation. △ Less

Submitted 19 December, 2023; originally announced December 2023.

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

Authors: Eduardo L. Brugnago, Enrique C. Gabrick, Kelly C. Iarosz, José D. Szezech Jr., Ricardo L. Viana, Antonio M. Batista, Iberê L. Caldas

Abstract: In this work, we study the dynamics of a SEIRS 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 strutures. We highlight the system multistability, iden… ▽ More In this work, we study the dynamics of a SEIRS 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 strutures. 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 has 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 reaching approximately $50\%$ higher than the chaotic one. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Fractional dynamics and recurrence analysis in cancer model

Authors: Enrique C. Gabrick, Matheus R. Sales, Elaheh Sayari, José Trobia, Ervin K. Lenzi, Fernando da S. Borges, José D. Szezech Jr., Kelly C. Iarosz, Ricardo L. Viana, Iberê L. Caldas, Antonio M. Batista

Abstract: In this work, we analyze the effects of fractional derivatives in the chaotic dynamics of a cancer model. We begin by studying the dynamics of a standard model, {\it i.e.}, with integer derivatives. We study the dynamical behavior by means of the bifurcation diagram, Lyapunov exponents, and recurrence quantification analysis (RQA), such as the recurrence rate (RR), the determinism (DET), and the r… ▽ More In this work, we analyze the effects of fractional derivatives in the chaotic dynamics of a cancer model. We begin by studying the dynamics of a standard model, {\it i.e.}, with integer derivatives. We study the dynamical behavior by means of the bifurcation diagram, Lyapunov exponents, and recurrence quantification analysis (RQA), such as the recurrence rate (RR), the determinism (DET), and the recurrence time entropy (RTE). We find a high correlation coefficient between the Lyapunov exponents and RTE. Our simulations suggest that the tumor growth parameter ($ρ_1$) is associated with a chaotic regime. Our results suggest a high correlation between the largest Lyapunov exponents and RTE. After understanding the dynamics of the model in the standard formulation, we extend our results by considering fractional operators. We fix the parameters in the chaotic regime and investigate the effects of the fractional order. We demonstrate how fractional dynamics can be properly characterized using RQA measures, which offer the advantage of not requiring knowledge of the fractional Jacobian matrix. We find that the chaotic motion is suppressed as $α$ decreases, and the system becomes periodic for $α\lessapprox 0.9966$. We observe limit cycles for $α\in (0.9966,0.899)$ and fixed points for $α<0.899$. The fixed point is determined analytically for the considered parameters. Finally, we discover that these dynamics are separated by an exponential relationship between $α$ and $ρ_1$. Also, the transition depends on a supper transient which obeys the same relationship. △ Less

Submitted 25 August, 2023; originally announced September 2023.

A nontwist field line mapping in a tokamak with ergodic magnetic limiter

Authors: Michele Mugnaine, José D. Szezech Jr., Ricardo L. Viana, Iberê L. Caldas

Abstract: For tokamaks with uniform magnetic shear, Martin and Taylor have proposed a symplectic map 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 non-monotonic profile, we obtain a nont… ▽ More For tokamaks with uniform magnetic shear, Martin and Taylor have proposed a symplectic map 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 non-monotonic profile, we obtain a nontwist map which presents the characteristic properties of degenerate systems, as the twin islands scenario, the 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. △ Less

Submitted 19 December, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Hamiltonian description for magnetic field lines: a tutorial

Authors: Ricardo L. Viana, Michele Mugnaine, Iberê L. Caldas

Abstract: Under certain circumstances, the equations for the magnetic field lines can be recast in a canonical form, after defining a suitable field line Hamiltonian. This analogy is extremely useful for dealing with a variety of problems involving magnetically confined plasmas, like in tokamaks and other toroidal devices, where there is usually one symmetric coordinate which plays the role of time in the c… ▽ More Under certain circumstances, the equations for the magnetic field lines can be recast in a canonical form, after defining a suitable field line Hamiltonian. This analogy is extremely useful for dealing with a variety of problems involving magnetically confined plasmas, like in tokamaks and other toroidal devices, where there is usually one symmetric coordinate which plays the role of time in the canonical equations. In this tutorial paper we review the basics of the Hamiltonian description for magnetic field lines, emphasizing the role of a variational principle and gauge invariance. We present representative applications of the formalism, using cylindrical and magnetic flux coordinates in tokamak plasmas. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Isochronous island bifurcations driven by resonant magnetic perturbations in Tokamaks

Authors: B. B. Leal, I. L. Caldas, M. C. de Sousa, R. L. Viana, A. M. Ozorio de Almeida

Abstract: Recent evidences show that heteroclinic bifurcations in magnetic islands may be caused by the amplitude variation of resonant magnetic perturbations in tokamaks. To investigate the onset of these bifurcations, we consider a large aspect ratio tokamak with an ergodic limiter composed of two pairs of rings that create external primary perturbations with two sets of wave numbers. An individual pair p… ▽ More Recent evidences show that heteroclinic bifurcations in magnetic islands may be caused by the amplitude variation of resonant magnetic perturbations in tokamaks. To investigate the onset of these bifurcations, we consider a large aspect ratio tokamak with an ergodic limiter composed of two pairs of rings that create external primary perturbations with two sets of wave numbers. An individual pair produces hyperbolic and elliptic periodic points, and its associated islands, that are consistent with the Poincaré-Birkhoff fixed point theorem. However, for two pairs producing external perturbations resonant on the same rational surface, we show that different configurations of isochronous island chains may appear on phase space according to the amplitude of the electric currents in each pair of the ergodic limiter. When one of the electric currents increases, isochronous bifurcations take place and new islands are created with the same winding number as the preceding islands. We present examples of bifurcation sequences displaying (a) direct transitions from the island chain configuration generated by one of the pairs to the configuration produced by the other pair, and (b) transitions with intermediate configurations produced by the limiter pairs coupling. Furthermore, we identify shearless bifurcations inside some isochronous islands, originating nonmonotonic local winding number profiles with associated shearless invariant curves. △ Less

Submitted 2 January, 2024; v1 submitted 1 August, 2023; originally announced August 2023.

ExB drift particle transport in tokamaks

Authors: L. A. Osorio-Quiroga, G. C. Grime, M. Roberto, R. L. Viana, Y. Elskens, I. L. Caldas

Abstract: In tokamaks, modification of the plasma profiles can reduce plasma transport, improving particle confinement. However, this improvement is still not completely understood. In this work, we consider a drift wave test particle model to investigate the influence of the electric and magnetic field profiles on plasma transport. Test particle orbits subjected to ExB drift are numerically integrated and… ▽ More In tokamaks, modification of the plasma profiles can reduce plasma transport, improving particle confinement. However, this improvement is still not completely understood. In this work, we consider a drift wave test particle model to investigate the influence of the electric and magnetic field profiles on plasma transport. Test particle orbits subjected to ExB drift are numerically integrated and their transport coefficient is obtained. We conclude that sheared profiles reduce particle transport, even for high amplitude perturbations. In particular, nonmonotonic electric and magnetic fields produce shearless transport barriers, which are particularly resistant to perturbations and reduce even more the transport coefficient. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Comments: 10 pages, 9 figures, 1 table. Published in Brazilian Journal of Physics

Journal ref: Braz J Phys, 53, 96 (2023)

Shearless curve breakup in the biquadratic nontwist map

Authors: Gabriel C. Grime, Marisa Roberto, Ricardo L. Viana, Yves Elskens, Iberê L. Caldas

Abstract: Nontwist area-preserving maps violate the twist condition along shearless invariant curves, which act as transport barriers in phase space. Recently, some plasma models have presented multiple shearless curves in phase space and these curves can break up independently. In this paper, we describe the different shearless curve breakup scenarios of the so-called biquadratic nontwist map, a recently p… ▽ More Nontwist area-preserving maps violate the twist condition along shearless invariant curves, which act as transport barriers in phase space. Recently, some plasma models have presented multiple shearless curves in phase space and these curves can break up independently. In this paper, we describe the different shearless curve breakup scenarios of the so-called biquadratic nontwist map, a recently proposed area-preserving map derived from a plasma model, that captures the essential behavior of systems with multiple shearless curves. Three different scenarios are found and their dependence on the system parameters is analyzed. The results indicate a relation between shearless curve breakup and periodic orbit reconnection-collision sequences. In addition, even after a shearless curve breakup, the remaining curves inhibit global transport. △ Less

Submitted 27 May, 2023; originally announced May 2023.

Comments: 12 pages, 7 figures

Confining and escaping magnetic field lines in Tokamaks: Analysis via symplectic map

Authors: Matheus S. Palmero, Iberê L. Caldas

Abstract: In magnetically confined plasma, it is possible to qualitatively describe the magnetic field configuration via phase spaces of suitable symplectic maps. These phase spaces are of mixed type, where chaos coexists with regular motion, and the complete understanding of the chaotic transport is a challenge that, when overcome, may provide further knowledge into the behaviour of confined fusion plasma.… ▽ More In magnetically confined plasma, it is possible to qualitatively describe the magnetic field configuration via phase spaces of suitable symplectic maps. These phase spaces are of mixed type, where chaos coexists with regular motion, and the complete understanding of the chaotic transport is a challenge that, when overcome, may provide further knowledge into the behaviour of confined fusion plasma. This work presents two numerical investigations into characteristics of mixed phase spaces which model distinct magnetic configurations in tokamaks under different perturbation regimes. The first approach relies on a recurrence-based analysis of ensembles of chaotic trajectories to detect open field lines that widely differ from the average. The second focuses on the transient dynamical behaviour of field lines before they escape the systems. These two methods provide insights into the influence of stickiness and invariant manifolds on the evolution of chaotic trajectories, improving our understanding of how these features affect the chaotic transport and diffusion properties in mixed phase spaces. These theoretical and numerical approaches may enhance our comprehension of confined plasma behaviour and plasma-wall interactions. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: Preprint submitted to Fundamental Plasma Physics

Bidimensional Symplectic Maps

Authors: Felipe G. Souza, Gabriel C. Grime, Iberê L. Caldas

Abstract: Symplectic maps can provide a straightforward and accurate way to visualize and quantify the dynamics of conservative systems with two degrees of freedom. These maps can be easily iterated from the simplest computers to obtain trajectories with great accuracy. Their usage arises in many fields, including celeste mechanics, plasma physics, chemistry, and so on. In this paper we introduce two exampl… ▽ More Symplectic maps can provide a straightforward and accurate way to visualize and quantify the dynamics of conservative systems with two degrees of freedom. These maps can be easily iterated from the simplest computers to obtain trajectories with great accuracy. Their usage arises in many fields, including celeste mechanics, plasma physics, chemistry, and so on. In this paper we introduce two examples of symplectic maps, the standard and the standard non-twist map, exploring the phase space transformation as their control parameters are varied. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: 15 pages, 9 figures

Stickiness and recurrence plots: an entropy-based approach

Authors: Matheus R. Sales, Michele Mugnaine, José D. Szezech Jr., Ricardo L. Viana, Iberê L. Caldas, Norbert Marwan, Jürgen Kurths

Abstract: 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 (RP), 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 recurr… ▽ More 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 (RP), 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. △ Less

Submitted 27 February, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: 16 pages, 7 figures

Fractal and Wada escape basins in the chaotic particle drift motion in tokamaks

Authors: Leonardo C. Souza, Amanda C. Mathias, Iberê L. Caldas, Yves Elskens, Ricardo L. Viana

Abstract: The ${\bf E}\times{\bf B}$ drift motion of particles in tokamaks provides valuable information on the turbulence-driven anomalous transport. One of the characteristic features of the drift motion dynamics is the presence of chaotic orbits for which the guiding center can experience large-scale drifts. If one or more exits are placed within that chaotic orbit, the corresponding escape basins struct… ▽ More The ${\bf E}\times{\bf B}$ drift motion of particles in tokamaks provides valuable information on the turbulence-driven anomalous transport. One of the characteristic features of the drift motion dynamics is the presence of chaotic orbits for which the guiding center can experience large-scale drifts. If one or more exits are placed within that chaotic orbit, the corresponding escape basins structure is complicated and, indeed, exhibits fractal structures. We investigate those structures through a number of numerical diagnostics, tailored to quantify the final-state uncertainty related to the fractal escape basins. We estimate the escape basin boundary dimension through the uncertainty exponent method, and quantify final-state uncertainty by the basin entropy and the basin boundary entropy. Finally, we describe the so-called Wada property, for the case of three or more escape basins. This property is verified both qualitatively and quantitatively, using a grid approach. △ Less

Submitted 16 February, 2024; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 15 pages, 10 figures

Unpredictability in seasonal infectious diseases spread

Authors: Enrique C. Gabrick, Elaheh Sayari, Paulo R. Protachevicz, José D. Szezech Jr., Kelly C. Iarosz, Silvio L. T. de Souza, Alexandre C. L. Almeida, Ricardo L. Viana, Iberê L. Caldas, Antonio M. Batista

Abstract: In this work, we study the unpredictability of seasonal infectious diseases considering a SEIRS model with seasonal forcing. To investigate the dynamical behaviour, we compute bifurcation diagrams type hysteresis and their respective Lyapunov exponents. Our results from bifurcations and the largest Lyapunov exponent show bistable dynamics for all the parameters of the model. Choosing the inverse o… ▽ More In this work, we study the unpredictability of seasonal infectious diseases considering a SEIRS model with seasonal forcing. To investigate the dynamical behaviour, we compute bifurcation diagrams type hysteresis and their respective Lyapunov exponents. Our results from bifurcations and the largest Lyapunov exponent show bistable dynamics for all the parameters of the model. Choosing the inverse of latent period as control parameter, over 70% of the interval comprises the coexistence of periodic and chaotic attractors, bistable dynamics. Despite the competition between these attractors, the chaotic ones are preferred. The bistability occurs in two wide regions. One of these regions is limited by periodic attractors, while periodic and chaotic attractors bound the other. As the boundary of the second bistable region is composed of periodic and chaotic attractors, it is possible to interpret these critical points as tipping points. In other words, depending on the latent period, a periodic attractor (predictability) can evolve to a chaotic attractor (unpredictability). Therefore, we show that unpredictability is associated with bistable dynamics preferably chaotic, and, furthermore, there is a tipping point associated with unpredictable dynamics. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Chaotic saddles and interior crises in a dissipative nontwist system

Authors: Rodrigo Simile Baroni, Ricardo Egydio de Carvalho, Iberê Luiz Caldas, Ricardo Luiz Viana, Philip J Morrison

Abstract: We consider a dissipative version of the standard nontwist map. Nontwist systems present a robust transport barrier, called the shearless curve, that becomes the shearless attractor when dissipation is introduced. This attractor can be regular or chaotic depending on the control parameters. Chaotic attractors can undergo sudden and qualitative changes as a parameter is varied. These changes are ca… ▽ More We consider a dissipative version of the standard nontwist map. Nontwist systems present a robust transport barrier, called the shearless curve, that becomes the shearless attractor when dissipation is introduced. This attractor can be regular or chaotic depending on the control parameters. Chaotic attractors can undergo sudden and qualitative changes as a parameter is varied. These changes are called crises, and at an interior crisis the attractor suddenly expands. Chaotic saddles are nonattracting chaotic sets that play a fundamental role in the dynamics of nonlinear systems, they are responsible for chaotic transients, fractal basin boundaries, chaotic scattering and they mediate interior crises. In this work we discuss the creation of chaotic saddles in a dissipative nontwist system and the interior crises they generate. We show how the presence of two saddles increase the transient times and analyze the phenomenon of crisis induced intermittency. △ Less

Submitted 24 February, 2023; v1 submitted 13 November, 2022; originally announced November 2022.

Plastic neural network with transmission delays promotes equivalence between function and structure

Authors: P. R. Protachevicz, F. S. Borges, A. M. Batista, M. S. Baptista, I. L. Caldas, E. E. N. Macau, E. L. Lameu

Abstract: The brain is formed by cortical regions that are associated with different cognitive functions. Neurons within the same region are more likely to connect than neurons in distinct regions, making the brain network to have characteristics of a network of subnetworks. The values of synaptic delays between neurons of different subnetworks are greater than those of the same subnetworks. This difference… ▽ More The brain is formed by cortical regions that are associated with different cognitive functions. Neurons within the same region are more likely to connect than neurons in distinct regions, making the brain network to have characteristics of a network of subnetworks. The values of synaptic delays between neurons of different subnetworks are greater than those of the same subnetworks. This difference in communication time between neurons has consequences on the firing patterns observed in the brain, which is directly related to changes in neural connectivity, known as synaptic plasticity. In this work, we build a plastic network of Hodgkin-Huxley neurons in which the connectivity modifications follow a spike-time dependent rule. We define an internal-delay among neurons communicating within the same subnetwork, an external-delay for neurons belonging to distinct subnetworks, and study how these communicating delays affect the entire network dynamics. We observe that the neuronal network exhibits a specific connectivity configuration for each synchronised pattern. Our results show how synaptic delays and plasticity work together to promote the formation of structural coupling among the neuronal subnetworks. We conclude that plastic neuronal networks are able to promote equivalence between function and structure meaning that topology emerges from behaviour and behaviour emerges from topology, creating a complex dynamical process where topology adapts to conform with the plastic rules and firing patterns reflect the changes on the synaptic weights. △ Less

Submitted 17 April, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Biquadratic Nontwist Map: a model for shearless bifurcations

Authors: G. C. Grime, M. Roberto, R. L. Viana, Y. Elskens, I. L. Caldas

Abstract: Area-preserving nontwist maps are used to describe a broad range of physical systems. In those systems, the violation of the twist condition leads to nontwist characteristic phenomena, such as reconnection-collision sequences and shearless invariant curves that act as transport barriers in the phase space. Although reported in numerical investigations, the shearless bifurcation, i.e., the emergenc… ▽ More Area-preserving nontwist maps are used to describe a broad range of physical systems. In those systems, the violation of the twist condition leads to nontwist characteristic phenomena, such as reconnection-collision sequences and shearless invariant curves that act as transport barriers in the phase space. Although reported in numerical investigations, the shearless bifurcation, i.e., the emergence scenario of multiple shearless curves, is not well understood. In this work, we derive an area-preserving map as a local approximation of a particle transport model for confined plasmas. Multiple shearless curves are found in this area-preserving map, with the same shearless bifurcation scenario numerically observed in the original model. Due to its symmetry properties and simple functional form, this map is proposed as a model to study shearless bifurcations. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 16 pages, 10 figures

Shaping the edge radial electric field to create shearless transport barriers in tokamaks

Authors: L. A. Osorio-Quiroga, M. Roberto, I. L. Caldas, R. L. Viana, Y. Elskens

Abstract: In tokamak-confined plasmas, particle transport can be reduced by modifying the radial electric field. In this paper, we investigate the influence of both a well-like and a hill-like shaped radial electric field profile on the creation of shearless transport barriers (STBs) at the plasma edge, which are a type of barrier that can prevent chaotic transport and are related to the presence of extreme… ▽ More In tokamak-confined plasmas, particle transport can be reduced by modifying the radial electric field. In this paper, we investigate the influence of both a well-like and a hill-like shaped radial electric field profile on the creation of shearless transport barriers (STBs) at the plasma edge, which are a type of barrier that can prevent chaotic transport and are related to the presence of extreme values in the rotation number profile. For that, we apply an ExB drift model to describe test particle orbits in large aspect-ratio tokamaks. We show how these barriers depend on the electrostatic fluctuation amplitudes and on the width and depth (height) of the radial electric field well-like (hill-like) profile. We find that, as the depth (height) increases, the STB at the plasma edge becomes more resistant to fluctuations, enabling access to an improved confinement regime that prevents chaotic transport. We also present parameter spaces with the radial electric field parameters, indicating the STB existence for several electric field configurations at the plasma edge, for which we obtain a fractal structure at the barrier/non-barrier frontier, typical of quasi-integrable Hamiltonian systems. △ Less

Submitted 2 June, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 12 pages and 8 figures

Journal ref: Fundamental Plasma Physics, 6, 100023(2023)

Effect of two vaccine doses in the SEIR epidemic model using a stochastic cellular automaton

Authors: Enrique C. Gabrick, Paulo R. Protachevicz, Antonio M. Batista, Kelly C. Iarosz, Silvio L. T. de Souza, Alexandre C. L. Almeida, José D. Szezech Jr, Michele Mugnaine, Iberê L. Caldas

Abstract: In this work, to support decision making of immunisation strategies, we propose the inclusion of two vaccination doses in the SEIR model considering a stochastic cellular automaton. We analyse three different scenarios of vaccination:… ▽ More In this work, to support decision making of immunisation strategies, we propose the inclusion of two vaccination doses in the SEIR model considering a stochastic cellular automaton. We analyse three different scenarios of vaccination: $i) unlimited doses, (ii) limited doses into susceptible individuals, and (iii) limited doses randomly distributed overall individuals. Our results suggest that the number of vaccinations and time to start the vaccination is more relevant than the vaccine efficacy, delay between the first and second doses, and delay between vaccinated groups. The scenario (i) shows that the solution can converge early to a disease-free equilibrium for a fraction of individuals vaccinated with the first dose. In the scenario (ii), few two vaccination doses divided into a small number of applications reduce the number of infected people more than into many applications. In addition, there is a low waste of doses for the first application and an increase of the waste in the second dose. The scenario (iii) presents an increase in the waste of doses from the first to second applications more than the scenario $(ii)$. In the scenario (iii), the total of wasted doses increases linearly with the number of applications. Furthermore, the number of effective doses in the application of consecutive groups decays exponentially overtime. △ Less

Submitted 31 August, 2022; originally announced September 2022.

Internal energy exchanges and chaotic dynamics in an intrinsically coupled system

Authors: M. C. de Sousa, A. B. Schelin, F. A. Marcus, R. L. Viana, I. L. Caldas

Abstract: Intrinsically coupled nonlinear systems present different oscillating components that exchange energy among themselves. A paradigmatic example is the spring pendulum, which displays spring, pendulum, and coupled oscillations. We analyze the energy exchanges among the oscillations, and obtain that it is enhanced for chaotic orbits. Moreover, the highest rates of energy exchange for the coupling occ… ▽ More Intrinsically coupled nonlinear systems present different oscillating components that exchange energy among themselves. A paradigmatic example is the spring pendulum, which displays spring, pendulum, and coupled oscillations. We analyze the energy exchanges among the oscillations, and obtain that it is enhanced for chaotic orbits. Moreover, the highest rates of energy exchange for the coupling occur along the homoclinic tangle of the primary hyperbolic point embedded in a chaotic sea. The results show a clear relation between internal energy exchanges and the dynamics of a coupled system. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: Submitted for publication

Journal ref: Physics Letters A, 453, 128481 (2022)

Shearless bifurcations in particle transport for reversed shear tokamaks

Authors: G. C. Grime, M. Roberto, R. L. Viana, Y. Elskens, I. L. Caldas

Abstract: Some internal transport barriers in tokamaks have been related to the vicinity of extrema of the plasma equilibrium profiles. This effect is numerically investigated by considering the guiding-center trajectories of plasma particles undergoing ExB drift motion, considering that the electric field has a stationary nonmonotonic radial profile and an electrostatic fluctuation. In addition, the equili… ▽ More Some internal transport barriers in tokamaks have been related to the vicinity of extrema of the plasma equilibrium profiles. This effect is numerically investigated by considering the guiding-center trajectories of plasma particles undergoing ExB drift motion, considering that the electric field has a stationary nonmonotonic radial profile and an electrostatic fluctuation. In addition, the equilibrium configuration has a nonmonotonic safety factor profile. The numerical integration of the equations of motion yields a symplectic map with shearless barriers. By changing the parameters of the safety factor profile, the appearance, and breakup of these shearless curves are observed. The successive shearless curves breakup and recovering is explained using concepts from bifurcation theory. We also present bifurcation sequences associated to the creation of multiple shearless curves. Physical consequences of scenarios with multiple shearless curves are discussed. △ Less

Submitted 23 September, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

Comments: 18 pages, 9 figures. Replacement improved the text

The effect of time delay for synchronisation suppression in neuronal networks

Authors: Matheus Hansen, Paulo R. Protachevicz, Kelly C. Iarosz, Ibere L. Caldas, Antonio M. Batista, Elbert E. N. Macau

Abstract: We study the time delay in the synaptic conductance for suppression of spike synchronisation in a random network of Hodgkin Huxley neurons coupled by means of chemical synapses. In the first part, we examine in detail how the time delay acts over the network during the synchronised and desynchronised neuronal activities. We observe a relation between the neuronal dynamics and the syaptic conductan… ▽ More We study the time delay in the synaptic conductance for suppression of spike synchronisation in a random network of Hodgkin Huxley neurons coupled by means of chemical synapses. In the first part, we examine in detail how the time delay acts over the network during the synchronised and desynchronised neuronal activities. We observe a relation between the neuronal dynamics and the syaptic conductance distributions. We find parameter values in which the time delay has high effectiveness in promoting the suppression of spike synchronisation. In the second part, we analyse how the delayed neuronal networks react when pulsed inputs with different profiles (periodic, random, and mixed) are applied on the neurons. We show the main parameters responsible for inducing or not synchronous neuronal oscillations in delayed networks. △ Less

Submitted 30 April, 2022; originally announced May 2022.

Diffusion transitions in a 2D periodic lattice

Authors: Matheus J. Lazarotto, Iberê L. Caldas, Yves Elskens

Abstract: Spatial diffusion of particles in periodic potential models has provided a good framework for studying the role of chaos in global properties of classical systems. Here a bidimensional "soft" billiard, classically modeled from an optical lattice hamiltonian system, is used to study diffusion transitions with the control parameters variation. Sudden transitions between normal and ballistic regimes… ▽ More Spatial diffusion of particles in periodic potential models has provided a good framework for studying the role of chaos in global properties of classical systems. Here a bidimensional "soft" billiard, classically modeled from an optical lattice hamiltonian system, is used to study diffusion transitions with the control parameters variation. Sudden transitions between normal and ballistic regimes are found and characterized by inspection of the topological changes in phase-space. Transitions correlated with increases in global stability area are shown to occur for energy levels where local maxima points become accessible, deviating trajectories approaching it. These instabilities promote a slowing down of the dynamics and an island myriad bifurcation phenomenon, along with the suppression of long flights within the lattice. Other diffusion regime variations occurring during small intervals of control parameters are shown to be related to the emergence of a set of orbits with long flights, thus altering the total average displacement for long integration times but without global changes in phase-space. △ Less

Submitted 13 May, 2022; v1 submitted 26 October, 2021; originally announced November 2021.

Comments: 14 pages, 14 figures, submitted to Communications in Nonlinear Sciences and Numerical Simulation

Journal ref: Communications in Nonlinear Science and Numerical Simulation, Volume 112, September 2022, 106525

Onset of internal transport barriers in tokamaks

Authors: L. A. Osorio, M. Roberto, I. L. Caldas, R. L. Viana, Y. Elskens

Abstract: Barriers have been identified in magnetically confined plasmas reducing the particle transport and improving the confinement. One of them, the primary shearless barriers are associated to extrema of non-monotonic plasma profiles. Previously, we identified these barriers in a model described by a map that allows the integration of charged particles motion in drift waves for a long time scale. In th… ▽ More Barriers have been identified in magnetically confined plasmas reducing the particle transport and improving the confinement. One of them, the primary shearless barriers are associated to extrema of non-monotonic plasma profiles. Previously, we identified these barriers in a model described by a map that allows the integration of charged particles motion in drift waves for a long time scale. In this work, we show how the existence of these robust barriers depends on the fluctuation amplitude and on the electric shear. Moreover, we also find control parameter intervals for which these primary barriers onset and break-up are recurrent. Another noticeable feature, in these transitions, is the appearance of a layer of particle trajectory stickiness after the shearless barrier break-up or before its onset. Besides the mentioned primary barriers, we also observe sequences of secondary shearless barriers, not reported before, created and destroyed by a sequence of bifurcations as the main control parameters, the fluctuation amplitude and electric shear, are varied. Furthermore, in these bifurcations, we also find hitherto unknown double and triple secondary shearless barriers which constitute a noticeable obstacle to the chaotic transport. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 12 pages, 28 figures. This article has been accepted by the Journal Physics of Plasmas. It can be found at: https://doi.org/10.1063/5.0056428

Journal ref: Physics of Plasmas, 28(8), 082305(2021)

Simulation of deterministic compartmental models for infectious diseases dynamics

Authors: Antonio M Batista, Silvio L T Souza, Kelly C Iarosz, Alexandre C L Almeida, José D Szezech, Enrique C Gabrick, Michele Mugnaine, Gefferson L dos Santos, Iberê L Caldas

Abstract: Infectious diseases are caused by pathogenic microorganisms and can spread through different ways. Mathematical models and computational simulation have been used extensively to investigate the transmission and spread of infectious diseases. In other words, mathematical model simulation can be used to analyse the dynamics of infectious diseases, aiming to understand the effects and how to control… ▽ More Infectious diseases are caused by pathogenic microorganisms and can spread through different ways. Mathematical models and computational simulation have been used extensively to investigate the transmission and spread of infectious diseases. In other words, mathematical model simulation can be used to analyse the dynamics of infectious diseases, aiming to understand the effects and how to control the spread. In general, these models are based on compartments, where each compartment contains individuals with the same characteristics, such as susceptible, exposed, infected, and recovered. In this paper, we cast further light on some classical epidemic models, reporting possible outcomes from numerical simulation. Furthermore, we provide routines in a repository for simulations. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Slippery-Sticky Transition of Interfacial Fluid Slip

Authors: Thiago F. Viscondi, Adriano Grigolo, Iberê L. Caldas, Julio R. Meneghini

Abstract: The influence of temperature on interfacial fluid slip, as measured by molecular-dynamics simulations of a Couette flow comprising a Lennard-Jones fluid and rigid crystalline walls, is examined as a function of the fluid-solid interaction strength. Two different types of thermal behavior are observed, namely, the slippery and sticky cases. The first is characterized by a steep and unlimited increa… ▽ More The influence of temperature on interfacial fluid slip, as measured by molecular-dynamics simulations of a Couette flow comprising a Lennard-Jones fluid and rigid crystalline walls, is examined as a function of the fluid-solid interaction strength. Two different types of thermal behavior are observed, namely, the slippery and sticky cases. The first is characterized by a steep and unlimited increase of the slip length at low temperatures, while the second presents a vanishing slip length in this regime. As the temperature increases in relation to a characteristic value, both cases converge to finite slip lengths. A recently proposed analytical model is found to well describe both thermal behaviors, also predicting the slippery-sticky transition that occurs at a critical value of the fluid-solid interaction parameter, for which, according to the model, fluid particles experience a smooth average energy landscape at the interface. △ Less

Submitted 29 June, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 10 pages, 4 figures

Journal ref: T. F. Viscondi, A. Grigolo, I. L. Caldas, and J. R. Meneghini, Physics of Fluids 33, 062012 (2021)

Transport barriers in symplectic maps

Authors: R. L. Viana, I. L. Caldas, J. D. Szezech Jr., A. M. Batista, C. V. Abud, A. B. Schelin, M. Mugnaine, M. S. Santos, B. B. Leal, B. Bartoloni, A. C. Mathias, J. V. Gomes, P. J. Morrison

Abstract: Chaotic transport is a subject of paramount importance in a variety of problems in plasma physics, specially those related to anomalous transport and turbulence. On the other hand, a great deal of information on chaotic transport can be obtained from simple dynamical systems like two-dimensional area-preserving (symplectic) maps, where powerful mathematical results like KAM theory are available. I… ▽ More Chaotic transport is a subject of paramount importance in a variety of problems in plasma physics, specially those related to anomalous transport and turbulence. On the other hand, a great deal of information on chaotic transport can be obtained from simple dynamical systems like two-dimensional area-preserving (symplectic) maps, where powerful mathematical results like KAM theory are available. In this work we review recent works on transport barriers in area-preserving maps, focusing on systems which do not obey the so-called twist property. For such systems KAM theory no longer holds everywhere and novel dynamical features show up as non-resistive reconnection, shearless curves and shearless bifurcations. After presenting some general features using a standard nontwist mapping, we consider magnetic field line maps for magnetically confined plasmas in tokamaks. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Growth and performance of the periodic orbits of a nonlinear driven oscillator

Authors: Danilo Rodrigues de Lima, Iberê Luiz Caldas

Abstract: Periodic orbits are fundamental to understand the dynamics of nonlinear systems. In this work, we focus on two aspects of interest regarding periodic orbits, in the context of a dissipative mapping, derived from a prototype model of a non-linear driven oscillator with fast relaxation and a limit cycle. For this map, we show numerically the exponential growth of periodic orbits quantity in certain… ▽ More Periodic orbits are fundamental to understand the dynamics of nonlinear systems. In this work, we focus on two aspects of interest regarding periodic orbits, in the context of a dissipative mapping, derived from a prototype model of a non-linear driven oscillator with fast relaxation and a limit cycle. For this map, we show numerically the exponential growth of periodic orbits quantity in certain regions of the parameter space and provide an analytical bound for such growth rate, by making use of the transition matrix associated with a given periodic orbit. Furthermore, we give numerical evidence to support that optimal orbits, those that maximize time averages, are often unstable periodic orbits with low period, by numerically comparing their performance under a family of sinusoidal functions. △ Less

Submitted 20 December, 2020; originally announced December 2020.

Comments: 15 pages, 6 figures, 3 appendices

Synchronization and attractors in a model simulating social jetlag

Authors: Flávia M. Ruziska, Iberê L. Caldas

Abstract: Much work has been done to investigate social jetlag, a misalignment between the biological clock and the social agenda caused by exposition to different light inputs, that causes several health issues. To investigate synchronization and attractors due to a sequence of light inputs, we introduce an extension of a model, previously used to describe jetlag caused by a single change in the light inpu… ▽ More Much work has been done to investigate social jetlag, a misalignment between the biological clock and the social agenda caused by exposition to different light inputs, that causes several health issues. To investigate synchronization and attractors due to a sequence of light inputs, we introduce an extension of a model, previously used to describe jetlag caused by a single change in the light input. The synchronization to the light input is sensitive to the control parameters of the system and to the light input periods. Depending on the parameter set, the observed synchronization is with one or another successive light input. Most of the solutions have the period of a light input. However, for some parameters, we also observed higher period, chaotic solutions, and bistability. △ Less

Submitted 23 December, 2020; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: article 21 pages 8 figures, supplementary material 10 pages, 7 figures

Mathematical model of brain tumour growth with drug resistance

Authors: José Trobia, Kun Tian, Antonio Marcos Batista, Celso Grebogi, Hai-Peng Ren, Moises Souza Santos, Paulo Ricardo Protachevicz, Fernando da Silva Borges, José Danilo Szezech Jr, Ricardo Luiz Viana, Iberê Luiz Caldas, Kelly Cristiane Iarosz

Abstract: Brain tumours are masses of abnormal cells that can grow in an uncontrolled way in the brain. There are different types of malignant brain tumours. Gliomas are malignant brain tumours that grow from glial cells and are identified as astrocytoma, oligodendroglioma, and ependymoma. We study a mathematical model that describes glia-neuron interaction, glioma, and chemotherapeutic agent. In this work,… ▽ More Brain tumours are masses of abnormal cells that can grow in an uncontrolled way in the brain. There are different types of malignant brain tumours. Gliomas are malignant brain tumours that grow from glial cells and are identified as astrocytoma, oligodendroglioma, and ependymoma. We study a mathematical model that describes glia-neuron interaction, glioma, and chemotherapeutic agent. In this work, we consider drug sensitive and resistant glioma cells. We show that continuous and pulsed chemotherapy can kill glioma cells with a minimal loss of neurons. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Low-dimensional chaos in the single wave model for self-consistent wave-particle Hamiltonian

Authors: J. V. Gomes, M. C. de Sousa, R. L. Viana, I. L. Caldas, Y. Elskens

Abstract: We analyze nonlinear aspects of the self-consistent wave-particle interaction using Hamiltonian dynamics in the single wave model, where the wave is modified due to the particle dynamics. This interaction plays an important role in the emergence of plasma instabilities and turbulence. The simplest case, where one particle (N = 1) is coupled with one wave (M = 1), is completely integrable, and the… ▽ More We analyze nonlinear aspects of the self-consistent wave-particle interaction using Hamiltonian dynamics in the single wave model, where the wave is modified due to the particle dynamics. This interaction plays an important role in the emergence of plasma instabilities and turbulence. The simplest case, where one particle (N = 1) is coupled with one wave (M = 1), is completely integrable, and the nonlinear effects reduce to the wave potential pulsating while the particle either remains trapped or circulates forever. On increasing the number of particles (N = 2, M = 1), integrability is lost and chaos develops. Our analyses identify the two standard ways for chaos to appear and grow (the homoclinic tangle born from a separatrix, and the resonance overlap near an elliptic fixed point). Moreover, a strong form of chaos occurs when the energy is high enough for the wave amplitude to vanish occasionally. △ Less

Submitted 23 July, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Sub-diffusive behavior in the Standard Map

Authors: Matheus S. Palmero, Gabriel I. Díaz, Iberê L. Caldas, Igor. M. Sokolov

Abstract: In this work, we investigate the presence of sub-diffusive behavior in the Chirikov-Taylor Standard Map. We show that the stickiness phenomena, present in the mixed phase space of the map setup, can be characterized as a Continuous Time Random Walk model and connected to the theoretical background for anomalous diffusion. Additionally, we choose a variant of the Ulam method to numerically approxim… ▽ More In this work, we investigate the presence of sub-diffusive behavior in the Chirikov-Taylor Standard Map. We show that the stickiness phenomena, present in the mixed phase space of the map setup, can be characterized as a Continuous Time Random Walk model and connected to the theoretical background for anomalous diffusion. Additionally, we choose a variant of the Ulam method to numerically approximate the Perron-Frobenius operator for the map, allowing us to calculate the anomalous diffusion exponent via an eigenvalue problem, compared to the solution of the Fractional Diffusion Equation. The results here corroborate other findings in the literature of anomalous transport in Hamiltonian maps and can be suitable to describe transport properties of other dynamical systems. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: Submitted to EPJ-ST

Influence of autapses on synchronisation in neural networks with chemical synapses

Authors: P. R. Protachevicz, K. C. Iarosz, I. L. Caldas, C. G. Antonopoulos, A. M. Batista, J. Kurths

Abstract: A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronisation. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autap… ▽ More A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronisation. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autapses, to study the effect of the latter on synchronous behaviour. We consider time delay in the conductance of the pre-synaptic neuron for excitatory and inhibitory connections. Interestingly, in neural networks consisting of both excitatory and inhibitory neurons, we uncover that synchronous behaviour depends on their synapse type. Our results provide evidence on the synchronous and desynchronous activities that emerge in random neural networks with chemical, inhibitory and excitatory synapses where neurons are equipped with autapses. △ Less

Submitted 6 November, 2020; v1 submitted 5 November, 2020; originally announced November 2020.

Measure, dimension, and complexity of the transient motion in Hamiltonian systems

Authors: Vitor M. de Oliveira, Matheus S. Palmero, Iberê L. Caldas

Abstract: Hamiltonian systems that are either open, leaking, or contain holes in the phase space possess solutions that eventually escape the system's domain. The motion described by such escape orbits before crossing the escape threshold can be understood as a transient behavior. In this work, we introduce a numerical method to visually illustrate and quantify the transient motion in Hamiltonian systems ba… ▽ More Hamiltonian systems that are either open, leaking, or contain holes in the phase space possess solutions that eventually escape the system's domain. The motion described by such escape orbits before crossing the escape threshold can be understood as a transient behavior. In this work, we introduce a numerical method to visually illustrate and quantify the transient motion in Hamiltonian systems based on the transient measure, a finite-time version of the natural measure. We apply this method to two physical systems: the single-null divertor tokamak, described by a symplectic map; and the Earth-Moon system, as modeled by the planar circular restricted three-body problem. Our results portray how different locations for the ensemble of initial conditions may lead to different transient dynamical scenarios in both systems. We show that these scenarios can be properly quantified from a geometrical aspect, the transient correlation dimension, and a dynamical aspect, the transient complexity coefficient. △ Less

Submitted 6 October, 2021; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 12 pages, 12 figures, 2 tables

Controlling Escape in the Standard Map

Authors: Gabriel I. Díaz, Matheus S. Palmero, Iberê Luiz Caldas, Edson D. Leonel

Abstract: We investigate how the diffusion exponent is affected by controlling small domains in the phase space.The main Kolomogorov-Arnold-Moser - KAM island of the Standard Map is considered to validate the investigation. The bifurcation scenario where the periodic island emits smaller resonance regions is considered and we show how closing paths escape from the island shore by controlling points and henc… ▽ More We investigate how the diffusion exponent is affected by controlling small domains in the phase space.The main Kolomogorov-Arnold-Moser - KAM island of the Standard Map is considered to validate the investigation. The bifurcation scenario where the periodic island emits smaller resonance regions is considered and we show how closing paths escape from the island shore by controlling points and hence making the diffusion exponent smaller. We notice the bigger controlled area the smaller the diffusion exponent. We show that controlling around the hyperbolic points associated to the bifurcation is better than a random control to reduce the diffusion exponent. The recurrence plot shows us channels of escape and a control applied there reduces the diffusion exponent. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Effects of drug resistance in the tumour-immune system with chemotherapy treatment

Authors: José Trobia, Enrique C Gabrick, Evandro G Seifert, Fernando S Borges, Paulo R Protachevicz, José D Szezech Jr, Kelly C Iarosz, Moises S Santos, Iberê L Caldas, Kun Tian, Hai-Peng Ren, Celso Grebogi, Antonio M Batista

Abstract: Cancer is a term used to refer to a large set of diseases. The cancerous cells grow and divide and, as a result, they form tumours that grow in size. The immune system recognise the cancerous cells and attack them, though, it can be weakened by the cancer. One type of cancer treatment is chemotherapy, which uses drugs to kill cancer cells. Clinical, experimental, and theoretical research has been… ▽ More Cancer is a term used to refer to a large set of diseases. The cancerous cells grow and divide and, as a result, they form tumours that grow in size. The immune system recognise the cancerous cells and attack them, though, it can be weakened by the cancer. One type of cancer treatment is chemotherapy, which uses drugs to kill cancer cells. Clinical, experimental, and theoretical research has been developed to understand the dynamics of cancerous cells with chemotherapy treatment, as well as the interaction between tumour growth and immune system. We study a mathematical model that describes the cancer growth, immune system response, and chemotherapeutic agents. The immune system is composed of resting cells that are converted to hunting cells to combat the cancer. In this work, we consider drug sensitive and resistant cancer cells. We show that the tumour growth can be controlled not only by means of different chemotherapy protocols, but also by the immune system that attacks both sensitive and resistant cancer cells. Furthermore, for all considered protocols, we demonstrate that the time delay from resting to hunting cells plays a crucial role in the combat against cancer cells. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Influence of inhibitory synapses on the criticality of excitable neuronal networks

Authors: F S Borges, P R Protachevicz, V Santos, M S Santos, E C Gabrick, K C Iarosz, E L Lameu, M S Baptista, I L Caldas, A M Batista

Abstract: In this work, we study the dynamic range of a neuronal network of excitable neurons with excitatory and inhibitory synapses. We obtain an analytical expression for the critical point as a function of the excitatory and inhibitory synaptic intensities. We also determine an analytical expression that gives the critical point value in which the maximal dynamic range occurs. Depending on the mean conn… ▽ More In this work, we study the dynamic range of a neuronal network of excitable neurons with excitatory and inhibitory synapses. We obtain an analytical expression for the critical point as a function of the excitatory and inhibitory synaptic intensities. We also determine an analytical expression that gives the critical point value in which the maximal dynamic range occurs. Depending on the mean connection degree and coupl\-ing weights, the critical points can exhibit ceasing or ceaseless dynamics. However, the dynamic range is equal in both cases. We observe that the external stimulus mask some effects of self-sustained activity (ceaseless dynamic) in the region where the dynamic range is calculated. In these regions, the firing rate is the same for ceaseless dynamics and ceasing activity. Furthermore, we verify that excitatory and inhibitory inputs are approximately equal for a network with a large number of connections, showing excitatory-inhibitory balance as reported experimentally. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Order-chaos-order and invariant manifolds in the bounded planar Earth-Moon system

Authors: Vitor M. de Oliveira, Priscilla A. Sousa-Silva, Iberê L. Caldas

Abstract: In this work, we investigate the Earth-Moon system, as modeled by the planar circular restricted three-body problem, and relate its dynamical properties to the underlying structure associated with specific invariant manifolds. We consider a range of Jacobi constant values for which the neck around the Lagrangian point $L_1$ is always open but the orbits are bounded due to Hill stability. First, we… ▽ More In this work, we investigate the Earth-Moon system, as modeled by the planar circular restricted three-body problem, and relate its dynamical properties to the underlying structure associated with specific invariant manifolds. We consider a range of Jacobi constant values for which the neck around the Lagrangian point $L_1$ is always open but the orbits are bounded due to Hill stability. First, we show that the system displays three different dynamical scenarios in the neighborhood of the Moon: two mixed ones, with regular and chaotic orbits, and an almost entirely chaotic one in between. We then analyze the transitions between these scenarios using the Monodromy matrix theory and determine that they are given by two specific types of bifurcations. After that, we illustrate how the phase space configurations, particularly the shapes of stability regions and stickiness, are intrinsically related to the hyperbolic invariant manifolds of the Lyapunov orbits around $L_1$ and also to the ones of some particular unstable periodic orbits. Lastly, we define transit time in a manner that is useful to depict dynamical trapping and show that the traced geometrical structures are also connected to the transport properties of the system. △ Less

Submitted 25 September, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 22 pages, 12 figures

Journal ref: Celest Mech Dyn Astr 132, 51 (2020)

Wave-particle interactions in a long traveling wave tube with upgraded helix

Authors: M. C. de Sousa, F. Doveil, Y. Elskens, I. L. Caldas

Abstract: We investigate the interaction of electromagnetic waves and electron beams in a 4 meters long traveling wave tube (TWT). The device is specially designed to simulate beam-plasma experiments without appreciable noise. This TWT presents an upgraded slow wave structure (SWS) that results in more precise measurements and makes new experiments possible. We introduce a theoretical model describing wave… ▽ More We investigate the interaction of electromagnetic waves and electron beams in a 4 meters long traveling wave tube (TWT). The device is specially designed to simulate beam-plasma experiments without appreciable noise. This TWT presents an upgraded slow wave structure (SWS) that results in more precise measurements and makes new experiments possible. We introduce a theoretical model describing wave propagation through the SWS and validated by the experimental dispersion relation, impedance, phase and group velocities. We analyze nonlinear effects arising from the beam-wave interaction, such as the modulation of the electron beam and the wave growth and saturation process. When the beam current is low, the wave growth coefficient and saturation amplitude follow the linear theory predictions. However, for high values of current, nonlinear space charge effects become important and these parameters deviate from the linear predictions, tending to a constant value. After saturation, we also observe trapping of the beam electrons, which alters the wave amplitude along the TWT. △ Less

Submitted 4 August, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: Submitted for publication

Journal ref: Physics of Plasmas 27, 093108 (2020)