In 2015, the New Horizons spacecraft flew past Pluto and its moon Charon, providing the first clear look at the surface of Charon. New Horizons images revealed an ancient surface, a large, intricate canyon system, and many fractures, among other geologic features. Here, we assess whether tidal stresses played a significant role in the formation of tensile fractures on Charon. Although presently in a circular orbit, most scenarios for the orbital evolution of Charon include an eccentric orbit for some period of time and possibly an internal ocean. Past work has shown that these conditions could have generated stresses comparable in magnitude to other tidally fractured moons, such as Europa and Enceladus. However, we find no correlation between observed fracture orientations and those predicted to form due to eccentricity-driven tidal stress. It thus seems more likely that the orbit of Charon circularized before its ocean froze, and that either tidal stresses alone were insufficient to fracture the surface or subsequent resurfacing remove these ancient fractures.

Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family members could contribute to the observed population of JFCs. Numerical integrations show that such dynamical evolution is indeed possible via a combination of eccentricity excitation apparently driven by the nearby 2:1 mean-motion resonance with Jupiter, gravitational interactions with planets other than Jupiter, and the Yarkovsky effect. We estimate that, at any given time, there may be tens of objects from the Themis family on JFC-like orbits with the potential to mimic active JFCs from the outer solar system, although not all, or even any, may necessarily be observably active. We find that dynamically evolved Themis family objects on JFC-like orbits have semimajor axes between 3.15 au and 3.40 au for the vast majority of their time on such orbits, consistent with the strong role that the 2:1 mean-motion resonance with Jupiter likely plays in their dynamical evolution. We conclude that a contribution from the Themis family to the active JFC population is plausible, although further work is needed to better characterize this contribution.

We present numerical simulations of terrestrial planet formation that examine the growth continuously from planetesimals to planets in the inner Solar System. Previous studies show that the growth will be inside-out, but it is still common practice to assume that the entire inner disk will eventually reach a bi-modal distribution of embryos and planetesimals. For the combinations of disk mass, initial planetesimal radius and gas disk lifetime explored in this work the entire disk never reaches a simple bi-modal mass distribution. We find that the inside-out growth is amplified by the combined effects of collisional evolution of solid bodies and interactions with a dissipating gas disk. This leads to oligarchic growth never being achieved in different places of the disk at the same time, where in some cases the disk can simultaneoulsy support chaotic growth and giant impacts inside 1 au and runaway growth beyond 2 au. The planetesimal population is efficiently depleted in the inner disk where embryo growth primarily advances in the presence of a significant gas disk. Further out in the disk growth is slower relative to the gas disk dissipation, resulting in more excited planetesimals at the same stage of growth and less efficient accretion. This same effect drives mass loss due to collisional grinding strongly altering the surface density of the accreted planets relative to the initial mass distribution. This effect decreases the Mars-to-Earth mass ratios compared to previous works with no collisional grinding. Similar to some previous findings utilizing vastly different growth scenarios these simulations produce a first generation of planetary embryos that are stable for 10-20 Myr, or 5-10 e-folding times of the gas dissipation timescale, before having an instability and entering the chaotic growth stage.

We measure the void probability function (VPF) of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). The VPF provides complementary information to standard two-point statistics in that it is sensitive to galaxy bias in the most extreme underdensities in the cosmic web. Thus the VPF is ideal for testing whether halo occupation of galaxies depends on large-scale density, an effect known as galaxy assembly bias. We find that standard HOD model---one parameterized by halo mass only---fit only to the two-point function, accurately predicts the VPF. Additionally, for HOD models where density dependence is explicitly incorporated, the best-fit models fit to the combination of the correlation function and the VPF have zero density dependence. Thus galaxy assembly bias is not a strong source of systematic uncertaintiy when modeling the clustering of massive galaxies.

The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the "Early Instability" scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets' orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar system's terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system.

The OSIRIS-REx Asteroid Sample Return Mission is the third mission in NASA's New Frontiers Program and is the first U.S. mission to return samples from an asteroid to Earth. The most important decision ahead of the OSIRIS-REx team is the selection of a prime sample-site on the surface of asteroid (101955) Bennu. Mission success hinges on identifying a site that is safe and has regolith that can readily be ingested by the spacecraft's sampling mechanism. To inform this mission-critical decision, the surface of Bennu is mapped using the OSIRIS-REx Camera Suite and the images are used to develop several foundational data products. Acquiring the necessary inputs to these data products requires observational strategies that are defined specifically to overcome the challenges associated with mapping a small irregular body. We present these strategies in the context of assessing candidate sample-sites at Bennu according to a framework of decisions regarding the relative safety, sampleability, and scientific value across the asteroid's surface. To create data products that aid these assessments, we describe the best practices developed by the OSIRIS-REx team for image-based mapping of irregular small bodies. We emphasize the importance of using 3D shape models and the ability to work in body-fixed rectangular coordinates when dealing with planetary surfaces that cannot be uniquely addressed by body-fixed latitude and longitude.

The moniker rubble pile is typically applied to all solar system bodies with Diameter between 200m and 10km - where in this size range there is an abundance of evidence that nearly every object is bound primarily by self-gravity with significant void space or bulk porosity between irregularly shaped constituent particles. The understanding of this population is derived from wide-ranging population studies of derived shape and spin, decades of observational studies in numerous wavelengths, evidence left behind from impacts on planets and moons and the in situ study of a few objects via spacecraft flyby or rendezvous. The internal structure, however, which is responsible for the name rubble pile, is never directly observed, but belies a violent history. Many or most of the asteroids on near-Earth orbits, and the ones most accessible for rendezvous and in situ study, are likely byproducts of the continued collisional evolution of the Main Asteroid Belt.

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

Many dynamical aspects of the solar system can be explained by the outer planets experiencing a period of orbital instability sometimes called the Nice Model. Though often correlated with a perceived delayed spike in the lunar cratering record known as the Late Heavy Bombardment (LHB), recent work suggests that this event may have occurred much earlier; perhaps during the epoch of terrestrial planet formation. While current simulations of terrestrial accretion can reproduce many observed qualities of the solar system, replicating the small mass of Mars requires modification to standard planet formation models. Here we use 800 dynamical simulations to show that an early instability in the outer solar system strongly influences terrestrial planet formation and regularly yields properly sized Mars analogs. Our most successful outcomes occur when the terrestrial planets evolve an additional 1-10 million years (Myr) following the dispersal of the gas disk, before the onset of the giant planet instability. In these simulations, accretion has begun in the Mars region before the instability, but the dynamical perturbation induced by the giant planets' scattering removes large embryos from Mars' vicinity. Large embryos are either ejected or scattered inward toward Earth and Venus (in some cases to deliver water), and Mars is left behind as a stranded embryo. An early giant planet instability can thus replicate both the inner and outer solar system in a single model.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

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

The thermal properties of the surfaces of asteroids determine the magnitude of the drift rate cause by the Yarkovsky force. In the general case of Main Belt asteroids, the Yarkovsky force is indirectly proportional to the thermal inertia, $\Gamma$. Following the proposed relationship between $\Gamma$ and asteroid diameter $D$, we find that asteroids' Yarkovsky drift rates might have a more complex size dependence than previous thought, leading to a curved family V-shape boundary in semi-major axis, a, vs. 1/$D$ space. This implies that asteroids are drifting faster at larger sizes than previously considered decreasing on average the known ages of asteroid families. The V-Shape curvature is determined for $>$25 families located throughout the Main Belt to quantify the Yarkovsky size-dependent drift rate. We find that there is no correlation between family age and V-shape curvature. In addition, the V-shape curvature decreases for asteroid families with larger heliocentric distances suggesting that the relationship between $\Gamma$ and $D$ is weaker in the outer MB possibly due to homogenous surface roughness among family members.

Ejection velocity fields of asteroid families are largely unconstrained due to the fact that members disperse relatively quickly on Myr time-scales by secular resonances and the Yarkovsky effect. The spreading of fragments in $a$ by the Yarkovsky effect is indistinguishable from the spreading caused by the initial ejection of fragments. By examining families $<$20 Myrs-old, we can use the V-shape identification technique to separate family shapes that are due to the initial ejection velocity field and those that are due to the Yarkovsky effect. $<$20 Myr-old asteroid families provide an opportunity to study the velocity field of family fragments before they become too dispersed. Only the Karin family's initial velocity field has been determined and scales inversely with diameter, $D^{-1}$. We have applied the V-shape identification technique to constrain young families' initial ejection velocity fields by measuring the curvature of their fragments' V-shape correlation in semi-major axis, $a$, vs. $D^{-1}$ space. Curvature from a straight line implies a deviation from a scaling of $D^{-1}$. We measure the V-shape curvature of 11 young asteroid families including the \FYnospace, Aeolia, Brangane, Brasilia, Clarissa, Iannini, Karin, Konig, Koronis(2), Theobalda and Veritas asteroid families. We find that the majority of asteroid families have initial ejection velocity fields consistent with $\sim D^{-1}$ supporting laboratory impact experiments and computer simulations of disrupting asteroid parent bodies.

Until 2013, only the giant planets were known to host ring systems. In June 2013, a stellar occulation revealed the presence of narrow and dense rings around Chariklo, a small Centaur object that orbits between Saturn and Uranus. Meanwhile, the Cassini spacecraft revealed evidence for the possible past presence of rings around the Saturnian satellites Rhea and Iapetus. Mars and Pluto are expected to have tenuous dusty rings, though they have so far evaded detection. More remotely, transit events observed around a star in 2007 may have revealed for the first time exoplanetary rings around a giant planet orbiting that star. So, evidence is building to show that rings are more common features in the universe than previously thought. Several interesting issues arising from the discovery (or suspicion) of new ring systems are described in this chapter.

Standard approaches to Bayesian parameter inference in large scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as Approximate Bayesian Computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter halos with galaxies. Using specific implementation of ABC supplemented with Population Monte Carlo importance sampling, a generative forward model using HOD, and a distance metric based on galaxy number density, two-point correlation function, and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected "true" HOD parameters. The parameter constraints we obtain from ABC are consistent with the "true" HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately our results suggest that ABC can and should be applied in parameter inference for LSS analyses.

In May of 2011, NASA selected the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) asteroid sample return mission as the third mission in the New Frontiers program. The other two New Frontiers missions are New Horizons, which explored Pluto during a flyby in July 2015 and is on its way for a flyby of Kuiper Belt object 2014 MU69 on Jan. 1, 2019, and Juno, an orbiting mission that is studying the origin, evolution, and internal structure of Jupiter. The spacecraft departed for near-Earth asteroid (101955) Bennu aboard an United Launch Alliance Atlas V 411 evolved expendable launch vehicle at 7:05 p.m. EDT on September 8, 2016, on a seven-year journey to return samples from Bennu. The spacecraft is on an outbound-cruise trajectory that will result in a rendezvous with Bennu in August 2018. The science instruments on the spacecraft will survey Bennu to measure its physical, geological, and chemical properties, and the team will use these data to select a site on the surface to collect at least 60 g of asteroid regolith. The team will also analyze the remote-sensing data to perform a detailed study of the sample site for context, assess Bennus resource potential, refine estimates of its impact probability with Earth, and provide ground-truth data for the extensive astronomical data set collected on this asteroid. The spacecraft will leave Bennu in 2021 and return the sample to the Utah Test and Training Range (UTTR) on September 24, 2023.

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.

It has been shown that some aspects of the terrestrial planets can be explained, particularly the Earth/Mars mass ratio, when they form from a truncated disk with an outer edge near 1.0 au (Hansen 2009). This has been previously modeled starting from an intermediate stage of growth utilizing pre-formed planetary embryos. We present simulations that were designed to test this idea by following the growth process from km-sized objects located between 0.7 to 1.0 au up to terrestrial planets. The simulations explore initial conditions where the solids in the disk are planetesimals with radii initially between 3 and 300 km, alternately including effects from a dissipating gaseous solar nebula and collisional fragmentation. We use a new Lagrangian code known as LIPAD (Levison et al. 2012), which is a particle-based code that models the fragmentation, accretion and dynamical evolution of a large number of planetesimals, and can model the entire growth process from km-sizes up to planets. A suite of large (Mars mass) planetary embryos is complete in only 1 Myr, containing most of the system mass. A quiescent period then persists for 10-20 Myr characterized by slow diffusion of the orbits and continued accretion of the remaining planetesimals. This is interupted by an instability that leads to embryos crossing orbits and embyro-embryo impacts that eventually produce the final set of planets. While this evolution is different than that found in other works exploring an annulus, the final planetary systems are similar, with roughly the correct number of planets and good Mars-analogs.

There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Bro\v{z} et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity, $e$, and inclination, $i$, due to secular resonances, close encounters with massive asteroids and the non-gravitational Yarkovsky force. This causes the family fragments to be indistinguishable from the background of the main belt making them more difficult to identify with the hierarchical clustering method (HCM) with increasing family age. The discovery of the Eulalia and new Polana families in the inner belt relied on new techniques because Yarkovsky spreading made them too disperse to be identified using the classical HCM. The technique used to discover the new Polana and Eulalia families is modified to identify asteroid families by searching for correlations between $a$ and asteroid diameter, $D$, or absolute magnitude, $H$. A group of asteroids is identified as a collisional family if its boundary in the $a$ vs. 1/$D$ or $a$ vs. $H$ planes has a characteristic V-shape which is due to the size dependent Yarkovsky spreading. The V-shape boundary is identified with two separate techniques. The first identifies a border by measuring the drop between the number of objects inside and outside of the border. The second identifies the V-shape border by measuring the number density of objects in $a$ vs. 1/$D$,$H$ space. The V-shape techniques are demonstrated on known families Erigone, Vesta, Koronis, and families difficult to identify by HCM such as Flora, Baptistina, new Polana, Eulalia and Karin. Application of the technique in a search for >2 Gyr-old families throughout the MB is discussed.

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

Directed energy for planetary defense is now a viable option and is superior in many ways to other proposed technologies, being able to defend the Earth against all known threats. This paper presents basic ideas behind a directed energy planetary defense system that utilizes laser ablation of an asteroid to impart a deflecting force on the target. A conceptual philosophy called DE-STAR, which stands for Directed Energy System for Targeting of Asteroids and exploRation, is an orbiting stand-off system, which has been described in other papers. This paper describes a smaller, stand-on system known as DE-STARLITE as a reduced-scale version of DE-STAR. Both share the same basic heritage of a directed energy array that heats the surface of the target to the point of high surface vapor pressure that causes significant mass ejection thus forming an ejection plume of material from the target that acts as a rocket to deflect the object. This is generally classified as laser ablation. DE-STARLITE uses conventional propellant for launch to LEO and then ion engines to propel the spacecraft from LEO to the near-Earth asteroid (NEA). During laser ablation, the asteroid itself provides the propellant source material; thus a very modest spacecraft can deflect an asteroid much larger than would be possible with a system of similar mission mass using ion beam deflection (IBD) or a gravity tractor. DE- STARLITE is capable of deflecting an Apophis-class (325 m diameter) asteroid with a 1- to 15-year targeting time (laser on time) depending on the system design. The mission fits within the rough mission parameters of the Asteroid Redirect Mission (ARM) program in terms of mass and size. DE-STARLITE also has much greater capability for planetary defense than current proposals and is readily scalable to match the threat. It can deflect all known threats with sufficient warning.

Laser ablation of a Near-Earth Object (NEO) on a collision course with Earth produces a cloud of ejecta which exerts a thrust on the NEO, deflecting it from its original trajectory. Ablation may be performed from afar by illuminating an Earth-targeting asteroid or comet with a stand-off "DE- STAR" system consisting of a large phased-array laser in Earth orbit. Alternatively, a much smaller stand-on "DE-STARLITE" system may travel alongside the target, slowly deflecting it from nearby over a long period. This paper presents orbital simulations comparing the effectiveness of both systems across a range of laser and NEO parameters. Simulated parameters include magnitude, duration and, for the stand-on system, direction of the thrust, as well as the type, size and orbital characteristics of the target NEO. These simulations indicate that deflection distance is approximately proportional to the magnitude of thrust and to the square of the duration of ablation, and is inversely proportional to the mass. Furthermore, deflection distance shows strong dependence on thrust direction with the optimal direction of thrust varying with the duration of laser activity. As one example, consider a typical 325 m asteroid: beginning 15 yr in advance, just 2 N of thrust from a ~20 kW stand-on DE-STARLITE system is sufficient to deflect the asteroid by 2 R_e. Numerous scenarios are discussed as is a practical implementation of such a system consistent with current launch vehicle capabilities.

Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 AU should roughly be the same mass as the Earth. Recently, a new model called Viscous Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants. This model envisions that the cores of the giant planets formed from 100 to 1000 km bodies that directly accreted a population of pebbles --- sub-meter sized objects that slowly grew in the protoplanetary disk. Here we apply this model to the terrestrial planet region and find that it can reproduce the basic structure of the inner Solar System, including a small Mars and a low-mass asteroid belt. Our models show that for an initial population of planetesimals with sizes similar to those of the main belt asteroids, VSPA becomes inefficient beyond $\sim\!$1.5 AU. As a result, Mars's growth is stunted and nothing large in the asteroid belt can accumulate.

Satellites of asteroids have been discovered in nearly every known small body population, and a remarkable aspect of the known satellites is the diversity of their properties. They tell a story of vast differences in formation and evolution mechanisms that act as a function of size, distance from the Sun, and the properties of their nebular environment at the beginning of Solar System history and their dynamical environment over the next 4.5 Gyr. The mere existence of these systems provides a laboratory to study numerous types of physical processes acting on asteroids and their dynamics provide a valuable probe of their physical properties otherwise possible only with spacecraft. Advances in understanding the formation and evolution of binary systems have been assisted by: 1) the growing catalog of known systems, increasing from 33 to nearly 250 between the Merline et al. (2002) Asteroids III chapter and now, 2) the detailed study and long-term monitoring of individual systems such as 1999 KW4 and 1996 FG3, 3) the discovery of new binary system morphologies and triple systems, 4) and the discovery of unbound systems that appear to be end-states of binary dynamical evolutionary paths. Specifically for small bodies (diameter smaller than 10 km), these observations and discoveries have motivated theoretical work finding that thermal forces can efficiently drive the rotational disruption of small asteroids. Long-term monitoring has allowed studies to constrain the system's dynamical evolution by the combination of tides, thermal forces and rigid body physics. The outliers and split pairs have pushed the theoretical work to explore a wide range of evolutionary end-states.

The past decade has brought major improvements in large-scale asteroid discovery and characterization with over half a million known asteroids and over 100,000 with some measurement of physical characterization. This explosion of data has allowed us to create a new global picture of the Main Asteroid Belt. Put in context with meteorite measurements and dynamical models, a new and more complete picture of Solar System evolution has emerged. The question has changed from "What was the original compositional gradient of the Asteroid Belt?" to "What was the original compositional gradient of small bodies across the entire Solar System?" No longer is the leading theory that two belts of planetesimals are primordial, but instead those belts were formed and sculpted through evolutionary processes after Solar System formation. This article reviews the advancements on the fronts of asteroid compositional characterization, meteorite measurements, and dynamical theories in the context of the heliocentric distribution of asteroid compositions seen in the Main Belt today. This chapter also reviews the major outstanding questions relating to asteroid compositions and distributions and summarizes the progress and current state of understanding of these questions to form the big picture of the formation and evolution of asteroids in the Main Belt. Finally, we briefly review the relevance of asteroids and their compositions in their greater context within our Solar System and beyond.

Pluto's system of 5 known satellites are in a puzzling orbital configuration. Each of the four small satellites are on low-eccentricity and low-inclination orbits situated near a mean motion resonance with the largest satellite Charon. The Pluto-Charon binary likely formed as a result of a giant impact and so the simplest explanation for the small satellites is that they accreted from debris of that collision. The Pluto-Charon binary has evolved outward since its formation due to tidal forces, which drove them into their current doubly synchronous state. Meanwhile, leftover debris from the formation of Charon was not initially distant enough from Pluto-Charon to explain the orbits of the current small satellites. The outstanding problems of the system are the movement of debris outward and the small satellites location near mean motion resonances with Charon. This work explores the dynamical behavior of collisionally interacting debris orbiting the Pluto-Charon system. While this work specifically tests initial disk and ring configurations designed to mimic the aftermath of the disruption of satellites by heliocentric impactors, we generally find that collisional interactions can help move material outwards and keep otherwise unstable material dynamically bound to the Pluto-Charon system. These processes can produce rings of debris whose orbits evolve rapidly due to collisional processes, with increasing pericenters and decreasing semimajor axes. While these rings and disks of debris eventually build satellites significantly further out than the initial locations of a disrupted satellite, they do not show a strong preference for building satellites in or near mean motion resonances with Charon under a wide array of tested conditions.

According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and initiated the final phase of core formation by melting the Earth's mantle. A key goal of geochemistry is to date this event, but different ages have been proposed. Some argue for an early Moon-forming event, approximately 30 million years (Myr) after the condensation of the first solids in the Solar System, whereas others claim a date later than 50 Myr (and possibly as late as around 100 My) after condensation. Here we show that a Moon-forming event at 40 Myr after condensation, or earlier, is ruled out at a 99.9 per cent confidence level. We use a large number of N-body simulations to demonstrate a relationship between the time of the last giant impact on an Earth-like planet and the amount of mass subsequently added during the era known as Late Accretion. As the last giant impact is delayed, the late-accreted mass decreases in a predictable fashion. This relationship exists within both the classical scenario and the Grand Tack scenario of terrestrial planet formation, and it holds across a wide range of disk conditions. The concentration of highly siderophile elements (HSEs) in Earth's mantle constrains the mass of chondritic material added to Earth during Late Accretion. Using HSE abundance measurements, we determine a Moon-formation age of 95 +/- 32 Myr since condensation. The possibility exists that some late projectiles were differentiated and left an incomplete HSE record in Earth's mantle. Even in this case, various isotopic constraints strongly suggest that the late-accreted mass did not exceed 1 per cent of Earth's mass, and so the HSE clock still robustly limits the timing of the Moon-forming event to significantly later than 40 My after condensation.

The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with the nebula gas lead to remarkable new accretion mechanisms for the formation of planetesimals and the growth of planetary embryos. Many of the popular models for the later stages of planet formation are presented. The classical models with the giant planets on fixed orbits are not consistent with the known history of the Solar System, fail to create a high Earth/Mars mass ratio, and, in many cases, are also internally inconsistent. The successful Grand Tack model creates a small Mars, a wet Earth, a realistic asteroid belt and the mass-orbit structure of the terrestrial planets. In the Grand Tack scenario, growth curves for Earth most closely match a Weibull model. The feeding zones, which determine the compositions of Earth and Venus follow a particular pattern determined by Jupiter, while the feeding zones of Mars and Theia, the last giant impactor on Earth, appear to randomly sample the terrestrial disk. The late accreted mass samples the disk nearly evenly.

The asteroid belt is the leftover of the original planetesimal population in the inner solar system. However, currently the asteroids have orbits with all possible values of eccentricities and inclinations compatible with long-term dynamical stability, whereas the initial planetesimal orbits should have been quasi-circular and almost co-planar. The total mass in the asteroid population is a small fraction of that existing primordially. Also, asteroids with different chemical/mineralogical properties are not ranked in an orderly manner with mean heliocentric distance as one could expect from the existence of a radial gradient of the temperature in the proto-planetary disk, but they are partially mixed. These properties show that the asteroid belt has been severely sculpted by one or a series of processes during its lifetime. This paper reviews the processes that have been proposed so far, discussing the properties that they explain and the problems that they are confronted with. Emphasis is paid to the interplay between the dynamical and the collisional evolution of the asteroid population, which allows the use of the size distribution to constrain the dynamical models. We divide the asteroid belt evolution into three phases. The first phase started during the lifetime of the gaseous proto-planetary disk, when the giant planets formed and presumably experienced large-scale migrations, and continued after the removal of the gas, during the build-up of the terrestrial planets. The second phase occurred after the removal of the gaseous proto-planetary disk and it became particularly lively for the asteroid belt when the giant planets suddenly changed their orbits, as a result of a mutual dynamical instability and the interaction with the trans-Neptunian planetesimal disk. The third phase covers the aftermath of the giant planet instability, until today.

We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations.

A new model for terrestrial planet formation (Hansen 2009, Walsh et al. 2011) has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanism to truncate the disk--a two-stage, inward-then-outward migration of Jupiter and Saturn, as found in numerous hydrodynamical simulations of giant planet formation. In addition to truncating the disk and producing a more realistic Earth/Mars mass ratio, the migration of the giant planets also populates the asteroid belt with two distinct populations of bodies--the inner belt is filled by bodies originating inside of 3 AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as `primitive' bodies). We find here that the planets will accrete on order 1-2% of their total mass from primitive planetesimals scattered onto planet-crossing orbits during the formation of the planets. For an assumed value of 10% for the water mass fraction of the primitive planetesimals, this model delivers a total amount of water comparable to that estimated to be on the Earth today. While the radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system, we find that the last giant impact is typically earlier than 20 Myr, and a substantial amount of mass is accreted after that event. However, 5 of the 27 planets larger than half an Earth mass formed in all simulations do experience large late impacts and subsequent accretion consistent with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following that event.

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.

We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\bar\tau_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\bar\tau_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population.

On 2013 November 28 Comet ISON (C/2012 S1) will pass by the Sun with a perihelion distance of 2.7 solar radii. Understanding the possible outcomes for the comet's response to such a close passage by the Sun is important for planning observational campaigns and for inferring ISON's physical properties. We present new numerical simulations and interpret them in context with the historical track record of comet disruptions and of sungrazing comet behavior. Historical data suggest that sizes below ~200 m are susceptible to destruction by sublimation driven mass loss, while we find that for ISON's perihelion distance, densities lower than 0.1 g cm^-3 are required to tidally disrupt a retrograde or non-spinning body. Such low densities are substantially below the range of the best-determined comet nucleus densities, though dynamically new comets such as ISON have few measurements of physical properties. Disruption may occur for prograde rotation at densities up to 0.7 g cm^-3, with the chances of disruption increasing for lower density, faster prograde rotation, and increasing elongation of the nucleus. Given current constraints on ISON's nucleus properties and the typically determined values for these properties among all comets, we find tidal disruption to be unlikely unless other factors (e.g., spin-up via torquing) affect ISON substantially. Whether or not disruption occurs, the largest remnant must be big enough to survive subsequent mass loss due to sublimation in order for ISON to remain a viable comet well after perihelion.

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

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 so-called Nysa-Polana complex of asteroids is a diverse and widespread group. It appears to be two overlapping families of different asteroid taxonomies: (44) Nysa is an E-type asteroid with the lowest number in the midst of a predominantly S-type cluster and (142) Polana is a B-type asteroid near the low-albedo B- and C-type cluster. Using the data from the Wide-field Infrared Survey Explorer (WISE) mission we have re-analyzed the region around the Nysa-Polana complex in the inner Main Belt, focusing on the low-albedo population. (142) Polana does not appear to be a member of the family of low-albedo asteroids in the Nysa-Polana complex. Rather, the largest is asteroid (495) Eulalia. This asteroid has never before been linked to this complex for an important dynamical reason: it currently has a proper eccentricity slightly below the range of most of the family members. However, its orbit is very close to the 3:1 mean motion resonance with Jupiter and is in a weak secular resonance. We show that its osculating eccentricity varies widely on short timescales and the averaged value diffuses over long timescales. The diffusive orbit, low-albedo, taxonomic similarity and semimajor axis strongly suggests that despite its current proper eccentricity, (495) Eulalia could have recently been at an orbit very central to the family. Hierarchical Clustering Method tests confirm that at an eccentricity of e=0.15, (495) Eulalia could be the parent of the family. The ``Eulalia family'' was formed between 900--1500 Myr ago, and likely resulted from the breakup of a 100--160 km parent body. There is also compelling evidence for an older and more widespread primitive family in the same region of the asteroid belt parented by asteroid (142) Polana.

This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth.

The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interesting spacecraft target: a binary asteroid with a low-DeltaV heliocentric orbit. The orbit of its satellite has provided valuable information about its mass density while its albedo and colors suggest it is primitive or part of the C-complex taxonomic grouping. We extend the physical characterization of this object with new observations of its emission at mid-Infrared (IR) wavelengths and with near-IR reflection spectroscopy. We derive an area-equivalent system diameter of 1.90 \pm 0.28 km (corresponding to approximate component diameters of 1.83 km and 0.51 km, respectively) and a geometric albedo of 0.039 \pm 0.012. 1996 FG3 was previously classified as a C-type asteroid, though the combined 0.4--2.5 micron spectrum with thermal correction indicates classification as B-type; both are consistent with the low measured albedo. Dynamical studies show that 1996 FG3 has most probably originated in the inner main asteroid belt. Recent work has suggested the inner Main Belt (142) Polana family as the possible origin of another low-DeltaV B-type NEA, (101955) 1999 RQ36. A similar origin for 1996 FG3 would require delivery by the overlapping Jupiter 7:2 and Mars 5:9 mean motion resonances rather than the nu-6 resonance, and we find this to be a low probability, but possible, origin.

Jupiter and Saturn formed in a few million years (Haisch et al. 2001) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ~100,000 years (Armitage 2007). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration (Masset & Snellgrove 2001, Morbidelli & Crida 2007, Pierens & Nelson 2008). The terrestrial planets finished accreting much later (Klein et al. 2009), and their characteristics, including Mars' small mass, are best reproduced by starting from a planetesimal disk with an outer edge at about one astronomical unit from the Sun (Wetherill 1978, Hansen 2009) (1 AU is the Earth-Sun distance). Here we report simulations of the early Solar System that show how the inward migration of Jupiter to 1.5 AU, and its subsequent outward migration, lead to a planetesimal disk truncated at 1 AU; the terrestrial planets then form from this disk over the next 30-50 million years, with an Earth/Mars mass ratio consistent with observations. Scattering by Jupiter initially empties but then repopulates the asteroid belt, with inner-belt bodies originating between 1 and 3 AU and outer-belt bodies originating between and beyond the giant planets. This explains the significant compositional differences across the asteroid belt. The key aspect missing from previous models of terrestrial planet formation is the substantial radial migration of the giant planets, which suggests that their behaviour is more similar to that inferred for extrasolar planets than previously thought.

It has long been suspected that metal polluted white dwarfs (types DAZ, DBZ, and DZ) and white dwarfs with dusty disks possess planetary systems, but a specific physical mechanism by which planetesimals are perturbed close to a white dwarf has not yet been fully posited. In this paper we demonstrate that mass loss from a central star during post main sequence evolution can sweep planetesimals into interior mean motion resonances with a single giant planet. These planetesimals are slowly removed through chaotic excursions of eccentricity that in time create radial orbits capable of tidally disrupting the planetesimal. Numerical N-body simulations of the Solar System show that a sufficient number of planetesimals are perturbed to explain white dwarfs with both dust and metal pollution, provided other white dwarfs have more massive relic asteroid belts. Our scenario requires only one Jupiter-sized planet and a sufficient number of asteroids near its 2:1 interior mean motion resonance. Finally, we show that once a planetesimal is perturbed into a tidal crossing orbit, it will become disrupted after the first pass of the white dwarf, where a highly eccentric stream of debris forms the main reservoir for dust producing collisions. These simulations, in concert with observations of white dwarfs, place interesting limits on the frequency of planetary systems around main sequence stars, the frequency of planetesimal belts, and the probability that dust may obscure future terrestrial planet finding missions.

We present a scenario for building the equatorial ridge and de-spinning Iapetus through an impact-generated disk and satellite. This impact puts debris into orbit, forming a ring inside the Roche limit and a satellite outside. This satellite rapidly pushes the ring material down to the surface of Iapetus, and then itself tidally evolves outward, thereby helping to de-spin Iapetus. This scenario can de-spin Iapetus an order of magnitude faster than when tides due to Saturn act alone, almost independently of its interior geophysical evolution. Eventually, the satellite is stripped from its orbit by Saturn. The range of satellite and impactor masses required is compatible with the estimated impact history of Iapetus.

We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the Martian obliquity is stable and oscillates with an amplitude lower than 20$^\circ$. These zones occur at obliquities below 30$^\circ$ and above 60$^\circ$; intermediate values show either resonant or chaotic behaviour depending on the primordial orbits of the terrestrial planets.

The migration of the giant planets due to the scattering of planetesimals causes powerful resonances to move through the asteroid belt and the terrestrial planet region. Exactly when and how the giant planets migrated is not well known. In this paper we present results of an investigation of the formation of the terrestrial planets during and after the migration of the giant planets. The latter is assumed to have occurred immediately after the dissipation of the nebular disk -- i.e. "early" with respect to the timing of the Late Heavy Bombardment (LHB). The presumed cause of our modeled early migration of the giant planets is angular mometum transfer between the planets and scattered planetesimals.

Asteroids are leftover pieces from the era of planet formation that help us understand conditions in the early Solar System. Unlike larger planetary bodies that were subject to global thermal modification during and subsequent to their formation, these small bodies have kept at least some unmodified primordial material from the solar nebula. However, the structural properties of asteroids have been modified considerably since their formation. Thus, we can find among them a great variety of physical configurations and dynamical histories. In fact, with only a few possible exceptions, all asteroids have been modified or completely disrupted many times during the age of the Solar System. This picture is supported by data from space mission encounters with asteroids that show much diversity of shape, bulk density, surface morphology, and other features. Moreover, the gravitational attraction of these bodies is so small that some physical processes occur in a manner far removed from our common experience on Earth. Thus, each visit to a small body has generated as many questions as it has answered. In this review we discuss the current state of research into asteroid disruption processes, focusing on collisional and rotational mechanisms. We find that recent advances in modeling catastrophic disruption by collisions have provided important insights into asteroid internal structures and a deeper understanding of asteroid families. Rotational disruption, by tidal encounters or thermal effects, is responsible for altering many smaller asteroids, and is at the origin of many binary asteroids and oddly shaped bodies.