In a previous paper (Voglis et al. 2006a, paper I) we demonstrated that, in a rotating galaxy with a strong bar, the unstable asymptotic manifolds of the short period family of unstable periodic orbits around the Lagrangian points L$_1$ or L$_2$ create correlations among the apocentric positions of many chaotic orbits, thus supporting a {\it spiral} structure beyond the bar. In the present paper we present evidence that the unstable manifolds of {\it all} the families of unstable periodic orbits near and beyond corotation contribute to the same phenomenon. Our results refer to a N-Body simulation, a number of drawbacks of which, as well as the reasons why these do not significantly affect the main results, are discussed. We explain the dynamical importance of the invariant manifolds as due to the fact that they produce a phenomenon of `stickiness' slowing down the rate of chaotic escape in an otherwise non-compact region of the phase space. We find a stickiness time of order 100 dynamical periods, which is sufficient to support a long-living spiral structure. Manifolds of different families become important at different ranges of values of the Jacobi constant. The projections of the manifolds of all the different families in the configuration space produce a pattern due to the `coalescence' of the invariant manifolds. This follows closely the maxima of the observed $m=2$ component near and beyond corotation. Thus, the manifolds support both the outer edge of the bar and the spiral arms.

We address the question of an appropriate choice of basis functions for the self-consistent field (SCF) method of simulation of the N-body problem. Our criterion is based on a comparison of the orbits found in N-body realizations of analytical potential-density models of triaxial galaxies, in which the potential is fitted by the SCF method using a variety of basis sets, with those of the original models. Our tests refer to maximally triaxial Dehnen gamma-models for values of $\gamma$ in the range 0<=gamma<=1. When an N-body realization of a model is fitted by the SCF method, the choice of radial basis functions affects significantly the way the potential, forces, or derivatives of the forces are reproduced, especially in the central regions of the system. We find that this results in serious discrepancies in the relative amounts of chaotic versus regular orbits, or in the distributions of the Lyapunov characteristic exponents, as found by different basis sets. Numerical tests include the Clutton-Brock and the Hernquist-Ostriker (HO) basis sets, as well as a family of numerical basis sets which are `close' to the HO basis set. The family of numerical basis sets is parametrized in terms of a quantity $\epsilon$ which appears in the kernel functions of the Sturm-Liouville (SL) equation defining each basis set. The HO basis set is the $\epsilon=0$ member of the family. We demonstrate that grid solutions of the SL equation yielding numerical basis sets introduce large errors in the variational equations of motion. We propose a quantum-mechanical method of solution of the SL equation which overcomes these errors. We finally give criteria for a choice of optimal value of $\epsilon$ and calculate the latter as a function of the value of gamma.

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.

We estimate the two-point correlation function in redshift space of the recently compiled HIPASS neutral hydrogen (HI) sources catalogue, which if modeled as a power law, $\xi(r)=(r_{0}/r)^{\gamma}$, the best-fitting parameters for the HI selected galaxies are found to be $r_{0}=3.3 \pm 0.3 h^{-1}$ Mpc with $\gamma=1.38 \pm 0.24$. Fixing the slope to its universal value $\gamma=1.8$, we obtain $r_{0}= 3.2\pm 0.2 h^{-1}$ Mpc. Comparing the measured two point correlation function with the predictions of the concordance cosmological model, we find that at the present epoch the HI selected galaxies are anti-biased with respect to the underlying matter fluctuation field with their bias value being $b_{0}\simeq 0.68$. Furthermore, dividing the HI galaxies into two richness subsamples we find that the low mass HI galaxies have a very low present bias factor ($b_{0}\simeq 0.48$), while the high mass HI galaxies trace the underlying matter distribution as the optical galaxies ($b_{0}\simeq 1$). Using our derived present-day HI galaxy bias we estimate their redshift space distortion parameter, and correct accordingly the correlation function for peculiar motions. The resulting real-space correlation length is $r^{\rm re}_{0}=1.8 \pm 0.2 h^{-1}$Mpc and $r^{\rm re}_{0}=3.9 \pm 0.6 h^{-1}$Mpc for the low and high mass HI galaxies, respectively. The low-mass HI galaxies appear to have the lowest correlation length among all extragalactic populations studied to-date. Also, we have correlated the IRAS-PSCz reconstructed density field, smoothed over scales of 5$h^{-1}$ Mpc, with the positions of the HI galaxies, to find that indeed the HI galaxies are typically found in negative overdensity regions ($\delta\rho/\rho_{\rm PSCz} \mincir 0$).

In the presence of a strong $m=2$ component in a rotating galaxy, the phase space structure near corotation is shaped to a large extent by the {\it invariant manifolds} of the short period family of unstable periodic orbits terminating at L$_1$ or L$_2$. The main effect of these manifolds is to create robust {\it phase correlations} among a number of chaotic orbits large enough to support a {\it spiral} density wave outside corotation. The phenomenon is described theoretically by soliton-like solutions of a Sine-Gordon equation. Numerical examples are given in an N-Body simulation of a barred spiral galaxy. In these examples, we demonstrate how the projection of unstable manifolds in configuration space reproduces essentially the entire observed bar-spiral pattern.

The present lecture notes are an introduction to selected topics of {\it Galactic Dynamics}. The focus is on topics that we consider more relevant to the main theme of this workshop, {\it Celestial Mechanics}. This is not intended to be a review article. In fact, any of the topics below could be the subject of a separate review. Only the main ideas and notions are introduced, as well as some important currently open problems in each topic. Some relevant results from our own research are also presented. We discuss topics related mostly to the so-called {\it ellipsoidal components} of galaxies. These are \textbf{a)} the dark halos of both elliptical and disk galaxies, \textbf{b)} the luminous matter in elliptical galaxies, and \textbf{c)} the bulges of disk galaxies. We shall only occasionally refer to the dynamics of disks, bars or spiral structure.

Dissipationless N-body models of rotating galaxies, iso-energetic to a non-rotating model, are examined as regards the mass in regular and in chaotic motion. The values of their spin parameters $\lambda$ are near the value $\lambda=0.22$ of our Galaxy. We obtain the distinction between the sets of particles moving in regular and in chaotic orbits and we show that the spatial distribution of these two sets of particles is much different. The rotating models are characterized by larger fractions of mass in chaotic motion ($\thickapprox 65%$) compared with the fraction of mass in chaotic motion in the non-rotating iso-energetic model ($\thickapprox 32%$). Furthermore, the Lyapunov numbers of the chaotic orbits in the rotating models become by about one order of magnitude larger than in the non-rotating model. Chaotic orbits are concentrated preferably in values of the Jacobi integral around the value of the effective potential at the corotation radius. We find that density waves form a central rotating bar embedded in a thin and a thick disc with exponential surface density profile. A surprising new result is that long living spiral arms are exited on the disc, composed almost completely by chaotic orbits. The bar excites an $m=2$ mode of spiral waves on the surface density of the disc, emanating from the corotation radius. These spiral waves are deformed, fade, or disappear temporarily, but they grow again re-forming a well developed spiral pattern. Spiral arms are discernible up to 20 or 30 rotations of the bar (lasting for about a Hubble time).

We study the virialization of the cosmic structures in the framework of flat cosmological models where the dark energy component plays an important role in the global dynamics of the universe. In particular, our analysis focuses on the study of the spherical matter perturbations, as the latter decouple from the background expansion and start to ``turn around'' and finally collapse. We generalize this procedure taking into account models with an equation of state which vary with time, and provide a complete formulation of the cluster virialization attempting to address the nonlinear regime of structure formation. In particular, assuming that clusters have collapsed prior to the epoch of $z_{\rm f}\simeq 1.4$, in which the most distant cluster has been found, we show that the behavior of the spherical collapse model depends on the functional form of the equation of state.