Discrete Symmetry and Stability in Hamiltonian Dynamics

In this tutorial we address the existence and stability of periodic and quasiperiodic orbits in N degree of freedom Hamiltonian systems and their connection with discrete symmetries. Of primary importance in our study are the nonlinear normal modes (NNMs), i.e periodic solutions which represent continuations of the system's linear normal modes in the nonlinear regime. We examine the existence of such solutions and discuss different methods for constructing them and studying their stability under fixed and periodic boundary conditions. In the periodic case, we employ group theoretical concepts to identify a special type of NNMs called one-dimensional "bushes". We describe how to use linear combinations such NNMs to construct s(>1)-dimensional bushes of quasiperiodic orbits, for a wide variety of Hamiltonian systems and exploit the symmetries of the linearized equations to simplify the study of their destabilization. Applying this theory to the Fermi Pasta Ulam (FPU) chain, we review a number of interesting results, which have appeared in the recent literature. We then turn to an analytical and numerical construction of quasiperiodic orbits, which does not depend on the symmetries or boundary conditions. We demonstrate that the well-known "paradox" of FPU recurrences may be explained in terms of the exponential localization of the energies Eq of NNM's excited at the low part of the frequency spectrum, i.e. q=1,2,3,.... Thus, we show that the stability of these low-dimensional manifolds called q-tori is related to the persistence or FPU recurrences at low energies. Finally, we discuss a novel approach to the stability of orbits of conservative systems, the GALIk, k=2,...,2N, by means of which one can determine accurately and efficiently the destabilization of q-tori, leading to the breakdown of recurrences and the equipartition of energy, at high values of the total energy E.Comment: 50 pages, 13 figure

Chaotic Dynamics of N-degree of Freedom Hamiltonian Systems

We investigate the connection between local and global dynamics of two N-degree of freedom Hamiltonian systems with different origins describing one-dimensional nonlinear lattices: The Fermi-Pasta-Ulam (FPU) model and a discretized version of the nonlinear Schrodinger equation related to Bose-Einstein Condensation (BEC). We study solutions starting in the vicinity of simple periodic orbits (SPOs) representing in-phase (IPM) and out-of-phase motion (OPM), which are known in closed form and whose linear stability can be analyzed exactly. Our results verify that as the energy E increases for fixed N, beyond the destabilization threshold of these orbits, all positive Lyapunov exponents exhibit a transition between two power laws, occurring at the same value of E. The destabilization energy E_c per particle goes to zero as N goes to infinity following a simple power-law. However, using SALI, a very efficient indicator we have recently introduced for distinguishing order from chaos, we find that the two Hamiltonians have very different dynamics near their stable SPOs: For example, in the case of the FPU system, as the energy increases for fixed N, the islands of stability around the OPM decrease in size, the orbit destabilizes through period-doubling bifurcation and its eigenvalues move steadily away from -1, while for the BEC model the OPM has islands around it which grow in size before it bifurcates through symmetry breaking, while its real eigenvalues return to +1 at very high energies. Still, when calculating Lyapunov spectra, we find for the OPMs of both Hamiltonians that the Lyapunov exponents decrease following an exponential law and yield extensive Kolmogorov--Sinai entropies per particle, in the thermodynamic limit of fixed energy density E/N with E and N arbitrarily large.Comment: 29 pages, 10 figures, published at International Journal of Bifurcation and Chaos (IJBC

Computation and stability of periodic orbits of nonlinear mappings

Abstract. In this paper, we first present numerical methods that allow us to compute accurately periodic orbits in high dimensional mappings and demonstrate the effectiveness of our methods by computing orbits of various stability types. We then use a terminology for the different stability types, which is perfectly suited for systems with many degrees of freedom, since it clearly reflects the configuration of the eigenvalues of the corresponding monodromy matrix on the complex plane. Studying the distribution of these eigenvalues over the points of an unstable periodic orbit, we attempt to find connections between local dynamics and the global morphology of the orbit.

© IC-SCCE NUMERICAL SOLUTION OF THE BOUSSINESQ EQUATION USING SPECTRAL METHODS AND STABILITY OF SOLITARY WAVE PROPAGATION

Abstract. We study numerically the propagation and stability properties of solitary waves (solitons) of the Boussinesq equation in one space dimension, using a combination of finite differences in time and spectral methods in space. Our schemes follow very accurately these solutions, which are given by simple closed formulas and are known to be stable under small perturbations, for small enough velocities. Studying the interaction of two such solitons, we determine in the velocity parameter plane a sharp curve beyond which they become unstable. This is achieved by applying a precise criterion, which predicts when the observed amplitude growth of the waves is caused by a dynamical instability rather than the accumulation of numerical errors.