Renato Calleja

A numerically accessible criterion for the breakdown of quasi-periodic solutions and its rigorous justification

We formulate and justify rigorously a numerically efficient criterion for the computation of the analyticity breakdown of quasi-periodic solutions in symplectic maps (any dimension) and 1D statistical mechanics models. Depending on the physical interpretation of the model, the analyticity breakdown may correspond to the onset of mobility of dislocations, or of spin waves (in the 1D models) and to the onset of global transport in symplectic twist maps in 2D.The criterion proposed here is based on the blow-up of Sobolev norms of the hull functions. We prove theorems that justify the criterion. These theorems are based on an abstract implicit function theorem, which unifies several results in the literature. The proofs also lead to fast algorithms, which have been implemented and used elsewhere. The method can be adapted to other contexts.

Collision of invariant bundles of quasi-periodic attractors in the dissipative standard map

We perform a numerical study of the breakdown of hyperbolicity of quasi-periodic attractors in the dissipative standard map. In this study, we compute the quasi-periodic attractors together with their stable and tangent bundles. We observe that the loss of normal hyperbolicity comes from the collision of the stable and tangent bundles of the quasi-periodic attractor. We provide numerical evidence that, close to the breakdown, the angle between the invariant bundles has a linear behavior with respect to the perturbing parameter. This linear behavior agrees with the universal asymptotics of the general framework of breakdown of hyperbolic quasi-periodic tori in skew product systems.

Breakdown of invariant attractors for the dissipative standard map

We implement different methods for the computation of the breakdown threshold of invariant attractors in the dissipative standard mapping. A first approach is based on the computation of the Sobolev norms of the function parametrizing the solution. Then we look for the approximating periodic orbits and we analyze their stability in order to compute the critical threshold at which an invariant attractor breaks down. We also determine the domain of convergence of the dissipative standard mapping by extending the computations to the complex parameter space as well as by investigating a two-frequency model.

Resonance Phenomena in a Scalar Delay Differential Equation with Two State-Dependent Delays

We study a scalar delay differential equation (DDE) with two delayed feedback terms that depend linearly on the state. The associated constant-delay DDE, obtained by freezing the state dependence, is linear and without recurrent dynamics. With state-dependent delay terms, on the other hand, the DDE shows very complicated dynamics. To investigate this, we perform a bifurcation analysis of the system and present its bifurcation diagram in the plane of the two feedback strengths. It is organized by Hopf-Hopf bifurcation points that give rise to curves of torus bifurcation and associated two-frequency dynamics in the form of invariant tori and resonance tongues. We numerically determine the type of the Hopf-Hopf bifurcation points by computing the normal form on the center manifold; this requires the expansion of the functional defining the state-dependent DDE in a power series whose terms up to order three contain only constant delays. We implemented this expansion and the computation of the normal form coef...

Boundary-value problem formulations for computing invariant manifolds and connecting orbits in the circular restricted three body problem

We demonstrate the remarkable effectiveness of boundary value formulations coupled to numerical continuation for the computation of stable and unstable manifolds in systems of ordinary differential equations. Specifically, we consider the circular restricted three-body problem (CR3BP), which models the motion of a satellite in an Earth–Moon-like system. The CR3BP has many well-known families of periodic orbits, such as the planar Lyapunov orbits and the non-planar vertical and halo orbits. We compute the unstable manifolds of selected vertical and halo orbits, which in several cases leads to the detection of heteroclinic connections from such a periodic orbit to invariant tori. Subsequent continuation of these connecting orbits with a suitable end point condition and allowing the energy level to vary leads to the further detection of apparent homoclinic connections from the base periodic orbit to itself, or the detection of heteroclinic connections from the base periodic orbit to other periodic orbits. Some of these connecting orbits are of potential interest in space mission design.

Heteroclinic bifurcations and chaotic transport in the two-harmonic standard map

We study a two-parameter family of standard maps: the so-called two-harmonic family. In particular, we study the areas of lobes formed by the stable and unstable manifolds. Variational methods are used to find heteroclinic orbits and their action. A specific pair of heteroclinic orbits is used to define a difference in action function and to study bifurcations in the stable and unstable manifolds. Using this idea, two phenomena are studied: the change of orientation of lobes and tangential intersections of stable and unstable manifolds.

Symmetries and choreographies in families that bifurcate from the polygonal relative equilibrium of the n-body problem

We use numerical continuation and bifurcation techniques in a boundary value setting to follow Lyapunov families of periodic orbits and subsequently bifurcating families. The Lyapunov families arise from the polygonal equilibrium of n bodies in a rotating frame of reference. When the frequency of a Lyapunov orbit and the frequency of the rotating frame have a rational relationship, then the orbit is also periodic in the inertial frame. We prove that a dense set of Lyapunov orbits, with frequencies satisfying a diophantine equation, correspond to choreographies. We present a sample of the many choreographies that we have determined numerically along the Lyapunov families and along bifurcating families, namely for the cases

Periodic Solutions of a Singularly Perturbed Delay Differential Equation with Two State-Dependent Delays

n=3