It is known that the asymptotic invariant manifolds around an unstable periodic orbit in conservative systems can be represented by convergent series (Cherry 1926, Moser 1956, 1958, Giorgilli 2001). The unstable and stable manifolds intersect at an infinity of homoclinic points, generating a complicated homoclinic tangle. In the case of simple mappings it was found (Da Silva Ritter et al. 1987) that the domain of convergence of the formal series extends to infinity along the invariant manifolds. This allows in practice to study the homoclinic tangle using only series. However in the case of Hamiltonian systems, or mappings with a finite analyticity domain,the convergence of the series along the asymptotic manifolds is also finite. Here, we provide numerical indications that the convergence does not reach any homoclinic points. We discuss in detail the convergence problem in various cases and we find the degree of approximation of the analytical invariant manifolds to the real (numerical) manifolds as i) the order of truncation of the series increases, and ii) we use higher numerical precision in computing the coefficients of the series. Then we introduce a new method of series composition, by using action-angle variables, that allows the calculation of the asymptotic manifolds up to an a arbitrarily large extent. This is the first case of an analytic development that allows the computation of the invariant manifolds and their intersections in a Hamiltonian system for an extent long enough to allow the study of homoclinic chaos by analytical means.

In this article we discuss a possible way of testing the Kerr black hole hypothesis by taking advantage of phenomena correlated with chaotic motion in the final stages of an accretion disk around a bumpy black hole. We anticipate that these phenomena should have an imprint in the electromagnetic spectrum coming from the accretion disk.

We study the dynamics in the neighborhood of simple and double unstable periodic orbits in a rotating 3D autonomous Hamiltonian system of galactic type. In order to visualize the four dimensional spaces of section we use the method of color and rotation. We investigate the structure of the invariant manifolds that we found in the neighborhood of simple and double unstable periodic orbits in the 4D spaces of section. We consider orbits in the neighborhood of the families x1v2, belonging to the x1 tree, and the z-axis (the rotational axis of our system). Close to the transition points from stability to simple instability, in the neighborhood of the bifurcated simple unstable x1v2 periodic orbits we encounter the phenomenon of stickiness as the asymptotic curves of the unstable manifold surround regions of the phase space occupied by rotational tori existing in the region. For larger energies, away from the bifurcating point, the consequents of the chaotic orbits form clouds of points with mixing of color in their 4D representations. In the case of double instability, close to x1v2 orbits, we find clouds of points in the four dimensional spaces of section. However, in some cases of double unstable periodic orbits belonging to the z-axis family we can visualize the associated unstable eigensurface. Chaotic orbits close to the periodic orbit remain sticky to this surface for long times (of the order of a Hubble time or more). Among the orbits we studied we found those close to the double unstable orbits of the x1v2 family having the largest diffusion speed.

We study the role of asymptotic curves in supporting the spiral structure of a N-body model simulating a barred spiral galaxy. Chaotic orbits with initial conditions on the unstable asymptotic curves of the main unstable periodic orbits follow the shape of the periodic orbits for an initial interval of time and then they are diffused outwards supporting the spiral structure of the galaxy. Chaotic orbits having small deviations from the unstable periodic orbits, stay close and along the corresponding unstable asymptotic manifolds, supporting the spiral structure for more than 10 rotations of the bar. Chaotic orbits of different Jacobi constants support different parts of the spiral structure. We also study the diffusion rate of chaotic orbits outwards and find that chaotic orbits that support the outer parts of the galaxy are diffused outwards more slowly than the orbits supporting the inner parts of the spiral structure.

We study the orbital behavior at the neighborhood of complex unstable periodic orbits in a 3D autonomous Hamiltonian system of galactic type. At a transition of a family of periodic orbits from stability to complex instability (also known as Hamiltonian Hopf Bifurcation) the four eigenvalues of the stable periodic orbits move out of the unit circle. Then the periodic orbits become complex unstable. In this paper we first integrate initial conditions close to the ones of a complex unstable periodic orbit, which is close to the transition point. Then, we plot the consequents of the corresponding orbit in a 4D surface of section. To visualize this surface of section we use the method of color and rotation [Patsis and Zachilas 1994]. We find that the consequents are contained in 2D "confined tori". Then, we investigate the structure of the phase space in the neighborhood of complex unstable periodic orbits, which are further away from the transition point. In these cases we observe clouds of points in the 4D surfaces of section. The transition between the two types of orbital behavior is abrupt.

We study the formation of the spiral structure of barred spiral galaxies, using an $N$-body model. The evolution of this $N$-body model in the adiabatic approximation maintains a strong spiral pattern for more than 10 bar rotations. We find that this longevity of the spiral arms is mainly due to the phenomenon of stickiness of chaotic orbits close to the unstable asymptotic manifolds originated from the main unstable periodic orbits, both inside and outside corotation. The stickiness along the manifolds corresponding to different energy levels supports parts of the spiral structure. The loci of the disc velocity minima (where the particles spend most of their time, in the configuration space) reveal the density maxima and therefore the main morphological structures of the system. We study the relation of these loci with those of the apocentres and pericentres at different energy levels. The diffusion of the sticky chaotic orbits outwards is slow and depends on the initial conditions and the corresponding Jacobi constant.

The unstable invariant manifolds of the short-period family of periodic orbits around the unstable Lagrangian points $L_1$ and $L_2$ of a barred galaxy define loci in the configuration space which take the form of a trailing spiral pattern. In the present paper we investigate this association in the case of the self-consistent models of Kaufmann & Contopoulos (1996) which provide an approximation of real barred-spiral galaxies. We also examine the relation of `response' models of barred-spiral galaxies with the theory of the invariant manifolds. Our main results are the following: The invariant manifolds yield the correct form of the imposed spiral pattern provided that their calculation is done with the spiral potential term turned on. We provide a theoretical model explaining the form of the invariant manifolds that supports the spiral structure. The azimuthal displacement of the Lagrangian points with respect to the bar's major axis is a crucial parameter in this modeling. When this is taken into account, the manifolds necessarily develop in a spiral-like domain of the configuration space, delimited from below by the boundary of a banana-like non-permitted domain, and from above either by rotational KAM tori or by cantori forming a stickiness zone. We construct `spiral response' models on the basis of the theory of the invariant manifolds and examine the connection of the latter to the `response' models (Patsis 2006) used to fit real barred-spiral galaxies, explaining how are the manifolds related to a number of morphological features seen in such models.

We distinguish two types of stickiness in systems of two degrees of freedom (a) stickiness around an island of stability and (b) stickiness in chaos, along the unstable asymptotic curves of unstable periodic orbits. We studied these effects in the standard map with a rather large nonlinearity K=5, and we emphasized the role of the asymptotic curves U, S from the central orbit O and the asymptotic curves U+U-S+S- from the simplest unstable orbit around the island O1. We calculated the escape times (initial stickiness times) for many initial points outside but close to the island O1. The lines that separate the regions of the fast from the slow escape time follow the shape of the asymptotic curves S+,S-. We explained this phenomenon by noting that lines close to S+ on its inner side (closer to O1) approach a point of the orbit 4/9, say P1, and then follow the oscillations of the asymptotic curve U+, and escape after a rather long time, while the curves outside S+ after their approach to P1 follow the shape of the asymptotic curves U- and escape fast into the chaotic sea. All these curves return near the original arcs of U+,U- and contribute to the overall stickiness close to U+,U-. The isodensity curves follow the shape of the curves U+,U- and the maxima of density are along U+,U-. For a rather long time the stickiness effects along U+,U- are very pronounced. However after much longer times (about 1000 iterations) the overall stickiness effects are reduced and the distribution of points in the chaotic sea outside the islands tends to be uniform.

We study the dynamics of the closed scalar field FRW cosmological models in the framework of the so called Unified Dark Matter (UDM) scenario. Performing a theoretical as well as a numerical analysis we find that there is a strong indication of chaos in agreement with previous studies. We find that a positive value of the spatial curvature is essential for the appearance of chaoticity, though the Lyapunov number seems to be independent of the curvature value. Models that are close to flat exhibit a chaotic behavior after a long time while pure flat models do not exhibit any chaos. Moreover, we find that some of the semiflat models in the UDM scenario exhibit similar dynamical behavior with the Lambda cosmology despite their chaoticity. Finally, we compare the measured evolution of the Hubble parameter derived from the differential ages of passively evolving galaxies with that expected in the semiflat unified scalar field cosmology. Based on a specific set of initial conditions we find that the UDM scalar field model matches well the observational data.

We study the orbital structure of a self-consistent N-body equilibrium configuration of a barred galaxy constructed from cosmological initial conditions. The value of its spin parameter L is near the observed value of our Galaxy L=0.22. We classify the orbits in regular and chaotic using a combination of two different methods and find 60% of them to be chaotic. We examine the phase space using projections of the 4D surfaces of section for test particles as well as for real N-body particles. The real particles are not uniformly distributed in the whole phase space but they avoid orbits that do not support the bar. We use frequency analysis for the regular orbits as well as for the chaotic ones to classify certain types of orbits of our self-consistent system. We find the main resonant orbits and their statistical weight in supporting the shape of the bar and we emphasize the role of weakly chaotic orbits in supporting the boxiness at the end of the bar.

This paper summarises an investigation of the statistical properties of orbits escaping from three different two-degree-of-freedom Hamiltonian systems which exhibit global stochasticity. Each H=H_{0}+eH', with H_{0} integrable and eH' a nonintegrable correction, not necessarily small. For e below a critical e_{0} escapes are impossible energetically. For somewhat higher values, escape is allowed energetically but many orbits never escape and the escape probability P for a generic orbit ensemble decays exponentially. At or near a critical e_{1}>e_{0}, there is an abrupt qualitative change in behaviour. Above e_{1}, P typically exhibits (1) a rapid evolution towards a nonzero P_{0}(e) followed by (2) a much slower subsequent decay towards zero which, in at least one case, is well fit by a power law P=const x t^{m}, with m=0.35-0.40. In all three cases, P_{0} and the time T required to converge towards P_{0} scales in e-e_{1}, i.e., P_{0}=const x (e-e_{1})^{a} and T=const x (e-e_{1})^{b}, and T also scales in the size r of the region sampled for initial conditions, i.e., T=const x r^{-d}. To within statistical uncertainties, the best fit values of the critical exponents are the same for all three potentials, namely: a=0.5, b=0.4, and d=0.1, and satisfy a-b-d=0. The transitional behaviour observed near e_{1} is attributed to the breakdown of some especially significant KAM tori or cantori. The power law behaviour at late times is interpreted as reflecting intrinsic diffusion of chaotic orbits through cantori surrounding islands of regular orbits.