Resonant Normal Form Research Articles

Parametric approximations of fast close encounters of the planar three-body problem as arcs of a focus-focus dynamics

Abstract A gravitational close encounter of a small body with a planet may produce a substantial change of its orbital parameters which can be studied using the circular restricted three-body problem. In this paper we provide parametric representations of the fast close encounters with the secondary body of the planar CRTBP as arcs of non-linear focus-focus dynamics. The result is the consequence of a remarkable factorisation of the Birkhoff normal forms of the Hamiltonian of the problem represented with the Levi–Civita regularisation. The parameterisations are computed using two different sequences of Birkhoff normalisations of given order N. For each value of N, the Birkhoff normalisations and the parameters of the focus-focus dynamics are represented by polynomials whose coefficients can be computed iteratively with a computer algebra system; no quadratures, such as those needed to compute action-angle variables of resonant normal forms, are needed. We also provide some numerical demonstrations of the method for values of the mass parameter representative of the Sun–Earth and the Sun–Jupiter cases.

Free propagation of resonant waves in nonlinear dissipative metamaterials

This paper deals with the free propagation problem of resonant and close-to-resonance waves in one-dimensional lattice metamaterials endowed with nonlinearly viscoelastic resonators. The resonators' constitutive and geometric nonlinearities imply a cubic coupling with the lattice. The analytical treatment of the nonlinear wave propagation equations is carried out via a perturbation approach. In particular, after a suitable reformulation of the problem in the Hamiltonian setting, the approach relies on the well-known resonant normal form techniques from Hamiltonian perturbation theory. It is shown how the constructive features of the Lie Series formalism can be exploited in the explicit computation of the approximations of the invariant manifolds. A discussion of the metamaterial dynamic stability, either in the general or in the weak dissipation case, is presented.

On the semi-analytical construction of halo orbits and halo tubes in the elliptic restricted three-body problem

The halo orbits of the spatial circular restricted three-body problem are largely considered in space-flight dynamics to design low-energy transfers between celestial bodies. A very efficient analytical method for the computation of halo orbits, and the related transfers, has been obtained from the high-order resonant Birkhoff normal forms defined at the Lagrangian points L1−L2. In this paper, by implementing a non-linear Floquet–Birkhoff resonant normal form, we provide the definition of orbits, as well as their manifold tubes, which exist in a large order approximation of the elliptic three-body problem and generalize the halo orbits of the circular problem. Since the libration amplitude of such halo orbits is large (comparable to the distance of L1−L2 from the secondary body), and the Birkhoff normal forms are obtained through series expansions at the Lagrangian points, we provide also an error analysis of the method with respect to the orbits of the genuine elliptic restricted three-body problem.

Metastability phenomena in two-dimensional rectangular lattices with nearest-neighbour interaction

We study analytically the dynamics of two-dimensional rectangular lattices with periodic boundary conditions. We consider anisotropic initial data supported on one low-frequency Fourier mode. We show that, in the continuous approximation, the resonant normal form of the system is given by integrable PDEs. We exploit the normal form in order to prove the existence of metastability phenomena for the lattices. More precisely, we show that the energy spectrum of the normal modes attains a distribution in which the energy is shared among a packet of low-frequencies modes; such distribution remains unchanged up to the time-scale of validity of the continuous approximation.

Normal forms for the Laplace resonance

We describe a comprehensive model for systems locked in the Laplace resonance. The framework is based on the simplest possible dynamical structure provided by the Keplerian problem perturbed by the resonant coupling truncated at second order in the eccentricities. The reduced Hamiltonian, constructed by a transformation to resonant coordinates, is then submitted to a suitable ordering of the terms and to the study of its equilibria. Henceforth, resonant normal forms are computed. The main result is the identification of two different classes of equilibria. In the first class, only one kind of stable equilibrium is present: the paradigmatic case is that of the Galilean system. In the second class, three kinds of stable equilibria are possible and at least one of them is characterised by a high value of the forced eccentricity for the ‘first planet’: here, the paradigmatic case is the exo-planetary system GJ-876, in which the combination of libration centres admits triple conjunctions otherwise not possible in the Galilean case. The normal form obtained by averaging with respect to the free eccentricity oscillations describes the libration of the Laplace argument for arbitrary amplitudes and allows us to determine the libration width of the resonance. The agreement of the analytic predictions with the numerical integration of the toy models is very good.

Lissajous and halo orbits in the restricted three-body problem by normalization method

We perform an analytical study of the Lissajous and halo orbits around collinear points L1 and L2 in a spatial circular restricted three-body problem of an arbitrary value of the mass ratio. Using a canonical transformation procedure, we generate complete and resonant normal forms through reduction to center manifolds. The coefficients in the normal forms are explicitly expressed as functions of mass ratio for the first time so that one can evaluate the energy level at which bifurcation of halo orbit takes place. Another contribution of this paper is giving the analytical solutions of Lissajous and halo orbits in the initial synodic reference system through the inverse transformation of normalization. The analytical results are the series form of normalized action-angle variables, and their coefficients are also explicitly expressed as functions of mass ratio. Finally, comparison results demonstrate that the solutions for a Lissajous orbit derived through normalization method and Lindstedt–Poincare method are completely the same, while the solutions for a halo orbit derived through these two methods are different but have the roughly equal accuracy.

Computer-assisted estimates for birkhoff normal forms

Birkhoff normal forms are commonly used in order to ensure the so called "effective stability" in the neighborhood of elliptic equilibrium points for Hamiltonian systems. From a theoretical point of view, this means that the eventual diffusion can be bounded for time intervals that are exponentially large with respect to the inverse of the distance of the initial conditions from such equilibrium points. Here, we focus on an approach that is suitable for practical applications: we extend a rather classical scheme of estimates for both the Birkhoff normal forms to any finite order and their remainders. This is made for providing explicit lower bounds of the stability time (that are valid for initial conditions in a fixed open ball), by using a fully rigorous computer-assisted procedure. We apply our approach in two simple contexts that are widely studied in Celestial Mechanics: the H\'enon-Heiles model and the Circular Planar Restricted Three-Body Problem. In the latter case, we adapt our scheme of estimates for covering also the case of resonant Birkhoff normal forms and, in some concrete models about the motion of the Trojan asteroids, we show that it can be more advantageous with respect to the usual non-resonant ones.

Nonlinear Oscillations of a Spring Pendulum at the 1:1:2 Resonance by Normal Form Methods

The object of this research is the process of small oscillations of a three-dimensional elastic pendulum, tuned to a 1:1:2 resonance to vertical and horizontal oscillations. The purpose is to develop a symbolic algorithm for calculating small oscillations of the pendulum. The primary efforts are aimed at building a software package that generates formulas that approximate the movements of the pendulum with sufficient accuracy. The algorithms of this work are developed based on the resonant normal form method. The importance of the study is due to the wide range of applicability of the method of normal forms for constructing approximations of periodic and conditionally periodic local families of solutions of ODEs. For high-dimensional resonance systems, the technique is a generalization of the Poincare–Linstedt method, and for coarse systems, the Carleman linearization method. The reliability of the results is confirmed by comparisons with the results of numerical solutions. Results can be useful for specialists working at the interface of computational mathematics and continuum mechanics. Approaches and methodologies developed in this paper can be applied to solve various modeling problems.

Structure of the centre manifold of the $$L_1,L_2$$ collinear libration points in the restricted three-body problem

We present a global analysis of the center manifold of the collinear points in the circular restricted three-body problem. The phase-space structure is provided by a family of resonant 2-DOF Hamiltonian normal forms. The near 1:1 commensurability leads to the construction of a detuned Birkhoff-Gustavson normal form. The bifurcation sequences of the main orbit families are investigated by a geometric theory based on the reduction of the symmetries of the normal form, invariant under spatial mirror symmetries and time reversion. This global picture applies to any values of the mass parameter.

Normal form method in search for periodic solutions of ordinary differential equations. Case of the fourth order

Introduction:The method of resonant normal form is based on reducing a system of nonlinear ordinary differential equations to a simpler form, easier to explore. Moreover, for a number of autonomous nonlinear problems, it is possible to obtain explicit formulas which approximate numerical calculations of families of their periodic solutions. Replacing numerical calculations with their precalculated formulas leads to significant savings in computational time. Similar calculations were made earlier, but their accuracy was insufficient, and their complexity was very high.Purpose:Application of the resonant normal form method and a software package developed for these purposes to fourth-order systems in order to increase the calculation speed.Results:It has been shown that with the help of a single algorithm it is possible to study equations of high orders (4th and higher). Comparing the tabulation of the obtained formulas with the numerical solutions of the corresponding equations shows good quantitative agreement. Moreover, the speed of calculation by prepared approximating formulas is orders of magnitude greater than the numerical calculation speed. The obtained approximations can also be successfully applied to unstable solutions. For example, in the Henon — Heyles system, periodic solutions are surrounded by chaotic solutions and, when numerically integrated, the algorithms are often unstable on them.Practical relevance:The developed approach can be used in the simulation of physical and biological systems.

Hopf bifurcation for the hydrogen atom in a circularly polarized microwave field

We consider the Rydberg electron in a circularly polarized microwave field, whose dynamics is described by a Hamiltonian depending on one parameter, K > 0. The corresponding Hamiltonian system of ODE has two equilibrium points L1 (unstable for all K and energy value h(L1)) and L2 (a center for K < Kcrit and a complex saddle for K > Kcrit, with energy value h(L2)). We study the Hamiltonian-Hopf bifurcation phenomena that take place for K close to Kcrit around L2. First, a local analysis based on the computation of the integrable normal form up to a finite order is carried out and the steps for the computation of this (resonant) normal form are explained in a constructive manner. The analysis of the normal form obtained allows: to claim the type of the Hopf bifurcation –supercritical–; to study the local behavior of the electron in a neighborhood of the equilibrium L2 for the original non integrable Hamiltonian (as a perturbative approach from the integrable normal form); to obtain (approximations for) the parametrizations of the relevant invariant objects that take place due to the bifurcation (periodic orbits and invariant manifolds of L2). We compute numerically such objects and analyse not only the local picture of the dynamics close to L2, but also a global description of the dynamics and the effect of the Hopf bifurcation as well as other objects that organize the dynamics are discussed. We conclude that, for K close to Kcrit and the energy level h(L2), the Hopf bifurcation has essentially no effect on the dynamics from a physical point of view. However, for bigger values of K > Kcrit, the Hopf bifurcation has a dramatic effect: different kind of orbits coexist, mostly chaotic. Such orbits provide a ionization mechanism with several passages far from and close to L2 before ionizing. Surprisingly enough, also robust confinement regions (where the electron remains confined for ever), exist in the middle of chaotic areas.

Secondary resonances and the boundary of effective stability of Trojan motions

One of the most interesting features in the libration domain of co-orbital motions is the existence of secondary resonances. For some combinations of physical parameters, these resonances occupy a large fraction of the domain of stability and rule the dynamics within the stable tadpole region. In this work, we present an application of a recently introduced `basic Hamiltonian model' Hb for Trojan dynamics, in Paez and Efthymiopoulos (2015), Paez, Locatelli and Efthymiopoulos (2016): we show that the inner border of the secondary resonance of lowermost order, as defined by Hb, provides a good estimation of the region in phase-space for which the orbits remain regular regardless the orbital parameters of the system. The computation of this boundary is straightforward by combining a resonant normal form calculation in conjunction with an `asymmetric expansion' of the Hamiltonian around the libration points, which speeds up convergence. Applications to the determination of the effective stability domain for exoplanetary Trojans (planet-sized objects or asteroids) which may accompany giant exoplanets are discussed.

Long time stability of small-amplitude Breathers in a mixed FPU-KG model

In the limit of small couplings in the nearest neighbor interaction, and small total energy, we apply the resonant normal form result of a previous paper of ours to a finite but arbitrarily large mixed Fermi–Pasta–Ulam Klein–Gordon chain, i.e., with both linear and nonlinear terms in both the on-site and interaction potential, with periodic boundary conditions. An existence and orbital stability result for Breathers of such a normal form, which turns out to be a generalized discrete nonlinear Schrodinger model with exponentially decaying all neighbor interactions, is first proved. Exploiting such a result as an intermediate step, a long time stability theorem for the true Breathers of the KG and FPU–KG models, in the anti-continuous limit, is proven.

Approximation of small-amplitude weakly coupled oscillators by discrete nonlinear Schrödinger equations

Small-amplitude weakly coupled oscillators of the Klein–Gordon lattices are approximated by equations of the discrete nonlinear Schrödinger type. We show how to justify this approximation by two methods, which have been very popular in the recent literature. The first method relies on a priori energy estimates and multi-scale decompositions. The second method is based on a resonant normal form theorem. We show that although the two methods are different in the implementation, they produce equivalent results as the end product. We also discuss the applications of the discrete nonlinear Schrödinger equation in the context of existence and stability of breathers of the Klein–Gordon lattice.

Bifurcation Sequences in the Symmetric 1:1 Hamiltonian Resonance

We present a general review of the bifurcation sequences of periodic orbits in general position of a family of resonant Hamiltonian normal forms with nearly equal unperturbed frequencies, invariant under [Formula: see text] symmetry. The rich structure of these classical systems is investigated with geometric methods and the relation with the singularity theory approach is also highlighted. The geometric approach is the most straightforward way to obtain a general picture of the phase-space dynamics of the family as is defined by a complete subset in the space of control parameters complying with the symmetry constraint. It is shown how to find an energy-momentum map describing the phase-space structure of each member of the family, a catastrophe map that captures its global features and formal expressions for action-angle variables. Several examples, mainly taken from astrodynamics, are used as applications.

Birth of periodic and artificial halo orbits in the restricted three-body problem

We investigate the bifurcation of artificial halo orbits from the Lyapunov planar family of periodic orbits around the collinear libration points of the circular, spatial, restricted three-body problem. Beside the gravitational forces, our model includes also the effect of the Solar Radiation Pressure (SRP) and this motivates the use of the term ‘artificial’ halo orbits. Indeed, as a typical problem, one may think of a solar sail, which is characterized by a performance parameter measuring the strength of the effect of the SRP on the spacecraft.To settle the model, we determine the position of the collinear points as a function of the mass and performance parameters and the energy values at which Hill׳s surfaces allow for transit orbits between the primaries. To analyze the dynamics we use a consolidated procedure which consists in the computation of a resonant normal form, allowing the reduction to the center manifold and providing an integrable approximation of the Hamiltonian dynamical system. Finally, we compute the bifurcation thresholds of the 1:1 resonant periodic orbit families (which have the standard ‘halo’ orbits as their first member) as a function of the performance and mass parameters.The results show that SRP is indeed a relevant ingredient for new dynamical features and must definitely be considered when planning a mission of a solar sail with trajectories in the neighborhoods of collinear points.

Existence of pearled patterns in the planar functionalized Cahn–Hilliard equation

The functionalized Cahn–Hilliard (FCH) equation supports planar and circular bilayer interfaces as equilibria which may lose their stability through the pearling bifurcation: a periodic, high-frequency, in-plane modulation of the bilayer thickness. In two spatial dimensions we employ spatial dynamics and a center manifold reduction to reduce the FCH equation to an 8th order ODE system. A normal form analysis and a fixed-point-theorem argument show that the reduced system admits a degenerate 1:1 resonant normal form, from which we deduce that the onset of the pearling bifurcation coincides with the creation of a two-parameter family of pearled equilibria which are periodic in the in-plane direction and exponentially localized in the transverse direction.

Resonant normal form and asymptotic normal form behaviour in magnetic bottle Hamiltonians

We consider normal forms in ‘magnetic bottle’ type Hamiltonians of the form (second frequency ω2 equal to zero in the lowest order). Our main results are: (i) a novel method to construct the normal form in cases of resonance, and (ii) a study of the asymptotic behaviour of both the non-resonant and the resonant series. We find that, if we truncate the normal form series at order r, the series remainder in both constructions decreases with increasing r down to a minimum, and then it increases with r. The computed minimum remainder turns to be exponentially small in , where ΔE is the mirror oscillation energy, while the optimal order scales as an inverse power of ΔE. We estimate numerically the exponents associated with the optimal order and the remainder's exponential asymptotic behaviour. In the resonant case, our novel method allows to compute a ‘quasi-integral’ (i.e. truncated formal integral) valid both for each particular resonance as well as away from all resonances. We applied these results to a specific magnetic bottle Hamiltonian. The non-resonant normal form yields theoretical invariant curves on a surface of section which fit well the empirical curves away from resonances. On the other hand the resonant normal form fits very well both the invariant curves inside the islands of a particular resonance as well as the non-resonant invariant curves. Finally, we discuss how normal forms allow to compute a critical threshold for the onset of global chaos in the magnetic bottle.

On classification of the Poincaré type maps on [formula omitted

In this paper we study polynomial normal forms for maps on R3. We obtain a complete classification of all diffeomorphisms F(x) on R3 having the Poincaré type fixed point O. The simplest resonant polynomial normal forms, namely, normal forms with the possible lowest degree together with the minimal number of parameters, are explicitly given.

Qualitative and analytical results of the bifurcation thresholds to halo orbits

We study the dynamics in the neighborhood of the collinear Lagrangian points in the spatial, circular, restricted three--body problem. We consider the case in which one of the primaries is a radiating body and the other is oblate (although the latter is a minor effect). Beside having an intrinsic mathematical interest, this model is particularly suited for the description of a mission of a spacecraft (e.g., a solar sail) to an asteroid. The aim of our study is to investigate the occurrence of bifurcations to halo orbits, which take place as the energy level is varied. The estimate of the bifurcation thresholds is performed by analytical and numerical methods: we find a remarkable agreement between the two approaches. As a side result, we also evaluate the influence of the different parameters, most notably the solar radiation pressure coefficient, on the dynamical behavior of the model. To perform the analytical and numerical computations, we start by implementing a center manifold reduction. Next, we estimate the bifurcation values using qualitative techniques (e.g. Poincar\'e surfaces, frequency analysis, FLIs). Concerning the analytical approach, following \cite{CPS} we implement a resonant normal form, we transform to suitable action-angle variables and we introduce a detuning parameter measuring the displacement from the synchronous resonance. The bifurcation thresholds are then determined as series expansions in the detuning. Three concrete examples are considered and we find in all cases a very good agreement between the analytical and numerical results.