We investigate regular and chaotic two-dimensional (2D) and three-dimensional (3D) orbits of stars in models of a galactic potential consisting in a disk, a halo and a bar, to find the origin of boxy components, which are part of the bar or (almost) the bar itself. Our models originate in snapshots of an N-body simulation, which develops a strong bar. We consider three snapshots of the simulation and for the orbital study we treat each snapshot independently, as an autonomous Hamiltonian system. The calculated corotation-to-bar-length ratios indicate that in all three cases the bar rotates slowly, while the orientation of the orbits of the main family of periodic orbits changes along its characteristic. We characterize the orbits as regular, sticky, or chaotic after integrating them for a 10 Gyr period by using the GALI$_2$ index. Boxiness in the equatorial plane is associated either with quasi-periodic orbits in the outer parts of stability islands, or with sticky orbits around them, which can be found in a large range of energies. We indicate the location of such orbits in diagrams, which include the characteristic of the main family. They are always found about the transition region from order to chaos. By perturbing such orbits in the vertical direction we find a class of 3D non-periodic orbits, which have boxy projections both in their face-on and side-on views.

We provide a concise presentation of the Smaller (SALI) and the Generalized Alignment Index (GALI) methods of chaos detection. These are efficient chaos indicators based on the evolution of two or more, initially distinct, deviation vectors from the studied orbit. After explaining the motivation behind the introduction of these indices, we sum up the behaviors they exhibit for regular and chaotic motion, as well as for stable and unstable periodic orbits, focusing mainly on finite-dimensional conservative systems: autonomous Hamiltonian models and symplectic maps. We emphasize the advantages of these methods in studying the global dynamics of a system, as well as their ability to identify regular motion on low dimensional tori. Finally we discuss several applications of these indices to problems originating from different scientific fields like celestial mechanics, galactic dynamics, accelerator physics and condensed matter physics.

We study the distinction and quantification of chaotic and regular motion in a time-dependent Hamiltonian barred galaxy model. Recently, a strong correlation was found between the strength of the bar and the presence of chaotic motion in this system, as models with relatively strong bars were shown to exhibit stronger chaotic behavior compared to those having a weaker bar component. Here, we attempt to further explore this connection by studying the interplay between chaotic and regular behavior of star orbits when the parameters of the model evolve in time. This happens for example when one introduces linear time dependence in the mass parameters of the model to mimic, in some general sense, the effect of self-consistent interactions of the actual N-body problem. We thus observe, in this simple time-dependent model also, that the increase of the bar's mass leads to an increase of the system's chaoticity. We propose a new way of using the Generalized Alignment Index (GALI) method as a reliable criterion to estimate the relative fraction of chaotic vs. regular orbits in such time-dependent potentials, which proves to be much more efficient than the computation of Lyapunov exponents. In particular, GALI is able to capture subtle changes in the nature of an orbit (or ensemble of orbits) even for relatively small time intervals, which makes it ideal for detecting dynamical transitions in time-dependent systems.

As originally formulated, the Generalized Alignment Index (GALI) method of chaos detection has so far been applied to distinguish quasiperiodic from chaotic motion in conservative nonlinear dynamical systems. In this paper we extend its realm of applicability by using it to investigate the local dynamics of periodic orbits. We show theoretically and verify numerically that for stable periodic orbits the GALIs tend to zero following particular power laws for Hamiltonian flows, while they fluctuate around non-zero values for symplectic maps. By comparison, the GALIs of unstable periodic orbits tend exponentially to zero, both for flows and maps. We also apply the GALIs for investigating the dynamics in the neighborhood of periodic orbits, and show that for chaotic solutions influenced by the homoclinic tangle of unstable periodic orbits, the GALIs can exhibit a remarkable oscillatory behavior during which their amplitudes change by many orders of magnitude. Finally, we use the GALI method to elucidate further the connection between the dynamics of Hamiltonian flows and symplectic maps. In particular, we show that, using for the computation of GALIs the components of deviation vectors orthogonal to the direction of motion, the indices of stable periodic orbits behave for flows as they do for maps.

We present a comparison of different numerical techniques for the integration of variational equations. The methods presented can be applied to any autonomous Hamiltonian system whose kinetic energy is quadratic in the generalized momenta, and whose potential is a function of the generalized positions. We apply the various techniques to the well-known H\'enon-Heiles system, and use the Smaller Alignment Index (SALI) method of chaos detection to evaluate the percentage of its chaotic orbits. The accuracy and the speed of the integration schemes in evaluating this percentage are used to investigate the numerical efficiency of the various techniques.

We present and compare different numerical schemes for the integration of the variational equations of autonomous Hamiltonian systems whose kinetic energy is quadratic in the generalized momenta and whose potential is a function of the generalized positions. We apply these techniques to Hamiltonian systems of various degrees of freedom, and investigate their efficiency in accurately reproducing well-known properties of chaos indicators like the Lyapunov Characteristic Exponents (LCEs) and the Generalized Alignment Indices (GALIs). We find that the best numerical performance is exhibited by the \textit{`tangent map (TM) method'}, a scheme based on symplectic integration techniques which proves to be optimal in speed and accuracy. According to this method, a symplectic integrator is used to approximate the solution of the Hamilton's equations of motion by the repeated action of a symplectic map $S$, while the corresponding tangent map $TS$, is used for the integration of the variational equations. A simple and systematic technique to construct $TS$ is also presented.

We present a survey of the theory of the Lyapunov Characteristic Exponents (LCEs) for dynamical systems, as well as of the numerical techniques developed for the computation of the maximal, of few and of all of them. After some historical notes on the first attempts for the numerical evaluation of LCEs, we discuss in detail the multiplicative ergodic theorem of Oseledec \cite{O_68}, which provides the theoretical basis for the computation of the LCEs. Then, we analyze the algorithm for the computation of the maximal LCE, whose value has been extensively used as an indicator of chaos, and the algorithm of the so--called `standard method', developed by Benettin et al. \cite{BGGS_80b}, for the computation of many LCEs. We also consider different discrete and continuous methods for computing the LCEs based on the QR or the singular value decomposition techniques. Although, we are mainly interested in finite--dimensional conservative systems, i. e. autonomous Hamiltonian systems and symplectic maps, we also briefly refer to the evaluation of LCEs of dissipative systems and time series. The relation of two chaos detection techniques, namely the fast Lyapunov indicator (FLI) and the generalized alignment index (GALI), to the computation of the LCEs is also discussed.

We use the Smaller ALignment Index (SALI) method of chaos detection, to study the global dynamics of conservative dynamical systems described by differential or difference equations. In particular, we consider the well--known 2 and 4--dimensional symplectic standard map, as well as an autonomous Hamiltonian system of 2 and 3 degrees of freedom describing the motion of a star in models of barred galaxies. The application of SALI helped us to compute rapidly and accurately the percentage of regular and chaotic motion for particular values of the parameters of these systems. We were also able to perform a computationally efficient determination of the dependence of these percentages on the variation of several parameters of the studied models.

We apply the Smaller ALignment Index (SALI) method to a 4--dimensional mapping of accelerator dynamics in order to distinguish rapidly, reliably and accurately between ordered and chaotic motion. The main advantage of this index is that it tends {\it exponentially} to zero in the case of chaotic orbits, while it fluctuates around non--zero values in the case of quasiperiodic trajectories. Thus, it avoids the notorious ambiguities concerning the eventual convergence of (maximum) Lyapunov exponents to (positive) non-zero values. Exploiting the different behavior of SALI in these two cases we produce phase space `charts' where regions of chaos and order are clearly identified. Evaluating the percentage of chaotic and escaping orbits as a function of the distance from the origin we are able to estimate rapidly and accurately the boundaries of the {\it dynamical aperture} of a proton beam, passing repeatedly through an array of magnetic focusing elements.

We propose the Particle Swarm Optimization (PSO) as an alternative method for locating periodic orbits in a three--dimensional (3D) model of barred galaxies. We develop an appropriate scheme that transforms the problem of finding periodic orbits into the problem of detecting global minimizers of a function, which is defined on the Poincar\'{e} Surface of Section (PSS) of the Hamiltonian system. By combining the PSO method with deflection techniques, we succeeded in tracing systematically several periodic orbits of the system. The method succeeded in tracing the initial conditions of periodic orbits in cases where Newton iterative techniques had difficulties. In particular, we found families of 2D and 3D periodic orbits associated with the inner 8:1 to 12:1 resonances, between the radial 4:1 and corotation resonances of our 3D Ferrers bar model. The main advantages of the proposed algorithm is its simplicity, its ability to work using function values solely, as well as its ability to locate many periodic orbits per run at a given Jacobian constant.

We use the Smaller Alignment Index (SALI) to distinguish rapidly and with certainty between ordered and chaotic motion in Hamiltonian flows. This distinction is based on the different behavior of the SALI for the two cases: the index fluctuates around non--zero values for ordered orbits, while it tends rapidly to zero for chaotic orbits. We present a detailed study of SALI's behavior for chaotic orbits and show that in this case the SALI exponentially converges to zero, following a time rate depending on the difference of the two largest Lyapunov exponents $\sigma_1$, $\sigma_2$ i.e. $SALI \propto e^{-(\sigma_1-\sigma_2)t}$. Exploiting the advantages of the SALI method, we demonstrate how one can rapidly identify even tiny regions of order or chaos in the phase space of Hamiltonian systems of 2 and 3 degrees of freedom.

We argue that inner rings in barred spiral galaxies are associated with specific 2D and 3D families of periodic orbits located just beyond the end of the bar. These are families located between the inner radial ultraharmonic 4:1 resonance and corotation. They are found in the upper part of a type-2 gap of the x1 characteristic, and can account for the observed ring morphologies without any help from families of the x1-tree. Due to the evolution of the stability of all these families, the ring shapes that are favored are mainly ovals, as well as polygons with `corners' on the minor axis, on the sides of the bar. On the other hand pentagonal rings, or rings of the NGC 7020 type hexagon, should be less probable. The orbits that make the rings belong in their vast majority to 3D families of periodic orbits and orbits trapped around them.

We study the conditions that favour boxiness of isodensities in the face-on views of orbital 3D models for barred galaxies. Using orbital weighted profiles we show that boxiness is in general a composite effect that appears when one considers stable orbits belonging to several families of periodic orbits. 3D orbits that are introduced due to vertical instabilities, play a crucial role in the face-on profiles and enhance their rectangularity. This happens because at the 4:1 radial resonance region we have several orbits with boxy face-on projections, instead of few rectangular-like x1 orbits, which, in a fair fraction of the models studied so far, are unstable at this region. Massive bars are characterized by rectangular-like orbits. However, we find that it is the pattern speed that affects most the elongation of the boxy feature, in the sense that fast bars are more elongated than slow ones. Boxiness in intermediate distances between the center of the model and the end of the bar can be attributed to x1v1 orbits, or to a combination of families related to the radial 3:1 resonance.

We present families, and sets of families, of periodic orbits that provide building blocks for boxy and peanut (hereafter b/p) edge-on profiles. We find cases where the b/p profile is confined to the central parts of the model and cases where a major fraction of the bar participates in this morphology. A b/p feature can be built either by 3D families associated with 3D bifurcations of the x1 family, or, in some models, even by families related with the z-axis orbits and existing over large energy intervals. The {\sf `X'} feature observed inside the boxy bulges of several edge-on galaxies can be attributed to the peaks of successive x1v1 orbits (Skokos et al. 2002a, hereafter paper I), provided their stability allows it. However in general, the x1v1 family has to overcome the obstacle of a S\ar\D\ar S transition in order to support the structure of a b/p feature. Other families that can be the backbones of b/p features are x1v4 and z3.1s. The morphology and the size of the boxy or peanut-shaped structures we find in our models is determined by the presence and stability of the families that support b/p features. The present study favours the idea that the observed edge-on profiles are the imprints of families of periodic orbits that can be found in appropriately chosen Hamiltonian systems, describing the potential of the bar.

Boxy edge-on profiles can be accounted for not only in models of barred galaxies, but also in models of normal (non-barred) galaxies. Thus, the presence of a bar is not a sine qua non condition for the appearance of this feature, as often assumed. We show that a `boxy' or a `peanut' structure in the central parts of a model is due to the presence of vertical resonances at which stable families of periodic orbits bifurcate from the planar x1 family. The orbits of these families reach in their projections on the equatorial plane a maximum distance from the center, beyond which they increase their mean radii by increasing only their deviations from the equatorial plane. The resulting orbital profiles are `stair-type' and constitute the backbone for the observed boxy structures in edge-on views of $N$-body models and, we believe, in edge-on views of disc galaxies. Since the existence of vertical resonances is independent of barred or spiral perturbations in the disc, `boxy' profiles may appear also in almost axisymmetric cases.

We investigate the orbital structure in a class of 3D models of barred galaxies. We consider different values of the pattern speed, of the strength of the bar and of the parameters of the central bulge of the galactic model. The morphology of the stable orbits in the bar region is associated with the degree of folding of the x1-characteristic. This folding is larger for lower values of the pattern speed. The elongation of rectangular-like orbits belonging to x1 and to x1-originated families depends mainly on the pattern speed. The detailed investigation of the trees of bifurcating families in the various models shows that major building blocks of 3D bars can be supplied by families initially introduced as unstable in the system, but becoming stable at another energy interval. In some models without radial and vertical 2:1 resonances we find, except for the x1 and x1-originated families, also families related to the z-axis orbits, which support the bar. Bifurcations of the x2 family can build a secondary 3D bar along the minor axis of the main bar. This is favoured in the slow rotating bar case.

In this series of papers we investigate the orbital structure of 3D models representing barred galaxies. In the present introductory paper we use a fiducial case to describe all families of periodic orbits that may play a role in the morphology of three-dimensional bars. We show that, in a 3D bar, the backbone of the orbital structure is not just the x1 family, as in 2D models, but a tree of 2D and 3D families bifurcating from x1. Besides the main tree we have also found another group of families of lesser importance around the radial 3:1 resonance. The families of this group bifurcate from x1 and influence the dynamics of the system only locally. We also find that 3D orbits elongated along the bar minor axis can be formed by bifurcations of the planar x2 family. They can support 3D bar-like structures along the minor axis of the main bar. Banana-like orbits around the stable Lagrangian points build a forest of 2D and 3D families as well. The importance of the 3D x1-tree families at the outer parts of the bar depends critically on whether they are introduced in the system as bifurcations in $z$ or in $\dot{z}$.

We study the spatial circular restricted problem of three bodies in the light of Nekhoroshev theory of stability over large time intervals. We consider in particular the Sun-Jupiter model and the Trojan asteroids in the neighborhood of the Lagrangian point $L_4$. We find a region of effective stability around the point $L_4$ such that if the initial point of an orbit is inside this region the orbit is confined in a slightly larger neighborhood of the equilibrium (in phase space) for a very long time interval. By combining analytical methods and numerical approximations we are able to prove that stability over the age of the universe is guaranteed on a realistic region, big enough to include one real asteroid. By comparing this result with the one obtained for the planar problem we see that the regions of stability in the two cases are of the same magnitude.