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.

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 present new Hubble Space Telescope and ground-based Keck observations and new Keplerian orbit solutions for the mutual orbit of binary Jupiter Trojan asteroid (617) Patroclus and Menoetius, targets of NASA's Lucy mission. We predict event times for the upcoming mutual event season, which is anticipated to run from late 2017 through mid 2019.

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 Trojan clouds contain a surprisingly diverse population of objects. It is only by understanding this diversity can we unravel their history. The more objects Lucy visits the better we can constrain the formation and evolution of the giant planets. While Lucy will visit four Trojans in the L$_4$ cloud, its trajectory only passes one known Trojan system, the PM binary, in L$_5$. There is ample time in the Lucy timeline to possibly add additional L$_5$ targets, thereby maximizing the mission's science return. In the next decade, the 8.4-m Large Synoptic Survey Telescope (LSST), is expected to detect nearly 300,000 Jupiter Trojans between approximately 16 and 24.5 r-band magnitudes in the $\sim$18,000 deg$^2$ main survey and the expected $\sim$4,000 deg$^2$ northern ecliptic extension. We explore potential opportunities for LSST to find new Lucy L$_5$ fly-by targets.

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.

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.

It is widely held that the first step in forming the gas giant planets, such as Jupiter and Saturn, is to form solid `cores' of roughly 10 M$_\oplus$. Getting the cores to form before the solar nebula dissipates ($\sim\!1-10\,$Myr) has been a major challenge for planet formation models. Recently models have emerged in which `pebbles' (centimeter- to meter-size objects) are first concentrated by aerodynamic drag and then gravitationally collapse to form 100 --- 1000 km objects. These `planetesimals' can then efficiently accrete leftover pebbles and directly form the cores of giant planets. This model known as `pebble accretion', theoretically, can produce 10 M$_\oplus$ cores in only a few thousand years. Unfortunately, full simulations of this process show that, rather than creating a few 10 M$_\oplus$ cores, it produces a population of hundreds of Earth-mass objects that are inconsistent with the structure of the Solar System. Here we report that this difficulty can be overcome if pebbles form slowly enough to allow the planetesimals to gravitationally interact with one another. In this situation the largest planetesimals have time to scatter their smaller siblings out of the disk of pebbles, thereby stifling their growth. Our models show that, for a large, and physically reasonable region of parameter space, this typically leads to the formation of one to four gas giants between 5 and 15 AU in agreement with the observed structure of the Solar System.

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.

When the EPOXI spacecraft flew by Comet 103P/Hartley 2, it observed large particles floating around the comet nucleus. These particles are likely low-density, centimeter- to decimeter-sized clumps of ice and dust. While the origin of these objects remains somewhat mysterious, it is possible that they are giving us important information about the earliest stages of our Solar System's formation. Recent advancements in planet formation theory suggest that planetesimals (or cometestimals) may grow directly from the gravitational collapse of aerodynamically concentrated small particles, often referred to as "pebbles." Here we show that the particles observed in the coma of 103P are consistent with the sizes of pebbles expected to efficiently form planetesimals in the region that this comet likely formed, while smaller pebbles are may be expected in the majority of comets, whose chemistry is often indicative of formation in the colder, outer regions of the protoplanetary disk.

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.

From observations by the Hubble Space Telescope, Keck II Telescope, and Gemini North Telescope, we have determined the mutual orbit of the large transneptunian object (174567) Varda and its satellite Ilmar\"e. These two objects orbit one another in a highly inclined, circular or near-circular orbit with a period of 5.75 days and a semimajor axis of 4810 km. This orbit reveals the system mass to be (2.664 +/- 0.064) x 10^20 kg, slightly greater than the mass of the second most massive main-belt asteroid (4) Vesta. The dynamical mass can in turn be combined with estimates of the surface area of the system from Herschel Space Telescope thermal observations to estimate a bulk density of 1.24 +0.50 -0.35 g cm^-3. Varda and Ilmar\"e both have colors similar to the combined colors of the system, B-V = 0.886 +/- 0.025 and V-I = 1.156 +/- 0.029.

Though ~10 Earth mass rocky/icy cores are commonly held as a prerequisite for the formation of gas giants, theoretical models still struggle to explain how these embryos can form within the lifetimes of gaseous circumstellar disks. In recent years, aerodynamic-aided accretion of "pebbles," objects ranging from centimeters to meters in size, has been suggested as a potential solution to this long-standing problem. While pebble accretion has been demonstrated to be extremely effective in local simulations that look at the detailed behavior of these pebbles in the vicinity of a single planetary embryo, to date there have been no global simulations demonstrating the effectiveness of pebble accretion in a more complicated, multi-planet environment. Therefore, we have incorporated the aerodynamic-aided accretion physics into LIPAD, a Lagrangian code that can follow the collisional / accretional / dynamical evolution of a protoplanetary system, to investigate how pebble accretion manifests itself in the larger planet formation picture. We find that under generic circumstances, pebble accretion naturally leads to an "oligarchic" type of growth in which a large number of planetesimals grow to similar-sized planets. In particular, our simulations tend to form hundreds of Mars- and Earth-mass objects between 4 and 10 AU. While the merging of some oligarchs may allow them to grow massive enough to form giant planet cores, leftover oligarchs lead to planetary systems that cannot be consistent with our own solar system. We investigate various ideas presented in the literature (including evaporation fronts and planet traps) and find that none easily overcome this tendency toward oligarchic growth.

We presented the first particle based, Lagrangian code that can follow the collisional/accretional/dynamical evolution of a large number of km-sized planetesimals through the entire growth process to become planets. We refer to it as the 'Lagrangian Integrator for Planetary Accretion and Dynamics' or LIPAD. LIPAD is built on top of SyMBA, which is a symplectic $N$-body integrator. In order to handle the very large number of planetesimals required by planet formation simulations, we introduce the concept of a `tracer' particle. Each tracer is intended to represent a large number of disk particles on roughly the same orbit and size as one another, and is characterized by three numbers: the physical radius, the bulk density, and the total mass of the disk particles represented by the tracer. We developed statistical algorithms that follow the dynamical and collisional evolution of the tracers due to the presence of one another. The tracers mainly dynamically interact with the larger objects (`planetary embryos') in the normal N-body way. LIPAD's greatest strength is that it can accurately model the wholesale redistribution of planetesimals due to gravitational interaction with the embryos, which has recently been shown to significantly affect the growth rate of planetary embryos . We verify the code via a comprehensive set of tests which compare our results with those of Eulerian and/or direct N-body codes.

Hubble Space Telescope observations between 2001 and 2010 resolved the binary components of the Cold Classical transneptunian object (79360) Sila-Nunam (provisionally designated 1997 CS29). From these observations we have determined the circular, retrograde mutual orbit of Nunam relative to Sila with a period of 12.50995 \pm 0.00036 days and a semimajor axis of 2777 \pm 19 km. A multi-year season of mutual events, in which the two near-equal brightness bodies alternate in passing in front of one another as seen from Earth, is in progress right now, and on 2011 Feb. 1 UT, one such event was observed from two different telescopes. The mutual event season offers a rich opportunity to learn much more about this barely-resolvable binary system, potentially including component sizes, colors, shapes, and albedo patterns. The low eccentricity of the orbit and a photometric lightcurve that appears to coincide with the orbital period are consistent with a system that is tidally locked and synchronized, like the Pluto-Charon system. The orbital period and semimajor axis imply a system mass of (10.84 \pm 0.22) \times 10^18 kg, which can be combined with a size estimate based on Spitzer and Herschel thermal infrared observations to infer an average bulk density of 0.72 +0.37 -0.23 g cm^-3, comparable to the very low bulk densities estimated for small transneptunian binaries of other dynamical classes.

Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own protosolar nebula using the properties of the small body reservoirs in the solar system. In standard planet formation theory, after Jupiter and Saturn formed they scattered a significant number of remnant planetesimals into highly eccentric orbits. In this paper, we show that if there had been a massive, extended protoplanetary disk at that time, then the disk would have excited Kozai oscillations in some of the scattered objects, driving them into high-inclination (i > 50 deg), low-eccentricity orbits (q > 30 AU). The dissipation of the gaseous disk would strand a subset of objects in these high-inclination orbits; orbits that are stable on Gyr time scales. To date, surveys have not detected any Kuiper Belt Objects with orbits consistent with this dynamical mechanism. Using these non-detections by the Deep Ecliptic Survey (DES) and the Palomar Distant Solar System Survey we are able to rule out an extended gaseous protoplanetary disk (R_D > 80 AU) in our solar system at the time of Jupiter's formation. Future deep all sky surveys such as the Large Synoptic Survey Telescope (LSST) will all us to further constrain the size of the protoplanetary disk.

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser et al., 2006) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster -- either a Hernquist or Plummer potential -- and chose other parameters based on the latest observations of embedded clusters from the literature. {These other parameters result in the stars being on radial orbits and the cluster collapses.} Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100\,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a `concentration radius' for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters...

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 present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK183, (149780) Altjira, 2001 QY297, 2003 QW111, and 2003 QY90 based on Hubble Space Telescope and Keck 2 laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits yields a sample of 22 mutual orbits for which the period P, semimajor axis a, and eccentricity e have been determined. These orbits have mutual periods ranging from 5 to over 800 days, semimajor axes ranging from 1,600 to 37,000 km, eccentricities ranging from 0 to 0.8, and system masses ranging from 2 x 10^17 to 2 x 10^22 kg. Based on the relative brightnesses of primaries and secondaries, most of these systems consist of near equal-sized pairs, although a few of the most massive systems are more lopsided. The observed distribution of orbital properties suggests that the most loosely-bound transneptunian binary systems are only found on dynamically cold heliocentric orbits. Of the 22 known binary mutual or-bits, orientation ambiguities are now resolved for 9, of which 7 are prograde and 2 are retro-grade, consistent with a random distribution of orbital orientations, but not with models predicting a strong preference for retrograde orbits. To the extent that other perturbations are not dominant, the binary systems undergo Kozai oscillations of their eccentricities and inclinations with periods of the order of tens of thousands to millions of years, some with strikingly high amplitudes.

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.

We report here on an extension of a previous study by Kirsh et al. (2009) of planetesimal-driven migration using our N-body code SyMBA (Duncan et al., 1998). The previous work focused on the case of a single planet of mass Mem, immersed in a planetesimal disk with a power-law surface density distribution and Rayleigh distributed eccentricities and inclinations. Typically 10^4-10^5 equal-mass planetesimals were used, where the gravitational force (and the back-reaction) on each planetesimal by the Sun and planetwere included, while planetesimal-planetesimal interactions were neglected. The runs reported on here incorporate the dynamical effects of a gas disk, where the Adachi et al. (1976) prescription of aerodynamic gas drag is implemented for all bodies. In some cases the Papaloizou and Larwood (2000) prescription of Type-I migration for the planet are implemented, as well as a mass distribution. In the gas-free cases, rapid planet migration was observed - at a rate independent of the planet's mass - provided the planet's mass was not large compared to the mass in planetesimals capable of entering its Hill sphere. In such cases, both inward and outward migrations can be self-sustaining, but there is a strong propensity for inward migration. When a gas disk is present, aerodynamic drag can substantially modify the dynamics of scattered planetesimals. For sufficiently large or small mono-dispersed planetesimals, the planet typically migrates inward. However, for a range of plausible planetesimal sizes (i.e. 0.5-5.0 km at 5.0 AU in a minimum mass Hayashi disk) outward migration is usually triggered, often accompanied by substantial planetary mass accretion. The origins of this behaviour are explained in terms of a toy model. The effects of including a size distribution and torques associated with Type-I migration are also discussed.

We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale of tau ~ 0.5My. However, self-consistent numerical simulations show that the planetesimal-driven migration of the giant planets is inconsistent with an exponential change in their semi major axes on such a short time scale (Hahn & Malhotra, 1999). In fact, the typical time scale is tau > 5My. When giant planet migration on this time scale is applied to the asteroid belt, the resulting orbital distribution is incompatible with the observed one. However, the planet migration can be significantly sped up by planet-planet encounters. Consider an evolution where both Jupiter and Saturn have close encounters with a Neptune-mass planet (presumably Uranus or Neptune themselves) and where this third planet, after being scattered inwards by Saturn, is scattered outwards by Jupiter. This scenario leads to a very rapid increase in the orbital separation between Jupiter and Saturn that we show here to have only mild effects on the structure of asteroid belt. This type of evolution is called a jumping-Jupiter case. Our results suggest that the total mass and dynamical excitation of the asteroid belt before migration were comparable to those currently observed. Moreover, they imply that, before migration, the orbits of Jupiter and Saturn were much less eccentric than the current ones.

One of the most challenging problems we face in our understanding of planet formation is how Jupiter and Saturn could have formed before the the solar nebula dispersed. The most popular model of giant planet formation is the so-called 'core accretion' model. In this model a large planetary embryo formed first, mainly by two-body accretion. This is then followed by a period of inflow of nebular gas directly onto the growing planet. The core accretion model has an Achilles heel, namely the very first step. We have undertaken the most comprehensive study of this process to date. In this study we numerically integrate the orbits of a number of planetary embryos embedded in a swarm of planetesimals. In these experiments we have included: 1) aerodynamic gas drag, 2) collisional damping between planetesimals, 3) enhanced embryo cross-sections due to their atmospheres, 4) planetesimal fragmentation, and 5) planetesimal driven migration. We find that the gravitational interaction between the embryos and the planetesimals lead to the wholesale redistribution of material - regions are cleared of material and gaps open near the embryos. Indeed, in 90% of our simulations without fragmentation, the region near that embryos is cleared of planetesimals before much growth can occur. The remaining 10%, however, the embryos undergo a burst of outward migration that significantly increases growth. On timescales of ~100,000 years, the outer embryo can migrate ~6 AU and grow to roughly 30 Earth-masses. We also find that the inclusion of planetesimal fragmentation tends to inhibit growth.

We present new analysis of HST images of (47171) 1999 TC36 that confirm it as a triple system. Fits to the point-spread function consistently show that the apparent primary is itself composed of two similar-sized components. The two central components, A1 and A2, can be consistently identified in each of nine epochs spread over seven years of time. In each instance the component separation, ranging from 0.023+/-0.002 to 0.031+/-0.003 arcsec, is roughly one half of the Hubble Space Telescope's diffraction limit at 606 nm. The orbit of the central pair has a semi-major axis of a~867 km with a period of P~1.9 days. These orbital parameters yield a system mass that is consistent with Msys = 12.75+/-0.06 10^18 kg derived from the orbit of the more distant secondary, component B. The diameters of the three components are dA1= 286(+45,-38) km, dA2= 265(+41,-35 km and dB= 139(+22,-18) km. The relative sizes of these components are more similar than in any other known multiple in the solar system. Taken together, the diameters and system mass yield a bulk density of p=542(+317,-211) kg m^-3. HST Photometry shows that component B is variable with an amplitude of >=0.17+/-0.05 magnitudes. Components A1 and A2 do not show variability larger than 0.08+/-0.03 magnitudes approximately consistent with the orientation of the mutual orbit plane and tidally-distorted equilibrium shapes. The system has high specific angular momentum of J/J'=0.93, comparable to most of the known Transneptunian binaries.

The habitability of planets is strongly affected by impacts from comets and asteroids. Indications from the ages of Moon rocks suggest that the inner Solar System experienced an increased rate of impacts roughly 3.8 Gya known as the Late Heavy Bombardment (LHB). Here we develop a model of how the Solar System would have appeared to a distant observer during its history based on the Nice model of Gomes et al. (2005). We compare our results with observed debris discs. We show that the Solar System would have been amongst the brightest of these systems before the LHB. Comparison with the statistics of debris disc evolution shows that such heavy bombardment events must be rare occurring around less than 12% of Sun-like stars.

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.

Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations. Moreover, the amplitudes of the eigenmodes characterising the current secular evolution of the eccentricities of Jupiter and Saturn would not be reproduced correctly; only one eigenmode is excited by resonance-crossing. We show that, at the very least, encounters between Saturn and one of the ice giants (Uranus or Neptune) need to have occurred, in order to reproduce the current secular properties of the giant planets, in particular the amplitude of the two strongest eigenmodes in the eccentricities of Jupiter and Saturn.

We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes the g5 eigenfrequency to cross resonances of the form g5=gk with k ranging from 1 to 4. Consequently these secular resonances cause large-amplitude responses in the eccentricities of the terrestrial planets. We show that the resonances g5=g_4 and g5=g3 do not pose a problem if Jupiter and Saturn have a fast approach and departure from their mutual 2:1 mean motion resonance. On the other hand, the resonance crossings g5=g2 and g5=g1 are more of a concern as they tend to yield a terrestrial system incompatible with the current one. We offer two solutions to this problem. The first uses the fact that a secular resonance crossing can also damp the amplitude of a Fourier mode if the latter is large originally. A second scenario involves a 'jumping Jupiter' in which encounters between an ice giant and Jupiter, without ejection of the former, cause the latter to migrate away from Saturn much faster than if migration is driven solely by encounters with planetesimals. In this case, the g5=g2 and g5=g1 resonances can be jumped over, or occur very briefly.

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.

The Nice model of Gomes et al. (2005) suggests that the migration of the giant planets caused a planetesimal clearing event which led to the Late Heavy Bombardment (LHB) at 880 Myr. Here we investigate the IR emission from the Kuiper belt during the history of the Solar System as described by the Nice model. We describe a method for easily converting the results of n-body planetesimal simulations into observational properties (assuming black-body grains and a single size distribution) and further modify this method to improve its realism (using realistic grain properties and a three-phase size distribution). We compare our results with observed debris discs and evaluate the plausibility of detecting an LHB-like process in extrasolar systems. Recent surveys have shown that 4% of stars exhibit 24 um excess and 16% exhibit 70 um excess. We show that the Solar System would have been amongst the brightest of these systems before the LHB at both 24 and 70 um. We find a significant increase in 24 um emission during the LHB, which rapidly drops off and becomes undetectable within 30 Myr, whereas the 70 um emission remains detectable until 360 Myr after the LHB. Comparison with the statistics of debris disc evolution shows that such depletion events must be rare occurring around less than 12% of Sun-like stars and with this level of incidence we would expect approximately 1 of the 413 Sun-like, field stars so far detected to have a 24 um excess to be currently going through an LHB. We also find that collisional processes are important in the Solar System before the LHB and that parameters for weak Kuiper belt objects are inconsistent with the Nice model interpretation of the LHB.

We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question.

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.

Extrasolar multiple-planet systems provide valuable opportunities for testing theories of planet formation and evolution. The architectures of the known multiple-planet systems demonstrate a fascinating level of diversity, which motivates the search for additional examples of such systems in order to better constrain their formation and dynamical histories. Here we describe a comprehensive investigation of 22 planetary systems in an effort to answer three questions: 1) Are there additional planets? 2) Where could additional planets reside in stable orbits? and 3) What limits can these observations place on such objects? We find no evidence for additional bodies in any of these systems; indeed, these new data do not support three previously announced planets (HD 20367b: Udry et al. 2003, HD 74156d: Bean et al. 2008, and 47 UMa c: Fischer et al. 2002). The dynamical simulations show that nearly all of the 22 systems have large regions in which additional planets could exist in stable orbits. The detection-limit computations indicate that this study is sensitive to close-in Neptune-mass planets for most of the systems targeted. We conclude with a discussion on the implications of these non-detections.

We present Hubble Space Telescope observations of six binary transneptunian systems: 2000 QL251, 2003 TJ58, 2001 XR254, 1999 OJ4, (134860) 2000 OJ67, and 2004 PB108. The mutual orbits of these systems are found to have periods ranging from 22 to 137 days, semimajor axes ranging from 2360 to 10500 km, and eccentricities ranging from 0.09 to 0.55. These orbital parameters enable estimation of system masses ranging from 0.2 to 9.7 x 10+18 kg. For reasonable assumptions of bulk density (0.5 to 2.0 g cm-3), the masses can be combined with visible photometry to constrain sizes and albedos. The resulting albedos are consistent with an emerging picture of the dynamically "Cold" Classical subpopulation having relatively high albedos, compared with comparablysized objects on more dynamically excited orbits.

We report resolved photometry of the primary and secondary components of 23 transneptunian binaries obtained with the Hubble Space Telescope. V-I colors of the components range from 0.7 to 1.5 with a median uncertainty of 0.06 magnitudes. The colors of the primaries and secondaries are correlated with a Spearman rank correlation probability of 99.99991%, 5 sigma for a normal distribution. Fits to the primary vs. secondary colors are identical to within measurement uncertainties. The color range of binaries as a group is indistinguishable from that of the larger population of apparently single transneptunian objects. Whatever mechanism produced the colors of apparently single TNOs acted equally on binary systems. The most likely explanation is that the colors of transneptunian objects and binaries alike are primordial and indicative of their origin in a locally homogeneous, globally heterogeneous protoplanetary disk.

The recent discovery of the 2003EL61 collisional family in the Kuiper belt (Brown et al., 2007) is surprising because the formation of such a family is a highly improbable event in today's belt. Assuming Brown et al.'s estimate of the size of progenitors, we find that the probability that a Kuiper belt object was involved in such a collision since primordial times is less than roughly 0.001. In addition, it is not possible for the collision to have occurred in a massive primordial Kuiper belt because the dynamical coherence of the family would not have survived whatever event produced the currently observed orbital excitation. Here we suggest that the family is the result of a collision between two scattered disk objects. We show that the probability that a collision occurred between two such objects with sizes similar to those advocated in Brown et al. (2007) and that the center of mass of the resulting family is on an orbit typical of the Kuiper belt can be as large as 47%. Given the large uncertainties involved in this estimate, this result is consistent with the existence of one such family. If true, this result has implications far beyond the origin of a single collisional family, because it shows that collisions played an important role in shaping the dynamical structure of the small body populations that we see today.

We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971 +/- 0.006 days and a semimajor axis of 1628 +/- 29 km, implying a system mass of (9.49 +/- 0.52) x 10^17 kg. The eccentricity of the orbit is 0.526 +/- 0.015. Combined with a radiometric size determined from Spitzer Space Telescope data and the assumption that Typhon and Echidna both have the same albedo, we estimate that their radii are 76 +14/-16 and 42 +8/-9 km, respectively. These numbers give an average bulk density of only 0.44 +0.44/-0.17 g cm^-3, consistent with very low bulk densities recently reported for two other small transneptunian binaries.

We have searched 101 Classical transneptunian objects for companions with the Hubble Space Telescope. Of these, at least 21 are binary. The heliocentric inclinations of the objects we observed range from 0.6-34 degrees. We find a very strong anticorrelation of binaries with inclination. Of the 58 targets that have inclinations of less than 5.5 degrees, 17 are binary, a binary fraction of 29 +7/-6%. All 17 are similar-brightness systems. On the contrary, only 4 of the 42 objects with inclinations greater than 5.5 degrees have satellites and only 1 of these is a similar-brightness binary. This striking dichotomy appears to agree with other indications that the low eccentricity, non-resonant Classical transneptunian objects include two overlapping populations with significantly different physical properties and dynamical histories.

We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ~35 AU to the Kuiper belt without excessive inclination excitation. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ~27 AU and an eccentricity of ~0.3, and that subsequently Neptune's eccentricity damped in ~1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level.

We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonance (i.e. each planet is in resonance with its neighbor). In total we found 6 such configurations. For the gas disk parameters found in MC07, these configurations are characterized by a negligible migration rate. After the disappearance of the gas, and in absence of planetesimals, only two of these six configurations (the least compact ones) are stable for a time of hundreds of millions of years or more. The others become unstable on a timescale of a few My. Our preliminary simulations show that, when a planetesimal disk is added beyond the orbit of the outermost planet, the planets can evolve from the most stable of these configurations to their current orbits in a fashion qualitatively similar to that described in Tsiganis et al. (2005).

Long time coverage and high radial velocity precision have allowed for the discovery of additional objects in known planetary systems. Many of the extrasolar planets detected have highly eccentric orbits, which raises the question of how likely those systems are to host additional planets. We investigate six systems which contain a very eccentric (e>0.6) planet: HD 3651, HD 37605, HD 45350, HD 80606, HD 89744, and 16 Cyg B. We present updated radial-velocity observations and orbital solutions, search for additional planets, and perform test particle simulations to find regions of dynamical stability. The dynamical simulations show that short-period planets could exist in the HD 45350 and 16 Cyg B systems, and we use the observational data to set tight detection limits, which rule out additional planets down to a few Neptune masses in the HD 3651, HD 45350, and 16 Cyg B systems.

Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX128) reveal it to be a close binary system. The mutual orbit has a period of 9.554 +/- 0.011 days and a semimajor axis of 1840 +/- 48 km. These values enable computation of a system mass of (5.41 +/- 0.42) 10^18 kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 microns are combined with visible photometry to constrain the system's effective radius (109 +10/-11 km) and geometric albedo (0.084 +0.021/-0.014). We estimate the average bulk density to be 1.37 +0.66/-0.32 g cm^-3, consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM165, and (47171) 1999 TC36, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity <= 0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the solar system.

This chapter discusses the dynamical properties of the Kuiper belt population. Then, it focuses on the characteristics of the Kuiper belt that cannot be explained by its evolution in the framework of the current solar system. We review models of primordial solar system evolution that have been proposed to reproduce the Kuiper belt features, outlining advantages and problems of each of them.

Chiang et al. 2006, hereafter C06 have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ~5 primordial ice giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich et al. (2004a,b, hereafter G04). Here we test this hypothesis with a series of numerical simulations using a new code designed to incorporate the dynamical effects of collisions. We find that we cannot reproduce the observed Solar System. In particular, G04 and C06 argue that during the instability, all but two of the ice giants would be ejected from the Solar System by Jupiter and Saturn, leaving Uranus and Neptune behind. We find that ejections are actually rare and that instead the systems spread outward. This always leads to a configuration with too many planets that are too far from the Sun. Thus, we conclude that both G04's scheme for the formation of Uranus and Neptune and C06's Kuiper belt formation scenario are not viable in their current forms.

We have identified a binary companion to (42355) 2002 CR46 in our ongoing deep survey using the Hubble Space Telescope's High Resolution Camera. It is the first companion to be found around an object in a non-resonant orbit that crosses the orbits of giant planets. Objects in orbits of this kind, the Centaurs, have experienced repeated strong scattering with one or more giant planets and therefore the survival of binaries in this transient population has been in question. Monte Carlo simulations suggest, however, that binaries in (42355) 2002 CR46 -like heliocentric orbits have a high probability of survival for reasonable estimates of the binary's still-unknown system mass and separation. Because Centaurs are thought to be precursors to short period comets, the question of the existence of binary comets naturally arises; none has yet been definitively identified. The discovery of one binary in a sample of eight observed by HST suggests that binaries in this population may not be uncommon.

Planetary embryos embedded in a gas disc suffer a decay in semimajor axis -- type I migration -- due to the asymmetric torques produced by the interior and exterior wakes raised by the body (Goldreich & Tremaine 1980; Ward 1986). This presents a challenge for standard oligarchic approaches to forming the terrestrial planets (Kokubo & Ida 1998) as the timescale to grow the progenitor objects near 1 AU is longer than that for them to decay into the Sun. In this paper we investigate the middle and late stages of oligarchic growth using both semi-analytic methods (based upon Thommes et al. 2003) and N-body integrations, and vary gas properties such as dissipation timescale in different models of the protoplanetary disc. We conclude that even for near-nominal migration efficiencies and gas dissipation timescales of ~1 Myr it is possible to maintain sufficient mass in the terrestrial region to form Earth and Venus if the disc mass is enhanced by factors of ~2-4 over the minimum mass model. The resulting configurations differ in several ways from the initial conditions used in previous simulations of the final stages of terrestrial accretion (e.g. Chambers 2001), chiefly in (1) larger inter-embryo spacings, (2) larger embryo masses, and (3) up to ~0.4 Earth masses of material left in the form of planetesimals when the gas vanishes. The systems we produce are reasonably stable for ~100 Myr and therefore require an external source to stir up the embryos sufficiently to produce final systems resembling the terrestrial planets.

Explaining the origin of the orbit of 2000 CR105 (a ~ 230AU, q ~ 45AU) is a major test for our understanding of the primordial evolution of the outer Solar System. Gladman et al. (2001) showed that this objects could not have been a normal member of the scattered disk that had its perihelion distance increased by chaotic diffusion. In this paper we explore four seemingly promising mechanisms for explaining the origin of the orbit of this peculiar object: (i) the passage of Neptune through a high-eccentricity phase, (ii) the past existence of massive planetary embryos in the Kuiper belt or the scattered disk, (iii) the presence of a massive trans-Neptunian disk at early epochs which exerted tides on scattered disk objects, and (iv) encounters with other stars. Of all these mechanisms, the only one giving satisfactory results is the passage of a star. Indeed, our simulations show that the passage of a solar mass star at about 800 AU only perturbs objects with semi-major axes larger than roughly 200 AU to large perihelion distances. This is in good agreement with the fact that 2000 CR105 has a semi-major axis of 230AU and no other bodies with similar perihelion distances but smaller semi-major axes have yet been discovered. The discovery of 2003 VB12, (a=450AU, q=75AU) announced a few days before the submission of this paper, strengthen our conclusions.

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as ``oligarchic growth.'' Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk ($\sim$ 10$\times$ minimum-mass) is required to produce giant planet core-sized bodies ($\sim$ 10 M$_{\oplus}$) within the lifetime of the nebular gas ($\la$ 10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of non-interacting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System's two outer giant planets at their current locations on {\it any} timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas, and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.

USNO-B is an all-sky catalog that presents positions, proper motions, magnitudes in various optical passbands, and star/galaxy estimators for 1,042,618,261 objects derived from 3,643,201,733 separate observations. The data were obtained from scans of 7,435 Schmidt plates taken for the various sky surveys during the last 50 years. USNO-B1.0 is believed to provide all-sky coverage, completeness down to V = 21, 0.2 arcsecond astrometric accuracy at J2000, 0.3 magnitude photometric accuracy in up to five colors, and 85% accuracy for distinguishing stars from non-stellar objects. A brief discussion of various issues is given here, but the actual data are available from http://www.nofs.navy.mil and other sites.

The outer giant planets, Uranus and Neptune, pose a challenge to theories of planet formation. They exist in a region of the Solar System where long dynamical timescales and a low primordial density of material would have conspired to make the formation of such large bodies ($\sim$ 15 and 17 times as massive as the Earth, respectively) very difficult. Previously, we proposed a model which addresses this problem: Instead of forming in the trans-Saturnian region, Uranus and Neptune underwent most of their growth among proto-Jupiter and -Saturn, were scattered outward when Jupiter acquired its massive gas envelope, and subsequently evolved toward their present orbits. We present the results of additional numerical simulations, which further demonstrate that the model readily produces analogues to our Solar System for a wide range of initial conditions. We also find that this mechanism may partly account for the high orbital inclinations observed in the Kuiper belt.

We present a new analysis of the currently available orbital elements for the known Kuiper belt objects. In the non-resonant, main Kuiper belt we find a statistically significant relationship between an object's absolute magnitude (H) and its inclination (i). Objects with H<6.5 (i.e. radii >~170 km for a 4% albedo) have higher inclinations than those with H>6.5 (radii >~ 170 km). We have shown that this relationship is not caused by any obvious observational bias. We argue that the main Kuiper belt consists of the superposition of two distinct distributions. One is dynamically hot with inclinations as large as \~35 deg and absolute magnitudes as bright as 4.5; the other is dynamically cold with i<~5 deg and H>6.5. The dynamically cold population is most likely dynamically primordial. We speculate on the potential causes of this relationship.

Our current knowledge of the dynamical structure of the Kuiper Belt is reviewed here. Numerical results on long term orbital evolution and dynamical mechanisms underlying the transport of objects out of the Kuiper Belt are discussed. Scenarios about the origin of the highly non-uniform orbital distribution of Kuiper Belt objects are described, as well as the constraints these provide on the formation and long term dynamical evolution of the outer Solar system. Possible mechanisms include an early history of orbital migration of the outer planets, a mass loss phase in the outer Solar system and scattering by large planetesimals. The origin and dynamics of the scattered component of the Kuiper Belt is discussed. Inferences about the primordial mass distribution in the trans-Neptune region are reviewed. Outstanding questions about Kuiper Belt dynamics are listed.

The limiting magnitude of the HST data set used by Cochran et al. (1995) to detect small objects in the Kuiper belt is reevaluated, and the methods used are described in detail. It is shown, by implanting artificial objects in the original HST images, and re-reducing the images using our original algorithm, that the limiting magnitude of our images (as defined by the 50% detectability limit) is $V=28.4$. This value is statistically the same as the value found in the original analysis. We find that $\sim50%$ of the moving Kuiper belt objects with $V=27.9$ are detected when trailing losses are included. In the same data in which these faint objects are detected, we find that the number of false detections brighter than $V=28.8$ is less than one per WFPC2 image. We show that, primarily due to a zero-point calibration error, but partly due to inadequacies in modeling the HST'S data noise characteristics and Cochran et al.'s reduction techniques, Brown et al. 1997 underestimate the SNR of objects in the HST dataset by over a factor of 2, and their conclusions are therefore invalid.

Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations in timescale, it is desirable to use a variable timestep. However, naively varying the timestep destroys the desirable properties of symplectic integrators. We discuss briefly the idea that choosing the timestep in a time symmetric manner can improve the performance of variable timestep integrators. Then we present a symplectic integrator which is based on decomposing the force into components and applying the component forces with different timesteps. This multiple timescale symplectic integrator has all the desirable properties of the constant timestep symplectic integrators.

We report the statistical detection (at the >99% confidence level) of a population of 28th magnitude objects exhibiting proper motions of ~1 arcsecond per hour at quadrature in deep HST/WFPC2 images. The drift directions imply a preponderance of objects on prograde orbits concentrated near the ecliptic plane. We interpret this as the detection of objects which reside in the Kuiper belt near or beyond the orbit of Neptune, comparable in size to comet 1P/Halley (radii ~10 km for an albedo of 4%). Our observations imply a population of ~2.5x10^4 objects deg^-2 with V<28.6. This must be viewed as a lower limit until our search of parameter space is complete. Our observations imply that there are >2x10^8 objects of this size in the Kuiper belt with inclinations <12 degrees and within ~40 AU of the Sun.