Wada basins and chaotic invariant sets in the Hénon-Heiles system

Relativistic chaotic scattering: Unveiling scaling laws for trapped trajectories

In this paper we study different types of phase space structures which appear in the context of relativistic chaotic scattering. By using the relativistic version of the Hénon-Heiles Hamiltonian, we numerically study the topology of different kind of exit basins and compare it with the case of low velocities in which the Newtonian version of the system is valid. Specifically, we numerically study the escapes in the phase space, in the energy plane, and in the β plane, which richly characterize the dynamics of the system. In all cases, fractal structures are present, and the escaping dynamics is characterized. In every case a scaling law is numerically obtained in which the percentage of the trapped trajectories as a function of the relativistic parameter β and the energy is obtained. Our work could be useful in the context of charged particles which eventually can be trapped in the magnetosphere, where the analysis of these structures can be relevant.

Deep learning-based analysis of basins of attraction

This research addresses the challenge of characterizing the complexity and unpredictability of basins within various dynamical systems. The main focus is on demonstrating the efficiency of convolutional neural networks (CNNs) in this field. Conventional methods become computationally demanding when analyzing multiple basins of attraction across different parameters of dynamical systems. Our research presents an innovative approach that employs CNN architectures for this purpose, showcasing their superior performance in comparison to conventional methods. We conduct a comparative analysis of various CNN models, highlighting the effectiveness of our proposed characterization method while acknowledging the validity of prior approaches. The findings not only showcase the potential of CNNs but also emphasize their significance in advancing the exploration of diverse behaviors within dynamical systems.

Probabilistic description of dissipative chaotic scattering

We investigate the extent to which the probabilistic properties of chaotic scattering systems with dissipation can be understood from the properties of the dissipation-free system. For large energies, a fully chaotic scattering leads to an exponential decay of the survival probability P(t)∼e^{-κt}, with an escape rate κ that decreases with energy. Dissipation leads to the appearance of different finite-time regimes in P(t). We show how these different regimes can be understood for small dissipations and long times from the (effective) escape rate κ (including the nonhyperbolic regime) of the conservative system, until the energy reaches a critical value at which no escape is possible. More generally, we argue that for small dissipation and long times the surviving trajectories in the dissipative system are distributed according to the conditionally invariant measure of the conservative system at the corresponding energy. Quantitative predictions of our general theory are compared with numerical simulations in the Hénon-Heiles model.

Period-doubling bifurcations and islets of stability in two-degree-of-freedom Hamiltonian systems

In this paper, we show that the destruction of the main Kolmogorov-Arnold-Moser (KAM) islands in two-degree-of-freedom Hamiltonian systems occurs through a cascade of period-doubling bifurcations. We calculate the corresponding Feigenbaum constant and the accumulation point of the period-doubling sequence. By means of a systematic grid search on exit basin diagrams, we find the existence of numerous very small KAM islands ("islets") for values below and above the aforementioned accumulation point. We study the bifurcations involving the formation of islets and we classify them in three different types. Finally, we show that the same types of islets appear in generic two-degree-of-freedom Hamiltonian systems and in area-preserving maps.

Control of escapes in two-degree-of-freedom open Hamiltonian systems

We investigate the possibility of avoiding the escape of chaotic scattering trajectories in two-degree-of-freedom Hamiltonian systems. We develop a continuous control technique based on the introduction of coupling forces between the chaotic trajectories and some periodic orbits of the system. The main results are shown through numerical simulations, which confirm that all trajectories starting near the stable manifold of the chaotic saddle can be controlled. We also show that it is possible to jump between different unstable periodic orbits until reaching a stable periodic orbit belonging to a Kolmogorov-Arnold-Moser island.

Dependence of the escape from an axially symmetric galaxy on the energy

The escape of a particle from a dynamical system depends on the intersection between the ingoing and outgoing asymptotic trajectories to certain periodic orbits placed at the openings of the curves of zero velocity of the system. Although many efforts have been devoted to the analysis of the escape from potentials presenting multiple openings, there are still few studies on potentials with only one opening. In this article, we clarify the way in which the energy affects the escape in this type of systems, showing that, contrary to what one could expect, there are several bifurcations for certain values of the energy.

Transient chaos under coordinate transformations in relativistic systems

We use the Hénon-Heiles system as a paradigmatic model for chaotic scattering to study the Lorentz factor effects on its transient chaotic dynamics. In particular, we focus on how time dilation occurs within the scattering region by measuring the time with a clock attached to the particle. We observe that the several events of time dilation that the particle undergoes exhibit sensitivity to the initial conditions. However, the structure of the singularities appearing in the escape time function remains invariant under coordinate transformations. This occurs because the singularities are closely related to the chaotic saddle. We then demonstrate using a Cantor-like set approach that the fractal dimension of the escape time function is relativistic invariant. In order to verify this result, we compute by means of the uncertainty dimension algorithm the fractal dimensions of the escape time functions as measured with an inertial frame and a frame comoving with the particle. We conclude that, from a mathematical point of view, chaotic transient phenomena are equally predictable in any reference frame and that transient chaos is coordinate invariant.

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.

Finding normally hyperbolic invariant manifolds in two and three degrees of freedom with Hénon-Heiles-type potential

We present a method based on a Lagrangian descriptor for revealing the high-dimensional phase space structures that are of interest in nonlinear Hamiltonian systems with index-1 saddle. These phase space structures include a normally hyperbolic invariant manifold and its stable and unstable manifolds, which act as codimension-1 barriers to phase space transport. In this article, finding the invariant manifolds in high-dimensional phase space will constitute identifying coordinates on these invariant manifolds. The method of Lagrangian descriptor is demonstrated by applying to classical two and three degrees of freedom Hamiltonian systems which have implications for myriad applications in chemistry, engineering, and physics.

Chaotic transport of navigation satellites

Navigation satellites are known from numerical studies to reside in a dynamically sensitive environment, which may be of profound importance for their long-term sustainability. We derive the fundamental Hamiltonian of Global Navigation Satellite System dynamics and show analytically that near-circular trajectories lie in the neighborhood of a Normally Hyperbolic Invariant Manifold (NHIM), which is the primary source of hyperbolicity. Quasicircular orbits escape through chaotic transport, regulated by NHIM's stable and unstable manifolds, following a power-law escape time distribution P(t)∼t^-α, with α∼0.8-1.5. Our study is highly relevant for the design of satellite disposal trajectories, using manifold dynamics.

On the escape from potentials with two exit channels

The aim of this paper is to investigate the escape dynamics in a Hamiltonian system describing the motion of stars in a galaxy with two exit channels through the analysis of the successive intersections of the stable and unstable manifolds to the main unstable periodic orbits with an adequate surface of section. We describe in detail the origin of the spirals shapes of the windows through which stars escape.

Ascertaining when a basin is Wada: the merging method

Trying to imagine three regions separated by a unique boundary seems a difficult task. However, this is exactly what happens in many dynamical systems showing Wada basins. Here, we present a new perspective on the Wada property: A Wada boundary is the only one that remains unaltered under the action of merging the basins. This observation allows to develop a new method to test the Wada property, which is much faster than the previous ones. Furthermore, another major advantage of the merging method is that a detailed knowledge of the dynamical system is not required.

Escape time, relaxation, and sticky states of a softened Henon-Heiles model: Low-frequency vibrational mode effects and glass relaxation

Here we study the relaxation of a chain consisting of three masses joined by nonlinear springs and periodic conditions when the stiffness is weakened. This system, when expressed in their normal coordinates, yields a softened Henon-Heiles system. By reducing the stiffness of one low-frequency vibrational mode, a faster relaxation is enabled. This is due to a reduction of the energy barrier heights along the softened normal mode as well as for a widening of the opening channels of the energy landscape in configurational space. The relaxation is for the most part exponential, and can be explained by a simple flux equation. Yet, for some initial conditions the relaxation follows as a power law, and in many cases there is a regime change from exponential to power-law decay. We pinpoint the initial conditions for the power-law decay, finding two regions of sticky states. For such states, quasiperiodic orbits are found since almost for all components of the initial momentum orientation, the system is trapped inside two pockets of configurational space. The softened Henon-Heiles model presented here is intended as the simplest model in order to understand the interplay of rigidity, nonlinear interactions and relaxation for nonequilibrium systems such as glass-forming melts or soft matter. Our softened system can be applied to model β relaxation in glasses and suggest that local reorientational jumps can have an exponential and a nonexponential contribution for relaxation, the latter due to asymmetric molecules sticking in cages for certain orientations.

From local uncertainty to global predictions: Making predictions on fractal basins

In nonlinear systems long term dynamics is governed by the attractors present in phase space. The presence of a chaotic saddle gives rise to basins of attraction with fractal boundaries and sometimes even to Wada boundaries. These two phenomena involve extreme difficulties in the prediction of the future state of the system. However, we show here that it is possible to make statistical predictions even if we do not have any previous knowledge of the initial conditions or the time series of the system until it reaches its final state. In this work, we develop a general method to make statistical predictions in systems with fractal basins. In particular, we have applied this new method to the Duffing oscillator for a choice of parameters where the system possesses the Wada property. We have computed the statistical properties of the Duffing oscillator for different phase space resolutions, to obtain information about the global dynamics of the system. The key idea is that the fraction of initial conditions that evolve towards each attractor is scale free-which we illustrate numerically. We have also shown numerically how having partial information about the initial conditions of the system does not improve in general the predictions in the Wada regions.

Uncertainty dimension and basin entropy in relativistic chaotic scattering

Chaotic scattering is an important topic in nonlinear dynamics and chaos with applications in several fields in physics and engineering. The study of this phenomenon in relativistic systems has received little attention as compared to the Newtonian case. Here we focus our work on the study of some relevant characteristics of the exit basin topology in the relativistic Hénon-Heiles system: the uncertainty dimension, the Wada property, and the basin entropy. Our main findings for the uncertainty dimension show two different behaviors insofar as we change the relativistic parameter β, in which a crossover behavior is uncovered. This crossover point is related with the disappearance of KAM islands in phase space, which happens for velocity values above the ultrarelativistic limit, v>0.1c. This result is supported by numerical simulations and by qualitative analysis, which are in good agreement. On the other hand, the computation of the exit basins in the phase space suggests the existence of Wada basins for a range of β<0.625. We also studied the evolution of the exit basins in a quantitative manner by computing the basin entropy, which shows a maximum value for β≈0.2. This last quantity is related to the uncertainty in the prediction of the final fate of the system. Finally, our work is relevant in galactic dynamics, and it also has important implications in other topics in physics such as as in the Störmer problem, among others.

Basin entropy: a new tool to analyze uncertainty in dynamical systems

In nonlinear dynamics, basins of attraction link a given set of initial conditions to its corresponding final states. This notion appears in a broad range of applications where several outcomes are possible, which is a common situation in neuroscience, economy, astronomy, ecology and many other disciplines. Depending on the nature of the basins, prediction can be difficult even in systems that evolve under deterministic rules. From this respect, a proper classification of this unpredictability is clearly required. To address this issue, we introduce the basin entropy, a measure to quantify this uncertainty. Its application is illustrated with several paradigmatic examples that allow us to identify the ingredients that hinder the prediction of the final state. The basin entropy provides an efficient method to probe the behavior of a system when different parameters are varied. Additionally, we provide a sufficient condition for the existence of fractal basin boundaries: when the basin entropy of the boundaries is larger than log2, the basin is fractal.

Testing for Basins of Wada

Nonlinear systems often give rise to fractal boundaries in phase space, hindering predictability. When a single boundary separates three or more different basins of attraction, we say that the set of basins has theWada property and initial conditions near that boundary are even more unpredictable. Many physical systems of interest with this topological property appear in the literature. However, so far the only approach to study Wada basins has been restricted to two-dimensional phase spaces. Here we report a simple algorithm whose purpose is to look for the Wada property in a given dynamical system. Another benefit of this procedure is the possibility to classify and study intermediate situations known as partially Wada boundaries.

Effects of periodic forcing in chaotic scattering

The effects of a periodic forcing on chaotic scattering are relevant in certain situations of physical interest. We investigate the effects of the forcing amplitude and the external frequency in both the survival probability of the particles in the scattering region and the exit basins associated to phase space. We have found an exponential decay law for the survival probability of the particles in the scattering region. A resonant-like behavior is uncovered where the critical values of the frequencies ω≃1 and ω≃2 permit the particles to escape faster than for other different values. On the other hand, the computation of the exit basins in phase space reveals the existence of Wada basins depending of the frequency values. We provide some heuristic arguments that are in good agreement with the numerical results. Our results are expected to be relevant for physical phenomena such as the effect of companion galaxies, among others.

Chaotic vibrations of the duffing system with fractional damping

We examined the Duffing system with a fractional damping term. Calculating the basins of attraction, we demonstrate a broad spectrum of non-linear behaviour connected with sensitivity to the initial conditions and chaos. To quantify dynamical response of the system, we propose the statistical 0-1 test as well as the maximal Lyapunov exponent; the application of the latter encounter a few difficulties because of the memory effect due to the fractional derivative. The results are confirmed by bifurcation diagrams, phase portraits, and Poincaré sections.

Weakly noisy chaotic scattering

The effect of a weak source of noise on the chaotic scattering is relevant to situations of physical interest. We investigate how a weak source of additive uncorrelated Gaussian noise affects both the dynamics and the topology of a paradigmatic chaotic scattering problem as the one taking place in the open nonhyperbolic regime of the Hénon-Heiles Hamiltonian system. We have found long transients for the time escape distributions for critical values of the noise intensity for which the particles escape slower as compared with the noiseless case. An analysis of the survival probability of the scattering function versus the Gaussian noise intensity shows a smooth curve with one local maximum and with one local minimum which are related to those long transients and with the basin structure in phase space. On the other hand, the computation of the exit basins in phase space shows a quadratic curve for which the basin boundaries lose their fractal-like structure as noise turned on.

New developments in classical chaotic scattering

Classical chaotic scattering is a topic of fundamental interest in nonlinear physics due to the numerous existing applications in fields such as celestial mechanics, atomic and nuclear physics and fluid mechanics, among others. Many new advances in chaotic scattering have been achieved in the last few decades. This work provides a current overview of the field, where our attention has been mainly focused on the most important contributions related to the theoretical framework of chaotic scattering, the fractal dimension, the basins boundaries and new applications, among others. Numerical techniques and algorithms, as well as analytical tools used for its analysis, are also included. We also show some of the experimental setups that have been implemented to study diverse manifestations of chaotic scattering. Furthermore, new theoretical aspects such as the study of this phenomenon in time-dependent systems, different transitions and bifurcations to chaotic scattering and a classification of boundaries in different types according to symbolic dynamics are also shown. Finally, some recent progress on chaotic scattering in higher dimensions is also described.

Community resilience and decision theory challenges for catastrophic events

Extreme and catastrophic events pose challenges for normative models of risk management decision making. They invite development of new methods and principles to complement existing normative decision and risk analysis. Because such events are rare, it is difficult to learn about them from experience. They can prompt both too little concern before the fact, and too much after. Emotionally charged and vivid outcomes promote probability neglect and distort risk perceptions. Aversion to acting on uncertain probabilities saps precautionary action; moral hazard distorts incentives to take care; imperfect learning and social adaptation (e.g., herd-following, group-think) complicate forecasting and coordination of individual behaviors and undermine prediction, preparation, and insurance of catastrophic events. Such difficulties raise substantial challenges for normative decision theories prescribing how catastrophe risks should be managed. This article summarizes challenges for catastrophic hazards with uncertain or unpredictable frequencies and severities, hard-to-envision and incompletely described decision alternatives and consequences, and individual responses that influence each other. Conceptual models and examples clarify where and why new methods are needed to complement traditional normative decision theories for individuals and groups. For example, prospective and retrospective preferences for risk management alternatives may conflict; procedures for combining individual beliefs or preferences can produce collective decisions that no one favors; and individual choices or behaviors in preparing for possible disasters may have no equilibrium. Recent ideas for building "disaster-resilient" communities can complement traditional normative decision theories, helping to meet the practical need for better ways to manage risks of extreme and catastrophic events.

Fractal structures in nonlinear plasma physics

Fractal structures appear in many situations related to the dynamics of conservative as well as dissipative dynamical systems, being a manifestation of chaotic behaviour. In open area-preserving discrete dynamical systems we can find fractal structures in the form of fractal boundaries, associated to escape basins, and even possessing the more general property of Wada. Such systems appear in certain applications in plasma physics, like the magnetic field line behaviour in tokamaks with ergodic limiters. The main purpose of this paper is to show how such fractal structures have observable consequences in terms of the transport properties in the plasma edge of tokamaks, some of which have been experimentally verified. We emphasize the role of the fractal structures in the understanding of mesoscale phenomena in plasmas, such as electromagnetic turbulence.

Noisy scattering dynamics in the randomly driven Hénon-Heiles oscillator

Noisy scattering dynamics in the randomly driven Hénon-Heiles oscillator is investigated when the energy is above the threshold to permit particles to escape from the scattering region. First, some basic simulation procedures are briefly introduced and the fractal exit basins appear to be robust when the bounded noisy excitation is imposed on the oscillator. Second, several key fractal characteristics of the sample basin boundaries, such as the delay-time function and the uncertainty dimension, are estimated from which this oscillator is found to be structurally unstable against the bounded noisy excitation. Moreover, the stable and unstable manifolds of some sample chaotic invariant sets are estimated and illustrated in a special two-dimensional Poincaré section. Lastly, several previous methods are developed to identify three arbitrarily chosen noisy scattering time series of the randomly driven Hénon-Heiles oscillator, from which the quasiperiodic-dominant and the chaotic-dominant dynamical behaviors are distinguished.

Effect of noise on chaotic scattering

When noise is present in a scattering system, particles tend to escape faster from the scattering region as compared with the noiseless case. For chaotic scattering, noise can render particle-decay exponential, and the decay rate typically increases with the noise intensity. We uncover a scaling law between the exponential decay rate and the noise intensity. The finding is substantiated by a heuristic argument and numerical results from both discrete-time and continuous-time models.

Fractal dimension in dissipative chaotic scattering

The effect of weak dissipation on chaotic scattering is relevant to situations of physical interest. We investigate how the fractal dimension of the set of singularities in a scattering function varies as the system becomes progressively more dissipative. A crossover phenomenon is uncovered where the dimension decreases relatively more rapidly as a dissipation parameter is increased from zero and then exhibits a much slower rate of decrease. We provide a heuristic theory and numerical support from both discrete-time and continuous-time scattering systems to establish the generality of this phenomenon. Our result is expected to be important for physical phenomena such as the advection of inertial particles in open chaotic flows, among others.

Basin topology in dissipative chaotic scattering

Chaotic scattering in open Hamiltonian systems under weak dissipation is not only of fundamental interest but also important for problems of current concern such as the advection and transport of inertial particles in fluid flows. Previous work using discrete maps demonstrated that nonhyperbolic chaotic scattering is structurally unstable in the sense that the algebraic decay of scattering particles immediately becomes exponential in the presence of weak dissipation. Here we extend the result to continuous-time Hamiltonian systems by using the Henon-Heiles system as a prototype model. More importantly, we go beyond to investigate the basin structure of scattering dynamics. A surprising finding is that, in the common case where multiple destinations exist for scattering trajectories, Wada basin boundaries are common and they appear to be structurally stable under weak dissipation, even when other characteristics of the nonhyperbolic scattering dynamics are not. We provide numerical evidence and a geometric theory for the structural stability of the complex basin topology.

Fractional kinetic model for chaotic transport in nonintegrable Hamiltonian systems

We propose a kinetic model of transport in nonintegrable Hamiltonian systems, based on a fractional kinetic equation with spatially dependent diffusion coefficient. The diffusion coefficient is estimated from the remainder of the optimal normal form for the given region of the phase space. After partitioning the phase space into building blocks, a separate equation can be constructed for each block. Solving the kinetic equations approximately and estimating the diffusion time scales, we convolve the solutions to get the description of the macroscopic behavior. We show that, in the limit of infinitely many blocks, one can expect an approximate scaling relation between the Lyapunov time and the diffusion (or escape) time, which is either an exponential or a power law. We check our results numerically on over a dozen Hamiltonians and find a good agreement.

Limit of small exits in open Hamiltonian systems

The nature of open Hamiltonian systems is analyzed, when the size of the exits decreases and tends to zero. Fractal basins appear typically in open Hamiltonian systems, but we claim that in the limit of small exits, the invariant sets tend to fill up the whole phase space with the strong consequence that a new kind of basin appears, where the unpredictability grows indefinitely. This means that for finite, arbitrarily small accuracy, we can find uncertain basins, where any information about the future of the system is lost. This total indeterminism had only been reported in dissipative systems, in particular in the so-called intermingled riddled basins, as well as in the riddledlike basins. We show that this peculiar, behavior is a general feature of open Hamiltonian systems.

Approximating stable and unstable manifolds in experiments

We introduce a procedure to reveal invariant stable and unstable manifolds, given only experimental data. We assume a model is not available and show how coordinate delay embedding coupled with invariant phase space regions can be used to construct stable and unstable manifolds of an embedded saddle. We show that the method is able to capture the fine structure of the manifold, is independent of dimension, and is efficient relative to previous techniques.

Deterministic diffusion in flower-shaped billiards

We propose a flower-shaped billiard in order to study the irregular parameter dependence of chaotic normal diffusion. Our model is an open system consisting of periodically distributed obstacles in the shape of a flower, and it is strongly chaotic for almost all parameter values. We compute the parameter dependent diffusion coefficient of this model from computer simulations and analyze its functional form using different schemes, all generalizing the simple random walk approximation of Machta and Zwanzig. The improved methods we use are based either on heuristic higher-order corrections to the simple random walk model, on lattice gas simulation methods, or they start from a suitable Green-Kubo formula for diffusion. We show that dynamical correlations, or memory effects, are of crucial importance in reproducing the precise parameter dependence of the diffusion coefficent.