We characterize the radii and masses of the star and planets in the Kepler-59 system, as well as their orbital parameters. The star parameters are determined through a standard spectroscopic analysis, resulting in a mass of $1.359\pm 0.155\,M_\odot$ and a radius of $1.367\pm 0.078\,R_\odot$. The planetary radii obtained are $1.5\pm 0.1\,R_\oplus$ for the inner and $2.2\pm 0.1\,R_\oplus$ for the outer planet. The orbital parameters and the planetary masses are determined by the inversion of Transit Timing Variations (TTV) signals. For this, we consider two different data sets, one provided by Holczer et al. 2016, with TTVs only for the planet Kepler-59c, and the other provided by Rowe et al. 2015, with TTVs signals for both planets. The inversion method is carried out by applying an algorithm of Bayesian inference (MultiNest) combined with an efficient N-body integrator (Swift). For each of the data sets, two possible solutions are found, both having the same probability according to their corresponding Bayesian evidences. All four solutions appear to be indistinguishable within their 2-$\sigma$ uncertainties. Nevertheless, statistical analyses show that the solutions from Rowe et al. 2015 data better characterize the data. The first and second solutions identify masses of $5_{-2}^{+4}~M_{\mathrm{\oplus}}$ and $4.6_{-2.0}^{+3.6}~M_{\mathrm{\oplus}}$, and $3.0^{+0.8}_{-0.8}~M_{\mathrm{\oplus}}$ and $2.6^{+1.9}_{-0.8}~M_{\mathrm{\oplus}}$ for the inner and outer planet, respectively. This points to a system with an inner super-Earth and an outer mini-Neptune. Dynamical studies show the planets have almost co-planar orbits with small eccentricities ($e<0.1$), close but not into the 3:2 mean motion resonance. Stability analysis indicates that this configuration is stable over million years of evolution.

We present the discovery of a pair of transiting giant planets, TOI-216b and c, using four sectors of TESS photometry. TOI-216 is a $0.87 M_{\odot}$ dwarf orbited by two transiters with radii of $8.2 R_{\oplus}$ and $11.3 R_{\oplus}$, and periods of $17.01$d and $34.57$d, respectively. Anti-correlated TTVs are clearly evident indicating that the transiters orbit the same star and interact via a near 2:1 mean motion resonance. By fitting the TTVs with a dynamical model, we infer masses of $26_{-11}^{+24} M_{\oplus}$ and $190_{-80}^{+220} M_{\oplus}$, establishing that the objects are planetary in nature and have likely sub-Kronian and Kronian densities. TOI-216 lies close to the southern ecliptic pole and thus will be observed by TESS throughout the first year, providing an opportunity for continuous dynamical monitoring and considerable refinement of the dynamical masses presented here. TOI-216 closely resembles Kepler-9 in architecture, and we hypothesize that in such systems these Saturn-analogs failed to fully open a gap and thus migrated far deeper into the system before becoming trapped into resonance, which would imply that future detections of new analogs may also have sub-Jupiter masses.

We develop an evolutionary model of the long-period comet (LPC) population, starting from their birthplace in a massive trans-Neptunian disk that was dispersed by migrating giant planets. Most comets that remain bound to the Solar system are stored in the Oort cloud. Galactic tides and passing stars make some of these bodies evolve into observable comets in the inner Solar system. Our approach models each step in a full-fledged numerical framework. Subsequent analysis consists of applying plausible fading models and computing the original orbits to compare with observations. Our results match the observed semimajor axis distribution of LPCs when Whipple's power-law fading scheme with an exponent $\kappa= 0.6^{+0.1}_{-0.2}$ is adopted. The cumulative perihelion ($q$) distribution is fit well by a linear increase plus a weak quadratic term. Beyond $q = 15$~au, however, the population increases steeply and the isotropy of LPC orbital planes breaks. We find tentative evidence from the perihelion distribution of LPCs that the returning comets are depleted in supervolatiles and become active due to water ice sublimation for $q\leq 3$ au. Using an independent calibration of the population of the initial disk, our predicted LPC flux is smaller than observations suggest by a factor of $\simeq 2$. Current data only characterize comets from the outer Oort cloud (semimajor axes $\gtrsim 10^4$~au). A true boost in understanding the Oort cloud's structure should result from future surveys when they detect LPCs with perihelia beyond $15$~au. Our results provide observational predictions of what can be expected from these new data.

Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydrocode. We find that embryos in a disk with high surface density ($\Sigma_0 = 990\,{\rm g}\,{\rm cm}^{-2}$) undergo `unsuccessful' two-body encounters which do not lead to a merger. Only when a 3rd protoplanet arrives to the convergence zone, three-body encounters lead to mergers. For a low-viscosity disk ($\nu = 5\times10^{13}\,{\rm cm}^2\,{\rm s}^{-1}$) a massive coorbital is a possible outcome, for which a pebble isolation develops and the coorbital is stabilised. For more massive protoplanets ($5\,M_\oplus$), the convergence radius is located further out, in the ice-giant zone. After a series of encounters, there is an evolution driven by a dynamical torque of a tadpole region, which is systematically repeated several times, until the coorbital configuration is disrupted and planets merge. This may be a pathway how to solve the problem that coorbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets ($0.1\,M_\oplus$) reveals completely different dynamics. The evolution is no longer smooth, but rather a random walk. This is because the spiral arms, developed in the gas disk due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surroundings. Our hydrodynamical simulations may have important implications for N-body simulations of planetary migration that use simplified torque prescriptions and are thus unable to capture protoplanet dynamics in its full glory.

We present radiation hydrodynamic simulations in which binary planets form by close encounters in a system of several super-Earth embryos. The embryos are embedded in a protoplanetary disk consisting of gas and pebbles and evolve in a region where the disk structure supports convergent migration due to Type I torques. As the embryos accrete pebbles, they become heated and thus affected by the thermal torque (Ben\'{i}tez-Llambay et al. 2015) and the hot-trail effect (Chrenko et al. 2017) which excites orbital eccentricities. Motivated by findings of Eklund & Masset (2017), we assume the hot-trail effect operates also vertically and reduces the efficiency of inclination damping. Non-zero inclinations allow the embryos to become closely packed and also vertically stirred within the convergence zone. Subsequently, close encounters of two embryos assisted by the disk gravity can form transient binary planets which quickly dissolve. Binary planets with a longer lifetime $\sim$$10^{4}$ yr form in 3-body interactions of a transient pair with one of the remaining embryos. The separation of binary components generally decreases in subsequent encounters and due to pebble accretion until the binary merges, forming a giant planet core. We provide an order-of-magnitude estimate of the expected occurrence rate of binary planets, yielding one binary planet per $\simeq$$2$--$5\times10^{4}$ planetary systems. Therefore, although rare, the binary planets may exist in exoplanetary systems and they should be systematically searched for.

Using numerical integrations, we find that the orbital eccentricity of Saturn's moon Iapetus undergoes prominent multi-Myr oscillations. We identify the responsible resonant argument to be $\varpi-\varpi_{g5}+\Omega-\Omega_{eq}$, with the terms being the longitudes of pericenter of Iapetus and planetary secular mode $g_5$, Iapetus's longitude of the node and Saturn's equinox. We find that this argument currently (on a $10^7$~yr timescale) appears to librate with a very large amplitude. On longer timescales, the behavior of this resonant angle is strongly dependent on the resonant interaction between Saturn's spin axis and the planetary mode $f_8$, with long-term secular resonance being possible if Saturn's equinox is librating relative to the node of the $f_8$ eigenmode. We present analytical estimates of the dependence of the resonant argument on the orbital elements of Iapetus. We find that this Iapetus-$g_5$ secular resonance could have been established only after the passage of Iapetus through the 5:1 mean-motion resonance with Titan, possibly in the last Gyr. Using numerical simulations, we show that the capture into the secular resonance appears to be a low-probability event. While the Iapetus-$g_5$ secular resonance can potentially help us put new constraints on the past dynamics of the Saturnian system, uncertainties in both the spin axis dynamics of Saturn and the tidal evolution rate of Titan make it impossible to make any firm conclusions about the resonance's longevity and origin.

Jupiter and Saturn each have complex systems of satellites and rings. These satellites can be classified into dynamical groups, implying similar formation scenarios. Recently, a larger number of additional irregular satellites have been discovered around both gas giants that have yet to be classified. The aim of this paper is to examine the relationships between the satellites and rings of the gas giants, using an analytical technique called cladistics. Cladistics is traditionally used to examine relationships between living organisms, the `tree of life'. In this work, we perform the first cladistical study of objects in a planetary science context. Our method uses the orbital, physical and compositional characteristics of satellites to classify the objects in the Jovian and Saturnian systems. We find that the major relationships between the satellites in the two systems, such as families, as presented in previous studies, are broadly preserved. In addition, based on our analysis of the Jovian system, we identify a new retrograde irregular family, the Iocaste family, and suggest that the Phoebe family of the Saturnian system can be further divided into two subfamilies. We also propose that the Saturnian irregular families be renamed, to be consistent with the convention used in Jovian families. Using cladistics, we are also able to assign the new unclassified irregular satellites into families. Taken together, the results of this study demonstrate the potential use of the cladistical technique in the investigation of relationships between orbital bodies.

We report on our study of asteroidal breakups, i.e. fragmentations of targets, subsequent gravitational reaccumulation and formation of small asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$ km. Simulations were performed with a smoothed-particle hydrodynamics (SPH) code combined with an efficient N-body integrator. We assumed various projectile sizes, impact velocities and impact angles (125 runs in total). Resulting size-frequency distributions are significantly different from scaled-down simulations with $D_{\rm pb} = 100$ km targets (Durda et al., 2007). We derive new parametric relations describing fragment distributions, suitable for Monte-Carlo collisional models. We also characterize velocity fields and angular distributions of fragments, which can be used as initial conditions for N-body simulations of small asteroid families. Finally, we discuss a number of uncertainties related to SPH simulations.

The age of a young asteroid family can be determined by tracking the orbits of family members backward in time and showing that they converge at some time in the past. Here we consider the Veritas family. We find that the membership of the Veritas family increased enormously since the last detailed analysis of the family. Using backward integration, we confirm the convergence of nodal longitudes $\Omega$, and, for the first time, also obtain a simultaneous convergence of pericenter longitudes $\varpi$. The Veritas family is found to be $8.23^{+0.37}_{-0.31}$~Myr old. To obtain a tight convergence of $\Omega$ and $\varpi$, as expected from low ejection speeds of fragments, the Yarkovsky effect needs to be included in the modeling of the past orbital histories of Veritas family members. Using this method, we compute the Yarkovsky semi-major axis drift rates, ${\rm d}a/{\rm d}t$, for 274 member asteroids. The distribution of ${\rm d}a/{\rm d}t$ values is consistent with a population of C-type objects with low densities and low thermal conductivities. The accuracy of individual ${\rm d}a/{\rm d}t$ measurements is limited by the effect of close encounters of member asteroids to (1) Ceres and other massive asteroids, which cannot be evaluated with confidence.

The Hungaria asteroid group is located interior to the main asteroid belt, with semimajor axes between 1.8 and 2 AU, low eccentricities and inclinations of 16-35 degrees. Recently, it has been proposed that Hungaria asteroids are a secularly declining population that may be related to the Late Heavy Bombardment (LHB) impactors (\'Cuk et al. 2012, Bottke et al. 2012). While \'Cuk et al. (2012) and Bottke et al. (2012) have reproduced a Hungaria-like population that declined exponentially, the real Hungarias were never confirmed to be unstable to the same degree. Here we find that the stability of Hungarias is strongly dependent on the evolution of the eccentricity of Mars, which is chaotic and unpredictable on Gyr timescales. We find that the high Martian eccentricity chiefly affects Hungarias through close approaches with Mars, rather than planetary secular modes. However, current minimum perihelia of Hungarias (over Myr timescales) are not diagnostic of their long-term stability due to a number of secular and mean motion resonances affecting the Hungaria region (Milani et al. 2010). We conclude that planetary chaos makes it impossible to determine the effective lifetimes of observed Hungarias. Furthermore, long-term changes of Martian eccentricity could lead to variable Hungaria loss over time. We speculate that some of the most stable Hungarias may have been placed in their present orbit when the eccentricity of Mars was significantly higher than today.

We use 16 quarters of the \textit{Kepler} mission data to analyze the transit timing variations (TTVs) of the extrasolar planet Kepler-46b (KOI-872). Our dynamical fits confirm that the TTVs of this planet (period $P=33.648^{+0.004}_{-0.005}$ days) are produced by a non-transiting planet Kepler-46c ($P=57.325^{+0.116}_{-0.098}$ days). The Bayesian inference tool \texttt{MultiNest} is used to infer the dynamical parameters of Kepler-46b and Kepler-46c. We find that the two planets have nearly coplanar and circular orbits, with eccentricities $\simeq 0.03$ somewhat higher than previously estimated. The masses of the two planets are found to be $M_{b}=0.885^{+0.374}_{-0.343}$ and $M_{c}=0.362^{+0.016}_{-0.016}$ Jupiter masses, with $M_{b}$ being determined here from TTVs for the first time. Due to the precession of its orbital plane, Kepler-46c should start transiting its host star in a few decades from now.

V-type asteroids are a taxonomic class whose surface is associated to a basaltic composition. The only known source of V-type asteroids in the Main Asteroid Belt is (4) Vesta, that is located in the inner part of the belt. However, many V-type asteroids cannot be dynamically linked to Vesta., in particular, those asteroids located in the middle and outer parts of the Belt. Previous works have failed to find mechanisms to transport V-type asteroids from the inner to the middle and outer belt. In this work we propose a dynamical mechanism that could have acted on primordial asteroid families. We consider a model of the giant planets migration known as the jumping Jupiter model with five planets. Our study is focused on the period of 10 Myr that encompasses the instability phase of the giant planets. We show that, for different hypothetical Vesta-like paleo-families in the inner belt, the perturbations caused by the ice giant that is scattered into the asteroid belt before being ejected from the solar system, are able to scatter V-type asteroids to the middle and outer belt. Based on the orbital distribution of V-type candidates identified from the Sloan Digital Sky Survey and the VISTA Survey colours, we show that this mechanism is efficient enough provided that the hypothetical paleo-family originated from a 100 to 500 km crater excavated on the surface of (4) Vesta. This mechanism is able to explain the currently observed V-type asteroids in the middle and outer belt, with the exception of (1459) Magnya.

In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Bro\v{z} and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body.

The Karin cluster is a young asteroid family thought to have formed only $\simeq 5.75$~My ago. The young age can be demonstrated by numerically integrating the orbits of Karin cluster members backward in time and showing the convergence of the perihelion and nodal longitudes (as well as other orbital elements). Previous work has pointed out that the convergence is not ideal if the backward integration only accounts for the gravitational perturbations from the Solar System planets. It improves when the thermal radiation force known as the Yarkovsky effect it is accounted for. This argument can be used to estimate the spin obliquities of the Karin cluster members. Here we take advantage of the fast growing membership of the Karin cluster and show that the obliquity distribution of diameter $D\simeq 1-2$ km Karin asteroids is bimodal, as expected if the YORP effect acted to move obliquities toward the extreme values ($0^\circ$ or $180^\circ$). The measured magnitude of the effect is consistent with the standard YORP model. The surface thermal conductivity is inferred to be $0.07$-0.2 W m$^{-1}$ K$^{-1}$ (thermal inertia $\simeq 300-500$ J m$^{-2}$ K$^{-1}$s$^{-1/2}$). We find that the strength of the YORP effect is roughly $\simeq 0.7$ of the nominal strength obtained for a collection of random Gaussian spheroids. These results are consistent with a surface composed of rough, rocky regolith. The obliquity values predicted here for 480 members of the Karin cluster can be validated by the light-curve inversion method.

We explore the past evolution of Saturn's moons using direct numerical integrations. We find that the past Tethys-Dione 3:2 orbital resonance predicted in standard models likely did not occur, implying that the system is less evolved than previously thought. On the other hand, the orbital inclinations of Tethys, Dione and Rhea suggest that the system did cross the Dione-Rhea 5:3 resonance, which is closely followed by a Tethys-Dione secular resonance. A clear implication is that either the moons are significantly younger than the planet, or that their tidal evolution must be extremely slow (Q > 80,000). As an extremely slow-evolving system is incompatible with intense tidal heating of Enceladus, we conclude that the moons interior to Titan are not primordial, and we present a plausible scenario for the system's recent formation. We propose that the mid-sized moons re-accreted from a disk about 100 Myr ago, during which time Titan acquired its significant orbital eccentricity. We speculate that this disk has formed through orbital instability and massive collisions involving the previous generation of Saturn's mid-sized moons. We identify the solar evection resonance perturbing a pair of mid-sized moons as the most likely trigger of such an instability. This scenario implies that most craters on the moons interior to Titan must have been formed by planetocentric impactors.

The orbit of Mercury has large values of eccentricity and inclination that cannot be easily explained if this planet formed on a circular and coplanar orbit. Here, we study the evolution of Mercury's orbit during the instability related to the migration of the giant planets in the framework of the jumping Jupiter model. We found that some instability models are able to produce the correct values of Mercury's eccentricity and inclination, provided that relativistic effects are included in the precession of Mercury's perihelion. The orbital excitation is driven by the fast change of the normal oscillation modes of the system corresponding to the perihelion precession of Jupiter (for the eccentricity), and the nodal regression of Uranus (for the inclination).

Asteroid families are groups of minor bodies produced by high-velocity collisions. After the initial dispersions of the parent bodies fragments, their orbits evolve because of several gravitational and non-gravitational effects,such as diffusion in mean-motion resonances, Yarkovsky and YORP effects, close encounters of collisions, etc. The subsequent dynamical evolution of asteroid family members may cause some of the original fragments to travel beyond the conventional limits of the asteroid family. Eventually, the whole family will dynamically disperse and no longer be recognizable. A natural question that may arise concerns the timescales for dispersion of large families. In particular, what is the oldest still recognizable family in the main belt? Are there any families that may date from the late stages of the Late Heavy Bombardment and that could provide clues on our understanding of the primitive Solar System? In this work, we investigate the dynamical stability of seven of the allegedly oldest families in the asteroid main belt. Our results show that none of the seven studied families has a nominally mean estimated age older than 2.7 Gyr, assuming standard values for the parameters describing the strength of the Yarkovsky force. Most "paleo-families" that formed between 2.7 and 3.8 Gyr would be characterized by a very shallow size-frequency distribution, and could be recognizable only if located in a dynamically less active region (such as that of the Koronis family). V-type asteroids in the central main belt could be compatible with a formation from a paleo-Eunomia family.

Ceres is the largest and most massive body in the asteroid main belt. Observational data from the Dawn spacecraft reveal the presence of at least two impact craters about 280~km in diameter on the Ceres surface, that could have expelled a significant number of fragments. Yet, standard techniques for identifying dynamical asteroid families have not detected any Ceres family. In this work, we argue that linear secular resonances with Ceres deplete the population of objects near Ceres. Also, because of the high escape velocity from Ceres, family members are expected to be very dispersed, with a considerable fraction of km-sized fragments that should be able to reach the pristine region of the main belt, the area between the 5J:-2A and 7J:-3A mean-motion resonances, where the observed number of asteroids is low. Rather than looking for possible Ceres family members near Ceres, here we propose to search in the pristine region. We identified 156 asteroids whose taxonomy, colors, albedo could be compatible with being fragments from Ceres. Remarkably, most of these objects have inclinations near that of Ceres itself.

An asteroid family forms as a result of a collision between an impactor and a parent body. The fragments with ejection speeds higher than the escape velocity from the parent body can escape its gravitational pull. The cloud of escaping debris can be identified by the proximity of orbits in proper element, or frequency, domains. Obtaining estimates of the original ejection speed can provide valuable constraints on the physical processes occurring during collision, and used to calibrate impact simulations. Unfortunately, proper elements of asteroids families are modified by gravitational and non-gravitational effects, such as resonant dynamics, encounters with massive bodies, and the Yarkovsky effect, such that information on the original ejection speeds is often lost, especially for older, more evolved families. It has been recently suggested that the distribution in proper inclination of the Koronis family may have not been significantly perturbed by local dynamics, and that information on the component of the ejection velocity that is perpendicular to the orbital plane ($v_W$), may still be available, at least in part. In this work we estimate the magnitude of the original ejection velocity speeds of Koronis members using the observed distribution in proper eccentricity and inclination, and accounting for the spread caused by dynamical effects. Our results show that i) the spread in the original ejection speeds is, to within a 15% error, inversely proportional to the fragment size, and ii) the minimum ejection velocity is of the order of 50 m/s, with larger values possible depending on the orbital configuration at the break-up.

Asteroid families form as a result of large-scale collisions among main belt asteroids. The orbital distribution of fragments after a family-forming impact could inform us about their ejection velocities. Unfortunately, however, orbits dynamically evolve by a number of effects, including the Yarkovsky drift, chaotic diffusion, and gravitational encounters with massive asteroids, such that it is difficult to infer the ejection velocities eons after each family's formation. Here we analyze the inclination distribution of asteroid families, because proper inclination can remain constant over long time intervals, and could help us to understand the distribution of the component of the ejection velocity that is perpendicular to the orbital plane ($v_{W}$). From modeling the initial breakup, we find that the distribution of $v_{W}$ of the fragments, which manage to escape the parent body's gravity, should be more peaked than a Gaussian distribution (i.e., be leptokurtic) even if the initial distribution was Gaussian. We surveyed known asteroid families for signs of a peaked distribution of $v_{W}$ using a statistical measure of the distribution peakedness or flatness known as kurtosis. We identified eight families whose $v_{W}$ distribution is significantly leptokurtic. These cases (e.g. the Koronis family) are located in dynamically quiet regions of the main belt, where, presumably, the initial distribution of $v_{W}$ was not modified by subsequent orbital evolution. We suggest that, in these cases, the inclination distribution can be used to obtain interesting information about the original ejection velocity field.

In this work, we investigate the evolution of a primordial belt of asteroids, represented by a large number of massless test particles, under the gravitational effect of migrating Jovian planets in the framework of the jumping-Jupiter model. We perform several simulations considering test particles distributed in the Main Belt, as well as in the Hilda and Trojan groups. The simulations start with Jupiter and Saturn locked in the mutual 3:2 mean motion resonance plus 3 Neptune-mass planets in a compact orbital configuration. Mutual planetary interactions during migration led one of the Neptunes to be ejected in less than 10 Myr of evolution, causing Jupiter to jump by about 0.3 au in semi-major axis. This introduces a large scale instability in the studied populations of small bodies. After the migration phase, the simulations are extended over 4 Gyr, and we compare the final orbital structure of the simulated test particles to the current Main Belt of asteroids with absolute magnitude $H<9.7$. The results indicate that, in order to reproduce the present Main Belt, the primordial belt should have had a distribution peaked at $\sim10^{\circ}$ in inclination and at $\sim0.1$ in eccentricity. We discuss the implications of this for the Grand Tack model. The results also indicate that neither primordial Hildas, nor Trojans, survive the instability, confirming the idea that such populations must have been implanted from other sources. In particular, we address the possibility of implantation of Hildas and Trojans from the Main Belt population, but find that this contribution should be minor.

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.

The Cybele region, located between the 2J:-1A and 5J:-3A mean-motion resonances, is adjacent and exterior to the asteroid main belt. An increasing density of three-body resonances makes the region between the Cybele and Hilda populations dynamically unstable, so that the Cybele zone could be considered the last outpost of an extended main belt. The presence of binary asteroids with large primaries and small secondaries suggested that asteroid families should be found in this region, but only relatively recently the first dynamical groups were identified in this area. Among these, the Sylvia group has been proposed to be one of the oldest families in the extended main belt. In this work we identify families in the Cybele region in the context of the local dynamics and non-gravitational forces such as the Yarkovsky and stochastic YORP effects. We confirm the detection of the new Helga group at $\simeq$3.65~AU, that could extend the outer boundary of the Cybele region up to the 5J:-3A mean-motion resonance. We obtain age estimates for the four families, Sylvia, Huberta, Ulla and Helga, currently detectable in the Cybele region, using Monte Carlo methods that include the effects of stochastic YORP and variability of the Solar luminosity. The Sylvia family should be $T = 1220 \pm 40$ Myr old, with a possible older secondary solution. Any collisional Cybele group formed prior to the late heavy bombardment would have been most likely completely dispersed in the jumping Jupiter scenario of planetary migration.

The Hungaria asteroids are interior to the main asteroid belt, with semimajor axes between 1.8 and 2 AU, low eccentricities and inclinations of 16-35 degrees. Small asteroids in the Hungaria region are dominated by a collisional family associated with (434) Hungaria. The dominant spectral type of the Hungaria group is the E or X-type (Warner et al, 2009), mostly due to the E-type composition of Hungaria and its genetic family. It is widely believed the E-type asteroids are related to the aubrite meteorites, also known as enstatite achondrites (Gaffey et al, 1992). Here we explore the hypothesis that aubrites originate in the Hungaria family. In order to test this connection, we compare model Cosmic Ray Exposure ages from orbital integrations of model meteoroids with those of aubrites. We show that long CRE ages of aubrites (longest among stony meteorite groups) reflect the delivery route of meteoroids from Hungarias to Earth being different than those from main-belt asteroids. We find that the meteoroids from Hungarias predominantly reach Earth by Yarkovsky-drifting across the orbit of Mars, with no assistance from orbital resonances. We conclude that the CRE ages of aubrites are fully consistent with a dominant source at the inner boundary of the Hungaria family at 1.7 AU. From here, meteoroids reach Earth through the Mars-crossing region, with relatively quick delivery times favored due to collisions (with Hungarias and the inner main-belt objects). We find that, after Vesta, (434) Hungaria is the best candidate for an asteroidal source of an achondrite group.

Lainey et al. (2012), by re-analyzing long-baseline astrometry of Saturn's moons, have found that the moons' tidal evolution is much faster than previously thought, implying an order of magnitude stronger tidal dissipation within Saturn. This result is controversial and implies recent formation of at least some of the mid-sized icy moons of Saturn. Here we show that this more intensive tidal dissipation is in full agreement with the evolved state of the Titan-Hyperion resonance. This resonance was previously thought to be non-tidal in origin, as the amount of tidal evolution required for its assembly is beyond what is possible in models that assume that all the major moons are primordial. We find that the survival of the Titan-Hyperion resonance is in agreement with a past Titan-Iapetus 5:1 resonance, but not with unbroken tidal evolution of Rhea from the rings to its current distance.

Olivine-dominated asteroids are a rare type of objects formed either in nebular processes or through magmatic differentiation. The analysis of meteorite samples suggest that at least 100 parent bodies in the main belt experienced partial or complete melting and differentiation before being disrupted. However, only a few olivine-dominated asteroids, representative of the mantle of disrupted differentiated bodies, are known to exist. Due to the paucity of these objects in the main belt their origin and evolution have been a matter of great debate over the years. In this work we present a detailed mineralogical analysis of twelve olivine-dominated asteroids. Within our sample we distinguish two classes, one that we call monomineralic-olivine asteroids and another referred to as olivine-rich asteroids. For the monomineralic-olivine asteroids the olivine chemistry was found to range from ~ Fo49 to Fo70, consistent with the values measured for brachinites and R chondrites. In the case of the olivine-rich asteroids we determined their olivine and low-Ca pyroxene abundance using a new set of spectral calibrations derived from the analysis of R chondrites spectra. We found that the olivine abundance for these asteroids varies from 0.68 to 0.93, while the fraction of low-Ca pyroxene to total pyroxene ranges from 0.6 to 0.9. A search for dynamical connections between the olivine-dominated asteroids and asteroid families found no genetic link (of the type core-mantel-crust) between these objects.

The dwarf planet Haumea has two large satellites, Namaka and Hi'iaka, which orbit at relatively large separations. Both moons have significant eccentricities and inclinations, in a pattern that is consistent with a past orbital resonance (Ragozzine and Brown, 2009). Based on our analysis, we find that the present system is not consistent with satellite formation close to the primary and tidal evolution though mean-motion resonances. We propose that Namaka experienced only limited tidal evolution, leading to the mutual 8:3 mean-motion resonance which redistributed eccentricities and inclinations between the moons. This scenario requires that the original orbit of Hi'iaka was mildly eccentric; we propose that this eccentricity was either primordial or acquired though encounters with other TNOs. Both dynamical stability and our preferred tidal evolution model imply that the moons' masses are only about one half of previously estimated values, suggesting high albedos and low densities. As the present orbits of the moons strongly suggest formation from a flat disk close to their present locations, we conclude that Hi'iaka and Namaka may be second-generation moons, formed after the breakup of a past large moon, previously proposed as the parent body of the Haumea family (Schlichting and Sari, 2009). We derive plausible parameters of that moon, consistent with the current models of Haumea's formation (Leinhardt et al, 2010). An interesting implication of this hypothesis is that Hi'iaka and Namaka may orbit retrograde with respect to Haumea's spin. Retrograde orbits of Haumea's moons would be in full agreement with available observations and our dynamical analysis, and could provide a unique confirmation the "disrupted satellite" scenario for the origin of the family.

Previous works have identified families halos by an analysis in proper elements domains, or by using Sloan Digital Sky Survey-Moving Object Catalog data, fourth release (SDSS-MOC4) multi-band photometry to infer the asteroid taxonomy, or by a combination of the two methods. The limited number of asteroids for which geometric albedo was known until recently discouraged in the past the extensive use of this additional parameter, which is however of great importance in identifying an asteroid taxonomy. The new availability of geometric albedo data from the Wide-field Infrared Survey Explorer (WISE) mission for about 100,000 asteroids significantly increased the sample of objects for which such information, with some errors, is now known. In this work we proposed a new method to identify families halos in a multi-domain space composed by proper elements, SDSS-MOC4 (a*,i-z) colors, and WISE geometric albedo for the whole main belt (and the Hungaria and Cybele orbital regions). Assuming that most families were created by the breakup of an undifferentiated parent body, they are expected to be homogeneous in colors and albedo. The new method is quite effective in determining objects belonging to a family halo, with low percentages of likely interlopers, and results that are quite consistent in term of taxonomy and geometric albedo of the halo members.

The observed wide eccentricity distribution of extrasolar giant planets is thought to be the result of dynamical instabilities and gravitational scattering among planets. Previously, it has been assumed that the orbits in giant planet systems become gravitationally unstable after the gas nebula dispersal. It was not well understood, however, how these unstable conditions were established in the first place. In this work we numerically simulate the evolution of systems of three planets as the planets sequentially grow to Jupiter's mass, and dynamically interact among themselves and with the gas disk. We use the hydro-dynamical code FARGO that we modified by implementing the $N$-body integrator SyMBA. The new code can handle close encounters and collisions between planets. To test their stability, the planetary systems were followed with SyMBA for up to $10^8$ yr after the gas disk dispersal. We find that dynamics of the growing planets is complex, because migration and resonances raise their orbital eccentricities, and cause dynamical instabilities when gas is still around. If the dynamical instabilities occur early, planets can be removed by collisions and ejections, and the system rearranges into a new, more stable configuration. In this case, the planetary systems emerging from the gas disks are expected to be stable, and would need to be destabilized by other means (low-mass planets, planetesimal disks, etc.). Alternatively, for the giant planet system to be intrinsically unstable upon the gas disk dispersal, a special timing would be required with the growth of (at least some of) the giant planets having to occur near the end of the gas disk lifetime.

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB.

This paper explores the possibility that the progenitors of the small satellites of Pluto got captured in the Pluto-Charon system from the massive heliocentric planetesimal disk in which Pluto was originally embedded into. We find that, if the dynamical excitation of the disk is small, temporary capture in the Pluto-Charon system can occur with non-negligible probability, due to the dynamical perturbations exerted by the binary nature of the Pluto-Charon pair. However, the captured objects remain on very elliptic orbits and the typical capture time is only 100 years. In order to explain the origin of the small satellites of Pluto, we conjecture that some of these objects got disrupted during their Pluto-bound phase by a collision with a planetesimal of the disk. This could have generated a debris disk, which damped under internal collisional evolution, until turning itself into an accretional disk that could form small satellites on circular orbits, co-planar with Charon. Unfortunately, we find that objects large enough to carry a sufficient amount of mass to generate the small satellites of Pluto have collisional lifetimes orders of magnitude longer than the capture time. Thus, this scenario cannot explain the origin of the small satellites of Pluto, which remains elusive.

Aims: We present 11 high-precision photometric transit observations of the transiting super-Earth planet GJ1214b. Combining these data with observations from other authors, we investigate the ephemeris for possible signs of transit timing variations (TTVs) using a Bayesian approach. Methods: The observations were obtained using telescope-defocusing techniques, and achieve a high precision with random errors in the photometry as low as 1mmag per point. To investigate the possibility of TTVs in the light curve, we calculate the overall probability of a TTV signal using Bayesian methods. Results: The observations are used to determine the photometric parameters and the physical properties of the GJ1214 system. Our results are in good agreement with published values. Individual times of mid-transit are measured with uncertainties as low as 10s, allowing us to reduce the uncertainty in the orbital period by a factor of two. Conclusions: A Bayesian analysis reveals that it is highly improbable that the observed transit times is explained by TTV, when compared with the simpler alternative of a linear ephemeris.

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

This document covers five broad topics from theory of small bodies: planetesimal formation, cosmochemistry, thermal evolution, collisions and dynamics. Each of these topics is described in a separate section, where we prioritize selected main issues over completeness. The text points out the principal unresolved problems in each area, and suggests ways how progress can be made. This includes support for new code development, observations and experimental work that can be used to constrain theory, new research directions, and studies of cross-over regimes where the system's behavior is determined by several competing processes. The suggested development areas are placed in the context of NASA space exploration.

We have extended our earlier work on space weathering of the youngest S-complex asteroid families to include results from asteroid clusters with ages <10^6 years and to newly identified asteroid pairs with ages <5x10^5 years. We have identified three S-complex asteroid clusters with ages in the range 10^{5-6} years. The average color of the objects in these clusters agree with the prediction of Willman et al., 2008. SDSS photometry of the members of very young asteroid pairs with ages <10^5 years was used to determine their taxonomy. The average color of the S-complex pairs is PC_1=0.49+/-0.03, over 5-sigma redder than predicted by Willman et al., 2008. Therefore, the most likely pair formation mechanism is gentle separation due to YORP spin-up leaving much of the aged and reddened surface undisturbed. In this case our color measurement allows us to set an upper limit of ~64% on the disturbed surface portion. Using pre-existing color data and our new results for the youngest S-complex asteroid clusters we have extended our space weather model to explicitly include the effects of regolith gardening and fit separate weathering and gardening characteristic timescales of tau_w=960+/-160My and tau_g=2000+/-290My respectively. The first principal component color for fresh S-complex material is 0.37+/-0.01 while the maximum amount of local reddening is 0.33+/-0.06. Our first-ever determination of the gardening time is in stark contrast to our calculated gardening time of tau_g~270My based on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. A possible resolution for the discrepancy is through a `honeycomb' mechanism in which the surface regolith structure absorbs small impactors without producing significant ejecta. This mechanism could also account for the paucity of small craters on (433) Eros.

Rewriting a shorter abstract.

We report a transit timing study of the transiting exoplanetary system HD 189733. In total we observed ten transits in 2006 and 2008 with the 2.6-m Nordic Optical Telescope, and two transits in 2007 with the 4.2-m William Herschel Telescope. We used Markov-Chain Monte Carlo simulations to derive the system parameters and their uncertainties, and our results are in a good agreement with previously published values. We performed two independent analyses of transit timing residuals to place upper mass limits on putative perturbing planets. The results show no evidence for the presence of planets down to 1 Earth mass near the 1:2 and 2:1 resonance orbits, and planets down to 2.2 Earth masses near the 3:5 and 5:3 resonance orbits with HD 189733b. These are the strongest limits to date on the presence of other planets in this system.