Search

The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary. DART revealed a largely oblate and ellipsoidal shape of Dimorphos before the impact, while the post-impact observations suggest that Dimorphos now has a prolate shape. Here we add those data points to the known properties of small binary asteroids and propose new paradigms of the radiative binary YORP (BYORP) effect as well as tidal dissipation in small binaries. We find that relatively spheroidal bodies like Dimorphos made of small debris may experience a weaker and more size-dependent BYORP effect than previously thought. This could explain the observed values of period drift in several well-characterized binaries. We also propose that energy dissipation in small binaries is dominated by relatively brief episodes of large-scale movement of (likely surface) materials, rather than long-term steady-state tidal dissipation. We propose that one such episode was triggered on Dimorphos by the DART impact. Depending on the longevity of this high-dissipation regime, it is possible that Dimorphos will be more dynamically relaxed in time for the Hera mission than it was in the weeks following the impact., Comment: Accepted for PSJ

Mercury is notoriously difficult to form in solar system simulations, due to its small mass and iron-rich composition. Smooth particle hydrodynamics simulations of collisions have found that a Mercury-like body could be formed by one or multiple giant impacts, but due to the chaotic nature of collisions it is difficult to create a scenario where such impacts will take place. Recent work has found more success forming Mercury analogues by adding additional embryos near Mercury's orbit. In this work, we aim to form Mercury by simulating the formation of the solar system in the presence of the giant planets Jupiter and Saturn. We test out the effect of an inner disk of embryos added on to the commonly-used narrow annulus of initial material. We form Mercury analogues with core-mass fractions (CMF) $> 0.4$ in $\sim 10\%$ of our simulations, and twice that number of Mercury analogues form during the formation process, but are unstable and do not last to the end of the simulations. Mercury analogues form at similar rates for both disks with and without an inner component, and most of our Mercury analogues have lower CMF than that of Mercury, $\sim 0.7$, due to significant accretion of debris material. We suggest that a more in-depth understanding of the fraction of debris mass that is lost to collisional grinding is necessary to understand Mercury's formation, or some additional mechanism is required to stop this debris from accreting., Comment: 15 pages, accepted for publication in ApJ

We explore binary asteroid formation by spin-up and rotational disruption considering the NASA DART mission's encounter with the Didymos-Dimorphos binary, which was the first small binary visited by a spacecraft. Using a suite of $N$-body simulations, we follow the gravitational accumulation of a satellite from meter-sized particles following a mass-shedding event from a rapidly rotating primary. The satellite's formation is chaotic, as it undergoes a series of collisions, mergers, and close gravitational encounters with other moonlets, leading to a wide range of outcomes in terms of the satellite's mass, shape, orbit, and rotation state. We find that a Dimorphos-like satellite can form rapidly, in a matter of days, following a realistic mass-shedding event in which only ${\sim}2-3\%$ of the primary's mass is shed. Satellites can form in synchronous rotation due to their formation near the Roche limit. There is a strong preference for forming prolate (elongated) satellites, although some simulations result in oblate spheroids like Dimorphos. The distribution of simulated secondary shapes is broadly consistent with other binary systems, measured through radar or lightcurves. Unless Dimorphos's shape is an outlier, and considering the observational bias against lightcurve-based determination of secondary elongations for oblate bodies, we suggest there could be a significant population of oblate secondaries. If these satellites initially form with elongated shapes, a yet-unidentified pathway is needed to explain how they become oblate. Finally, we show that this chaotic formation pathway occasionally forms asteroid pairs and stable triples, including co-orbital satellites and satellites in mean motion resonances., Comment: 41 pages, 23 figures, accepted to PSJ, movies available at https://zenodo.org/doi/10.5281/zenodo.8387043

The impact of the DART spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos' orbit substantially, largely from the ejection of material. We present results from twelve Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ~1.4 magnitudes, we find consistent dimming rates of 0.11-0.12 magnitudes/day in the first week, and 0.08-0.09 magnitudes/day over the entire study period. The system returned to its pre-impact brightness 24.3-25.3 days after impact through the primary ejecta tail remained. The dimming paused briefly eight days after impact, near in time to the appearance of the second tail. This was likely due to a secondary release of material after re-impact of a boulder released in the initial impact, through movement of the primary ejecta through the aperture likely played a role., Comment: 16 pages, 5 Figures, accepted in the Astrophysical Journal Letters (ApJL) on October 16, 2023

It remains an elusive goal to simultaneously model the astrophysics of Solar System accretion while reproducing the mantle chemistry of more than one inner terrestrial planet. Here, we used a multistage core-mantle differentiation model based on Rubie et al. (2011,2015) to track the formation and composition of Earth and Mars in various Grand Tack formation simulations. Prior studies showed that in order to recreate Earth's mantle composition, it must grow first from reduced (Fe-metal rich and O-poor) building blocks and then from increasingly oxidized (FeO rich) material. This accretion chemistry occurs when an oxidation gradient exists across the disk so that the innermost solids are reduced and increasingly oxidized material is found at greater heliocentric distances. For a suite of Grand Tack simulations, we investigated whether Earth and Mars can be simultaneously produced by the same oxidation gradient. Our model did not find an oxidation gradient that simultaneously reproduces the mantle composition of Earth and Mars. Due to its small mass and rapid formation, the formation history of Mars-like planets is very stochastic which decreases the likelihood of compatibility with an Earth-producing oxidation gradient in any given realization. To reconcile the accretion history and ideal chemistry of the Mars-like planet with the oxidation gradient of an Earth-producing disk, we determined where in the Earth-producing disk Mars must have formed. We find that the FeO-rich composition of the Martian mantle requires that Mars' building blocks must originate exterior to 1.0 astronomical units (AU)., Comment: 20 pages, 7 figures, Accepted for publication at ICARUS 04/2023

The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on September 26, 2022 as a planetary defense test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defense, intended to validate kinetic impact as a means of asteroid deflection. Here we report the first determination of the momentum transferred to an asteroid by kinetic impact. Based on the change in the binary orbit period, we find an instantaneous reduction in Dimorphos's along-track orbital velocity component of 2.70 +/- 0.10 mm/s, indicating enhanced momentum transfer due to recoil from ejecta streams produced by the impact. For a Dimorphos bulk density range of 1,500 to 3,300 kg/m$^3$, we find that the expected value of the momentum enhancement factor, $\beta$, ranges between 2.2 and 4.9, depending on the mass of Dimorphos. If Dimorphos and Didymos are assumed to have equal densities of 2,400 kg/m$^3$, $\beta$= 3.61 +0.19/-0.25 (1 $\sigma$). These $\beta$ values indicate that significantly more momentum was transferred to Dimorphos from the escaping impact ejecta than was incident with DART. Therefore, the DART kinetic impact was highly effective in deflecting the asteroid Dimorphos., Comment: accepted by Nature

While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority space mission related to asteroid mitigation. NASA's Double Asteroid Redirection Test (DART) mission is the first full-scale test of kinetic impact technology. The mission's target asteroid was Dimorphos, the secondary member of the S-type binary near-Earth asteroid (65803) Didymos. This binary asteroid system was chosen to enable ground-based telescopes to quantify the asteroid deflection caused by DART's impact. While past missions have utilized impactors to investigate the properties of small bodies those earlier missions were not intended to deflect their targets and did not achieve measurable deflections. Here we report the DART spacecraft's autonomous kinetic impact into Dimorphos and reconstruct the impact event, including the timeline leading to impact, the location and nature of the DART impact site, and the size and shape of Dimorphos. The successful impact of the DART spacecraft with Dimorphos and the resulting change in Dimorphos's orbit demonstrates that kinetic impactor technology is a viable technique to potentially defend Earth if necessary., Comment: Accepted by Nature

Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstrated the activation process of an asteroid from an impact under precisely known impact conditions. Here we report the observations of the DART impact ejecta with the Hubble Space Telescope (HST) from impact time T+15 minutes to T+18.5 days at spatial resolutions of ~2.1 km per pixel. Our observations reveal a complex evolution of ejecta, which is first dominated by the gravitational interaction between the Didymos binary system and the ejected dust and later by solar radiation pressure. The lowest-speed ejecta dispersed via a sustained tail that displayed a consistent morphology with previously observed asteroid tails thought to be produced by impact. The ejecta evolution following DART's controlled impact experiment thus provides a framework for understanding the fundamental mechanisms acting on asteroids disrupted by natural impact., Comment: accepted by Nature

The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, Grand Tack, and Early Instability) for which we compile a total of 433 N-body simulations. We relate the outputs of these simulations to the chemistry of Earth's core and mantle using a melt-scaling law combined with a multi-stage model of core formation. We find the compositions of Earth analogs to be largely governed by the fraction of equilibrating embryo cores and the initial embryo masses in N-body simulations. Simulation type may be important when considering magma ocean lifetimes, where Grand Tack simulations have the largest amounts of material accreted after the last giant impact. However, we cannot rule out any accretion scenarios or initial embryo masses due to the sensitivity of Earth's mantle composition to different parameters and the stochastic nature of N-body simulations. Comparing the last embryo impacts experienced by Earth analogs to specific Moon-forming scenarios, we find the characteristics of the Moon-forming impact are dependent on the initial conditions in N-body simulations where larger initial embryo masses promote larger and slower Moon-forming impactors. Mars-sized initial embryos are most consistent with the canonical hit-and-run scenario onto a solid mantle. Our results suggest that constraining the fraction of equilibrating impactor core and the initial embryo masses in N-body simulations could be significant for understanding both Earth's accretion history and characteristics of the Moon-forming impact.

Rocky planets both in and outside of our solar system are observed to have a range of core-mass fractions (CMFs). Imperfect collisions can preferentially strip mantle material from a planet, changing its CMF, and are therefore thought to be the most likely cause of this observed CMF variation. However, previous work that implements these collisions into N-body simulations of planet formation has struggled to reliably form high CMF super-Earths. In this work, we specify our initial conditions and simulation parameters to maximize the prevalence of high-energy, CMF-changing collisions in order to form planets with highly diverse CMFs. High-energy collisions have a large $v_{imp}/v_{esc}$ ratio, so we maximize this ratio by starting simulations with high-eccentricity and inclination disks to increase the difference in their orbital velocities, maximizing $v_{imp}$. Additionally, we minimize $v_{esc}$ by starting with small embryos. The final planets undergo more high-energy, debris-producing collisions, and experience significant CMF change over their formation. However, we find that a number of processes work together to average out the CMF of a planet over time, therefore we do not consistently form high-CMF, high mass planets. We do form high-CMF planets below 0.5 $M_{\oplus}$. Additionally, we find in these highly eccentric environments, loss of debris mass due to collisional grinding has a significant effect on final planet masses and CMFs, resulting in smaller planets and a higher average planet CMF. This work highlights the importance of improving measurements of high-density planets to better constrain their CMFs., Comment: 19 pages. This paper has been accepted to ApJ

The composition the atmosphere of Venus results from the integration of many processes entering into play over the entire geological history of the planet. Determining the elemental abundances and isotopic ratios of noble gases (He, Ne, Ar, Kr, Xe) and stable isotopes (H, C, N, O, S) in the Venus atmosphere is a high priority scientific target since it could open a window on the origin and early evolution of the entire planet. This chapter provides an overview of the existing dataset on noble gases and stable isotopes in the Venus atmosphere. The current state of knowledge on the origin and early and long-term evolution of the Venus atmosphere deduced from this dataset is summarized. A list of persistent and new unsolved scientific questions stemming from recent studies of planetary atmospheres (Venus, Earth and Mars) are described. Important mission requirements pertaining to the measurement of volatile elements in the atmosphere of Venus as well as potential technical difficulties are outlined., Comment: 40 pages, 10 figures, 1 table

Conventional planet formation theory suggests that chondritic materials have delivered crucial atmospheric and hydrospheric elements such as carbon (C), nitrogen (N), and hydrogen (H) onto primitive Earth. However, recent measurements highlight the significant elemental ratio discrepancies between terrestrial parent bodies and the supposed planet building blocks. Here we present a volatile evolution model during the assembly of Earth and Earth-like planets. Our model includes impact losses, atmosphere-mantle exchange, and time dependent effects of accretion and outgassing calculated from dynamical modeling outcomes. Exploring a wide range of planetesimal properties (i.e., size and composition) as well as impact history informed by N-body accretion simulations, we find that while the degree of CNH fractionation has inherent stochasticity, the evolution of C/N and C/H ratios can be traced back to certain properties of the protoplanet and projectiles. Interestingly, the majority of our Earth-like planets acquire superchondritic final C/N ratios, implying that the volatile elemental ratios on terrestrial planets are driven by the complex interplay between delivery, atmospheric ablation, and mantle degassing., Comment: 24 pages, 10 figures, 3 tables; accepted to Earth and Planetary Science Letters

Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to that expected from loss by partial vaporisation/condensation. Here we reconcile these observations by showing, using N-body simulations, that Earth accreted stochastically from many precursor bodies whose variable compositions reflect the temperatures at which they formed. Impact-induced atmospheric loss was efficient only when the proto-Earth was small, and elements that accreted thereafter retain near-chondritic isotope ratios. Earth's composition is reproduced when initial temperatures of planetesimal- to embryo-sized bodies are set by disk accretion rates of (1.08 $\pm$ 0.17) $\times$ 10$^{-7}$ solar masses/yr, although they may be perturbed by $^{26}$Al heating on bodies formed at different times. The model implies a heliocentric gradient in composition and rapid planetesimal formation within $\sim$ 1 Myr, in accord with radiometric volatile depletion ages of Earth., Comment: 13 pages, 4 figures. Nat Astron (2022)

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $\beta$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $\beta$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous., Comment: 23 pages, 13 figures, published in PSJ

The Solar System's orbital structure is thought to have been sculpted by an episode of dynamical instability among the giant planets. However, the instability trigger and timing have not been clearly established. Hydrodynamical modeling has shown that while the Sun's gaseous protoplanetary disk was present the giant planets migrated into a compact orbital configuration in a chain of resonances. Here we use dynamical simulations to show that the giant planets' instability was likely triggered by the dispersal of the gaseous disk. As the disk evaporated from the inside-out, its inner edge swept successively across and dynamically perturbed each planet's orbit in turn. The associated orbital shift caused a dynamical compression of the exterior part of the system, ultimately triggering instability. The final orbits of our simulated systems match those of the Solar System for a viable range of astrophysical parameters. The giant planet instability therefore took place as the gaseous disk dissipated, constrained by astronomical observations to be a few to ten million years after the birth of the Solar System. Terrestrial planet formation would not complete until after such an early giant planet instability; the growing terrestrial planets may even have been sculpted by its perturbations, explaining the small mass of Mars relative to Earth., Comment: Authors' version, 12 pages, 4 figures

Planets between 1-4 Earth radii with orbital periods <100 days are strikingly common. The migration model proposes that super-Earths migrate inwards and pile up at the disk inner edge in chains of mean motion resonances. After gas disk dispersal, simulations show that super-Earth's gravitational interactions can naturally break their resonant configuration leading to a late phase of giant impacts. The instability phase is key to matching the orbital spacing of observed systems. Yet, most previous simulations have modelled collisions as perfect accretion events, ignoring fragmentation. In this work, we investigate the impact of imperfect accretion on the breaking the chains scenario. We performed N-body simulations starting from distributions of planetary embryos and modelling the effects of pebble accretion and migration in the gas disk. Our simulations also follow the long-term dynamical evolution of super-Earths after the gas disk dissipation. We compared the results of simulations where collisions are treated as perfect merging events with those where imperfect accretion and fragmentation are allowed. We concluded that the perfect accretion is a suitable approximation in this regime, from a dynamical point of view. Although fragmentation events are common, only ~10% of the system mass is fragmented during a typical "late instability phase", with fragments being mostly reacreted by surviving planets. This limited total mass in fragments proved to be insufficient to alter qualitatively the final system dynamical configuration -- e.g. promote strong dynamical friction or residual migration -- compared to simulations where fragmentation is neglected., Comment: 13 pages, 8 figures and 1 table. Accepted for publication in MNRAS

Most close-in planetary satellites are in synchronous rotation, which is usually the stable end-point of tidal despinning. Saturn's moon Hyperion is a notable exception by having a chaotic rotation. Hyperion's dynamical state is a consequence of its high eccentricity and its highly prolate shape (Wisdom et al. 1984). As many binary asteroids also have elongated secondaries, chaotic rotation is expected for moons in eccentric binaries (\'Cuk and Nesvorn\'y, 2010), and a minority of asteroidal secondaries may be in that state (Pravec et al. 2016). The question of secondary rotation is also important for the action of the BYORP effect, which can quickly evolve orbits of synchronous (but not non-synchronous) secondaries (\'Cuk and Burns, 2005). Here we report results of a large set of short numerical simulations which indicate that, apart from synchronous and classic chaotic rotation, close-in irregularly-shaped asteroidal secondaries can occupy an additional, intermediate rotational state. In this ``barrel instability'' the secondary slowly rolls along its long axis, while the longest axis is staying largely aligned with the primary-secondary line. This behavior may be more difficult to detect through lightcurves than a fully chaotic rotation, but would likewise shut down BYORP. We show that the binary's eccentricity, separation measured in secondary's radii and the secondary's shape are all important for determining whether the system settles in synchronous rotation, chaotic tumbling, or barrel instability. We compare our results for synthetic asteroids with known binary pairs to determine which of these behaviors may be present in the Near-Earth Asteroid binary population., Comment: Accepted for PSJ

The NASA Double Asteroid Redirection Test (DART) mission is a planetary defense-driven test of a kinetic impactor on Dimorphos, the satellite of the binary asteroid 65803 Didymos. DART will intercept Dimorphos at a relative speed of ${\sim}6.5 \text{ km s}^{-1}$, perturbing Dimorphos's orbital velocity and changing the binary orbital period. We present three independent methods (one analytic and two numerical) to investigate the post-impact attitude stability of Dimorphos as a function of its axial ratios, $a/b$ and $b/c$ ($a \ge b \ge c$), and the momentum transfer efficiency $\beta$. The first method uses a novel analytic approach in which we assume a circular orbit and a point-mass primary that identifies four fundamental frequencies of motion corresponding to the secondary's mean motion, libration, precession, and nutation frequencies. At resonance locations among these four frequencies, we find that attitude instabilities are possible. Using two independent numerical codes, we recover many of the resonances predicted by the analytic model and indeed show attitude instability. With one code, we use fast Lyapunov indicators to show that the secondary's attitude can evolve chaotically near the resonance locations. Then, using a high-fidelity numerical model, we find that Dimorphos enters a chaotic tumbling state near the resonance locations and is especially prone to unstable rotation about its long axis, which can be confirmed by ESA's Hera mission arriving at Didymos in late 2026. We also show that a fully coupled treatment of the spin and orbital evolution of both bodies is crucial to accurately model the long-term evolution of the secondary's spin state and libration amplitude. Finally, we discuss the implications of a post-impact tumbling or rolling state, including the possibility of terminating BYORP evolution if Dimorphos is no longer in synchronous rotation., Comment: 38 pages, 13 figures, Accepted for publication in Icarus

Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core-mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (``perfect merging'') and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of non-accreted projectiles in hit-and-run collisions. We find that the results of the neural-network model feature a wider range of final core mass fractions and metal-silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than 0.1 Earth's masses. For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies., Comment: 21 pages, 2 tables, 7 figures. Published open access on PSJ: https://iopscience.iop.org/article/10.3847/PSJ/abf0ad

Although numerous white dwarf stars host dusty debris disks, the temperature distribution of these stars differs significantly from the white dwarf population as a whole. Dusty debris disks exist exclusively around white dwarfs cooler than 27,000 K. This is all the more enigmatic given that the formation processes of dusty debris disks should favor younger, hotter white dwarfs, which likely host more dynamically unstable planetary systems. Here we apply a sophisticated material sublimation model to white dwarf systems to show that these statistics are actually a natural result of the interplay of thermal and tidal forces, and show how they define the circumstellar regions where dusty debris disks can form. We demonstrate that these processes tend to prevent stability against both sublimative destruction and reaccretion into planetesimals for rocky materials until white dwarfs cool to below ~25,000-32,000 K, in agreement with the observed limit of ~27,000 K. For pure water ice, this critical temperature is less than 2,700 K (requiring a cooling age older the universe); this precludes pure water ice-rich debris disks forming through the accepted two-step mechanism. The critical temperature is size-dependent; more massive white dwarfs could potentially host dusty debris disks at warmer temperatures. Our model suggests that the location of the disks within the PG 0010+280, GD 56, GD 362, and PG 1541+651 systems are consistent with a forsterite-dominated olivine composition. We also find that very cool white dwarfs may simultaneously host multiple, independently formed dusty debris disks, consistent with observations of the LSPM J0207+3331 system., Comment: Astrophysical Journal Letters. Accepted for Publication

Dynamical scenarios of terrestrial planets formation involve strong perturbations of the inner part of the solar system by the giant-planets, leading to enhanced impact velocities and subsequent collisional erosion. We quantitatively estimate the effect of collisional erosion on the resulting composition of Earth, and estimate how it may provide information on the dynamical context of its formation. We simulate and quantify the erosion of Earth's crust in the context of Solar System formation scenarios, including the classical model and Grand Tack scenario that invokes orbital migration of Jupiter during the gaseous disk phase (Walsh et al., 2011; Raymond et al., 2018). We find that collisional erosion of the early crust is unlikely to produce an excess of about 6% of the Sm/Nd ratio in terrestrial rock samples compared to chondrites for most simulations. Only Grand Tack simulations in which the last giant impact on Earth occurred later than 50 million years after the start of Solar System formation can account for such an offset. However, this time frame is consistent with current cosmochemical and dynamical estimates of the Moon forming impact (Chyba, 1991; Walker, 2009; Touboul et al.,2007, 2009, 2015; Pepin and Porcelli, 2006; Norman et al., 2003; Nyquist et al., 2006; Boyet et al.,2015). Such a late fractionation in the Sm/Nd ratio is unlikely to be responsible for a 20-ppm $^{142}$Nd excess in terrestrial rocks due to the half life of the radiogenic system. Additionally, such a large and late fractionation in the Sm/Nd ratio would accordingly induce non-observed anomalies in the $^{143}$Nd/$^{144}$Nd ratio. Considering our results, the Grand Tack model with a late Moon-forming impact cannot be easily reconciled with the Nd isotopic Earth contents., Comment: 22 pages, 3 figures, supplementary, to be published

Small bodies are time-capsules of different eras of solar system history from the most primitive materials within the solar system to evolved pieces of larger bodies. A small body sample return program is an essential component of small body exploration, and such a program should include opportunities for both missions and laboratory analysis., Comment: A White Paper submitted to the Planetary Science and Astrobiology Decadal Survey 2023-2032

Here, we develop scaling laws for (1) the distribution of impact-induced heat within the mantle and (2) shape of the impact-induced melt 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 $\theta$ ($0^{\circ}, 30^{\circ}, 60^{\circ}$, and $90^{\circ}$), total mass $M_T$ ($1M_{\rm Mars}-53M_{\rm Mars}$, where $M_{\rm Mars}$ is the mass of Mars), impact velocity $v_{\rm imp}$ ($v_{\rm esc} - 2v_{\rm esc}$, where $v_{\rm esc}$ is the mutual escape velocity), and impactor-to-total mass ratio $\gamma$ ($0.03-0.5$). We find that the equilibrium pressure at the base of a melt pool can be higher (up to $\approx 80 \%$) than those of radially-uniform global magma ocean models. This could have a significant impact on element partitioning. These melt scaling laws are publicly available on GitHub ($\href{https://github.com/mikinakajima/MeltScalingLaw}{https://github.com/mikinakajima/MeltScalingLaw}$)., Comment: Published in Earth and Planetary Science Letters

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., Comment: 15 pages, 14 figures. Accepted for publication in MNRAS

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., Comment: 29 pages, 32 figures

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., Comment: 34 pages, 5 figures, published in EPSL

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., Comment: 34 pages, 9 figures, to be published in Icarus

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., Comment: Astro2020 APC White Paper: State of the Profession Consideration

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)., Comment: 25 pages, accepted by A&A, see companion papers by Izidoro et al. (2019) and Lambrechts et al. (2019)

At least 30\% of main sequence stars host planets with sizes of 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 inner edge of the disk 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 original resonant configurations. Our simulations naturally match observations when they produce 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 (abridged)., Comment: Accepted in A&A, version including language editing

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., Comment: Accepted in A&A

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., Comment: Accepted Jan 17, 2018 for publication in Space Science Reviews under the topical collection The Delivery of Water to Protoplanets, Planets and Satellites and in an ISSI Space Science Series book of the same title

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., Comment: Supplementary material included

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., Comment: 10 text pages, 3 figures, 3 SI pages

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., Comment: 12 pages, 6 figures. Accepted to MNRAS. Blog post about the paper at http://planetplanet.net/2016/02/21/did-the-solar-system-form-inside-0ut/

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., Comment: 11 pages, 3 figures, submitted to IAU Proceedings

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)., Comment: 56 pages, 6 figures, 1 table. Accepted for publication in Icarus

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., Comment: 42 pages, 4 figures, contribution to the Asteroids 4 book