We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the project. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 $\mu$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $\mu$m to constrain the brightness temperatures of their daysides. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We estimate for TRAPPIST-1 transit depth measurements mean noise floors of $\sim$35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 $\mu$m, and can only set 3$\sigma$ upper limits on their dayside brightness temperatures (611K for b 586K for c).

In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not representative of the past evolution of real bodies. Instead, these simulations are only useful as a means to quantify the short dynamical lifetime of the considered bodies and the fast decay of their population. In light of this fast decay, if the observed bodies were the survivors of populations of objects captured from interstellar space in the early Solar System, these populations should have been implausibly large (e.g. about 10 times the current main asteroid belt population for the retrograde coorbital of Jupiter). More likely, the observed objects are just transient members of a population that is maintained in quasi-steady state by a continuous flux of objects from some parent reservoir in the distant Solar System. We identify in the Halley type comets and the Oort cloud the most likely sources of retrograde coorbitals and highly-inclined Centaurs.

Models of planet formation are built on underlying physical processes. In order to make sense of the origin of the planets we must first understand the origin of their building blocks. This review comes in two parts. The first part presents a detailed description of six key mechanisms of planet formation: 1) The structure and evolution of protoplanetary disks 2) The formation of planetesimals 3) Accretion of protoplanets 4) Orbital migration of growing planets 5) Gas accretion and giant planet migration 6) Resonance trapping during planet migration. While this is not a comprehensive list, it includes processes for which our understanding has changed in recent years or for which key uncertainties remain. The second part of this review shows how global models are built out of planet formation processes. We present global models to explain different populations of known planetary systems, including close-in small/low-mass planets (i.e., super-Earths), giant exoplanets, and the Solar System's planets. We discuss the different sources of water on rocky exoplanets, and use cosmochemical measurements to constrain the origin of Earth's water. We point out the successes and failings of different models and how they may be falsified. Finally, we lay out a path for the future trajectory of planet formation studies.

The dynamical structure of the Solar System can be explained by a period of orbital instability experienced by the giant planets. While a late instability was originally proposed to explain the Late Heavy Bombardment, recent work favors an early instability. We model the early dynamical evolution of the outer Solar System to self-consistently constrain the most likely timing of the instability. We first simulate the dynamical sculpting of the primordial outer planetesimal disk during the accretion of Uranus and Neptune from migrating planetary embryos during the gas disk phase, and determine the separation between Neptune and the inner edge of the planetesimal disk. We performed simulations with a range of migration histories for Jupiter. We find that, unless Jupiter migrated inwards by 10 AU or more, the instability almost certainly happened within 100 Myr of the start of Solar System formation. There are two distinct possible instability triggers. The first is an instability that is triggered by the planets themselves, with no appreciable influence from the planetesimal disk. Of those, the median instability time is $\sim4$Myr. Among self-stable systems -- where the planets are locked in a resonant chain that remains stable in the absence of a planetesimal's disk-- our self-consistently sculpted planetesimal disks nonetheless trigger a giant planet instability with a median instability time of 37-62 Myr for a reasonable range of migration histories of Jupiter. The simulations that give the latest instability times are those that invoked long-range inward migration of Jupiter from 15 AU or beyond; however these simulations over-excited the inclinations of Kuiper belt objects and are inconsistent with the present-day Solar System. We conclude on dynamical grounds that the giant planet instability is likely to have occurred early in Solar System history.

We review the state of knowledge on the origin of Earth's water. Empirical constraints come from chemical and isotopic measurements of solar system bodies and of Earth itself. Dynamical models have revealed pathways for water delivery to Earth during its formation; most are anchored to specific models for terrestrial planet formation. Meanwhile, disk chemical models focus on determining how the isotopic ratios of the building blocks of planets varied as a function of radial distance and time, defining markers of material transported along those pathways. Carbonaceous chondrite meteorites -- representative of the outer asteroid belt -- match Earth's water isotopes (although mantle plumes have been measured at lower D/H). But how was this connection established -- did Earth's water originate among the asteroids (as in the classical model of terrestrial planet formation)? Or, more likely, was Earth's water delivered from the same parent population as the hydrated asteroids (e.g., external pollution, as in the Grand Tack model)? We argue that the outer asteroid belt -- the boundary between the inner and outer solar system -- is the next frontier for new discoveries. The outer asteroid belt contains a population of icy bodies with volatile-driven activity seen on twelve main belt comets (MBCs); seven of which exhibit sublimation-driven activity on repeated perihelion passages. Measurements of the isotopic characteristics of MBCs would provide essential missing links in the chain between disk models and dynamical models. Finally, we extrapolate to rocky exoplanets. Migration is the only mechanism likely to produce very water-rich planets with more than a few percent water by mass (and even with migration, some planets are pure rock). While water loss mechanisms remain to be studied in more detail, we expect that water should be delivered to the vast majority of rocky exoplanets.

The asteroid belt is characterized by an extreme low total mass of material on dynamically excited orbits. The Nice Model explains many peculiar qualities of the solar system, including the belt's excited state, by invoking an orbital instability between the outer planets. However, previous studies of the Nice Model's effect on the belt's structure struggle to reproduce the innermost asteroids' orbital inclination distribution. Here, we show how the final phase of giant planet migration sculpts the asteroid belt, in particular its inclination distribution. As interactions with leftover planetesimals cause Saturn to move away from Jupiter, its rate of orbital precession slows as the two planets' mutual interactions weaken. When the planets approach their modern separation, where Jupiter completes just short of five orbits for every two of Saturn's, Jupiter's eccentric forcing on Saturn strengthens. We use numerical simulations to show that the absence of asteroids with orbits that precess between 24-28 arcsec/yr is related to the inclination problem. As Saturn's precession speeds back up, high inclination asteroids are excited on to planet crossing orbits and removed from the inner main belt. Through this process, the asteroid belt's orbital structure is reshaped, leading to markedly improved simulation outcomes.

Most giant exoplanets discovered by radial velocity surveys have much higher eccentricities than those in the solar system. The planet--planet scattering mechanism has been shown to match the broad eccentricity distribution, but the highest-eccentricity planets are often attributed to Kozai-Lidov oscillations induced by a stellar companion. Here we investigate whether the highly eccentric exoplanet population can be produced entirely by scattering. We ran 500 N-body simulations of closely packed giant-planet systems that became unstable under their own mutual perturbations. We find that the surviving bound planets can have eccentricities up to $e > 0.99$, with a maximum of 0.999017 in our simulations. This suggests that there is no maximum eccentricity that can be produced by planet--planet scattering. Importantly, we find that extreme eccentricities are not extremely rare; the eccentricity distribution for all giant exoplanets with $e > 0.3$ is consistent with all planets concerned being generated by scattering. Our results show that the discovery of planets with extremely high eccentricities does not necessarily signal the action of the Kozai-Lidov mechanism.

In the innermost regions of protoplanerary discs, the solid-to-gas ratio can be increased considerably by a number of processes, including photoevaporative and particle drift. MHD disc models also suggest the existence of a dead-zone at $R\lesssim 10$ AU, where the regions close to the midplane remain laminar. In this context, we use two-fluid hydrodynamical simulations to study the interaction between a low-mass planet ($\sim 1.7 \;{\rm M_\oplus}$) on a fixed orbit and an inviscid pebble-rich disc with solid-to-gas ratio $\epsilon\ge 0.5$. For pebbles with Stokes numbers St=0.1, 0.5, multiple dusty vortices are formed through the Rossby Wave Instability at the planet separatrix. Effects due to gas drag then lead to a strong enhancement in the solid-to-gas ratio, which can increase by a factor of $\sim 10^3$ for marginally coupled particles with St=0.5. As in streaming instabilities, pebble clumps reorganize into filaments that may plausibly collapse to form planetesimals. When the planet is allowed to migrate in a MMSN disc, the vortex instability is delayed due to migration but sets in once inward migration stops due a strong positive pebble torque. Again, particle filaments evolving in a gap are formed in the disc while the planet undergoes an episode of outward migration. Our results suggest that vortex instabilities triggered by low-mass planets could play an important role in forming planetesimals in pebble-rich, inviscid discs, and may significantly modify the migration of low-mass planets. They also imply that planetary dust gaps may not necessarily contain planets if these migrated away.

The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

Recent observations have found a valley in the size distribution of close-in super-Earths that is interpreted as a signpost that close-in super-Earths are mostly rocky in composition. However, new models predict that planetesimals should first form at the water ice line such that close-in planets are expected to have a significant water ice component. Here we investigate the water contents of super-Earths by studying the interplay between pebble accretion, planet migration and disc evolution. Planets' compositions are determined by their position relative to different condensation fronts (ice lines) throughout their growth. Migration plays a key role. Assuming that planetesimals start at or exterior to the water ice line ($r>r_{\rm H_2O}$), inward migration causes planets to leave the source region of icy pebbles and therefore to have lower final water contents than in discs with either outward migration or no migration. The water ice line itself moves inward as the disc evolves, and delivers water as it sweeps across planets that formed dry. The relative speed and direction of planet migration and inward drift of the water ice line is thus central in determining planets' water contents. If planet formation starts at the water ice line, this implies that hot close-in super-Earths (r<0.3 AU) with water contents of a few percent are a signpost of inward planet migration during the early gas phase. Hot super-Earths with larger water ice contents on the other hand, experienced outward migration at the water ice line and only migrated inwards after their formation was complete either because they become too massive to be contained in the region of outward migration or in chains of resonant planets. Measuring the water ice content of hot super-Earths may thus constrain their migration history.

Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nominal type-I and type-II migration rates and a pebble flux of 100-200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30AU migrate into the inner system (r<1AU) and form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explore slower migration rates which allow the formation of gas giants from embryos originating from the 5-10AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. Giant planet growth requires a sufficiently-high pebble flux to enable growth to out-compete migration. Even higher pebble fluxes produce systems with multiple gas giants. We show that planetary embryos starting interior to 5AU do not grow into gas giants, even if migration is slow and the pebble flux is large. Instead they grow to the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is independent of the pebble flux. Additionally we show that the long term evolution of our formed planetary systems can produce systems with hot super-Earths and outer gas giants as well as systems of giants on eccentric orbits (abridged).

At least 30% of main sequence stars host planets with sizes between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk-planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes ($\lesssim15M_{\oplus}$) and super-massive planetary cores potentially able to become gas giant planets ($\gtrsim15M_{\oplus}$). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the disk's inner edge forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to $\sim$95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the resonant configuration. Our simulations match observations if we combine a dominant fraction ($\gtrsim95\%$) of unstable systems with a sprinkling ($\lesssim5\%$) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (${\rm N\geq2}$) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet-planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice-rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths (abridged)

Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation pathways that are regulated by the disc's pebble reservoir. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass-flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass-flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of Mars-mass embryos when the gas disc dissipates. Without gas being present, the embryos become unstable and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most 5 Earth masses. Instead, if the pebble mass-flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the 5 to 20 Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. The pebble flux - which controls the transition between the two growth modes - may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.

Each of the giant planets within the Solar System has large moons but none of these moons have their own moons (which we call ${\it submoons}$). By analogy with studies of moons around short-period exoplanets, we investigate the tidal-dynamical stability of submoons. We find that 10 km-scale submoons can only survive around large (1000 km-scale) moons on wide-separation orbits. Tidal dissipation destabilizes the orbits of submoons around moons that are small or too close to their host planet; this is the case for most of the Solar System's moons. A handful of known moons are, however, capable of hosting long-lived submoons: Saturn's moons Titan and Iapetus, Jupiter's moon Callisto, and Earth's Moon. Based on its inferred mass and orbital separation, the newly-discovered exomoon candidate Kepler-1625b-I can in principle host a large submoon, although its stability depends on a number of unknown parameters. We discuss the possible habitability of submoons and the potential for subsubmoons. The existence, or lack thereof, of submoons, may yield important constraints on satellite formation and evolution in planetary systems.

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

Exoplanet surveys have confirmed one of humanity's (and all teenagers') worst fears: we are weird. If our Solar System were observed with present-day Earth technology -- to put our system and exoplanets on the same footing -- Jupiter is the only planet that would be detectable. The statistics of exo-Jupiters indicate that the Solar System is unusual at the ~1% level among Sun-like stars (or ~0.1% among all stars). But why are we different? Successful formation models for both the Solar System and exoplanet systems rely on two key processes: orbital migration and dynamical instability. Systems of close-in super-Earths or sub-Neptunes require substantial radial inward motion of solids either as drifting mm- to cm-sized pebbles or migrating Earth-mass or larger planetary embryos. We argue that, regardless of their formation mode, the late evolution of super-Earth systems involves migration into chains of mean motion resonances, generally followed by instability when the disk dissipates. This pattern is likely also ubiquitous in giant planet systems. We present three models for inner Solar System formation -- the low-mass asteroid belt, Grand Tack, and Early Instability models -- each invoking a combination of migration and instability. We identify bifurcation points in planetary system formation. We present a series of events to explain why our Solar System is so weird. Jupiter's core must have formed fast enough to quench the growth of Earth's building blocks by blocking the flux of inward-drifting pebbles. The large Jupiter/Saturn mass ratio is rare among giant exoplanets but may be required to maintain Jupiter's wide orbit. The giant planets' instability must have been gentle, with no close encounters between Jupiter and Saturn, also unusual in the larger (exoplanet) context. Our Solar System system is thus the outcome of multiple unusual, but not unheard of, events.

Containing only a few percent the mass of the moon, the current asteroid belt is around three to four orders of magnitude smaller that its primordial mass inferred from disk models. Yet dynamical studies have shown that the asteroid belt could not have been depleted by more than about an order of magnitude over the past ~4 Gyr. The remainder of the mass loss must have taken place during an earlier phase of the solar system's evolution. An orbital instability in the outer solar system occurring during the process of terrestrial planet formation can reproduce the broad characteristics of the inner solar system. Here, we test the viability of this model within the constraints of the main belt's low present-day mass and orbital structure. While previous studies modeled asteroids as massless test particles because of limited computing power, our work uses GPU (Graphics Processing Unit) acceleration to model a fully self-gravitating asteroid belt. We find that depletion in the main belt is related to the giant planets' exact evolution within the orbital instability. Simulations that produce the closest matches to the giant planets' current orbits deplete the main belt by two to three orders of magnitude. These simulated asteroid belts are also good matches to the actual asteroid belt in terms of their radial mixing and broad orbital structure.

We present empirical evidence, supported by a planet formation model, to show that the curve $R/R_\oplus = 1.05\,(F/F_\oplus)^{0.11}$ approximates the location of the so-called photo-evaporation valley. Planets below that curve are likely to have experienced complete photo-evaporation, and planets just above it appear to have inflated radii; thus we identify a new population of inflated super-Earths and mini-Neptunes. Our N-body simulations are set within an evolving protoplanetary disk and include prescriptions for orbital migration, gas accretion, and atmospheric loss due to giant impacts. Our simulated systems broadly match the sizes and periods of super-Earths in the Kepler catalog. They also reproduce the relative sizes of adjacent planets in the same system, with the exception of planet pairs that straddle the photo-evaporation valley. This latter group is populated by planet pairs with either very large or very small size ratios ($R_{\rm out} / R_{\rm in} \gg 1$ or $R_{\rm out} / R_{\rm in} \ll 1$) and a dearth of size ratios near unity. It appears that this feature could be reproduced if the planet outside the photo-evaporation valley (typically the outer planet, but some times not) has its atmosphere significantly expanded by stellar irradiation. This new population of planets may be ideal targets for future transit spectroscopy observations with the upcoming James Webb Space Telescope.

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

A leading model for the origin of super-Earths proposes that planetary embryos migrate inward and pile up on close-in orbits. As large embryos are thought to preferentially form beyond the snow line, this naively predicts that most super-Earths should be very water-rich. Here we show that the shortest-period planets formed in the migration model are often purely rocky. The inward migration of icy embryos through the terrestrial zone accelerates the growth of rocky planets via resonant shepherding. We illustrate this process with a simulation that provided a match to the Kepler-36 system of two planets on close orbits with very different densities. In the simulation, two super-Earths formed in a Kepler-36-like configuration; the inner planet was pure rock while the outer one was ice-rich. We conclude from a suite of simulations that the feeding zones of close-in super-Earths are likely to be broad and disconnected from their final orbital radii.

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

The past decade has seen major progress in our understanding of terrestrial planet formation. Yet key questions remain. In this review we first address the growth of 100 km-scale planetesimals as a consequence of dust coagulation and concentration, with current models favoring the streaming instability. Planetesimals grow into Mars-sized (or larger) planetary embryos by a combination of pebble- and planetesimal accretion. Models for the final assembly of the inner Solar System must match constraints related to the terrestrial planets and asteroids including their orbital and compositional distributions and inferred growth timescales. Two current models -- the Grand-Tack and low-mass (or empty) primordial asteroid belt scenarios -- can each match the empirical constraints but both have key uncertainties that require further study. We present formation models for close-in super-Earths -- the closest current analogs to our own terrestrial planets despite their very different formation histories -- and for terrestrial exoplanets in gas giant systems. We explain why super-Earth systems cannot form in-situ but rather may be the result of inward gas-driven migration followed by the disruption of compact resonant chains. The Solar System is unlikely to have harbored an early system of super-Earths; rather, Jupiter's early formation may have blocked the ice giants' inward migration. Finally, we present a chain of events that may explain why our Solar System looks different than more than 99\% of exoplanet systems.

'Oumuamua was discovered passing through our Solar System on a hyperbolic orbit. It presents an apparent contradiction, with colors similar to those of volatile-rich Solar System bodies but with no visible outgassing or activity during its close approach to the Sun. Here we show that this contradiction can be explained by the dynamics of planetesimal ejection by giant planets. We propose that 'Oumuamua is an extinct fragment of a comet-like planetesimal born in a planet-forming disk that also formed Neptune- to Jupiter-mass giant planets. On its pathway to ejection 'Oumuamua's parent body underwent a close encounter with a giant planet and was tidally disrupted into small pieces, similar to comet Shoemaker-Levy 9's disruption after passing close to Jupiter. We use dynamical simulations to show that 0.1-1% of cometary planetesimals undergo disruptive encounters prior to ejection. Rocky asteroidal planetesimals are unlikely to disrupt due to their higher densities. After disruption, the bulk of fragments undergo enough close passages to their host stars to lose their surface volatiles and become extinct. Planetesimal fragments such as 'Oumuamua contain little of the mass in the population of interstellar objects but dominate by number. Our model makes predictions that will be tested in the coming decade by LSST.

'Oumuamua, the first bona-fide interstellar planetesimal, was discovered passing through our Solar System on a hyperbolic orbit. This object was likely dynamically ejected from an extrasolar planetary system after a series of close encounters with gas giant planets. To account for 'Oumuamua's detection, simple arguments suggest that ~1 Earth mass of planetesimals are ejected per Solar mass of Galactic stars. However, that value assumes mono-sized planetesimals. If the planetesimal mass distribution is instead top-heavy the inferred mass in interstellar planetesimals increases to an implausibly high value. The tension between theoretical expectations for the planetesimal mass function and the observation of 'Oumuamua can be relieved if a small fraction (~0.1-1%) of planetesimals are tidally disrupted on the pathway to ejection into 'Oumuamua-sized fragments. Using a large suite of simulations of giant planet dynamics including planetesimals, we confirm that 0.1-1% of planetesimals pass within the tidal disruption radius of a gas giant on their pathway to ejection. 'Oumuamua may thus represent a surviving fragment of a disrupted planetesimal. Finally, we argue that an asteroidal composition is dynamically disfavoured for 'Oumuamua, as asteroidal planetesimals are both less abundant and ejected at a lower efficiency than cometary planetesimals.

The planetary building blocks that formed in the terrestrial planet region were likely very dry, yet water is comparatively abundant on Earth. We review the various mechanisms proposed for the origin of water on the terrestrial planets. Various in-situ mechanisms have been suggested, which allow for the incorporation of water into the local planetesimals in the terrestrial planet region or into the planets themselves from local sources, although all of those mechanisms have difficulties. Comets have also been proposed as a source, although there may be problems fitting isotopic constraints, and the delivery efficiency is very low, such that it may be difficult to deliver even a single Earth ocean of water this way. The most promising route for water delivery is the accretion of material from beyond the snow line, similar to carbonaceous chondrites, that is scattered into the terrestrial planet region as the planets are growing. Two main scenarios are discussed in detail. First is the classical scenario in which the giant planets begin roughly in their final locations and the disk of planetesimals and embryos in the terrestrial planet region extends all the way into the outer asteroid belt region. Second is the Grand Tack scenario, where early inward and outward migration of the giant planets implants material from beyond the snow line into the asteroid belt and terrestrial planet region, where it can be accreted by the growing planets. Sufficient water is delivered to the terrestrial planets in both scenarios. While the Grand Tack scenario provides a better fit to most constraints, namely the small mass of Mars, planets may form too fast in the nominal case discussed here. This discrepancy may be reduced as a wider range of initial conditions is explored. Finally, we discuss several more recent models that may have important implications for water delivery to the terrestrial planets.

Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting to study the formation and evolution of terrestrial planets that formed in the same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all known to better than 5% precision, their densities have significant uncertainties (between 28% and 95%) because of poor constraints on the planet's masses. Aims.The goal of this paper is to improve our knowledge of the TRAPPIST-1 planetary masses and densities using transit-timing variations (TTV). The complexity of the TTV inversion problem is known to be particularly acute in multi-planetary systems (convergence issues, degeneracies and size of the parameter space), especially for resonant chain systems such as TRAPPIST-1. Methods. To overcome these challenges, we have used a novel method that employs a genetic algorithm coupled to a full N-body integrator that we applied to a set of 284 individual transit timings. This approach enables us to efficiently explore the parameter space and to derive reliable masses and densities from TTVs for all seven planets. Results. Our new masses result in a five- to eight-fold improvement on the planetary density uncertainties, with precisions ranging from 5% to 12%. These updated values provide new insights into the bulk structure of the TRAPPIST-1 planets. We find that TRAPPIST-1\,c and e likely have largely rocky interiors, while planets b, d, f, g, and h require envelopes of volatiles in the form of thick atmospheres, oceans, or ice, in most cases with water mass fractions less than 5%.

Evidence of mutually inclined planetary orbits has been reported for giant planets these last years. Here we aim to study the impact of eccentric and inclined massive giant planets on the terrestrial planet formation process, and investigate whether it can possibly lead to the existence of inclined terrestrial planets. We have performed 126 simulations of the late-stage planetary accretion in eccentric and inclined giant planet systems. The physical and orbital parameters of the giant planet systems result from n-body simulations of three giant planets in the late stage of the gas disc, under the combined action of Type II migration and planet-planet scattering. Fourteen two- and three-planet configurations have been selected, with diversified masses, semi-major axes (resonant configurations or not), eccentricities and inclinations (including coplanar systems) at the dispersal of the gas disc. We have then followed the gravitational interactions of these systems with an inner disc of planetesimals and embryos (9 runs per system), studying in detail the final configurations of the formed terrestrial planets. While coplanar giant systems harbour several terrestrial planets, generally as massive as the Earth and mainly on low eccentric and low inclined orbits, terrestrial planets formed in systems with mutually inclined giant planets are usually fewer, less massive (<0.5 M_{\Earth}) and with larger eccentricities and inclinations. This work shows that terrestrial planets can form on stable inclined orbits through the classical accretion theory, even in coplanar giant planet systems emerging from the disc phase.

Planet-planet scattering is the leading mechanism to explain the broad eccentricity distribution of observed giant exoplanets. Here we study the orbital stability of primordial giant planet moons in this scenario. We use N-body simulations including realistic oblateness and evolving spin evolution for the giant planets. We find that the vast majority (80-90% across all our simulations) of orbital parameter space for moons is destabilized. There is a strong radial dependence, as moons past 0.1 Hill radii are systematically removed. Closer-in moons on Galilean-moon-like orbits (< 0.04 Hill radii) have a good ( 20-40%) chance of survival. Destabilized moons may undergo a collision with the star or a planet, be ejected from the system, be captured by another planet, be ejected but still orbiting its free-floating host planet, or survive on heliocentric orbits as "planets." The survival rate of moons increases with the host planet mass but is independent of the planet's final (post-scattering) orbits. Based on our simulations we predict the existence of an abundant galactic population of free- floating (former) moons.

The asteroid belt contains less than a thousandth of Earth's mass and is radially segregated, with S-types dominating the inner belt and C-types the outer belt. It is generally assumed that the belt formed with far more mass and was later strongly depleted. Here we show that the present-day asteroid belt is consistent with having formed empty, without any planetesimals between Mars and Jupiter's present-day orbits. This is consistent with models in which drifting dust is concentrated into an isolated annulus of terrestrial planetesimals. Gravitational scattering during terrestrial planet formation causes radial spreading, transporting planetesimals from inside 1-1.5 AU out to the belt. Several times the total current mass in S-types is implanted, with a preference for the inner main belt. C-types are implanted from the outside, as the giant planets' gas accretion destabilizes nearby planetesimals and injects a fraction into the asteroid belt, preferentially in the outer main belt. These implantation mechanisms are simple byproducts of terrestrial- and giant planet formation. The asteroid belt may thus represent a repository for planetary leftovers that accreted across the Solar System but not in the belt itself.

There is a long-standing debate regarding the origin of the terrestrial planets' water as well as the hydrated C-type asteroids. Here we show that the inner Solar System's water is a simple byproduct of the giant planets' formation. Giant planet cores accrete gas slowly until the conditions are met for a rapid phase of runaway growth. As a gas giant's mass rapidly increases, the orbits of nearby planetesimals are destabilized and gravitationally scattered in all directions. Under the action of aerodynamic gas drag, a fraction of scattered planetesimals are deposited onto stable orbits interior to Jupiter's. This process is effective in populating the outer main belt with C-type asteroids that originated from a broad (5-20 AU-wide) region of the disk. As the disk starts to dissipate, scattered planetesimals reach sufficiently eccentric orbits to cross the terrestrial planet region and deliver water to the growing Earth. This mechanism does not depend strongly on the giant planets' orbital migration history and is generic: whenever a giant planet forms it invariably pollutes its inner planetary system with water-rich bodies.

"Hot super-Earths" (or "Mini-Neptunes") between 1 and 4 times Earth's size with period shorter than 100 days orbit 30-50\% of Sun-like type stars. Their orbital configuration -- measured as the period ratio distribution of adjacent planets in multi-planet systems -- is a strong constraint for formation models. Here we use N-body simulations with synthetic forces from an underlying evolving gaseous disk to model the formation and long-term dynamical evolution of super-Earth systems. While the gas disk is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disk. These resonant chains are far more compact than the observed super-Earth systems. Once the gas dissipates resonant chains may become dynamically unstable. They undergo a phase of giant impacts that spreads the systems out. Disk turbulence has no measurable effect on the outcome. Our simulations match observations if a small fraction of resonant chains remain stable, while most super-Earths undergo a late dynamical instability. Our statistical analysis restricts the contribution of stable systems to less than $25\%$. Our results also suggest that the large fraction of observed single planet systems does not necessarily imply any dichotomy in the architecture of planetary systems. Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in $\sim 5\%$ of the cases. This leads to a mystery: in our simulations only 50-60\% of resonant chains became unstable whereas at least 75\% (and probably 90-95\%) must be unstable to match observations.

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf.

Planets are thought to form via accretion from a remnant disk of gas and solids around a newly formed star. During this process material in the disk either remains bound to the star as part of either a planet, a smaller celestial body, or makes up part of the interplanetary medium; falls into the star; or is ejected from the system. Herein we use dynamical models to probe the abundance and properties of ejected material during late-stage planet formation and estimate their contribution to the free-floating planet population. We present 300 N-body simulations of terrestrial planet formation around a solar-type star, with and without giant planets present, using a model that accounts for collisional fragmentation. In simulations with Jupiter and Saturn analogs, about one-third of the initial (~5 Mearth) disk mass is ejected, about half in planets more massive than Mercury but less than 0.3 Mearth, and the remainder in smaller bodies. Most ejections occur within 25 Myr, which is shorter than the timescale typically required for Earth-mass planets to grow (30-100 Myr). When giant planets are omitted from our simulations, almost no material is ejected within 200 Myr and only about 1% of the initial disk is ejected by 2 Gyr. We show that about 2.5 terrestrial-mass planets are ejected per star in the Galaxy. We predict that the space-borne microlensing search for free-floating planets from the Wide-Field Infra-Red Space Telescope (WFIRST) will discover up to 15 Mars-mass planets, but few free-floating Earth-mass planets.

CoRoT-9b is one of the rare long-period (P=95.3 days) transiting giant planets with a measured mass known to date. We present a new analysis of the CoRoT-9 system based on five years of radial-velocity (RV) monitoring with HARPS and three new space-based transits observed with CoRoT and Spitzer. Combining our new data with already published measurements we redetermine the CoRoT-9 system parameters and find good agreement with the published values. We uncover a higher significance for the small but non-zero eccentricity of CoRoT-9b ($e=0.133^{+0.042}_{-0.037}$) and find no evidence for additional planets in the system. We use simulations of planet-planet scattering to show that the eccentricity of CoRoT-9b may have been generated by an instability in which a $\sim 50~M_\oplus$ planet was ejected from the system. This scattering would not have produced a spin-orbit misalignment, so we predict that CoRoT-9b orbit should lie within a few degrees of the initial plane of the protoplanetary disk. As a consequence, any significant stellar obliquity would indicate that the disk was primordially tilted.

One focus of modern astronomy is to detect temperate terrestrial exoplanets well-suited for atmospheric characterisation. A milestone was recently achieved with the detection of three Earth-sized planets transiting (i.e. passing in front of) a star just 8% the mass of the Sun 12 parsecs away. Indeed, the transiting configuration of these planets with the Jupiter-like size of their host star - named TRAPPIST-1 - makes possible in-depth studies of their atmospheric properties with current and future astronomical facilities. Here we report the results of an intensive photometric monitoring campaign of that star from the ground and with the Spitzer Space Telescope. Our observations reveal that at least seven planets with sizes and masses similar to the Earth revolve around TRAPPIST-1. The six inner planets form a near-resonant chain such that their orbital periods (1.51, 2.42, 4.04, 6.06, 9.21, 12.35 days) are near ratios of small integers. This architecture suggests that the planets formed farther from the star and migrated inward. The seven planets have equilibrium temperatures low enough to make possible liquid water on their surfaces.

Over the past two decades, large strides have been made in the field of planet formation. Yet fundamental questions remain. Here we review our state of understanding of five fundamental bottlenecks in planet formation. These are: 1) the structure and evolution of protoplanetary disks; 2) the growth of the first planetesimals; 3) orbital migration driven by interactions between proto-planets and gaseous disk; 4) the origin of the Solar System's orbital architecture; and 5) the relationship between observed super-Earths and our own terrestrial planets. Given our lack of understanding of these issues, even the most successful formation models remain on shaky ground.

Ultracool dwarfs (UCD; $T_{\rm eff}<\sim3000~$K) cool to settle on the main sequence after $\sim$1 Gyr. For brown dwarfs, this cooling never stops. Their habitable zone (HZ) thus sweeps inward at least during the first Gyr of their lives. Assuming they possess water, planets found in the HZ of UCDs have experienced a runaway greenhouse phase too hot for liquid water prior to entering the HZ. It has been proposed that such planets are desiccated by this hot early phase and enter the HZ as dry worlds. Here we model the water loss during this pre-HZ hot phase taking into account recent upper limits on the XUV emission of UCDs and using 1D radiation-hydrodynamic simulations. We address the whole range of UCDs but also focus on the planets recently found around the $0.08~M_\odot$ dwarf TRAPPIST-1. Despite assumptions maximizing the FUV-photolysis of water and the XUV-driven escape of hydrogen, we find that planets can retain significant amounts of water in the HZ of UCDs, with a sweet spot in the $0.04$-$0.06~M_\odot$ range. We also studied the TRAPPIST-1 system using observed constraints on the XUV-flux. We find that TRAPPIST-1b and c may have lost as much as 15 Earth Oceans and planet d -- which might be inside the HZ -- may have lost less than 1 Earth Ocean. Depending on their initial water contents, they could have enough water to remain habitable. TRAPPIST-1 planets are key targets for atmospheric characterization and could provide strong constraints on the water erosion around UCDs.

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

Radial velocity monitoring has found the signature of a $M \sin i = 1.3$~M$_\oplus$ planet located within the Habitable Zone (HZ) of Proxima Centauri \citep{Anglada16}. Despite a hotter past and an active host star the planet Proxima~b could have retained enough volatiles to sustain surface habitability \citep{Ribas2016}. Here we use a 3D Global Climate Model (GCM) to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect. We find that a broad range of atmospheric compositions allow surface liquid water. On a tidally-locked planet with sufficient surface water inventory, liquid water is always present, at least in the substellar region. With a non-synchronous rotation, this requires a minimum greenhouse warming ($\sim$10~mbar of CO$_2$ and 1~bar of N$_2$). If the planet is dryer, $\sim$0.5~bar/1.5~bars of CO$_2$ (respectively for asynchronous/synchronous rotation) suffice to prevent the trapping of any arbitrary small water inventory into polar/nightside ice caps. We produce reflection/emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible by direct imaging with forthcoming large telescopes. The angular separation of $7 \lambda/D$ at 1~$\mu$m (with the E-ELT) and a contrast of $\sim$10$^{-7}$ will enable high-resolution spectroscopy and the search for molecular signatures, including H$_2$O, O$_2$, and CO$_2$. The observation of thermal phase curves can be attempted with JWST, thanks to a contrast of $2\times10^{-5}$ at 10~$\mu$m. Proxima~b will also be an exceptional target for future IR interferometers. Within a decade it will be possible to image Proxima~b and possibly determine whether this exoplanet's surface is habitable.

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

We present the results of hydrodynamical simulations of the orbital evolution of planets undergoing runaway gas accretion in radiative discs. We consider accreting disc models with constant mass flux through the disc, and where radiative cooling balances the effect of viscous heating and stellar irradiation. We assume that 20-30 $M_\oplus$ giant planet cores are formed in the region where viscous heating dominates and migrate outward under the action of a strong corotation torque. In the case where gas accretion is neglected, we find evidence for strong dynamical torques in accreting discs with accretion rates ${\dot M}\gtrsim 7\times 10^{-8} \;M_\odot/yr$. Their main effect is to increase outward migration rates by a factor of $\sim 2$ typically. In the presence of gas accretion, however, runaway outward migration is observed with the planet passing through the zero-torque radius and the transition between the viscous heating and stellar heating dominated regimes. The ability for an accreting planet to enter a fast migration regime is found to depend strongly on the planet growth rate, but can occur for values of the mass flux through the disc of ${\dot M}\gtrsim 5\times 10^{-8} \;M_\odot/yr$. We find that an episode of runaway outward migration can cause an accreting planet formed in the 5-10 AU region to temporarily orbit at star-planet separations as large as $\sim$60-70 AU. However, increase in the amplitude of the Lindblad torque associated with planet growth plus change in the streamline topology near the planet systematically cause the direction of migration to be reversed. Our results indicate that a planet can reach large orbital distances under the combined effect of dynamical torques and gas accretion, but an alternative mechanism is required to explain the presence of massive planets on wide orbits.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

We investigate the prospects for the capture of the proposed Planet 9 from other stars in the Sun's birth cluster. Any capture scenario must satisfy three conditions: the encounter must be more distant than ~150 au to avoid perturbing the Kuiper belt; the other star must have a wide-orbit planet (a>~100au); the planet must be captured onto an appropriate orbit to sculpt the orbital distribution of wide-orbit Solar System bodies. Here we use N-body simulations to show that these criteria may be simultaneously satisfied. In a few percent of slow close encounters in a cluster, bodies are captured onto heliocentric, Planet 9-like orbits. During the ~100 Myr cluster phase, many stars are likely to host planets on highly-eccentric orbits with apastron distances beyond 100 au if Neptune-sized planets are common and susceptible to planet--planet scattering. While the existence of Planet 9 remains unproven, we consider capture from one of the Sun's young brethren a plausible route to explain such an object's orbit. Capture appears to predict a large population of Trans-Neptunian Objects (TNOs) whose orbits are aligned with the captured planet, and we propose that different formation mechanisms will be distinguishable based on their imprint on the distribution of TNOs.

Jupiter's core is generally assumed to have formed beyond the snow line. Here we consider an alternative scenario, that Jupiter's core may have accumulated in the innermost parts of the protoplanetary disk. A growing body of research suggests that small particles ("pebbles") continually drift inward through the disk. If a fraction of drifting pebbles is trapped at the inner edge of the disk a several Earth-mass core can quickly grow. Subsequently, the core may migrate outward beyond the snow line via planet-disk interactions. Of course, to reach the outer Solar System Jupiter's core must traverse the terrestrial planet-forming region. We use N-body simulations including synthetic forces from an underlying gaseous disk to study how the outward migration of Jupiter's core sculpts the terrestrial zone. If the outward migration is fast (Tmig~10^4 years), the core simply migrates past resident planetesimals and planetary embryos. However, if its migration is slower (Tmig~10^5 years) the core removes solids from the inner disk by shepherding objects in mean motion resonances. In many cases the disk interior to 0.5-1 AU is cleared of embryos and most planetesimals. By generating a mass deficit close to the Sun, the outward migration of Jupiter's core may thus explain the absence of terrestrial planets closer than Mercury. Jupiter's migrating core often stimulates the growth of another large (~Earth-mass) core -- that may provide a seed for Saturn's core -- trapped in exterior resonance. The migrating core also may transport a fraction of terrestrial planetesimals, such as the putative parent bodies of iron meteorites, to the asteroid belt.

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass $0.1~M_\odot$, Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water.

Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disk between 1 and 3 AU. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 AU and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas-drag, as has been proposed to explain the formation of close-in super Earths. Here we show that the asteroid belt's orbital excitation provides a crucial constraint against this scenario for the Solar System. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from disks of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Disks with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogs significantly more massive than the real planet. In contrast, a disk with a surface density gradient proportional to $r^{-5.5}$ reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disk profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.

Reproducing Uranus and Neptune remains a challenge for simulations of solar system formation. The ice giants' peculiar obliquities suggest that they both suffered giant collisions during their formation. Thus, there must have been an epoch of accretion dominated by collisions among large planetary embryos in the primordial outer solar system. We test this idea using N-body numerical simulations including the effects of a gaseous protoplanetary disk. One strong constraint is that the masses of the ice giants are very similar -- the Neptune/Uranus mass ratio is $\sim1.18$. We show that similar-size ice giants do indeed form by collisions between planetary embryos beyond Saturn. The fraction of successful simulations varies depending on the initial number of planetary embryos in the system, their individual and total masses. Similar-sized ice giants are consistently reproduced in simulations starting with 5-10 planetary embryos with initial masses of $\sim$3-6 ${\rm M_\oplus}$. We conclude that accretion from a population of planetary embryos is a plausible scenario for the origin of Uranus and Neptune.

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

Most of the properties of the Earth-Moon system can be explained by a collision between a planetary embryo and the growing Earth late in the accretion process. Simulations show that most of the material that eventually aggregates to form the Moon originates from the impactor. However, analysis of the terrestrial and lunar isotopic composition show them to be highly similar. In contrast, the compositions of other solar system bodies are significantly different than the Earth and Moon. This poses a major challenge to the giant impact scenario since the Moon-forming impactor is then thought to also have differed in composition from the proto-Earth. Here we track the feeding zones of growing planets in a suite of simulations of planetary accretion, in order to measure the composition of Moon-forming impactors. We find that different planets formed in the same simulation have distinct compositions, but the compositions of giant impactors are systematically more similar to the planets they impact. A significant fraction of planet-impactor pairs have virtually identical compositions. Thus, the similarity in composition between the Earth and Moon could be a natural consequence of a late giant impact.

Planets of 1-4 times Earth's size on orbits shorter than 100 days exist around 30-50% of all Sun-like stars. In fact, the Solar System is particularly outstanding in its lack of "hot super-Earths" (or "mini-Neptunes"). These planets -- or their building blocks -- may have formed on wider orbits and migrated inward due to interactions with the gaseous protoplanetary disk. Here, we use a suite of dynamical simulations to show that gas giant planets act as barriers to the inward migration of super-Earths initially placed on more distant orbits. Jupiter's early formation may have prevented Uranus and Neptune (and perhaps Saturn's core) from becoming hot super-Earths. Our model predicts that the populations of hot super-Earth systems and Jupiter-like planets should be anti-correlated: gas giants (especially if they form early) should be rare in systems with many hot super-Earths. Testing this prediction will constitute a crucial assessment of the validity of the migration hypothesis for the origin of close-in super-Earths.

We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.

Planetary embryos embedded in gaseous protoplanetary disks undergo Type I orbital migration. Migration can be inward or outward depending on the local disk properties but, in general, only planets more massive than several $M_\oplus$ can migrate outward. Here we propose that an embryo's migration history determines whether it becomes a hot super-Earth or the core of a giant planet. Systems of hot super-Earths (or mini-Neptunes) form when embryos migrate inward and pile up at the inner edge of the disk. Giant planet cores form when inward-migrating embryos become massive enough to switch direction and migrate outward. We present simulations of this process using a modified N-body code, starting from a swarm of planetary embryos. Systems of hot super-Earths form in resonant chains with the innermost planet at or interior to the disk inner edge. Resonant chains are disrupted by late dynamical instabilities triggered by the dispersal of the gaseous disk. Giant planet cores migrate outward toward zero-torque zones, which move inward and eventually disappear as the disk disperses. Giant planet cores migrate inward with these zones and are stranded at ~1-5 AU. Our model reproduces several properties of the observed extra-solar planet populations. The frequency of giant planet cores increases strongly when the mass in solids is increased, consistent with the observed giant exoplanet - stellar metallicity correlation. The frequency of hot super-Earths is not a function of stellar metallicity, also in agreement with observations. Our simulations can reproduce the broad characteristics of the observed super-Earth population.

Our knowledge of planets' orbital dynamics, which was based on Solar System studies, has been challenged by the diversity of exoplanetary systems. Around cool and ultra cool dwarfs, the influence of tides on the orbital and spin evolution of planets can strongly affect their climate and their capacity to host surface liquid water. We illustrate the role of tides and dynamics with the extreme case of planets orbiting around brown dwarfs. In multiple planet systems, the eccentricity is excited by planet-planet interactions. Planets are therefore heated up from the inside by the tidally-induced friction. This process can heat a habitable zone planet to such a level that surface liquid water cannot exist. We also talk about the newly discovered potentially habitable Earth-sized planet Kepler-186f. Given the poorly estimated age of the system, the planet could still be evolving towards synchronization and have a high obliquity or be pseudo-synchronized with a zero obliquity. These two configurations would have a different effect on the climate of this planet.

The `Grand Tack' model proposes that the inner Solar System was sculpted by the giant planets' orbital migration in the gaseous protoplanetary disk. Jupiter first migrated inward then Jupiter and Saturn migrated back outward together. If Jupiter's turnaround or "tack" point was at ~1.5 AU the inner disk of terrestrial building blocks would have been truncated at ~1 AU, naturally producing the terrestrial planets' masses and spacing. During the gas giants' migration the asteroid belt is severely depleted but repopulated by distinct planetesimal reservoirs that can be associated with the present-day S and C types. The giant planets' orbits are consistent with the later evolution of the outer Solar System. Here we confront common criticisms of the Grand Tack model. We show that some uncertainties remain regarding the Tack mechanism itself; the most critical unknown is the timing and rate of gas accretion onto Saturn and Jupiter. Current isotopic and compositional measurements of Solar System bodies -- including the D/H ratios of Saturn's satellites -- do not refute the model. We discuss how alternate models for the formation of the terrestrial planets each suffer from an internal inconsistency and/or place a strong and very specific requirement on the properties of the protoplanetary disk. We conclude that the Grand Tack model remains viable and consistent with our current understanding of planet formation. Nonetheless, we encourage additional tests of the Grand Tack as well as the construction of alternate models.

A large number of observed exoplanets are part of multiple planet systems. Most of these systems are sufficiently close-in to be tidally evolving. In such systems, there is a competition between the excitation caused by planet-planet interactions and tidal damping. Using as an example two multiple planet systems, which host planets in the surface liquid water habitable zone (HZ): Kepler-62 and Kepler-186, we show the importance and effect of both planetary and stellar tides on the dynamical evolution of planets and on the climate of the HZ planets.

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

Super-Earths with orbital periods less than 100 days are extremely abundant around Sun-like stars. It is unlikely that these planets formed at their current locations. Rather, they likely formed at large distances from the star and subsequently migrated inward. Here we use N-body simulations to study the effect of super-Earths on the accretion of rocky planets. In our simulations, one or more super-Earths migrates inward through a disk of planetary embryos and planetesimals embedded in a gaseous disk. We tested a wide range of migration speeds and configurations. Fast-migrating super-Earths ($\tau_{mig} \sim$0.01-0.1 Myr) only have a modest effect on the protoplanetary embryos and planetesimals. Sufficient material survives to form rocky, Earth-like planets on orbits exterior to the super-Earths'. In contrast, slowly migrating super-Earths shepherd rocky material interior to their orbits and strongly deplete the terrestrial planet-forming zone. In this situation any Earth-sized planets in the habitable zone are extremely volatile-rich and are therefore probably not Earth-like.

The Kepler-186 system consists of five planets orbiting an early-M dwarf. The planets have physical radii of 1.0-1.50 R$_\oplus$ and orbital periods of 4 to 130 days. The $1.1~$R$_\oplus$ Kepler-186f with a period of 130 days is of particular interest. Its insolation of roughly $0.32~S_\odot$places it within the liquid water habitable zone. We present a multi-faceted study of the Kepler-186 system. First, we show that the distribution of planet masses can be roughly reproduced if the planets accreted from a high-surface density disk presumably sculpted by an earlier phase of migration. However, our simulations predict the existence of 1-2 undetected planets between planets e and f. Next, we present a dynamical analysis of the system including the effect of tides. The timescale for tidal evolution is short enough that the four inner planets must have small obliquities and near-synchronous rotation rates. Tidal evolution of Kepler-186f is slow enough that its current spin state depends on a combination of its dissipation rate and the stellar age. Finally, we study the habitability of Kepler-186f with a 1-D climate model. The planet's surface temperature can be raised above 273 K with 0.5-5 bars of CO$_2$, depending on the amount of N$_2$ present. Kepler-186f represents a case study of an Earth-sized planet in the cooler regions of the habitable zone of a cool star.

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

We present Rossiter-McLaughlin observations of the transiting super-Earth 55 Cnc e collected during six transit events between January 2012 and November 2013 with HARPS and HARPS-N. We detect no radial-velocity signal above 35 cm/s (3-sigma) and confine the stellar v sin i to 0.2 +/- 0.5 km/s. The star appears to be a very slow rotator, producing a very low amplitude Rossiter-McLaughlin effect. Given such a low amplitude, the Rossiter-McLaughlin effect of 55 Cnc e is undetected in our data, and any spin-orbit angle of the system remains possible. We also performed Doppler tomography and reach a similar conclusion. Our results offer a glimpse of the capacity of future instrumentation to study low amplitude Rossiter-McLaughlin effects produced by super-Earths.

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

Planets close to their stars are thought to form farther out and migrate inward due to angular momentum exchange with gaseous protoplanetary disks. This process can produce systems of planets in co-orbital (Trojan or 1:1) resonance, in which two planets share the same orbit, usually separated by 60 degrees. Co-orbital systems are detectable among the planetary systems found by the Kepler mission either directly or by transit timing variations. However, no co-orbital systems have been found within the thousands of Kepler planets and candidates. Here we study the orbital evolution of co-orbital planets embedded in a protoplanetary disk using a grid-based hydrodynamics code. We show that pairs of similar-mass planets in co-orbital resonance are disrupted during large-scale orbital migration. Destabilization occurs when one or both planets is near the critical mass needed to open a gap in the gaseous disk. A confined gap is opened that spans the 60 degree azimuthal separation between planets. This alters the torques imparted by the disk on each planet -- pushing the leading planet outward and the trailing planet inward -- and disrupts the resonance. The mechanism applies to systems in which the two planets' masses differ by a factor of two or less. In a simple flared disk model the critical mass for gap opening varies from a few Earth masses at the inner edge of the disk to 1 Saturn-mass at 5 AU. A pair of co-orbital planets with masses in this range that migrates will enter a region where the planets are at the gap-opening limit. At that point the resonance is disrupted. We therefore predict an absence of planets on co-orbital configurations with masses in the super-Earth to Saturn mass range with similar masses.

Vega has been shown to host multiple dust populations, including both hot exo-zodiacal dust at sub-AU radii and a cold debris disk extending beyond 100 AU. We use dynamical simulations to show how Vega's hot dust can be created by long-range gravitational scattering of planetesimals from its cold outer regions. Planetesimals are scattered progressively inward by a system of 5-7 planets from 30-60 AU to very close-in. In successful simulations the outermost planets are typically Neptune-mass. The back-reaction of planetesimal scattering causes these planets to migrate outward and continually interact with fresh planetesimals, replenishing the source of scattered bodies. The most favorable cases for producing Vega's exo-zodi have negative radial mass gradients, with sub-Saturn- to Jupiter-mass inner planets at 5-10 AU and outer planets of 2.5 to 20 Earth masses. The mechanism fails if a Jupiter-sized planet exists beyond ~15 AU because the planet preferentially ejects planetesimals before they can reach the inner system. Direct-imaging planet searches can therefore directly test this mechanism.

The quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-sized - that transit a 0.47+\-0.05 Rsun star. The intensity and spectrum of the star's radiation places Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and H2O at its surface, then some of this H2O is likely to be in liquid form.

High levels of exozodiacal dust have been observed in the inner regions of a large fraction of main sequence stars. Given the short lifetime of the observed small dust grains, these 'exozodis' are difficult to explain, especially for old (>100 Myr) stars. The exozodiacal dust may be observed as excess emission in the mid-infrared, or using interferometry. We hypothesise that exozodi are sustained by planetesimals scattered by planets inwards from an outer planetesimal belt, where collision timescales are long. In this work, we use N-body simulations to show that the outwards migration of a planet into a belt, driven by the scattering of planetesimals, can increase, or sustain, the rate at which planetesimals are scattered from the outer belt to the exozodi region. We hypothesise that this increase is sufficient to sustain the observed exozodi on Gyr timescales. No correlation between observations of an outer belt and an exozodi is required for this scenario to work, as the outer belt may be too faint to detect. If planetesimal driven migration does explain the observed exozodi, this work suggests that the presence of an exozodi indicates the presence of outer planets and a planetesimal belt.

It has been proposed that the observed systems of hot super-Earths formed in situ from high-mass disks. By fitting a disk profile to the entire population of Kepler planet candidates, Chiang & Laughlin (2013) constructed a "minimum-mass extrasolar nebula" with surface density profile Sigma r^-1.6. Here we use multiple-planet systems to show that it is inconsistent to assume a universal disk profile. Systems with 3-6 low-mass planets (or planet candidates) produce a diversity of minimum-mass disks with surface density profiles ranging from Sigma r^-3.2 to Sigma r^0.5 (5th-95th percentile). By simulating the transit detection of populations of synthetic planetary systems designed to match the properties of observed super-Earth systems, we show that a universal disk profile is statistically excluded at high confidence. Rather, the underlying distribution of minimum-mass disks is characterized by a broad range of surface density slopes. Models of gaseous disks can only explain a narrow range of slopes (roughly between r^0 and r^-1.5). Yet accretion of terrestrial planets in a gas-free environment preserves the initial radial distribution of building blocks. The known systems of hot super-Earths must therefore not represent the structure of their parent gas disks and can not have predominantly formed in situ. We instead interpret the diversity of disk slopes as the imprint of a process that re-arranged the solids relative to the gas in the inner parts of protoplanetary disks. A plausible mechanism is inward type 1 migration of Mars- to Earth-mass planetary embryos, perhaps followed by a final assembly phase.

Since exoplanets were detected using the radial velocity method, they have revealed a diverse distribution of orbital configurations. Amongst these are planets in highly eccentric orbits (e > 0.5). Most of these systems consist of a single planet but several have been found to also contain a longer period planet in a near-circular orbit. Here we use the latest Keplerian orbital solutions to investigate four known systems which exhibit this extreme eccentricity diversity; HD 37605, HD 74156, HD 163607, and HD 168443. We place limits on the presence of additional planets in these systems based on the radial velocity residuals. We show that the two known planets in each system exchange angular momentum through secular oscillations of their eccentricities. We calculate the amplitude and timescale for these eccentricity oscillations and associated periastron precession. We further demonstrate the effect of mutual orbital inclinations on the amplitude of high-frequency eccentricity oscillations. Finally, we discuss the implications of these oscillations in the context of possible origin scenarios for unequal eccentricities.

We present numerical simulations modeling the orbital evolution of very wide binaries, pairs of stars separated by over ~1000 AU. Due to perturbations from other passing stars and the Milky Way's tide, the orbits of very wide binary stars occasionally become extremely eccentric, which forces close encounters between the companion stars (Kaib et al. 2013). We show that this process causes a stellar collision between very wide binary companion stars once every 1000-7500 years on average in the Milky Way. One of the main uncertainties in this collision rate is the amount of energy dissipated by dynamic tides during close (but not collisional) periastron passages. This dissipation presents a dynamical barrier to stellar collisions and can instead transform very wide binaries into close or contact binaries. However, for any plausible tidal dissipation model, very wide binary stars are an unrealized, and potentially the dominant, source of stellar collisions in our Galaxy. Such collisions should occur throughout the thin disk of the Milky Way. Stellar collisions within very wide binaries should yield a small population of single, Li-depleted, rapidly rotating massive stars.

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

The field of exoplanetary science has experienced a recent surge of new systems that is largely due to the precision photometry provided by the Kepler mission. The latest discoveries have included compact planetary systems in which the orbits of the planets all lie relatively close to the host star, which presents interesting challenges in terms of formation and dynamical evolution. The compact exoplanetary systems are analogous to the moons orbiting the giant planets in our Solar System, in terms of their relative sizes and semi-major axes. We present a study that quantifies the scaled sizes and separations of the Solar System moons with respect to their hosts. We perform a similar study for a large sample of confirmed Kepler planets in multi-planet systems. We show that a comparison between the two samples leads to a similar correlation between their scaled sizes and separation distributions. The different gradients of the correlations may be indicative of differences in the formation and/or long-term dynamics of moon and planetary systems.

Earth-mass bodies are expected to undergo Type I migration directed either inward or outward depending on the thermodynamical state of the protoplanetary disc. Zones of convergent migration exist where the Type I torque cancels out. We study the evolution of multiple protoplanets of a few Earth masses embedded in a non-isothermal protoplanetary disc. The protoplanets are located in the vicinity of a convergence zone located at the transition between two different opacity regimes. Inside the convergence zone, Type I migration is directed outward and outside the zone migration is directed inward. We used a grid-based hydrodynamical code that includes radiative effects. We performed simulations varying the initial number of embryos and tested the effect of including stochastic forces to mimic the effects resulting from turbulence. We also performed N-body runs calibrated on hydrodynamical calculations to follow the evolution on Myr timescales. For a small number of initial embryos (N = 5-7) and in the absence of stochastic forcing, the population of protoplanets migrates convergently toward the zero-torque radius and forms a stable resonant chain that protects embryos from close encounters. In systems with a larger initial number of embryos, or in which stochastic forces were included, these resonant configurations are disrupted. This in turn leads to the growth of larger cores via a phase of giant impacts, after which the system settles to a new stable resonant configuration. Giant planets cores with masses of 10 Earth masses formed in about half of the simulations with initial protoplanet masses of m_p = 3 Earth masses but in only 15% of simulations with m_p = 1 Earth mass. This suggests that if ~2-3 Earth mass protoplanets can form in less than ~1 Myr, convergent migration and giant collisions can grow giant planet cores at Type I migration convergence zones.

If mutual gravitational scattering among exoplanets occurs, then it may produce unique orbital properties. For example, two-planet systems that lie near the boundary between circulation and libration of their periapses could result if planet-planet scattering ejected a former third planet quickly, leaving one planet on an eccentric orbit and the other on a circular orbit. We first improve upon previous work that examined the apsidal behavior of known multiplanet systems by doubling the sample size and including observational uncertainties. This analysis recovers previous results that demonstrated that many systems lay on the apsidal boundary between libration and circulation. We then performed over 12,000 three-dimensional N-body simulations of hypothetical three-body systems that are unstable, but stabilize to two-body systems after an ejection. Using these synthetic two-planet systems, we test the planet-planet scattering hypothesis by comparing their apsidal behavior, over a range of viewing angles, to that of the observed systems and find that they are statistically consistent regardless of the multiplicity of the observed systems. Finally, we combine our results with previous studies to show that, from the sampled cases, the most likely planetary mass function prior to planet-planet scattering follows a power law with index -1.1. We find that this pre-scattering mass function predicts a mutual inclination frequency distribution that follows an exponential function with an index between -0.06 and -0.1.

Using N-body simulations with planet-disk interactions, we present a mechanism capable of forming compact systems of hot super Earths such as Kepler 11. Recent studies show that outward migration is common in the inner parts of radiative disks. However we show that two processes naturally tip the balance in favor of inward migration. First the corotation torque is too weak to generate outward migration for planetary embryos less massive than $4\mearth$. Second, system of multiple embryos generate sustained non-zero eccentricities that damp the corotation torque and again favor inward migration. Migration and accretion of planetary embryos in realistic disks naturally produce super Earths in resonant chains near the disk inner edge. Their compact configuration is similar to the observed systems.

Given their tendency to be incorporated into the core during differentiation, the highly-siderophile elements (HSEs) in Earth's mantle are thought to have been accreted as a `late veneer' after the end of the giant impact phase. Bottke et al (2010) proposed that the large Earth-to-Moon HSE abundance ratio can be explained if the late veneer was characterized by large (D = 1000-4000km) impactors. Here we simulate the evolution of the terrestrial planets during a stochastic late veneer phase from the end of accretion until the start of the late heavy bombardment ~500 Myr later. We show that a late veneer population of 0.05 Earth masses dominated by large (D > 1000km) bodies naturally delivers a ~0.01 Earth mass veneer to Earth, consistent with constraints. The eccentricities and inclinations of the terrestrial planets are excited by close encounters with the largest late veneer bodies. We find the best agreement with their post-veneer orbits if either a) the terrestrial planets' pre-veneer angular momentum deficit AMD_0 was less than half of the current one AMD_now, or b) AMD_0 <= AMD_now and the veneer was limited to smaller (D_max <= 2000km) bodies. Veneer impacts on Venus, Earth and Mars were mostly accretionary but on Mercury and the Moon they were mostly erosive. In ~20% of simulations an energetic impact occurred that could have removed >25% of Mercury's mass, thereby increasing its iron mass fraction. We show that, due to the erosive nature of larger impacts, the Moon cannot accrete any material from objects larger than 500-1000km. The large Earth-to-Moon HSE abundance ratio is explained if the late veneer included large impactors (D >= 500-1000km) regardless of their size distribution, as long as most of Earth's veneer came from large bodies. The spin angular momenta imparted by stochastic late veneer impacts was far in excess of the current ones for Mercury and Venus.

Excess emission, associated with warm, dust belts, commonly known as exozodis, has been observed around a third of nearby stars. The high levels of dust required to explain the observations are not generally consistent with steady-state evolution. A common suggestion is that the dust results from the aftermath of a dynamical instability, an event akin to the Solar System's Late Heavy Bombardment. In this work, we use a database of N-body simulations to investigate the aftermath of dynamical instabilities between giant planets in systems with outer planetesimal belts. We find that, whilst there is a significant increase in the mass of material scattered into the inner regions of the planetary system following an instability, this is a short-lived effect. Using the maximum lifetime of this material, we determine that even if every star has a planetary system that goes unstable, there is a very low probability that we observe more than a maximum of 1% of sun-like stars in the aftermath of an instability, and that the fraction of planetary systems currently in the aftermath of an instability is more likely to be limited to <0.06. This probability increases marginally for younger or higher mass stars. We conclude that the production of warm dust in the aftermath of dynamical instabilities is too short-lived to be the dominant source of the abundantly observed exozodiacal dust.

Short period (<50 days) low-mass (<10Mearth) exoplanets are abundant and the few of them whose radius and mass have been measured already reveal a diversity in composition. Some of these exoplanets are found on eccentric orbits and are subjected to strong tides affecting their rotation and resulting in significant tidal heating. Within this population, some planets are likely to be depleted in volatiles and have no atmosphere. We model the thermal emission of these "Super Mercuries" to study the signatures of rotation and tidal dissipation on their infrared light curve. We compute the time-dependent temperature map at the surface and in the subsurface of the planet and the resulting disk-integrated emission spectrum received by a distant observer for any observation geometry. We calculate the illumination of the planetary surface for any Keplerian orbit and rotation. We include the internal tidal heat flow, vertical heat diffusion in the subsurface and generate synthetic light curves. We show that the different rotation periods predicted by tidal models (spin-orbit resonances, pseudo-synchronization) produce different photometric signatures, which are observable provided that the thermal inertia of the surface is high, like that of solid or melted rocks (but not regolith). Tidal dissipation can also directly affect the light curves and make the inference of the rotation more difficult or easier depending on the existence of hot spots on the surface. Infrared light curve measurement with the James Webb Space Telescope and EChO can be used to infer exoplanets' rotation periods and dissipation rates and thus to test tidal models. This data will also constrain the nature of the (sub)surface by constraining the thermal inertia.

Given the very close proximity of their habitable zones, brown dwarfs represent high-value targets in the search for nearby transiting habitable planets that may be suitable for follow-up occultation spectroscopy. In this paper we develop search strategies to find habitable planets transiting brown dwarfs depending on their maximum habitable orbital period (PHZ out). Habitable planets with PHZ out shorter than the useful duration of a night (e.g. 8-10 hrs) can be screened with 100 percent completeness from a single location and in a single night (near-IR). More luminous brown dwarfs require continuous monitoring for longer duration, e.g. from space or from a longitude-distributed network (one test scheduling achieved - 3 telescopes, 13.5 contiguous hours). Using a simulated survey of the 21 closest known brown dwarfs (within 7 pc) we find that the probability of detecting at least one transiting habitable planet is between 4.5 +5.6-1.4 and 56 +31-13 percent, depending on our assumptions. We calculate that brown dwarfs within 5-10 pc are characterizable for potential biosignatures with a 6.5 m space telescope using approx. 1 percent of a 5-year mission's life-time spread over a contiguous segment only one fifth to one tenth of this duration.

Exoplanetary science has reached a historic moment. The James Webb Space Telescope will be capable of probing the atmospheres of rocky planets, and perhaps even search for biologically produced gases. However this is contingent on identifying suitable targets before the end of the mission. A race therefore, is on, to find transiting planets with the most favorable properties, in time for the launch. Here, we describe a realistic opportunity to discover extremely favorable targets - rocky planets transiting nearby brown dwarfs - using the Spitzer Space Telescope as a survey instrument. Harnessing the continuous time coverage and the exquisite precision of Spitzer in a 5,400 hour campaign monitoring nearby brown dwarfs, we will detect a handful of planetary systems with planets as small as Mars. The survey we envision is a logical extension of the immense progress that has been realized in the field of exoplanets and a natural outcome of the exploration of the solar neighborhood to map where the nearest habitable rocky planets are located (as advocated by the 2010 Decadal Survey). Our program represents an essential step towards the atmospheric characterization of terrestrial planets and carries the compelling promise of studying the concept of habitability beyond Earth-like conditions. In addition, our photometric monitoring will provide invaluable observations of a large sample of nearby brown dwarfs situated close to the M/L transition. This is why, we also advocate an immediate public release of the survey data, to guarantee rapid progress on the planet search and provide a treasure trove of data for brown dwarf science.

We use the super Earth 55 Cnc e as a case study to address an observable effect of tidal heating. We investigate whether planet-planet interactions can force the eccentricity of this planet to a level affecting the eclipse depth observed with Spitzer. Using the constant time lag tidal model, we first calculate the observed planet flux as a function of albedo and eccentricity, for different tidal dissipation constants and for two extreme cases: a planet with no heat redistribution and a planet with full heat redistribution. We derive the values of albedo and eccentricity that match the observed transit depth. We then perform N-body simulations of the planetary system including tides and General Relativity to follow the evolution of the eccentricity of planet e. We compare the range of eccentricities given by the simulations with the eccentricities required to alter the eclipse depth. We find that the eccentricity of planet e can be large enough to contribute at a measurable level to the thermal emission measured with Spitzer. This affects the constraints on the albedo of the planet, which can be as high as 0.9 (instead of 0.55 when ignoring tidal heating). We also derive a maximum value for the eccentricity of planet e directly from the eclipse depth: e<0.015 assuming Earth's dissipation constant. Transiting exoplanets in multiple planet systems are exceptional targets for testing tidal models because their tidal luminosity may be observable. Future multi-wavelengths observations of eclipse depth and phase curves (for instance with EChO and JWST) should allow us to better resolve the temperature map of these planets and break the degeneracy between albedo and tidal heating that remains for single band observations. In addition, an accurate determination of the eccentricity will make it possible to constrain the dissipation rate of the planet and to probe its internal structure.

Earth-mass planets embedded in gaseous protoplanetary disks undergo Type I orbital migration. In radiative disks an additional component of the corotation torque scaling with the entropy gradient across the horseshoe region can counteract the general inward migration, Type I migration can then be directed either inward or outward depending on the local disk properties. Thus, special locations exist in the disk toward which planets migrate in a convergent way. Here we present N-body simulations of the convergent migration of systems of low-mass (M=1-10 m_earth) planets. We show that planets do not actually converge in convergence zones. Rather, they become trapped in chains of mean motion resonances (MMRs). This causes the planets' eccentricities to increase to high enough values to affect the structure of the horseshoe region and weaken the positive corotation torque. The zero-torque equilibrium point of the resonant chain of planets is determined by the sum of the attenuated corotation torques and unattenuated differential Lindblad torques acting on each planet. The effective convergence zone is shifted inward. Systems with several planets can experience stochastic migration as a whole due to continuous perturbations from planets entering and leaving resonances.

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.

Starting from planetary systems with three giant planets and an outer disk of planetesimals, we use dynamical simulations to show how dynamical instabilities can transform planetesimal disks into 100-1000 AU-scale isotropic clouds. The instabilities involve a phase of planet-planet scattering that concludes with the ejection of one or more planets and the inward-scattering of the surviving gas giant(s) to remove them from direct dynamical contact with the planetesimals. "Mini-Oort clouds" are thus formed from scattered planetesimals whose orbits are frozen by the abrupt disappearance of the perturbing giant planet. Although the planetesimal orbits are virtually isotropic, the surviving giant planets tend to have modest inclinations (typically ~10 degrees) with respect to the initial orbital plane. The collisional lifetimes of mini-Oort clouds are long (10 Myr to >10 Gyr) and there is a window of ~100 Myr or longer during which they produce spherical clouds of potentially observable dust at 70 microns. If the formation channel for hot Jupiters commonly involves planetary close encounters, we predict a correlation between this subset of extrasolar planetary systems and mini-Oort clouds.

In this chapter we review the astrophysical origins of Earth's carbon, starting from the products of the Big Bang and culminating with the Earth's formation. We review the measured compositions of different primitive objects including comets, various classes of meteorites and interstellar dust particles. We discuss the composition of the Solar Nebula, especially with regards to the distribution of volatiles such as carbon. We discuss dynamical models of planetary formation from planetesimals and planetary embryos, and the timescale for volatile delivery to the growing Earth from different sources. Finally, we review Earth's carbon reservoirs. Throughout the chapter we highlight open questions related to planet formation, meteoritics, and geochemistry.

Prior work has found that a variety of terrestrial planetary compositions are expected to occur within known extrasolar planetary systems. However, such studies ignored the effects of giant planet migration, which is thought to be very common in extra-solar systems. Here we present calculations of the compositions of terrestrial planets that formed in dynamical simulations incorporating varying degrees of giant planet migration. We used chemical equilibrium models of the solid material present in the disks of five known planetary host stars: the Sun, GJ 777, HD4203, HD19994 and HD213240. Giant planet migration has a strong effect on the compositions of simulated terrestrial planets as the migration results large-scale mixing between terrestrial planet building blocks that condensed at a range of temperatures. This mixing acts to 1) increase the typical abundance of Mg-rich silicates in the terrestrial planets feeding zones and thus increase the frequency of planets with Earth-like compositions compared with simulations with static giant planet orbits; and 2) drastically increase the efficiency of the delivery of hydrous phases (water and serpentine) to terrestrial planets and thus produce water worlds and/or wet Earths. Our results demonstrate that although a wide variety of terrestrial planet compositions can still be produced, planets with Earth-like compositions should be common within extrasolar planetary systems.

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

We investigate how the evolution of the stellar spin rate affects, and is affected by, planets in close orbits, via star-planet tidal interactions. To do this, we used a standard equilibrium tidal model to compute the orbital evolution of single planets orbiting both Sun-like stars and 0.1 M\odot M-dwarfs. We tested two stellar spin evolution profiles, one with fast initial rotation (P=1.2 day) and one with slow initial rotation (P=8 day). We tested the effect of varying the stellar and planetary dissipation and the planet's mass and initial orbital radius. Conclusions: Tidal evolution allows to differentiate the early behaviors of extremely close-in planets orbiting either a rapidly rotating star or a slowly rotating star. The early spin-up of the star allows the close-in planets around fast rotators to survive the early evolution. For planets around M-dwarfs, surviving the early evolution means surviving on Gyr timescales whereas for Sun-like stars the spin-down brings about late mergers of Jupiter planets. In light of this study, we can say that differentiating between one spin evolution from another given the present position of planets can be very tricky. Unless we can observe some markers of former evolution it is nearly impossible to distinguish the two very different spin profiles, let alone intermediate spin profiles. Though some conclusions can still be drawn from statistical distributions of planets around fully convective M-dwarfs. However, if the tidal evolution brings about a merger late in its history it can also entail a noticeable acceleration of the star in late ages, so that it is possible to have old stars that spin rapidly. This raises the question of better constraining the age of stars.

We present models for the formation of terrestrial planets, and the collisional evolution of debris disks, in planetary systems that contain multiple unstable gas giants. We previously showed that the dynamics of the giant planets introduces a correlation between the presence of terrestrial planets and debris disks. Here we present new simulations that show that this connection is qualitatively robust to changes in: the mass distribution of the giant planets, the width and mass distribution of the outer planetesimal disk, and the presence of gas in the disk. We discuss how variations in these parameters affect the evolution. Systems with equal-mass giant planets undergo the most violent instabilities, and these destroy both terrestrial planets and the outer planetesimal disks that produce debris disks. In contrast, systems with low-mass giant planets efficiently produce both terrestrial planets and debris disks. A large fraction of systems with low-mass outermost giant planets have stable gaps between these planets that are frequently populated by planetesimals. Planetesimal belts between outer giant planets may affect debris disk SEDs. If Earth-mass seeds are present in outer planetesimal disks, the disks radially spread to colder temperatures. We argue that this may explain the very low frequency of > 1 Gyr-old solar-type stars with observed 24 micron excesses. Among the (limited) set of configurations explored, the best candidates for hosting terrestrial planets at ~1 AU are stars older than 0.1-1 Gyr with bright debris disks at 70 micron but with no currently-known giant planets. These systems combine evidence for rocky building blocks, with giant planet properties least likely to undergo destructive dynamical evolution. We predict an anti-correlation between debris disks and eccentric giant planets, and a positive correlation between debris disks and terrestrial planets.

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

Recent gravitational microlensing observations predict a vast population of free-floating giant planets that outnumbers main sequence stars almost twofold. A frequently-invoked mechanism for generating this population is a dynamical instability that incites planet-planet scattering and the ejection of one or more planets in isolated main sequence planetary systems. Here, we demonstrate that this process alone probably cannot represent the sole source of these galactic wanderers. By using straightforward quantitative arguments and N-body simulations, we argue that the observed number of exoplanets exceeds the plausible number of ejected planets per system from scattering. Thus, other potential sources of free-floaters, such as planetary stripping in stellar clusters and post-main-sequence ejection, must be considered.

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.

The tidal evolution of planets orbiting brown dwarfs (BDs) presents an interesting case study because BDs' terrestrial planet forming region is located extremely close-in. In fact, the habitable zones of BDs range from roughly 0.001 to 0.03 AU and for the lowest-mass BDs are located interior to the Roche limit. In contrast with stars, BDs spin up as they age. Thus, the corotation distance moves inward. This has important implications for the tidal evolution of planets around BDs. We used a standard equilibrium tidal model to compute the orbital evolution of a large ensemble of planet-BD systems. We tested the effect of numerous parameters such as the initial semi-major axis and eccentricity, the rotation period of the BD, the masses of both the BD and planet, and the tidal dissipation factors. We find that all planets that form at or beyond the corotation distance and with initial eccentricities smaller than \sim 0.1 are repelled from the BD. Some planets initially interior to corotation can survive if their inward tidal evolution is slower than the BD's spin evolution, but most initially close-in planets fall onto the BD. We find that the most important parameter for the tidal evolution is the initial orbital distance with respect to the corotation distance. Some planets can survive in the habitable zone for Gyr timescales, although in many cases the habitable zone moves inward past the planet's orbit in just tens to hundreds of Myr. Surviving planets can have orbital periods of less than 10 days (as small as 10 hrs), so they could be observable by transit.

(abridged) A complete census of planetary systems around a volume-limited sample of solar-type stars (FGK dwarfs) in the Solar neighborhood with uniform sensitivity down to Earth-mass planets within their Habitable Zones out to several AUs would be a major milestone in extrasolar planets astrophysics. This fundamental goal can be achieved with a mission concept such as NEAT - the Nearby Earth Astrometric Telescope. NEAT is designed to carry out space-borne extremely-high-precision astrometric measurements sufficient to detect dynamical effects due to orbiting planets of mass even lower than Earth's around the nearest stars. Such a survey mission would provide the actual planetary masses and the full orbital geometry for all the components of the detected planetary systems down to the Earth-mass limit. The NEAT performance limits can be achieved by carrying out differential astrometry between the targets and a set of suitable reference stars in the field. The NEAT instrument design consists of an off-axis parabola single-mirror telescope, a detector with a large field of view made of small movable CCDs located around a fixed central CCD, and an interferometric calibration system originating from metrology fibers located at the primary mirror. The proposed mission architecture relies on the use of two satellites operating at L2 for 5 years, flying in formation and offering a capability of more than 20,000 reconfigurations (alternative option uses deployable boom). The NEAT primary science program will encompass an astrometric survey of our 200 closest F-, G- and K-type stellar neighbors, with an average of 50 visits. The remaining time might be allocated to improve the characterization of the architecture of selected planetary systems around nearby targets of specific interest (low-mass stars, young stars, etc.) discovered by Gaia, ground-based high-precision radial-velocity surveys.

It has recently been shown that the terrestrial planets and asteroid belt can be reproduced if the giant planets underwent an inward-then-outward migration (the "Grand Tack"; Walsh et al 2011). Inward migration occurs when Jupiter opens a gap and type II migrates inward. The planets "tack" and migrate outward when Saturn reaches the gap-opening mass and is caught in the 3:2 resonance with Jupiter. The aim is to test the viability of the Grand Tack model and to study the dynamical evolution of Jupiter and Saturn during their growth from 10 Earth masses cores. We have performed numerical simulations using a grid-based hydrodynamical code. Most of our simulations assume an isothermal equation of state for the disk but a subset use a fully-radiative version of the code. For an isothermal disk the two phase migration of Jupiter and Saturn is very robust and independent of the mass-growth history of these planets provided the disk is cool enough. For a radiative disk the we find some outcomes with two phase migrations and others with more complicated behavior. We construct a simple, 1-D model of an evolving viscous disk to calculate the evolution of the disk's radiative properties: the disk transitions from radiative to isothermal from its outermost regions inward in time. We show that a two-phase migration is a natural outcome at late times even under the limiting assumption that isothermal conditions are required. Thus, our simulations provide strong support for the Grand Tack scenario.

The bright star 55 Cancri is known to host five planets, including a transiting super-Earth. The study presented here yields directly determined values for 55 Cnc's stellar astrophysical parameters based on improved interferometry: $R=0.943 \pm 0.010 R_{\odot}$, $T_{\rm EFF} = 5196 \pm 24$ K. We use isochrone fitting to determine 55 Cnc's age to be 10.2 $\pm$ 2.5 Gyr, implying a stellar mass of $0.905 \pm 0.015 M_{\odot}$. Our analysis of the location and extent of the system's habitable zone (0.67--1.32 AU) shows that planet f, with period $\sim$ 260 days and $M \sin i = 0.155 M_{Jupiter}$, spends the majority of the duration of its elliptical orbit in the circumstellar habitable zone. Though planet f is too massive to harbor liquid water on any planetary surface, we elaborate on the potential of alternative low-mass objects in planet f's vicinity: a large moon, and a low-mass planet on a dynamically stable orbit within the habitable zone. Finally, our direct value for 55 Cancri's stellar radius allows for a model-independent calculation of the physical diameter of the transiting super-Earth 55 Cnc e ($\sim 2.05 \pm 0.15 R_{\earth}$), which, depending on the planetary mass assumed, implies a bulk density of 0.76 $\rho_{\earth}$ or 1.07 $\rho_{\earth}$.

The eccentric orbits of the known extrasolar giant planets provide evidence that most planet-forming environments undergo violent dynamical instabilities. Here, we numerically simulate the impact of giant planet instabilities on planetary systems as a whole. We find that populations of inner rocky and outer icy bodies are both shaped by the giant planet dynamics and are naturally correlated. Strong instabilities -- those with very eccentric surviving giant planets -- completely clear out their inner and outer regions. In contrast, systems with stable or low-mass giant planets form terrestrial planets in their inner regions and outer icy bodies produce dust that is observable as debris disks at mid-infrared wavelengths. Fifteen to twenty percent of old stars are observed to have bright debris disks (at wavelengths of ~70 microns) and we predict that these signpost dynamically calm environments that should contain terrestrial planets.

Circumstantial evidence suggests that most known extra-solar planetary systems are survivors of violent dynamical instabilities. Here we explore how giant planet instabilities affect the formation and survival of terrestrial planets. We simulate planetary system evolution around Sun-like stars from initial conditions that comprise: an inner disk of planetesimals and planetary embryos, three giant planets at Jupiter-Saturn distances, and a massive outer planetesimal disk. We then calculate dust production rates and debris disk SEDs assuming that each planetesimal particle represents an ensemble of smaller bodies in collisional equilibrium. We predict a strong correlation between the presence of terrestrial planets and debris disks, mediated by the giant planets. Strong giant planet instabilities destroy all rocky material - including fully-formed terrestrial planets if the instabilities occur late - along with the icy planetesimals. Stable or weakly unstable systems allow terrestrial planets to accrete and significant dust to be produced in their outer regions. Stars older than ~100 Myr with bright cold dust emission (at ~70 microns) signpost the dynamically calm environments conducive to efficient terrestrial accretion. We predict that while the typical eccentricities of terrestrial planets are small, there should exist a novel class of terrestrial planet system whose single planet undergoes large amplitude oscillations in eccentricity and inclination. By scaling to the observed semimajor axis distribution of giant exoplanets, we estimate that terrestrial exoplanets in the same systems should be a few times more abundant at 0.5 AU than giant or terrestrial exoplanets at 1 AU. Finally, we discuss the Solar System, which appears to be unusual in combining a rich terrestrial planet system with a low dust content.

The search for rocky exoplanets plays an important role in our quest for extra-terrestrial life. Here, we discuss the extreme physical properties possible for the first characterized rocky super-Earth, CoRoT-7b (R_pl = 1.58 \pm 0.10 R_Earth, Mpl = 6.9 \pm 1.2 M_Earth). It is extremely close to its star (a = 0.0171 AU = 4.48 R_st), with its spin and orbital rotation likely synchronized. The comparison of its location in the (Mpl, Rpl) plane with the predictions of planetary models for different compositions points to an Earth-like composition, even if the error bars of the measured quantities and the partial degeneracy of the models prevent a definitive conclusion. The proximity to its star provides an additional constraint on the model. It implies a high extreme-UV flux and particle wind, and the corresponding efficient erosion of the planetary atmosphere especially for volatile species including water. Consequently, we make the working hypothesis that the planet is rocky with no volatiles in its atmosphere, and derive the physical properties that result. As a consequence, the atmosphere is made of rocky vapours with a very low pressure (P \leq 1.5 Pa), no cloud can be sustained, and no thermalisation of the planetary is expected. The dayside is very hot (2474 \leq 71 K at the sub-stellar point) while the nightside is very cold (50 to 75 K). The sub-stellar point is as hot as the tungsten filament of an incandescent bulb, resulting in the melting and distillation of silicate rocks and the formation of a lava ocean. These possible features of CoRoT-7b could be common to many small and hot planets, including the recently discovered Kepler-10b. They define a new class of objects that we propose to name "Lava-ocean planets".

In the outer regions of the habitable zone, the risk of transitioning into a globally frozen "snowball" state poses a threat to the habitability of planets with the capacity to host water-based life. We use a one-dimensional energy balance climate model (EBM) to examine how obliquity, spin rate, orbital eccentricity, and ocean coverage might influence the onset of such a snowball state. For an exoplanet, these parameters may be strikingly different from the values observed for Earth. Since, for constant semimajor axis, the annual mean stellar irradiation scales with (1-e^2)^(-1/2), one might expect the greatest habitable semimajor axis (for fixed atmospheric composition) to scale as (1-e^2)^(-1/4). We find that this standard ansatz provides a reasonable lower bound on the outer boundary of the habitable zone, but the influence of obliquity and ocean fraction can be profound in the context of planets on eccentric orbits. For planets with eccentricity 0.5, our EBM suggests that the greatest habitable semimajor axis can vary by more than 0.8 AU (78%!) depending on obliquity, with higher obliquity worlds generally more stable against snowball transitions. One might also expect that the long winter at an eccentric planet's apoastron would render it more susceptible to global freezing. Our models suggest that this is not a significant risk for Earth-like planets around Sun-like stars since such planets are buffered by the thermal inertia provided by oceans covering at least 10% of their surface. Since planets on eccentric orbits spend much of their year particularly far from the star, such worlds might turn out to be especially good targets for direct observations with missions such as TPF-Darwin. Nevertheless, the extreme temperature variations achieved on highly eccentric exo-Earths raise questions about the adaptability of life to marginally or transiently habitable conditions.

Tides may be crucial to the habitability of exoplanets. If such planets form around low-mass stars, then those in the circumstellar habitable zone will be close enough to their host stars to experience strong tidal forces. Tides may result in orbital decay and circularization, evolution toward zero obliquity, a fixed rotation rate (not necessarily synchronous), and substantial internal heating. Due to tidal effects, the range of habitable orbital locations may be quite different from that defined by the traditional concept of a habitable zone (HZ) based on stellar insolation, atmospheric effects, and liquid water on a planet's surface. Tidal heating may make locations within the traditional HZ too hot, while planets outside the traditional zone could be rendered quite habitable due to tides. Here we consider these effects on the exoplanet GJ 581 d.

Although the Earth's orbit is never far from circular, terrestrial planets around other stars might experience substantial changes in eccentricity that could lead to climate changes, including possible "phase transitions" such as the snowball transition (or its opposite). There is evidence that Earth has gone through at least one globally frozen, "snowball" state in the last billion years, which it is thought to have exited after several million years because global ice-cover shut off the carbonate-silicate cycle, thereby allowing greenhouse gases to build up to sufficient concentration to melt the ice. Due to the positive feedback caused by the high albedo of snow and ice, susceptibility to falling into snowball states might be a generic feature of water-rich planets with the capacity to host life. This paper has two main thrusts. First, we revisit one-dimensional energy balance climate models as tools for probing possible climates of exoplanets, investigate the dimensional scaling of such models, and introduce a simple algorithm to treat the melting of the ice layer on a globally-frozen planet. We show that if a terrestrial planet undergoes Milankovitch-like oscillations of eccentricity that are of great enough magnitude, it could melt out of a snowball state. Second, we examine the kinds of variations of eccentricity that a terrestrial planet might experience due to the gravitational influence of a giant companion. We show that a giant planet on a sufficiently eccentric orbit can excite extreme eccentricity oscillations in the orbit of a habitable terrestrial planet. More generally, these two results demonstrate that the longterm habitability (and astronomical observables) of a terrestrial planet can depend on the detailed architecture of the planetary system in which it resides.

CoRoT-7 b is the first confirmed rocky exoplanet, but, with an orbital semi-major axis of 0.0172 AU, its origins may be unlike any rocky planet in our solar system. In this study, we consider the roles of tidal evolution and evaporative mass loss in CoRoT-7 b's history, which together have modified the planet's mass and orbit. If CoRoT-7 b has always been a rocky body, evaporation may have driven off almost half its original mass, but the mass loss may depend sensitively on the extent of tidal decay of its orbit. As tides caused CoRoT-7 b's orbit to decay, they brought the planet closer to its host star, thereby enhancing the mass loss rate. Such a large mass loss also suggests the possibility that CoRoT-7 b began as a gas giant planet and had its original atmosphere completely evaporated. In this case, we find that CoRoT-7 b's original mass probably didn't exceed 200 Earth masses (about 2/3 of a Jupiter mass). Tides raised on the host star by the planet may have significantly reduced the orbital semi-major axis, perhaps causing the planet to migrate through mean-motion resonances with the other planet in the system, CoRoT-7 c. The coupling between tidal evolution and mass loss may be important not only for CoRoT-7 b but also for other close-in exoplanets, and future studies of mass loss and orbital evolution may provide insight into the origin and fate of close-in planets, both rocky and gaseous.

We develop an idealized dynamical model to predict the typical properties of outer extrasolar planetary systems, at radii beyond 5 AU. Our hypothesis is that dynamical evolution in outer planetary systems is controlled by a combination of planet-planet scattering and planetary interactions with an exterior disk of small bodies ("planetesimals"). Using 5,000 long duration N-body simulations, we follow the evolution of three planets surrounded by a 50 Earth mass primordial planetesimal disk. For large planet masses (above that of Saturn) the influence of the disk is modest, and we recover the observed eccentricity distribution of extrasolar planets (observed primarily at smaller radii). We explain the observed mass dependence of the eccentricity by invoking strong correlations between planet masses in the same system. For lower mass planets we observe diverse dynamical behavior: strong scattering events, sudden jumps in eccentricity due to resonance crossings, and re-circularization of scattered low-mass planets in the disk. We present distributions of the final eccentricity and inclination, and discuss how they vary with planet mass and initial system architecture. We predict a transition to lower eccentricities for low mass planets at radii where disks influence the dynamics. Radial velocity measurements capable of detecting planets with K~5 m/s and periods in excess of 10 years will constrain this regime. We also study the population of resonant and non-resonant multiple planet systems. We show that, among systems with Jupiter-mass planets that avoid close encounters, the planetesimal disk acts as a damping mechanism that frequently populates mean motion resonances. Resonant chains ought to be common among massive planets in outer planetary systems.

We review how tides may impact the habitability of terrestrial-like planets. If such planets form around low-mass stars, then planets in the circumstellar habitable zone will be close enough to their host stars to experience strong tidal forces. We discuss 1) decay of semi-major axis, 2) circularization of eccentric orbits, 3) evolution toward zero obliquity, 4) fixed rotation rates (not necessarily synchronous), and 5) internal heating. We briefly describe these effects using the example of a 0.25 solar mass star with a 10 Earth-mass companion. We suggest that the concept of a habitable zone should be modified to include the effects of tides.

CoRoT-7 b, a planet about 70% larger than the Earth orbiting a Sun-like star, is the first-discovered rocky exoplanet, and hence has been dubbed a "super-Earth". Some initial studies suggested that since the planet is so close to its host star, it receives enough insolation to partially melt its surface. However, these past studies failed to take into consideration the role that tides may play in this system. Even if the planet's eccentricity has always been zero, we show that tidal decay of semi-major axis could have been large enough that the planet formed on a wider orbit which received less insolation. Moreover, CoRoT-7 b could be tidally heated at a rate that dominates its geophysics and drives extreme volcanism. In this case, CoRoT-7 b is a "super-Io" that, like Jupiter's volcanic moon, is dominated by volcanism and rapid resurfacing. Such heating could occur with an eccentricity of just 10^-5. This small value could be driven by CoRoT-7 c if its own eccentricity is larger than ~10^-4. CoRoT-7 b may be the first of a class of planetary super-Ios likely to be revealed by the CoRoT and Kepler spacecraft.

We review the problem of the formation of terrestrial planets, with particular emphasis on the interaction of dynamical and geochemical models. The lifetime of gas around stars in the process of formation is limited to a few million years based on astronomical observations, while isotopic dating of meteorites and the Earth-Moon system suggest that perhaps 50-100 million years were required for the assembly of the Earth. Therefore, much of the growth of the terrestrial planets in our own system is presumed to have taken place under largely gas-free conditions, and the physics of terrestrial planet formation is dominated by gravitational interactions and collisions. The earliest phase of terrestrial-planet formation involve the growth of km-sized or larger planetesimals from dust grains, followed by the accumulations of these planetesimals into ~100 lunar- to Mars-mass bodies that are initially gravitationally isolated from one-another in a swarm of smaller planetesimals, but eventually grow to the point of significantly perturbing one-another. The mutual perturbations between the embryos, combined with gravitational stirring by Jupiter, lead to orbital crossings and collisions that drive the growth to Earth-sized planets on a timescale of 10-100 million years. Numerical treatment of this process has focussed on the use of symplectic integrators which can rapidly integrate the thousands of gravitationally-interacting bodies necessary to accurately model planetary growth. While the general nature of the terrestrial planets--their sizes and orbital parameters--seem to be broadly reproduced by the models, there are still some outstanding dynamical issues. One of these is the presence of an embryo-sized body, Mars, in our system in place of the more massive objects that simulations tend to yield. [Abridged]

The habitable zones of main sequence stars have traditionally been defined as the range of orbits that intercept the appropriate amount of stellar flux to permit surface water on a planet. Terrestrial exoplanets discovered to orbit M stars in these zones, which are close-in due to decreased stellar luminosity, may also undergo significant tidal heating. Tidal heating may span a wide range for terrestrial exoplanets and may significantly affect conditions near the surface. For example, if heating rates on an exoplanet are near or greater than that on Io (where tides drive volcanism that resurface the planet at least every 1 Myr) and produce similar surface conditions, then the development of life seems unlikely. On the other hand, if the tidal heating rate is less than the minimum to initiate plate tectonics, then CO_2 may not be recycled through subduction, leading to a runaway greenhouse that sterilizes the planet. These two cases represent potential boundaries to habitability and are presented along with the range of the traditional habitable zone for main sequence, low-mass stars. We propose a revised habitable zone that incorporates both stellar insolation and tidal heating. We apply these criteria to GJ 581 d and find that it is in the traditional habitable zone, but its tidal heating alone may be insufficient for plate tectonics.

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include 1) the orbits, in particular the small eccentricities, and 2) the masses of the terrestrial planets -- Mars' relatively small mass in particular has not been adequately reproduced in previous simulations; 3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; 4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and 5) Earth's relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N=1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the "Nice model." We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars' mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. (abridged)

We study the final architecture of planetary systems that evolve under the combined effects of planet-planet and planetesimal scattering. Using N-body simulations we investigate the dynamics of marginally unstable systems of gas and ice giants both in isolation and when the planets form interior to a planetesimal belt. The unstable isolated systems evolve under planet-planet scattering to yield an eccentricity distribution that matches that observed for extrasolar planets. When planetesimals are included the outcome depends upon the total mass of the planets. For system masses exceeding about one Jupiter mass the final eccentricity distribution remains broad, whereas for lower mass planetary systems a combination of divergent orbital evolution and recircularization of scattered planets results in a preponderance of nearly circular final orbits. We also study the fate of marginally stable multiple planet systems in the presence of planetesimal disks, and find that for high planet masses the majority of such systems evolve into resonance. A significant fraction lead to resonant chains that are planetary analogs of Jupiter's Galilean satellites. We predict that a transition from eccentric to near-circular orbits will be observed once extrasolar planet surveys detect sub-Jovian mass planets at orbital radii of 5-10 AU.

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

The known extrasolar multiple-planet systems share a surprising dynamical attribute: they cluster just beyond the Hill stability boundary. Here we show that the planet-planet scattering model, which naturally explains the observed exoplanet eccentricity distribution, can reproduce the observed distribution of dynamical configurations. We calculated how each of our scattered systems would appear over an appropriate range of viewing geometries; as Hill stability is weakly dependent on the masses, the mass-inclination degeneracy does not significantly affect our results. We consider a wide range of initial planetary mass distributions and find that some are poor fits to the observed systems. In fact, many of our scattering experiments overproduce systems very close to the stability boundary. The distribution of dynamical configurations of two-planet systems actually may provide better discrimination between scattering models than the distribution of eccentricity. Our results imply that, at least in their inner regions which are weakly affected by gas or planetesimal disks, planetary systems should be "packed", with no large gaps between planets.

We study the dynamical stability of an additional, potentially habitable planet in the HD 47186 planetary system. Two planets are currently known in this system: a "hot Neptune" with a period of 4.08 days and a Saturn-mass planet with a period of 3.7 years. Here we consider the possibility that one or more undetected planets exist between the two known planets and possibly within the habitable zone in this system. Given the relatively low masses of the known planets, additional planets could have masses $\lsim 10 \mearth$, and hence be terrestrial-like and further improving potential habitability. We perform $N$-body simulations to identify the stable zone between planets $b$ and $c$ and find that much of the inner habitable zone can harbor a $10 \mearth$ planet. With the current radial-velocity threshold of $\sim 1$ m/s, an additional planet should be detectable if it lies at the inner edge of the habitable zone at 0.8 AU. We also show that the stable zone could contain two additional planets of $10 \mearth$ each if their eccentricities are lower than $\sim 0.3$.

In agreement with previous work, we show that the presence of the short-lived radionuclide Al-26 in the early Solar System was unlikely (<2% a priori probability) to be the result of direct introduction of supernova ejecta into the gaseous disk during the Class II stage of protosolar evolution. We also show that any Bondi-Hoyle accretion of contaminated residual gas from the natal star cluster made a negligible contribution to the primordial Al-26 inventory of the Solar System. These results are consistent with the absence of the oxygen isotopic signature expected with any late introduction of supernova ejecta into the protoplanetary disk. Instead, the presence of Al-26 in the oldest Solar System solids (calcium-aluminum-rich inclusions or CAIs) and its apparent uniform distribution with the inferred canonical Al-26/Al-27 ratio of (4.5-5)E-5 support the inheritance of Al-26 from the parent giant molecular cloud. We propose that this radionuclide originated in a prior generation of massive stars that formed in the same host molecular cloud as the Sun and contaminated that cloud by Wolf-Rayet winds. We calculated the Galactic distribution of Al-26/Al-27 ratios that arise from such contamination using the established embedded cluster mass and stellar initial mass functions, published nucleosynthetic yields from the winds of massive stars, and by assuming rapid and uniform mixing into the cloud. Although our model predicts that the majority of stellar systems contain no Al-26 from massive stars, and that the a priori probability that the Al-26/Al-27 ratio will reach or exceed the canonical Solar System value is only ~6%, the maximum in the distribution of non-zero values is close to the canonical ratio.

Three planets with minimum masses less than 10 Earth masses orbit the star HD 40307, suggesting these planets may be rocky. However, with only radial velocity data, it is impossible to determine if these planets are rocky or gaseous. Here we exploit various dynamical features of the system in order to assess the physical properties of the planets. Observations allow for circular orbits, but a numerical integration shows that the eccentricities must be at least 0.0001. Also, planets b and c are so close to the star that tidal effects are significant. If planet b has tidal parameters similar to the terrestrial planets in the Solar System and a remnant eccentricity larger than 0.001, then, going back in time, the system would have been unstable within the lifetime of the star (which we estimate to be 6.1 +/- 1.6 Gyr). Moreover, if the eccentricities are that large and the inner planet is rocky, then its tidal heating may be an order of magnitude greater than extremely volcanic Io, on a per unit surface area basis. If planet b is not terrestrial, e.g. Neptune-like, these physical constraints would not apply. This analysis suggests the planets are not terrestrial-like, and are more like our giant planets. In either case, we find that the planets probably formed at larger radii and migrated early-on (via disk interactions) into their current orbits. This study demonstrates how the orbital and dynamical properties of exoplanet systems may be used to constrain the planets' physical properties.

Planet-planet scattering is the leading mechanism to explain the large eccentricities of the observed exoplanet population. However, scattering has not been considered important to the production of pairs of planets in mean motion resonances (MMRs). We present results from a large number of numerical simulations of dynamical instabilities in 3-planet systems. We show that MMRs arise naturally in about five percent of cases. The most common resonances we populate are the 2:1 and 3:1 MMRs, although a wide variety of MMRs can occur, including high-order MMRs (up to eleventh order). MMRs are generated preferentially in systems with uneven mass distributions: the smallest planet is typically ejected after a series of close encounters, leaving the remaining, more massive planets in resonance. The distribution of resonant planets is consistent with the phase-space density of resonant orbits, meaning that planets are randomly thrown into MMRs rather than being slowly pulled into them. It may be possible to distinguish between MMRs created by scattering vs. convergent migration in a gaseous disk by considering planetary mass ratios: resonant pairs of planets beyond ~1 AU with more massive outer planets are likely to have formed by scattering. In addition, scattering may be responsible for pairs of planets in high-order MMRs (3:1 and higher) that are not easily populated by migration. The frequency of MMRs from scattering is comparable to the expected survival rate of MMRs in turbulent disks. Thus, planet-planet scattering is likely to be a major contributor to the population of resonant planets.

Five planets are known to orbit the star 55 Cancri. The recently-discovered planet f at 0.78 AU (Fischer et al. 2008) is located at the inner edge of a previously-identified stable zone that separates the three close-in planets from planet d at 5.9 AU. Here we map the stability of the orbital space between planets f and d using a suite of n-body integrations that include an additional, yet-to-be-discovered planet g with a radial velocity amplitude of 5 m/s (planet mass = 0.5-1.2 Saturn masses). We find a large stable zone extending from 0.9 to 3.8 AU at eccentricities below 0.4. For each system we quantify the probability of detecting planets b-f on their current orbits given perturbations from hypothetical planet g, in order to further constrain the mass and orbit of an additional planet. We find that large perturbations are associated with specific mean motion resonances (MMRs) with planets f and d. We show that two MMRs, 3f:1g (the 1:3 MMR between planets g and f) and 4g:1d cannot contain a planet g. The 2f:1g MMR is unlikely to contain a planet more massive than about 20 Earth masses. The 3g:1d and 5g:2d MMRs could contain a resonant planet but the resonant location is strongly confined. The 3f:2g, 2g:1d and 3g:2d MMRs exert a stabilizing influence and could contain a resonant planet. Furthermore, we show that the stable zone may in fact contain 2-3 additional planets, if they are ~50 Earth masses each. Finally, we show that any planets exterior to planet d must reside beyond 10 AU.

Gravitational scattering between massive planets has been invoked to explain the eccentricity distribution of extrasolar planets. For scattering to occur, the planets must either form in -- or migrate into -- an unstable configuration. In either case, it is likely that a residual gas disk, with a mass comparable to that of the planets, will be present when scattering occurs. Using explicit hydrodynamic simulations, we study the impact of gas disks on the outcome of two-planet scattering. We assume a specific model in which the planets are driven toward instability by gravitational torques from an outer low mass disk. We find that the accretion of mass and angular momentum that occurs when a scattered planet impacts the disk can strongly influence the subsequent dynamics by reducing the number of close encounters. The eccentricity of the innermost surviving planet at the epoch when the system becomes Hill stable is not substantially altered from the gas-free case, but the outer planet is circularized by its interaction with the disk. The signature of scattering initiated by gas disk migration is thus a high fraction of low eccentricity planets at larger radii accompanying known eccentric planets. Subsequent secular evolution of the two planets in the presence of damping can further damp both eccentricities, and tends to push systems away from apsidal alignment and toward anti-alignment. We note that the late burst of accretion when the outer planet impacts the disk is in principle observable, probably via detection of a strong near-IR excess in systems with otherwise weak disk and stellar accretion signatures.

Tides raised on a planet by its host star's gravity can reduce a planet's orbital semi-major axis and eccentricity. This effect is only relevant for planets orbiting very close to their host stars. The habitable zones of low-mass stars are also close-in and tides can alter the orbits of planets in these locations. We calculate the tidal evolution of hypothetical terrestrial planets around low-mass stars and show that tides can evolve planets past the inner edge of the habitable zone, sometimes in less than 1 billion years. This migration requires large eccentricities (>0.5) and low-mass stars (<0.35 M_Sun). Such migration may have important implications for the evolution of the atmosphere, internal heating and the Gaia hypothesis. Similarly, a planet detected interior to the habitable zone could have been habitable in the past. We consider the past habitability of the recently-discovered, ~5 M_Earth planet, Gliese 581 c. We find that it could have been habitable for reasonable choices of orbital and physical properties as recently as 2 Gyr ago. However, when we include constraints derived from the additional companions, we see that most parameter choices that predict past habitability require the two inner planets of the system to have crossed their mutual 3:1 mean motion resonance. As this crossing would likely have resulted in resonance capture, which is not observed, we conclude that Gl 581 c was probably never habitable.

Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that would place these systems near instability. The hypothetical planet in HD 74156 was expected to have a mass about equal to that of Saturn, a semi-major axis between 0.9 and 1.4 AU, and an eccentricity less than 0.2. HD 74156 d, a planet with a mass of 1.3 Saturn masses at 1.04 AU with an eccentricity of 0.25, was recently announced. We have reanalyzed all published data on this system in order to place tighter constraints on the properties of the new planet. We find two possible orbits for this planet, one close to that already identified and another (with a slightly better fit to the data) at ~0.89 AU. We also review the current status of other planet predictions, discuss the observed single planet systems, and suggest other systems which may contain planets in between those that are already known. The confirmation of the existence of HD 74156 d suggests that planet formation is an efficient process, and planetary systems should typically contain many planets.

Terrestrial planets form in a series of dynamical steps from the solid component of circumstellar disks. First, km-sized planetesimals form likely via a combination of sticky collisions, turbulent concentration of solids, and gravitational collapse from micron-sized dust grains in the thin disk midplane. Second, planetesimals coalesce to form Moon- to Mars-sized protoplanets, also called "planetary embryos". Finally, full-sized terrestrial planets accrete from protoplanets and planetesimals. This final stage of accretion lasts about 10-100 Myr and is strongly affected by gravitational perturbations from any gas giant planets, which are constrained to form more quickly, during the 1-10 Myr lifetime of the gaseous component of the disk. It is during this final stage that the bulk compositions and volatile (e.g., water) contents of terrestrial planets are set, depending on their feeding zones and the amount of radial mixing that occurs. The main factors that influence terrestrial planet formation are the mass and surface density profile of the disk, and the perturbations from giant planets and binary companions if they exist. Simple accretion models predicts that low-mass stars should form small, dry planets in their habitable zones. The migration of a giant planet through a disk of rocky bodies does not completely impede terrestrial planet growth. Rather, "hot Jupiter" systems are likely to also contain exterior, very water-rich Earth-like planets, and also "hot Earths", very close-in rocky planets. Roughly one third of the known systems of extra-solar (giant) planets could allow a terrestrial planet to form in the habitable zone.

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets (< 5-10 Earth masses). Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form "hot Earths" and their observable predictions. Models include: 1) in situ accretion; 2) formation at larger orbital distance followed by inward "type 1" migration; 3) formation from material being "shepherded" inward by a migrating gas giant planet; 4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disk; 5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances; and 6) photo-evaporative mass loss of a close-in giant planet. Models 1-4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets (> 5 Earth masses) with very close-in perihelion distances (< 0.025 AU), and even then the net inward movement in orbital distance is at most only 0.1-0.15 AU. For planets of less than about 70 Earth masses, photo-evaporation can remove the planet's envelope and leave behind the solid core on a Gyr timescale, but only for planets inside 0.025-0.05 AU. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit-timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).

The search for habitable planets like Earth around other stars fulfils an ancient imperative to understand our origins and place in the cosmos. The past decade has seen the discovery of hundreds of planets, but nearly all are gas giants like Jupiter and Saturn. Recent advances in instrumentation and new missions are extending searches to planets the size of the Earth, but closer to their host stars. There are several possible ways such planets could form, and future observations will soon test those theories. Many of these planets we discover may be quite unlike Earth in their surface temperature and composition, but their study will nonetheless inform us about the process of planet formation and the frequency of Earth-like planets around other stars.

Smaller terrestrial planets (< 0.3 Earth masses) are less likely to retain the substantial atmospheres and ongoing tectonic activity probably required to support life. A key element in determining if sufficiently massive "sustainably habitable" planets can form is the availability of solid planet-forming material. We use dynamical simulations of terrestrial planet formation from planetary embryos and simple scaling arguments to explore the implications of correlations between terrestrial planet mass, disk mass, and the mass of the parent star. We assume that the protoplanetary disk mass scales with stellar mass as Mdisk ~ f Mstar^h, where f measures the relative disk mass, and 1/2 < h < 2, so that disk mass decreases with decreasing stellar mass. We consider systems without Jovian planets, based on current models and observations for M stars. We assume the mass of a planet formed in some annulus of a disk with given parameters is proportional to the disk mass in that annulus, and show with a suite of simulations of late-stage accretion that the adopted prescription is surprisingly accurate. Our results suggest that the fraction of systems with sufficient disk mass to form > 0.3 Earth mass habitable planets decreases for low-mass stars for every realistic combination of parameters. This "habitable fraction" is small for stellar masses below a mass in the interval 0.5 to 0.8 Solar masses, depending on disk parameters, an interval that excludes most M stars. Radial mixing and therefore water delivery are inefficient in lower-mass disks commonly found around low-mass stars, such that terrestrial planets in the habitable zones of most low-mass stars are likely to be small and dry.

One of the most surprising discoveries of extrasolar planets is the detection of planets in moderately close binary star systems. The Jovian-type planets in the two binaries of Gamma Cephei and GJ 86 have brought to the forefront questions on the formation of giant planets and the possibility of the existence of smaller bodies in such dynamically complex environments. The diverse dynamical characteristics of these objects have made scientists wonder to what extent the current theories of planet formation can be applied to binaries and multiple star systems. At present, the sensitivity of the detection techniques does not allow routine discovery of Earth-sized bodies in binary systems. However, with the advancement of new techniques, and with the recent launch of CoRoT and the launch of Kepler in late 2008, the detection of more planets (possibly terrestrial-class objects) in such systems is on the horizon. Theoretical studies and numerical modeling of terrestrial and habitable planet formation are, therefore, necessary to gain fundamental insights into the prospects for life in such systems and have great strategic impact on NASA science and missions.

Recent radial velocity observations have indicated that Jovian-type planets can exist in moderately close binary star systems. Numerical simulations of the dynamical stability of terrestrial-class planets in such environments have shown that, in addition to their giant planets, these systems can also harbor Earth-like objects. In this paper, we study the late stage of terrestrial planet formation in such binary-planetary systems, and present the results of the simulations of the formation of Earth-like bodies in their habitable zones. We consider a circumprimary disk of Moon- to Mars-sized objects and numerically integrate the orbits of these bodies at the presence of the Jovian-type planet of the system and for different values of the mass, semimajor axis, and orbital eccentricity of the secondary star. Results indicate that, Earth-like objects, with substantial amounts of water, can form in the habitable zone of the primary star. Simulations also indicate that, by transferring angular momentum from the secondary star to protoplanetary objects, the giant planet of the system plays a key role in the radial mixing of these bodies and the water contents of the final terrestrial planets. We will discuss the results of our simulation and show that the formation of habitable planets in binary-planetary systems is more probable in binaries with moderate to large perihelia.

The water content and habitability of terrestrial planets are determined during their final assembly, from perhaps a hundred 1000-km "planetary embryos" and a swarm of billions of 1-10 km "planetesimals." During this process, we assume that water-rich material is accreted by terrestrial planets via impacts of water-rich bodies that originate in the outer asteroid region. We present analysis of water delivery and planetary habitability in five high-resolution simulations containing about ten times more particles than in previous simulations (Raymond et al 2006a, Icarus, 183, 265-282). These simulations formed 15 terrestrial planets from 0.4 to 2.6 Earth masses, including five planets in the habitable zone. Every planet from each simulation accreted at least the Earth's current water budget; most accreted several times that amount (assuming no impact depletion). Each planet accreted at least five water-rich embryos and planetesimals from past 2.5 AU; most accreted 10-20 water-rich bodies. We present a new model for water delivery to terrestrial planets in dynamically calm systems, with low-eccentricity or low-mass giant planets -- such systems may be very common in the Galaxy. We suggest that water is accreted in comparable amounts from a few planetary embryos in a "hit or miss" way and from millions of planetesimals in a statistically robust process. Variations in water content are likely to be caused by fluctuations in the number of water-rich embryos accreted, as well as from systematic effects such as planetary mass and location, and giant planet properties.

Close-in giant planets are thought to have formed in the cold outer regions of planetary systems and migrated inward, passing through the orbital parameter space occupied by the terrestrial planets in our own Solar System. We present dynamical simulations of the effects of a migrating giant planet on a disk of protoplanetary material and the subsequent evolution of the planetary system. We numerically investigate the dynamics of post-migration planetary systems over 200 million years using models with a single migrating giant planet, one migrating and one non-migrating giant planet, and excluding the effects of a gas disk. Material that is shepherded in front of the migrating giant planet by moving mean motion resonances accretes into "hot Earths", but survival of these bodies is strongly dependent on dynamical damping. Furthermore, a significant amount of material scattered outward by the giant planet survives in highly excited orbits; the orbits of these scattered bodies are then damped by gas drag and dynamical friction over the remaining accretion time. In all simulations Earth-mass planets accrete on approximately 100 Myr timescales, often with orbits in the Habitable Zone. These planets range in mass and water content, with both quantities increasing with the presence of a gas disk and decreasing with the presence of an outer giant planet. We use scaling arguments and previous results to derive a simple recipe that constrains which giant planet systems are able to form and harbor Earth-like planets in the Habitable Zone, demonstrating that roughly one third of the known planetary systems are potentially habitable.

Astrobiology, the study of life as a planetary phenomenon, aims to understand the fundamental nature of life on earth and the possibility of life elsewhere. To achieve this goal, astrobiologists have initiated unprecedented communication between the disciplines of astronomy, biology, chemistry, and geology. The Astrobiology Primer has been created as a reference tool for those who are interested in the interdisciplinary field of astrobiology. The field incorporates many diverse research endeavors, but it is our hope that this slim volume will present the reader with all he or she needs to know to become involved and to understand, at least at a fundamental level, the state of the art. Because of the great diversity of material, each section was written by a different author with a different expertise. The Primer was constructed collaboratively. Ninety researchers from around the world contributed information with regard to what they expected from other astrobiologists and what they would like to know themselves but still had difficulty understanding (see Contributors). Those submissions were read and considered by the Editors who produced a list of seven general categories of knowledge, represented by the seven chapters in the Primer: 1) Stellar Formation and Evolution, 2) Planetary Formation and Evolution, 3) Astrobiogeochemistry and the Origin of Life, 4) Evolution of Life through Time, 5) Planet Detection & Characterization, 6) Diversity of Life, and 7) Science in Space. No one volume, of course, can contain the vast amount of information brought to play in astrobiology, but we believe that the Primer will provide a forum and a language around which the community will have the opportunity to develop a consensus about central issues.

Close-in giant planets (e.g. ``Hot Jupiters'') are thought to form far from their host stars and migrate inward, through the terrestrial planet zone, via torques with a massive gaseous disk. Here we simulate terrestrial planet growth during and after giant planet migration. Several-Earth mass planets also form interior to the migrating Jovian planet, analogous to recently-discovered ``Hot Earths''. Very water-rich, Earth-mass planets form from surviving material outside the giant planet's orbit, often in the Habitable Zone and with low orbital eccentricities. More than a third of the known systems of giant planets may harbor Earth-like planets.

Recent results have shown that many of the known extrasolar planetary systems contain regions which are stable for both Earth-mass and Saturn-mass planets. Here we simulate the formation of terrestrial planets in four planetary systems -- 55 Cancri, HD 38529, HD 37124, and HD 74156 -- under the assumption that these systems of giant planets are complete and that their orbits are well-determined. Assuming the giant planets formed and migrated quickly, then terrestrial planets may form from a second generation of planetesimals. In each case, Moon- to Mars-sized planetary embryos are placed in between the giant planets and evolved for 100 Myr. We find that planets form relatively easily in 55 Cnc, with masses up to 0.6 Earth masses and in some cases substantial water contents and orbits in the habitable zone. HD 38529 is likely to support an asteroid belt but no terrestrial planets of significant mass. No terrestrial planets form in HD 37124 and HD 74156, although in some cases 1-2 lone embryos survive for 100 Myr. If migration occurred later, depleting the planetesimal disk, then massive terrestrial planets are unlikely to form in any of these systems.

Gas giant planets are far easier than terrestrial planets to detect around other stars, and are thought to form much more quickly than terrestrial planets. Thus, in systems with giant planets, the late stages of terrestrial planet formation are strongly affected by the giant planets' dynamical presence. Observations of giant planet orbits may therefore constrain the systems that can harbor potentially habitable, Earth-like planets. We present results of 460 N-body simulations of terrestrial accretion from a disk of Moon- to Mars-sized planetary embryos. We systematically vary the orbital semimajor axis of a Jupiter-mass giant planet between 1.6 and 6 AU, and eccentricity between 0 and 0.4. We find that for Sun-like stars, giant planets inside roughly 2.5 AU inhibit the growth of 0.3 Earth-mass planets in the habitable zone. If planets accrete water from volatile-rich embryos past 2-2.5 AU, then water-rich habitable planets can only form in systems with giant planets beyond 3.5 AU. Giant planets with significant orbital eccentricities inhibit both accretion and water delivery. The majority of the current sample of extra-solar giant planets appears unlikely to form habitable planets.

The final stage in the formation of terrestrial planets consists of the accumulation of ~1000-km ``planetary embryos'' and a swarm of billions of 1-10 km ``planetesimals.'' During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ~0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: 1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. 2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. 3) If planetary embryos form relatively slowly, following the models of Kokubo & Ida, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. 4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.

The ``minimum-mass solar nebula'' (MMSN) model estimates the surface density distribution of the protoplanetary disk by assuming the planets to have formed in situ. However, significant radial migration of the giant planets likely occurred in the Solar system, implying a distortion in the values derived by the MMSN method. The true density profiles of protoplanetary disks is therefore uncertain. Here we present results of simulations of late-stage terrestrial accretion, each starting from a disk of planetary embryos. We assume a power-law surface density profile that varies with heliocentric distance r as r^-alpha, and vary alpha between 1/2 and 5/2 (alpha = 3/2 for the MMSN model). We find that for steeper profiles (higher values of alpha), the terrestrial planets (i) are more numerous, (ii) form more quickly, (iii) form closer to the star, (iv) are more massive, (v) have higher iron contents, and (vi) have lower water contents. However, the possibility of forming potentially habitable planets does not appear to vary strongly with alpha.

`Hot jupiters,' giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form (Armitage 2003). We present results of simulations of terrestrial planet formation in the presence of hot jupiters, broadly defined as having orbital radii <= 0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g. 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with hot jupiters.

Recent results have shown that many of the known extrasolar planetary systems contain regions which are stable for massless test particles. We examine the possibility that Saturn-mass planets exist in these systems, just below the detection threshold, and attempt to predict likely orbital parameters for such unseen planets. To do this, we insert a Saturn-mass planet into the stable regions of these systems and integrate its orbit for 100 million years. We conduct 200-600 of these experiments to test parameter space in HD37124, HD38529, 55Cnc, and HD74156. In HD37124 the global maximum of the survival rate of Saturns in parameter space is at semimajor axis a = 1.03 AU, eccentricity e=0.1. In HD38529, only 5% of Saturns are unstable, and the region in which a Saturn could survive is very broad, centered on 0.5<a<0.6, e<0.2. In 55Cnc we find three maxima at (a,e) = (1.0 AU, 0.02), (2.0 AU, 0.08), and (3.0 AU, 0.17). In HD74156 we find a broad maximum with $a$ = 0.9-1.2 AU, e<=0.15. Several of these maxima are located in the habitable zones of their parent stars and are therefore of astrobiological interest. We suggest the possibility that companions may lie in these locations of parameter space, and encourage further observational investigation of these systems.

We present a spectroscopic analysis of nearly 8000 late-type dwarfs in the Sloan Digital Sky Survey. Using the Halpha emission line as an activity indicator, we investigate the fraction of active stars as a function of spectral type and find a peak near type M8, confirming previous results. In contrast to past findings, we find that not all M7-M8 stars are active. We show that this may be a selection effect of the distance distributions of previous samples, as the active stars appear to be concentrated near the Galactic Plane. We also examine the activity strength (ratio of the luminosity emitted in Halpha to the bolometric luminosity) for each star, and find that the mean activity strength is constant over the range M0-M5 and declines at later types. The decline begins at a slightly earlier spectral type than previously found. We explore the effect that activity has on the broadband photometric colors and find no significant differences between active and inactive stars. We also carry out a search for subdwarfs using spectroscopic metallicity indicators, and find 60 subdwarf candidates. Several of these candidates are near the extreme subdwarf boundary. The spectroscopic subdwarf candidates are redder by \~0.2 magnitudes in g-r compared to disk dwarfs at the same r-i color.

Recent work has suggested that many planetary systems lie near instability. If all systems are near instability, an additional planet must exist in stable regions of well-separated extra-solar planetary systems to push these systems to the edge of stability. We examine the known systems by placing massless test particles in between the planets and integrating for 1-10 million years. We find that some systems, HD168443 and HD74156, eject nearly all test particles within 2 million years. However we find that HD37124, HD38529, and 55Cnc have large contiguous regions in which particles survive for 10 million years. These three systems, therefore, seem the most likely candidates for additional companions. Furthermore HD74156 and HD168443 must be complete and therefore radial velocity surveys should only focus on detecting more distant companions. We also find that several systems show stable regions that only exist at nonzero eccentricities.

We present a detailed observational strategy for finding Near Earth Objects (NEOs) with the Sloan Digital Sky Survey (SDSS) telescope. We investigate strategies in normal, unbinned mode as well as binning the CCDs 2x2 or 3x3, which affects the sky coverage rate and the limiting apparent magnitude. We present results from 1 month, 3 year and 10 year simulations of such surveys. For each cadence and binning mode, we evaluate the possibility of achieving the Spaceguard goal of detecting 90% of 1 km NEOs (absolute magnitude H <= 18 for an albedo of 0.1). We find that an unbinned survey is most effective at detecting H <= 20 NEOs in our sample. However, a 3x3 binned survey reaches the Spaceguard Goal after only seven years of operation. As the proposed large survey telescopes (PanStarss; LSST) are at least 5-10 years from operation, an SDSS NEO survey could make a significant contribution to the detection and photometric characterization of the NEO population.

The first full year of operation following the commissioning year of the Sloan Digital Sky Survey has revealed a wide variety of newly discovered cataclysmic variables. We show the SDSS spectra of forty-two cataclysmic variables observed in 2002, of which thirty-five are new classifications, four are known dwarf novae (CT Hya, RZ Leo, T Leo and BZ UMa), one is a known CV identified from a previous quasar survey (Aqr1) and two are known ROSAT or FIRST discovered CVs (RX J09445+0357, FIRST J102347.6+003841). The SDSS positions, colors and spectra of all forty-two systems are presented. In addition, the results of follow-up studies of several of these objects identify the orbital periods, velocity curves and polarization that provide the system geometry and accretion properties. While most of the SDSS discovered systems are faint (>18th mag) with low accretion rates (as implied from their spectral characteristics), there are also a few bright objects which may have escaped previous surveys due to changes in the mass transfer rate.

We present results from 42 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary "embryos" (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 42 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 "habitable" planets between 0.9 and 1.1 AU. These planets range from dry worlds to "water worlds" with 100+ oceans of water (1 ocean = 1.5x10^24 g), and vary in mass between 0.23 and 3.85 Earth masses. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call "super embryos," 1-2 Earth mass protoplanets which can serve as the accretion seeds for 2+ Earth mass planets with large water contents.

We have identified 109 White Dwarf (WD) - M dwarf pairs in the Sloan Digital Sky Survey (SDSS) with g < 20th magnitude. For each system we determined the temperature of the WD primary and the spectral type of the M dwarf secondary. Using H-alpha emission as a proxy for the chromospheric activity level of the M dwarf, we investigated correlations between the activity level and properties of the system. Compared with field M dwarfs (Hawley et al. 1996), we see a slightly higher active fraction of early-type stars, with activity levels similar to the field. We have conducted followup observations at the ARC 3.5m telescope to obtain radial velocity information and to search for short period binaries, which may be on the verge of interacting. We report on one system with a 4.1 hour period, and several additional systems with significant velocity variations.