We obtained spectra of the pre-main sequence star AU Microscopii during a transit of its Neptune-sized planet to investigate its orbit and atmosphere. We used the high-dispersion near-infrared spectrograph IRD on the Subaru telescope to detect the Doppler "shadow" from the planet and constrain the projected stellar obliquity; we find that the stellar spin and orbit are aligned ($\lambda=-4.7_{-6.4}^{+6.8}$ degrees), suggesting that the planet formed and possibly migrated within the protoplanetary disk. We use both the IRD data and spectra obtained with NIRSPEC on Keck-2 to search for absorption in the 1083\,nm line of metastable triplet He\,I by the planet's atmosphere and place an upper limit for the equivalent width of 3.7 m\AA{} at 99\,\% confidence. From this limit, we constrain the mass loss rate of the planetary atmosphere to $<0.15-0.45$ \mearth\,Gyr$^{-1}$, comparable to the rate predicted by an energy-limited escape calculation and in the range required to test models of XUV-driven atmospheric escape from close-in young "Neptunes".

Growing protoplanets experience a number of impacts during the accretion stage. A large impactor hits the surface of a protoplanet and produces a magma ocean, where the impactor's iron emulsifies and experiences metal-silicate equilibration with the mantle of the protoplanet while it descends towards the base of the magma ocean. This process repeatedly occurs and determines the chemical compositions of both mantle and core. The partitioning is controlled by parameters such as the equilibration pressure and temperature, which are often associated with or assumed to be proportional to the pressure and temperature at the base of the magma ocean. The pressure and temperature depend on both the depth and shape of a magma ocean because a spatially confined magma ocean, namely, a melt pool, can have a larger equilibrium pressure than a radially uniform (global) magma ocean even if their melt volumes are the same. Here, we develop scaling laws for (1) the total internal energy gain due to an impact, and (2) the heat distribution within the mantle based on more than 100 smoothed particle hydrodynamic (SPH) simulations. We use Legendre polynomials to describe these scaling laws and determine their coefficients by linear regression, minimizing the error between our model and SPH simulations. The input parameters are the impact angle, total mass, impact velocity, and impactor-to-total mass ratio. We find that the internal energy gain by a large impact is well characterized by the summation of the kinetic energy and accretional potential energy release as a function of the impact angle. We determine that the equilibrium pressure at the base of a melt pool can be higher (by $10-50 \%$) than those obtained from conventional radially-uniform global magma ocean models. These melt scaling laws are publicly available on GitHub (https://github.com/mikinakajima/MeltScalingLaw).

We obtained high-resolution infrared spectroscopy and short-cadence photometry of the 600-800 Myr Praesepe star K2-100 during transits of its 1.67-day planet. This Neptune-size object, discovered by the NASA K2 mission, is an interloper in the "desert" of planets with similar radii on short period orbits. Our observations can be used to understand its origin and evolution by constraining the orbital eccentricity by transit fitting, measuring the spin-orbit obliquity by the Rossiter-McLaughlin effect, and detecting any extended, escaping hydrogen-helium envelope with the 10830A line of neutral helium in the 2s3S triplet state. Transit photometry with 1-min cadence was obtained by the K2 satellite during Campaign 18 and transit spectra were obtained with the IRD spectrograph on the Subaru telescope. While the elevated activity of K2-100 prevented us from detecting the Rossiter-McLaughlin effect, the new photometry combined with revised stellar parameters allowed us to constrain the eccentricity to e < 0.15/0.28 with 90%/99% confidence. We modeled atmospheric escape as an isothermal, spherically symmetric Parker wind, with photochemistry driven by UV radiation that we estimate by combining the observed spectrum of the active Sun with calibrations from observations of K2-100 and similar young stars in the nearby Hyades cluster. Our non-detection (<5.7mA) of a transit-associated He I line limits mass loss of a solar-composition atmosphere through a T<10000K wind to <0.3Me/Gyr. Either K2-100b is an exceptional desert-dwelling planet, or its mass loss is occurring at a lower rate over a longer interval, consistent with a core accretion-powered scenario for escape.

The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the subsequent formation of residual condensates with strongly sub-solar Al/Si and Mg/Si ratios. The mixing of residual condensates with different amounts of material with solar refractory/Si element ratios explains the Al/Si and Mg/Si values of non-carbonaceous chondrites. To match quantitatively the observed ratios, we find that the first-planetesimals should have accreted when the disk temperature was ~1,330-1,400 K depending on pressure and assuming a solar C/O ratio of the disk. We discuss how this model relates to our current understanding of disk evolution, grain dynamics, and planetesimal formation. We also extend the discussion to moderately volatile elements (e.g., Na), explaining how it may be possible that the depletion of these elements in non-carbonaceous chondrites is correlated with the depletion of refractory elements (e.g., Al). Extending the analysis to Cr, we find evidence for a higher than solar C/O ratio in the protosolar disk's gas when/where condensation from a fractionated gas occurred. Finally, we discuss the possibility that the supra-solar Al/Si and Mg/Si ratios of the Earth are due to the accretion of ~40% of the mass of our planet from the first-generation of refractory-rich planetesimals.

Recent observations of rocky super-Earths have revealed an apparent wider distribution of Fe/Mg ratios, or core to mantle ratios, than the planets in our Solar System. This study aims to understand how much of the chemical diversity in the super-Earth population can arise from giant impacts during planetary formation. Planet formation simulations have only recently begun to treat collisions more realistically in an attempt to replicate the planets in our Solar System. We investigate planet formation more generally by simulating the formation of rocky super-Earths with varying initial conditions using a version of SyMBA, a gravitational N-body code, that incorporates realistic collisions. We track the maximum plausible change in composition after each impact. The final planets span a range of Fe/Mg ratios similar to the Solar System planets, but do not completely match the distribution in super-Earth data. We only form a few planets with minor iron-depletion, suggesting other mechanisms are at work. The most iron-rich planets have a lower Fe/Mg ratio than Mercury, and are less enriched than planets such as Kepler-100b. This indicates that further work on our understanding of planet formation and further improvement of precision of mass and radius measurements are required to explain planets at the extremes of this Fe/Mg distribution.

In an effort to measure the Rossiter-McLaughlin effect for the TRAPPIST-1 system, we performed high-resolution spectroscopy during transits of planets e, f, and b. The spectra were obtained with the InfraRed Doppler spectrograph on the Subaru 8.2-m telescope, and were supplemented with simultaneous photometry obtained with a 1-m telescope of the Las Cumbres Observatory Global Telescope. By analyzing the anomalous radial velocities, we found the projected stellar obliquity to be $\lambda=1\pm 28$ degrees under the assumption that the three planets have coplanar orbits, although we caution that the radial-velocity data show correlated noise of unknown origin. We also sought evidence for the expected deformations of the stellar absorption lines, and thereby detected the "Doppler shadow" of planet b with a false alarm probability of $1.7\,\%$. The joint analysis of the observed residual cross-correlation map including the three transits gave $\lambda=19_{-15}^{+13}$ degrees. These results indicate that the the TRAPPIST-1 star is not strongly misaligned with the common orbital plane of the planets, although further observations are encouraged to verify this conclusion.

NASA's Kepler mission discovered $\sim700$ planets in multi-planet systems containing 3 or more transiting bodies, many of which are super-Earths and mini-Neptunes in compact configurations. Using $N$-body simulations, we examine the in situ, final stage assembly of multi-planet systems via the collisional accretion of protoplanets. Our initial conditions are constructed using a subset of the Kepler 5-planet systems as templates. Two different prescriptions for treating planetary collisions are adopted. The simulations address numerous questions: do the results depend on the accretion prescription?; do the resulting systems resemble the Kepler systems, and do they reproduce the observed distribution of planetary multiplicities when synthetically observed?; do collisions lead to significant modification of protoplanet compositions, or to stripping of gaseous envelopes?; do the eccentricity distributions agree with those inferred for the Kepler planets? We find the accretion prescription is unimportant in determining the outcomes. The final planetary systems look broadly similar to the Kepler templates adopted, but the observed distributions of planetary multiplicities or eccentricities are not reproduced, because scattering does not excite the systems sufficiently. In addition, we find that $\sim 1$% of our final systems contain a co-orbital planet pair in horseshoe or tadpole orbits. Post-processing the collision outcomes suggests they would not significantly change the ice fractions of initially ice-rich protoplanets, but significant stripping of gaseous envelopes appears likely. Hence, it may be difficult to reconcile the observation that many low mass Kepler planets have H/He envelopes with an in situ formation scenario that involves giant impacts after dispersal of the gas disc.

We examine 141 N-body simulations of terrestrial planet late-stage accretion that use the Grand Tack scenario, coupling the collisional results with a hafnium-tungsten (Hf-W) isotopic evolution model. Accretion in the Grand Tack scenario results in faster planet formation than classical accretion models because of higher planetesimal surface density induced by a migrating Jupiter. Planetary embryos that grow rapidly experience radiogenic ingrowth of mantle $^{182}$W that is inconsistent with the measured terrestrial composition, unless much of the tungsten is removed by an impactor core that mixes thoroughly with the target mantle. For physically Earth-like surviving planets, we find that the fraction of equilibrating impactor core $k_\text{core} \geq 0.6$ is required to produce results agreeing with observed terrestrial tungsten anomalies (assuming equilibration with relatively large volumes of target mantle material; smaller equilibrating mantle volumes would require even larger $k_\text{core}$). This requirement of substantial core re-equilibration may be difficult to reconcile with fluid dynamical predictions and hydrocode simulations of mixing during large impacts, and hence this result does not favor the rapid planet building that results from Grand Tack accretion.

The Keck Planet Imager and Characterizer (KPIC) is an upgrade to the Keck II adaptive optics system enabling high contrast imaging and high-resolution spectroscopic characterization of giant exoplanets in the mid-infrared (2-5 microns). The KPIC instrument will be developed in phases. Phase I entails the installation of an infrared pyramid wavefront sensor (PyWFS) based on a fast, low-noise SAPHIRA IR-APD array. The ultra-sensitive infrared PyWFS will enable high contrast studies of infant exoplanets around cool, red, and/or obscured targets in star forming regions. In addition, the light downstream of the PyWFS will be coupled into an array of single-mode fibers with the aid of an active fiber injection unit (FIU). In turn, these fibers route light to Keck's high-resolution infrared spectrograph NIRSPEC, so that high dispersion coronagraphy (HDC) can be implemented for the first time. HDC optimally pairs high contrast imaging and high-resolution spectroscopy allowing detailed characterization of exoplanet atmospheres, including molecular composition, spin measurements, and Doppler imaging. Here we provide an overview of the instrument, its science scope, and report on recent results from on-sky commissioning of Phase I. The instrument design and techniques developed will be key for more advanced instrument concepts needed for the extremely large telescopes of the future.

We explore the hypothesis that the population of Martian Trojans is the result of a balance between the production of new asteroids ("YORPlets") through the YORP effect and their eventual escape from the Trojan clouds through Yarkovsky-driven orbital evolution. Our principal observables are: (5261) Eureka, its family of 8 asteroids and the family-less Trojans (101429) 1998 VF31 & (121514) 1999 UJ7. We model the population evolution as a birth-death random process and assume it is in a steady state. We then simulate the discovery of Trojans to-date and find that family members of 101429 and 121514, if they exist, are intrinsically more difficult to detect than Eureka's. Their non-discovery can be used as evidence of their non-existence only under the assumption that their brightness relative to the parent asteroid is similar to that in the Eureka family. To find out how efficiently Mars Trojans are lost from the Trojan clouds due to the Yarkovsky effect, we carry out dynamical simulations of test particles originating from these parent bodies. We find that objects originating from Eureka and 121514 begin escaping after ~1 Gyr, but that those from 101429 are already lost by that time, probably due to that asteroid's proximity to an eccentricity-type secular resonance within Mars's co-orbital region. This is the likely cause behind the absence of Trojans in the orbital vicinity of 101429. In contrast, the solitary status of 121514 points to an intrinsic inability of the asteroid to produce YORPlets during the most recent ~20% of the solar system's history, a finding potentially related to 121514's present, low angular momentum rotational state, unless the Eureka family formed rapidly during a single fission event.

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

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.

Wavefront sensing and control are important for enabling one of the key advantages of using large apertures, namely higher angular resolutions. Pyramid wavefront sensors are becoming commonplace in new instrument designs owing to their superior sensitivity. However, one remaining roadblock to their widespread use is the fabrication of the pyramidal optic. This complex optic is challenging to fabricate due to the pyramid tip, where four planes need to intersect in a single point. Thus far, only a handful of these have been produced due to the low yields and long lead times. To address this, we present an alternative implementation of the pyramid wavefront sensor that relies on two roof prisms instead. Such prisms are easy and inexpensive to source. We demonstrate the successful operation of the roof prism pyramid wavefront sensor on a 8-m class telescope, at visible and near infrared wavelengths ---for the first time using a SAPHIRA HgCdTe detector without modulation for a laboratory demonstration---, and elucidate how this sensor can be used more widely on wavefront control test benches and instruments.

We report on initial results from the largest infrared AO direct imaging survey searching for wide orbit (>100 AU) massive exoplanets and brown dwarfs as companions around young nearby stars using Robo-AO at the 2.1-m telescope on Kitt Peak, Arizona. The occurrence rates of these rare substellar companions are critical to furthering our understanding of the origin of planetary-mass companions on wide orbits. The observing efficiency of Robo-AO allows us to conduct a survey an order of magnitude larger than previously possible. We commissioned a low-noise high-speed SAPHIRA near-infrared camera to conduct this survey and report on its sensitivity, performance, and data reduction process.

We discuss some of the unique details of the operation and behavior of Leonardo SAPHIRA detectors, particularly in relation to their usage for adaptive optics wavefront sensing. SAPHIRA detectors are 320$\times$256@24 $\mu$m pixel HgCdTe linear avalanche photodiode arrays and are sensitive to 0.8-2.5 $\mu m$ light. SAPHIRA arrays permit global or line-by-line resets, of the entire detector or just subarrays of it, and the order in which pixels are reset and read enable several readout schemes. We discuss three readout modes, the benefits, drawbacks, and noise sources of each, and the observational modes for which each is optimal. We describe the ability of the detector to read subarrays for increased frame rates, and finally clarify the differences between the avalanche gain (which is user-adjustable) and the charge gain (which is not).

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.

The timeline of the lunar bombardment in the first Gy of the Solar System remains unclear. Some basin-forming impacts occurred 3.9-3.7Gy ago. Many other basins formed before, but their exact ages are not precisely known. There are two possible interpretations of the data: in the cataclysm scenario there was a surge in the impact rate approximately 3.9Gy ago, while in the accretion tail scenario the lunar bombardment declined since the era of planet formation and the latest basins formed in its tail-end. Here, we revisit the work of Morbidelli et al.(2012) that examined which scenario could be compatible with both the lunar crater record in the 3-4Gy period and the abundance of highly siderophile elements (HSE) in the lunar mantle. We use updated numerical simulations of the fluxes of impactors. Under the traditional assumption that the HSEs track the total amount of material accreted by the Moon since its formation, we conclude that only the cataclysm scenario can explain the data. The cataclysm should have started ~3.95Gy ago. However we show that HSEs could have been sequestered from the lunar mantle due to iron sulfide exsolution during magma ocean crystallization, followed by mantle overturn. Based on the hypothesis that the lunar magma ocean crystallized about 100-150My after Moon formation, and therefore that HSEs accumulated in the lunar mantle only after this time, we show that the bombardment in the 3-4Gy period can be explained in the accretion tail scenario. This hypothesis would also explain why the Moon appears so depleted in HSEs relative to the Earth. We also extend our analysis of the cataclysm and accretion tail scenarios to the case of Mars. The accretion tail scenario requires a global resurfacing event on Mars ~4.4Gy ago, possibly associated with the formation of the Borealis basin, and it is consistent with the HSE budget of the planet.

A suite of science instruments is critical to any high contrast imaging facility, as it defines the science capabilities and observing modes available. SCExAO uses a modular approach which allows for state-of-the-art visitor modules to be tested within an observatory environment on an 8-m class telescope. This allows for rapid prototyping of new and innovative imaging techniques that otherwise take much longer in traditional instrument design. With the aim of maturing science modules for an advanced high contrast imager on an giant segmented mirror telescopes (GSMTs) that will be capable of imaging terrestrial planets, we offer an overview and status update on the various science modules currently under test within the SCExAO instrument.

Earth possesses a persistent, internally-generated magnetic field, whereas no trace of a dynamo has been detected on Venus, at present or in the past, although a high surface temperature and recent resurfacing events may have removed paleomagnetic evidence. Whether or not a terrestrial body can sustain an internally generated magnetic field by convection inside its metallic fluid core is determined in part by its initial thermodynamic state and its compositional structure, both of which are in turn set by the processes of accretion and differentiation. Here we show that the cores of Earth- and Venus-like planets should grow with stable compositional stratification unless disturbed by late energetic impacts. They do so because higher abundances of light elements are incorporated into the liquid metal that sinks to form the core as the temperatures and pressures of metal-silicate equilibration increase during accretion. We model this process and determine that this establishes a stable stratification that resists convection and inhibits the onset of a geodynamo. However, if a late energetic impact occurs, it could mechanically stir the core creating a single homogenous region within which a long-lasting geodynamo would operate. While Earth's accretion has been punctuated by a late giant impact with likely enough energy to mix the core (e.g. the impact that formed the Moon), we hypothesize that the accretion of Venus is characterized by the absence of such energetic giant impacts and the preservation of its primordial stratifications.

Seven of the nine known Mars Trojan asteroids belong to an orbital cluster named after its largest member 5261 Eureka. Eureka is likely the progenitor of the whole cluster, which formed at least 1 Gyr ago. It was suggested that the thermal YORP effect spun-up Eureka resulting with fragments being ejected by the rotational-fission mechanism. Eureka's spectrum exhibits a broad and deep absorption band around 1 {\mu}m, indicating an olivine-rich composition. Here we show evidence that the Trojan Eureka cluster progenitor could have originated as impact debris excavated from the Martian mantle. We present new near-infrared observations of two Trojans (311999 2007 NS2 and 385250 2001 DH47) and find that both exhibit an olivine-rich reflectance spectrum similar to Eureka's. These measurements confirm that the progenitor of the cluster has an achondritic composition. Olivine-rich reflectance spectra are rare amongst asteroids but are seen around the largest basins on Mars. They are also consistent with some Martian meteorites (e.g. Chassigny), and with the material comprising much of the Martian mantle. Using numerical simulations, we show that the Mars Trojans are more likely to be impact ejecta from Mars than captured olivine-rich asteroids transported from the main belt. This result links directly specific asteroids to debris from the forming planets.

Highly siderophile elements (HSEs) are strongly depleted in the bulk silicate Earth (BSE) but are present in near-chondritic relative abundances. The conventional explanation is that the HSEs were stripped from the mantle by the segregation of metal during core formation but were added back in near-chondritic proportions by late accretion, after core formation had ceased. Here we show that metal-silicate equilibration and segregation during Earth's core formation actually increased HSE mantle concentrations because HSE partition coefficients are relatively low at the high pressures of core formation within Earth. The pervasive exsolution and segregation of iron sulfide liquid from silicate liquid (the "Hadean matte") stripped magma oceans of HSEs during cooling and crystallization, before late accretion, and resulted in slightly suprachondritic palladium/iridium and ruthenium/iridium ratios.

From the results of a comprehensive asteroid population evolution model, we conclude that the YORP-induced rotational fission hypothesis can be consistent with the observed population statistics of small asteroids in the main belt including binaries and contact binaries. The foundation of this model is the asteroid rotation model of Marzari et al. (2011), which incorporates both the YORP effect and collisional evolution. This work adds to that model the rotational fission hypothesis and the binary evolution model of Jacobson & Scheeres (2011). The asteroid population evolution model is highly constrained by these and other previous works, and therefore it has only two significant free parameters: the ratio of low to high mass ratio binaries formed after rotational fission events and the mean strength of the binary YORP (BYORP) effect. We successfully reproduce characteristic statistics of the small asteroid population: the binary fraction, the fast binary fraction, steady-state mass ratio fraction and the contact binary fraction. We find that in order for the model to best match observations, rotational fission produces high mass ratio ($> 0.2$) binary components with four to eight times the frequency as low mass ratio ($< 0.2$) components, where the mass ratio is the mass of the secondary component divided by the mass of the primary component. This is consistent with post-rotational fission binary system mass ratio being drawn from either a flat or a positive and shallow distribution, since the high mass ratio bin is four times the size of the low mass ratio bin; this is in contrast to the observed steady-state binary mass ratio, which has a negative and steep distribution. This can be understood in the context of the BYORP-tidal equilibrium hypothesis, which predicts that low mass ratio binaries survive for a significantly longer period of time than high mass ratio systems.

Because of the high energies involved, giant impacts that occur during planetary accretion cause large degrees of melting. The depth of melting in the target body after each collision determines the pressure and temperature conditions of metal-silicate equilibration and thus geochemical fractionation that results from core-mantle differentiation. The accretional collisions involved in forming the terrestrial planets of the inner Solar System have been calculated by previous studies using N-body accretion simulations. Here we use the output from such simulations to determine the volumes of melt produced and thus the pressure and temperature conditions of metal-silicate equilibration, after each impact, as Earth-like planets accrete. For these calculations a parametrised melting model is used that takes impact velocity, impact angle and the respective masses of the impacting bodies into account. The evolution of metal-silicate equilibration pressures (as defined by evolving magma ocean depths) during Earth's accretion depends strongly on the lifetime of impact-generated magma oceans compared to the time interval between large impacts. In addition, such results depend on starting parameters in the N-body simulations, such as the number and initial mass of embryos. Thus, there is the potential for combining the results, such as those presented here, with multistage core formation models to better constrain the accretional history of the Earth.

Earth and Moon are shown here to be composed of oxygen isotope reservoirs that are indistinguishable, with a difference in {\Delta}"17O of -1 +/- 5ppm (2se). Based on these data and our new planet formation simulations that include a realistic model for oxygen isotopic reservoirs, our results favor vigorous mixing during the giant impact and therefore a high-energy high- angular-momentum impact. The results indicate that the late veneer impactors had an average {\Delta}"17O within approximately 1 per mil of the terrestrial value, suggesting that these impactors were water rich.

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.

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

Rotationally fissioned asteroids produce unbound daughter asteroids that have very similar heliocentric orbits. Backward integration of their current heliocentric orbits provides an age of closest proximity that can be used to date the rotational fission event. Most asteroid pairs follow a predicted theoretical relationship between the primary spin period and the mass ratio of the two pair members that is a direct consequence of the YORP-induced rotational fission hypothesis. If the progenitor asteroid has strength, asteroid pairs may have high mass ratios with possibly fast rotating primaries. However, secondary fission leaves the originally predicted trend unaltered. We also describe the characteristics of pair members produced by four alternative routes from a rotational fission event to an asteroid pair. Unlike direct formation from the event itself, the age of closest proximity of these pairs cannot generally be used to date the rotational fission event since considerable time may have passed.

Sublimating gas molecules scatter off of the surface of an icy body in the same manner as photons. This means that for every photon-driven body force, there should be a sublimation-driven analogue that affects icy bodies. Thermal photons emitted from the surfaces of asymmetrically shaped bodies in the Solar System generate net torques that change the spin rates of these bodies over time. The long-term averaging of this torque is called the YORP effect. Here we propose a sublimation-driven analogue to the YORP effect (SYORP), in which sublimating gas molecules emitted from the surfaces of icy bodies also generate net torques on the bodies. However, sublimating molecules carry momentum away from the body at a rate ~10^4-10^5 greater than thermal photons, resulting in much greater body torques. While previous studies of sublimative torques focused on emissions from highly localized sources on the surfaces of Jupiter Family Comet nuclei, SYORP applies to non-localized emissions across the entire body, which likely dominates sublimation-drive torques on small icy chunks and Dynamically Young Comets, and can therefore be applied without high-resolution spacecraft observations of their surfaces. Instead, we repurpose the equations of the YORP effect to account for sublimation-driven torques. We show how an SYORP-driven mechanism best matches observations of the Sun-oriented lineaments (striae) of comet tails, whose formation mechanism has remained enigmatic for decades. The SYORP effect explains why striae are observed between near-perihelion and ~1 AU from the Sun for comets with perihelia less than 0.6 AU, and solves longstanding problems with moving enough material into the cometary tail to form visible striae. We show that SYORP can form striae that match those of Comet West, and produce a power-law size frequency distribution of their parent chunks with an index of -1.4 (-1.1- 2).

In the past decade, the number of known binary near-Earth asteroids has more than quadrupled and the number of known large main belt asteroids with satellites has doubled. Half a dozen triple asteroids have been discovered, and the previously unrecognized populations of asteroid pairs and small main belt binaries have been identified. The current observational evidence confirms that small (<20 km) binaries form by rotational fission and establishes that the YORP effect powers the spin-up process. A unifying paradigm based on rotational fission and post-fission dynamics can explain the formation of small binaries, triples, and pairs. Large (>20 km) binaries with small satellites are most likely created during large collisions.

Infrared avalanche photodiode arrays represent a panacea for many branches of astronomy by enabling extremely low-noise, high-speed and even photon-counting measurements at near-infrared wavelengths. We recently demonstrated the use of an early engineering-grade infrared avalanche photodiode array that achieves a correlated double sampling read noise of 0.73 e- in the lab, and a total noise of 2.52 e- on sky, and supports simultaneous high-speed imaging and tip-tilt wavefront sensing with the Robo-AO visible-light laser adaptive optics system at the Palomar Observatory 1.5-m telescope. We report here on the improved image quality achieved simultaneously at visible and infrared wavelengths by using the array as part of an image stabilization control-loop with adaptive-optics sharpened guide stars. We also discuss a newly enabled survey of nearby late M-dwarf multiplicity as well as future uses of this technology in other adaptive optics and high-contrast imaging applications.

Satellites of asteroids have been discovered in nearly every known small body population, and a remarkable aspect of the known satellites is the diversity of their properties. They tell a story of vast differences in formation and evolution mechanisms that act as a function of size, distance from the Sun, and the properties of their nebular environment at the beginning of Solar System history and their dynamical environment over the next 4.5 Gyr. The mere existence of these systems provides a laboratory to study numerous types of physical processes acting on asteroids and their dynamics provide a valuable probe of their physical properties otherwise possible only with spacecraft. Advances in understanding the formation and evolution of binary systems have been assisted by: 1) the growing catalog of known systems, increasing from 33 to nearly 250 between the Merline et al. (2002) Asteroids III chapter and now, 2) the detailed study and long-term monitoring of individual systems such as 1999 KW4 and 1996 FG3, 3) the discovery of new binary system morphologies and triple systems, 4) and the discovery of unbound systems that appear to be end-states of binary dynamical evolutionary paths. Specifically for small bodies (diameter smaller than 10 km), these observations and discoveries have motivated theoretical work finding that thermal forces can efficiently drive the rotational disruption of small asteroids. Long-term monitoring has allowed studies to constrain the system's dynamical evolution by the combination of tides, thermal forces and rigid body physics. The outliers and split pairs have pushed the theoretical work to explore a wide range of evolutionary end-states.

The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk. What was the origin of this dichotomy in the mass distribution of embryos/cores? We show in this paper that the classic processes of runaway and oligarchic growth from a disk of planetesimals cannot explain this dichotomy, even if the original surface density of solids increased at the snowline. Instead, the accretion of drifting pebbles by embryos and cores can explain the dichotomy, provided that some assumptions hold true. We propose that the mass-flow of pebbles is two-times lower and the characteristic size of the pebbles is approximately ten times smaller within the snowline than beyond the snowline (respectively at heliocentric distance $r<r_{ice}$ and $r>r_{ice}$, where $r_{ice}$ is the snowline heliocentric distance), due to ice sublimation and the splitting of icy pebbles into a collection of chondrule-size silicate grains. In this case, objects of original sub-lunar mass would grow at drastically different rates in the two regions of the disk. Within the snowline these bodies would reach approximately the mass of Mars while beyond the snowline they would grow to $\sim 20$ Earth masses. The results may change quantitatively with changes to the assumed parameters, but the establishment of a clear dichotomy in the mass distribution of protoplanets appears robust, provided that there is enough turbulence in the disk to prevent the sedimentation of the silicate grains into a very thin layer.

The formation of the Earth's core is a consequence of planetary accretion and processes in the Earth's interior. The mechanical process of planetary differentiation is likely to occur in large, if not global, magma oceans created by the collisions of planetary embryos. Metal-silicate segregation in magma oceans occurs rapidly and efficiently unlike grain scale percolation according to laboratory experiments and calculations. Geochemical models of the core formation process as planetary accretion proceeds are becoming increasingly realistic. Single stage and continuous core formation models have evolved into multi-stage models that are couple to the output of dynamical models of the giant impact phase of planet formation. The models that are most successful in matching the chemical composition of the Earth's mantle, based on experimentally-derived element partition coefficients, show that the temperature and pressure of metal-silicate equilibration must increase as a function of time and mass accreted and so must the oxygen fugacity of the equilibrating material. The latter can occur if silicon partitions into the core and through the late delivery of oxidized material. Coupled dynamical accretion and multi-stage core formation models predict the evolving mantle and core compositions of all the terrestrial planets simultaneously and also place strong constraints on the bulk compositions and oxidation states of primitive bodies in the protoplanetary disk.

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.

The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with the nebula gas lead to remarkable new accretion mechanisms for the formation of planetesimals and the growth of planetary embryos. Many of the popular models for the later stages of planet formation are presented. The classical models with the giant planets on fixed orbits are not consistent with the known history of the Solar System, fail to create a high Earth/Mars mass ratio, and, in many cases, are also internally inconsistent. The successful Grand Tack model creates a small Mars, a wet Earth, a realistic asteroid belt and the mass-orbit structure of the terrestrial planets. In the Grand Tack scenario, growth curves for Earth most closely match a Weibull model. The feeding zones, which determine the compositions of Earth and Venus follow a particular pattern determined by Jupiter, while the feeding zones of Mars and Theia, the last giant impactor on Earth, appear to randomly sample the terrestrial disk. The late accreted mass samples the disk nearly evenly.

In order to test planetary accretion and differentiation scenarios, we integrated a multistage core-mantle differentiation model with N-body accretion simulations. Impacts between embryos and planetesimals result in magma ocean formation and episodes of core formation. The core formation model combines rigorous chemical mass balance with metal-silicate element partitioning data. The primary constraint on the combined model is the composition of the Earth's primitive mantle, the composition of the Martian mantle, and the mass fractions of the metallic cores of Earth and Mars. The model is refined by least squares minimization with up to five fitting parameters that consist of the metal-silicate equilibrium pressure and 1-4 parameters that define the starting compositions of primitive bodies. This integrated model has been applied to 6 Grand Tack simulations. Investigations of a broad parameter space indicate that: accretion of Earth was heterogeneous, metal-silicate equilibration pressures increase as accretion progresses and are 60-70% of core-mantle boundary pressures at the time of each impact, and a large fraction (70-100%) of the metal of impactor cores equilibrates with a small fraction of the silicate mantles of protoplanets during each core formation event. Acceptable fits to the Earth's mantle composition are obtained only when bodies that originated close to the Sun, at <0.9-1.2 AU, are highly reduced and those beyond this distance are increasingly oxidized. The FeO content of the Martian mantle depends critically on the heliocentric distance at which the Mars-forming embryo originated. Finally, the Earth's core is predicted to contain 8-9 wt% silicon, 2-4 wt% oxygen and 10-60 ppm hydrogen, whereas the Martian core is predicted to contain low concentrations (<1 wt%) of Si and O.

Observations with radar, photometric and direct imaging techniques have discovered that multiple asteroid systems can be divided clearly into a handful of different morphologies, and recently, the discovery of small unbound asteroid systems called asteroid pairs have revolutionized the study of small asteroid systems. Simultaneously, new theoretical advances have demonstrated that solar radiation dictates the evolution of small asteroids with strong implications for asteroid internal structure. We review our current understanding of how small asteroid systems evolve and point to the future.

Using our photometric observations taken between 1996 and 2013 and other published data, we derived properties of the binary near-Earth asteroid (175706) 1996 FG3 including new measurements constraining evolution of the mutual orbit with potential consequences for the entire binary asteroid population. We also refined previously determined values of parameters of both components, making 1996 FG3 one of the most well understood binary asteroid systems. We determined the orbital vector with a substantially greater accuracy than before and we also placed constraints on a stability of the orbit. Specifically, the ecliptic longitude and latitude of the orbital pole are 266{\deg} and -83{\deg}, respectively, with the mean radius of the uncertainty area of 4{\deg}, and the orbital period is 16.1508 +/- 0.0002 h (all quoted uncertainties correspond to 3sigma). We looked for a quadratic drift of the mean anomaly of the satellite and obtained a value of 0.04 +/- 0.20 deg/yr^2, i.e., consistent with zero. The drift is substantially lower than predicted by the pure binary YORP (BYORP) theory of McMahon and Scheeres (McMahon, J., Scheeres, D. [2010]. Icarus 209, 494-509) and it is consistent with the theory of an equilibrium between BYORP and tidal torques for synchronous binary asteroids as proposed by Jacobson and Scheeres (Jacobson, S.A., Scheeres, D. [2011]. ApJ Letters, 736, L19). Based on the assumption of equilibrium, we derived a ratio of the quality factor and tidal Love number of Q/k = 2.4 x 10^5 uncertain by a factor of five. We also derived a product of the rigidity and quality factor of mu Q = 1.3 x 10^7 Pa using the theory that assumes an elastic response of the asteroid material to the tidal forces. This very low value indicates that the primary of 1996 FG3 is a 'rubble pile', and it also calls for a re-thinking of the tidal energy dissipation in close asteroid binary systems.

Although discs of dust and gas have been observed orbiting white dwarfs, the origin of this circumstellar matter is uncertain. We hypothesize that the in-situ breakup of small bodies such as asteroids spun to fission during the giant branch phases of stellar evolution provides an important contribution to this debris. The YORP effect, which arises from radiation pressure, accelerates the spin rate of asymmetric asteroids, which can eventually shear themselves apart. This pressure is maintained and enhanced around dying stars because the outward push of an asteroid due to stellar mass loss is insignificant compared to the resulting stellar luminosity increase. Consequently, giant star radiation will destroy nearly all bodies with radii in the range 100 m - 10 km that survive their parent star's main sequence lifetime within a distance of about 7 au; smaller bodies are spun apart to their strongest, competent components. This estimate is conservative, and would increase for highly asymmetric shapes or incorporation of the inward drag due to giant star stellar wind. The resulting debris field, which could extend to thousands of au, may be perturbed by remnant planetary systems to reproduce the observed dusty and gaseous discs which accompany polluted white dwarfs.

The prototype Robo-AO system at the Palomar Observatory 1.5-m telescope is the world's first fully automated laser adaptive optics instrument. Scientific operations commenced in June 2012 and more than 12,000 observations have since been performed at the ~0.12" visible-light diffraction limit. Two new infrared cameras providing high-speed tip-tilt sensing and a 2' field-of-view will be integrated in 2014. In addition to a Robo-AO clone for the 2-m IGO and the natural guide star variant KAPAO at the 1-m Table Mountain telescope, a second generation of facility-class Robo-AO systems are in development for the 2.2-m University of Hawai'i and 3-m IRTF telescopes which will provide higher Strehl ratios, sharper imaging, ~0.07", and correction to {\lambda} = 400 nm.

We present conclusions from a large number of N-body simulations of the giant impact phase of terrestrial planet formation. We focus on new results obtained from the recently proposed Grand Tack model, which couples the gas-driven migration of giant planets to the accretion of the terrestrial planets. The giant impact phase follows the oligarchic growth phase, which builds a bi-modal mass distribution within the disc of embryos and planetesimals. By varying the ratio of the total mass in the embryo population to the total mass in the planetesimal population and the mass of the individual embryos, we explore how different disc conditions control the final planets. The total mass ratio of embryos to planetesimals controls the timing of the last giant (Moon forming) impact and its violence. The initial embryo mass sets the size of the lunar impactor and the growth rate of Mars. After comparing our simulated outcomes with the actual orbits of the terrestrial planets (angular momentum deficit, mass concentration) and taking into account independent geochemical constraints on the mass accreted by the Earth after the Moon forming event and on the timescale for the growth of Mars, we conclude that the protoplanetary disc at the beginning of the giant impact phase must have had most of its mass in Mars-sized embryos and only a small fraction of the total disc mass in the planetesimal population. From this, we infer that the Moon forming event occurred between $\sim$60 and $\sim$130 My after the formation of the first solids, and was caused most likely by an object with a mass similar to that of Mars.

We present a model of near-Earth asteroid (NEA) rotational fission and ensuing dynamics that describes the creation of synchronous binaries and all other observed NEA systems including: doubly synchronous binaries, high- e binaries, ternary systems, and contact binaries. Our model only presupposes the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, "rubble pile" asteroid geophysics, and gravitational interactions. The YORP effect torques a "rubble pile" asteroid until the asteroid reaches its fission spin limit and the components enter orbit about each other (Scheeres, D.J. [2007]. Icarus 189, 370-385). Non-spherical gravitational potentials couple the spin states to the orbit state and chaotically drive the system towards the observed asteroid classes along two evolutionary tracks primarily distinguished by mass ratio. Related to this is a new binary process termed secondary fission - the secondary asteroid of the binary system is rotationally accelerated via gravitational torques until it fissions, thus creating a chaotic ternary system. The initially chaotic binary can be stabilized to create a synchronous binary by components of the fissioned secondary asteroid impacting the primary asteroid, solar gravitational perturbations, and mutual body tides. These results emphasize the importance of the initial component size distribution and configuration within the parent asteroid. NEAs may go through multiple binary cycles and many YORP-induced rotational fissions during their approximately 10 Myr lifetime in the inner Solar System. Rotational fission and the ensuing dynamics are responsible for all NEA systems including the most commonly observed synchronous binaries.

The size distribution of small asteroids in the Main Belt is assumed to be determined by an equilibrium between the creation of new bodies out of the impact debris of larger asteroids and the destruction of small asteroids by collisions with smaller projectiles. However, for a diameter less than 6 km we find that YORP-induced rotational disruption significantly contributes to the erosion even exceeding the effects of collisional fragmentation. Including this additional grinding mechanism in a collision evolution model for the asteroid belt, we generate size-frequency distributions from either an accretional (Weidenschilling, 2011) or an "Asteroids were born big" (Morbidelli, 2009) initial size-frequency distribution that are consistent with observations reported in Gladman et al. (2009). Rotational disruption is a new mechanism that must be included in all future collisional evolution models of asteroids.

We propose and analyze a new mechanism for the formation of the wide asynchronous binary population. These binary asteroids have wide semi-major axes relative to most near-Earth and Main Belt asteroid systems. Confirmed members have rapidly rotating primaries and satellites that are not tidally locked. Previously suggested formation mechanisms from impact ejecta, planetary flybys and directly from rotational fission events cannot satisfy all of the observations. The newly hypothesized mechanism works as follows: (i) these systems are formed from rotational fission, (ii) their satellites are tidally locked, (iii) their orbits are expanded by the BYORP effect, (iv) their satellites de-synchronize due to the adiabatic invariance between the libration of the secondary and the mutual orbit, and (v) the secondary avoids resynchronization due to the the YORP effect. This seemingly complex chain of events is a natural pathway for binaries with satellites that have particular shapes, which define the BYORP effect torque that acts on the system. After detailing the theory, we analyze each of the wide asynchronous binary members and candidates to assess their most likely formation mechanism. Finally, we suggest possible future observations to check and constrain our hypothesis.

Synchronous binary asteroids may exist in a long-term stable equilibrium, where the opposing torques from mutual body tides and the binary YORP (BYORP) effect cancel. Interior of this equilibrium, mutual body tides are stronger than the BYORP effect and the mutual orbit semi-major axis expands to the equilibrium; outside of the equilibrium, the BYORP effect dominates the evolution and the system semi-major axis will contract to the equilibrium. If the observed population of small (0.1 - 10 km diameter) synchronous binaries are in static configurations that are no longer evolving, then this would be confirmed by a null result in the observational tests for the BYORP effect. The confirmed existence of this equilibrium combined with a shape model of the secondary of the system enables the direct study of asteroid geophysics through the tidal theory. The observed synchronous asteroid population cannot exist in this equilibrium if described by the canonical "monolithic" geophysical model. The "rubble pile" geophysical model proposed by \citet{Goldreich2009} is sufficient, however it predicts a tidal Love number directly proportional to the radius of the asteroid, while the best fit to the data predicts a tidal Love number inversely proportional to the radius. This deviation from the canonical and \citet{Goldreich2009} models motivates future study of asteroid geophysics. Ongoing BYORP detection campaigns will determine whether these systems are in an equilibrium, and future determination of secondary shapes will allow direct determination of asteroid geophysical parameters.

Asteroid pairs sharing similar heliocentric orbits were found recently. Backward integrations of their orbits indicated that they separated gently with low relative velocities, but did not provide additional insight into their formation mechanism. A previously hypothesized rotational fission process4 may explain their formation - critical predictions are that the mass ratios are less than about 0.2 and, as the mass ratio approaches this upper limit, the spin period of the larger body becomes long. Here we report photometric observations of a sample of asteroid pairs revealing that primaries of pairs with mass ratios much less than 0.2 rotate rapidly, near their critical fission frequency. As the mass ratio approaches 0.2, the primary period grows long. This occurs as the total energy of the system approaches zero requiring the asteroid pair to extract an increasing fraction of energy from the primary's spin in order to escape. We do not find asteroid pairs with mass ratios larger than 0.2. Rotationally fissioned systems beyond this limit have insufficient energy to disrupt. We conclude that asteroid pairs are formed by the rotational fission of a parent asteroid into a proto-binary system which subsequently disrupts under its own internal system dynamics soon after formation.

Microlensing events observed in the line of sight toward the LMC indicate that a significant fraction of the mass of the dark halo of the Galaxy is probably composed of white dwarfs. In addition, white dwarf sequences have now be observed in the HR diagrams of several globular clusters. Because of the unavailability of white dwarf atmospheres for Teff < 4000K, cooling time scales for white dwarfs older than ~ 10 Gyr are very uncertain. Moreover, the identification of a MACHO white dwarf population by direct observation depends on a knowledge of the colors and bolometric corrections of very-cool white dwarfs. In this paper we present the first detailed model atmospheres and spectra of very cool hydrogen white dwarfs for Teff < 4000K. We include the latest description of the opacities of hydrogen and significantly, we introduce a non-ideal equation of state in the atmosphere calculation. We find that due to strong absorption from H_2 in the infrared, very old white dwarfs are brightest in the V, R, and I bands, and we confirm that they become bluer in most color indices as they cool below Teff ~ 3500K.