We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in $\sim20\%$ of the cases on a wide range of timescales (up to $5\times10^9$ orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multi-planet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet-planet collisions in $\sim80\%$ of Kepler-102 simulations. Kepler-102 b likely has a mass $>\sim0.1M_{\oplus}$, hence a bulk density exceeding about half Earth's, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system's AMD evolution with a "spectral fraction" calculated from the power spectrum of short integrations ($\sim5\times10^6$ orbits). We find that small spectral fractions ($\lesssim0.01$) are strongly associated with dynamical stability on long timescales ($5\times10^9$ orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many non-resonant multi-planet systems, and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long term (in)stability of multi-planet systems from short simulations.

Orbital resonances play an important role in the dynamics of planetary systems. Classical theoretical analyses found in textbooks report that libration widths of first order mean motion resonances diverge for nearly circular orbits. Here we examine the nature of this divergence with a non-perturbative analysis of a few first order resonances interior to a Jupiter-mass planet. We show that a first order resonance has two branches, the pericentric and the apocentric resonance zone. As the eccentricity approaches zero, the centers of these zones diverge away from the nominal resonance location but their widths shrink. We also report a novel finding of "bridges" between adjacent first order resonances: at low eccentricities, the apocentric libration zone of a first order resonance smoothly connects with the pericentric libration zone of the neighboring first order resonance. These bridges may facilitate resonant migration across large radial distances in planetary systems, entirely in the low eccentricity regime.

Luna 3 (or Lunik 3 in Russian sources) was the first spacecraft to perform a flyby of the Moon. Launched in October 1959 on a translunar trajectory with large semi-major axis and eccentricity, it collided with the Earth in late March 1960. The short, 6-month dynamical lifetime has often been explained through an increase in eccentricity due to the Lidov-Kozai effect. However, the classical Lidov-Kozai solution is only valid in the limit of small semi-major axis ratio, a condition that is satisfied only for solar (but not for lunar) perturbations. We undertook a study of the dynamics of Luna 3 with the aim of assessing the principal mechanisms affecting its evolution. We analyze the Luna 3 trajectory by generating accurate osculating solutions, and by comparing them to integrations of singly- and doubly-averaged equations of motion in vectorial form. Lunar close encounters, which cannot be reproduced in an averaging approach, decisively affect the trajectory and break the doubly-averaged dynamics. Solar perturbations induce oscillations of intermediate period that affect the geometry of the close encounters and cause the singly-averaged and osculating inclinations to change quadrants (the orbital plane ``flips''). We find that the peculiar evolution of Luna 3 can only be explained by taking into account lunar close encounters and intermediate-period terms; such terms are averaged out in the Lidov-Kozai solution, which is not adequate to describe translunar or cislunar trajectories. Understanding the limits of the Lidov-Kozai solution is of particular significance for the motion of objects in the Earth-Moon environment and of exoplanetary systems.

We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first establish a baseline for evaluating various techniques for close-encounter prediction and collision-probability calculation (Hoots et al. 1984, Gronchi 2005, JeongAhn and Malhotra 2015) by carrying out brute-force numerical simulations and using a sequence of filters to greatly reduce the computational expense of the algorithm. Next, we estimate conjunction events in the orbital environment following the anticipated deployments of the OneWeb LEO and SpaceX Starlink mega-constellations. As a final step, we investigate Minimum Space Occupancy (MiSO) orbits (Bombardelli et al. 2018), a generalization of the well-known frozen orbits that account for the perturbed-Keplerian dynamics of the Earth-Moon-Sun-satellite system. We evaluate the ability of MiSO configurations of the proposed mega-constellations, as suggested by Bombardelli et al. 2018, to reduce the risk of endogenous (intra-constellation) collisions. The results indicate that the adoption of the MiSO orbital configuration can significantly reduce risk with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites.

Most of the major planets in the Solar System support populations of co-orbiting bodies, known as Trojans, at their L4 and L5 Lagrange points. In contrast, Earth has only one known co-orbiting companion. This paper presents the results from a search for Earth Trojans using the DECam instrument on the Blanco Telescope at CTIO. This search found no additional Trojans in spite of greater coverage compared to previous surveys of the L5 point. Therefore, the main result of this work is to place the most stringent constraints to date on the population of Earth Trojans. These constraints depend on assumptions regarding the underlying population properties, especially the slope of the magnitude distribution (which in turn depends on the size and albedo distributions of the objects). For standard assumptions, we calculate upper limits to a 90% confidence limit on the L5 population of $N_{ET}<1$ for magnitude $H<15.5$, $N_{ET}=60-85$ for $H<19.7$, and $N_{ET}\ $= 97 for $H=20.4$. This latter magnitude limit corresponds to Trojans $\sim$300 m in size for albedo $0.15$. At H=19.7, these upper limits are consistent with previous L4 Earth Trojan constraints and significantly improve L5 constraints.

We utilize ALMA archival data to estimate the dust disk size of 152 protoplanetary disks in Lupus (1-3 Myr), Chamaeleon I (2-3 Myr), and Upper-Sco (5-11 Myr). We combine our sample with 47 disks from Tau/Aur and Oph whose dust disk radii were estimated, as here, through fitting radial profile models to visibility data. We use these 199 homogeneously derived disk sizes to identify empirical disk-disk and disk-host property relations as well as to search for evolutionary trends. In agreement with previous studies, we find that dust disk sizes and millimeter luminosities are correlated, but show for the first time that the relationship is not universal between regions. We find that disks in the 2-3 Myr-old Cha I are not smaller than disks in other regions of similar age, and confirm the Barenfeld et al. (2017) finding that the 5-10 Myr USco disks are smaller than disks belonging to younger regions. Finally, we find that the outer edge of the Solar System, as defined by the Kuiper Belt, is consistent with a population of dust disk sizes which have not experienced significant truncation.

Our understanding of the history of the solar system has undergone a revolution in recent years, owing to new theoretical insights into the origin of Pluto and the discovery of the Kuiper belt and its rich dynamical structure. The emerging picture of dramatic orbital migration of the planets driven by interaction with the primordial Kuiper belt is thought to have produced the final solar system architecture that we live in today. This paper gives a brief summary of this new view of our solar system's history, and reviews the astronomical evidence in the resonant populations of the Kuiper belt.

The Scattered Disk Objects (SDOs) are thought to be a small fraction of the ancient population of leftover planetesimals in the outer solar system that were gravitationally scattered by the giant planets and have managed to survive primarily by capture and sticking in Neptune's exterior mean motion resonances (MMRs). In order to advance understanding of the role of MMRs in the dynamics of the SDOs, we investigate the phase space structure of a large number of Neptune's MMRs in the semi-major axis range 33--140~au by use of Poincar\'e sections of the circular planar restricted three body model for the full range of particle eccentricity pertinent to SDOs. We find that, for eccentricities corresponding to perihelion distances near Neptune's orbit, distant MMRs have stable regions with widths that are surprisingly large and of similar size to those of the closer-in MMRs. We identify a phase-shifted second resonance zone that exists in the phase space at planet-crossing eccentricities but not at lower eccentricities; this second resonance zone plays an important role in the dynamics of SDOs in lengthening their dynamical lifetimes. Our non-perturbative measurements of the sizes of the stable resonance zones confirm previous results and provide an additional explanation for the prominence of the $N$:1 sequence of MMRs over the $N$:2, $N$:3 sequences and other MMRs in the population statistics of SDOs; our results also provide a tool to more easily identify resonant objects.

We present numerical simulations of giant planet migration in our solar system and examine how the speed of planetary migration affects inclinations in the resulting population of small bodies (test particles) scattered outward and subsequently captured into Neptune's 3:2 mean motion resonance (the Plutinos) as well as the hot classical Kuiper belt population. We do not find a consistent relationship between the degree of test particle inclination excitation and e-folding planet migration timescales in the range 5-50 Myr. Our results present a counter-example to Nesvorny 2015's finding that the Plutino and hot classical inclinations showed a marked increase with increasing e-folding timescales for Neptune's migration. We argue that these differing results are likely due to differing secular architectures of the giant planets during and after migration. Small changes in the planets' initial conditions and differences in the numerical implementation of planet migration can result in different amplitudes of the planets' inclination secular modes, and this can lead to different final inclination distributions for test particles in the simulations. We conclude that the observed large inclination dispersion of Kuiper belt objects does not require Neptune's migration to be slow; planetary migration with e-folding timescales of 5, 10, 30, and 50 Myr can all yield inclination dispersions similar to the observed Plutino and hot classical populations, with no correlation between the degree of inclination excitation and migration speed.

The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-size bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system argue for two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation.

The existence of Earth's Trojan asteroids is not well constrained and represents a major gap in our inventory of small bodies in near-Earth space. Their discovery would be of high scientific and human interest.

We propose a LSST Solar System near-Sun Survey, to be implemented during twilight hours, that extends the seasonal reach of LSST to its maximum as fresh sky is uncovered at about 50 square degrees per night (1500 sq. deg. per lunation) in the morning eastern sky, and surveyable sky is lost at the same rate to the western evening sky due to the Earth's synodic motion. By establishing near-horizon fence post picket lines to the far west and far east we address Solar System science use cases (including Near Earth Objects, Interior Earth Objects, Potentially Hazardous Asteroids, Earth Trojans, near-Sun asteroids, sun-grazing comets, and dormant comets) as well as provide the first look and last look that LSST will have at the transient and variable objects within each survey field. This proposed near-Sun Survey will also maximize the overlap with the field of regard of the proposed NEOCam spacecraft that will be stationed at the Earth's L1 Lagrange point and survey near quadrature with the Sun. This will allow LSST to incidently follow-up NEOCam targets and vice-versa (as well as targets from missions such as Euclid), and will roughly correspond to the Earth's L4 and L5 regions.

Two planetary mass objects in the far outer Solar System --- collectively referred to here as Planet X --- have recently been hypothesized to explain the orbital distribution of distant Kuiper Belt Objects. Neither planet is thought to be exceptionally faint, but the sky locations of these putative planets are poorly constrained. Therefore, a wide area survey is needed to detect these possible planets. The Large Synoptic Survey Telescope (LSST) will carry out an unbiased, large area (around 18,000 deg$^2$), deep (limiting magnitude of individual frames of 24.5) survey (the "wide-fast-deep" survey) of the southern sky beginning in 2022, and is therefore an important tool to search for these hypothesized planets. Here we explore the effectiveness of LSST as a search platform for these possible planets. Assuming the current baseline cadence (which includes the wide-fast-deep survey plus additional coverage) we estimate that LSST will confidently detect or rule out the existence of Planet X in 61\% of the entire sky. At orbital distances up to $\sim$75 au, Planet X could simply be found in the normal nightly moving object processing; at larger distances, it will require custom data processing. We also discuss the implications of a non-detection of Planet X in LSST data.

Observations of Kuiper belt objects (KBOs) in Neptune's 5:2 resonance present two puzzles: this third order resonance hosts a surprisingly large population, comparable to the prominent populations of Plutinos and Twotinos in the first order 3:2 and 2:1 resonances, respectively; secondly, their eccentricities are concentrated near $0.4$. To shed light on these puzzles, we investigate the phase space near this resonance with use of Poincar\'e sections of the circular planar restricted three body model. We find several transitions in the phase space structure with increasing eccentricity, which we explain with the properties of the resonant orbit relative to Neptune's. The resonance width is narrow for very small eccentricities, but widens dramatically for $e\gtrsim0.2$, reaching a maximum near $e\approx0.4$, where it is similar to the maximum widths of the 2:1 and 3:2 resonances. We confirm these results with N-body numerical simulations, including the effects of all four giant planets and a wide range of orbital inclinations of the KBOs. We find that the boundaries of the stable resonance zone are not strongly sensitive to inclination and remain very similar to those found with the simplified three body model, with the caveat that orbits of eccentricity above $\sim0.53$ are unstable; higher eccentricity orbits are phase-protected from destabilizing encounters with Neptune but not with Uranus. These results show that the 5:2 resonant KBOs are not more puzzling than the Plutinos and Twotinos; however, detailed understanding of the origins of eccentric, inclined resonant KBOs remains a challenge.

Dwarf planet Ceres, the largest object in the Main Asteroid Belt, has a surface that exhibits a range of crater densities for a crater diameter range of 5-300 km. In all areas the shape of the craters' size-frequency distribution is very similar to those of the most ancient heavily cratered surfaces on the terrestrial planets. The most heavily cratered terrain on Ceres covers ~15% of its surface and has a crater density similar to the highest crater density on <1% of the lunar highlands. This region of higher crater density on Ceres probably records the high impact rate at early times and indicates that the other 85% of Ceres was partly resurfaced after the Late Heavy Bombardment (LHB) at ~4 Ga. The Ceres cratering record strongly indicates that the period of Late Heavy Bombardment originated from an impactor population whose size-frequency distribution resembles that of the Main Belt Asteroids.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

We measure the mid-plane of the main asteroid belt by using the observational data of a nearly complete and unbiased sample of asteroids, and find that it has inclination $\bar{I}=0.93\pm0.04$ degrees and longitude of ascending node $\bar{\Omega}=87.6\pm2.6$ degrees (in J2000 ecliptic-equinox coordinate system). This plane differs significantly from previously published measurements, and it is also distinctly different than the solar system's invariable plane as well as Jupiter's orbit plane. The mid-plane of the asteroid belt is theoretically expected to be a slightly warped sheet whose local normal is controlled by the gravity of the major planets. Specifically, its inclination and longitude of ascending node varies with semi-major axis and time (on secular timescales), and is defined by the forced solution of secular perturbation theory; the $\nu_{16}$ nodal secular resonance is predicted to cause a significant warp of the mid-plane in the inner asteroid belt. We test the secular theory by measuring the current location of the asteroids' mid-plane in finer semi-major axis bins. We find that the measured mid-plane in the middle and outer asteroid belt is consistent, within 3--$\sigma$ confidence level, with the prediction of secular perturbation theory, but a notable discrepancy is present in the inner asteroid belt near $\sim2$ AU.

Several anomalous features in the orbital distribution of distant Kuiper Belt objects have recently sparked suggestions of an unseen large planet (or two) orbiting beyond Neptune. I review the theoretical arguments and the prospects for the existence of unseen massive bodies in the distant solar system.

We measured the mean plane of the Kuiper belt as a function of semi-major axis. For the classical Kuiper belt as a whole (the non-resonant objects in the semi-major axis range 42--48~au), we find a mean plane of inclination $i_m=1.8^{\circ}$$^{+0.7^{\circ}}_{-0.4^{\circ}}$ and longitude of ascending node $\Omega_m=77^{\circ}$$^{+18^{\circ}}_{-14^{\circ}}$ (in the J2000 ecliptic-equinox coordinate system), in accord with theoretical expectations of the secular effects of the known planets. With finer semi-major axis bins, we detect a statistically significant warp in the mean plane near semi-major axes 40--42~au. Linear secular theory predicts a warp near this location due to the $\nu_{18}$ nodal secular resonance, however the measured mean plane for the 40.3-42~au semi-major axis bin (just outside the $\nu_{18}$) is inclined $\sim13^{\circ}$ to the predicted plane, a nearly 3-$\sigma$ discrepancy. For the more distant Kuiper belt objects of semi-major axes in the range 50--80~au, the expected mean plane is close to the invariable plane of the solar system, but the measured mean plane deviates greatly from this: it has inclination $i_m=9.1^{\circ}$$^{+6.6^{\circ}}_{-3.8^{\circ}}$ and longitude of ascending node $\Omega_m=227^{\circ}$$^{+18^{\circ}}_{-44^{\circ}}$. We estimate this deviation from the expected mean plane to be statistically significant at the $\sim97-99\%$ confidence level. We discuss several possible explanations for this deviation, including the possibility that a relatively close-in ($a\lesssim100$~au), unseen small planetary-mass object in the outer solar system is responsible for the warping.

Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay close to their parent bodies. Here we present ALMA 1.3 mm observations of HD 95086, a young early-type star that hosts a directly imaged giant planet b and a massive debris disk with both asteroid- and Kuiper-belt analogs. The location of the Kuiper-belt analog is resolved for the first time. The system can be depicted as a broad ($\Delta R/R \sim$0.84), inclined (30\arcdeg$\pm$3\arcdeg) ring with millimeter emission peaked at 200$\pm$6 au from the star. The 1.3 mm disk emission is consistent with a broad disk with sharp boundaries from 106$\pm$6 to 320$\pm$20 au with a surface density distribution described by a power law with an index of --0.5$\pm$0.2. Our deep ALMA map also reveals a bright source located near the edge of the ring, whose brightness at 1.3 mm and potential spectral energy distribution are consistent with it being a luminous star-forming galaxy at high redshift. We set constraints on the orbital properties of planet b assuming co-planarity with the observed disk.

Many topics in planetary studies demand an estimate of the collision probability of two objects moving on nearly Keplerian orbits. In the classic works of \"Opik (1951) and Wetherill (1967), the collision probability was derived by linearizing the motion near the collision points, and there is now a vast literature using their method. We present here a simpler and more physically motivated derivation for non-tangential collisions in Keplerian orbits, as well as for tangential collisions that were not previously considered. Our formulas have the added advantage of being manifestly symmetric in the parameters of the two colliding bodies. In common with the \"Opik-Wetherill treatments, we linearize the motion of the bodies in the vicinity of the point of orbit intersection (or near the points of minimum distance between the two orbits) and assume a uniform distribution of impact parameter within the collision radius. We point out that the linear approximation leads to singular results for the case of tangential encounters. We regularize this singularity by use of a parabolic approximation of the motion in the vicinity of a tangential encounter.

Mean motion resonances [MMRs] play an important role in the formation and evolution of planetary systems and have significantly influenced the orbital properties and distribution of planets and minor planets in the solar system as well as exo-planetary systems. Most previous theoretical analyses have focused on the low-to-moderate eccentricity regime, but with new discoveries of high eccentricity resonant minor planets and even exoplanets, there is increasing motivation to examine MMRs in the high eccentricity regime. Here we report on a study of the high eccentricity regime of MMRs in the circular planar restricted three-body problem. Non-perturbative numerical analyses of the 2:1 and the 3:2 interior resonances are carried out for a wide range of secondary-to-primary mass ratio, and for a wide range of eccentricity of the test particle. The surface-of-section technique is used to study the phase space structure near resonances. We identify transitions in phase space at certain critical eccentricities related to the geometry of resonant orbits; new stable libration zones appear at high eccentricity at libration centers shifted from those at low eccentricities. We present novel results on the mass and eccentricity dependence of the resonance libration centers and their widths in semi-major axis. Our results show that MMRs have sizable libration zones at high eccentricities, comparable to those at lower eccentricities.

The observationally complete sample of the main belt asteroids now spans more than two orders of magnitude in size and numbers more than 64,000 (excluding collisional family members). We undertook an analysis of asteroids' eccentricities and their interpretation with simple physical models. We find that Plummer's (1916) conclusion that the asteroids' eccentricities follow a Rayleigh distribution holds for the osculating eccentricities of large asteroids, but the proper eccentricities deviate from a Rayleigh distribution: there is a deficit of eccentricities smaller than $\sim0.1$ and an excess of larger eccentricities. We further find that the proper eccentricities do not depend significantly on asteroid size but have strong dependence on heliocentric distance: the outer asteroid belt follows a Rayleigh distribution, but the inner belt is strikingly different. Eccentricities in the inner belt can be modeled as a vector sum of a primordial eccentricity vector of random orientation and magnitude drawn from a Rayleigh distribution of parameter $\sim0.06$, and an excitation of random phase and magnitude $\sim0.13$. These results imply that a late dynamical excitation of the asteroids occurred, it was independent of asteroid size, it was stronger in the inner belt than in the outer belt. We discuss implications for the primordial asteroid belt and suggest that the observationally complete sample size of main belt asteroids is large enough that more sophisticated model-fitting of the eccentricities is warranted and could serve to test alternative theoretical models of the dynamical excitation history of asteroids and its links to the migration history of the giant planets.

The four longest period Kuiper belt objects have orbital periods close to integer ratios with each other. A hypothetical planet with orbital period $\sim$17,117 years, semimajor axis $\sim$665 AU, would have N/1 and N/2 period ratios with these four objects. The orbital geometries and dynamics of resonant orbits constrain the orbital plane, the orbital eccentricity and the mass of such a planet, as well as its current location in its orbital path.

The Karin family is a very young asteroid family created by an asteroid breakup 5.8 Myr ago. Since the members of this family probably have not experienced significant orbital or collisional evolution yet, it is possible that they still preserve properties of the original family-forming event in terms of their spin state. As we carried out a series of photometric observations of the Karin family asteroids, here we report an analysis result of lightcurves including the rotation period of eleven members. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs or smaller MBAs (diameter D<12 km), and even lower than that of larger MBAs (D>130 km). We investigated a correlation between the peak-to-peak variation magnitude reduced to zero solar phase angle and the rotation period of the eleven Karin family asteroids, and found a possible trend that elongated members have lower spin rate, and less elongated members have higher spin rate. However, this has to be confirmed by another series of future observations.

The distribution of orbital period ratios of adjacent planets in extra-solar planetary systems discovered by the {\it Kepler} space telescope exhibits a peak near $\sim1.5$--$2$, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below $\sim1.5$. We find from this data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill's criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed $10^{-3.2}m_*$, or about two-thirds Jupiter mass for solar mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill's criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of $\log m/M_\oplus$ of $\sim(0.6-1.0)$ and $\sim(1.1-1.2)$, respectively.

We calculate the present-day impact flux on Mars and its variation over the Martian year, using the current data on the orbital distribution of known Mars-crossing minor planets. We adapt the {\"O}pik-Wetherill formulation for calculating collision probabilities, paying careful attention to the non-uniform distribution of the perihelion longitude and the argument of perihelion owed to secular planetary perturbations. We find that these previously neglected non-uniformities have a significant effect on the mean annual impact flux as well as its seasonal variation. The impact flux peaks when Mars is at aphelion, but the near-alignment of Mars' eccentricity vector with the mean direction of the eccentricity vectors of Mars-crossers causes the mean annual impact flux as well as the amplitude of the seasonal variation to be significantly lower than the estimate based on a uniform random distribution of perihelion longitudes of Mars-crossers. We estimate that the flux of large impactors (of absolute magnitude $H<16$) within $\pm30^\circ$ of Mars' aphelion is about four times larger than when the planet is near perihelion. Extrapolation of our results to a model population of meter-size Mars-crossers shows that if these small impactors have a uniform distribution of their angular elements, then their aphelion-to-perihelion impact flux ratio would be as large as 25. These theoretical predictions can be tested with observational data of contemporary impacts that is becoming available from spacecraft currently in orbit about Mars.

HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing 5% of the total dust mass. We also present detailed analysis of the disk SED and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (175 K) belt, a cold (55 K) disk, and an extended disk halo (up to 800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass 5 MJ, could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system.

We investigate the orbital evolution of particles in a planet's chaotic zone to determine their final destinations and their timescales of clearing. There are four possible final states of chaotic particles: collision with the planet, collision with the star, escape, or bounded but non-collision orbits. In our investigations, within the framework of the planar circular restricted three body problem for planet-star mass ratio $\mu$ in the range $10^{-9}$ to $10^{-1.5}$, we find no particles hitting the star. The relative frequencies of escape and collision with the planet are not scale-free, as they depend upon the size of the planet. For planet radius $R_p\ge0.001R_H$ where $R_H$ is the planet's Hill radius, we find that most chaotic zone particles collide with the planet for $\mu\lesssim10^{-5}$; particle scattering to large distances is significant only for higher mass planets. For fixed ratio $R_p/R_H$, the particle clearing timescale, $T_{cl}$, has a broken power-law dependence on $\mu$. A shallower power-law, $T_{cl}\sim \mu^{-{1/3}}$, prevails at small $\mu$ where particles are cleared primarily by collisions with the planet; a steeper power law, $T_{cl}\sim\mu^{-{3/2}}$, prevails at larger $\mu$ where scattering dominates the particle loss. In the limit of vanishing planet radius, we find $T_{cl}\approx0.024\mu^{-{3\over2}}$. The interior and exterior boundaries of the annular zone in which chaotic particles are cleared are increasingly asymmetric about the planet's orbit for larger planet masses; the inner boundary coincides well with the classical first order resonance overlap zone, $\Delta a_{cl,int}\simeq1.2\mu^{0.28}a_p$; the outer boundary is better described by $\Delta a_{cl,ext}\simeq1.7\mu^{0.31}a_p$, where $a_p$ is the planet-star separation.

We review previously published and newly obtained crater size-frequency distributions in the inner solar system. These data indicate that the Moon and the terrestrial planets have been bombarded by two populations of objects. Population 1, dominating at early times, had nearly the same size distribution as the present-day asteroid belt, and produced the heavily cratered surfaces with a complex, multi-sloped crater size-frequency distribution. Population 2, dominating since about 3.8-3.7 Ga, has the same size distribution as near-Earth objects (NEOs), had a much lower impact flux, and produced a crater size distribution characterized by a differential -3 single-slope power law in the crater diameter range 0.02 km to 100 km. Taken together with the results from a large body of work on age-dating of lunar and meteorite samples and theoretical work in solar system dynamics, a plausible interpretation of these data is as follows. The NEO population is the source of Population 2 and it has been in near-steady state over the past ~3.7-3.8 gigayears; these objects are derived from the main asteroid belt by size-dependent non-gravitational effects that favor the ejection of smaller asteroids. However, Population 1 were main belt asteroids ejected from their source region in a size-independent manner, possibly by means of gravitational resonance sweeping during giant planet orbit migration; this caused the so-called Late Heavy Bombardment (LHB). The LHB began some time before ~3.9 Ga, peaked and declined rapidly over the next ~100 to 300 megayears, and possibly more slowly from about 3.8-3.7 Ga to ~2 Ga. A third crater population (Population S) consists of secondary impact craters that can dominate the cratering record at small diameters.

Planets can affect debris disk structure by creating gaps, sharp edges, warps, and other potentially observable signatures. However, there is currently no simple way for observers to deduce a disk-shepherding planet's properties from the observed features of the disk. Here we present a single equation that relates a shepherding planet's maximum mass to the debris ring's observed width in scattered light, along with a procedure to estimate the planet's eccentricity and minimum semimajor axis. We accomplish this by performing dynamical N-body simulations of model systems containing a star, a single planet, and a disk of parent bodies and dust grains to determine the resulting debris disk properties over a wide range of input parameters. We find that the relationship between planet mass and debris disk width is linear, with increasing planet mass producing broader debris rings. We apply our methods to five imaged debris rings to constrain the putative planet masses and orbits in each system. Observers can use our empirically-derived equation as a guide for future direct imaging searches for planets in debris disk systems. In the fortuitous case of an imaged planet orbiting interior to an imaged disk, the planet's maximum mass can be estimated independent of atmospheric models.

We examine the angular distributions of near-Earth objects (NEOs) which are often regarded as uniform. The apparent distribution of the longitude of ascending node, $\Omega$, is strongly affected by well-known seasonal effects in the discovery rate of NEOs. The deviation from the expected $\pi$-periodicity in the apparent distribution of $\Omega$ indicates that its intrinsic distribution is slightly enhanced along a mean direction, $\bar\Omega=111^\circ$; approximately 53% of NEOs have $\Omega$ values within $\pm90^\circ$ of $\bar\Omega$. We also find that each sub-group of NEOs (Amors, Apollos and Atens) has different observational selection effects which cause different non-uniformities in the apparent distributions of their arguments of perihelion $\omega$, and longitudes of perihelion $\varpi$. For their intrinsic distributions, our analysis reveals that the Apollo asteroids have non-uniform $\omega$ due to secular dynamics associated with inclination-eccentricity-$\omega$ coupling, and the Amors' $\varpi$ distribution is peaked towards the secularly forced eccentricity vector. The Apollos' $\omega$ distribution is axial, favoring values near 0 and $180^\circ$; the two quadrants centered at 0 and $180^\circ$ account for 55% of the Apollos' $\omega$ values. The Amors' $\varpi$ distribution peaks near $\bar\varpi=4^\circ$; sixty-one percent of Amors have $\varpi$ within $\pm90^\circ$ of this peak. We show that these modest but statistically significant deviations from uniform random distributions of angular elements are owed to planetary perturbations, primarily Jupiter's. It is remarkable that this strongly chaotic population of minor planets reveals the presence of Jupiter in its angular distributions.

The Centaurs are a population of small, planet-crossing objects in the outer solar system. They are dynamically short-lived and represent the transition population between the Kuiper belt and the Jupiter family short-period comets. Dynamical models and observations of the physical properties of the Centaurs indicate that they may have multiple source populations in the trans-Neptunian region. It has been suggested that the inclination distribution of the Centaurs may be useful in distinguishing amongst these source regions. The Centaurs, however, undergo many close encounters with the giant planets during their orbital evolution; here we show that these encounters can substantially determine the inclination distribution of the Centaurs. Almost any plausible initial inclination distribution of a Kuiper belt source results in Centaurs having inclinations peaked near 10-20 degrees. Our studies also find that the Kuiper belt is an extremely unlikely source of the retrograde Centaur that has been observed.

Vega and Fomalhaut, are similar in terms of mass, ages, and global debris disk properties; therefore, they are often referred as "debris disk twins". We present Spitzer 10-35 um spectroscopic data centered at both stars, and identify warm, unresolved excess emission in the close vicinity of Vega for the first time. The properties of the warm excess in Vega are further characterized with ancillary photometry in the mid infrared and resolved images in the far-infrared and submillimeter wavelengths. The Vega warm excess shares many similar properties with the one found around Fomalhaut. The emission shortward of ~30 um from both warm components is well described as a blackbody emission of ~170 K. Interestingly, two other systems, eps Eri and HR 8799, also show such an unresolved warm dust using the same approach. These warm components may be analogous to the solar system's zodiacal dust cloud, but of far greater. The dust temperature and tentative detections in the submillimeter suggest the warm excess arises from dust associated with a planetesimal ring located near the water-frost line and presumably created by processes occurring at similar locations in other debris systems as well. We also review the properties of the 2 um hot excess around Vega and Fomalhaut, showing that the dust responsible for the hot excess is not spatially associated with the dust we detected in the warm belt. We suggest it may arise from hot nano grains trapped in the magnetic field of the star. Finally, the separation between the warm and cold belt is rather large with an orbital ratio >~10 in all four systems. In light of the current upper limits on the masses of planetary objects and the large gap, we discuss the possible implications for their underlying planetary architecture, and suggest that multiple, low-mass planets likely reside between the two belts in Vega and Fomalhaut.

The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the restricted three-body problem. We show that for the 3:2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100M_{\oplus}. This mass range is higher than expected, by at least a factor of three, from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family's age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50-400 m/s, approximately 20-45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m/s of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ~50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1 sigma level, however both models predict more objects at larger relative velocities (>150 m/s) than have been identified.

We examine a low energy mechanism for the transfer of meteoroids between two planetary systems embedded in a star cluster using quasi-parabolic orbits of minimal energy. Using Monte Carlo simulations, we find that the exchange of meteoroids could have been significantly more efficient than previously estimated. Our study is relevant to astrobiology as it addresses whether life on Earth could have been transferred to other planetary systems in the solar system's birth cluster and whether life on Earth could have been transferred here from beyond the solar system. In the solar system, the timescale over which solid material was delivered to the region from where it could be transferred via this mechanism likely extended to several hundred million years (as indicated by the 3.8-4.0 Ga epoch of the Late Heavy Bombardment). This timescale could have overlapped with the lifetime of the Solar birth cluster (~100-500 Myr). Therefore, we conclude that lithopanspermia is an open possibility if life had an early start. Adopting parameters from the minimum mass solar nebula, considering a range of planetesimal size distributions derived from observations of asteroids and Kuiper Belt Objects and theoretical coagulation models, and taking into account Oort Cloud formation models, the expected number of bodies with mass > 10 kg that could have been transferred between the Sun and its nearest cluster neighbor could be of the order of 1E14-3E16, with transfer timescales of 10s Myr. We estimate that of the order of 3E8 x l(km) could potentially be life-bearing, where l(km) is the depth of the Earth crust in km that was ejected as the result of the early bombardment.

Kozai-Lidov oscillations of Jupiter-mass planets, excited by comparable planetary or brown dwarf mass perturbers were recently shown in numerical experiments to be slowly modulated and to exhibit striking features, including extremely high eccentricities and the generation of retrograde orbits with respect to the perturber. Here we solve this problem analytically for the case of a test particle orbiting a host star and perturbed by a distant companion whose orbit is eccentric and highly inclined. We give analytic expressions for the conditions that produce retrograde orbits and high eccentricities. This mechanism likely operates in various systems thought to involve Kozai-Lidov oscillations such as tight binaries, mergers of compact objects, irregular moons of planets and many others. In particular, it could be responsible for exciting eccentricities and inclinations of exo-planetary orbits and be important for understanding the spin-orbit (mis)alignment of hot Jupiters.

We investigate the long-term evolution of the inclinations of the known classical and resonant Kuiper belt objects (KBOs). This is partially motivated by the observed bimodal inclination distribution and by the putative physical differences between the low- and high-inclination populations. We find that some classical KBOs undergo large changes in inclination over gigayear timescales, which means that a current member of the low-inclination population may have been in the high-inclination population in the past, and vice versa. The dynamical mechanisms responsible for the time-variability of inclinations are predominantly distant encounters with Neptune and chaotic diffusion near the boundaries of mean motion resonances. We reassess the correlations between inclination and physical properties including inclination time-variability. We find that the size-inclination and color-inclination correlations are less statistically significant than previously reported (mostly due to the increased size of the data set since previous works with some contribution from inclination variability). The time-variability of inclinations does not change the previous finding that binary classical KBOs have lower inclinations than non-binary objects. Our study of resonant objects in the classical Kuiper belt region includes objects in the 3:2, 7:4, 2:1, and eight higher-order mean motion resonances. We find that these objects (some of which were previously classified as non-resonant) undergo larger changes in inclination compared to the non-resonant population, indicating that their current inclinations are not generally representative of their original inclinations. They are also less stable on gigayear timescales.

We calculate the eccentricity excitation of asteroids produced by the sweeping $\nu_6$ secular resonance during the epoch of planetesimal-driven giant planet migration in the early history of the solar system. We derive analytical expressions for the magnitude of the eccentricity change and its dependence on the sweep rate and on planetary parameters; the $\nu_6$ sweeping leads to either an increase or a decrease of eccentricity depending on an asteroid's initial orbit. Based on the slowest rate of $\nu_6$ sweeping that allows a remnant asteroid belt to survive, we derive a lower limit on Saturn's migration speed of $\sim0.15\AU\My^{-1}$ during the era that the $\nu_6$ resonance swept through the inner asteroid belt (semimajor axis range 2.1--$2.8\AU$). This rate limit is for Saturn's current eccentricity, and scales with the square of Saturn's eccentricity; the limit on Saturn's migration rate could be lower if Saturn's eccentricity were lower during its migration. Applied to an ensemble of fictitious asteroids, our calculations show that a prior single-peaked distribution of asteroid eccentricities would be transformed into a double-peaked distribution due to the sweeping of the $\nu_6$. Examination of the orbital data of main belt asteroids reveals that the proper eccentricities of the known bright ($H \leq10.8$) asteroids may be consistent with a double-peaked distribution. If so, our theoretical analysis then yields two possible solutions for the migration rate of Saturn and for the dynamical states of the pre-migration asteroid belt: a dynamically cold state (single-peaked eccentricity distribution with mean of $\sim0.05$) linked with Saturn's migration speed $\sim 4\AU\My^{-1}$, or a dynamically hot state (single-peaked eccentricity distribution with mean of $\sim0.3$) linked with Saturn's migration speed $\sim 0.8\AU\My^{-1}$.

In addition to the Sun, six other stars are known to harbor multiple planets and debris disks: HD 69830, HD 38529, HD 128311, HD 202206, HD 82943 and HR 8799. In this paper we set constraints on the location of the dust-producing planetesimals around the latter four systems. We use a radiative transfer model to analyze the spectral energy distributions of the dust disks (including two new Spitzer IRS spectra presented in this paper), and a dynamical model to assess the long-term stability of the planetesimals' orbits. As members of a small group of stars that show evidence of harboring a multiple planets and planetesimals, their study can help us learn about the diversity of planetary systems.

Recent lunar crater studies have revealed an asymmetric distribution of rayed craters on the lunar surface. The asymmetry is related to the synchronous rotation of the Moon: there is a higher density of rayed craters on the leading hemisphere compared with the trailing hemisphere. Rayed craters represent generally the youngest impacts. The purpose of this paper is to test the hypotheses that (i) the population of Near-Earth asteroids (NEAs) is the source of the impactors that have made the rayed craters, and (ii) that impacts by this projectile population account quantitatively for the observed asymmetry. We carried out numerical simulations of the orbital evolution of a large number of test particles representing NEAs in order to determine directly their impact flux on the Moon. The simulations were done in two stages. In the first stage we obtained encounter statistics of NEAs on the Earth's activity sphere. In the second stage we calculated the direct impact flux of the encountering particles on the surface of the Moon; the latter calculations were confined within the activity sphere of the Earth. A steady-state synthetic population of NEAs was generated from a debiased orbital distribution of the known NEAs. We find that the near-Earth asteroids do have an asymmetry in their impact flux on the Moon: apex-to-antapex ratio of 1.32 +/- 0.01. However, the observed rayed crater distribution's asymmetry is significantly more pronounced: apex-to-antapex ratio of 1.65 +/- 0.16. Our results suggest the existence of an undetected population of slower (low impact velocity) projectiles, such as a population of objects nearly coorbiting with Earth; more observational study of young lunar craters is needed to secure this conclusion.

The cumulative effects of weak resonant and secular perturbations by the major planets produce chaotic behavior of asteroids on long timescales. Dynamical chaos is the dominant loss mechanism for asteroids with diameters D > 10 km in the current asteroid belt. In a numerical analysis of the long term evolution of test particles in the main asteroid belt region, we find that the dynamical loss history of test particles from this region is well described with a logarithmic decay law. In our simulations the loss rate function that is established at t = 1 My persists with little deviation to at least t = 4 Gy. Our study indicates that the asteroid belt region has experienced a significant amount of depletion due to this dynamical erosion - having lost as much as ~50% of the large asteroids - since 1 My after the establishment of the current dynamical structure of the asteroid belt. Because the dynamical depletion of asteroids from the main belt is approximately logarithmic, an equal amount of depletion occurred in the time interval 10-200 My as in 0.2-4 Gy, roughly ~30% of the current number of large asteroids in the main belt over each interval. We find that asteroids escaping from the main belt due to dynamical chaos have an Earth impact probability of ~0.3%. Our model suggests that the rate of impacts from large asteroids has declined by a factor of 3 over the last 3 Gy, and that the present-day impact flux of D > 10 km objects on the terrestrial planets is roughly an order of magnitude less than estimates currently in use in crater chronologies and impact hazard risk assessments.

We have obtained a full suite of Spitzer observations to characterize the debris disk around HR 8799 and to explore how its properties are related to the recently discovered set of three massive planets orbiting the star. We distinguish three components to the debris system: (1) warm dust (T ~150 K) orbiting within the innermost planet; (2) a broad zone of cold dust (T ~45 K) with a sharp inner edge, orbiting just outside the outermost planet and presumably sculpted by it; and (3) a dramatic halo of small grains originating in the cold dust component. The high level of dynamical activity implied by this halo may arise due to enhanced gravitational stirring by the massive planets. The relatively young age of HR 8799 places it in an important early stage of development and may provide some help in understanding the interaction of planets and planetary debris, an important process in the evolution of our own solar system.

We carried out extensive numerical orbit integrations to probe the long-term chaotic dynamics of the two strongest mean motion resonances of Neptune in the Kuiper belt, the 3:2 (Plutinos) and 2:1 (Twotinos). Our primary results include a computation of the relative volumes of phase space characterized by large- and small-resonance libration amplitudes, and maps of resonance stability measured by mean chaotic diffusion rate. We find that Neptune's 2:1 resonance has weaker overall long-term stability than the 3:2 -- only ~15% of Twotinos are projected to survive for 4 Gyr, compared to ~27% of Plutinos, based on an extrapolation from our 1-Gyr integrations. We find that Pluto has only a modest effect, causing a ~4% decrease in the Plutino population that survives to 4 Gyr. Given current observational estimates, and assuming an initial distribution of particles proportional to the local phase space volume in the resonance, we conclude that the primordial populations of Plutinos and Twotinos formerly made up more than half the population of the classical and resonant Kuiper Belt. We also conclude that Twotinos were originally nearly as numerous as Plutinos; this is consistent with predictions from early models of smooth giant planet migration and resonance sweeping of the Kuiper Belt, and provides a useful constraint for more detailed models.

The Centaurs are a transient population of small bodies in the outer solar system whose orbits are strongly chaotic. These objects typically suffer significant changes of orbital parameters on timescales of a few thousand years, and their orbital evolution exhibits two types of behaviors described qualitatively as random-walk and resonance-sticking. We have analyzed the chaotic behavior of the known Centaurs. Our analysis has revealed that the two types of chaotic evolution are quantitatively distinguishable: (1) the random walk-type behavior is well described by so-called generalized diffusion in which the rms deviation of the semimajor axis grows with time t as ~t^H, with Hurst exponent H in the range 0.22--0.95, however (2) orbital evolution dominated by intermittent resonance sticking, with sudden jumps from one mean motion resonance to another, has poorly defined H. We further find that these two types of behavior are correlated with Centaur dynamical lifetime: most Centaurs whose dynamical lifetime is less than ~22 Myr exhibit generalized diffusion, whereas most Centaurs of longer dynamical lifetimes exhibit intermittent resonance sticking. We also find that Centaurs in the diffusing class are likely to evolve into Jupiter-family comets during their dynamical lifetimes, while those in the resonance-hopping class do not.

The main asteroid belt lies between the orbits of Mars and Jupiter, but the region is not uniformly filled with asteroids. There are gaps, known as the Kirkwood gaps, in the asteroid distribution in distinct locations that are associated with orbital resonances with the giant planets; asteroids placed in these locations follow chaotic orbits and escape from the asteroid belt. Here we show that the observed distribution of main belt asteroids does not fill uniformly even those regions that are dynamically stable over the age of the solar system. We find a pattern of excess depletion of asteroids, particularly just outward of the Kirkwood Gaps associated with the 5:2, the 7:3, and the 2:1 jovian resonances. These features are not accounted for by planetary perturbations in the current structure of the solar system, but are consistent with dynamical ejection of asteroids by the sweeping of gravitational resonances during the migration of Jupiter and Saturn ~4 gigayears ago.

The exchange of meteorites among the terrestrial planets of our Solar System is a well established phenomenon that has triggered discussion of lithopanspermia within the Solar System. Similarly, could solid material be transferred across planetary systems? To address this question, we explore the dynamics of the transfer of small bodies between planetary systems. In particular, we examine a dynamical process that yields very low escape velocities using nearly parabolic trajectories, and the reverse process that allows for low velocity capture. These processes are chaotic and provide a mechanism for minimal energy transfer that yield an increased transfer probability compared to that of previously studied mechanisms that have invoked hyperbolic trajectories. We estimate the transfer probability in a stellar cluster as a function of stellar mass and cluster size. We find that significant amounts of solid material could potentially have been transferred from the early Solar System to our nearest neighbor stars. While this low velocity mechanism improves the odds for interstellar lithopanspermia, the exchange of biologically active materials across stellar systems depends greatly upon the highly uncertain viability of organisms over the timescales for transfer, typically millions of years.

The short period Jupiter family comets (JFCs) are thought to originate in the Kuiper Belt; specifically, a dynamical subclass of the Kuiper Belt known as the `scattered disk' is argued to be the dominant source of JFCs. However, the best estimates from observational surveys indicate that this source may fall short by more than two orders of magnitude the estimates obtained from theoretical models of the dynamical evolution of Kuiper belt objects into JFCs. We re-examine the scattered disk as a source of the JFCs and make a rigorous estimate of the discrepancy. We find that the uncertainties in the dynamical models combined with a change in the size distribution function of the scattered disk at faint magnitudes (small sizes) beyond the current observational limit offer a possible but problematic resolution to the discrepancy. We discuss several other possibilities: that the present population of JFCs is a large fluctuation above their long term average, that larger scattered disk objects tidally break-up into multiple fragments during close planetary encounters as their orbits evolve from the trans-Neptune zone to near Jupiter, or that there are alternative source populations that contribute significantly to the JFCs. Well-characterized observational investigations of the Centaurs, objects that are transitioning between the trans-Neptune Kuiper belt region and the inner solar system, can test the predictions of the non-steady state and the tidal break-up hypotheses. The classical and resonant classes of the Kuiper belt are worth re-consideration as significant additional or alternate sources of the JFCs.

The extrasolar system OGLE-2006-BLG-109L is the first multiple-planet system to be discovered by gravitational microlensing (Gaudi et al., 2008); the two large planets that have been detected have mass ratios, semimajor axis ratios, and equilibrium temperatures that are similar to those of Jupiter and Saturn; the mass of the host star is only 0.5 M_sun, and the system is more compact than our own Solar system. We find that in the habitable zone of the host star, the two detected planets resonantly excite large orbital eccentricities on a putative earth-mass planet, driving such a planet out of the habitable zone. We show that an additional inner planet of ~>0.3M_earth at <~0.1 AU would suppress the eccentricity perturbation and greatly improve the prospects for habitability of the system. Thus, the planetary architecture of a potentially habitable OGLE-2006-BLG-109L planetary system -- with two ``terrestrial'' planets and two jovian planets -- could bear very close resemblance to our own Solar system.

HD 38529 is a post-main sequence G8III/IV star (3.5 Gyr old) with a planetary system consisting of at least two planets having Msin(i) of 0.8 MJup and 12.2 MJup, semimajor axes of 0.13 AU and 3.74 AU, and eccentricities of 0.25 and 0.35, respectively. Spitzer observations show that HD 38529 has an excess emission above the stellar photosphere, with a signal-to-noise ratio (S/N) at 70 micron of 4.7, a small excess at 33 micron (S/N=2.6) and no excess <30 micron. We discuss the distribution of the potential dust-producing planetesimals from the study of the dynamical perturbations of the two known planets, considering in particular the effect of secular resonances. We identify three dynamically stable niches at 0.4-0.8 AU, 20-50 AU and beyond 60 AU. We model the spectral energy distribution of HD 38529 to find out which of these niches show signs of harboring dust-producing plantesimals. The secular analysis, together with the SED modeling resuls, suggest that the planetesimals responsible for most of the dust emission are likely located within 20-50 AU, a configuration that resembles that of the Jovian planets + Kuiper Belt in our Solar System. Finally, we place upper limits (8E-6 lunar masses of 10 micron particles) to the amount of dust that could be located in the dynamically stable region that exists between the two planets (0.25--0.75 AU).

(Abridged) The discovery of extrasolar planetary systems revealed an unexpected diversity of planetary systems that has revolutionized planet formation theory. A strong program of theoretical research is essential to maximize both the discovery potential and the scientific returns of future observational programs, so as to achieve a deeper understanding of the formation and evolution of planetary systems. We outline three broad categories of theoretical research: detailed studies of specific planetary systems, testing planet formation models by comparing their predictions to the observed exoplanet population, and detailed modeling of specific physical processes. We describe how such theoretical research plays an important role in analyzing observations for a wide range detection methods and contributes to understanding the Earth's place in the universe and the potential for Earth-like life beyond our solar system. In this white paper, we suggest how to maximize the scientific return of future exoplanet observations. Our recommendations include a strong theory program, support for multiple observational programs that will study a diverse set of planets and stars, significant observing time devoted to follow-up observations, and healthy collaboration between observers and theorists.

The dust disks observed around mature stars are evidence that plantesimals are present in these systems on spatial scales that are similar to that of the asteroids and the KBOs in the Solar System. These dust disks (a.k.a. ``debris disks'') present a wide range of sizes, morphologies and properties. It is inferred that their dust mass declines with time as the dust-producing planetesimals get depleted, and that this decline can be punctuated by large spikes that are produced as a result of individual collisional events. The lack of solid state features indicate that, generally, the dust in these disks have sizes larger than approximately 10 microns, but exceptionally, strong silicate features in some disks suggest the presence of large quantities of small grains, thought to be the result of recent collisions. Spatially resolved observations of debris disks show a diversity of structural features, such as inner cavities, warps, offsets, brightness asymmetries, spirals, rings and clumps. There is growing evidence that, in some cases, these structures are the result of the dynamical perturbations of a massive planet. Our Solar System also harbors a debris disk and some of its properties resemble those of extra-solar debris disks. From the cratering record, we can infer that its dust mass has decayed with time, and that there was at least one major ``spike'' in the past during the Late Heavy Bombardment. This offers a unique opportunity to use extra-solar debris disks to shed some light in how the Solar System might have looked in the past. Similarly, our knowledge of the Solar System is influencing our understanding of the types of processes which might be at play in the extra-solar debris disks.

We provide an overview of the Spitzer Legacy Program ``Formation and Evolution of Planetary Systems" (FEPS) which was proposed in 2000, begun in 2001, and executed aboard the Spitzer Space Telescope between 2003 and 2006. This program exploits the sensitivity of Spitzer to carry out mid-infrared spectrophotometric observations of solar-type stars. With a sample of ~ 328 stars ranging in age from ~ 3 Myr to ~ 3 Gyr, we trace the evolution of circumstellar gas and dust from primordial planet-building stages in young circumstellar disks through to older collisionally generated debris disks. When completed, our program will help define the time scales over which terrestrial and gas giant planets are built, constrain the frequency of planetesimal collisions as a function of time, and establish the diversity of mature planetary architectures. In addition to the observational program, we have coordinated a concomitant theoretical effort aimed at understanding the dynamics of circumstellar dust with and without the effects of embedded planets, dust spectral energy distributions, and atomic and molecular gas line emission. Together with the observations, these efforts will provide astronomical context for understanding whether our Solar System - and its habitable planet - is a common or a rare circumstance. Additional information about the FEPS project can be found on the team website: feps.as.arizona.edu

A moderately massive early Sun has been proposed to resolve the so-called faint early Sun paradox. We calculate the time-evolution of the solar mass that would be required by this hypothesis, using a simple parametrized energy-balance model for Earth's climate. Our calculations show that the solar mass loss rate would need to have been 2-3 orders of magnitude higher than present for a time on the order of ~2 Gy. Such a mass loss history is significantly at variance (both in timescale and in the magnitude of the mass loss rates) with that inferred from astronomical observations of mass loss in younger solar analogues. While suggestive, the astronomical data cannot completely rule out the possibility that the Sun had the required mass loss history; therefore, we also examine the effects of the hypothetical historical solar mass loss on orbital dynamics in the solar system, with a view to identifying additional tests of the hypothesis. Planetary and satellite orbits provide a few tests, but these are weak or non-unique.

Using data from the Spitzer Space Telescope Legacy Science Program ``Formation and Evolution of Planetary Systems'' (FEPS), we have searched for debris disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to have one or more massive planets. Only one of the sources, HD 38529, has excess emission above the stellar photosphere; at 70 micron the signal-to-noise ratio in the excess is 4.7 while at wavelengths < 30 micron there is no evidence of excess. The remaining sources show no excesses at any Spitzer wavelengths. Applying survival tests to the FEPS sample and the results for the FGK survey published in Bryden et al. (2006), we do not find a significant correlation between the frequency and properties of debris disks and the presence of close-in planets. We discuss possible reasons for the lack of a correlation.

A large disruption in the main asteroid belt can cause a large flux, an "asteroid shower", on the terrestrial planets. We quantitatively examine the hypothesis that such an event was the cause of the lunar late heavy bombardment (LHB). We performed numerical integrations of about 20000 test particles starting in the vicinity of the nu6 secular resonance in the main asteroid belt. The purpose of these integrations is to calculate, for each of the terrestrial planets, the collision probability of asteroids coming from an asteroid break-up event in the inner part of the main belt. Compared with previous studies, we simulate nearly two orders of magnitude larger number of particles, and we include the orbital effects of the eight planets, Mercury to Neptune. We also examined in detail the orbital evolution of asteroid fragments once they enter the Earth's activity sphere, including the effect of the Earth-Moon orbit. We obtained the collision probability, the distributions of impact velocities, impact positions, and impact angles of asteroid fragments on the Moon. The collision probability with the Moon (~0.1%) suggests that a fairly large parent body, 1000-1500 km in diameter, is required if the LHB event is to be ascribed to a single asteroid disruption. An even larger parent body is required for less favorable initial conditions than we investigated here. We conclude that an asteroid disruption event is not a viable explanation for the LHB.

We have carried out a sensitive search for gas emission lines at infrared and millimeter wavelengths for a sample of 15 young sun-like stars selected from our dust disk survey with the Spitzer Space Telescope. We have used mid-infrared lines to trace the warm (300-100 K) gas in the inner disk and millimeter transitions of 12CO to probe the cold (~20 K) outer disk. We report no gas line detections from our sample. Line flux upper limits are first converted to warm and cold gas mass limits using simple approximations allowing a direct comparison with values from the literature. We also present results from more sophisticated models following Gorti and Hollenbach (2004) which confirm and extend our simple analysis. These models show that the SI line at 25.23 micron can set constraining limits on the gas surface density at the disk inner radius and traces disk regions up to a few AU. We find that none of the 15 systems have more than 0.04 MJ of gas within a few AU from the disk inner radius for disk radii from 1 AU up to ~40 AU. These gas mass upper limits even in the 8 systems younger than ~30 Myr suggest that most of the gas is dispersed early. The gas mass upper limits in the 10-40 AU region, that is mainly traced by our CO data, are <2 Mearth. If these systems are analogs of the Solar System, either they have already formed Uranus- and Neptune-like planets or they will not form them beyond 100 Myr. Finally, the gas surface density upper limits at 1 AU are smaller than 0.01% of the minimum mass solar nebula for most of the sources. If terrestrial planets form frequently and their orbits are circularized by gas, then circularization occurs early.

New insights into the history of the inner solar system are derived from the impact cratering record of the Moon, Mars, Venus and Mercury, and from the size distributions of asteroid populations. Old craters from a unique period of heavy bombardment that ended $\sim$3.8 billion years ago were made by asteroids that were dynamically ejected from the main asteroid belt, possibly due to the orbital migration of the giant planets. The impactors of the past $\sim$3.8 billion years have a size distribution quite different from the main belt asteroids, but very similar to the population of near-Earth asteroids.

Nbody simulations are used to examine the consequences of Neptune's outward migration into the Kuiper Belt, with the simulated endstates being compared rigorously and quantitatively to the observations. These simulations confirm the findings of Chiang et al. (2003), who showed that Neptune's migration into a previously stirred-up Kuiper Belt can account for the Kuiper Belt Objects (KBOs) known to librate at Neptune's 5:2 resonance. We also find that capture is possible at many other weak, high-order mean motion resonances, such as the 11:6, 13:7, 13:6, 9:4, 7:3, 12:5, 8:3, 3:1, 7:2, and the 4:1. The more distant of these resonances, such as the 9:4, 7:3, 5:2, and the 3:1, can also capture particles in stable, eccentric orbits beyond 50 AU, in the region of phase space conventionally known as the Scattered Disk. Indeed, 90% of the simulated particles that persist over the age of the Solar System in the so-called Scattered Disk zone never had a close encounter with Neptune, but instead were promoted into these eccentric orbits by Neptune's resonances during the migration epoch. This indicates that the observed Scattered Disk might not be so scattered. This model also produced only a handful of Centaurs, all of which originated at Neptune's mean motion resonances in the Kuiper Belt. We also report estimates of the abundances and masses of the Belt's various subpopulations (e.g., the resonant KBOs, the Main Belt, and the so-called Scattered Disk), and also provide upper limits on the abundance of Centaurs and Neptune's Trojans, as well as upper limits on the sizes and abundances of hypothetical KBOs that might inhabit the a>50 AU zone.

We report several results related to the dynamical evolution of dust produced in the Kuiper Belt (KB). We show that its particle size frequency distribution in space is greatly changed from the distribution at production, as a results of the combined effects of radiation forces and the perturbations of the planets. We estimate the contribution of KB dust to the zodiacal cloud by calculating the radial profile of its number density near the ecliptic. We also study the contribution of KB dust to the population of interplanetary dust particles (IDPs) collected at Earth, by calculating geocentric encounter velocities and capture rates. Our models show, in contrast with previous studies, that KB dust grains on Earth-crossing orbits have high eccentricities and inclinations and, therefore, their encounter velocities are similar to those of cometary grains and not to asteroidal grains. We estimate that at most 25% in number of captured IDPs have cometary or KB origin; the KB contribution may be as low as 1%-2%. We present the velocity field of KB dust throughout the solar system; this, together with the number density radial profile, is potentially useful for planning spacecraft missions to the outer solar system.

Main sequence stars are commonly surrounded by debris disks, composed of cold dust continuously replenished by a reservoir of undetected dust-producing planetesimals. In a planetary system with a belt of planetesimals (like the Solar System's Kuiper Belt) and one or more interior giant planets, the trapping of dust particles in the mean motion resonances with the planets can create structure in the dust disk, as the particles accumulate at certain semimajor axes. Sufficiently massive planets may also scatter and eject dust particles out of a planetary system, creating a dust depleted region inside the orbit of the planet. In anticipation of future observations of spatially unresolved debris disks with the Spitzer Space Telescope, we are interested in studying how the structure carved by planets affects the shape of the disk's spectral energy distribution (SED), and consequently if the SED can be used to infer the presence of planets. We numerically calculate the equilibrium spatial density distributions and SEDs of dust disks originated by a belt of planetesimals in the presence of interior giant planets in different planetary configurations, and for a representative sample of chemical compositions. The dynamical models are necessary to estimate the enhancement of particles near the mean motion resonances with the planets, and to determine how many particles drift inside the planet's orbit. Based on the SEDs and predicted $\it{Spitzer}$ colors we discuss what types of planetary systems can be distinguishable from one another and the main parameter degeneracies in the model SEDs.

Main sequence stars are commonly surrounded by debris disks, formed by cold far-IR-emitting dust that is thought to be continuously replenished by a reservoir of undetected dust-producing planetesimals. We have investigated the orbital evolution of dust particles in debris disks harboring massive planets. Small dust grains are blown out by radiation pressure, as is well known; in addition, gravitational scattering by the giant planets also creates an outflow of large grains. We describe the characteristics of this large-particle outflow in different planetary architectures and for different particle sizes. In addition, the ejection of particles is responsible for the clearing of dust inside the orbit of the planet. We study the efficiency of particle ejection and the resulting dust density contrast inside and outside the orbit of the planet, as a function of the planet's mass and orbital elements and the particle size. We discuss its implications for exo-planetary debris disks and for the interpretation of in-situ dust detection experiments on space probes traveling in the outer solar system.

We present the discovery of debris systems around three solar mass stars based upon observations performed with the Spitzer Space Telescope as part of a Legacy Science Program, ``the Formation and Evolution of Planetary Systems'' (FEPS). We also confirm the presence of debris around two other stars. All the stars exhibit infrared emission in excess of the expected photospheres in the 70 micron band, but are consistent with photospheric emission at <= 33 micron. This restricts the maximum temperature of debris in equilibrium with the stellar radiation to T < 70 K. We find that these sources are relatively old in the FEPS sample, in the age range 0.7 - 3 Gyr. Based on models of the spectral energy distributions, we suggest that these debris systems represent materials generated by collisions of planetesimal belts. We speculate on the nature of these systems through comparisons to our own Kuiper Belt, and on the likely planet(s) responsible for stirring the system and ultimately releasing dust through collisions. We further report observations of a nearby star HD 13974 (d =11 pc) that is indistinguishable from a bare photosphere at both 24 micron and 70 micron. The observations place strong upper limits on the presence of any cold dust in this nearby system (L_IR/L_* < 10^{-5.2}).

[Condensed] We search 0.02 deg^2 for trans-Neptunian objects (TNOs) with m<=29.2 (diameter ~15 km) using the ACS on HST. Three new objects are discovered, roughly 25 times fewer than expected from extrapolation of the differential sky density Sigma(m) of brighter objects. The ACS and other recent TNO surveys show departures from a power law size distribution. Division of the TNO sample into ``classical Kuiper belt'' (CKB) and ``Excited'' samples reveals that Sigma(m) differs for the two populations at 96% confidence. A double power law adequately fits all data. Implications include: The total mass of the CKB is ~0.010 M_Earth, only a few times Pluto's mass, and is predominately in the form of ~100 km bodies. The mass of Excited objects is perhaps a few times larger. The Excited class has a shallower bright-end size distribution; the largest objects, including Pluto, comprise tens of percent of the total mass whereas the largest CKBOs are only ~2% of its mass. The predicted mass of the largest Excited body is close to the Pluto mass; the largest CKBO is ~60 times less massive. The deficit of small TNOs occurs for sizes subject to disruption by present-day collisions, suggesting extensive depletion by collisions. Both accretion and erosion appearing to have proceeded to more advanced stages in the Excited class than the CKB. The absence of distant TNOs implies that any distant (60 AU) population must have less than the CKB mass in the form of objects 40 km or larger. The CKB population is sparser than theoretical estimates of the required precursor population for short period comets, but the Excited population could be a viable precursor population.

We present 3-160 micron photometry obtained with the IRAC and MIPS instruments for the first five targets from the Spitzer Legacy Science Program "Formation and Evolution of Planetary Systems" and 4-35 micron spectro-photometry obtained with the IRS for two sources. We discuss in detail our observations of the debris disks surrounding HD 105 (G0V, 30 +- 10 Myr) and HD 150706 (G3V, ~ 700 +- 300 Myr). For HD 105, possible interpretations include large bodies clearing the dust inside of 45 AU or a reservoir of gas capable of sculpting the dust distribution. The disk surrounding HD 150706 also exhibits evidence of a large inner hole in its dust distribution. Of the four survey targets without previously detected IR excess, spanning ages 30 Myr to 3 Gyr, the new detection of excess in just one system of intermediate age suggests a variety of initial conditions or divergent evolutionary paths for debris disk systems orbiting solar-type stars.

(abbreviated) We use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in eccentric co-planar orbits. We show that the secular dynamics of this system can be described using only two parameters, the ratios of the semimajor axes and the planetary masses. The main dynamical features of the system are presented in geometrical pictures that allows us to investigate a large domain of the secular phase space of this three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. We show that three different regimes of secular motion with respect to the secular angle $\Delta\varpi$ are possible: they are circulation, oscillation (around 0 and $180^\circ$), and high eccentricity libration in a nonlinear secular resonance. The first two regimes are a continuous extension of the classical linear secular perturbation theory; the last is a new feature, hitherto unknown, in the secular dynamics of the three-body problem.

We have numerically investigated the long term dynamical behavior of known Centaurs. This class of objects is thought to constitute the transitional population between the Kuiper Belt and the Jupiter-family comets (JFCs). In our study, we find that over their dynamical lifetimes, these objects diffuse into the JFCs and other sinks, and also make excursions into the Scattered Disk, but (not surprisingly) do not diffuse into the parameter space representing the main Kuiper Belt. These Centaurs spend most of their dynamical lifetimes in orbits of eccentricity 0.2-to-0.6 and perihelion distance 12-to-30 AU. Their orbital evolution is characterized by frequent close encounters with the giant planets. Most of these Centaurs will escape from the solar system (or enter the Oort Cloud), while a fraction will enter the JFC population and a few percent will impact a giant planet. Their median dynamical lifetime is 9 Myr, although there is a wide dispersion in lifetimes, ranging from less than 1 Myr to more than 100 Myr. We find the dynamical evolution of this sample of Centaurs to be less orderly than the planet-to-planet "hand-off" described in previous investigations. We discuss the implications of our study for the spatial distribution of the Centaurs as a whole.

We investigate the survivability of Trojan-type companions of Neptune during primordial radial migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. Loss of Neptune Trojans during planetary migration is not a random diffusion process. Rather, losses occur almost exclusively during discrete episodes when Trojan particles are swept by secondary resonances associated with mean-motion commensurabilities of Uranus with Neptune. The single greatest episode of loss ejects nearly 75% of existing Neptune Trojans and occurs just prior to Neptune reaching its final orbit.

Almost all of the $>600$ known Kuiper belt objects (KBOs) have been discovered within 50 AU of the Sun. One possible explanation for the observed lack of KBOs beyond 50 AU is that the distant Kuiper belt is dynamically very cold, and thus thin enough on the sky to have slipped between previous deep survey fields. We have completed a survey designed to search for a dynamically cold distant Kuiper belt near the invariable plane of the Solar system. In $2.3 {\rm deg}^2$ we have discovered a total of 33 KBOs and 1 Centaur, but no objects in circular orbits beyond 50 AU. We find that we can exclude at 95% CL the existence of a distant disk inclined by $i\le1\arcdeg$ to the invariable plane and containing more than 1.2 times as many $D>185$ km KBOs between 50 and 60 AU as the observed inner Kuiper belt, if the distant disk is thinner than $\sigma=1\fdg75$.

The dust produced in the Kuiper Belt (KB) spreads throughout the Solar System forming a dust disk. We numerically model the orbital evolution of KB dust and estimate its equilibrium spatial distribution and its brightness and spectral energy distributions (SED), assuming greybody absorption and emission by the dust grains. We show that the planets modify the KB disk SED, so potentially we can infer the presence of planets in spatially unresolved debris disks by studying the shape of their SEDs. We point out that there are inherent uncertainties in the prediction of structure in the dust disk, owing to the chaotic dynamics of dust orbital evolution imposed by resonant gravitational perturbations of the planets.

We show that in a system of two planets initially in nearly circular orbits, an impulse perturbation that imparts a finite eccentricity to one planet's orbit causes the other planet's orbit to become eccentric as well, and also naturally results in a libration of their relative apsidal longitudes for a wide range of initial conditions. We suggest that such a mechanism may explain orbital eccentricities and apsidal resonance in some exo-planetary systems. The eccentricity impulse could be caused by the ejection of a planet from these systems, or by torques from a primordial gas disk. The amplitude of secular variations provides an observational constraint on the dynamical history of such systems.

We will utilize the sensitivity of SIRTF through the Legacy Science Program to carry out spectrophotometric observations of solar-type stars aimed at (1) defining the timescales over which terrestrial and gas giant planets are built, from measurements diagnostic of dust/gas masses and radial distributions; and (2) establishing the diversity of planetary architectures and the frequency of planetesimal collisions as a function of time through observations of circumstellar debris disks. Together, these observations will provide an astronomical context for understanding whether our solar system - and its habitable planet - is a common or a rare circumstance. Achieving our science goals requires measuring precise spectral energy distributions for a statistically robust sample capable of revealing evolutionary trends and the diversity of system outcomes. Our targets have been selected from two carefully assembled databases of solar-like stars: (1) a sample located within 50 pc of the Sun spanning an age range from 100-3000 Myr for which a rich set of ancillary measurements (e.g. metallicity, stellar activity, kinematics) are available; and (2) a selection located between 15 and 180 pc and spanning ages from 3 to 100 Myr. For stars at these distances SIRTF is capable of detecting stellar photospheres with SNR >30 at lambda < 24 microns for our entire sample, as well as achieving SNR >5 at the photospheric limit for over 50% of our sample at lambda=70 microns. Thus we will provide a complete census of stars with excess emission down to the level produced by the dust in our present-day solar system. More information concerning our program can be found at: http://gould.as.arizona.edu/feps

The population of the Kuiper Belt within 50 AU of the Sun has likely been severely depleted by gravitational perturbations from the giant planets, particularly Neptune. The density of Kuiper Belt objects is expected to be two orders of magnitude higher just beyond 50 AU, where planetary perturbations are insignificant. In 1998 and 1999, we surveyed for Kuiper Belt Objects (KBOs) in 6 fields of the ecliptic (total sky area 1.5 deg^2) to limiting magnitudes between R=24.9 and R=25.9. This is deep enough to detect KBOs of diameter >~ 160 km at a distance of 65 AU. We detected 24 objects. None of these objects, however, is beyond 53 AU. Our survey places a 95% CL upper limit of Sigma<5 deg^-2 on the surface density of KBOs larger than ~160 km beyond 55 AU. This can be compared to the surface density of ~6 deg^-2 of 160 km KBOs at distances 30--50 AU determined from this survey and previous shallower surveys. The mean volume density of D>160 km KBOs in the 55--65 AU region is, at >95% confidence, less than the mean density in the 30--50 AU region, and at most twothirds of the mean density from 40--50 AU. Thus, a substantial density increase beyond 50 AU is excluded in this model-independent estimate, implying that some process or event in the history of the Solar System has truncated the distribution of 160-km planetesimals at ~50 AU. A dense primordial disk could be present beyond 50 AU if it contains only smaller objects, or is sufficiently thin and inclined to have escaped detection in our 6 survey fields.

Radial velocity observations of the F8 V star Upsilon Andromedae taken at Lick and at Whipple Observatories have revealed evidence of three periodicities in the line-of-sight velocity of the star. These periodicities have been interpreted as evidence for at least three low mass companions (LMCs) revolving around Upsilon Andromedae. The mass and orbital parameters inferred for these companions raise questions about the dynamical stability of the system. We report here results from our independent analysis of the published radial velocity data as well as new unpublished data taken at Lick Observatory. Our results confirm the finding of three periods in the data. Our best fits to the data, on the assumption that these periods arise from the gravitational perturbations of companions in keplerian orbits, is also generally in agreement, but with some differences, from the earlier findings. We find that the available data do not constrain well the orbital eccentricity of the middle companion in a three-companion model of the data. We also find that in order for our best-fit model to the Lick data to be dynamically stable over the lifetime of the star (~2 billion years), the system must have a mean inclination to the plane of the sky greater than 13 degrees. The corresponding minimum inclination for the best fit to the Whipple data set is 19 degrees. These values imply that the maximum mass for the outer companion can be no greater than about 20 Jupiter masses. Our analysis of the stability of the putative systems also places constraints on the relative inclinations of the orbital planes of the companions. We comment on global versus local (i.e., method of steepest descent) means of finding best-fit orbits from radial velocity data sets.

The existence of the Oort Comet Cloud, the Kuiper Belt, and plausible inefficiencies in planetary core formation, all suggest that there was once a residual planetesimal disk of mass 10-100 Earth-masses in the vicinity of the giant planets following their formation. Since removal of this disk requires an exchange of orbital energy and angular momentum with the planets, significant planetary migration can ensue. The planet migration phenomenon is examined numerically by evolving the orbits of the giant planets while they are embedded in a planetesimal disk having a mass of M_d=10 to 200 Earth-masses. We find that Saturn, Uranus, and Neptune evolve radially outwards as they scatter the planetesimals, while Jupiter's orbit shrinks as it ejects mass. Higher-mass disks result in more rapid and extensive planet migration. If orbit expansion and resonance trapping by Neptune is invoked to explain the eccentricities of Pluto and its cohort of Kuiper Belt Objects at Neptune's 3:2 mean-motion resonance, then our simulations suggest that a disk mass of order M_d~50 Earth-masses is required to expand Neptune's orbit by ~7 AU in order to pump up Plutino eccentricities to e~0.3. Such planet migration implies that the initial Solar System was more compact in the past, with the Jupiter-Neptune separation having been smaller by about 30%. The planetesimal disk is also the source of the Oort Cloud of comets. Using the results of our simulations together with a simple treatment of Oort Cloud dynamics, we estimate that ~12 Earth-masses of disk material was initially deposited in the Oort Cloud, of which ~4 Earth-masses will persist over the age of the Solar System. The majority of these comets originated from the Saturn-Neptune region of the solar nebula.

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 Solar system beyond Neptune is believed to house a population of small primordial bodies left over from the planet formation process. The region up to heliocentric distance $\sim 50 AU$ (a.k.a. the Kuiper Belt) may be the source of the observed short period comets. In this region, the phase space structure near orbital resonances with Neptune is of special interest for the long term stability of orbits. There is reason to believe that a significant fraction (perhaps most) of the Kuiper Belt objects reside preferentially in these resonance locations. This paper describes the dynamics of small objects near the major orbital resonances with Neptune. Estimates of the widths of stable resonance zones as well as the properties of resonant orbits are obtained from the circular, planar restricted three-body model. Although this model does not contain the full complexity of the long term orbital dynamics of Kuiper Belt objects subject to the full N-body perturbations of all the planets, it does provide a baseline for the phase space structure and properties of resonant orbits in the trans-Neptunian Solar system.

The origin of the highly eccentric, inclined, and resonance-locked orbit of Pluto has long been a puzzle. A possible explanation has been proposed recently [Malhotra, R., {\it Nature} 365:819-21 (1993)] which suggests that these extraordinary orbital properties may be a natural consequence of the formation and early dynamical evolution of the outer Solar system. A resonance capture mechanism is possible during the clearing of the residual planetesimal debris and the formation of the Oort Cloud of comets by planetesimal mass loss from the vicinity of the giant planets. If this mechanism were in operation during the early history of the planetary system, the entire region between the orbit of Neptune and approximately 50 AU would have been swept by first order mean motion resonances. Thus, resonance capture could occur not only for Pluto, but quite generally for other trans-Neptunian small bodies. Some consequences of this evolution for the present-day dynamical structure of the trans-Neptunian region are: (i) most of the objects in the region beyond Neptune and up to $\sim\!50$ AU exist in very narrow zones located at orbital resonances with Neptune (particularly the 3:2 and the 2:1 resonances), and (ii) these resonant objects would have significantly large eccentricities. The distribution of objects in the Kuiper Belt as predicted by this theory is presented here. Keywords{solar system: origins, dynamics --- Pluto --- comets --- Kuiper Belt}

Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ODE integration methods. Here I present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.

Orbital resonances are ubiquitous in the Solar system. They play a decisive role in the long term dynamics, and in some cases the physical evolution, of the planets and of their natural satellites, as well as the evolution of small bodies (including dust) in the planetary system. The few-body gravitational problem of hierarchical planetary-type systems allows for a complex range of dynamical timescales, from the fast orbital periods to the very slow orbit precession rates. The interaction of fast and slow degrees of freedom produces a rich diversity of resonance phenomena. Weak dissipative effects --- such as tides or radiation drag forces --- also produce unexpectedly rich dynamical behaviors. This paper provides a mostly qualitative discussion of simple dynamical models for the commonly encountered orbital resonance phenomena in the Solar system.

About two years ago, Wolszczan and Frail announced the detection of a possible planetary system consisting of two Earth-mass planets around a millisecond pulsar. It was pointed out shortly thereafter that the mutual gravitational interaction of the planets leads to predictable evolution of their orbits. Because millisecond pulsars are very stable ``clocks'' and allow pulse timing measurements with microsecond precision, the orbit evolution would lead to a measurable modulation of the arrival times of the pulses. Detection of this effect would serve to corroborate the planetary interpretation of the data and provide additional constraints on the system. The theory of planetary perturbations in the context of pulsar timing observations is summarized here.