Jupiter Trojan asteroids are located around L4 and L5 Lagrangian points on relatively stable orbits, in 1:1 MMR with Jupiter. However, not all of them lie in orbits that remain stable over the age of the Solar System. Unstable zones allow some Trojans to escape in time scales shorter than the Solar System age. This may contribute to populate other small body populations. In this paper, we study this process by performing long-term numerical simulations of the observed Trojans, focusing on the trajectories of those that leave the resonance. The orbits of current Trojans are taken as initial conditions and their evolution is followed under the gravitational action of the Sun and the planets. We find the rate of escape of Trojans from L5, ~1.1 times greater than from L4. The majority of escaped Trojans have encounters with Jupiter although they have encounters with the other planets too. Almost all escaped Trojans reach the comet zone, ~90% cross the Centaur zone and only L4 Trojans reach the transneptunian zone. Considering the real asymmetry between L4 and L5, we show that 18 L4 Trojans and 14 L5 Trojans with diameter D > 1 km are ejected from the resonance every Myr. The contribution of the escaped Trojans to other minor body populations would be negligible, being the contribution from L4 and L5 to JFCs and no-JFCs almost the same, and the L4 contribution to Centaurs and TNOs, orders of magnitude greater than that of L5. Considering the collisional removal, besides the dynamical one, and assuming that Trojans that escape due to collisions follow the same dynamical behavior that the ones removed by dynamics, we would have a minor contribution of Trojans to comets and Centaurs. However, there would be some specific regions were escaped Trojans could be important such as ACOs, Encke-type comets, S-L 9-type impacts on Jupiter and NEOs.

Based on the model described in Ramos et al., 2017, we present an analytical+numerical study of the resonance capture under Type-I migration for the Kepler-25 (Marcy et al., 2014) and K2-24 (Petigura et al., 2016) Kepler systems, both close to a 2/1 mean-motion resonance. We find that, depending on the flare index and the proximity to the central star, the average value of the period-ratio between two consecutive planets show a significant deviation with respect to the resonant nominal value, up to values well in agreement with the observations.

We present a series of dynamical maps for fictitious 3-planets systems in initially circular coplanar orbits. These maps have unveiled a rich resonant structure involving two or three planets, as well as indicating possible migration routes from secular to double resonances or pure 3-planet commensurabilities. These structures are then compared to the present-day orbital architecture of observed resonant chains. In a second part of the paper we describe N-body simulations of type-I migration. Depending on the orbital decay timescale, we show that 3-planet systems may be trapped in different combinations of independent commensurabilities: (i) double resonances, (ii) intersection between a 2-planet and a first-order 3-planet resonance, and (iii) simultaneous libration in two first-order 3-planet resonances. These latter outcomes are found for slow migrations, while double resonances are almost always the final outcome in high-density disks. Finally, we discuss an application to the TRAPPIST-1 system. We find that, for low migration rates and planetary masses of the order of the estimated values, most 3-planet sub-systems are able to reach the observed double resonances after following evolutionary routes defined by pure 3-planet resonances. The final orbital configuration shows resonance offsets comparable with present-day values without the need of tidal dissipation. For the 8/5 resonance proposed to dominate the dynamics of the two inner planets, we find little evidence of its dynamical significance; instead, we propose that this relation between mean motions could be a consequence of the interaction between a pure 3-planet resonance and a 2-planet commensurability between planets c and d.

We present an analytical and numerical study of the orbital migration and resonance capture of fictitious two-planet systems with masses in the super-Earth range undergoing Type-I migration. We find that, depending on the flare index and proximity to the central star, the average value of the period ratio, $P_2/P_1$, between both planets may show a significant deviation with respect to the nominal value. For planets trapped in the 2:1 commensurability, offsets may reach values on the order of $0.1$ for orbital periods on the order of $1$ day, while systems in the 3:2 mean-motion resonance (MMR) show much smaller offsets for all values of the semimajor axis. These properties are in good agreement with the observed distribution of near-resonant exoplanets, independent of their detection method. We show that 2:1-resonant systems far from the star, such as HD82943 and HR8799, are characterized by very small resonant offsets, while higher values are typical of systems discovered by Kepler with orbital periods approximately a few days. Conversely, planetary systems in the vicinity of the 3:2 MMR show little offset with no significant dependence on the orbital distance. In conclusion, our results indicate that the distribution of Kepler planetary systems around the 2:1 and 3:2 MMR are consistent with resonant configurations obtained as a consequence of a smooth migration in a laminar flared disk, and no external forces are required to induce the observed offset or its dependence with the commensurability or orbital distance from the star.

We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.

We present a series of numerical integrations of observed and fictitious Jupiter Trojan asteroids, under the gravitational effects of the four outer planets, for time-spans comparable with the age of the Solar System. From these results we calculate the escape rate from each Lagrange point, and construct dynamical maps of "permanence" time in different regions of the phase space. Fictitious asteroids in L4 and L5 show no significant difference, showing almost identical dynamical maps and escape rates. For real Trojans, however, we found that approximately 23% o f the members of the leading swarm escaped after 4.5 Gyrs, while this number increased to 28.3% for L5. This implies that the asymmetry between the two populations increases with time, indicating that it may have been smaller at the time of formation/capture of these asteroids. Nevertheless, the difference in chaotic diffusion cannot, in itself, account for the current observed asymmetry (~40%), and must be primarily primordial and characteristic of the capture mechanism of the Trojans. Finally, we calculate new proper elements for all the numbered Trojans using the semi-analytical approach of Beaug\'e and Roig (2001), and compare the results with the numerical estimations by Bro\v{z} and Rosehnal (2011). For asteroids that were already numbered in 2011, both methods yield very similar results, while significant differences were found for those bodies that became numbered after 2011.

We review the orbital stability of the planar circular restricted three-body problem, in the case of massless particles initially located between both massive bodies. We present new estimates of the resonance overlap criterion and the Hill stability limit, and compare their predictions with detailed dynamical maps constructed with N-body simulations. We show that the boundary between (Hill) stable and unstable orbits is not smooth but characterized by a rich structure generated by the superposition of different mean-motion resonances which does not allow for a simple global expression for stability. We propose that, for a given perturbing mass $m_1$ and initial eccentricity $e$, there are actually two critical values of the semimajor axis. All values $a < a_{\rm Hill}$ are Hill-stable, while all values $a > a_{\rm unstable}$ are unstable in the Hill sense. The first limit is given by the Hill-stability criterion and is a function of the eccentricity. The second limit is virtually insensitive to the initial eccentricity, and closely resembles a new resonance overlap condition (for circular orbits) developed in terms of the intersection between first and second-order mean-motion resonances.