Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $\Delta a$ of a celestial body on its size, spinning obliquity, initial orbit and thermal parameters on its surface. With appropriate simplification and approximation, we obtain the analytical solutions to the perturbation equations for the motion of asteroids influenced by the Yarkovsky effect, and they are then verified by numerical simulations of the full equations of motion. These solutions present explicitly the dependencies of $\Delta a$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $\Delta a$, we investigate the combined seasonal and diurnal Yarkovsky effects. The critical points where the migration direction reverses are calculated and the consequent selective effects according to the size and rotation state of asteroids are discussed. %Solely the Yarkovsky effect is found to be able to produce some ring structure in the aged circumstellar debris disk. Finally, we apply the analytical formulae to calculate the migration of Eos family members. The space distribution of asteroids is well reproduced. Our calculations suggest that statistically the orientations of spin axes of family members satisfy a random-obliquity distribution, and the rotation rate $\omega_{\rm rot}$ of asteroid depends on its size $R$ by $\omega_{\rm rot}\propto R^{-1}$.

We study the presence of resonances in the region of space around the Earth. We consider a massless body (e.g, a dust particle or a small space debris) subject to different forces: the gravitational attraction of the geopotential, the effects of Sun and Moon. We distinguish different types of resonances: tesseral resonances are due to a commensurability involving the revolution of the particle and the rotation of the Earth, semi-secular resonances include the rates of variation of the mean anomalies of Moon and Sun, while secular resonances just depend on the rates of variation of the arguments of perigee and the longitudes of the ascending nodes of the perturbing bodies. We characterize such resonances, giving precise statements on the regions where the resonances can be found and provide examples of some specific commensurability relations.

We investigate the dynamics of charged dust close to outer mean-motion resonances with planet Jupiter. The importance of the interplanetary magnetic field on the orbital evolution of dust is clearly demonstrated. New dynamical phenomena are found that do not exist in the classical problem of uncharged dust. We find changes in the orientation of the orbital planes of dust particles, an increased amount of chaotic orbital motions, sudden 'jumps' in the resonant argument, and a decrease in time of temporary capture due to the Lorentz force. Variations in the orbital planes of dust grain orbits are found to be related to the angle between the orbital angular momentum and magnetic axes of the heliospheric field and the rotation rate of the Sun. These variations are bound using a simplified model derived from the full dynamical problem using first order averaging theory. It is found that the interplanetary magnetic field does not affect the capture process, that is still dominated by the other non-gravitational forces. Our study is based on a dynamical model in the framework of the inclined circular restricted three-body problem. Additional forces include solar radiation pressure, solar wind drag, the Poynting-Robertson effect, and the influence of a Parker spiral type interplanetary magnetic field model. The analytical estimates are derived on the basis of Gauss' form of planetary equations of motion. Numerical results are obtained by simulations of dust grain orbits together with the system of variational equations. Chaotic regions in phase space are revealed by means of Fast Lyapunov Chaos Indicators.

We investigate the collective behaviour of particle orbits in the vicinity of magnetic reconnection in Earth's magneto-tail. Various regions of different kinds of orbital stability of particle motions are found. We locate regimes of temporary capture of particle orbits in configuration space as well as locations, where strong particle accelerations take place. With this study we are able to provide a detailed map, i.e. the topology, of high and low acceleration centers close to the reconnection site. Quasi-regular and chaotic kinds of motions of elementary particles can be determined as well. The orbital stability of particle orbits is obtained by a statistical analysis of the outcome of the system of variational equations of particle orbits within the framework of particle-in-cell simulations. Using the concept of Lyapunov Characteristic Numbers to ensembles of particle orbits we introduce Lyapunov Ensemble Averages to describe the response of particle orbits to local perturbations induced by the electro-magnetic field.

The aim of this study is first to determine the gravity field of the comet 67P Churyumov-Gerasimenko and second to derive the solar component of the precession rate and nutation coefficients of the spin axis of the comet nucleus, i.e. without the direct, usually larger, effect of outgassing. The gravity field, and related moments of inertia, are obtained from two polyhedra, that are provided by the OSIRIS and NAV-CAM experiments on Rosetta, and are based on the assumption of uniform density for the comet nucleus. We also calculate the forced precession rate as well as the nutation coefficients on the basis of Kinoshita's theory of rotation of the rigid Earth and adapted it to be able to indirectly include the effect of outgassing on the rotational parameters. The 2nd degree denormalized Stokes coefficients of comet 67P C-G turn out to be (bracketed numbers refer to second shape model) C20 = -6.74 [-7.93] * 10^-2, C22 = 2.60 [2.71] * 10^-2 consistent with normalized principal moments of inertia A/MR^2 = 0.13 [0.11], B/MR^2 = 0.23 [0.22], with polar moment c = C/MR^2 = 0.25, depending on the choice of the polyhedron model. The obliquity between the rotation axis and the mean orbit normal is 52 degree, and the precession rate only due to solar torques becomes 20-30 "/y. Oscillations in longitude caused by the gravitational pull of the Sun turn out to be of the order of 1', oscillations in obliquity can be estimated to be of the order of 0.5'.

We investigate the combined effect of solar wind, Poynting-Robertson drag, and the frozen-in interplanetary magnetic field on the motion of charged dust grains in our solar system. For this reason we derive a secular theory of motion by the means of averaging method and validate it with numerical simulations of the un-averaged equations of motions. The theory predicts that the secular motion of charged particles is mainly affected by the z-component of the solar magnetic axis, or the normal component of the interplanetary magnetic field. The normal component of the interplanetary magnetic field leads to an increase or decrease of semi-major axis depending on its functional form and sign of charge of the dust grain. It is generally accepted that the combined effects of solar wind and photon absorption and re-emmision (Poynting-Robertson drag) lead to a decrease in semi-major axis on secular time scales. On the contrary, we demonstrate that the interplanetary magnetic field may counteract these drag forces under certain circumstances. We derive a simple relation between the parameters of the magnetic field, the physical properties of the dust grain as well as the shape and orientation of the orbital ellipse of the particle, which is a necessary conditions for the stabilization in semi-major axis.

We analyze the combined effect of Poynting-Robertson and solar wind drag on space debris. We derive a model within Cartesian, Gaussian and Hamiltonian frameworks. We focus on the geosynchronous resonance, although the results can be easily generalized to any resonance. By numerical and analytical techniques, we compute the drift in semi-major axis due to Poynting-Robertson and solar wind drag. After a linear stability analysis of the equilibria, we combine a careful investigation of the regular, resonant, chaotic behavior of the phase space with a long-term propagation of a sample of initial conditions. The results strongly depend on the value of the area-to-mass ratio of the debris, which might show different dynamical behaviors: temporary capture or escape from the geosynchronous resonance, as well as temporary capture or escape from secondary resonances involving the rate of variation of the longitude of the Sun. Such analysis shows that Poynting-Robertson and solar wind drag must be taken into account, when looking at the long-term behavior of space debris. Trapping or escape from the resonance can be used to place the debris in convenient regions of the phase space.

Mercury is the unique known planet that is situated in a 3:2 spin-orbit resonance nowadays. Observations and models converge to the same conclusion: the planet is presently deeply trapped in the resonance and situated at the Cassini state $1$, or very close to it. We investigate the complete non-linear stability of this equilibrium, with respect to several physical parameters, in the framework of Birkhoff normal form and Nekhoroshev stability theory. We use the same approach adopted for the 1:1 spin-orbit case with a peculiar attention to the role of Mercury's non negligible eccentricity. The selected parameters are the polar moment of inertia, the Mercury's inclination and eccentricity and the precession rates of the perihelion and node. Our study produces a bound to both the latitudinal and longitudinal librations (of 0.1 radians) for a long but finite time (greatly exceeding the age of the solar system). This is the so-called effective stability time. Our conclusion is that Mercury, placed inside the 3:2 spin-orbit resonance, occupies a very stable position in the space of these physical parameters, but not the most stable possible one.

We investigate the stability of motion close to the Lagrangian equilibrium points L4 and L5 in the framework of the spatial, elliptic, restricted three- body problem, subject to the radial component of Poynting-Robertson drag. For this reason we develop a simplified resonant model, that is based on averaging theory, i.e. averaged over the mean anomaly of the perturbing planet. We find temporary stability of particles displaying a tadpole motion in the 1:1 resonance. From the linear stability study of the averaged simplified resonant model, we find that the time of temporary stability is proportional to beta a1 n1 , where beta is the ratio of the solar radiation over the gravitational force, and a1, n1 are the semi-major axis and the mean motion of the perturbing planet, respectively. We extend previous results (Murray (1994)) on the asymmetry of the stability indices of L4 and L5 to a more realistic force model. Our analytical results are supported by means of numerical simulations. We implement our study to Jupiter-like perturbing planets, that are also found in extra-solar planetary systems.

In a dissipative system the time to reach an attractor is often influenced by the peculiarities of the model and in particular by the strength of the dissipation. In particular, as a dissipative model we consider the spin-orbit problem providing the dynamics of a triaxial satellite orbiting around a central planet and affected by tidal torques. The model is ruled by the oblateness parameter of the satellite, the orbital eccentricity, the dissipative parameter and the drift term. We devise a method which provides a reliable indication on the transient time which is needed to reach an attractor in the spin-orbit model; the method is based on an analytical result, precisely a suitable normal form construction. This method provides also information about the frequency of motion. A variant of such normal form used to parametrize invariant attractors provides a specific formula for the drift parameter, which in turn yields a constraint - which might be of interest in astronomical problems - between the oblateness of the satellite and its orbital eccentricity.

Earth-based radar observations of the rotational dynamics of Mercury (Margot et al. 2012) combined with the determination of its gravity field by MESSENGER (Smith et al. 2012) give clues on the internal structure of Mercury, in particular its polar moment of inertia C, deduced from the obliquity (2.04 +/- 0.08) arcmin. The dynamics of the obliquity of Mercury is a very-long term motion (a few hundreds of kyrs), based on the regressional motion of Mercury's orbital ascending node. This paper, following the study of Noyelles & D'Hoedt (2012), aims at first giving initial conditions at any time and for any values of the internal structure parameters for numerical simulations, and at using them to estimate the influence of usually neglected parameters on the obliquity, like J3, the Love number k2 and the secular variations of the orbital elements. We use for that averaged representations of the orbital and rotational motions of Mercury, suitable for long-term studies. We find that J3 should alter the obliquity by 250 milli-arcsec, the tides by 100 milli-arcsec, and the secular variations of the orbital elements by 10 milli-arcsec over 20 years. The resulting value of C could be at the most changed from 0.346mR^2 to 0.345mR^2.

Recently the first Earth Trojan has been observed (Mainzer et al., ApJ 731) and found to be on an interesting orbit close to the Lagrange point L4 (Connors et al., Nature 475). In the present study we therefore perform a detailed investigation on the stability of its orbit and moreover extend the study to give an idea of the probability to find additional Earth-Trojans. Our results are derived using different approaches: a) we derive an analytical mapping in the spatial elliptic restricted three-body problem to find the phase space structure of the dynamical problem. We explore the stability of the asteroid in the context of the phase space geometry, including the indirect influence of the additional planets of our Solar system. b) We use precise numerical methods to integrate the orbit forward and backward in time in different dynamical models. Based on a set of 400 clone orbits we derive the probability of capture and escape of the Earth Trojan asteroids 2010 TK7. c) To this end we perform an extensive numerical investigation of the stability region of the Earth's Lagrangian points. We present a detailed parameter study in the regime of possible stable tadpole and horseshoe orbits of additional Earth-Trojans, i.e. with respect to the semi-major axes and inclinations of thousands of fictitious Trojans. All three approaches underline that the Earth Trojan asteroid 2010 TK7 finds himself in an unstable region on the edge of a stable zone; additional Earth-Trojan asteroids may be found in this regime of stability.

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore we did numerical experiments in different dynamical models with different initial conditions -- on one hand the couple Venus-Earth was set close to different mean motion resonances (MMR), and on the other hand Venus' orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40deg). The couple Venus-Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth's semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in eccentricity and/or inclination of Venus in these resonances. Even escapes for the innermost planet are possible which may happen quite rapidly.