We present a new measurement of the Lense-Thirring effect on the orbits of the geodetic satellites LAGEOS, LAGEOS II and LARES. This secular precession is a general relativity effect produced by the gravitomagnetic field of the Earth generated by its rotation. The effect is a manifestation of spacetime curvature generated by mass-currents, a peculiarity of Einstein's theory of gravitation. This measurement stands out, compared to previous measurements in the same context, for its precision ($\simeq7.4\times10^{-3}$) and accuracy ($\simeq16\times10^{-3}$), i.e. for a reliable and robust evaluation of the systematic sources of error due to both gravitational and non-gravitational perturbations. For this new measurement, we have largely exploited the results of GRACE mission to significantly improve the description of the gravitational field of the Earth, by also modeling its time dependence. In this way, we strongly reduced the systematic errors due to the uncertainty in the knowledge of the Earth even zonal harmonics and, at the same time, avoided a possible bias of the final result and, consequently, of the precision of the measurement, linked to a non-reliable handling of the unmodeled and mismodeled periodic effects.

We present a global analysis of the center manifold of the collinear points in the circular restricted three-body problem. The phase-space structure is provided by a family of resonant 2-DOF Hamiltonian normal forms. The near 1:1 commensurability leads to the construction of a detuned Birkhoff-Gustavson normal form. The bifurcation sequences of the main orbit families are investigated by a geometric theory based on the reduction of the symmetries of the normal form, invariant under spatial mirror symmetries and time reversion. This global picture applies to any values of the mass parameter.

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.

We consider the three-body mean motion resonance defined by the Jovian moons Io, Europa, and Ganymede, which is commonly known as the Laplace resonance. In particular, we construct approximate models for the evolution of the librating argument over the period of 100 years, focusing on its principal amplitude and frequency, and on the observed mean motion combinations associated with the quasi-resonant interactions. First, we numerically propagated the Cartesian equations of motion of the Jovian system for the period under examination, and by comparing the results with a suitable set of ephemerides, we derived the main dynamical effects on the target quantities. Using these effects, we built an alternative Hamiltonian formulation and used the normal forms theory to locate the resonance and to compute its main amplitude and frequency. From the Cartesian model we observe that on the timescale considered and with ephemerides as initial conditions, both the librating argument and the diagnostics are well approximated by considering the mutual gravitational interactions of Jupiter and the Galilean moons (including Callisto), and the effect of Jupiter's J2 harmonic. Under the same initial conditions, the Hamiltonian formulation in which Callisto and J2 are reduced to their secular contributions achieves larger errors for the quantities above, particularly for the librating argument. By introducing appropriate resonant variables, we show that these errors can be reduced by moving in a certain action-angle phase plane, which in turn implies the necessity of a tradeoff in the selection of the initial conditions. In addition to being a good starting point for a deeper understanding of the Laplace resonance, the models and methods described are easily generalizable to different types of multi-body mean motion resonances. They are also prime tools for studying the dynamics of extrasolar systems.

We study the dynamics of the de Sitter resonance, namely the stable equilibrium configuration of the first three Galilean satellites. We clarify the relation between this family of configurations and the more general Laplace resonant states. In order to describe the dynamics around the de Sitter stable equilibrium, a one-degree of freedom Hamiltonian normal form is constructed and exploited to identify initial conditions leading to the two families. The normal form Hamiltonian is used to check the accuracy in the location of the equilibrium positions. Besides, it gives a measure of how sensitive it is with respect to the different perturbations acting on the system. By looking at the phase-plane of the normal form, we can identify a \sl Laplace-like \rm configuration, which highlights many substantial aspects of the observed one.

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.

We provide a detailed derivation of the analytical expansion of the lunar and solar disturbing functions. Although there exist several papers on this topic, many derivations contain mistakes in the final expansion or rather (just) in the proof, thereby necessitating a recasting and correction of the original derivation. In this work, we provide a self-consistent and definite form of the lunisolar expansion. We start with Kaula's expansion of the disturbing function in terms of the equatorial elements of both the perturbed and perturbing bodies. Then we give a detailed proof of Lane's expansion, in which the elements of the Moon are referred to the ecliptic plane. Using this approach the inclination of the Moon becomes nearly constant, while the argument of perihelion, the longitude of the ascending node, and the mean anomaly vary linearly with time. We make a comparison between the different expansions and we profit from such discussion to point out some mistakes in the existing literature, which might compromise the correctness of the results. As an application, we analyze the long--term motion of the highly elliptical and critically--inclined Molniya orbits subject to quadrupolar gravitational interactions. The analytical expansions presented herein are very powerful with respect to dynamical studies based on Cartesian equations, because they quickly allow for a more holistic and intuitively understandable picture of the dynamics.

Using bifurcation theory, we study the secular resonances induced by the Sun and Moon on space debris orbits around the Earth. In particular, we concentrate on a special class of secular resonances, which depend only on the debris' orbital inclination. This class is typically subdivided into three distinct types of secular resonances: those occurring at the critical inclination, those corresponding to polar orbits, and a third type resulting from a linear combination of the rates of variation of the argument of perigee and the longitude of the ascending node. The model describing the dynamics of space debris includes the effects of the geopotential, as well as the Sun's and Moon's attractions, and it is defined in terms of suitable action-angle variables. We consider the system averaged over both the mean anomaly of the debris and those of the Sun and Moon. Such multiply-averaged Hamiltonian is used to study the lunisolar resonances which depend just on the inclination. Borrowing the technique from the theory of bifurcations of Hamiltonian normal forms, we study the birth of periodic orbits and we determine the energy thresholds at which the bifurcations of lunisolar secular resonances take place. This approach gives us physically relevant information on the existence and location of the equilibria, which help us to identify stable and unstable regions in the phase space. Besides their physical interest, the study of inclination dependent resonances offers interesting insights from the dynamical point of view, since it sheds light on different phenomena related to bifurcation theory.

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.

We present a general review of the bifurcation sequences of periodic orbits in general position of a family of resonant Hamiltonian normal forms with nearly equal unperturbed frequencies, invariant under $Z_2 \times Z_2$ symmetry. The rich structure of these classical systems is investigated with geometric methods and the relation with the singularity theory approach is also highlighted. The geometric approach is the most straightforward way to obtain a general picture of the phase-space dynamics of the family as is defined by a complete subset in the space of control parameters complying with the symmetry constraint. It is shown how to find an energy-momentum map describing the phase space structure of each member of the family, a catastrophe map that captures its global features and formal expressions for action-angle variables. Several examples, mainly taken from astrodynamics, are used as applications.

We compute the normal forms for the Hamiltonian leading to the epicyclic approximations of the (perturbed) Kepler problem in the plane. The Hamiltonian setting corresponds to the dynamics in the Hill synodic system where, by means of the tidal expansion of the potential, the equations of motion take the form of perturbed harmonic oscillators in a rotating frame. In the unperturbed, purely Keplerian case, the post-epicyclic solutions produced with the normal form coincide with those obtained from the expansion of the solution of the Kepler equation. In all cases where the perturbed problem can be cast in autonomous form, the solution is easily obtained as a perturbation series. The generalization to the spatial problem and/or the non-autonomous case is straightforward.

This paper illustrates the application of Lie transform normal-form theory to the construction of the 1:2 resonant normal form corresponding to a wide class of natural Hamiltonian systems. We show how to compute the bifurcations of the main periodic orbits in a potential with double reflection symmetries. The stability analysis of the normal modes and of the periodic orbits in general position allows us to get overall informations on the phase-space structure of systems in which this resonance is dominating. As an example we apply these results to a class of models useful as galactic potentials.

We present a general analysis of the orbit structure of 2D potentials with self-similar elliptical equipotentials by applying the method of Lie transform normalization. We study the most relevant resonances and related bifurcations. We find that the 1:1 resonance is associated only to the appearance of the loops and leads to the destabilization of either one or the other normal modes, depending on the ellipticity of equipotentials. Inclined orbits are never present and may appear only when the equipotentials are heavily deformed. The 1:2 resonance determines the appearance of bananas and anti-banana orbits: the first family is stable and always appears at a lower energy than the second, which is unstable. The bifurcation sequence also produces the variations in the stability character of the major axis orbit and is modified only by very large deformations of the equipotentials. Higher-order resonances appear at intermediate or higher energies and can be described with good accuracy.

We investigate an analytical treatment of bifurcations of families of resonant `thin' tubes in axisymmetric galactic potentials. We verify that the most relevant bifurcations are due to the (1:1) resonance producing the `inclined' orbits through two different mechanisms: from the disk orbit and from the `thin' tube associated to the vertical oscillation. The closest resonances occurring after these are the (4:3) resonance in the oblate case and the (2:1) resonance in the prolate case. The (1:1) resonances are treated in a straightforward way using a 2nd-order truncated normal form. The higher-order resonances are instead cumbersome to investigate, because the normal form has to be truncated to a high degree and the number of terms grows very rapidly. We therefore adopt a further simplification giving analytic formulas for the values of the parameters at which bifurcations ensue and compare them with selected numerical results. Thanks to the asymptotic nature of the series involved, the predictions are reliable well beyond the convergence radius of the original series.

This paper aims to illustrate the applications of resonant Hamiltonian normal forms to some problems of galactic dynamics. We detail the construction of the 1:1 resonant normal form corresponding to a wide class of potentials with self-similar elliptical equi-potentials and apply it to investigate relevant features of the orbit structure of the system. We show how to compute the bifurcation of the main periodic orbits in the symmetry planes of a triaxial ellipsoid and in the meridional plane of an axi-symmetric spheroid and briefly discuss how to refine these results with higher-order approaches.

Analytic methods to investigate periodic orbits in galactic potentials. To evaluate the quality of the approximation of periodic orbits in the logarithmic potential constructed using perturbation theory based on Hamiltonian normal forms. The solutions of the equations of motion corresponding to periodic orbits are obtained as series expansions computed by inverting the normalizing canonical transformation. To improve the convergence of the series a resummation based on a continued fraction may be performed. This method is analogous to that looking for approximate rational solutions (Prendergast method). It is shown that with a normal form truncated at the lowest order incorporating the relevant resonance it is possible to construct quite accurate solutions both for normal modes and periodic orbits in general position.

The study of self-gravitating stellar systems has provided important hints to develop tools of analytical mechanics. In the present contribution we review how to exploit detuned resonant normal forms to extract information on several aspects of the dynamics in systems with self-similar elliptical equipotentials. In particular, using energy and ellipticity as parameters, we compute the instability thresholds of axial orbits, bifurcation values of low-order boxlets and phase-space fractions pertaining to the families around them. We also show how to infer something about the singular limit of the potential.

The phase-space structure of two families of galactic potentials is approximated with a resonant detuned normal form. The normal form series is obtained by a Lie transform of the series expansion around the minimum of the original Hamiltonian. Attention is focused on the quantitative predictive ability of the normal form. We find analytical expressions for bifurcations of periodic orbits and compare them with other analytical approaches and with numerical results. The predictions are quite reliable even outside the convergence radius of the perturbation and we analyze this result using resummation techniques of asymptotic series.

We investigate the dynamics in the logarithmic galactic potential with an analytical approach. The phase-space structure of the real system is approximated with resonant detuned normal forms constructed with the method based on the Lie transform. Attention is focused on the properties of the axial periodic orbits and of low order `boxlets' that play an important role in galactic models. Using energy and ellipticity as parameters, we find analytical expressions of several useful indicators, such as stability-instability thresholds, bifurcations and phase-space fractions of some orbit families and compare them with numerical results available in the literature.

We investigate the dynamics in a galactic potential with two reflection symmetries. The phase-space structure of the real system is approximated with a resonant detuned normal form constructed with the method based on the Lie transform. Attention is focused on the stability properties of the axial periodic orbits that play an important role in galactic models. Using energy and ellipticity as parameters, we find analytical expressions of bifurcations and compare them with numerical results available in the literature.

In this paper we investigate the gravothermal instability of spherical stellar systems endowed with a radially anisotropic velocity distribution. We focus our attention on the effects of anisotropy on the conditions for the onset of the instability and in particular we study the dependence of the spatial structure of critical models on the amount of anisotropy present in a system. The investigation has been carried out by the method of linear series which has already been used in the past to study the gravothermal instability of isotropic systems. We consider models described by King, Wilson and Woolley-Dickens distribution functions. In the case of King and Woolley-Dickens models, our results show that, for quite a wide range of amount of anisotropy in the system, the critical value of the concentration of the system (defined as the ratio of the tidal to the King core radius of the system) is approximately constant and equal to the corresponding value for isotropic systems. Only for very anisotropic systems the critical value of the concentration starts to change and it decreases significantly as the anisotropy increases and penetrates the inner parts of the system. For Wilson models the decrease of the concentration of critical models is preceded by an intermediate regime in which critical concentration increases, it reaches a maximum and then it starts to decrease. The critical value of the central potential always decreases as the anisotropy increases.