Nontwist Map Research Articles

Magnetic structures in the explicitly time-dependent nontwist map.

We investigate how the magnetic structures of the plasma change in a large aspect ratio tokamak perturbed by an ergodic magnetic limiter, when a system parameter is non-adiabatically varied in time. We model such a scenario by considering the Ullmann-Caldas nontwist map, where we introduce an explicit time-dependence to the ratio of the limiter and plasma currents. We apply the tools developed recently in the field of chaotic Hamiltonian systems subjected to parameter drift. Namely, we follow trajectory ensembles initially forming Kolmogorov Arnold Moser (KAM) tori and island chains in the autonomous configuration space. With a varying parameter, these ensembles, called snapshot tori, develop time-dependent shapes. An analysis of the time evolution of the average distance of point pairs in such an ensemble reveals that snapshot tori go through a transition to chaos, with a positive Lyapunov exponent. We find empirical power-law relationships between both the Lyapunov exponent and the beginning of the transition to chaos (the so-called critical instant), as a function of the rate of the parameter drift, with the former showing an increasing trend and the latter a decreasing trend. We conclude that, in general, coherent tori and magnetic islands tend to break up and become chaotic as the perturbation increases, similar to the case of subsequent constant perturbations. However, because of the continuous drift, some structures can persist longer and exist even at perturbation values where they would not be observable in the constant perturbation case.

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

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. Here we describe a 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 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 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 structure that we call 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 nonuniformly distributed points.

Lagrangian descriptors: The shearless curve and the shearless attractor.

Hamiltonian systems with a nonmonotonic frequency profile are called nontwist. One of the key properties of such systems, depending on adjustable parameters, is the presence of a robust transport barrier in the phase space called the shearless curve, which becomes the equally robust shearless attractor when dissipation is introduced. We consider the standard nontwist map with and without dissipation. We derive analytical expressions for the Lagrangian descriptor (LD) for the unperturbed map and show how they are related to the rotation number profile. We show how the LDs can reconstruct finite segments of the invariant manifolds for the perturbed map. In the conservative case, we demonstrate how the LDs distinguish the chaotic seas from regular structures. The LDs also provide a remarkable tool to identify when the shearless curve is destroyed: we present a fractal boundary, in the parameter space, for the existence or not of the shearless torus. In the dissipative case, we show how the LDs can be used to localize point attractors and the shearless attractor and distinguish their basins of attraction.

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

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

Shearless curve breakup in the biquadratic nontwist map

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.

Chaotic saddles and interior crises in a dissipative nontwist system.

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, and 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 increases the transient times and we analyze the phenomenon of crisis induced intermittency.

Mapas padrão twist e não-twist

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 map and the standard non-twist map, exploring the phase space transformation as their control parameters are varied. Furthermore, for the non-twist map, we also present, that its chaotic motion is separated in phase space by the non-twist transport barrier.

Dynamics, multistability, and crisis analysis of a sine-circle nontwist map.

We propose a one-dimensional dynamical system, the sine-circle nontwist map, that can be considered a local approximation of the standard nontwist map and an extension of the paradigmatic sine-circle map. The map depends on three parameters, exhibiting a simple mathematical form but with a rich dynamical behavior. We identify periodic, quasiperiodic, and chaotic solutions for different parameter sets with the Lyapunov exponent and Slater's theorem. From the bifurcation analysis, we determine two bifurcation lines, those that depend on just two of the control parameters, for which the bifurcation that occurs is of the saddle-node type. In order to investigate multistability, we analyze the bifurcation diagrams in the two directions of parameter variation and we observe some regions of hysteresis, representing the coexistence of different attractors. We also analyze different multistable scenarios, as single attractor, coexistence of periodic attractors, coexistence of chaotic and periodic attractors, chaotic behavior, and coexistence of different chaotic bands, by the Lyapunov exponent and the analysis of the domain occupied by the solutions. From the parameter spaces constructed, we observe the prevalence of single attractor and only chaotic behavior scenarios. The multistable scenario is, mostly, formed by different periodic attractors. Lastly, we analyze the crisis in chaotic attractors and we identify the interior and the boundary crisis. From our results, the boundary crisis plays a key role for the extinction of multistability.

Destruction and resurgence of the quasiperiodic shearless attractor.

We consider a dissipative version of the standard nontwist map. It is known that nontwist systems may present a robust transport barrier, called shearless curve, that gives rise to an attractor that retains some of its properties when dissipation is introduced. This attractor is known as shearless attractor, and it may be quasiperiodic or chaotic depending on the control parameters. We describe a route for the destruction and resurgence of the quasiperiodic shearless attractor by analyzing the manifolds of the unstable periodic orbits (UPOs) which are fixed points of the map. We show that the shearless attractor is destroyed by a collision with the UPOs and it resurges after the reconnection of the unstable manifolds of different UPOs.

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

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

Ratchet current in nontwist Hamiltonian systems.

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

Fractal structures in the parameter space of nontwist area-preserving maps.

Fractal structures are very common in the phase space of nonlinear dynamical systems, both dissipative and conservative, and can be related to the final state uncertainty with respect to small perturbations on initial conditions. Fractal structures may also appear in the parameter space, since parameter values are always known up to some uncertainty. This problem, however, has received less attention, and only for dissipative systems. In this work we investigate fractal structures in the parameter space of two conservative dynamical systems: the standard nontwist map and the quartic nontwist map. For both maps there is a shearless invariant curve in the phase space that acts as a transport barrier separating chaotic orbits. Depending on the values of the system parameters this barrier can break up. In the corresponding parameter space the set of parameter values leading to barrier breakup is separated from the set not leading to breakup by a curve whose properties are investigated in this work, using tools as the uncertainty exponent and basin entropies. We conclude that this frontier in parameter space is a complicated curve exhibiting both smooth and fractal properties, that are characterized using the uncertainty dimension and basin and basin boundary entropies.

Transport barriers with shearless attractors.

We present a mechanism to generate a sequence of shearless curves or attractors to form a band of transport barriers. We consider the labyrinthic nontwist standard map to prepare a scenario with three shearless curves. Dissipation is introduced and three shearless attractors coexist, very close to each other. In both cases a collective transport barrier is formed.

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

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

Dynamical characterization of transport barriers in nontwist Hamiltonian systems.

The turnstile provides us a useful tool to describe the flux in twist Hamiltonian systems. Thus, its determination allows us to find the areas where the trajectories flux through barriers. We show that the mechanism of the turnstile can increase the flux in nontwist Hamiltonian systems. A model which captures the essence of these systems is the standard nontwist map, introduced by del Castillo-Negrete and Morrison. For selected parameters of this map, we show that chaotic trajectories entering in resonances zones can be explained by turnstiles formed by a set of homoclinic points. We argue that for nontwist systems, if the heteroclinic points are sufficiently close, they can connect twin-islands chains. This provides us a scenario where the trajectories can cross the resonance zones and increase the flux. For these cases the escape basin boundaries are nontrivial, which demands the use of an appropriate characterization. We applied the uncertainty exponent and the entropies of the escape basin boundary in order to quantify the degree of unpredictability of the asymptotic trajectories.

Rigorous Computer-Assisted Application of KAM Theory: A Modern Approach

In this paper, we present and illustrate a general methodology to apply KAM theory in particular problems, based on an a posteriori approach. We focus on the existence of real analytic quasi-periodic Lagrangian invariant tori for symplectic maps. The purpose is to verify the hypotheses of a KAM theorem in an a posteriori format: Given a parameterization of an approximately invariant torus, we have to check non-resonance (Diophantine) conditions, non-degeneracy conditions and certain inequalities to hold. To check such inequalities, we require to control the analytic norm of some functions that depend on the map, the ambient structure and the parameterization. To this end, we propose an efficient computer-assisted methodology, using fast Fourier transform, having the same asymptotic cost of using the parameterization method for obtaining numerical approximations of invariant tori. We illustrate our methodology by proving the existence of invariant curves for the standard map (up to $$\varepsilon =0.9716$$ ), meandering curves for the non-twist standard map and 2-dimensional tori for the Froeschle map.

Area-preserving maps models of gyroaveraged E×B chaotic transport

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 E × 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 the counting of recurrences times are used to quantify the dependence on the Larmor radius of the threshold for the destruction of transport barriers.

Breakup of inverse golden mean shearless tori in the two-frequency standard nontwist map

The breakup of shearless invariant tori with winding number ω = ( 5 − 1 ) / 2 (inverse golden mean) is studied using Greeneʼs residue criterion in the recently derived two-frequency or extended standard nontwist map (ESNM). Depending on the frequency ratio, the ESNM has or does not have a particular spatial symmetry. If the symmetry is present, the breakup is shown to be the same as in the standard nontwist map; if not, the results are very different. ► We study the breakup of shearless invariant tori with inverse golden mean winding number in an area-preserving nontwist map perturbed by two harmonics. ► Torus breakup is studied using Greeneʼs residue criterion. ► Symmetry of the map (or lack of) is found to affect the details of the breakup.

Effective transport barriers in nontwist systems

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. These barriers have been identified in symplectic nontwist maps that model such zonal flows. We use the so-called standard nontwist map, a paradigmatic example of nontwist systems, to analyze the parameter dependence of the transport through a broken shearless barrier. On 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. Moreover, we show how these stickiness structures, and the concomitant transport reduction in the shearless region, are determined by a homoclinic tangle of the remaining dominant twin island chains. We use the finite-time rotation number, a recently proposed diagnostic, to identify transport barriers that separate different regions of stickiness. The identified barriers are comparable to those obtained by using finite-time Lyapunov exponents.

Breakup of shearless invariant tori in cubic and quartic nontwist maps

The effect of symmetry on invariant torus breakup in nontwist maps is investigated. In particular, the breakup of shearless invariant tori with winding number ω = ( 5 - 1 ) / 2 (inverse golden mean) and ω = 2 - 1 (an inverse silver mean) is studied numerically using Greene’s residue criterion in a cubic and a quartic nontwist map. The details of the breakup are compared to those previously obtained for the standard nontwist map, which has the same particular spatial symmetry as the quartic map. The cubic map lacks this symmetry. The results show that if the symmetry exists, the details of the breakup are the same as in the standard nontwist map. If the symmetry does not exist, the breakup is shown to be different.