Navigation satellites are known from numerical studies to reside in a dynamically sensitive environment, which may be of profound importance for their long-term sustainability. We derive the fundamental Hamiltonian of GNSS dynamics and show analytically that near-circular trajectories lie in the neighborhood of a Normally Hyperbolic Invariant Manifold (NHIM), which is the primary source of hyperbolicity. Quasi-circular orbits escape through chaotic transport, regulated by the NHIM's stable and unstable manifolds, following a power-law escape time distribution $P(t) \sim t^{-\alpha}$, with $\alpha \sim 0.8 - 1.5$. Our study is highly relevant for the design of satellite disposal trajectories, using manifold dynamics.

We provide an analytical approximation to the dynamics in each of the three most important low order secondary resonances (1:1, 2:1, and 3:1) bifurcating from the synchronous primary resonance in the gravitational spin-orbit problem. To this end we extend the perturbative approach introduced in Gkolias et. al. (2016), based on normal form series computations. This allows to recover analytically all non-trivial features of the phase space topology and bifurcations associated with these resonances. Applications include the characterization of spin states of irregular planetary satellites or double systems of minor bodies with irregular shapes. The key ingredients of our method are: i) the use of a detuning parameter measuring the distance from the exact resonance, and ii) an efficient scheme to `book-keep' the series terms, which allows to simultaneously treat all small parameters entering the problem. Explicit formulas are provided for each secondary resonance, yielding i) the time evolution of the spin state, ii) the form of phase portraits, iii) initial conditions and stability for periodic solutions, and iv) bifurcation diagrams associated with the periodic orbits. We give also error estimates of the method, based on analyzing the asymptotic behavior of the remainder of the normal form series.

The manifold theory of barred-spiral structure provides a dynamical mechanism explaining how spiral arms beyond the ends of galactic bars can be supported by chaotic flows extending beyond the bar's co-rotation zone. We discuss its applicability to N-body simulations of secularly evolving barred galaxies. In these simulations, we observe consecutive `incidents' of spiral activity, leading to a time-varying disc morphology. Besides disc self-excitations, we provide evidence of a newly noted excitation mechanism related to the `off-centering' effect: particles ejected in elongated orbits at major incidents cause the disc center-of-mass to recoil and be set in a wobble-type orbit with respect to the halo center of mass. The time-dependent m=1 perturbation on the disc by the above mechanism correlates with the excitation of new incidents of non-axisymmetric activity beyond the bar. At every new excitation, the manifolds act as dynamical avenues attracting particles which are directed far from corotation along chaotic orbits. The fact that the manifolds evolve morphologically in time, due to varying non-axisymmetric perturbations, allows to reconcile manifolds with the presence of multiple patterns and frequencies in the disc. We find a time-oscillating pattern speed profile $\Omega_p(R)$ at distances R between the bar's corotation, at resonance with the succession of minima and maxima of the non-axisymmetric activity beyond the bar. Finally, we discuss disc thermalization, i.e., the evolution of the disc velocity dispersion profile and its connection with disc responsiveness to manifold spirals.

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.

In a realistic scenario, the evolution of the rotational dynamics of a celestial or artificial body is subject to dissipative effects. Time-varying non-conservative forces can be due to, for example, a variation of the moments of inertia or to tidal interactions. In this work, we consider a simplified model describing the rotational dynamics, known as the spin-orbit problem, where we assume that the orbital motion is provided by a fixed Keplerian ellipse. We consider different examples in which a non-conservative force acts on the model and we propose an analytical method, which reduces the system to a Hamiltonian framework. In particular, we compute a time parametrisation in a series form, which allows us to transform the original system into a Hamiltonian one. We also provide applications of our method to study the rotational motion of a body with time-varying moments of inertia, e.g. an artificial satellite with flexible components, as well as subject to a tidal torque depending linearly on the velocity.

The increasing threat raised by space debris led to the development of different mathematical models and approaches to investigate the dynamics of small particles orbiting around the Earth. Such models and methods strongly depend on the altitude of the objects above Earth's surface, since the strength of the different forces acting on an Earth orbiting object (geopotential, atmospheric drag, lunar and solar attractions, solar radiation pressure, etc.) varies with the altitude of the debris. In this review, our focus is on presenting different analytical and numerical approaches employed in modern studies of the space debris problem. We start by considering a model including the geopotential, solar and lunar gravitational forces and the solar radiation pressure. We summarize the equations of motion using different formalisms: Cartesian coordinates, Hamiltonian formulation using Delaunay and epicyclic variables, Milankovitch elements. Some of these methods lead in a straightforward way to the analysis of resonant motions. In particular, we review results found recently about the dynamics near tesseral, secular and semi-secular resonances. As an application of the above methods, we proceed to analyze a timely subject namely the possible causes for the onset of chaos in space debris dynamics. Precisely, we discuss the phenomenon of overlapping of resonances, the effect of a large area-to-mass ratio, the influence of lunisolar secular resonances. We conclude with a short discussion about the effect of the dissipation due to the atmospheric drag and we provide a list of minor effects, which could influence the dynamics of space debris.

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we i) give theoretical estimates for the long-term behavior of such orbits, and ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average `forced eccentricity', as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant `Laplace plane'. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.

Using N-body simulations we study the structures induced on a galactic disc by repeated flybys of a companion in decaying eccentric orbit around the disc. Our system is composed by a stellar disc, bulge and live dark matter halo, and we study the system's dynamical response to a sequence of a companion's flybys, when we vary i) the disc's temperature (parameterized by Toomre's Q-parameter) and ii) the companion's mass and initial orbit. We use a new 3D Cartesian grid code: MAIN (Mesh-adaptive Approximate Inverse N-body solver). The main features of MAIN are reviewed, with emphasis on the use of a new Symmetric Factored Approximate Sparse Inverse (SFASI) matrix in conjunction with the multigrid method that allows the efficient solution of Poisson's equation in three space variables. We find that: i) companions need to be assigned initial masses in a rather narrow window of values in order to produce significant and more long-standing non-axisymmetric structures (bars and spirals) in the main galaxy's disc by the repeated flyby mechanism. ii) a crucial phenomenon is the antagonism between companion-excited and self-excited modes on the disc. Values of $Q >1.5$ are needed in order to allow for the growth of the companion-excited modes to prevail over the the growth of the disc's self-excited modes. iii) We give evidence that the companion-induced spiral structure is best represented by a density wave with pattern speed nearly constant in a region extending from the ILR to a radius close to, but inside, corotation.

A number of studies, referring to the observed Trojan asteroids of various planets in our Solar System, or to hypothetical Trojan bodies in extrasolar planetary systems, have emphasized the importance of so-called secondary resonances in the problem of the long term stability of Trojan motions. Such resonances describe commensurabilities between the fast, synodic, and secular frequency of the Trojan body, and, possibly, additional slow frequencies produced by more than one perturbing bodies. The presence of secondary resonances sculpts the dynamical structure of the phase space. Hence, identifying their location is a relevant task for theoretical studies. In the present paper we combine the methods introduced in two recent papers (Paez & Efthymiopoulos, 2015, Paez & Locatelli, 2015) in order to analytically predict the location of secondary resonances in the Trojan problem (SEE FILE FOR COMPLETE ABSTRACT)

We study the resonant dynamics in a simple one degree of freedom, time dependent Hamiltonian model describing spin-orbit interactions. The equations of motion admit periodic solutions associated with resonant motions, the most important being the synchronous one in which most evolved satellites of the Solar system, including the Moon, are observed. Such primary resonances can be surrounded by a chain of smaller islands which one refers to as secondary resonances. Here, we propose a novel canonical normalization procedure allowing to obtain a higher order normal form, by which we obtain analytical results on the stability of the primary resonances as well as on the bifurcation thresholds of the secondary resonances. The procedure makes use of the expansion in a parameter, called the detuning, measuring the shift from the exact secondary resonance. Also, we implement the so-called `book-keeping' method, i.e., the introduction of a suitable separation of the terms in orders of smallness in the normal form construction, which deals simultaneously with all the small parameters of the problem. Our analytical computation of the bifurcation curves is in excellent agreement with the results obtained by a numerical integration of the equations of motion, thus providing relevant information on the parameter regions where satellites can be found in a stable configuration.

In the framework of the ERTBP, we study an example of the influence of secondary resonances over the long term stability of Trojan motions. By the integration of ensembles of orbits, we find various types of chaotic diffusion, slow and fast. We show that the distribution of escape times is bi-modular, corresponding to two populations of short and long escape times. The objects with long escape times produce a power-law tail in the distribution.

Using both numerical and analytical approaches, we demonstrate the existence of an effective power-law relation $L\propto m^p$ between the mean Lyapunov exponent $L$ of stellar orbits chaotically scattered by a supermassive black hole in the center of a galaxy and the mass parameter $m$, i.e. ratio of the mass of the black hole over the mass of the galaxy. The exponent $p$ is found numerically to obtain values in the range $p \approx 0.3$--$0.5$. We propose a theoretical interpretation of these exponents, based on estimates of local `stretching numbers', i.e. local Lyapunov exponents at successive transits of the orbits through the black hole's sphere of influence. We thus predict $p=2/3-q$ with $q\approx 0.1$--$0.2$. Our basic model refers to elliptical galaxy models with a central core. However, we find numerically that an effective power law scaling of $L$ with $m$ holds also in models with central cusp, beyond a mass scale up to which chaos is dominated by the influence of the cusp itself. We finally show numerically that an analogous law exists also in disc galaxies with rotating bars. In the latter case, chaotic scattering by the black hole affects mainly populations of thick tube-like orbits surrounding some low-order branches of the $x_1$ family of periodic orbits, as well as its bifurcations at low-order resonances, mainly the Inner Lindbland resonance and the 4/1 resonance. Implications of the correlations between $L$ and $m$ to determining the rate of secular evolution of galaxies are discussed.

We investigate the dynamics of small trojan exoplanets in domains of secondary resonances within the tadpole domain of motion. We consider the limit of a massless trojan companion of a giant planet. Without other planets, this is a case of the elliptic restricted three body problem (ERTBP). The presence of more planets (the restricted multi-planet problem, RMPP) induces new direct and indirect secular effects on the trojan's dynamics. In the theoretical part of this paper, we develop a Hamiltonian formalism in action-angle variables, which allows to treat in a unified way resonant dynamics and secular effects on the trojan body in both the ERTBP or the RMPP. Our formalism leads to a decomposition of the Hamiltonian in two parts, $H=H_b+H_{sec}$. $H_b$, called the basic model, describes resonant dynamics in the short-period (epicyclic) and synodic (libration) degrees of freedom. $H_{sec}$ contains only terms depending on slow (secular) angles. $H_b$ is formally identical in the ERTBP and the RMPP, apart from a re-definition of angular variables. An important physical consequence is that the slow chaotic diffusion proceeds in both the ERTBP and the RMPP by a qualitatively similar dynamical mechanism better approximated by the paradigm of `modulational diffusion'. In the numerical part, we focus on the ERTBP for making a numerical demonstration of the chaotic diffusion process along resonances. Using color stability maps, we provide a survey of the resonant web for characteristic mass parameters of the primary, in which the secondary resonances from 1:5 to 1:12 (ratio of the short over the synodic period) and their resonant multiplets appear. We give numerical examples of diffusion of weakly chaotic orbits in the resonant web. We make a statistics of the escaping times in the resonant domain, and find power-law tails of the distribution of escaping times for slowly diffusing chaotic orbits.

We consider the dynamics of a small trojan companion of a hypothetical giant exoplanet under the secular perturbations of additional planets. By a suitable choice of action-angle variables, the problem is amenable to the study of the slow modulation, induced by secular perturbations, to the dynamics of an otherwise called `basic' Hamiltonian model of two degrees of freedom (planar case). We present this Hamiltonian decomposition, which implies that the slow chaotic diffusion at resonances is best described by the paradigm of modulational diffusion.

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.

We model nanoflare heating of extrapolated active-region coronal loops via the acceleration of electrons and protons in Harris-type current sheets. The kinetic energy of the accelerated particles is estimated using semi-analytical and test-particle-tracing approaches. Vector magnetograms and photospheric Doppler velocity maps of NOAA active region 09114, recorded by the Imaging Vector Magnetograph (IVM), were used for this analysis. A current-free field extrapolation of the active-region corona was first constructed. The corresponding Poynting fluxes at the footpoints of 5000 extrapolated coronal loops were then calculated. Assuming that reconnecting current sheets develop along these loops, we utilized previous results to estimate the kinetic-energy gain of the accelerated particles and we related this energy to nanoflare heating and macroscopic loop characteristics. Kinetic energies of 0.1 to 8 keV (for electrons) and 0.3 to 470 keV (for protons) were found to cause heating rates ranging from $10^{-6}$ to 1 $\mathrm{erg\, s^{-1} cm^{-3}}$. Hydrodynamic simulations show that such heating rates can sustain plasma in coronal conditions inside the loops and generate plasma thermal distributions which are consistent with active region observations. We concluded the analysis by computing the form of X-ray spectra generated by the accelerated electrons using the thick target approach that were found to be in agreement with observed X-ray spectra, thus supporting the plausibility of our nanoflare-heating scenario.

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.

The acceleration of charged particles (electrons and protons) in flaring solar active regions is analyzed by numerical experiments. The acceleration is modelled as a stochastic process taking place by the interaction of the particles with local magnetic reconnection sites via multiple steps. Two types of local reconnecting topologies are studied: the Harris-type and the X-point. A formula for the maximum kinetic energy gain in a Harris-type current sheet, found in a previous work of ours, fits well the numerical data for a single step of the process. A generalization is then given approximating the kinetic energy gain through an X-point. In the case of the multiple step process, in both topologies the particles' kinetic energy distribution is found to acquire a practically invariant form after a small number of steps. This tendency is interpreted theoretically. Other characteristics of the acceleration process are given, such as the mean acceleration time and the pitch angle distributions of the particles.

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.

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.