We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg$_2$SiO$_4$), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

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

A group of 72 members of the astrochemistry and spectroscopy communities assembled at Emory University in October, 2012, to discuss the state-of-the art in molecular laboratory spectroscopic techniques used to study molecules of astrophysical interest, and highlight recent successes in laboratory efforts that complement new astronomical observations. The format of the meeting facilitated detailed discussion of the needs for collecting and interpreting molecular spectroscopic information to support observational astrophysics in the era of Far-IR astronomy. We present the findings and recommendations from this workshop in this report. The participants in this workshop identified the development of new laboratory capabilities that offer rapid, broadband, high-resolution, high-sensitivity measurements of molecular spectra in the frequency ranges that overlap with new observational facilities as the top priority for far-infrared laboratory astrophysics initiatives over the next 5-10 years. An additional major need is the development of suitable analysis tools for broadband molecular spectra from both the laboratory and space; though such efforts should be conducted within the framework of larger collaborative funding structures, rather than through single-PI efforts, so as to provide tools that are generally applicable rather than highly specialized. Additional efforts in theory, computation, and modeling that support these goals are also important for advancing this field.

Nearly half of the exoplanets found within binary star systems reside in very wide binaries with average stellar separations beyond 1,000 AU (1 AU being the Earth-Sun distance), yet the influence of such distant binary companions on planetary evolution remains largely unstudied. Unlike their tighter counterparts, the stellar orbits of wide binaries continually change under the influence of the Galactic tide and impulses from other passing stars. Here we report numerical simulations demonstrating that the variable nature of wide binary star orbits dramatically reshapes the planetary systems they host, typically Gyrs after formation. Contrary to previous understanding, wide binary companions may often strongly perturb planetary systems, triggering planetary ejections and exciting orbital eccentricities of surviving planets. Indeed, observed exoplanet eccentricities offer evidence of this; giant exoplanet orbits within wide binaries are statistically more eccentric than those around isolated stars. Both eccentricity distributions are well-reproduced when we assume isolated stars and wide binaries host similar planetary systems whose outermost giant planets are scattered beyond ~10 AU from their parent stars via early internal instabilities. Consequently, our results suggest that although wide binaries eventually truncate their planetary systems, most isolated giant exoplanet systems harbor additional distant, still undetected planets.

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

Although the 55 Cnc system contains multiple, closely packed planets that are presumably in a coplanar configuration, we use numerical simulations to demonstrate that they are likely to be highly inclined to their parent star's spin axis. Due to perturbations from its distant binary companion, this planetary system precesses like a rigid body about its parent star. Consequently, the parent star's spin axis and the planetary orbit normal likely diverged long ago. Because only the projected separation of the binary is known, we study this effect statistically, assuming an isotropic distribution for wide binary orbits. We find that the most likely projected spin-orbit angle is ~50 degrees, with a ~30% chance of a retrograde configuration. Transit observations of the innermost planet - 55 Cnc e - may be used to verify these findings via the Rossiter-McLaughlin effect. 55 Cancri may thus represent a new class of planetary systems with well-ordered, coplanar orbits that are inclined with respect to the stellar equator.

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

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

Using orbital integration and analytical arguments, we have found a new mechanism (an "eccentricity trap") to halt type I migration of planets near the inner edge of a protoplanetary disk. Because asymmetric eccentricity damping due to disk-planet interaction on the innermost planet at the disk edge plays a crucial role in the trap, this mechanism requires continuous eccentricity excitation and hence works for a resonantly interacting convoy of planets. This trap is so strong that the edge torque exerted on the innermost planet can completely halt type I migrations of many outer planets through mutual resonant perturbations. Consequently, the convoy stays outside the disk edge, as a whole. We have derived semi-analytical formula for the condition for the eccentricity trap and predict how many planets are likely to be trapped. We found that several planets or more should be trapped by this mechanism in protoplanetary disks that have cavities. It can be responsible for the formation of non-resonant, multiple, close-in super-Earth systems extending beyond 0.1AU. Such systems are being revealed by radial velocity observations to be quite common around solar-type stars.

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

More than 100 extrasolar planets have been discovered since 1990s. Different from the solar system, these planets' orbital eccentricities cover a huge range from 0 to 0.7. Incidently, the first Kuiper Belt Object was discovered in 1992. Thus, an interesting and important question will be whether extrasolar planetary systems could have structures like Kuiper Belt or asteroid belt. We investigate the stability of these planetary systems with different orbital eccentricities by the similar procedures in Rabl & Dvorak (1988) and Holman & Wiegert (1999). We claim that most extrasolar planetary systems can have their own belts at the outer regions. However, we find that the orbits with high--eccentricity is very powerful in depletion of these populations.

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

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

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

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

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

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

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