The current giant planet region is a transitional zone where transneptunian objects (TNOs) cross in their way to becoming Jupiter Family Comets (JFCs). Their dynamical behavior is conditioned by the intrinsic dynamical features of TNOs and also by the encounters with the giant planets. We address the Giant Planet Crossing (GPC) population (those objects with $5.2$ au $ < q < 30$ au) studying their number and their evolution from their sources, considering the current configuration of the Solar System. This subject is reviewed from previous investigations and also addressed by new numerical simulations of the dynamical evolution of Scattered Disk Objects (SDOs). We obtain a model of the intrinsic orbital element distribution of GPCs. The Scattered Disk represents the main source of prograde GPCs and Centaurs, while the contribution from Plutinos lies between one and two orders of magnitude below that from the SD. We obtain the number and size distribution of GPCs from our model, computing 9600 GPCs from the SD with $D > 100$ km and $\sim 10^8$ with $D > 1$ km in the current population. The contribution from other sources is considered negligible. The mean lifetime in the Centaur zone is 7.2 Myr, while the mean lifetime of SDOs in the GPC zone is of 68 Myr. The latter is dependent on the initial inclination, being the ones with high inclinations the ones that survive the longest in the GPC zone. There is also a correlation of lifetime with perihelion distance, where greater perihelion leads to longer lifetime. The dynamical evolution of observed GPCs is different for prograde and retrograde objects. Retrograde GPCs have lower median lifetime than prograde ones, thus experiencing a comparatively faster evolution. However, it is probable that this faster evolution is due to the fact that the majority of retrograde GPCs have low perihelion values and then, lower lifetimes.

Aims. We developed a six-part collisional evolution model of the main asteroid belt (MB) and used it to study the contribution of the different regions of the MB to the near-earth asteroids (NEAs). Methods. We built a statistical code called ACDC that simulates the collisional evolution of the MB split into six regions (namely Inner, Middle, Pristine, Outer, Cybele and High-Inclination belts) according to the positions of the major resonances present there ($\nu_{6}$, 3:1J, 5:2J, 7:3J and 2:1J). We consider the Yarkovsky effect and the mentioned resonances as the main mechanism that removes asteroids from the different regions of the MB and delivers them to the NEA region. We calculated the evolution of the NEAs coming from the different source regions by considering the bodies delivered by the resonances and mean dynamical timescales in the NEA population. Results. Our model is in agreement with the major observational constraints associated with the MB, such as the size distributions of the different regions of the MB and the number of large asteroid families. It is also able to reproduce the observed NEAs with H < 16 and agrees with recent estimations for H < 20, but deviates for smaller sizes. We find that most sources make a significant contribution to the NEAs; however the Inner and Middle belts stand out as the most important source of NEAs followed by the Outer belt. The contributions of the Pristine and Cybele regions are minor. The High-Inclination belt is the source of only a fraction of the actual observed NEAs with high inclination, as there are dynamical processes in that region that enable asteroids to increase and decrease their inclinations.

We report follow-up observations of transiting exoplanets that have either large uncertainties (>10 minutes) in their transit times or have not been observed for over three years. A fully robotic ground-based telescope network, observations from citizen astronomers and data from TESS have been used to study eight planets, refining their ephemeris and orbital data. Such follow-up observations are key for ensuring accurate transit times for upcoming ground and space-based telescopes which may seek to characterise the atmospheres of these planets. We find deviations from the expected transit time for all planets, with transits occurring outside the 1 sigma uncertainties for seven planets. Using the newly acquired observations, we subsequently refine their periods and reduce the current predicted ephemeris uncertainties to 0.28 - 4.01 minutes. A significant portion of this work has been completed by students at two high schools in London as part of the Original Research By Young Twinkle Students (ORBYTS) programme.

From its discovery, the WASP-18 system with its massive transiting planet on a tight orbit was identified as a unique laboratory for studies on tidal planet-star interactions. In an analysis of Doppler data, which include five new measurements obtained with the HIRES/Keck-I instrument between 2012 and 2018, we show that the radial velocity signal of the photosphere following the planetary tidal potential can be distilled for the host star. Its amplitude is in agreement with both theoretical predictions of the equilibrium tide approximation and an ellipsoidal modulation observed in an orbital phase curve. Assuming a circular orbit, we refine system parameters using photometric time series from TESS. With a new ground-based photometric observation, we extend the span of transit timing observations to 28 years in order to probe the rate of the orbital period shortening. Since we found no departure from a constant-period model, we conclude that the modified tidal quality parameter of the host star must be greater than 3.9x10^6 with 95% confidence. This result is in line with conclusions drawn from studies of the population of hot Jupiters, predicting that the efficiency of tidal dissipation is 1 or 2 orders of magnitude weaker. As the WASP-18 system is one of the prime candidates for detection of orbital decay, further timing observations are expected to push the boundaries of our knowledge on stellar interiors.

The small ($\le$ 135 km mean radius) saturnian satellites are closely related to the rings and together they constitute a complex dynamical system where formation and destruction mechanisms compete against each other. We model the cratering process by Centaur objects on the small saturnian satellites, and compare our results with the observational crater counts obtained from the Voyager and Cassini missions. Using a theoretical model previously developed we calculate the crater production on these satellites considering two slopes of the size-frequency distribution (SFD) for the smaller objects of the Centaur population and compare our results with the available observations. In addition, we consider the case of catastrophic collisions between these satellites and Centaur objects and calculate the age of formation of those satellites that suffer one or more disruptions. In general we find that the observed crater distributions are best modeled by the crater size distribution corresponding to the $s_2 = 3.5$ index of the SFD of impactors with diameters smaller than 60 km. However, for crater diameters $D \lesssim 3-8$ km (which correspond to impactor diameters $d \sim 0.04 - 0.15$ km), the observed distributions become flatter and deviate from our results, which may evidence processes of erosion and/or crater saturation at small crater sizes or a possible break in the SFD of impactors at $d \sim 0.04 - 0.15$ km to a much shallower differential slope of $\sim -1.5$. Our results suggest that Pan, Daphnis, Atlas, Aegaeon, Methone, Anthe, Pallene, Calypso and Polydeuces suffered one or more catastrophic collisions over the age of the Solar System, the younger being associated to arcs with ages of $\sim10^8$ years. We have also calculated the satellites' surface ages, which indicate ongoing resurfacing processes.

We homogeneously analyse $\sim 3.2\times 10^5$ photometric measurements for $\sim 1100$ transit lightcurves belonging to $17$ exoplanet hosts. The photometric data cover $16$ years 2004--2019 and include amateur and professional observations. Old archival lightcurves were reprocessed using up-to-date exoplanetary parameters and empirically debiased limb-darkening models. We also derive self-consistent transit and radial-velocity fits for $13$ targets. We confirm the nonlinear TTV trend in the WASP-12 data at a high significance, and with a consistent magnitude. However, Doppler data reveal hints of a radial acceleration about $(-7.5\pm 2.2)$~m/s/yr, indicating the presence of unseen distant companions, and suggesting that roughly $10$ per cent of the observed TTV was induced via the light-travel (or Roemer) effect. For WASP-4, a similar TTV trend suspected after the recent TESS observations appears controversial and model-dependent. It is not supported by our homogeneus TTV sample, including $10$ ground-based EXPANSION lightcurves obtained in 2018 simultaneously with TESS. Even if the TTV trend itself does exist in WASP-4, its magnitude and tidal nature are uncertain. Doppler data cannot entirely rule out the Roemer effect induced by possible distant companions.

In 1977, while Apple II and Atari computers were being sold, a tiny dot was observed in an inconvenient orbit. The minor body 1977 UB, to be named (2060) Chiron, with an orbit between Saturn and Uranus, became the first Centaur, a new class of minor bodies orbiting roughly between Jupiter and Neptune. The observed overabundance of short-period comets lead to the downfall of the Oort Cloud as exclusive source of comets and to the rise of the need for a Trans-Neptunian comet belt. Centaurs were rapidly seen as the transition phase between Kuiper Belt Objects (KBOs), also known as Trans-Neptunian Objects (TNOs) and the Jupiter-Family Comets (JFCs). Since then, a lot more has been discovered about Centaurs: they can have cometary activity and outbursts, satellites, and even rings. Over the past four decades since the discovery of the first Centaur, rotation periods, surface colors, reflectivity spectra and albedos have been measured and analyzed. However, despite such a large number of studies and complementary techniques, the Centaur population remains a mystery as they are in so many ways different from the TNOs and even more so from the JFCs.

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.

Recent studies have analyzed the dynamical evolution of outer small body populations under the effects of an eccentric inner massive perturber. Such outer reservoirs are composed of particles on prograde and retrograde orbits, as well as particles whose orbit flips from prograde to retrograde and back again showing a coupling between the inclination i and the ascending node longitude \Omega (Type-F particles). We analyze the role of the General Relativity (GR) on the dynamics of outer particles under the influence of an inner eccentric Jupiter-mass planet produced by a planetary scattering event. In particular, we study how the GR affects the dynamical evolution of the outer Type-F particles, which experience an eccentric Lidov-Kozai mechanism. We carry out N-body simulations with and without GR effects. When the GR is included, the extreme values of \Omega are obtained for retrograde inclinations, while the minimum and maximum inclinations allowed for Type-F particles increase in comparison with that derived without GR effects. According to this, if the GR is included in the simulations, the range of prograde (retrograde) inclinations of the libration region is reduced (increased) respect to that obtained in absence of GR. We find two new class of particles when the GR effects are included in the simulations. On the one hand, particles whose orbital plane flips from prograde to retrograde and back again without experiencing a coupling between i and \Omega. On the other hand, retrograde particles that show a strong coupling between i and \Omega. We infer that the GR may significantly modify the dynamical properties of the outer reservoirs that evolve under the effects of an eccentric inner perturber.

We analyze the process of planetary scattering around M0-type stars. To do this, we carry out N-body simulations with three Jupiter-mass planets close to their instability limit together with an outer planetesimal disk. This paper focuses on the analysis of systems in which a single Jupiter-mass planet survives after the dynamical instability event. The small body reservoirs show different dynamical behaviors. In fact, our simulations produce particles on prograde and retrograde orbits, as well as particles whose orbital plane flips from prograde to retrograde and back again along their evolution. Such particles are called "Type-F particles". We find strong correlations between the inclination $i$ and the ascending node longitude $\Omega$ of such particles. First, $\Omega$ librates around 90$^{\circ}$ or/and 270$^{\circ}$. This property is very important since it represents a necessary and sufficient condition for the flipping of an orbit. Moreover, the libration periods of $i$ and $\Omega$ are equal and they are out to phase by a quarter period. We also remark that the larger the libration amplitude of $i$, the larger the libration amplitude of $\Omega$. Finally, we analyze the initial conditions of Type-F particles of all our simulations immediately after the dynamical instability event, when a single Jupiter-mass planet survives in the system. We carry out this study with the goal to determine the parameter space that lead to the flipping of an orbit. Our results suggest that the orbit of a test particle can flip for any value of its initial eccentricity, although we found only two Type-F particles with initial inclinations $i <$ 17$^{\circ}$. Moreover, our study indicates that the minimum value of the inclination of the Type-F particles in a given system decreases with an increase in the eccentricity of the giant planet.

The secular approximation of the hierarchical three body systems has been proven to be very useful in addressing many astrophysical systems, from planets, stars to black holes. In such a system two objects are on a tight orbit, and the tertiary is on a much wider orbit. Here we study the dynamics of a system by taking the tertiary mass to zero and solve the hierarchical three body system up to the octupole level of approximation. We find a rich dynamics that the outer orbit undergoes due to gravitational perturbations from the inner binary. The nominal result of the precession of the nodes is mostly limited for the lowest order of approximation, however, when the octupole-level of approximation is introduced the system becomes chaotic, as expected, and the tertiary oscillates below and above 90deg, similarly to the non-test particle flip behavior (e.g., Naoz 2016). We provide the Hamiltonian of the system and investigate the dynamics of the system from the quadrupole to the octupole level of approximations. We also analyze the chaotic and quasi-periodic orbital evolution by studying the surfaces of sections. Furthermore, including general relativity, we show case the long term evolution of individual debris disk particles under the influence of a far away interior eccentric planet. We show that this dynamics can naturally result in retrograde objects and a puffy disk after a long timescale evolution (few Gyr) for initially aligned configuration.

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.

The observations of the surfaces of the mid sized Saturnian satellites made by Cassini Huygens mission have shown a variety of features that allows study of the processes that took place and are taking place on those worlds. Research of the Saturnian satellite surfaces has clear implications for Saturn history and surroundings. In a recent paper, the production of craters on the mid sized Saturnian satellites by Centaur objects was calculated considering the current Solar System. We have compared our results with crater counts from Cassini images and we have noted that the number of observed small craters is less than our calculated number. In this paper we estimate the age of the surface for each observed terrain on each mid sized satellite of Saturn. We have noticed that since there are less observed small craters than calculated (except on Iapetus), this results in younger ages. This could be the result of efficient endogenous or exogenous process(es) for erasing small craters and or crater saturation at those sizes. The size limit from which the observed number of smaller craters is less than the calculated is different for each satellite, possibly indicating processes that are unique to each, but other potential common explanations would be crater saturation and or deposition of E ring particles. These processes are also suggested by the findings that the smaller craters are being preferentially removed, and the erasure process is gradual. On Enceladus, only mid and high latitude plains have remnants of old terrains; the other regions could be young; the regions near the South Polar Terrain could be as young as 50 Myr old. On the contrary for Iapetus, all the surface is old and it notably registers a primordial source of craters. As the crater size is decreased, it would be perceived to approach saturation until D less than 2 km craters, where saturation is complete.

We study the impactor flux and cratering on Ceres and Vesta caused by the collisional and dynamical evolution of the asteroid Main Belt. We develop a statistical code based on a well-tested model for the simultaneous evolution of the Main Belt and NEA size distributions. This code includes catastrophic collisions and noncollisional removal processes such as the Yarkovsky effect and the orbital resonances. The model assumes that the dynamical depletion of the early Main Belt was very strong, and owing to that, most Main Belt comminution occurred when its dynamical structure was similar to the present one. Our results indicate that the number of D > 1 km Main Belt asteroids striking Ceres and Vesta over the Solar System history are approximately 4 600 and 1 100 respectively. The largest Main Belt asteroids expected to have impacted Ceres and Vesta had diameters of 71.7 km and 21.1 km. The number of D > 0.1 km craters on Ceres is \sim 3.4 \times 10^8 and 6.2 \times 10^7 on Vesta. The number of craters with D > 100 km are 47 on Ceres and 8 on Vesta. Our study indicates that the D = 460 km crater observed on Vesta had to be formed by the impact of a D \sim 66.2 km projectile, which has a probability of occurr \sim 30% over the Solar System history. If significant discrepancies between our results about the cratering on Ceres and Vesta and data obtained from the Dawn Mission were found, they should be linked to a higher degree of collisional evolution during the early Main Belt and/or the existence of the late heavy bombardment. An increase in the collisional activity in the early phase may be provided for an initial configuration of the giant planets consistent with, for example, the Nice model. From this, the Dawn Mission would be able to give us clues about the initial configuration of the early Solar System and its subsequent dynamical evolution.

Phoebe is one of the irregular satellites of Saturn; the images taken by Cassini-Huygens spacecraft allowed us to analyze its surface and the craters on it. We study the craters on Phoebe produced by Centaur objects from the Scattered Disk (SD) and plutinos escaped from the 3:2 mean motion resonance with Neptune and compare our results with the observations by Cassini. We use previous simulations on trans-Neptunian Objects and a method that allows us to obtain the number of craters and the cratering rate on Phoebe. We obtain the number of craters and the greatest crater on Phoebe produced by Centaurs in the present configuration of the Solar System. Moreover, we obtain a present normalized rate of encounters of Centaurs with Saturn of $\dot F = 7.1 \times 10^{-11}$ per year, from which we can infer the current cratering rate on Phoebe for each crater diameter. Our study and the comparison with the observations suggest that the main crater features on Phoebe are unlikely to have been produced in the present configuration of the Solar System and that they must have been acquired when the SD were depleted in the early Solar System. If this is what happened and the craters were produced when Phoebe was a satellite of Saturn, then it had to be captured, very early in the evolution of the Solar System.

We study the impactor flux and cratering on Pluto and Charon due to the collisional evolution of Plutinos. Plutinos are those trans-Neptunian objects located at 39.5 AU, in the 3:2 mean motion resonance with Neptune. To do this, we develop a statistical code that includes catastrophic collisions and cratering events, and takes into account the stability and instability zones of the 3:2 mean motion resonance with Neptune. We proposes different initial populations that account for the uncertainty in the size distribution of Plutinos at small sizes. Depending on the initial population, our results indicate the following. The number of D > 1 km Plutinos streaking Pluto over 3.5 Gyr is between 1271 and 5552. For Charon, the number of D > 1 km Plutino impactors is between 354 and 1545. The number of D > 1 km craters on Pluto produced by Plutinos during the last 3.5 Gyr is between 43076 and 113879. For Charon, the number of D > 1 km craters is between 20351 and 50688.The largest Plutino impactor onto Pluto has a diameter between 17 and 23 km, which produces a craterwith a diameter of 31 - 39 km. The largest Plutino impactor onto Charon has a diameter between 10 and 15 km, which produces a crater with a diameter of 24 - 33 km. We test if 2 Pluto-sized objects are assumed in the 3:2 Neptune resonance, then the total number of Plutino impactors both onto Pluto as Charon with diameters D > 1 km is a factor of 1.6 - 1.8 larger if considering 1 Pluto-sized object. Given the dynamical structure of the trans-Neptunian region, it is necessary to study in detail the contribution of all the potential sources of impactors on the Pluto-Charon system, to obtain the main contributor and the whole production of craters. Then, we will be able to contrast those studies with observations which will help us to understand the geological processes and history of the surface of those worlds.

It was shown in previous works the existence of weakly chaotic orbits in the plutino population that diffuse very slowly. These orbits correspond to long-term plutino escapers and then represent the plutinos that are escaping from the resonance at present. In this paper we perform numerical simulations in order to explore the dynamical evolution of plutinos recently escaped from the resonance. The numerical simulations were divided in two parts. In the first one we evolved 20,000 test particles in the resonance in order to detect and select the long-term escapers. In the second one, we numerically integrate the selected escaped plutinos in order to study their dynamical post escaped behavior. Our main results include the characterization of the routes of escape of plutinos and their evolution in the Centaur zone. We obtained a present rate of escape of plutinos between 1 and 10 every 10 years. The escaped plutinos have a mean lifetime in the Centaur zone of 108 Myr and their contribution to the Centaur population would be a fraction of less than 6 % of the total Centaur population. In this way, escaped plutinos would be a secondary source of Centaurs.