The formation of the solar system's terrestrial planets concluded with a period of giant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20-30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-type wind. Here we show that H2O and other high-mean molecular weight outgassed species can be efficiently lost through this thermal wind if present in a hydrogen-dominated atmosphere, substantially altering the final volatile inventory of terrestrial planets. Specifically, we demonstrate that a giant impact of a Mars-sized embryo with a proto-Earth can remove several Earth oceans' worth of H2O, and other heavier volatile species, together with a primordial hydrogen-dominated atmosphere. These results may offer an explanation for the observed depletion in Earth's light noble gas budget and for its depleted xenon inventory, which suggest that Earth underwent significant atmospheric loss by the end of its accretion. Because planetary embryos are massive enough to accrete primordial hydrogen envelopes and because giant impacts are stochastic and occur concurrently with other early atmospheric evolutionary processes, our results suggest a wide diversity in terrestrial planet volatile budgets.

Recent studies have shown that atmospheric mass-loss powered by the cooling luminosity of a planet's core can explain the observed radius valley separating super-Earths and sub-Neptunes, even without photoevaporation. In this work, we investigate the dependence of this core-powered mass-loss mechanism on stellar mass ($M_\ast$), metallicity ($Z_\ast$) and age ($\tau_\ast$). Without making any changes to the underlying planet population, we find that the core-powered mass-loss model yields a shift in the radius valley to larger planet sizes around more massive stars with a slope given by $\text{d log}R_p/\text{d log}M_\ast \simeq 0.35$, in agreement with observations. To first order, this slope is driven by the dependence of core-powered mass-loss on the bolometric luminosity of the host star and is given by $\text{d log}R_p/\text{d log}M_\ast \simeq (3\alpha-2)/36 \simeq 0.33$, where $(L_\ast/L_\odot) = (M_\ast/M_\odot)^\alpha$ is the stellar mass-luminosity relation and $\alpha\simeq 4.6$ for the CKS dataset. We therefore find, in contrast to photoevaporation models, no evidence for a linear correlation between planet and stellar mass, but can't rule it out either. In addition, we show that the location of the radius valley is, to first order, independent of stellar age and metallicity. Since core-powered mass-loss proceeds over Gyr timescales, the abundance of super-Earths relative to sub-Neptunes increases with age but decreases with stellar metallicity. Finally, due the dependence of the envelope's cooling timescale on metallicity, we find that the radii of sub-Neptunes increase with metallicity and decrease with age with slopes given by $\text{d log}R_p/\text{d log}Z_\ast \simeq 0.1$ and $\text{d log}R_p/\text{d log}\tau_\ast \simeq -0.1$, respectively. We conclude with a series of observational tests that can differentiate between core-powered mass-loss and photoevaporation models.

Oxygen fugacity is a measure of rock oxidation that influences planetary structure and evolution. Most rocky bodies in the Solar System formed at oxygen fugacities approximately five orders of magnitude higher than a hydrogen-rich gas of solar composition. It is unclear whether this oxidation of rocks in the Solar System is typical among other planetary systems. We exploit the elemental abundances observed in six white dwarfs polluted by the accretion of rocky bodies to determine the fraction of oxidized iron in those extrasolar rocky bodies and therefore their oxygen fugacities. The results are consistent with the oxygen fugacities of Earth, Mars, and typical asteroids in the Solar System, suggesting that at least some rocky exoplanets are geophysically and geochemically similar to Earth.

Silicon and Mg in differentiated rocky bodies exhibit heavy isotope enrichments that have been attributed to evaporation of partially or entirely molten planetesimals. We evaluate the mechanisms of planetesimal evaporation in the early solar system and the conditions that controlled attendant isotope fractionations. Energy balance at the surface of a body accreted within ~1 Myr of CAI formation and heated from within by 26Al decay results in internal temperatures exceeding the silicate solidus, producing a transient magma ocean with a thin surface boundary layer of order < 1 meter that would be subject to foundering. Bodies that are massive enough to form magma oceans by radioisotope decay (ge 0.1%) can retain hot rock vapor even in the absence of ambient nebular gas. We find that a steady-state rock vapor forms within minutes to hours and results from a balance between rates of magma evaporation and atmospheric escape. Vapor pressure buildup adjacent to the surfaces of the evaporating magmas would have inevitably led to an approach to equilibrium isotope partitioning between the vapor phase and the silicate melt. Numerical simulations of this near-equilibrium evaporation process for a body with a radius of ~ 700 km yield a steady-state far-field vapor pressure corresponding to 95% saturation. Approaches to equilibrium isotope fractionation between vapor and melt should have been the norm during planet formation due to the formation of steady-state rock vapor atmospheres and/or the presence of protostellar gas. We model the Si and Mg isotopic composition of bulk Earth and show that the best fit is for a carbonaceous chondrite-like source material with about 12% loss of Mg and 15% loss of Si resulting from near-equilibrium evaporation into the solar protostellar disk of hydrogen gas on timescales of 10,000 to 100,000 years.

We present a novel method for determining the surface density of protoplanetary disks through consideration of disk 'dust lines' which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the disk is equal to both the drift timescale and growth timescale of the maximum particle size at a given dust line. We provide an initial proof of concept of our model through an application to the disk TW Hya and are able to estimate the disk dust-to-gas ratio, CO abundance, and accretion rate in addition to the total disk surface density. We find that our derived surface density profile and dust-to-gas ratio are consistent with the lower limits found through measurements of HD gas. The CO ice line also depends on surface density through grain adsorption rates and drift and we find that our theoretical CO ice line estimates have clear observational analogues. We further apply our model to a large parameter space of theoretical disks and find three observational diagnostics that may be used to test its validity. First we predict that the dust lines of disks other than TW Hya will be consistent with the normalized CO surface density profile shape for those disks. Second, surface density profiles that we derive from disk ice lines should match those derived from disk dust lines. Finally, we predict that disk dust and ice lines will scale oppositely, as a function of surface density, across a large sample of disks.

We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of their host stellar mass, substantially larger than the minimum mass solar nebula. All are stable to gravitational collapse except for one which approaches the limit of Toomre-Q stability. Our mass estimates are 2-15 times larger than estimates from integrated optically thin dust emission. We derive depleted dust-to-gas ratios with typical values of ~$10^{-3}$ in the outer disk. Using coagulation models we derive dust surface density profiles that are consistent with millimeter dust observations. In these models, the disks formed with an initial dust mass that is a factor of ~10 greater than is presently observed. Of the three disks in our sample with resolved CO line emission, the masses of HD 163296, AS 209, and TW Hya are roughly 3, 115, and 40 times more massive than estimates from CO respectively. This range indicates that CO depletion is not uniform across different disks and that dust is a more robust tracer of total disk mass. Our method of determining surface density using dust lines is robust even if particles form as aggregates and is useful even in the presence of dust substructure caused by pressure traps. The low Toomre-Q values observed in this sample indicate that at least some disks do not accrete efficiently.

Recent observations revealed a bimodal radius distribution of small, short-period exoplanets with a paucity in their occurrence, a radius `valley', around $1.5-2.0$ $R_\oplus$. In this work, we investigate the effect of a planet's own cooling luminosity on its thermal evolution and atmospheric mass-loss (core-powered mass-loss) and determine its observational consequences for the radius distribution of small, close-in exoplanets. Using simple analytical descriptions and numerical simulations, we demonstrate that planetary evolution based on the core-powered mass-loss mechanism alone (i.e., without any photoevaporation) can produce the observed valley in the radius distribution. Our results match the valley's location, shape and slope in planet radius-orbital period parameter space, and the relative magnitudes of the planet occurrence rate above and below the valley. We find that the slope of the valley is, to first order, dictated by the atmospheric mass-loss timescale at the Bondi radius and given by $\text{d log} R_p/ \text{d log} P \simeq 1/(3(1-\beta))$ which evaluates to $ -0.11$ for $ \beta \simeq 4$, where $M_c/M_\oplus = (R_c/R_{\oplus})^{\beta} (\rho_{c*}/\rho_{\oplus})^{\beta/3}$ is the mass-radius relation of the core. This choice for $\beta$ yields good agreement with observations and attests to the significance of internal compression for massive planetary cores. We further find that the location of the valley scales as $\rho_{c*}^{-4/9}$ and that the observed planet population must have predominantly rocky cores with typical water-ice fractions of less than $\sim 20\%$. Furthermore, we show that the relative magnitude of the planet occurrence rate above and below the valley is sensitive to the details of the planet-mass distribution but that the location of the valley is not.

Using moving mesh hydrodynamic simulations, we determine the shock propagation and resulting ground velocities for a planet hit by a high velocity impactor. We use our results to determine the atmospheric mass loss caused by the resulting ground motion due to the impact shock wave. We find that there are two distinct shock propagation regimes: In the limit in which the impactor is significantly smaller than the target ($R_i<< R_t$), the solutions are self-similar and the shock velocity at a fixed point on the target scale as $m_i^{2/3}$, where $m_i$ is the mass of the impactor. In addition, the ground velocities follow a universal profile given by $v_g/v_i=(14.2x^2-25.3x+11.3)/(x^2-2.5x+1.9) +2\ln{R_i/R_t}$, where $x=\sin\left(\theta/2\right)$, $\theta$ is the latitude on the target measured from the impact site, and $v_g$ and $v_i$ are the ground velocity and impact velocity, respectively. In contrast, in the limit in which the impactor is comparable to the size of the target ($R_i \sim R_t$), we find that shock velocities decline with the mass of the impactor significantly more weakly than $m_i^{2/3}$. We use the resulting surface velocity profiles to calculate the atmospheric mass loss for a large range of impactor masses and impact velocities and apply them to the Kepler-36 system and the Moon forming impact. Finally, we present and generalise our results in terms of the $v_g/v_i$ and the impactor to target size ratio ($R_i/R_t$) such that they can easily be applied to other collision scenarios.

The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution.

The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond.

Systems of close-in super-Earths display striking diversity in planetary bulk density and composition. Giant impacts are expected to play a role in the formation of many of these worlds. Previous works, focused on the mechanical shock caused by a giant impact, have shown that these impacts can eject large fractions of the planetary envelope, offering a partial explanation for the observed spread in exoplanet compositions. Here, we examine the thermal consequences of giant impacts, and show that the atmospheric loss caused by these effects can significantly exceed that caused by mechanical shocks for hydrogen-helium (H/He) envelopes. When a giant impact occurs, part of the impact energy is converted into thermal energy, heating the rocky core and the envelope. We find that the ensuing thermal expansion of the envelope can lead to a period of sustained, rapid mass loss through a Parker wind, resulting in the partial or complete erosion of the H/He envelope. The fraction of the envelope lost depends on the planet's orbital distance from its host star and its initial thermal state, and hence age. Planets closer to their host stars are more susceptible to thermal atmospheric loss triggered by impacts than ones on wider orbits. Similarly, younger planets, with rocky cores which are still hot and molten from formation, suffer greater atmospheric loss. This is especially interesting because giant impacts are expected to occur $10{-}100~\mathrm{Myr}$ after formation. For planets where the thermal energy of the core is much greater than the envelope energy, the impactor mass required for significant atmospheric removal is $M_\mathrm{imp} / M_p \sim \mu / \mu_c \sim 0.1$, approximately the ratio of the heat capacities of the envelope and core. When the envelope energy dominates the total energy budget, complete loss can occur when the impactor mass is comparable to the envelope mass.

Recent observations identify a valley in the radius distribution of small exoplanets, with planets in the range $1.5-2.0\,{\rm R}_{\oplus}$ significantly less common than somewhat smaller or larger planets. This valley may suggest a bimodal population of rocky planets that are either engulfed by massive gas envelopes that significantly enlarge their radius, or do not have detectable atmospheres at all. One explanation of such a bimodal distribution is atmospheric erosion by high-energy stellar photons. We investigate an alternative mechanism: the luminosity of the cooling rocky core, which can completely erode light envelopes while preserving heavy ones, produces a deficit of intermediate sized planets. We evolve planetary populations that are derived from observations using a simple analytical prescription, accounting self-consistently for envelope accretion, cooling and mass loss, and demonstrate that core-powered mass loss naturally reproduces the observed radius distribution, regardless of the high-energy incident flux. Observations of planets around different stellar types may distinguish between photoevaporation, which is powered by the high-energy tail of the stellar radiation, and core-powered mass loss, which depends on the bolometric flux through the planet's equilibrium temperature that sets both its cooling and mass-loss rates.

Super-Earths are the most abundant planets known to date and are characterized by having sizes between that of Earth and Neptune, typical orbital periods of less than 100 days and gaseous envelopes that are often massive enough to significantly contribute to the planet's overall radius. Furthermore, super-Earths regularly appear in tightly-packed multiple-planet systems, but resonant configurations in such systems are rare. This chapters summarizes current super-Earth formation theories. It starts from the formation of rocky cores and subsequent accretion of gaseous envelopes. We follow the thermal evolution of newly formed super-Earths and discuss their atmospheric mass loss due to disk dispersal, photoevaporation, core-cooling and collisions. We conclude with a comparison of observations and theoretical predictions, highlighting that even super-Earths that appear as barren rocky cores today likely formed with primordial hydrogen and helium envelopes and discuss some paths forward for the future.

The measurement of an exoplanet's oblateness and obliquity provides insights into the planet's internal structure and formation history. Previous work using small differences in the shape of the transit light curve has been moderately successful, but was hampered by the small signal and extreme photometric precision required. The measurement of changes in transit depth, caused by the spin precession of an oblate planet, was proposed as an alternative method. Here, we present the first attempt to measure these changes. Using Kepler photometry, we examined the brown dwarf Kepler-39b and the warm Saturn Kepler-427b. We could not reliably constrain the oblateness of Kepler-39b. We find transit depth variations for Kepler-427b at $90.1\%$ significance ($1.65\sigma$) consistent with a precession period of $P_\mathrm{prec} = 5.45^{+0.46}_{-0.37}~\mathrm{years}$ and an oblateness, $f~=~0.19^{+0.32}_{-0.16}$. This oblateness is comparable to Solar System gas giants, and would raise questions about the dynamics and tidal synchronization of Kepler-427b.

We present 1-second cadence observations of M32 (NGC221) with the CHIMERA instrument at the Hale 200-inch telescope of the Palomar Observatory. Using field stars as a baseline for relative photometry, we are able to construct a light curve of the nucleus in the g-prime and r-prime band with 1sigma=36 milli-mag photometric stability. We derive a temporal power spectrum for the nucleus and find no evidence for a time-variable signal above the noise as would be expected if the nuclear black hole were accreting gas. Thus, we are unable to constrain the spin of the black hole although future work will use this powerful instrument to target more actively accreting black holes. Given the black hole mass of (2.5+/-0.5)*10^6 Msun inferred from stellar kinematics, the absence of a contribution from a nuclear time-variable signal places an upper limit on the accretion rate which is 4.6*10^{-8} of the Eddington rate, a factor of two more stringent than past upper limits from HST. The low mass of the black hole despite the high stellar density suggests that the gas liberated by stellar interactions was primarily at early cosmic times when the low-mass black hole had a small Eddington luminosity. This is at least partly driven by a top-heavy stellar initial mass function at early cosmic times which is an efficient producer of stellar mass black holes. The implication is that supermassive black holes likely arise from seeds formed through the coalescence of 3-100 Msun mass black holes that then accrete gas produced through stellar interaction processes.

A commonly noted feature of the population of multi-planet extrasolar systems is the rarity of planet pairs in low-order mean-motion resonances. We revisit the physics of resonance capture via convergent disk-driven migration. We point out that for planet spacings typical of stable configurations for Kepler systems, the planets can routinely maintain a small but nonzero eccentricity due to gravitational perturbations from their neighbors. Together with the upper limit on the migration rate needed for capture, the finite eccentricity can make resonance capture difficult or impossible in Sun-like systems for planets smaller than ~Neptune-sized. This mass limit on efficient capture is broadly consistent with observed exoplanet pairs that have mass determinations: of pairs with the heavier planet exterior to the lighter planet -- which would have been undergoing convergent migration in their disks -- those in or nearly in resonance are much more likely to have total mass greater than two Neptune masses than to have smaller masses. The agreement suggests that the observed paucity of resonant pairs around sun-like stars may simply arise from a small resonance capture probability for lower-mass planets. Planet pairs that thereby avoid resonance capture are much less likely to collide in an eventual close approach than to simply migrate past one another to become a divergently migrating pair with the lighter planet exterior. For systems around M stars we expect resonant pairs to be much more common, since there the minimum mass threshhold for efficient capture is about an Earth mass.

The majority of discovered exoplanetary systems harbour a new class of planets, bodies typically several times more massive than Earth but orbiting their host stars well inside the orbit of Mercury. The origin of these close-in super-Earths and mini-Neptunes is a major unanswered question in planet formation. Unlike Earth, whose atmosphere contains $<10^{-6}$ its total mass, a large fraction of close-in planets have significant gaseous envelopes, containing $1 -10\%$ or more of their total mass. It has been proposed that these close-in planets formed in situ either by delivery of $50-100M_{\oplus}$ of rocky material to the inner disc, or in a disc enhanced relative to the MMSN. In both cases, final assembly of the planets occurs by giant impacts (GIs). Here we test the viability of these scenarios. We show that atmospheres accreted by isolation masses are small ($10^{-3}-10^{-2}$ the core mass) and that atmospheric mass-loss during GIs is significant, with typical post-GI atmospheres that are $8 \times 10^{-4}$ the core mass. Such values are consistent with terrestrial planet atmospheres but more than an order of magnitude below atmospheric masses of $1-10\%$ inferred for many close-in exoplanets. In the most optimistic scenario with no core luminosity, post-GI envelope accretion from a depleted gas disc yields atmospheric masses that are several per cent the core mass. If the gravitational potential energy due to the last mass doubling of the planet by GIs is released over the disc dissipation time-scale as core luminosity, then envelope masses are reduced by about an order of magnitude. Finally we show that radial drift time-scales due to gas drag for many isolation masses are shorter than typical disc lifetimes. Given these challenges, we conclude that most observed close-in planets with envelopes larger than several per cent likely formed at larger separations from their host stars.

Some recently discovered short-period Earth to Neptune sized exoplanets (super Earths) have low observed mean densities which can only be explained by voluminous gaseous atmospheres. Here, we study the conditions allowing the accretion and retention of such atmospheres. We self-consistently couple the nebular gas accretion onto rocky cores and the subsequent evolution of gas envelopes following the dispersal of the protoplanetary disk. Specifically, we address mass-loss due to both photo-evaporation and cooling of the planet. We find that planets shed their outer layers (dozens of percents in mass) following the disk's dispersal (even without photo-evaporation), and their atmospheres shrink in a few Myr to a thickness comparable to the radius of the underlying rocky core. At this stage, atmospheres containing less particles than the core (equivalently, lighter than a few % of the planet's mass) can be blown away by heat coming from the cooling core, while heavier atmospheres cool and contract on a timescale of Gyr at most. By relating the mass-loss timescale to the accretion time, we analytically identify a Goldilocks region in the mass-temperature plane in which low-density super Earths can be found: planets have to be massive and cold enough to accrete and retain their atmospheres, while not too massive or cold, such that they do not enter runaway accretion and become gas giants (Jupiters). We compare our results to the observed super-Earth population and find that low-density planets are indeed concentrated in the theoretically allowed region. Our analytical and intuitive model can be used to investigate possible super-Earth formation scenarios.

The Caltech HIgh-speed Multi-color camERA (CHIMERA) is a new instrument that has been developed for use at the prime focus of the Hale 200-inch telescope. Simultaneous optical imaging in two bands is enabled by a dichroic beam splitter centered at 567 nm, with Sloan u' and g' bands available on the blue arm and Sloan r', i' and z_s' bands available on the red arm. Additional narrow-band filters will also become available as required. An Electron Multiplying CCD (EMCCD) detector is employed for both optical channels, each capable of simultaneously delivering sub-electron effective read noise under multiplication gain and frame rates of up to 26 fps full frame (several 1000 fps windowed), over a fully corrected 5 x 5 arcmin field of view. CHIMERA was primarily developed to enable the characterization of the size distribution of sub-km Kuiper Belt Objects via stellar occultation, a science case that motivates the frame-rate, the simultaneous multi-color imaging and the wide field of view of the instrument. In addition, it also has unique capability in the detection of faint near-Earth asteroids and will be used for the monitoring of short duration transient and periodic sources, particularly those discovered by the intermediate Palomar Transient Factory (iPTF), and the upcoming Zwicky Transient Facility (ZTF).

Although current sensitivity limits are such that true Solar System analogs remain challenging to detect, numerous planetary systems have been discovered that are very different from our own Solar System. The majority of systems harbor a new class of planets, bodies that are typically several times more massive than the Earth but that orbit their host stars well inside the orbit of Mercury. These planets frequently show evidence for large Hydrogen and Helium envelopes containing several percent of the planet's mass and display a large diversity in mean densities. Here we show that this wide range can be achieved by one or two late giant impacts, which are frequently needed to achieve long-term orbital stability in multiple planet systems once the gas disk has disappeared. We demonstrate using hydrodynamical simulations that a single collision between similarly sized exoplanets can easily reduce the envelope-to-core-mass ratio by a factor of two and show that this leads to a corresponding increase in the observed mean density by factors of 2-3. In addition we investigate how envelope-mass-loss depends on envelope mass, planet radius, semi-major axis, and the mass distribution inside the envelope. We propose that a small number of giant impacts may be responsible for the large observed spread in mean densities, especially for multiple-planet systems containing planets with very different densities and which have not been significantly sculpted by photo evaporation.

Models are developed to simulate lightcurves of stars dimmed by transiting exoplanets with and without rings. These models are then applied to \textit{Kepler} photometry to search for planetary rings in a sample of 21 exoplanets, mostly hot Jupiters, chosen to offer the best observational opportunity for discovering potential rings. We also examine what kinds of rings might be expected for these planets, in terms of both size and orientation, based on arguments involving the host planet's equilibrium temperature, its likely obliquities, and the formation and stability of possible ring systems. Finding no evidence for rings, for each of the 21 studied planets it is determined on an observational basis which potential rings can be rejected out of a representative set of fiducial rings, varying in both size and orientation. For 12 of the 21 planets, we determined that Saturn-like rings could be ruled out for at least certain orientations. Additionally, the detectability of rings is studied, and it is found that ringed planets with small obliquities (roughly $5^{\circ}-10^{\circ}$) can yield large signals, which is encouraging for future work, since such small obliquities are expected for hot Jupiters.

Recent observations by the {\it Kepler} space telescope have led to the discovery of more than 4000 exoplanet candidates consisting of many systems with Earth- to Neptune-sized objects that reside well inside the orbit of Mercury, around their respective host stars. How and where these close-in planets formed is one of the major unanswered questions in planet formation. Here we calculate the required disk masses for {\it in situ} formation of the {\it Kepler} planets. We find that, if close-in planets formed as {\it isolation masses}, then standard gas-to-dust ratios yield corresponding gas disks that are gravitationally unstable for a significant fraction of systems, ruling out such a scenario. We show that the maximum width of a planet's accretion region in the absence of any migration is $2 v_{esc}/\Omega$, where $v_{esc}$ is the escape velocity of the planet and $\Omega$ the Keplerian frequency and use it to calculate the required disk masses for {\it in situ} formation with giant impacts. Even with giant impacts, formation without migration requires disk surface densities in solids at semi-major axes less than 0.1~AU of $10^3-10^5 \rm{~g~cm^{-2}}$ implying typical enhancements above the minimum-mass solar nebular (MMSN) by at least a factor of 20. Corresponding gas disks are below, but not far from, the gravitational stability limit. In contrast, formation beyond a few AU is consistent with MMSN disk masses. This suggests that migration of either solids or fully assembled planets is likely to have played a major role in the formation of close-in super-Earths and mini-Neptunes.

We quantify the atmospheric mass loss during planet formation by examining the contributions to atmospheric loss from both giant impacts and planetesimal accretion. Giant impacts cause global motion of the ground. Using analytic self-similar solutions and full numerical integrations we find (for isothermal atmospheres with adiabatic index ($\gamma=5/3$) that the local atmospheric mass loss fraction for ground velocities $v_g < 0.25 v_{esc}$ is given by $\chi_{loss}=(1.71 v_g/v_{esc})^{4.9}$, where $v_{esc}$ is the escape velocity from the target. Yet, the global atmospheric mass loss is a weaker function of the impactor velocity $v_{Imp}$ and mass $m_{Imp}$ and given by $X_{loss} ~ 0.4x+1.4x^2-0.8x^3$ (isothermal atmosphere) and $X_{loss} ~ 0.4x+1.8x^2-1.2x^3$ (adiabatic atmosphere), where $x=(v_{Imp}m/v_{esc}M)$. Atmospheric mass loss due to planetesimal impacts proceeds in two different regimes: 1) Large enough impactors $m > \sqrt{2} \rho_0 (\pi h R)^{3/2}$ (25~km for the current Earth), are able to eject all the atmosphere above the tangent plane of the impact site, which is $h/2R$ of the whole atmosphere, where $h$, $R$ and $\rho_0$ are the atmospheric scale height, radius of the target, and its atmospheric density at the ground. 2) Smaller impactors, but above $m>4 \pi \rho_0 h^3$ (1~km for the current Earth) are only able to eject a fraction of the atmospheric mass above the tangent plane. We find that the most efficient impactors (per unit impactor mass) for atmospheric loss are planetesimals just above that lower limit and that the current atmosphere of the Earth could have resulted from an equilibrium between atmospheric erosion and volatile delivery to the atmosphere from planetesimals. We conclude that planetesimal impacts are likely to have played a major role in atmospheric mass loss over the formation history of the terrestrial planets. (Abridged)

We assess the multi-planet systems discovered by the Kepler satellite in terms of current ideas about orbital migration and eccentricity damping due to planet-disk interactions. Our primary focus is on mean motion resonances. Only a few percent of planet pairs are in close proximity to a resonance. However, predicted migration rates (parameterized by $\tau_n=n/{|\dot n|}$) imply that during convergent migration most planets would have been captured into first order resonances. Eccentricity damping (parameterized by $\tau_e=e/{|\dot e|}$) offers a plausible resolution. Estimates suggest $\tau_e/\tau_n\sim (h/a)^2\sim 10^{-2}$, where $h/a$ is the ratio of disk thickness to radius. Together, eccentricity damping and orbital migration give rise to an equilibrium eccentricity, $e_{eq}\sim(\tau_e/\tau_n)^{1/2}$. Capture is permanent provided $e_{eq}\lesssim \mu^{1/3}$, where $\mu$ denotes the planet to star mass ratio. But for $e_{eq}\gtrsim \mu^{1/3}$, capture is only temporary because librations around equilibrium are overstable and lead to passage through resonance on timescale $\tau_e$. Most Kepler planet pairs have $e_{eq}>\mu^{1/3}$. Since $\tau_n>> \tau_e$ is the timescale for migration between neighboring resonances, only a modest percentage of pairs end up trapped in resonances after the disk disappears. Planet pairs close to a mean motion resonance typically exhibit period ratios 1-2% larger than those for exact resonance. The direction of this shift undoubtedly reflects the same asymmetry that requires convergent migration for resonance capture. Permanent resonance capture at these separations from exact resonance would demand $\mu (\tau_n/\tau_e)^{1/2}\gtrsim 0.01$, a value that estimates of $\mu$ from transit data and $(\tau_e/\tau_n)^{1/2}$ from theory are insufficient to match. Plausible alternatives involve eccentricity damping during or after disk dispersal. (Abridged)

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.

We show, by comparing observations with theoretical models, that the observed Kuiper Belt size distribution is well matched by coagulation models, which start from an initial planetesimal population with radii of about 1km, and subsequent collisional evolution. We find that the observed size distribution for R > 30km has not been modified by collisional evolution over the age of the solar system, and that the size distribution below R ~ 30km has been modified by collisions and that its slope is well matched by collisional evolution models that use published strength laws. We investigate in detail the resulting size distribution of bodies ranging from 0.01km to 30km and find that its slope changes several times as a function of radius before approaching the expected value for an equilibrium collisional cascade of material strength dominated bodies for R < 0.1km. Compared to a single power law size distribution that would span the whole range from 0.01km to 30km, we find a strong deficit of bodies around R ~10km and a strong excess of bodies around 2km. This deficit and excess are caused by an excess mass in small planetesimals in the km size range that was left over from the runaway growth phase and that leaves a signature in the size distribution that is not erased after 4.5 Gyrs of collisional evolution. Observations of the small KBO size distribution can therefore test if large KBOs grew as a result of runaway growth and constrain the initial planetesimal sizes. We find that results from recent KBO occultation surveys and the observed KBO size distribution can be best matched by an initial planetesimal population that contained about equal mass per logarithmic mass bin in bodies ranging from 0.4km to 4km in radius. In addition, we find that we cannot match the observed KBO size distribution if most of the planetesimal mass was contained in bodies that were 10km in radius or larger.

We present here the analysis of about 19,500 new star hours of low ecliptic latitude observations (|b| < 20 deg) obtained by the Hubble Space Telescope's FGS over a time span of more than nine years; which is an addition to the 12,000 star hours previously analyzed by Schlichting et al. (2009). Our search for stellar occultations by small Kuiper belt objects (KBOs) yielded one new candidate event corresponding to a body with a 530 +/-70m radius at a distance of about 40AU. Using bootstrap simulations, we estimate a probability of approx 5%, that this event is due to random statistical fluctuations within the new data set. Combining this new event with the single KBO occultation reported by Schlichting et al. (2009) we arrive at the following results: 1) The ecliptic latitudes of 6.6 deg and 14.4 deg of the two events are consistent with the observed inclination distribution of larger, 100km-sized KBOs. 2) Assuming that small, sub-km sized KBOs have the same ecliptic latitude distribution as their larger counterparts, we find an ecliptic surface density of KBOs with radii larger than 250m of N(r>250m) = 1.1^{+1.5}_{-0.7} x 10^7 deg^{-2}; if sub-km sized KBOs have instead a uniform ecliptic latitude distribution for -20 deg < b< 20 deg then N(r>250m) = 4.4^{+5.8}_{-2.8} x 10^6 deg^{-2}. This is the best measurement of the surface density of sub-km sized KBOs to date. 3) Assuming the KBO size distribution can be well described by a single power law given by N(>r) \propto r^{1-q}, where N(>r) is the number of KBOs with radii greater than r, and q is the power law index, we find q=3.8+/-0.2 for a KBO ecliptic latitude distribution that follows the observed distribution for larger, 100-km sized KBOs. 4) Regardless of the exact power law, our results suggest that small KBOs are numerous enough to satisfy the required supply rate for the Jupiter family comets. (Abridged)

The final stage of terrestrial planet formation consists of the cleanup of residual planetesimals after the giant impact phase. Dynamically, a residual planetesimal population is needed to damp the high eccentricities of the terrestrial planets after the giant impact stage. Geochemically, highly siderophile element (HSE) abundance patterns inferred for the terrestrial planets and the Moon suggest that a total of about 0.01 M_Earth of chondritic material was delivered as `late veneer' by planetesimals to the terrestrial planets after the end of giant impacts. Here we combine these two independent lines of evidence for a leftover population of planetesimals and show that: 1) A residual planetesimal population containing 0.01 M_Earth is able to damp the eccentricities of the terrestrial planets after giant impacts to their observed values. 2) At the same time, this planetesimal population can account for the observed relative amounts of late veneer added to the Earth, Moon and Mars provided that the majority of the late veneer was delivered by small planetesimals with radii <10m. These small planetesimal sizes are required to ensure efficient damping of the planetesimal's velocity dispersion by mutual collisions, which in turn ensures that the planets' accretion cross sections are significantly enhanced by gravitational focusing above their geometric values. Specifically we find, in the limit that the relative velocity between the terrestrial planets and the planetesimals is significantly less than the terrestrial planets' escape velocities, that gravitational focusing yields an accretion ratio Earth/Mars~17, which agrees well with the accretion ratio inferred from HSEs of 12-23. For the Earth-Moon system, we find an accretion ratio of ~200, which is consistent with estimates of 150-700 derived from HSE abundances that include the lunar crust as well as mantle component. (Abridged)

The standard theoretical treatment of collisional cascades derives a steady-state size distribution assuming a single constant velocity dispersion for all bodies regardless of size. Here we relax this assumption and solve self-consistently for the bodies' steady-state size and size-dependent velocity distributions. Specifically, we account for viscous stirring, dynamical friction, and collisional damping of the bodies' random velocities in addition to the mass conservation requirement typically applied to find the size distribution in a steady-state cascade. The resulting size distributions are significantly steeper than those derived without velocity evolution. For example, accounting self-consistently for the velocities can change the standard q=3.5 power-law index of the Dohnanyi (1969) differential size spectrum to an index as large as q=4. Similarly, for bodies held together by their own gravity, the corresponding power-law index range 2.88<q<3.14 of Pan & Sari (2005) can steepen to values as large as q=3.26. Our velocity results allow quantitative predictions of the bodies' scale heights as a function of size. Together with our predictions, observations of the scale heights for different sized bodies for the Kuiper belt, the asteroid belt, and extrasolar debris disks may constrain the mass and number of large bodies stirring the cascade as well as the colliding bodies' internal strengths.

We discuss the nature of rings that may exist around extrasolar planets. Taking the general properties of rings around the gas giants in the Solar System, we infer the likely properties of rings around exoplanets that reside inside the ice line. Due to their proximity to their host star, rings around such exoplanets must primarily consist of rocky materials. However, we find that despite the higher densities of rock compared to ice, most of the observed extrasolar planets with reliable radii measurements have sufficiently large Roche radii to support rings. For the currently known transiting extrasolar planets, Poynting-Robertson drag is not effective in significantly altering the dynamics of individual ring particles over a time span of $10^8$ years provided that they exceed about 1 m in size. In addition, we show that significantly smaller ring particles can exist in optically thick rings, for which we find typical ring lifetimes ranging from a few times $10^6$ to a few times $10^9$ years. Most interestingly, we find that many of the rings could have nontrivial Laplacian planes due to the increased effects of the orbital quadrupole caused by the exoplanets' proximity to their host star, allowing a constraint on the $J_2$ of extrasolar planets from ring observations. This is particular exciting, since a planet's $J_2$ reveals information about its interior structure. Furthermore, measurements of an exoplanet's oblateness and of its $J_2$, from warped rings, would together place limits on its spin period. Based on the constraints that we have derived for extrasolar rings, we anticipate that the best candidates for ring detections will come from transit observations by the Kepler spacecraft of extrasolar planets with semi-major axes $\sim 0.1$ AU and larger.

Approximately 10-20% of all Kuiper belt objects (KBOs) occupy mean-motion resonances with Neptune. This dynamical configuration likely resulted from resonance capture as Neptune migrated outward during the late stages of planet formation. The details of Neptune's planetesimal-driven migration, including its radial extent and the concurrent eccentricity evolution of the planet, are the subject of considerable debate. Two qualitatively different proposals for resonance capture have been proposed--migration-induced capture driven by smooth outward evolution of Neptune's orbit and chaotic capture driven by damping of the planet's eccentricity near its current semi-major axis. We demonstrate that the distribution of comparable-mass, wide-separation binaries occupying resonant orbits can differentiate between these two scenarios. If migration-induced capture occurred, this fraction records information about the formation locations of different populations of KBOs. Chaotic capture, in contrast, randomizes the orbits of bodies as they are placed in resonance. In particular, migration-induced capture produces the following signatures. The 2:1 resonance should contain a dynamically cold component, with inclinations less than 5-10 degrees, having a binary fraction comparable to that among cold classical KBOs. If the 3:2 resonance also hosts a cold component, its binary fraction should be 20-30% lower than in the cold classical belt. Among cold 2:1 (and if present 3:2) KBOs, objects with eccentricities e < 0.2 should have a binary fraction ~20% larger than those with e > 0.2. Searches for cold components in the binary fractions of resonant KBOs are currently practical. The additional migration-generated trends described here may be distinguished with objects discovered by LSST. (Abstract abridged.)

Runway growth is an important stage in planet formation during which large protoplanets form, while most of the initial mass remains in small planetesimals. The amount of mass converted into large protoplanets and their resulting size distribution are not well understood. Here, we use analytic work, that we confirm by coagulation simulations, to describe runaway growth and the corresponding evolution of the velocity dispersion. We find that runaway growth proceeds as follows: Initially all the mass resides in small planetesimals, with mass surface density \sigma, and large protoplanets start to form by accreting small planetesimals. This growth continues until growth by merging large protoplanets becomes comparable to growth by planetesimal accretion. This condition sets in when \Sigma/\sigma ~\alpha^{3/4} ~ 10^{-3}, where \Sigma is the mass surface density in protoplanets in a given logarithmic mass interval and \alpha is the ratio of the size of a body to its Hill radius. From then on, protoplanetary growth and the evolution of the velocity dispersion become self-similar and \Sigma remains roughly constant, since an increase in \Sigma by accretion of small planetesimals is balanced by a decrease due to merging with large protoplanets. We show that this growth leads to a protoplanet size distribution given by N(>R) \propto R^{-3} where N(>R) is the number of objects with radii greater than R (i.e., a differential power-law index of 4). We apply our results to the Kuiper Belt, which is a relic of runaway growth. Our results successfully match the observed Kuiper belt size distribution, they illuminate the physical processes that shaped it and explain the total mass that is present in large Kuiper belt objects (KBOs) today. This work suggests that the current mass in large KBOs is primordial and that it has not been significantly depleted. Abridged

The Kuiper belt is a remnant of the primordial Solar System. Measurements of its size distribution constrain its accretion and collisional history, and the importance of material strength of Kuiper belt objects (KBOs). Small, sub-km sized, KBOs elude direct detection, but the signature of their occultations of background stars should be detectable. Observations at both optical and X-ray wavelengths claim to have detected such occultations, but their implied KBO abundances are inconsistent with each other and far exceed theoretical expectations. Here, we report an analysis of archival data that reveals an occultation by a body with a 500 m radius at a distance of 45 AU. The probability of this event to occur due to random statistical fluctuations within our data set is about 2%. Our survey yields a surface density of KBOs with radii larger than 250 m of 2.1^{+4.8}_{-1.7} x 10^7 deg^{-2}, ruling out inferred surface densities from previous claimed detections by more than 5 sigma. The fact that we detected only one event, firmly shows a deficit of sub-km sized KBOs compared to a population extrapolated from objects with r>50 km. This implies that sub-km sized KBOs are undergoing collisional erosion, just like debris disks observed around other stars.

Recently, the first collisional family was discovered in the Kuiper belt. The parent body of this family, Haumea, is one of the largest objects in the Kuiper belt and is orbited by two satellites. It has been proposed that the Haumea family was created from dispersed fragments that resulted from a giant impact. This proposed origin of the Haumea family is however in conflict with the observed velocity dispersion between the family members (\sim 140 m/s) which is significantly less than the escape velocity from Haumea's surface (\sim 900 m/s). In this paper we propose a different formation scenario for Haumea's collisional family. In our scenario the family members are ejected while in orbit around Haumea. This scenario, therefore, gives naturally rise to a lower velocity dispersion among the family members than expected from direct ejection from Haumea's surface. In our scenario Haumea's giant impact forms a single moon that tidally evolves outward until it suffers a destructive collision from which the family is created. We show that this formation scenario yields a velocity dispersion of \sim 190m/s among the family members which is in good agreement with the observations. The probability for Haumea's initial giant impact in todays Kuiper belt is less than 10^{-3}. In our scenario, however, Haumea's giant impact can occur before the excitation of the Kuiper belt and the ejection of the family members afterwards. This has the advantage that one can preserve the dynamical coherence of the family and explain Haumea's original giant impact, which is several orders of magnitude more likely to have occurred in the primordial dynamically cold Kuiper belt compared to the dynamically excited Kuiper belt today. Abridged

The evidence is mounting that star formation necessarily involves planet formation. We clearly have a vested interest in finding other Earths but a true understanding of planet formation requires completing the census and mapping planetary architecture in all its grandeur and diversity. Here, we show that a 2000-star survey undertaken with SIM Lite will uniquely probe planets around B-A-F stars, bright and binary stars and white dwarfs. In addition, we show that the high precision of SIM Lite allows us to gain unique insights into planet formation via accurate measurements of mutual inclinations.

With the discovery of Kuiper Belt binaries that have wide separations and roughly equal masses new theories were proposed to explain their formation. Two formation scenarios were suggested by Goldreich and collaborators: In the first, dynamical friction that is generated by a sea of small bodies enables a transient binary to become bound ($L^2s$ mechanism); in the second, a transient binary gets bound by an encounter with a third body ($L^3$ mechanism). We show that these different binary formation scenarios leave their own unique signatures in the relative abundance of prograde to retrograde binary orbits. This signature is due to stable retrograde orbits that exist much further out in the Hill sphere than prograde orbits. It provides an excellent opportunity to distinguish between the different binary formation scenarios observationally. We predict that if binary formation proceeded while sub-Hill velocities prevailed, the vast majority of all comparable mass ratio binaries have retrograde orbits. This dominance of retrograde binary orbits is a result of binary formation via the $L^2s$ mechanism, or any other mechanism that dissipates energy in a smooth and gradual manner. For super-Hill velocities binary formation proceeds via the $L^3$ mechanism which produces a roughly equal number of prograde and retrograde binaries. These predictions assume that subsequent orbital evolution due to dynamical friction and dynamical stirring of the Kuiper belt did not alter the sense of the binary orbit after formation.

The discovery that a substantial fraction of Kuiper Belt objects (KBOs) exists in binaries with wide separations and roughly equal masses, has motivated a variety of new theories explaining their formation. Goldreich et al. (2002) proposed two formation scenarios: In the first, a transient binary is formed, which becomes bound with the aid of dynamical friction from the sea of small bodies (L^2s mechanism); in the second, a binary is formed by three body gravitational deflection (L^3 mechanism). Here, we accurately calculate the L^2s and L^3 formation rates for sub-Hill velocities. While the L^2s formation rate is close to previous order of magnitude estimates, the L^3 formation rate is about a factor of 4 smaller. For sub-Hill KBO velocities (v << v_H) the ratio of the L^3 to the L^2s formation rate is 0.05 (v/v_H) independent of the small bodies' velocity dispersion, their surface density or their mutual collisions. For Super-Hill velocities (v >> v_H) the L^3 mechanism dominates over the L^2s mechanism. Binary formation via the L^3 mechanism competes with binary destruction by passing bodies. Given sufficient time, a statistical equilibrium abundance of binaries forms. We show that the frequency of long-lived transient binaries drops exponentially with the system's lifetime and that such transient binaries are not important for binary formation via the L^3 mechanism, contrary to Lee et al. (2007). For the L^2s mechanism we find that the typical time, transient binaries must last, to form Kuiper Belt binaries (KBBs) for a given strength of dynamical friction, D, increases only logarithmically with D. Longevity of transient binaries only becomes important for very weak dynamical friction (i.e. D \lesssim 0.002) and is most likely not crucial for KBB formation.

Planetesimal accretion during planet formation is usually treated as collisionless. Such accretion from a uniform and dynamically cold disk predicts protoplanets with slow retrograde rotation. However, if the building blocks of protoplanets, planetesimals, are small, of order of a meter in size, then they are likely to collide within the protoplanet's sphere of gravitational influence, creating a prograde accretion disk around the protoplanet. The accretion of such a disk results in the formation of protoplanets spinning in the prograde sense with the maximal spin rate allowed before centrifugal forces break them apart. As a result of semi-collisional accretion, the final spin of a planet after giant impacts is not completely random but is biased toward prograde rotation. The eventual accretion of the remaining planetesimals in the post giant-impact phase might again be in the semi-collisional regime and delivers a significant amount of additional prograde angular momentum to the terrestrial planets. We suggest that in our Solar System, semi-collisional accretion gave rise to the preference for prograde rotation observed in the terrestrial planets and perhaps the largest asteroids.

We examine the growth of eccentricities of a population of particles with initially circular orbits around a central massive body. Successive encounters between pairs of particles increase the eccentricities in the disk on average. As long as the epicyclic motions of the particles are small compared to the shearing motion between Keplerian orbits, there is no preferred scale for the eccentricities. The simplification due to this self-similarity allows us to find an analytic form for the distribution function; full numerical integrations of a disk with 200 planetesimals verify our analytical self-similar distribution. The shape of this non-equilibrium profile is identical to the equilibrium profile of a shear-dominated population whose mutual excitations are balanced by dynamical friction or Epstein gas drag.