Seth A. Jacobson

Accretion and differentiation of the terrestrial planets with implications for the compositions of early-formed Solar System bodies and accretion of water

Abstract In order to test accretion simulations as well as planetary differentiation scenarios, we have integrated a multistage core–mantle differentiation model with N-body accretion simulations. Impacts between embryos and planetesimals are considered to 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 and requires that the bulk compositions of all starting embryos and planetesimals are defined as a function of their heliocentric distances of origin. To do this, we assume that non-volatile elements are present in Solar System (CI) relative abundances in all bodies and that oxygen and H2O contents are the main compositional variables. The primary constraint on the combined model is the composition of the Earth’s primitive mantle. In addition, we aim to reproduce 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 equilibration pressure and 1–4 parameters that define the starting compositions of primitive bodies. This integrated model has been applied to six Grand Tack N-body accretion simulations. Investigations of a broad parameter space indicate that: (1) accretion of Earth was heterogeneous, (2) metal–silicate equilibration pressures increase as accretion progresses and are, on average, 60–70% of core–mantle boundary pressures at the time of each impact, and (3) a large fraction (70–100%) of the metal of impactor cores equilibrates with a small fraction of the silicate mantles of proto-planets during each core formation event. Results are highly sensitive to the compositional model for the primitive starting bodies and several accretion/core-formation models can thus be excluded. Acceptable fits to the Earth’s mantle composition are obtained only when bodies that originated close to the Sun, at

Highly siderophile elements in Earth’s mantle as a clock for the Moon-forming impact

According to the generally accepted scenario, the last giant impact on Earth formed the Moon and initiated the final phase of core formation by melting 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 Myr) 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 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 after 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 Myr after condensation.

Formation of asteroid pairs by rotational fission

Pairs of asteroids sharing similar heliocentric orbits, but not bound together, 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 process 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 the 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 primarys 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.

Oxygen isotopic evidence for vigorous mixing during the Moon-forming giant impact

Rehomogenizing the Earth-Moon system A giant impact formed the Moon, and lunar rocks provide insight into that process. Young et al. found that rocks on Earth and the Moon have identical oxygen isotopes. This suggests that well-mixed material from the giant impact must have formed both the Moon and Earths mantle. The finding also constrains the composition of the “late veneer”: material sprinkled onto Earth after the Moon-forming impact. Science, this issue p. 493 Identical oxygen isotope composition requires homogenization after the Moon-forming giant impact. Earth and the Moon are shown here to have indistinguishable oxygen isotope ratios, with a difference in Δ′17O of −1 ± 5 parts per million (2 standard error). On the basis of these data and our new planet formation simulations that include a realistic model for primordial 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 Δ′17O within approximately 1 per mil of the terrestrial value, limiting possible sources for this late addition of mass to the Earth-Moon system.

The great dichotomy of the Solar System: Small terrestrial embryos and massive giant planet cores

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) 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 similar to 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

Lunar and terrestrial planet formation in the Grand Tack scenario

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 time scale 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 approximately 60 and approximately 130 Myr after the formation of the first solids and was caused most likely by an object with a mass similar to that of Mars.

Fossilized condensation lines in the Solar System protoplanetary disk

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplan-etary disk the temperature should have decreased below water-condensation level well before the disk was photo-evaporated. Thus, the global water depletion of the inner Solar System is puzzling. 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. Thus the gas in the vicinity of the snowline always comes from farther out, where it should have already condensed, and therefore it should be dry. 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 planet with a mass larger than 20 Earth mass is present, the radial drift of particles is interrupted, because such a planet gives the disk a super-Keplerian rotation just outside of its own orbit. From this result, we propose that the precursor of Jupiter achieved this threshold mass when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. In fact, even if it cooled later, the disk inside of Jupiters 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.7 AU, 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.7 AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

Highly siderophile elements were stripped from Earth’s mantle by iron sulfide segregation

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.

The binary near-Earth Asteroid (175706) 1996 FG3 - An observational constraint on its orbital evolution

Abstract Using our photometric observations taken between April 1996 and January 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. With our 17-year long dataset, 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 ° and - 83 ° , respectively, with the mean radius of the uncertainty area of 4 ° , and the orbital period is 16.1508 ± 0.0002 h (all quoted uncertainties correspond to 3σ). 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 tigidity and quality factor of μ Q = 1.3 × 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.

Superbend project at the Advanced Light Source

To satisfy the demand for more high energy, high brightness X-ray sources at the Advanced Light Source (ALS), a plan is in place to replace three 1.3 Tesla normal conducting bending magnets with three 5 Tesla superconducting magnets (Superbends) in the year 2001. This will result in 12 new X-ray beam lines (four from each superbend) for users. The Superbend sources will be an order of magnitude higher in X-ray brightness and flux at 12 keV than the conventional 1.3 Tesla bending magnets. The Superbend project is a major upgrade to the ALS where the 3 superconducting magnets will be an integral part of the machine lattice. In this paper we discuss the current status of the Superbend projects well as precomissioning studies prior to the 2001 installation for users.