Search

We investigate the influence of the finite Larmor radius on the dynamics of guiding-center test particles subjected to an $\mathbf{E} \times \mathbf{B}$ drift in a large aspect-ratio tokamak. For that, we adopt the drift-wave test particle transport model presented by W. Horton [Physics of Plasmas \textbf{5}, 3910 (1998)] and introduce a second-order gyro-averaged extension, which accounts for the finite Larmor radius effect that arises from a spatially varying electric field. Using this extended model, we numerically examine the influence of the finite Larmor radius on chaotic transport and the formation of transport barriers. For non-monotonic plasma profiles, we show that the twist condition of the dynamical system, i.e.,\ KAM theorem's non-degeneracy condition for the Hamiltonian, is violated along a special curve, which, under non-equilibrium conditions, exhibits significant resilience to destruction, thereby inhibiting chaotic transport. This curve acts as a robust barrier to transport and is usually called shearless transport barrier. While varying the amplitude of the electrostatic perturbations, we analyze bifurcation diagrams of the shearless barriers and escape rates of orbits to explore the impact of the finite Larmor radius on controlling chaotic transport. Our findings show that increasing the Larmor radius enhances the robustness of transport barriers, as larger electrostatic perturbation amplitudes are required to disrupt them. Additionally, as the Larmor radius increases, even in the absence of transport barriers, we observe a reduction in the escape rates, indicating a decrease in chaotic transport., Comment: Forthcoming in Physics of Plasmas

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.

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\'e-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.

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., Comment: 10 pages, 9 figures, 1 table. Published in Brazilian Journal of Physics

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.

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., Comment: 16 pages, 10 figures

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., Comment: 12 pages and 8 figures

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., Comment: Submitted for publication

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., Comment: 18 pages, 9 figures. Replacement improved the text

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., Comment: 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

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.

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.

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.

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.

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., Comment: Submitted for publication

Noise appears in the brain due to various sources, such as ionic channel fluctuations and synaptic events. They affect the activities of the brain and influence neuron action potentials. Stochastic differential equations have been used to model firing patterns of neurons subject to noise. In this work, we consider perturbing noise in the adaptive exponential integrate-and-fire (AEIF) neuron. The AEIF is a two-dimensional model that describes different neuronal firing patterns by varying its parameters. Noise is added in the equation related to the membrane potential. We show that a noise current can induce continuous and noncontinuous transitions in neuronal interspike intervals. Moreover, we show that the noncontinuous transition occurs mainly for parameters close to the border between tonic spiking and burst activities of the neuron without noise.

Numerical experiments of the statistical evolution of an ensemble of non-interacting particles in a time-dependent billiard with inelastic collisions, reveals the existence of three statistical regimes for the evolution of the speeds ensemble, namely, diffusion plateau, normal growth/exponential decay and stagnation. These regimes are linked numerically to the transition from Gauss-like to Boltzmann-like speed distributions. Further, the different evolution regimes are obtained analytically through velocity-space diffusion analysis. From these calculations the asymptotic root mean square of speed, initial plateau, and the growth/decay rates for intermediate number of collisions are determined in terms of the system parameters. The analytical calculations match the numerical experiments and point to a dynamical mechanism for thermalization, where inelastic collisions and a high-dimensional phase space lead to a bounded diffusion in the velocity space towards a stationary distribution function with a kind of reservoir temperature determined by the boundary oscillation amplitude and the restitution coefficient.

Excessively high, neural synchronisation has been associated with epileptic seizures, one of the most common brain diseases worldwide. A better understanding of neural synchronisation mechanisms can thus help control or even treat epilepsy. In this paper, we study neural synchronisation in a random network where nodes are neurons with excitatory and inhibitory synapses, and neural activity for each node is provided by the adaptive exponential integrate-and-fire model. In this framework, we verify that the decrease in the influence of inhibition can generate synchronisation originating from a pattern of desynchronised spikes. The transition from desynchronous spikes to synchronous bursts of activity, induced by varying the synaptic coupling, emerges in a hysteresis loop due to bistability where abnormal (excessively high synchronous) regimes exist. We verify that, for parameters in the bistability regime, a square current pulse can trigger excessively high (abnormal) synchronisation, a process that can reproduce features of epileptic seizures. Then, we show that it is possible to suppress such abnormal synchronisation by applying a small-amplitude external current on less than 10% of the neurons in the network. Our results demonstrate that external electrical stimulation not only can trigger synchronous behaviour, but more importantly, it can be used as a means to reduce abnormal synchronisation and thus, control or treat effectively epileptic seizures.

In tokamaks, internal transport barriers, produced by modifications of the plasma current profile, reduce particle transport and improve plasma confinement. The triggering of the internal transport barriers and their dependence on the plasma profiles is a key nonlinear dynamics problem still under investigation. We consider the onset of shearless invariant curves inside the plasma which create internal transport barriers. A non-integrable drift-kinetic model is used to describe particle transport driven by drift waves and to investigate these shearless barriers onset in tokamaks. We show that for some currently observed plasma profiles shearless particle transport barriers can be triggered by properly modifying the electric field profile and the influence of non-resonant modes in the barriers onset. In particular, we show that a broken barrier can be restored by enhancing non-resonant modes., Comment: 9 pages, 15 figures, original article

Brain plasticity refers to brain's ability to change neuronal connections, as a result of environmental stimuli, new experiences, or damage. In this work, we study the effects of the synaptic delay on both the coupling strengths and synchronisation in a neuronal network with synaptic plasticity. We build a network of Hodgkin-Huxley neurons, where the plasticity is given by the Hebbian rules. We verify that without time delay the excitatory synapses became stronger from the high frequency to low frequency neurons and the inhibitory synapses increases in the opposite way, when the delay is increased the network presents a non-trivial topology. Regarding the synchronisation, only for small values of the synaptic delay this phenomenon is observed.

Invariant manifolds of unstable periodic orbits organize the dynamics of chaotic orbits in phase space. They provide insight into the mechanisms of transport and chaotic advection and have important applications in physical situations involving three-dimensional flows. The numerical determination of invariant manifolds for planar maps is a problem on its own. Efficient and practical techniques to describe these curves must be made available to the scientific community. 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 in a refinement procedure that prevents redundant calculations occurring in non-refinement approaches and the approximated method is based on an intuitive geometrical decomposition of the manifold in bare and fine details. The resulting approximated manifold is precise when compared to the exact manifold, and its calculation is computationally more efficient, making it ideal for mappings involving intensive calculations like inverse functions or numerical integration of ODEs between crossings in a surface of section.

Nowadays, divertors are used in the main tokamaks to control the magnetic field and to improve the plasma confinement. In this article, we present analytical symplectic maps describing Poincar\'e maps of the magnetic field lines in confined plasmas with a single null poloidal divertor. Initially, we present a divertor map and the tokamap for a diverted configuration. We also introduce the Ullmann map for a diverted plasma, whose control parameters are determined from tokamak experiments. Finally, an explicit, area-preserving and integrable magnetic field line map for a single-null divertor tokamak is obtained using a trajectory integration method to represent toroidal equilibrium magnetic surfaces. In this method, we also give examples of onset of chaotic field lines at the plasma edge due to resonant perturbations.

The dynamics of wave-particle interactions in magnetized plasmas restricts the wave amplitude to moderate values for particle beam acceleration from rest energy. We analyze how a perturbing invariant robust barrier modifies the phase space of the system and enlarges the wave amplitude interval for particle acceleration. For low values of the wave amplitude, the acceleration becomes effective for particles with initial energy close to the rest energy. For higher values of the wave amplitude, the robust barrier controls chaos in the system and restores the acceleration process. We also determine the best position for the perturbing barrier in phase space in order to increase the final energy of the particles., Comment: Submitted for publication

The characterisation of neuronal connectivity is one of the most important matters in neuroscience. In this work, we show that a recently proposed informational quantity, the causal mutual information, employed with an appropriate methodology, can be used not only to correctly infer the direction of the underlying physical synapses, but also to identify their excitatory or inhibitory nature, considering easy to handle and measure bivariate time-series. The success of our approach relies on a surprising property found in neuronal networks by which non-adjacent neurons do "understand" each other (positive mutual information), however this exchange of information is not capable of causing effect (zero transfer entropy). Remarkably, inhibitory connections, responsible for enhancing synchronisation, transfer more information than excitatory connections, known to enhance entropy in the network. We also demonstrate that our methodology can be used to correctly infer directionality of synapses even in the presence of dynamic and observational Gaussian noise, and is also successful in providing the effective directionality of inter modular connectivity, when only mean fields can be measured.

In this work, we study the dynamic range in a neuronal network modelled by cellular automaton. We consider deterministic and non-deterministic rules to simulate electrical and chemical synapses. Chemical synapses have an intrinsic time-delay and are susceptible to parameter variations guided by learning Hebbian rules of behaviour. Our results show that chemical synapses can abruptly enhance sensibility of the neural network, a manifestation that can become even more predominant if learning rules of evolution are applied to the chemical synapses.

Intrinsically nonlinear coupled systems present different oscillating components that exchange energy among themselves. We present a new approach to deal with such energy exchanges and to investigate how it depends on the system control parameters. The method consists in writing the total energy of the system, and properly identifying the energy terms for each component and, especially, their coupling. To illustrate the proposed approach, we work with the bi-dimensional spring pendulum, which is a paradigm to study nonlinear coupled systems, and is used as a model for several systems. For the spring pendulum, we identify three energy components, resembling the spring and pendulum like motions, and the coupling between them. With these analytical expressions, we analyze the energy exchange for individual trajectories, and we also obtain global characteristics of the spring pendulum energy distribution by calculating spatial and time average energy components for a great number of trajectories (periodic, quasi-periodic and chaotic) throughout the phase space. Considering an energy term due to the nonlinear coupling, we identify regions in the parameter space that correspond to strong and weak coupling. The presented procedure can be applied to nonlinear coupled systems to reveal how the coupling mediates internal energy exchanges, and how the energy distribution varies according to the system parameters., Comment: Submitted for publication

Neuronal systems have been modeled by complex networks in different description levels. Recently, it has been verified that networks can simultaneously exhibit one coherent and other incoherent domain, known as chimera states. In this work, we study the existence of chimera states in a network considering the connectivity matrix based on the cat cerebral cortex. The cerebral cortex of the cat can be separated in 65 cortical areas organised into the four cognitive regions: visual, auditory, somatosensory-motor and frontolimbic. We consider a network where the local dynamics is given by the Hindmarsh-Rose model. The Hindmarsh-Rose equations are a well known model of neuronal activity that has been considered to simulate membrane potential in neuron. Here, we analyse under which conditions chimera states are present, as well as the affects induced by intensity of coupling on them. We observe the existence of chimera states in that incoherent structure can be composed of desynchronised spikes or desynchronised bursts. Moreover, we find that chimera states with desynchronised bursts are more robust to neuronal noise than with desynchronised spikes.

A statistical study of finite Larmor radius (FLR) effects on transport driven by electrostatic drift-waves is presented. The study is based on a reduced discrete Hamiltonian dynamical system known as the gyro-averaged standard map (GSM). In this system, FLR effects are incorporated through the gyro-averaging of a simplified weak-turbulence model of electrostatic fluctuations. Formally, the GSM is a modified version of the standard map in which the perturbation amplitude, $K_0$, becomes $K_0 J_0(\hat{\rho})$, where $J_0$ is the zeroth-order Bessel function and $\hat{\rho}$ is the Larmor radius. Assuming a Maxwellian probability density function (pdf) for $\hat{\rho}$, we compute analytically and numerically the pdf and the cumulative distribution function of the effective drift-wave perturbation amplitude $K_0 J_0(\hat{\rho})$. Using these results we compute the probability of loss of confinement (i.e., global chaos), $P_{c}$, and the probability of trapping in the main drift-wave resonance, $P_{t}$. It is shown that $P_{c}$ provides an upper bound for the escape rate, and that $P_{t}$ provides a good estimate of the particle trapping rate. The analytical results are compared with direct numerical Monte-Carlo simulations of particle transport.

Non-axisymmetric stationary magnetic perturbations lead to the formation of homoclinic tangles near the divertor magnetic saddle in tokamak discharges. These tangles intersect the divertor plates in static helical structures that delimit the regions reached by open magnetic field lines reaching the plasma column and leading the charged particles to the strike surfaces by parallel transport. In this article we introduce a non-axisymmetric rotating magnetic perturbation to model the time development of the three-dimensional magnetic field of a single-null DIII-D tokamak discharge developing a rotating tearing mode. The stable and unstable manifolds of the asymmetric magnetic saddle are calculated through an adaptive method providing the manifold cuts at a given poloidal plane and the strike surfaces. For the modeled shot, the experimental heat pattern and its time development are well described by the rotating unstable manifold, indicating the emergence of homoclinic lobes in a rotating frame due to the plasma instabilities. In the model it is assumed that the magnetic field is created by a stationary axisymmetric plasma current and a set of rotating internal helical filamentary currents. The currents in the filaments are adjusted to match the waveforms of the magnetic probes at the mid-plane and the rotating magnetic field is introduced as a perturbation to the axisymmetric field obtained from a Grad-Shafranov equilibrium reconstruction code.

We investigate the dynamical behavior of a non-ideal Duffing oscillator, a system composed of a mass-spring-pendulum driven by a DC motor with limited power supply. To identify new features on Duffing oscillator parameter space due to the limited power supply, we provide an extensive numerical characterization in the bi-parameter space by using Lyapunov exponents. Following this procedure, we identify remarkable new periodic windows, the ones known as Arnold tongues and also shrimp-shaped structures. Such windows appear in highly organized distribution with typical self-similar structures for the shrimps, and, surprisingly, codimension-2 bifurcation as a point of accumulations for the tongues.

In this paper, we study the effects of spike timing-dependent plasticity on synchronisation in a network of Hodgkin-Huxley neurons. Neuron plasticity is a flexible property of a neuron and its network to change temporarily or permanently their biochemical, physiological, and morphological characteristics, in order to adapt to the environment. Regarding the plasticity, we consider Hebbian rules, specifically for spike timing-dependent plasticity (STDP), and with regard to network, we consider that the connections are randomly distributed. We analyse the synchronisation and desynchronisation according to an input level and probability of connections. Moreover, we verify that the transition for synchronisation depends on the neuronal network architecture, and the external perturbation level.

Discrete maps have been extensively used to model 2-dimensional chaotic transport in plasmas and fluids. Here we focus on area-preserving maps describing finite Larmor radius (FLR) effects on ${\bf E} \times {\bf B}$ chaotic transport in magnetized plasmas with zonal flows perturbed by electrostatic drift waves. FLR effects are included by gyro-averaging the Hamiltonians of the maps which, depending on the zonal flow profile, can have monotonic or non-monotonic frequencies. In the limit of zero Larmor radius, the monotonic frequency map reduces to the standard Chirikov-Taylor map, and, in the case of non-monotonic frequency, the map reduces to the standard nontwist map. We show that in both cases FLR leads to chaos suppression, changes in the stability of fixed points, and robustness of transport barriers. FLR effects are also responsible for changes in the phase space topology and zonal flow bifurcations. Dynamical systems methods based on recurrence time statistics are used to quantify the dependence on the Larmor radius of the threshold for the destruction of transport barriers., Comment: Accepted for publication in Physics of Plasmas

In toroidally confined plasmas, the Grad-Shafranov equation, in general a non-linear PDE, describes the hydromagnetic equilibrium of the system. This equation becomes linear when the kinetic pressure is proportional to the poloidal magnetic flux and the squared poloidal current is a quadratic function of it. In this work, the eigenvalue of the associated homogeneous equation is related with the safety factor on the magnetic axis, the plasma beta and the Shafranov shift, then, the adjustable parameters of the particular solution are bounded through physical constrains. The poloidal magnetic flux becomes a linear superposition of independent solutions and its parameters are adjusted with a non-linear fitting algorithm. This method is used to find hydromagnetic equilibria with normal and reversed magnetic shear and defined values of the elongation, triangularity, aspect-ratio, and X-point(s). The resultant toroidal and poloidal beta, the safety factor at the $95\%$ flux surface and the plasma current are in agreement with usual experimental values for high beta discharges and the model can be used locally to describe reversed magnetic shear equilibria., Comment: 8 pages, 7 figures, 14 references

Master-slave systems have been intensively investigated to model the application of chaos to communications. We considered Colpitts oscillators coupled according to a master-slave configuration to study chaos synchronisation. We revealed the existence of super persistent transients in this coupled system. Moreover, we showed that an additive noise added to the slave system may suppress chaos synchronisation. When synchronisation is not suppressed, the noise induces longer transients.

Nonlinear dynamical systems, ranging from insect populations to lasers and chemical reactions, might exhibit sensitivity to small perturbations in their control parameters, resulting in uncertainties on the predictability of tunning parameters that lead these systems to either a chaotic or a periodic behavior. By quantifying such uncertainties in four different classes of nonlinear systems, we show that this sensitivity is to be expected because the boundary between the sets of parameters leading to chaos and to periodicity is fractal. Moreover, the dimension of this fractal boundary was shown to be roughly the same for these classes of systems. Using an heuristic model for the way periodic windows appear in parameter spaces, we provide an explanation for the universal character of this fractal boundary.

We analyze a Hamiltonian model with five wire loops that delineates the magnetic surfaces of the tokamak ITER, including a similar safety factor profile and the X-point related to the presence of a poloidal divertor. Non-axisymmetric magnetic perturbations are added by external coils, similar to the correction coils installed at the tokamak DIII-D and those that will be installed at ITER. To show the influence of magnetic perturbations on the field line escape, we integrate numerically the field line differential equations and obtain the footprints and deposition patterns on the divertor plate. Moreover, we show that the homoclinic tangle describes the deposition patterns in the divertor plate, agreeing with results observed in sophisticated simulation codes. Additionally, we show that while chaotic lines escape to the divertor plates, some of them are trapped, for many toroidal turns, in complex structures around magnetic islands, embedded in the chaotic region, giving rise to the so called stickiness effect characteristic of chaotic Hamiltonian systems. Finally, we introduce a random collisional term to the field line mapping to investigate stickiness alterations due to particle collisions. Within this model, we conclude that, even reduced by collisions, stickiness still influences the field line transport.

We present simulated figures of the diverted magnetic field lines of the tokamak ITER, obtained by numerically integrating a Hamiltonian model with electrical currents in five wire loops and control coils. We show evidences of a sticky island embedded in the chaotic region near the divertor plates, which traps magnetic field lines for many toroidal turns increasing their connection lengths to these plates.

We analyze the dynamics of a relativistic particle moving in a uniform magnetic field and perturbed by a standing electrostatic wave. We show that a pulsed wave produces an infinite number of perturbative terms with the same winding number, which may generate islands in the same region of phase space. As a consequence, the number of isochronous island chains varies as a function of the wave parameters. We observe that in all the resonances, the number of chains is related to the amplitude of the various resonant terms. We determine analytically the position of the periodic points and the number of island chains as a function of the wave number and wave period. Such information is very important when one is concerned with regular particle acceleration, since it is necessary to adjust the initial conditions of the particle to obtain the maximum acceleration., Comment: Submitted

In parameter space of nonlinear dynamical systems, windows of periodic states are aligned following routes of period-adding configuring periodic window sequences. In state space of driven nonlinear oscillators, we determine the torsion associated with the periodic states and identify regions of uniform torsion in the window sequences. Moreover, we find that the measured of torsion differs by a constant between successive windows in periodic window sequences. We call this phenomenon as torsion-adding. Finally, combining the torsion and the period adding rules, we deduce a general rule to obtain the asymptotic winding number in the accumulation limit of such periodic window sequences.

In fluids and plasmas with zonal flow reversed shear, a peculiar kind of transport barrier appears in the shearless region, one that is associated with a proper route of transition to chaos. Using a symplectic nontwist maps, which model such zonal flows, we investigate these barriers. In particular, the standard nontwist map, a paradigm of nontwist systems, is used to analyze the parameter dependence of the transport through a broken shearless barrier. Varying a proper control parameter, we identify the onset of structures with high stickiness that give rise to an effective barrier near the broken shearless curve. Several diagnostics are used: stickiness is measured by escape time and transmissivity plots, while barriers that separate different regions of stickiness are identified by employing two indicators in the phase space, the finite-time rotation number and the finite-time Lyapunov exponent., Comment: 20 pages, 9 figures

We analyze the behavior of a relativistic particle moving under the influence of a uniform magnetic field and a stationary electrostatic wave. We work with a set of pulsed waves that allows us to obtain an exact map for the system. We also use a method of control for near-integrable Hamiltonians that consists in the addition of a small and simple control term to the system. This control term creates invariant tori in phase space that prevent chaos from spreading to large regions and make the controlled dynamics more regular. We show numerically that the control term just slightly modifies the system but is able to drastically reduce chaos with a low additional cost of energy. Moreover, we discuss how the control of chaos and the consequent recovery of regular trajectories in phase space are useful to improve regular particle acceleration., Comment: 8 pages, 2 figures. Published in Physical Review E

In the bi-dimensional parameter space of driven oscillators, shrimp-shaped periodic windows are immersed in chaotic regions. For two of these oscillators, namely, Duffing and Josephson junction, we show that a weak harmonic perturbation replicates these periodic windows giving rise to parameter regions correspondent to periodic orbits. The new windows are composed of parameters whose periodic orbits have periodicity and pattern similar to stable and unstable periodic orbits already existent for the unperturbed oscillator. These features indicate that the reported replicate periodic windows are associated with chaos control of the considered oscillators.

Lagrangian coherent structures are effective barriers, sticky regions, that separate phase space regions of different dynamical behavior. The usual way to detect such structures is via finite-time Lyapunov exponents. We show that similar results can be obtained for single-frequency systems from finite-time rotation numbers, which are much faster to compute. We illustrate our claim by considering examples of continuous and discrete-time dynamical systems of physical interest., Comment: 4 pages, 3 figures

The existence of a special periodic window in the two-dimensional parameter space of an experimental Chua's circuit is reported. One of the main reasons that makes such a window special is that the observation of one implies that other similar periodic windows must exist for other parameter values. However, such a window has never been experimentally observed, since its size in parameter space decreases exponentially with the period of the periodic attractor. This property imposes clear limitations for its experimental detection., Comment: 4.2 pages, 4 figures

We show experimentally the scenario of a two-frequency torus $T^2$ breakdown, in which a global bifurcation occurs due to the collision of a torus with an unstable periodic orbit, creating a heteroclinic saddle connection, followed by an intermittent behavior.

We show that, in periodically perturbed chaotic systems, Phase Synchronization appears, associated to a special type of stroboscopic map, in which not only averages quantities are equal to invariants of the perturbation, the angular frequency, but also it exists a very large number of non-transient transformations, possibly infinity. In cases where there is not phase synchronization there is either only transitive transformations on the attractor, or a finite number of non-transitive transformations. We base our statements in experimental and numerical results from the sinusoidally perturbed Chua's circuit.

Concepts from the Ergodic Theory are used to describe the existence of non-transitive maps in attractors of phase synchronous chaotic systems. It is shown that for a class of phase-coherent systems, e.g. the sinusoidally forced Chua's circuit and two coupled R{\"o}ssler oscillators, phase synchronization implies that such maps exist. These ideas are also extended to other coupled chaotic systems. In addition, a phase for a chaotic attractor is defined from the tangent vector of the flow. Finally, it is discussed how these maps can be used to real time detection of phase synchronization in experimental systems.

We investigate the statistics of recurrences to finite size intervals for chaotic dynamical systems. We find that the typical distribution presents an exponential decay for almost all recurrence times except for a few short times affected by a kind of memory effect. We interpret this effect as being related to the unstable periodic orbits inside the interval. Although it is restricted to a few short times it changes the whole distribution of recurrences. We show that for systems with strong mixing properties the exponential decay converges to the Poissonian statistics when the width of the interval goes to zero. However, we alert that special attention to the size of the interval is required in order to guarantee that the short time memory effect is negligible when one is interested in numerically or experimentally calculated Poincar{\'e} recurrence time statistics., Comment: 8 pages and 12 figures.Corrected version