Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.

The Jovian Trojans are two swarms of objects located around the L$_4$ and L$_5$ Lagrange points. The population is thought to have been captured by Jupiter during the Solar system's youth. Within the swarms, six collisional families have been identified in previous work, with four in the L$_4$ swarm, and two in the L$_5$. Our aim is to investigate the stability of the two Trojan swarms, with a particular focus on these collisional families. We find that the members of Trojan swarms escape the population at a linear rate, with the primordial L$_4$ (23.35% escape) and L$_5$ (24.89% escape) population sizes likely 1.31 and 1.35 times larger than today. Given that the escape rates were approximately equal between the two Trojan swarms, our results do not explain the observed asymmetry between the two groups, suggesting that the numerical differences are primordial in nature, supporting previous studies. Upon leaving the Trojan population, the escaped objects move onto orbits that resemble those of the Centaur and short-period comet populations. Within the Trojan collisional families, the 1996 RJ and 2001 UV$_{209}$ families are found to be dynamically stable over the lifetime of the Solar system, whilst the Hektor, Arkesilos and Ennomos families exhibit various degrees of instability. The larger Eurybates family shows 18.81\% of simulated members escaping the Trojan population. Unlike the L4 swarm, the escape rate from the Eurybates family is found to increase as a function of time, allowing an age estimation of approximately $1.045\pm 0.364 \times 10^9$ years.

For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.

A recent examination of K2 lightcurves indicates that ~15% of Jupiter Trojans have very slow rotation (spin periods Ps>100 h). Here we consider the possibility that these bodies formed as equal-size binaries in the massive outer disk at ~20-30 au. Prior to their implantation as Jupiter Trojans, tight binaries tidally evolved toward a synchronous state with Ps~Pb, where Pb is the binary orbit period. They may have been subsequently dissociated by impacts and planetary encounters with at least one binary component retaining its slow rotation. Surviving binaries on Trojan orbits would continue to evolve by tides and spin-changing impacts over 4.5 Gyr. To explain the observed fraction of slow rotators, we find that at least ~15-20% of outer disk bodies with diameters 15<D<50 km would have to form as equal-size binaries with 12<ab/R<30, where ab is the binary semimajor axis and R=D/2. The mechanism proposed here could also explain very slow rotators found in other small body populations.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

The discovery and characteristics of transneptunian binaries are reviewed. In the 20 years since their first discovery, a wealth of detail has emerged including the frequency of binaries in different populations, their relative sizes and separations, and colors. Taken globally, these properties give strong clues to the origin and evolution of the populations where these binaries are found. In the last ten years and increasing number of binary orbits have been determined which yields a new trove of information on their masses and densities as well as details of their orbits including inclination, eccentricity and the timing of mutual events. In 2018, the study of transneptunian binaries remains one of the most active areas of progress in understanding the Solar System beyond Neptune.

We use published models of the early Solar System evolution with a slow, long-range and grainy migration of Neptune to predict the orbital element distributions and the number of modern-day Centaurs. The model distributions are biased by the OSSOS survey simulator and compared with the OSSOS Centaur detections. We find an excellent match to the observed orbital distribution, including the wide range of orbital inclinations which was the most troublesome characteristic to fit in previous models. A dynamical model, in which the original population of outer disk planetesimals was calibrated from Jupiter Trojans, is used to predict that OSSOS should detect 11+/-4 Centaurs with semimajor axis a<30 au, perihelion distance q>7.5 au and diameter D>10 km (absolute magnitude H_r<13.7 for a 6% albedo). This is consistent with 15 actual OSSOS Centaur detections with H_r<13.7. The population of Centaurs is estimated to be 21,000+/-8,000 for D>10 km. The inner scattered disk at 50<a<200 au should contain (2.0+/-0.8)x10^7 D>10 km bodies and the Oort cloud should contain (5.0+/-1.9)x10^8 D>10 km comets. Population estimates for different diameter cutoffs can be obtained from the size distribution of Jupiter Trojans (N(>D) proportional to D^(-2.1) for 5<D<100 km). We discuss model predictions for the Large Synoptic Survey Telescope observations of Centaurs.

A critical step toward the emergence of planets in a protoplanetary disk consists in accretion of planetesimals, bodies 1-1000 km in size, from smaller disk constituents. This process is poorly understood partly because we lack good observational constraints on the complex physical processes that contribute to planetesimal formation. In the outer solar system, the best place to look for clues is the Kuiper belt, where icy planetesimals survived to this day. Here we report evidence that Kuiper belt planetesimals formed by the streaming instability, a process in which aerodynamically concentrated clumps of pebbles gravitationally collapse into 100-km-class bodies. Gravitational collapse was previously suggested to explain the ubiquity of equal-size binaries in the Kuiper belt. We analyze new hydrodynamical simulations of the streaming instability to determine the model expectations for the spatial orientation of binary orbits. The predicted broad inclination distribution with 80% of prograde binary orbits matches the observations of trans-Neptunian binaries. The formation models which imply predominantly retrograde binary orbits can be ruled out. Given its applicability over a broad range of protoplanetary disk conditions, it is expected that the streaming instability seeded planetesimal formation also elsewhere in the solar system, and beyond.

Here we describe a story behind the discovery of Kepler-46, which was the first exoplanetary system detected and characterized from a method known as the transit timing variations (TTVs). The TTV method relies on the gravitational interaction between planets orbiting the same star. If transits of at least one of the planets are detected, precise measurements of its transit times can be used, at least in principle, to detect and characterize other non-transiting planets in the system. Kepler-46 was the first case for which this method was shown to work in practice. Other detections and characterizations followed (e.g., Kepler-88). The TTV method plays an important role in addressing the incompleteness of planetary systems detected from transits.

As indicated by their special characteristics, the cold classical Kuiper belt objects (KBOs) formed and survived at 42-47 au. Notably, they show a large fraction of equal-size binaries whose formation is probably related to the accretion of KBOs themselves. These binaries are uncommon in other --hot, resonant, scattered-- populations, which are thought to have been implanted from the massive disk below 30 au to >30 au during Neptune's migration. Here we highlight the possibility that equal-size binaries formed in the disk but were subsequently removed by impacts and/or dynamical effects (e.g., scattering encounters with Neptune). We determine the dependence of these processes on the size and separation of binary components. Our results indicate that tighter binaries, if they formed in the massive disk, have relatively good chances of survival (unless the disk was long-lived). In contrast, the widest binaries in the hot population, such as 2002 VF130, have a very low survival probability (<1%) even if the massive disk was short-lived. They may represent a trace of lucky survivors of a much larger population of the original disk binaries, or they formed at ~30-40~au and dodged the impact- and encounter-related perturbations that we studied here. We find that all known satellites of the largest KBOs would survive during the dynamical implantation of these bodies in the Kuiper belt. The low orbital eccentricities of Pluto's small moons may have been excited by impacts and/or encounters of the Pluto system to Neptune.

This chapter reviews accretion models for Kuiper belt objects (KBOs), discussing in particular the compatibility of the observed properties of the KBO population with the streaming instability paradigm. Then it discusses how the dynamical structure of the KBO population, including the formation of its 5 sub-components (cold, hot, resonant, scattered and fossilized), can be quantitatively understood in the framewok of the giant planet instability. We also establish the connections between the KBO population and the Trojans of Jupiter and Neptune, the irregular satellites of all giant planets, the Oort cloud and the D-type main belt asteroids. Finally, we discuss the collisional evolution of the KBO population, arguing that the current size-frequency distribution below 100~km in size has been achieved as a collisional equilibrium in a few tens of My inside the original massive trans-Neptunain disk, possibly with the exception of the cold population sub-component.

We perform dynamical fits to the Transit Timing Variations (TTVs) of Kepler-419b. The TTVs from 17 Kepler quarters are obtained from Holczer et al (2016, ApJS 225,9). The dynamical fits are performed using the MultiNest Bayesian inference tool, coupled to an efficient symplectic N-body integrator. We find that the existing TTV data alone are able to uniquely constrain the planetary masses of Kepler-419b and c. Our estimates are in a good agreement with previous mass determinations that combined different techniques and observations, such as TTVs, radial velocity measurements, and the photoeccentric effect. As expected, however, our mass estimates have larger uncertainty. We study the global stability of the system within the parameters uncertainties to discard possible unstable solutions. We conclude that our method applied to the Holczer et al. data can provide reliable determinations of the planetary parameters. In the fore-coming work, we will use Holczer et al. (2016b) data to determine planetary parameters for a large set of Kepler systems.

The orbital distribution of trans-Neptunian objects provides strong evidence for the radial migration of Neptune. The outer planets' orbits are thought to have become unstable during the early stages with Jupiter having scattering encounters with a Neptune-class planet. As a consequence, Jupiter jumped inward by a fraction of an au, as required from inner solar system constraints, and obtained its current orbital eccentricity. The timing of these events is often linked to the lunar Late Heavy Bombardment that ended ~700 Myr after the dispersal of the protosolar nebula (t_0). Here we show instead that planetary migration started shortly after t_0. Such early migration is inferred from the survival of the Patroclus-Menoetius binary Jupiter Trojan. The binary formed at t <~ t_0 within a massive planetesimal disk once located beyond Neptune. The longer the binary stayed in the disk, the greater the likelihood that collisions would strip its components from one another. The simulations of its survival indicate that the disk had to have been dispersed by migrating planets within <100 Myr of t_0. This constraint implies that the planetary migration is unrelated to the formation of the youngest lunar basins.

The main asteroid belt (MB) is low in mass but dynamically excited. Here we propose a new mechanism to excite the MB during the giant planet ('Nice model') instability, which is expected to have featured repeated close encounters between Jupiter and one or more ice giants ('Jumping Jupiter' -- JJ). We show that, when Jupiter temporarily reaches a high enough level of excitation, both in eccentricity and inclination it induces strong forced vectors of eccentricity and inclination across the MB region. Because during the JJ instability Jupiter's orbit `jumps' around, the forced vectors keep changing both in magnitude and phase throughout the whole MB region. The entire cold primordial MB is thus excited as a natural outcome of the JJ instability. The level of such an excitation, however, is typically larger than the current orbital excitation observed in the MB. We show that the subsequent evolution of the Solar System is capable of reshaping the resultant over-excited MB to its present day orbital state, and that a strong mass depletion ($\sim$90$\%$) is associated to the JJ instability phase and its subsequent evolution throughout the age of the Solar System.

Several properties of the Solar System, including the wide radial spacing of the giant planets, can be explained if planets radially migrated by exchanging orbital energy and momentum with outer disk planetesimals. Neptune's planetesimal-driven migration, in particular, has a strong advocate in the dynamical structure of the Kuiper belt. A dynamical instability is thought to have occurred during the early stages with Jupiter having close encounters with a Neptune-class planet. As a result of the encounters, Jupiter acquired its current orbital eccentricity and jumped inward by a fraction of an au, as required for the survival of the terrestrial planets and from asteroid belt constraints. Planetary encounters also contributed to capture of Jupiter Trojans and irregular satellites of the giant planets. Here we discuss the dynamical evolution of the early Solar System with an eye to determining how models of planetary migration/instability can be constrained from its present architecture.

The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $\nu_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $\nu_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component).

The Rosetta spacecraft observations revealed that the nucleus of comet 67P/Churyumov-Gerasimenko consists of two similarly sized lobes connected by a narrow neck. Here we evaluate the possibility that 67P is a collapsed binary. We assume that the progenitor of 67P was a binary and consider various physical mechanisms that could have brought the binary components together, including small-scale impacts and gravitational encounters with planets. We find that 67P could be a primordial body (i.e., not a collisional fragment) if the outer planetesimal disk lasted <10 Myr before it was dispersed by migrating Neptune. The probability of binary collapse by impact is 30% for tightly bound binaries. Most km-class binaries become collisionally dissolved. Roughly 10% of the surviving binaries later evolve to become contact binaries during the disk dispersal, when bodies suffer gravitational encounters to Neptune. Overall, the processes described in this work do not seem to be efficient enough to explain the large fraction (~67%) of bi-lobed cometary nuclei inferred from spacecraft imaging.

The timeline of the lunar bombardment in the first Gy of the Solar System remains unclear. Some basin-forming impacts occurred 3.9-3.7Gy ago. Many other basins formed before, but their exact ages are not precisely known. There are two possible interpretations of the data: in the cataclysm scenario there was a surge in the impact rate approximately 3.9Gy ago, while in the accretion tail scenario the lunar bombardment declined since the era of planet formation and the latest basins formed in its tail-end. Here, we revisit the work of Morbidelli et al.(2012) that examined which scenario could be compatible with both the lunar crater record in the 3-4Gy period and the abundance of highly siderophile elements (HSE) in the lunar mantle. We use updated numerical simulations of the fluxes of impactors. Under the traditional assumption that the HSEs track the total amount of material accreted by the Moon since its formation, we conclude that only the cataclysm scenario can explain the data. The cataclysm should have started ~3.95Gy ago. However we show that HSEs could have been sequestered from the lunar mantle due to iron sulfide exsolution during magma ocean crystallization, followed by mantle overturn. Based on the hypothesis that the lunar magma ocean crystallized about 100-150My after Moon formation, and therefore that HSEs accumulated in the lunar mantle only after this time, we show that the bombardment in the 3-4Gy period can be explained in the accretion tail scenario. This hypothesis would also explain why the Moon appears so depleted in HSEs relative to the Earth. We also extend our analysis of the cataclysm and accretion tail scenarios to the case of Mars. The accretion tail scenario requires a global resurfacing event on Mars ~4.4Gy ago, possibly associated with the formation of the Borealis basin, and it is consistent with the HSE budget of the planet.

Dynamical models of the asteroid delivery from the main belt suggest that the current impact flux of diameter D>10 km asteroids on the Earth is 0.5-1 per Gyr. Studies of the Near-Earth Asteroid (NEA) population find a much higher flux, with ~7 D>10-km asteroid impacts per Gyr. Here we show that this problem is rooted in the application of impact probability of small NEAs (1.5 per Gyr per object), whose population is well characterized, to large NEAs. In reality, large NEAs evolve from the main belt by different escape routes, have a different orbital distribution, and lower impact probabilities (0.8+/-0.3 per Gyr per object) than small NEAs. In addition, we find that the current population of two D>10 km NEAs (Ganymed and Eros) is a slight fluctuation over the long term average of 1.1+/-0.5 D>10 km NEAs in a steady state. These results have important implications for our understanding of the occurrence of the K/T-scale impacts on the terrestrial worlds.

The origin of the orbital structure of the cold component of the Kuiper belt is still a hot subject of investigation. Several features of the solar system suggest that the giant planets underwent a phase of global dynamical instability, but the actual dynamical evolution of the planets during the instability is still debated. To explain the structure of the cold Kuiper belt, Nesvorny (2015, AJ 150,68) argued for a soft instability, during which Neptune never achieved a very eccentric orbit. Here we investigate the possibility of a more violent instability, from an initially more compact fully resonant configuration of 5 giant planets. We show that the orbital structure of the cold Kuiper belt can be reproduced quite well provided that the cold population formed in situ, with an outer edge between 44 and 45 au and never had a large mass.

Comets are icy objects that orbitally evolve from the trans-Neptunian region (the Kuiper belt and beyond) into the inner Solar System, where they are heated by solar radiation and become active due to sublimation of water ice. Here we perform end-to-end simulations in which cometary reservoirs are formed in the early Solar System and evolved over 4.5 Gyr. The gravitational effects of Planet 9 (P9), hypothesized to circle the Sun on a wide orbit, are included in some of our simulations. Different models are considered for comets to be active, including a simple assumption that comets remain active for Np(q) perihelion passages with perihelion distance q<2.5 au. The orbital distribution and number of active comets produced in our model is compared to observations. The orbital distribution of ecliptic comets (ECs) is well reproduced in models with Np(2.5)=500 and without P9. With P9, the inclination distribution of model ECs is wider than the observed one. We find that the known Halley-type comets (HTCs) have a nearly isotropic inclination distribution (with only a slight preference for prograde orbits). In our model, the HTCs appear to be an extension of the population of returning Oort-cloud comets (OCCs) to shorter orbital periods. The inclination distribution of model HTCs becomes broader with increasing Np, but the existing observational data are not good enough to constrain Np from orbital fits. Np(2.5)>1000 is required to obtain a steady-state population of large active HTCs that is consistent with observations. To fit the ratio of the returning-to-new OCCs, by contrast, our model implies that Np(2.5)<10, possibly because the detected long-period comets are smaller and much easier to disrupt than observed HTCs.

Recent works on planetary migration show that the orbital structure of the Kuiper belt can be very well reproduced if before the onset of the planetary instability Neptune underwent a long-range planetesimal-driven migration up to $\sim$28 au. However, considering that all giant planets should have been captured in mean motion resonances among themselves during the gas-disk phase, it is not clear whether such a very specific evolution for Neptune is possible, nor whether the instability could have happened at late times. Here, we first investigate which initial resonant configuration of the giant planets can be compatible with Neptune being extracted from the resonant chain and migrating to $\sim$28 au before that the planetary instability happened. We address the late instability issue by investigating the conditions where the planets can stay in resonance for about 400 My. Our results indicate that this can happen only in the case where the planetesimal disk is beyond a specific minimum distance $\delta_{stab}$ from Neptune. Then, if there is a sufficient amount of dust produced in the planetesimal disk, that drifts inwards, Neptune can enter in a slow dust-driven migration phase for hundreds of Mys until it reaches a critical distance $\delta_{mig}$ from the disk. From that point, faster planetesimal-driven migration takes over and Neptune continues migrating outward until the instability happens. We conclude that, although an early instability reproduces more easily the evolution of Neptune required to explain the structure of the Kuiper belt, such evolution is also compatible with a late instability.

The Asteroid Belt is characterized by the radial mixing of bodies with different physical properties, a very low mass compared to Minimum Mass Solar Nebula expectations and has an excited orbital distribution. Models of the evolution of the Asteroid Belt show that the origin of its structure is strongly linked to the process of terrestrial planet formation. The Grand Tack model presents a possible solution to the conundrum of reconciling the small mass of Mars with the properties of the Asteroid Belt, including the mass depletion, radial mixing and orbital excitation. However, while the inclination distribution produced in the Grand Tack model is in good agreement with the one observed, the eccentricity distribution is skewed towards values larger than those found today. Here, we evaluate the evolution of the orbital properties of the Asteroid Belt from the end of the Grand Tack model (at the end of the gas nebula phase when planets emerge from the dispersing gas disk), throughout the subsequent evolution of the Solar System including an instability of the Giant Planets approximately 400 My later. Before the instability, the terrestrial planets were modeled on dynamically cold orbits with Jupiter and Saturn locked in a 3:2 mean motion resonance. The model continues for an additional 4.1 Gy after the giant planet instability. Our results show that the eccentricity distribution obtained in the Grand Tack model evolves towards one very similar to that currently observed, and the semimajor axis distribution does the same. The inclination distribution remains nearly unchanged with a slight preference for depletion at low inclination; this leads to the conclusion that the inclination distribution at the end of the Grand Tack is a bit over-excited.

The impact cratering record of the Moon and the terrestrial planets provides important clues about the formation and evolution of the Solar System. Especially intriguing is the epoch 3.8-3.9 Gyr ago (Ga), known as the Late Heavy Bombardment (LHB), when the youngest lunar basins such as Imbrium and Orientale formed. The LHB was suggested to originate from a slowly declining impactor flux or from a late dynamical instability. Here we develop a model for the historical flux of large asteroid impacts and discuss how it depends on various parameters, including the time and nature of the planetary migration/instability. We find that the asteroid impact flux dropped by 1 to 2 orders of magnitude during the first 1 Gyr and remained relatively unchanged over the last 3 Gyr. The early impacts were produced by asteroids whose orbits became excited during the planetary migration/instability, and by those originating from the inner extension of the main belt (E-belt; semimajor axis 1.6<a<2.1 au). The profiles obtained for the early and late versions of the planetary instability initially differ, but end up being similar after ~3 Ga. Thus, the time of the instability can only be determined by considering the cratering and other constraints during the first ~1.5 Gyr of the Solar System history. Our absolute calibration of the impact flux indicates that asteroids were probably not responsible for the LHB, independently of whether the instability happened early or late, because the calibrated flux is not large enough to explain Imbrium/Orientale and a significant share of large lunar craters. Comets and leftovers of the terrestrial planet formation provided additional, and probably dominant source of impacts during early epochs.

The orbital structure of the asteroid belt holds a record of the Solar System's dynamical history. The current belt only contains ${\rm \sim 10^{-3}}$ Earth masses yet the asteroids' orbits are dynamically excited, with a large spread in eccentricity and inclination. In the context of models of terrestrial planet formation, the belt may have been excited by Jupiter's orbital migration. The terrestrial planets can also be reproduced without invoking a migrating Jupiter; however, as it requires a severe mass deficit beyond Earth's orbit, this model systematically under-excites the asteroid belt. Here we show that the orbits of the asteroids may have been excited to their current state if Jupiter and Saturn's early orbits were chaotic. Stochastic variations in the gas giants' orbits cause resonances to continually jump across the main belt and excite the asteroids' orbits on a timescale of tens of millions of years. While hydrodynamical simulations show that the gas giants were likely in mean motion resonance at the end of the gaseous disk phase, small perturbations could have driven them into a chaotic but stable state. The gas giants' current orbits were achieved later, during an instability in the outer Solar System. Although it is well known that the present-day Solar System exhibits chaotic behavior, our results suggest that the early Solar System may also have been chaotic.

The dynamical structure of the Kuiper belt beyond 50 au is not well understood. Here we report results of a numerical model with long-range, slow and grainy migration of Neptune. The model implies that bodies scattered outward by Neptune to semimajor axes a>50 au often evolve into resonances which subsequently act to raise the perihelion distances of orbits to q>40 au. The implication of the model is that the orbits with 50<a<100 au and q>40 au should cluster near (but not in) the resonances with Neptune (3:1 at a=62.6 au, 4:1 at a=75.9 au, 5:1 at a=88.0 au, etc.). The recent detection of several distant Kuiper Belt Objects (KBOs) near resonances is consistent with this prediction, but it is not yet clear whether the orbits are really non-resonant as our model predicts. We estimate from the model that there should presently be ~1600-2400 bodies at the 3:1 resonance and ~1000-1400 bodies at the 4:1 resonance (for q>40 au and diameters D>100 km). These results favorably compare with the population census of distant KBOs inferred from existing observations.

The Transit Timing Variations (TTVs) are deviations of the measured mid-transit times from the exact periodicity. One of the most interesting causes of TTVs is the gravitational interaction between planets. Here we consider a case of two planets in a mean motion resonance (orbital periods in a ratio of small integers). This case is important because the resonant interaction can amplify the TTV effect and allow planets to be detected more easily. We develop an analytic model of the resonant dynamics valid for small orbital eccentricities and use it to derive the principal TTV terms. We find that a resonant system should show TTV terms with two basic periods (and their harmonics). The resonant TTV period is proportional (m/M_*)^(-2/3), where m and M_* are the planetary and stellar masses. For m=10^(-4) M_*, for example, the TTV period exceeds the orbital period by ~2 orders of magnitude. The amplitude of the resonant TTV terms scales linearly with the libration amplitude. The ratio of the TTV amplitudes of two resonant planets is inversely proportional to the ratio of their masses. These and other relationships discussed in the main text can be used to aid the interpretation of TTV observations.

The Kuiper belt is a population of icy bodies beyond the orbit of Neptune. The complex orbital structure of the Kuiper belt, including several categories of objects inside and outside of resonances with Neptune, emerged as a result of Neptune's migration into an outer planetesimal disk. An outstanding problem with the existing migration models is that they invariably predict excessively large resonant populations, while observations show that the non-resonant orbits are in fact common (e.g., the main belt population is 2-4 times larger than Plutinos in the 3:2 resonance). Here we show that this problem can resolved if it is assumed that Neptune's migration was grainy, as expected from scattering encounters of Neptune with massive planetesimals. The grainy migration acts to destabilize resonant bodies with large libration amplitudes, a fraction of which ends up on stable non-resonant orbits. Thus, the non-resonant--to--resonant ratio obtained with the grainy migration is higher, up to ~10 times higher for the range of parameters investigated here, than in a model with smooth migration. The grainy migration leads to a narrower distribution of the libration amplitudes in the 3:2 resonance. The best fit to observations is obtained when it is assumed that the outer planetesimal disk below 30 AU contained 1000-4000 Plutos, or ~1000 bodies twice as massive as Pluto. We find that the probability for an outer disk object to end up on a stable orbit in the Kuiper belt is ~10^{-3}. Together, these results are consistent with having only two (known) Pluto-mass bodies in the Kuiper belt today (Pluto and Eris). The combined mass of Pluto-class objects in the original disk was ~2-8 Earth masses (M_Earth), which represents 10-40% of the estimated disk mass (M_disk ~ 20 M_Earth).

Themis is an old and statistically robust asteroid family populating the outer main belt, and resulting from a catastrophic collision that took place 2.5$\pm$1.0 Gyr ago. Within the old Themis family a young sub-family, Beagle, formed less than 10 Myr ago, has been identified. We present the results of a spectroscopic survey in the visible and near infrared range of 22 Themis and 8 Beagle families members. The Themis members investigated exhibit a wide range of spectral behaviors, while the younger Beagle family members look spectrally bluer with a smaller spectral slope variability. The best meteorite spectral analogues found for both Themis and Beagle families members are carbonaceous chondrites having experienced different degrees of aqueous alteration, prevalently CM2 but also CV3 and CI, and some of them are chondrite samples being unusual or heated. We extended the spectral analysis including the data available in the literature on Themis and Beagle families members, and we looked for correlations between spectral behavior and physical parameters using the albedo and size values derived from the WISE data. The analysis of this larger sample confirm the spectral diversity within the Themis family and that Beagle members tend to be bluer and to have an higher albedo. The differences between the two family may be partially explained by space weathering processes, which act on these primitive surfaces in a similar way than on S-type asteroids, i.e. producing reddening and darkening. However we see several Themis members having albedos and spectral slopes similar to the young Beagle members. Alternative scenarios are proposed including heterogeneity in the parent body having a compositional gradient with depth, and/or the survival of projectile fragments having a different composition than the parent body.

We present a survey of 41 Kepler Objects of Interest (KOIs) for exomoons using Bayesian photodynamics, more than tripling the number of KOIs surveyed with this technique. We find no compelling evidence for exomoons although thirteen KOIs yield spurious detections driven by instrumental artifacts, stellar activity and/or perturbations from unseen bodies. Regarding the latter, we find seven KOIs exhibiting >5 sigma evidence of transit timing variations, including the 'mega-Earth' Kepler-10c, likely indicating an additional planet in that system. We exploit the moderately large sample of 57 unique KOIs surveyed to date to infer several useful statistics. For example, although there is a diverse range in sensitivities, we find that we are sensitive to Pluto-Charon mass-ratio systems for ~40% of KOIs studied and Earth-Moon mass-ratios for 1 in 8 cases. In terms of absolute mass, our limits probe down to 1.7 Ganymede masses, with a sensitivity to Earth-mass moons for 1 in 3 cases studied and to the smallest moons capable of sustaining an Earth-like atmosphere (0.3 Earth masses) for 1 in 4. Despite the lack of positive detections to date, we caution against drawing conclusions yet, since our most interesting objects remain under analysis. Finally, we point out that had we searched for the photometric transit signals of exomoons alone, rather than using photodynamics, we estimate that 1 in 4 KOIs would have erroneously been concluded to harbor exomoons due to residual time correlated noise in the Kepler data, posing a serious problem for alternative methods.

The solar system has changed dramatically since its birth, and so did our understanding of it. A considerable research effort has been invested in the past decade in an attempt to reconstruct the solar system history, including the earliest stages some 4.5 billion years ago. The results indicate how several processes, such as planetary migration and dynamical instabilities, acted to relax the orbital spacing of the outer planets, and provided the needed perturbation to explain the present planetary orbits that are not precisely circular and coplanar. Here we highlight this work and illustrate the key results in a computer simulation that unifies several recently developed theories. The emerging view represents another step away from the initial perception of the solar system as part of unchanging heavens.

The Kuiper belt is a population of icy bodies beyond the orbit of Neptune. A particularly puzzling and up-to-now unexplained feature of the Kuiper belt is the so-called `kernel', a concentration of orbits with semimajor axes a~44 AU, eccentricities e~0.05, and inclinations i<5 deg. Here we show that the Kuiper belt kernel can be explained if Neptune's otherwise smooth migration was interrupted by a discontinuous change of Neptune's semimajor axis when Neptune reached ~28 AU. Before the discontinuity happened, planetesimals located at ~40 AU were swept into Neptune's 2:1 resonance, and were carried with the migrating resonance outwards. The 2:1 resonance was at ~44 AU when Neptune reached ~28 AU. If Neptune's semimajor axis changed by fraction of AU at this point, perhaps because Neptune was scattered off of another planet, the 2:1 population would have been released at ~44 AU, and would remain there to this day. We show that the orbital distribution of bodies produced in this model provides a good match to the orbital properties of the kernel. If Neptune migration was conveniently slow after the jump, the sweeping 2:1 resonance would deplete the population of bodies at ~45-47 AU, thus contributing to the paucity of the low-inclination orbits in this region. Special provisions, probably related to inefficiencies in the accretional growth of sizable objects, are still needed to explain why only a few low-inclination bodies have been so far detected beyond ~47 AU.

We study the effects of planetary late migration on the gas giants obliquities. We consider the planetary instability models from Nesvorny & Morbidelli (2012), in which the obliquities of Jupiter and Saturn can be excited when the spin-orbit resonances occur. The most notable resonances occur when the $s_7$ and $s_8$ frequencies, changing as a result of planetary migration, become commensurate with the precession frequencies of Jupiter's and Saturn's spin vectors. We show that Jupiter may have obtained its present obliquity by crossing of the $s_8$ resonance. This would set strict constrains on the character of migration during the early stage. Additional effects on Jupiter's obliquity are expected during the last gasp of migration when the $s_7$ resonance was approached. The magnitude of these effects depends on the precise value of the Jupiter's precession constant. Saturn's large obliquity was likely excited by capture into the $s_8$ resonance. This probably happened during the late stage of planetary migration when the evolution of the $s_8$ frequency was very slow, and the conditions for capture into the spin-orbit resonance with $s_8$ were satisfied. However, whether or not Saturn is in the spin-orbit resonance with $s_8$ at the present time is not clear, because the existing observations of Saturn's spin precession and internal structure models have significant uncertainties.

Much of the dynamical structure of the Kuiper belt can be explained if Neptune migrated over several AU, and/or if Neptune was scattered to an eccentric orbit during planetary instability. An outstanding problem with the existing formation models is that the distribution of orbital inclinations they predicted is narrower than the one inferred from observations. Here we perform numerical simulations of Kuiper belt formation starting from an initial state with Neptune at 20 < a_{N,0} < 30 AU and a dynamically cold outer disk extending from beyond a_{N,0} to 30 AU. Neptune's orbit is migrated into the disk on an e-folding timescale 1 <= tau <= 100 Myr. A small fraction (~10^{-3}) of the disk planetesimals become implanted into the Kuiper belt in the simulations. By analyzing the orbital distribution of the implanted bodies in different cases we find that the inclination constraint implies that tau >= 10 Myr and a_{N,0} <= 25 AU. The models with tau < 10 Myr do not satisfy the inclination constraint, because there is not enough time for various dynamical processes to raise inclinations. The slow migration of Neptune is consistent with other Kuiper belt constraints, and with recently developed models of planetary instability/migration. Neptune's eccentricity and inclination are never large in these models (e_N<0.1, i_N<2 deg), as required to avoid excessive orbital excitation in the >40 AU region, where the Cold Classicals presumably formed.

Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large asteroids have been disrupted by impacts over the age of the Solar System, producing groups of fragments known as asteroid families. Here we explain how the asteroid families are identified, review their current inventory, and discuss how they can be used to get insights into long-term dynamics of main belt asteroids. Electronic tables of the membership for 122 notable families are reported on the Planetary Data System node.

The discovery of an exomoon would provide deep insights into planet formation and the habitability of planetary systems, with transiting examples being particularly sought after. Of the hundreds of Kepler planets now discovered, the seven-planet system Kepler-90 is unusual for exhibiting an unidentified transit-like signal in close proximity to one of the transits of the long-period gas-giant Kepler-90g, as noted by Cabrera et. al. (2014). As part of the 'Hunt for Exomoons with Kepler' (HEK) project, we investigate this possible exomoon signal and find it passes all conventional photometric, dynamical and centroid diagnostic tests. However, pixel-level light curves indicate that the moon-like signal occurs on nearly all of the target's pixels, which we confirm using a novel way of examining pixel-level data which we dub the 'transit centroid'. This test reveals that the possible exomoon to Kepler-90g is likely a false positive, perhaps due to a cosmic ray induced Sudden Pixel Sensitivity Dropout (SPSD). This work highlights the extreme care required for seeking non-periodic low-amplitude transit signals, such as exomoons.

Embedded in the gaseous protoplanetary disk, Jupiter and Saturn naturally become trapped in 3:2 resonance and migrate outward. This serves as the basis of the Grand Tack model. However, previous hydrodynamical simulations were restricted to isothermal disks, with moderate aspect ratio and viscosity. Here we simulate the orbital evolution of the gas giants in disks with viscous heating and radiative cooling. We find that Jupiter and Saturn migrate outward in 3:2 resonance in modest-mass ($M_{disk} \approx M_{MMSN}$, where MMSN is the "minimum-mass solar nebula") disks with viscous stress parameter $\alpha$ between $10^{-3}$ and $10^{-2} $. In disks with relatively low-mass ($M_{disk} \lesssim M_{MMSN}$) , Jupiter and Saturn get captured in 2:1 resonance and can even migrate outward in low-viscosity disks ($\alpha \le 10^{-4}$). Such disks have a very small aspect ratio ($h\sim 0.02-0.03$) that favors outward migration after capture in 2:1 resonance, as confirmed by isothermal runs which resulted in a similar outcome for $h \sim 0.02$ and $\alpha \le 10^{-4}$. We also performed N-body runs of the outer Solar System starting from the results of our hydrodynamical simulations and including 2-3 ice giants. After dispersal of the gaseous disk, a Nice model instability starting with Jupiter and Saturn in 2:1 resonance results in good Solar Systems analogs. We conclude that in a cold Solar Nebula, the 2:1 resonance between Jupiter and Saturn can lead to outward migration of the system, and this may represent an alternative scenario for the evolution of the Solar System.

Saturn's moon Iapetus has an orbit in a transition region where the Laplace surface is bending from the equator to the orbital plane of Saturn. The orbital inclination of Iapetus to the local Laplace plane is ~8 deg, which is unexpected, because the inclination should be ~0 if Iapetus formed from a circumplanetary disk on the Laplace surface. It thus appears that some process has pumped up Iapetus's inclination while leaving its eccentricity near zero (e=0.03 at present). Here we examined the possibility that Iapetus's inclination was excited during the early solar system instability when encounters between Saturn and ice giants occurred. We found that the dynamical effects of planetary encounters on Iapetus's orbit sensitively depend on the distance of the few closest encounters. In four out of ten instability cases studied here, the orbital perturbations were too large to be plausible. In one case, Iapetus's orbit was practically unneffected. In the remaining five cases, the perturbations of Iapetus's inclination were adequate to explain its present value. In three of these cases, however, Iapetus's eccentricity was excited to 0.1-0.25, and it is not clear whether it could have been damped to its present value (0.03) by some subsequent process (e.g., tides and dynamical friction from captured irregular satellites do not seem to be strong enough). Our results therefore imply that only 2 out of 10 instability cases (~20%) can excite Iapetus's inclination to its present value (~30% of trials lead to >5 deg) while leaving its orbital eccentricity low.

KOI-227, KOI-319 and KOI-884 are identified here as (at least) two planet systems. For KOI-319 and KOI-884, the observed Transit Timing Variations (TTVs) of the inner transiting planet are used to detect an outer non-transiting planet. The outer planet in KOI-884 is 2.6 Jupiter masses and has the orbital period just narrow of the 3:1 resonance with the inner planet (orbital period ratio 2.93). The distribution of parameters inferred from KOI-319.01's TTVs is bimodal with either a 1.6 Neptune-mass (M_N) planet wide of the 5:3 resonance (period 80.1 d) or a 1 Saturn-mass planet wide of the 7:3 resonance (period 109.2 d). The radial velocity measurements can be used in this case to determine which of these parameter modes is correct. KOI-227.01's TTVs with large 10 hour amplitude can be obtained for planetary-mass companions in various major resonances. Based on the Bayesian evidence, the current TTV data favor the outer 2:1 resonance with a companion mass 1.5 M_N, but this solution implies a very large density of KOI-227.01. The inner and outer 3:2 resonance solutions with sub-Neptune-mass companions are physically more plausible, but will need to be verified.

While not detected yet, pairs of exoplanets in the 1:1 mean motion resonance probably exist. Low eccentricity, near-planar orbits, which in the comoving frame follow the horseshoe trajectories, are one of the possible stable configurations. Here we study transit timing variations produced by mutual gravitational interaction of planets in this orbital architecture, with the goal to develop methods that can be used to recognize this case in observational data. In particular, we use a semi-analytic model to derive parametric constraints that should facilitate data analysis. We show that characteristic traits of the transit timing variations can directly constrain the (i) ratio of planetary masses, and (ii) their total mass (divided by that of the central star) as a function of the minimum angular separation as seen from the star. In an ideal case, when transits of both planets are observed and well characterized, the minimum angular separation can also be inferred from the data. As a result, parameters derived from the observed transit timing series alone can directly provide both planetary masses scaled to the central star mass.

We develop an analytic model for transit timing variations produced by orbital conjunctions between gravitationally interacting planets. If the planetary orbits have tight orbital spacing, which is a common case among the Kepler planets, the effect of a single conjunction can be best described as: (1) a step-like change of the transit timing ephemeris with subsequent transits of the inner planet being delayed and those of the outer planet being sped up, and (2) a discrete change in sampling of the underlying oscillations from eccentricity-related interaction terms. In the limit of small orbital eccentricities, our analytic model gives explicit equations for these effects as a function of the mass and orbital separation of planets. We point out that a detection of the conjunction effect in real data is of crucial importance for the physical characterization of planetary systems from transit timing variations.

The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present Solar System, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites, and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorny & Morbidelli (2012) that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorny & Morbidelli (2012) clearly fails this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early Solar System evolved.

Transit timing variations (TTVs) have proven to be a powerful technique for confirming Kepler planet candidates, for detecting non-transiting planets, and for constraining the masses and orbital elements of multi-planet systems. These TTV applications often require the numerical integration of orbits for computation of transit times (as well as impact parameters and durations); frequently tens of millions to billions of simulations are required when running statistical analyses of the planetary system properties. We have created a fast code for transit timing computation, TTVFast, which uses a symplectic integrator with a Keplerian interpolator for the calculation of transit times (Nesvorny et al. 2013). The speed comes at the expense of accuracy in the calculated times, but the accuracy lost is largely unnecessary, as transit times do not need to be calculated to accuracies significantly smaller than the measurement uncertainties on the times. The time step can be tuned to give sufficient precision for any particular system. We find a speed-up of at least an order of magnitude relative to dynamical integrations with high precision using a Bulirsch-Stoer integrator.

With their smaller radii and high cosmic abundance, transiting planets around cool stars hold a unique appeal. As part of our on-going project to measure the occurrence rate of extrasolar moons, we here present results from a survey focussing on eight Kepler planetary candidates associated with M-dwarfs. Using photodynamical modeling and Bayesian multimodal nested sampling, we find no compelling evidence for an exomoon in these eight systems. Upper limits on the presence of such bodies probe down to $\sim0.4M_{\oplus}$ in the best case. For KOI-314, we are able to confirm the planetary nature of two out of the three known transiting candidates using transit timing variations. Of particular interest is KOI-314c, which is found to have a mass of $1.0_{-0.3}^{+0.4}M_{\oplus}$, making it the lowest mass transiting planet discovered to date. With a radius of $1.61_{-0.15}^{+0.16}R_{\oplus}$, this Earth-mass world is likely enveloped by a significant gaseous envelope comprising $\geq17_{-13}^{+12}$% of the planet by radius. We find evidence to support the planetary nature of KOI-784 too via transit timing, but we advocate further observations to verify the signals. In both systems, we infer that the inner planet has a higher density than the outer world, which may be indicative of photo-evaporation. These results highlight both the ability of Kepler to search for sub-Earth mass moons and the exciting ancillary science which often results from such efforts.

The irregular satellites of outer planets are thought to have been captured from heliocentric orbits. The exact nature of the capture process, however, remains uncertain. We examine the possibility that irregular satellites were captured from the planetesimal disk during the early Solar System instability when encounters between the outer planets occurred (Nesvorny, Vokrouhlicky & Morbidelli 2007, AJ 133; hereafter NVM07). NVM07 already showed that the irregular satellites of Saturn, Uranus and Neptune were plausibly captured during planetary encounters. Here we find that the current instability models present favorable conditions for capture of irregular satellites at Jupiter as well, mainly because Jupiter undergoes a phase of close encounters with an ice giant. We show that the orbital distribution of bodies captured during planetary encounters provides a good match to the observed distribution of irregular satellites at Jupiter. The capture efficiency for each particle in the original transplanetary disk is found to be (1.3-3.6)x10^-8. This is roughly enough to explain the observed population of jovian irregular moons. We also confirm NVM07's results for the irregular satellites of Saturn, Uranus and Neptune.

Dawson and Murray-Clay (2012) pointed out that the inner part of the cold population in the Kuiper belt (that with semi major axis a<43.5 AU) has orbital eccentricities significantly smaller than the limit imposed by stability constraints. Here, we confirm their result by looking at the orbital distribution and stability properties in proper element space. We show that the observed distribution could have been produced by the slow sweeping of the 4/7 mean motion resonance with Neptune that accompanied the end of Neptune's migration process. The orbital distribution of the hot Kuiper belt is not significantly affected in this process, for the reasons discussed in the main text. Therefore, the peculiar eccentricity distribution of the inner cold population can not be unequivocally interpreted as evidence that the cold population formed in-situ and was only moderately excited in eccentricity; it can simply be the signature of Neptune's radial motion, starting from a moderately eccentric orbit. We discuss how this agrees with a scenario of giant planet evolution following a dynamical instability and, possibly, with the radial transport of the cold population.

Kepler-22b is the first transiting planet to have been detected in the habitable-zone of its host star. At 2.4 Earth radii, Kepler-22b is too large to be considered an Earth-analog, but should the planet host a moon large enough to maintain an atmosphere, then the Kepler-22 system may yet possess a telluric world. Aside from being within the habitable-zone, the target is attractive due to the availability of previously measured precise radial velocities and low intrinsic photometric noise, which has also enabled asteroseismology studies of the star. For these reasons, Kepler-22b was selected as a target-of-opportunity by the 'Hunt for Exomoons with Kepler' (HEK) project. In this work, we conduct a photodynamical search for an exomoon around Kepler-22b leveraging the transits, radial velocities and asteroseismology plus several new tools developed by the HEK project to improve exomoon searches. We find no evidence for an exomoon around the planet and exclude moons of mass >0.5 Earth masses to 95% confidence. By signal injection and blind retrieval, we demonstrate that an Earth-like moon is easily detected for this planet even when the time-correlated noise of the data set is taken into account. We provide updated parameters for the planet Kepler-22b including a revised mass of <53 Earth masses to 95% confidence and an eccentricity of 0.13(-0.13)(+0.36) by exploiting Single-body Asterodensity Profiling (SAP). Finally, we show that Kepler-22b has a >95% probability of being within the empirical habitable-zone but a <5% probability of being within the conservative habitable-zone.

Evidence in the Solar System suggests that the giant planets underwent an epoch of radial migration that was very rapid, with an e-folding timescale shorter than 1~Myr. It is probable that the cause of this migration was that the giant planets experienced an orbital instability that caused them to encounter each other, resulting in radial migration. Several works suggest that this dynamical instability occurred `late', long after all the planets had formed and the solar nebula had dissipated. Assuming that the terrestrial planets had already formed, then their orbits would have been affected by the migration of the giant planets. As a result, how did the orbits of the terrestrial planets change? And can we use this migration to obtain information on the primordial orbits of the terrestrial planets? We directly model a large number of terrestrial planet systems and their response to giant planet migration. We study the change in the Angular Momentum Deficit (AMD) of the terrstrials. We conclude that the primordial AMD should have been lower than ~70\% of the current value, but higher than 10\%. We find that a scenario with five giant planets better satisfies the orbital constraints of the terrestrial planets. We predict that Mars was initially on an eccentric and inclined orbit while the orbits of Mercury, Venus and Earth were more circular and coplanar. The lower primordial dynamical excitement and the peculiar partitioning between planets impose new constraints for terrestrial planet formation simulations.

The Transit Timing Variations (TTVs) can be used as a diagnostic of gravitational interactions between planets in a multi-planet system. Many Kepler Objects of Interest (KOIs) exhibit significant TTVs, but KOI-142.01 stands out among them with an unrivaled, 12-hour TTV amplitude. Here we report a thorough analysis of KOI-142.01's transits. We discover periodic Transit Duration Variations (TDVs) of KOI-142.01 that are nearly in phase with the observed TTVs. We show that KOI-142.01's TTVs and TDVs uniquely detect a non-transiting companion with a mass 0.7 that of Jupiter (KOI-142c). KOI-142.01's mass inferred from the transit variations is consistent with the measured transit depth, suggesting a Neptune class planet (KOI-142b). The orbital period ratio P_c/P_b=2.03 indicates that the two planets are just wide of the 2:1 resonance. The present dynamics of this system, characterized here in detail, can be used to test various formation theories that have been proposed to explain the near-resonant pairs of exoplanets.

Jupiter Trojans are thought to be survivors of a much larger population of planetesimals that existed in the planetary region when planets formed. They can provide important constraints on the mass and properties of the planetesimal disk, and its dispersal during planet migration. Here we tested a possibility that the Trojans were captured during the early dynamical instability among the outer planets (aka the Nice model), when the semimajor axis of Jupiter was changing as a result of scattering encounters with an ice giant. The capture occurs in this model when Jupiter's orbit and its Lagrange points become radially displaced in a scattering event and fall into a region populated by planetesimals (that previously evolved from their natal transplanetary disk to ~5 AU during the instability). Our numerical simulations of the new capture model, hereafter jump capture, satisfactorily reproduce the orbital distribution of the Trojans and their total mass. The jump capture is potentially capable of explaining the observed asymmetry in the number of leading and trailing Trojans. We find that the capture probability is (6-8)x10^-7 for each particle in the original transplanetary disk, implying that the disk contained (3-4)x10^7 planetesimals with absolute magnitude H<9 (corresponding to diameter D=80 km for a 7% albedo). The disk mass inferred from this work, M_disk=14-28 M_Earth, is consistent with the mass deduced from recent dynamical simulations of the planetary instability.

From the list of 2321 transiting planet candidates announced by the Kepler Mission, we select seven targets with favorable properties for the capacity to dynamically maintain an exomoon and present a detectable signal. These seven candidates were identified through our automatic target selection (TSA) algorithm and target selection prioritization (TSP) filtering, whereby we excluded systems exhibiting significant time-correlated noise and focussed on those with a single transiting planet candidate of radius less than 6 Earth radii. We find no compelling evidence for an exomoon around any of the seven KOIs but constrain the satellite-to-planet mass ratios for each. For four of the seven KOIs, we estimate a 95% upper quantile of M_S/M_P<0.04, which given the radii of the candidates, likely probes down to sub-Earth masses. We also derive precise transit times and durations for each candidate and find no evidence for dynamical variations in any of the KOIs. With just a few systems analyzed thus far in the in-going HEK project, projections on eta-moon would be premature, but a high frequency of large moons around Super-Earths/Mini-Neptunes would appear to be incommensurable with our results so far.

We analyze the distribution of extrasolar planets (both confirmed and Kepler candidates) according to their orbital periods P and planetary radii R. Among confirmed planets, we find compelling evidence for a paucity of bodies with 3 < R < 10 R_\oplus, where R_\oplus in the Earth's radius, and P < 2-3 days. We have christened this region a "sub-Jovian Pampas". The same trend is detected in multiplanet Kepler candidates. Although approximately 16 Kepler single-planet candidates inhabit this Pampas, at least 7 are probable false positives (FP). This last number could be significantly higher if the ratio of FP is higher than 10%, as suggested by recent studies. In a second part of the paper we analyze the distribution of planets in the (P,R) plane according to stellar metallicities. We find two interesting trends: (i) a lack of small planets (R < 4 R_\oplus) with orbital periods P < 5 days in metal-poor stars, and (ii) a paucity of sub-Jovian planets (4 R_\oplus < R < 8 R_\oplus) with P < 100 days, also around metal-poor stars. Although all these trends are preliminary, they appear statistically significant and deserve further scrutiny. If confirmed, they could represent important constraints on theories of planetary formation and dynamical evolution.

Understanding the growth of the cores of giant planets is a difficult problem. Recently, Lambrechts and Johansen (2012; LJ12) proposed a new model in which the cores grow by the accretion of pebble-size objects, as the latter drift towards the star due to gas drag. Here, we investigate the dynamics of pebble-size objects in the vicinity of planetary embryos of 1 and 5 Earth masses and the resulting accretion rates. We use hydrodynamical simulations, in which the embryo influences the dynamics of the gas and the pebbles suffer gas drag according to the local gas density and velocities. The pebble dynamics in the vicinity of the planetary embryo is non-trivial, and it changes significantly with the pebble size. Nevertheless, the accretion rate of the embryo that we measure is within an order of magnitude of the rate estimated in LJ12 and tends to their value with increasing pebble-size. We conclude that the model by LJ12 has the potential to explain the rapid growth of giant planet cores. The actual accretion rates however, depend on the surface density of pebble size objects in the disk, which is unknown to date.

Several properties of the Solar System, including the wide radial spacing and orbital eccentricities of giant planets, can be explained if the early Solar System evolved through a dynamical instability followed by migration of planets in the planetesimal disk. Here we report the results of a statistical study, in which we performed nearly 10^4 numerical simulations of planetary instability starting from hundreds of different initial conditions. We found that the dynamical evolution is typically too violent, if Jupiter and Saturn start in the 3:2 resonance, leading to ejection of at least one ice giant from the Solar System. Planet ejection can be avoided if the mass of the transplanetary disk of planetesimals was large (M_disk>50 M_Earth), but we found that a massive disk would lead to excessive dynamical damping (e.g., final e_55 < 0.01 compared to present e_55=0.044, where e_55 is the amplitude of the fifth eccentric mode in the Jupiter's orbit), and to smooth migration that violates constraints from the survival of the terrestrial planets. Better results were obtained when the Solar System was assumed to have five giant planets initially and one ice giant, with the mass comparable to that of Uranus and Neptune, was ejected into interstellar space by Jupiter. The best results were obtained when the ejected planet was placed into the external 3:2 or 4:3 resonance with Saturn and M_disk ~ 20 M_Earth. The range of possible outcomes is rather broad in this case, indicating that the present Solar System is neither a typical nor expected result for a given initial state, and occurs, in best cases, with only a ~5% probability (as defined by the success criteria described in the main text). The case with six giant planets shows interesting dynamics but does offer significant advantages relative to the five planet case.

NASA's Dawn spacecraft observations of asteroid (4) Vesta reveal a surface with the highest albedo and color variation of any asteroid we have observed so far. Terrains rich in low albedo dark material (DM) have been identified using Dawn Framing Camera (FC) 0.75 {\mu}m filter images in several geologic settings: associated with impact craters (in the ejecta blanket material and/or on the crater walls and rims); as flow-like deposits or rays commonly associated with topographic highs; and as dark spots (likely secondary impacts) nearby impact craters. This DM could be a relic of ancient volcanic activity or exogenic in origin. We report that the majority of the spectra of DM are similar to carbonaceous chondrite meteorites mixed with materials indigenous to Vesta. Using high-resolution seven color images we compared DM color properties (albedo, band depth) with laboratory measurements of possible analog materials. Band depth and albedo of DM are identical to those of carbonaceous chondrite xenolith-rich howardite Mt. Pratt (PRA) 04401. Laboratory mixtures of Murchison CM2 carbonaceous chondrite and basaltic eucrite Millbillillie also show band depth and albedo affinity to DM. Modeling of carbonaceous chondrite abundance in DM (1-6 vol%) is consistent with howardite meteorites. We find no evidence for large-scale volcanism (exposed dikes/pyroclastic falls) as the source of DM. Our modeling efforts using impact crater scaling laws and numerical models of ejecta reaccretion suggest the delivery and emplacement of this DM on Vesta during the formation of the ~400 km Veneneia basin by a low-velocity (<2 km/sec) carbonaceous impactor. This discovery is important because it strengthens the long-held idea that primitive bodies are the source of carbon and probably volatiles in the early Solar System.

The Kepler Mission is monitoring the brightness of ~150,000 stars searching for evidence of planetary transits. As part of the "Hunt for Exomoons with Kepler" (HEK) project, we report a planetary system with two confirmed planets and one candidate planet discovered using the publicly available data for KOI-872. Planet b transits the host star with a period P_b=33.6d and exhibits large transit timing variations indicative of a perturber. Dynamical modeling uniquely detects an outer nontransiting planet c near the 5:3 resonance (P_c=57.0d) of mass 0.37 times that of Jupiter. Transits of a third planetary candidate are also found: a 1.7-Earth radius super-Earth with a 6.8d period. Our analysis indicates a system with nearly coplanar and circular orbits, reminiscent of the orderly arrangement within the solar system.

Two decades ago, empirical evidence concerning the existence and frequency of planets around stars, other than our own, was absent. Since this time, the detection of extrasolar planets from Jupiter-sized to most recently Earth-sized worlds has blossomed and we are finally able to shed light on the plurality of Earth-like, habitable planets in the cosmos. Extrasolar moons may also be frequent habitable worlds but their detection or even systematic pursuit remains lacking in the current literature. Here, we present a description of the first systematic search for extrasolar moons as part of a new observational project called "The Hunt for Exomoons with Kepler" (HEK). The HEK project distills the entire list of known transiting planet candidates found by Kepler (2326 at the time of writing) down to the most promising candidates for hosting a moon. Selected targets are fitted using a multimodal nested sampling algorithm coupled with a planet-with-moon light curve modelling routine. By comparing the Bayesian evidence of a planet-only model to that of a planet-with-moon, the detection process is handled in a Bayesian framework. In the case of null detections, upper limits derived from posteriors marginalised over the entire prior volume will be provided to inform the frequency of large moons around viable planetary hosts, eta-moon. After discussing our methodologies for target selection, modelling, fitting and vetting, we provide two example analyses.

Exoplanets show a pile-up of Jupiter-size planets in orbits with a 3-day period. A fraction of these hot Jupiters have retrograde orbits with respect to the parent star's rotation. To explain these observations we performed a series of numerical integrations of planet scattering followed by the tidal circularization. We considered planetary systems having 3 and 4 planets initially. We found that the standard Kozai migration is an inefficient mechanism for the formation of hot Jupiters. Our results show the formation of two distinct populations of hot Jupiters. The inner population of hot Jupiters with semimajor axis a < 0.03 AU formed in the systems where no planetary ejections occurred. This group contained a significant fraction of highly inclined and retrograde orbits, with distributions largely independent of the initial setup. However, our follow-up integrations showed that this populations was transient with most planets falling inside the Roche radius of the star in <1 Gyr. The outer population of hot Jupiters formed in systems where at least one planet was ejected. This population survived the effects of tides over >1 Gyr. The semimajor axis distribution of Population II fits nicely the observed 3-day pile-up. The inclination distribution of the outer hot planets depends on the number of planets in the initial systems and the 4-planet case showed a larger proportion (up to 10%), and a wider spread in inclination values. As the later results roughly agrees with observations, this may suggest that the planetary systems with observed hot Jupiters were originally rich in the number of planets, some of which were ejected. In a broad perspective, our work therefore hints on an unexpected link between the hot Jupiters and recently discovered free floating planets.

Recent studies of solar system formation suggest that the solar system's giant planets formed and migrated in the protoplanetary disk to reach resonant orbits with all planets inside 15 AU from the Sun. After the gas disk's dispersal, Uranus and Neptune were likely scattered by gas giants, and approached their current orbits while dispersing the transplanetary disk of planetesimals, whose remains survived to this time in the region known as the Kuiper belt. Here we performed N-body integrations of the scattering phase between giant planets in an attempt to determine which initial states are plausible. We found that the dynamical simulations starting with a resonant system of four giant planets have a low success rate in matching the present orbits of giant planets, and various other constraints (e.g., survival of the terrestrial planets). The dynamical evolution is typically too violent, if Jupiter and Saturn start in the 3:2 resonance, and leads to final systems with fewer than four planets. Several initial states stand out in that they show a relatively large likelihood of success in matching the constraints. Some of the statistically best results were obtained when assuming that the solar system initially had five giant planets and one ice giant, with the mass comparable to that of Uranus and Neptune, was ejected to interstellar space by Jupiter. This possibility appears to be conceivable in view of the recent discovery of a large number free-floating planets in interstellar space, which indicates that planet ejection should be common.

The Oort Cloud Comets (OCCs), exemplified by the Great Comet of 1997 (Hale-Bopp), are occasional visitors from the heatless periphery of the solar system. Previous works hypothesized that a great majority of OCCs must physically disrupt after one or two passages through the inner solar system, where strong thermal gradients can cause phase transitions or volatile pressure buildup. Here we study the fate of small debris particles produced by OCC disruptions to determine whether the imprints of a hypothetical population of OCC meteoroids can be found in the existing meteor radar data. We find that OCC particles with diameters D<10 um are blown out from the solar system by radiation pressure, while those with D>1 um have a very low Earth-impact probability. The intermediate particle sizes, D=100 um, represent a sweet spot. About 1% of these particles orbitally evolve by Poynting-Robertson drag to reach orbits with semimajor axis a=1 AU. They are expected to produce meteors with radiants near the apex of the Earth's orbital motion. We find that the model distributions of their impact speeds and orbits provide a good match to radar observations of apex meteors, except for the eccentricity distribution, which is more skewed toward e=1 in our model. Finally, we propose an explanation for the long-standing problem in meteor science related to the relative strength of apex and helion/antihelion sources. As we show in detail, the observed trend, with the apex meteors being more prominent in observations of highly sensitive radars, can be related to orbital dynamics of particles released on the long-period orbits.

The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (>10^5 yr at 1 AU) than postulated in the standard collisional models (10^4 yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5)x10^11 km^2 and 4x10^19 g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be 10^4-10^5 kg/s. The input is up to 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate.

The size distribution in the cold classical Kuiper belt can be approximated by two idealized power laws: one with steep slope for radii R>R* and one with shallow slope for R<R*, where R*~25-50 km. Previous works suggested that the SFD roll-over at R* can be the result of extensive collisional grinding in the Kuiper belt that led to the catastrophic disruption of most bodies with R<R*. Here we use a new code to test the effect of collisions in the Kuiper belt. We find that the observed roll-over could indeed be explained by collisional grinding provided that the initial mass in large bodies was much larger than the one in the present Kuiper belt, and was dynamically depleted. In addition to the size distribution changes, our code also tracks the effects of collisions on binary systems. We find that it is generally easier to dissolve wide binary systems, such as the ones existing in the cold Kuiper belt today, than to catastrophically disrupt objects with R~R*. Thus, the binary survival sets important limits on the extent of collisional grinding in the Kuiper belt. We find that the extensive collisional grinding required to produce the SFD roll-over at R* would imply a strong gradient of the binary fraction with R and separation, because it is generally easier to dissolve binaries with small components and/or those with wide orbits. The expected binary fraction for R<R* is <0.1. The present observational data do not show such a gradient. Instead, they suggest a large binary fraction of ~0.4 for R=30-40 km. This may indicate that the roll-over was not produced by disruptive collisions, but is instead a fossil remnant of the KBO formation process.

A large fraction of 100-km-class low-inclination objects in the classical Kuiper Belt (KB) are binaries with comparable mass and wide separation of components. A favored model for their formation was capture during the coagulation growth of bodies in the early KB. Instead, recent studies suggested that large objects can rapidly form in the protoplanetary disks when swarms of locally concentrated solids collapse under their own gravity. Here we examine the possibility that KB binaries formed during gravitational collapse when the excess of angular momentum prevented the agglomeration of available mass into a solitary object. We find that this new mechanism provides a robust path toward the formation of KB binaries with observed properties, and can explain wide systems such as 2001 QW322 and multiples such as (47171) 1999 TC36. Notably, the gravitational collapse is capable of producing 100% binary fraction for a wide range of the swarm's initial angular momentum values. The binary components have similar masses (80% have the secondary-over-primary radius ratio >0.7) and their separation ranges from ~1,000 to ~100,000 km. The binary orbits have eccentricities from e=0 to ~1, with the majority having e<0.6. The binary orbit inclinations with respect to the initial angular momentum of the swarm range from i=0 to ~90 deg, with most cases having i<50 deg. Our binary formation mechanism implies that the primary and secondary components in each binary pair should have identical bulk composition, which is consistent with the current photometric data. We discuss the applicability of our results to the Pluto-Charon, Orcus-Vanth, (617) Patroclus-Menoetius and (90) Antiope binary systems.

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these `space weathering' processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r*, at which planetary encounters can reset the surface. We found that t_sw ~ 1e6 yr (at 1 AU) and r* ~ 5 R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 1e5 yr would require that r* > 20 R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 1e7 yr. We found that t_sw ~ q^2, where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.

The Transit Timing Variation (TTV) method relies on monitoring changes in timing of transits of known exoplanets. Non-transiting planets in the system can be inferred from TTVs by their gravitational interaction with the transiting planet. The TTV method is sensitive to low-mass planets that cannot be detected by other means. Here we describe a fast algorithm that can be used to determine the mass and orbit of the non-transiting planets from the TTV data. We apply our code, ttvim.f, to a wide variety of planetary systems to test the uniqueness of the TTV inversion problem and its dependence on the precision of TTV observations. We find that planetary parameters, including the mass and mutual orbital inclination of planets, can be determined from the TTV datasets that should become available in near future. Unlike the radial velocity technique, the TTV method can therefore be used to characterize the inclination distribution of multi-planet systems.

The zodiacal cloud is a thick circumsolar disk of small debris particles produced by asteroid collisions and comets. Here, we present a zodiacal cloud model based on the orbital properties and lifetimes of comets and asteroids, and on the dynamical evolution of dust after ejection. The model is quantitatively constrained by IRAS observations of thermal emission, but also qualitatively consistent with other zodiacal cloud observations. We find that 85-95% of the observed mid-infrared emission is produced by particles from the Jupiter-family comets (JFCs) and $<$10% by dust from long period comets. Asteroidal dust is found to be present at $<$10%. We suggest that spontaneous disruptions of JFCs, rather than the usual cometary activity driven by sublimating volatiles, is the main mechanism that librates cometary particles into the zodiacal cloud. Our results imply that JFC particles represent $\sim$85% of the total mass influx at Earth. Since their atmospheric entry speeds are typically low ($\approx$14.5 km s$^{-1}$ mean for D=100-200 $\mu$m with $\approx$12 km s$^{-1}$ being the most common case), many JFC grains should survive frictional heating and land on the Earth's surface. This explains why most micrometeorites collected in antarctic ice have primitive carbonaceous composition. The present mass of the inner zodiacal cloud at $<$5 AU is estimated to be 1-$2\times10^{19}$ g, mainly in D=100-200 $\mu$m particles. The inner zodiacal cloud should have been $>10^4$ times brighter during the Late Heavy Bombardment (LHB) epoch $\approx$3.8 Gyr ago, when the outer planets scattered numerous comets into the inner solar system. The bright debris disks with a large 24-$\mu$m excess observed around mature stars may be an indication of massive cometary populations existing in those systems.

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter D~100km in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ~100km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk ``jumped'' from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100km, probably even 1,000km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size-range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.

By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-Neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the `Nice model'. This leads us to predict that the HD of hot population should have the same slope of the HD of the cold population for 6.5 < H < 9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H~9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-Neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. (2004) and Fuentes and Holman (2008). Finally, we show that the constraint on the total mass of the primordial trans-Neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU.

We present the observational results of a survey designed to target and detect asteroids whose colors are similar to those of Vesta family members and thus may be considered as candidates for having a basaltic composition. Fifty basaltic candidates were selected with orbital elements that lie outside of the Vesta dynamical family. Optical and near-infrared spectra were used to assign a taxonomic type to 11 of the 50 candidates. Ten of these were spectroscopically confirmed as V-type asteroids, suggesting that most of the candidates are basaltic and can be used to constrain the distribution of basaltic material in the Main Belt. Using our catalog of V-type candidates and the success rate of the survey, we calculate unbiased size-frequency and semi-major axis distributions of V-type asteroids. These distributions, in addition to an estimate for the total mass of basaltic material, suggest that Vesta was the predominant contributor to the basaltic asteroid inventory of the Main Belt, however scattered planetesimals from the inner Solar System (a < 2.0 AU) and other partially/fully differentiated bodies likely contributed to this inventory. In particular, we infer the presence of basaltic fragments in the vicinity of asteroid 15 Eunomia, which may be derived from a differentiated parent body in the middle Main Belt (2.5 < a < 2.8). We find no asteroidal evidence for a large number of previously undiscovered basaltic asteroids, which agrees with previous theories suggesting that basaltic fragments from the ~100 differentiated parent bodies represented in meteorite collections have been "battered to bits" [Burbine, T.H., Meibom, A., Binzel, R.P., 1996. Mantle material in the Main Belt: Battered to bits? Met. & Planet. Sci. 31, 607].

Using 8m-10m class telescopes and their Adaptive Optics (AO) systems, we conducted a long-term adaptive optics campaign initiated in 2003 focusing on four binary asteroid systems: (130) Elektra, (283) Emma, (379) Huenna, and (3749) Balam. The analysis of these data confirms the presence of their asteroidal satellite. We did not detect any additional satellite around these systems even though we have the capability of detecting a loosely-bound fragment (located at 1/4 x RHill) ~40 times smaller in diameter than the primary. The orbits derived for their satellites display significant eccentricity, ranging from 0.1 to 0.9, suggesting a different origin. Based on AO size estimate, we show that (130) Elektra and (283) Emma, G-type and P-type asteroids respectively, have a significant porosity (30-60% considering CI-CO meteorites as analogs) and their satellite's eccentricities (e~0.1) are possibly due to excitation by tidal effects. (379) Huenna and (3749) Balam, two loosely bound binary systems, are most likely formed by mutual capture. (3749) Balam's possible high bulk density is similar to (433) Eros, another S-type asteroid, and should be poorly fractured as well. (379) Huenna seems to display both characteristics: the moonlet orbits far away from the primary in term of stability (20% x RHill), but the primary's porosity is significant (30-60%).

The recent discovery of the first V-type asteroid in the middle belt, (21238) 1995WV7, located at ~2.54 AU, raises the question of whether it came from (4) Vesta or not. In this paper, we present spectroscopic observations indicating the existence of another V-type asteroid at ~2.53 AU, (40521) 1999RL95, and we investigate the possibility that these two asteroids evolved from the Vesta family to their present orbits by drifting in semi-major axis due to the Yarkovsky effect. The main problem with this scenario is that the asteroids need to cross the 3/1 mean motion resonance with Jupiter, which is highly unstable. Combining numerical simulations of the orbital evolution, that include the Yarkovsky effect, with Monte Carlo models, we compute the probability of an asteroid of given diameter D to evolve from the Vesta family and to cross over the 3/1 resonance, reaching a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 1995WV7 has a low probability of having evolved through this mechanism due to its large size (~5 km). However, the mechanism might explain the orbit of smaller bodies like (40521) 1999RL95 (~3 km), provided that we assume that the Vesta family formed > 3.5 Gy ago. We estimate that about 10% or more of the V-type bodies with D>1 km may come from the Vesta family by crossing over the 3/1 resonance. The remaining 90% must have a different origin.

Given a primary interest in "mitigation of the potential hazard" of near-Earth objects impacting the Earth, the subject of characterization takes on an aspect not normally present when considering asteroids as abstract bodies. Many deflection concepts are interested in the classic geophysical characteristics of asteroids when considering the physical challenge of modifying their orbits in order to cause them to subsequently miss an impact with Earth. Yet for all deflection concepts there are characteristics of the threat which overwhelm these traditional factors. For example, a close gravitational encounter with Earth some years or decades prior to impact can reduce the velocity change necessary for deflection by several orders of magnitude if the deflection precedes the close encounter (or encounters). Conversely this "benefit" comes at a "price"; a corresponding increase in the accuracy of tracking required to determine the probability of impact. Societal issues, both national and international, also characterize the NEO deflection process and these may strongly contend with the purely technical issues normally considered. Therefore critical factors not normally considered must be brought into play as one characterizes the threat of NEO impacts.