We revisit the problem of secular resonance sweeping during the dissipation of a protoplanetary disk and its possible role in exciting the orbits of primordial asteroids, in the light of recent models of solar system evolution. We develop an integrator that incorporates the gravitational effect of a uniformly (or not) depleting, axisymmetric disk with arbitrary surface density profile; its performance is verified by analytical calculations. The secular response of fictitious asteroids, under perturbations from Jupiter, Saturn and a decaying disk, is thoroughly studied. Note that the existence of a symmetry plane induced by the disk, lifts the inherent degeneracy of the two-planet system, such that the '$s_5$' nodal frequency can also play a major role. We examine different resonant configurations for the planets (2:1, 3:2, 5:3), disk models and depletion scenarios. For every case we compute the corresponding time paths of secular resonances, which show when and where resonance crossing occurs. The excitation of asteroids, particularly in inclination, is studied in the various models and compared to analytical estimates. We find that inclination excitation in excess of $\sim 10^{\circ}$ is possible in the asteroid belt, but the nearly uniform spread of $\Delta i\sim 20^{\circ}$ observed calls for the combined action of secular resonance sweeping with other mechanisms (e.g. scattering by Mars-sized embryos) that would be operating during terrestrial planet formation. Our results are also applicable to extrasolar planetary systems.

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 have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term (120 years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes (400 km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded "ReDSHIFT" project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an "atlas" of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth's oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device.

The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbits). Recent studies have highlighted the complicated secular dynamics in the MEO region, but also the possibility of exploiting these dynamics, for designing removal strategies. In this paper, we present our numerical exploration of the long-term dynamics in MEO, performed with the purpose of unveiling the set of reentry and graveyard solutions that could be reached with maneuvers of reasonable DV cost. We simulated the dynamics over 120-200 years for an extended grid of millions of fictitious MEO satellites that covered all inclinations from 0 to 90deg, using non-averaged equations of motion and a suitable dynamical model that accounted for the principal geopotential terms, 3rd-body perturbations and solar radiation pressure (SRP). We found a sizeable set of usable solutions with reentry times that exceed ~40years, mainly around three specific inclination values: 46deg, 56deg, and 68deg; a result compatible with our understanding of MEO secular dynamics. For DV <= 300 m/s (i.e., achieved if you start from a typical GNSS orbit and target a disposal orbit with e<0.3), reentry times from GNSS altitudes exceed ~70 years, while low-cost (DV ~= 5-35 m/s) graveyard orbits, stable for at lest 200 years, are found for eccentricities up to e~0.018. This investigation was carried out in the framework of the EC-funded "ReDSHIFT" project.

In this paper, we study the long-term dynamical evolution of highly-elliptical orbits (HEOs) in the medium-Earth orbit (MEO) region around the Earth. The real population consists primarily of Geosynchronous Transfer Orbits (GTOs), launched at specific inclinations, Molniya-type satellites and related debris. We performed a suite of long-term numerical integrations (up to 200 years) within a realistic dynamical model, aimed primarily at recording the dynamical lifetime of such orbits (defined as the time needed for atmospheric reentry) and understanding its dependence on initial conditions and other parameters, such as the area-to-mass ratio (A/m). Our results are presented in the form of 2-D lifetime maps, for different values of inclination, A/m, and drag coefficient. We find that the majority of small debris (> 70%, depending on the inclination) can naturally reenter within 25-90 years, but these numbers are significantly less optimistic for large debris (e.g., upper stages), with the notable exception of those launched from high latitude (Baikonur). We estimate the reentry probability and mean dynamical lifetime for different classes of GTOs and we find that both quantities depend primarily and strongly on initial perigee altitude. Atmospheric drag and higher A/m values extend the reentry zones, especially at low inclinations. For high inclinations, this dependence is weakened, as the primary mechanisms leading to reentry are overlapping lunisolar resonances. This study forms part of the EC-funded (H2020) "ReDSHIFT" project.

We present the technical specifications and first results of the ESA-funded, lunar monitoring project "NELIOTA" (NEO Lunar Impacts and Optical TrAnsients) at the National Observatory of Athens, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. For the purposes of this project a twin camera instrument was specially designed and installed at the 1.2 m Kryoneri telescope utilizing the fast-frame capabilities of scientific Complementary Metal-Oxide Semiconductor detectors (sCMOS). The system provides a wide field-of-view (17.0' $\times$ 14.4') and simultaneous observations in two photometric bands (R and I), reaching limiting magnitudes of 18.7 mag in 10 sec in both bands at a 2.5 signal-to-noise level. This makes it a unique instrument that can be used for the detection of NEO impacts on the Moon, as well as for any astronomy projects that demand high-cadence multicolor observations. The wide field-of-view ensures that a large portion of the Moon is observed, while the simultaneous, high-cadence, monitoring in two photometric bands makes possible, for the first time, the determination of the temperatures of the impacts on the Moon's surface and the validation of the impact flashes from a single site. Considering the varying background level on the Moon's surface we demonstrate that the NELIOTA system can detect NEO impact flashes at a 2.5 signal-to-noise level of ~12.4 mag in the I-band and R-band for observations made at low lunar phases ~0.1. We report 31 NEO impact flashes detected during the first year of the NELIOTA campaign. The faintest flash was at 11.24 mag in the R-band (about two magnitudes fainter than ever observed before) at lunar phase 0.32. Our observations suggest a detection rate of $1.96 \times 10^{-7}$ events $km^{-2} h^{-1}$.

Librational motion in celestial mechanics is generally associated with the existence of stable resonant configurations and signified by the existence of stable periodic solutions and oscillation of critical (resonant) angles. When such an oscillation takes place around a value different than 0 or $\pi$, the libration is called asymmetric. In the context of the planar circular restricted three-body problem (CRTBP), asymmetric librations have been identified for the exterior mean-motion resonances (MMRs) 1:2, 1:3 etc. as well as for co-orbital motion (1:1). In exterior MMRs the massless body is the outer one. In this paper, we study asymmetric librations in the 3-dimensional space. We employ the computational approach of Markellos (1978) and compute families of asymmetric periodic orbits and their stability. Stable, asymmetric periodic orbits are surrounded in phase space by domains of initial conditions which correspond to stable evolution and librating resonant angles. Our computations were focused on the spatial circular restricted three-body model of the Sun-Neptune-TNO system (TNO= trans-Neptunian object). We compare our results with numerical integrations of observed TNOs, which reveal that some of them perform 1:2-resonant, inclined asymmetric librations. For the stable 1:2 TNOs librators, we find that their libration seems to be related with the vertically stable planar asymmetric orbits of our model, rather than the 3-dimensional ones found in the present study.

The rectilinear elliptic restricted Three Body Problem (TBP) is the limiting case of the elliptic restricted TBP when the motion of the primaries is described by a Keplerian ellipse with eccentricity $e'=1$, but the collision of the primaries is assumed to be a non-singular point. The rectilinear model has been proposed as a starting model for studying the dynamics of motion around highly eccentric binary systems. Broucke (1969) explored the rectilinear problem and obtained isolated periodic orbits for mass parameter $\mu=0.5$ (equal masses of the primaries). We found that all orbits obtained by Broucke are linearly unstable. We extend Broucke's computations by using a finer search for symmetric periodic orbits and computing their linear stability. We found a large number of periodic orbits, but only eight of them were found to be linearly stable and are associated with particular mean motion resonances. These stable orbits are used as generating orbits for continuation with respect to $\mu$ and $e'<1$. Also, continuation of periodic solutions with respect to the mass of the small body can be applied by using the general TBP. FLI maps of dynamical stability show that stable periodic orbits are surrounded in phase space with regions of regular orbits indicating that systems of very highly eccentric orbits can be found in stable resonant configurations. As an application we present a stability study for the planetary system HD7449.

In this work we aim to constrain the slope of the size distribution of main-belt asteroids, at their primordial state. To do so we turn out attention to the part of the main asteroid belt between 2.82 and 2.96~AU, the so-called "pristine zone", which has a low number density of asteroids and few, well separated asteroid families. Exploiting these unique characteristics, and using a modified version of the hierarchical clustering method we are able to remove the majority of asteroid family members from the region. The remaining, background asteroids should be of primordial origin, as the strong 5/2 and 7/3 mean-motion resonances with Jupiter inhibit transfer of asteroids to and from the neighboring regions. The size-frequency distribution of asteroids in the size range $17<D(\rm{km})<70$ has a slope $q\simeq-1$. Using Monte-Carlo methods, we are able to simulate, and compensate for the collisional and dynamical evolution of the asteroid population, and get an upper bound for its size distribution slope $q=-1.43$. In addition, applying the same 'family extraction' method to the neighboring regions, i.e. the middle and outer belts, and comparing the size distributions of the respective background populations, we find statistical evidence that no large asteroid families of primordial origin had formed in the middle or pristine zones.

It is generally accepted today that our solar system has undergone a phase during which the orbits of the giant planets became very unstable. In recent years, many studies have identified traces of this event and have provided reasonable justification for this occurrence. The magnitude (in terms of orbital variation) and the timing of the instability though (early or late with respect to the dispersal of the gas disk) still remains an open debate. The terrestrial planets seem to set a strict constraint: either the giant planet instability happened early, while the terrestrial planets were still forming, or the orbits of Jupiter and Saturn had to separate from each other impulsively, with a large enough `jump' in semimajor axis (Brasser et al. 2009; Kaib and Chambers 2016) for the terrestrial planets to remain stable. Because a large orbital jump is a low probability event, the early instability hypothesis seems to be favored. However, the asteroid belt would also evolve in a different way, assuming different instability amplitudes. These two constraints need to match each other in order to favor one scenario over the other. Considering an initially dynamically cold disk of asteroids, Morbidelli et al. (2010) concluded that a comparably large jump is needed to reconstruct the current asteroid belt. Here we confirm the same conclusion, but considering an asteroid population already strongly excited in eccentricity, such as that produced in the Grand Tack scenario (Walsh et al. 2011). Because the asteroids existed since the time of removal of the gas disk, unlike the terrestrial planets, this constraint on the width of the giant planet jump is valid regardless of whether the instability happened early or late. Hence, at this stage, assuming an early instability does not appear to provide any advantage in terms of the probabilistic reconstruction of the solar system structure.

The navigation satellite constellations in medium-Earth orbit exist in a background of third-body secular resonances stemming from the perturbing gravitational effects of the Moon and the Sun. The resulting chaotic motions, emanating from the overlapping of neighboring resonant harmonics, induce especially strong perturbations on the orbital eccentricity, which can be transported to large values, thereby increasing the collision risk to the constellations and possibly leading to a proliferation of space debris. We show here that this transport is of a diffusive nature and we present representative diffusion maps that are useful in obtaining a global comprehension of the dynamical structure of the navigation satellite orbits.

NELIOTA is a new ESA activity launched at the National Observatory of Athens in February 2015 aiming to determine the distribution and frequency of small near-earth objects (NEOs) via lunar monitoring. The project involves upgrading the 1.2m Kryoneri telescope at the National Observatory of Athens, procuring two fast-frame cameras, and developing a software system, which will control the telescope and the cameras, process the images and automatically detect NEO impacts. NELIOTA will provide a web-based user interface, where the impact events will be reported and made available to the scientific community and the general public. The objective of this 3.5 year activity is to design, develop and implement a highly automated lunar monitoring system, which will conduct an observing campaign for 2 years in search of NEO impact flashes on the Moon. The impact events will be verified, characterised and reported. The 1.2m telescope will be capable of detecting flashes much fainter than current, small-aperture, lunar monitoring telescopes. NELIOTA is therefore expected to characterise the frequency and distribution of NEOs weighing as little as a few grams.

The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of asteroids. The short-lived population is known to be a transient group sustained in steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, can exhibit dynamical lifetimes comparable to $4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space denoted $\mathrm{A}$ and $\mathrm{B}$, with an uneven population ratio $\mathrm{B}/\mathrm{A} \simeq 10$. The orbits of $\mathrm{A}$-island asteroids are predominantly highly inclined, compared to island $\mathrm{B}$. The size-frequency distribution is steep but the orbital distribution lacks any evidence of a collisional cluster. These observational constraints are somewhat puzzling and therefore the origin of the long-lived asteroids has not been explained so far. With the aim to provide a viable explanation, we first update the resonant population and revisit its physical properties. Using an $N$-body model with seven planets and the Yarkovsky effect included, we demonstrate that the dynamical depletion of island $\mathrm{A}$ is faster, in comparison with island $\mathrm{B}$. Then we investigate (i) the survivability of primordial resonant asteroids and (ii) capture of the population during planetary migration, following a recently described scenario with an escaping fifth giant planet and a jumping-Jupiter instability. We also model the collisional evolution of the resonant population over past $4\,\mathrm{Gyr}$. Our conclusion is that the long-lived group was created by resonant capture from a narrow part of hypothetical outer main-belt family during planetary migration. Primordial asteroids surviving the migration were probably not numerous enough to substantially contribute to the observed population.

We consider a two-planet system, which migrates under the influence of dissipative forces that mimic the effects of gas-driven (Type II) migration. It has been shown that, in the planar case, migration leads to resonant capture after an evolution that forces the system to follow families of periodic orbits. Starting with planets that differ slightly from a coplanar configuration, capture can, also, occur and, additionally, excitation of planetary inclinations has been observed in some cases. We show that excitation of inclinations occurs, when the planar families of periodic orbits, which are followed during the initial stages of planetary migration, become vertically unstable. At these points, {\em vertical critical orbits} may give rise to generating stable families of $3D$ periodic orbits, which drive the evolution of the migrating planets to non-coplanar motion. We have computed and present here the vertical critical orbits of the $2/1$ and $3/1$ resonances, for various values of the planetary mass ratio. Moreover, we determine the limiting values of eccentricity for which the "inclination resonance" occurs.

We show that the existence of prograde equatorial satellites is consistent with a collisional tilting scenario for Uranus. In fact, if the planet was surrounded by a proto-satellite disk at the time of the tilting and a massive ring of material was temporarily placed inside the Roche radius of the planet by the collision, the proto-satellite disk would have started to precess incoherently around the equator of the planet, up to a distance greater than that of Oberon. Collisional damping would then have collapsed it into a thin equatorial disk, from which the satellites eventually formed. The fact that the orbits of the satellites are prograde requires Uranus to have had a non-negligible initial obliquity (comparable to that of Neptune) before it was finally tilted to 98 degrees.

The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely, that some exoplanets are not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing planetary system, created in a tight cluster, and is captured as a new planet. In the present work we study this interaction process by integrating trajectories of planet-sized bodies, which encounter a binary system consisting of a Jupiter-sized planet revolving around a Sun-like star. To simplify the problem we assume coplanar orbits for the bound and the free-floating planet and an initially parabolic orbit for the free-floating planet. By calculating the uncertainty exponent, a quantity that measures the dependence of the final state of the system on small changes of the initial conditions, we show that the interaction process is a fractal classical scattering. The uncertainty exponent is in the range (0.2-0.3) and is a decreasing function of time. In this way we see that the statistical approach we follow to tackle the problem is justified. The possible final outcomes of this interaction are only four, namely flyby, planet exchange, capture or disruption. We give the probability of each outcome as a function of the incoming planet's mass. We find that the probability of exchange or capture (in prograde as well as retrograde orbits and for very long times) is non-negligible, a fact that might explain the possible future observations of planetary systems with orbits that are either retrograde or tight and highly eccentric.

We study the establishment of three-planet resonances -similar to the Laplace resonance in the Galilean satellites- and their effects on the mutual inclinations of the orbital planes of the planets, assuming that the latter undergo migration in a gaseous disc. In particular, we examine the resonance relations that occur, by varying the physical and initial orbital parameters of the planets (mass, initial semi-major axis and eccentricity) as well as the parameters of the migration forces (migration rate and eccentricity damping rate), which are modeled here through a simplified analytic prescription. We find that, in general, for planetary masses below 1.5 M_J, multiple-planet resonances of the form n3:n2:n1=1:2:4 and 1:3:6 are established, as the inner planets, m1 and m2, get trapped in a 1:2 resonance and the outer planet m3 subsequently is captured in a 1:2 or 1:3 resonance with m2. For mild eccentricity damping, the resonance pumps the eccentricities of all planets on a relatively short time-scale, to the point where they enter an inclination-type resonance (as in Libert & Tsiganis 2011); then mutual inclinations can grow to ~35{\deg}, thus forming a "3-D system". On the other hand, we find that trapping of m2 in a 2:3 resonance with m1 occurs very rarely, for the range of masses used here, so only two cases of capture in a respective three-planet resonance were found. Our results suggest that trapping in a three-planet resonance can be common in exoplanetary systems, provided that the planets are not very massive. Inclination pumping could then occur relatively fast, provided that eccentricity damping is not very efficient so that at least one of the inner planets acquires an orbital eccentricity higher than e=0.3.

Assuming that giant planets are formed in thin protoplanetary discs, a '3D' system can form, provided that the mutual inclination is excited by some dynamical mechanism. Resonant interactions and close planetary encounters are thought to be the primary inclination-excitation mechanisms, resulting in a resonant and non-resonant system, respectively. Here we propose an alternative formation scenario, starting from a system composed of three giant planets in a nearly coplanar configuration. As was recently shown for the case of the Solar system, planetary migration in the gas disc (Type II migration) can force the planets to become trapped in a multiply resonant state. We simulate this process, assuming different values for the planetary masses and mass ratios. We show that such a triple resonance generally becomes unstable as the resonance excites the eccentricities of all planets and planet-planet scattering sets in. One of the three planets is typically ejected from the system, leaving behind a dynamically 'hot' (but stable) two-planet configuration. The resulting two-planet systems typically have large values of semimajor axial ratios (a1/a2 < 0.3), while the mutual inclination can be as high as 70{\deg}, with a median of \sim30{\deg}. A small fraction of our two-planet systems (\sim5 per cent) ends up in the stability zone of the Kozai resonance. In a few cases, the triple resonance can remain stable for long times and a '3D' system can form by resonant excitation of the orbital inclinations; such a three-planet system could be stable if enough eccentricity damping is exerted on the planets. Finally, in the single-planet resulting systems, which are formed when two planets are ejected from the system, the inclination of the planet's orbital plane with respect to the initial invariant plane -presumably the plane perpendicular to the star's spin axis- can be as large as \sim40{\deg}.

The presence of highly eccentric extrasolar planets in binary stellar systems suggests that the Kozai effect has played an important role in shaping their dynamical architectures. However, the formation of planets in inclined binary systems poses a considerable theoretical challenge, as orbital excitation due to the Kozai resonance implies destructive, high-velocity collisions among planetesimals. To resolve the apparent difficulties posed by Kozai resonance, we seek to identify the primary physical processes responsible for inhibiting the action of Kozai cycles in protoplanetary disks. Subsequently, we seek to understand how newly-formed planetary systems transition to their observed, Kozai-dominated dynamical states. We find that theoretical difficulties in planet formation arising from the presence of a distant companion star, posed by the Kozai effect and other secular perturbations, can be overcome by a proper account of gravitational interactions within the protoplanetary disk. In particular, fast apsidal recession induced by disk self-gravity tends to erase the Kozai effect, and ensure that the disk's unwarped, rigid structure is maintained. Subsequently, once a planetary system has formed, the Kozai effect can continue to be wiped out as a result of apsidal precession, arising from planet-planet interactions. However, if such a system undergoes a dynamical instability, its architecture may change in such a way that the Kozai effect becomes operative. The results presented here suggest that planetary formation in highly inclined binary systems is not stalled by perturbations, arising from the stellar companion. Consequently, planet formation in binary stars is probably no different from that around single stars on a qualitative level. Furthermore, it is likely that systems where the Kozai effect operates, underwent a transient phase of dynamical instability in the past.

We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale of tau ~ 0.5My. However, self-consistent numerical simulations show that the planetesimal-driven migration of the giant planets is inconsistent with an exponential change in their semi major axes on such a short time scale (Hahn & Malhotra, 1999). In fact, the typical time scale is tau > 5My. When giant planet migration on this time scale is applied to the asteroid belt, the resulting orbital distribution is incompatible with the observed one. However, the planet migration can be significantly sped up by planet-planet encounters. Consider an evolution where both Jupiter and Saturn have close encounters with a Neptune-mass planet (presumably Uranus or Neptune themselves) and where this third planet, after being scattered inwards by Saturn, is scattered outwards by Jupiter. This scenario leads to a very rapid increase in the orbital separation between Jupiter and Saturn that we show here to have only mild effects on the structure of asteroid belt. This type of evolution is called a jumping-Jupiter case. Our results suggest that the total mass and dynamical excitation of the asteroid belt before migration were comparable to those currently observed. Moreover, they imply that, before migration, the orbits of Jupiter and Saturn were much less eccentric than the current ones.

We present a transport model that describes the orbital diffusion of asteroids in chaotic regions of the 3-D space of proper elements. Our goal is to use a simple random-walk model to study the evolution and derive accurate age estimates for dynamically complex asteroid families. To this purpose, we first compute local diffusion coefficients, which characterize chaotic diffusion in proper eccentricity (e_p) and inclination (I_p), in a selected phase-space region. Then, a Monte-Carlo-type code is constructed and used to track the evolution of random walkers (i.e. asteroids), by coupling diffusion in (e_p,I_p) with a drift in proper semi-major axis (a_p) induced by the Yarkovsky/YORP thermal effects. We validate our model by applying it to the family of (490) Veritas, for which we recover previous estimates of its age (~8.3 Myr). Moreover, we show that the spreading of chaotic family members in proper elements space is well reproduced in our random-walk simulations. Finally, we apply our model to the family of (3556) Lixiaohua, which is much older than Veritas and thus much more affected by thermal forces. We find the age of the Lixiaohua family to be 155+/-36 Myr.

We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes the g5 eigenfrequency to cross resonances of the form g5=gk with k ranging from 1 to 4. Consequently these secular resonances cause large-amplitude responses in the eccentricities of the terrestrial planets. We show that the resonances g5=g_4 and g5=g3 do not pose a problem if Jupiter and Saturn have a fast approach and departure from their mutual 2:1 mean motion resonance. On the other hand, the resonance crossings g5=g2 and g5=g1 are more of a concern as they tend to yield a terrestrial system incompatible with the current one. We offer two solutions to this problem. The first uses the fact that a secular resonance crossing can also damp the amplitude of a Fourier mode if the latter is large originally. A second scenario involves a 'jumping Jupiter' in which encounters between an ice giant and Jupiter, without ejection of the former, cause the latter to migrate away from Saturn much faster than if migration is driven solely by encounters with planetesimals. In this case, the g5=g2 and g5=g1 resonances can be jumped over, or occur very briefly.

Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations. Moreover, the amplitudes of the eigenmodes characterising the current secular evolution of the eccentricities of Jupiter and Saturn would not be reproduced correctly; only one eigenmode is excited by resonance-crossing. We show that, at the very least, encounters between Saturn and one of the ice giants (Uranus or Neptune) need to have occurred, in order to reproduce the current secular properties of the giant planets, in particular the amplitude of the two strongest eigenmodes in the eccentricities of Jupiter and Saturn.

We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ~35 AU to the Kuiper belt without excessive inclination excitation. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ~27 AU and an eccentricity of ~0.3, and that subsequently Neptune's eccentricity damped in ~1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level.

We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonance (i.e. each planet is in resonance with its neighbor). In total we found 6 such configurations. For the gas disk parameters found in MC07, these configurations are characterized by a negligible migration rate. After the disappearance of the gas, and in absence of planetesimals, only two of these six configurations (the least compact ones) are stable for a time of hundreds of millions of years or more. The others become unstable on a timescale of a few My. Our preliminary simulations show that, when a planetesimal disk is added beyond the orbit of the outermost planet, the planets can evolve from the most stable of these configurations to their current orbits in a fashion qualitatively similar to that described in Tsiganis et al. (2005).

Among the many new irregular satellites that have been discovered in the last five years, at least six are in the so-called Kozai resonance. Due to solar perturbations, the argument of pericenter of a satellite usually precesses from 0 to 360 degrees. However, at inclinations higher than 39.3 degrees and lower than 140.7 degrees a new kind of behavior occurs for which the argument of pericenter oscillates around +/-90 degrees. In this work we will concentrate on the orbital history of the saturnian satellite S/2000 S5 Kiviuq, one of the satellites currently known to be in such resonance Kiviuq's orbit is very close to the separatrix of the Kozai resonance. Due to perturbations from the other jovian planets, it is expected that orbits near the Kozai separatrix may show significant chaotic behavior. This is important because chaotic diffusion may transfer orbits from libration to circulation, and vice versa. To identify chaotic orbits we used two well-known methods: the Frequency Analysis Method (Laskar 1990) and Maximum Lyapunov Exponents (Benettin et al. 1980). Our results show that the Kozai resonance is crossed by a web of secondary resonances, whose arguments involve combinations of the argument of pericenter, the argument of the Great Inequality, longitude of the node, and other terms related to the secular frequencies g5, g6, and s6. Many test orbits whose precession period is close to the period of the Great Inequality (883 yrs), or some of its harmonics, are trapped by these secondary resonances, and show significant chaotic behavior. Planetary migration, by moving the locations of these secondary resonances, may have depleted an original population of Kozai resonators.