A numerical method for the approximation of stable and unstable manifolds of microscopic simulators

We address a numerical methodology for the approximation of coarse-grained stable and unstable manifolds of saddle equilibria/stationary states of multiscale/stochastic systems for which a macroscopic description does not exist analytically in a closed form. Thus, the underlying hypothesis is that we have a detailed microscopic simulator (Monte Carlo, molecular dynamics, agent-based model etc.) that describes the dynamics of the subunits of a complex system (or a black-box large-scale simulator) but we do not have explicitly available a dynamical model in a closed form that describes the emergent coarse-grained/macroscopic dynamics. Our numerical scheme is based on the equation-free multiscale framework, and it is a three-tier procedure including (a) the convergence on the coarse-grained saddle equilibrium, (b) its coarse-grained stability analysis, and (c) the approximation of the local invariant stable and unstable manifolds; the later task is achieved by the numerical solution of a set of homological/functional equations for the coefficients of a polynomial approximation of the manifolds.

Parity and time reversal elucidate both decision-making in empirical models and attractor scaling in critical Boolean networks

We present new applications of parity inversion and time reversal to the emergence of complex behavior from simple dynamical rules in stochastic discrete models. Our parity-based encoding of causal relationships and time-reversal construction efficiently reveal discrete analogs of stable and unstable manifolds. We demonstrate their predictive power by studying decision-making in systems biology and statistical physics models. These applications underpin a novel attractor identification algorithm implemented for Boolean networks under stochastic dynamics. Its speed enables resolving a long-standing open question of how attractor count in critical random Boolean networks scales with network size and whether the scaling matches biological observations. Via 80-fold improvement in probed network size (N = 16,384), we find the unexpectedly low scaling exponent of 0.12 ± 0.05, approximately one-tenth the analytical upper bound. We demonstrate a general principle: A system’s relationship to its time reversal and state-space inversion constrains its repertoire of emergent behaviors.

Topological and geometric hyperbolicity criteria for polynomial automorphisms of

We prove that uniform hyperbolicity is invariant under topological conjugacy for dissipative polynomial automorphisms of $\mathbb {C}^2$ . Along the way we also show that a sufficient condition for hyperbolicity is that local stable and unstable manifolds of saddle points have uniform geometry.

Bogdanov–Takens bifurcation of a Holling IV prey–predator model with constant-effort harvesting

A prey–predator model with constant-effort harvesting on the prey and predators is investigated in this paper. First, we discuss the number and type of the equilibria by analyzing the equations of equilibria and the distribution of eigenvalues. Second, with the rescaled harvesting efforts as bifurcation parameters, a subcritical Hopf bifurcation is exhibited near the multiple focus and a Bogdanov–Takens bifurcation is also displayed near the BT singularity by analyzing the versal unfolding of the model. With the variation of bifurcation parameters, the system shows multi-stable structure, and the attractive domains for different attractors are constituted by the stable and unstable manifolds of saddles and the limit cycles bifurcated from Hopf and Bogdanov–Takens bifurcations. Finally, a cusp point and two generalized Hopf points are found on the saddle-node bifurcation curve and the Hopf bifurcation curves, respectively. Several phase diagrams for parameters near one of the generalized Hopf points are exhibited through the generalized Hopf bifurcation.

Locating and Stabilizing Unstable Periodic Orbits Embedded in the Horseshoe Map

Based on the theory of symbolic dynamical systems, we propose a novel computation method to locate and stabilize the unstable periodic points (UPPs) in a two-dimensional dynamical system with a Smale horseshoe. This method directly implies a new framework for controlling chaos. By introducing the subset based correspondence between a planar dynamical system and a symbolic dynamical system, we locate regions sectioned by stable and unstable manifolds comprehensively and identify the specified region containing a UPP with the particular period. Then Newton’s method compensates the accurate location of the UPP with the regional information as an initial estimation. On the other hand, the external force control (EFC) is known as an effective method to stabilize the UPPs. By applying the EFC to the located UPPs, robust controlling chaos is realized. In this framework, we never use ad hoc approaches to find target UPPs in the given chaotic set. Moreover, the method can stabilize UPPs with the specified period regardless of the situation where the targeted chaotic set is attractive. As illustrative numerical experiments, we locate and stabilize UPPs and the corresponding unstable periodic orbits in a horseshoe structure of the Duffing equation. In spite of the strong instability of UPPs, the controlled orbit is robust and the control input retains being tiny in magnitude.

The Role of Time-Dependent Phase Space Structures in Reaction Dynamics and the No-Recrossing Property of Dividing Surfaces

We analyze benchmark models for reaction dynamics associated with a time-dependent saddle point. Our model allows us to incorporate time dependence of a general form, subject to an exponential growth restriction. Under these conditions, we analytically compute the time-dependent normally hyperbolic invariant manifold; its time-dependent stable and unstable manifolds; and a time-dependent dividing surface that has the no-recrossing property. Consideration of the time dependence of these phase space structures is necessary in order to precisely capture reacting and nonreacting trajectories. Moreover, we show that a time-dependent dividing surface is necessary in order to eliminate recrossing in the time-dependent setting. In other words, if the dividing surface is not time-dependent, recrossing may occur.

Dynamical Analysis and Simulation of a New Lorenz-Like Chaotic System

This work presents and investigates a new chaotic system with eight terms. By numerical simulation, the two-scroll chaotic attractor is found for some certain parameters. And, by theoretical analysis, we discuss the dynamical behavior of the new-type Lorenz-like chaotic system. Firstly, the local dynamical properties, such as the distribution and the local stability of all equilibrium points, the local stable and unstable manifolds, and the Hopf bifurcations, are all revealed as the parameters varying in the space of parameters. Secondly, by applying the way of Poincaré compactification in ℝ 3 , the dynamics at infinity are clearly analyzed. Thirdly, combining the dynamics at finity and those at infinity, the global dynamical behaviors are formulated. Especially, we have proved the existence of the infinite heteroclinic orbits. Furthermore, all obtained theoretical results in this paper are further verified by numerical simulations.

Rendezvous in Cis-Lunar Space near Rectilinear Halo Orbit: Dynamics and Control Issues

The paper presents the development of a fully-safe, automatic rendezvous strategy between a passive vehicle and an active one orbiting around the Earth–Moon L2 Lagrangian point. This is one of the critical phases of future missions to permanently return to the Moon, which are of interest to the majority of space organizations. The first step in the study is the derivation of a suitable full 6-DOF relative motion model in the Local Vertical Local Horizontal reference frame, most suitable for the design of the guidance. The main dynamic model is approximated using both the elliptic and circular three-body motion, due to the contribution of Earth and Moon gravity. A rather detailed set of sensors and actuator dynamics was also implemented in order to ensure the reliability of the guidance algorithms. The selection of guidance and control is presented, and evaluated using a sample scenario as described by ESA’s HERACLES program. The safety, in particular the passive safety, concept is introduced and different techniques to guarantee it are discussed that exploit the ideas of stable and unstable manifolds to intrinsically guarantee some properties at each hold-point, in which the rendezvous trajectory is divided. Finally, the rendezvous dynamics are validated using available Ephemeris models in order to verify the validity of the results and their limitations for future more detailed design.

Different types of stable periodic points of diffeomorphism of a plane with a homoclinic orbit

A diffeomorphism of the plane into itself with a fixed hyperbolic point is considered; the presence of a nontransverse homoclinic point is assumed. Stable and unstable manifolds touch each other at a homoclinic point; there are various ways of touching a stable and unstable manifold. In the works of Sh. Newhouse, L. P. Shilnikov and other authors, studied diffeomorphisms of the plane with a nontranverse homoclinic point, under the assumption that this point is a tangency point of finite order. It follows from the works of these authors that an infinite set of stable periodic points can lie in a neighborhood of a homoclinic point; the presence of such a set depends on the properties of the hyperbolic point. In this paper, it is assumed that a homoclinic point is not a point at which the tangency of a stable and unstable manifold is a tangency of finite order. Allocate a countable number of types of periodic points lying in the vicinity of a homoclinic point; points belonging to the same type are called n-pass (multi-pass), where n is a natural number. In the present paper, it is shown that if the tangency is not a tangency of finite order, the neighborhood of a nontransverse homolinic point can contain an infinite set of stable single-pass, double-pass, or three-pass periodic points with characteristic exponents separated from zero.

A Surface of Heteroclinic Connections Between Two Saddle Slow Manifolds in the Olsen Model

The Olsen model for the biochemical peroxidase-oxidase reaction has a parameter regime where one of its four variables evolves much slower than the other three. It is characterized by the existence of periodic orbits along which a large oscillation is followed by many much smaller oscillations before the process repeats. We are concerned here with a crucial ingredient for such mixed-mode oscillations (MMOs) in the Olsen model: a surface of connecting orbits that is followed closely by the MMO periodic orbit during its global, large-amplitude transition back to another onset of small oscillations. Importantly, orbits on this surface connect two one-dimensional saddle slow manifolds, which exist near curves of equilibria of the limit where the slow variable is frozen and acts as a parameter of the so-called fast subsystem. We present a numerical method, based on formulating suitable boundary value problems, to compute such a surface of connecting orbits. It involves a number of steps to compute the slow manifolds, certain submanifolds of their stable and unstable manifolds and, finally, a first connecting orbit that is then used to sweep out the surface by continuation. If it exists, such a surface of connecting orbits between two one-dimensional saddle slow manifolds is robust under parameter variations. We compute and visualize it in the Olsen model and show how this surface organizes the global return mechanism of MMO periodic orbits from the end of small oscillations back to a region of phase space where they start again.