Colors and collision rates within the Kuiper belt: problems with the collisional resurfacing scenario

We present a numerical check of the collisional resurfacing (CR) hypothesis proposed to explain the observed color diversity within the Kuiper Belt (where surface reddening due to space weathering is counteracted by regular resurfacing of neutral material after mutual collisions). Deterministic simulations are performed in order to estimate the relative spatial distribution of kinetic energy received by collisions, , for a population of target Kuiper Belt objects (KBOs) embedded in a swarm of impactors distributed within the belt. Four different impactor disks have been considered, depending on the excitation and the external limit of the belt and the density of the scattered KBOs (SKBOs) population. The obtained results are compared to the relative color index distribution within the observed Kuiper Belt, in order to derive possible similarities between the high vs low objects spatial distribution in our simulations and the bluer vs redder KBOs distribution in the “real” Kuiper Belt. Such similarities are found for several important features, in particular the general correlations between highly impacted objects and high rms excitation and low perihelion q values that are in good agreement with equivalent correlations found for the bluest objects of the observed belt. Nevertheless, simulations disagree with observations on two crucial points. (1) The plutinos are significantly more collisionally affected than the rest of our test KBO population, whereas there is no observed tendency toward bluer plutinos. (2) There is always a much stronger correlation between and eccentricities than inclinations, whereas observations show just the opposite feature. The presence of numerous SKBO impactors could significantly damp these problematic features, but cannot erase them. Whether these contradictions invalidate the whole CR scenario or not remains yet uncertain, since the physical processes at play are still far from being fully understood and the sample of available observational data is still relatively limited. But it seems nevertheless that the scenario might not hold in its simple present form.     

Origin of the structure of the Kuiper belt during a dynamical instability in the orbits of Uranus and Neptune

We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. It successfully explains, for the first time, the current orbital architecture of the giant planets [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461], the existence of the Trojans populations of Jupiter and Neptune [Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465], and the origin of the late heavy bombardment of the terrestrial planets [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469]. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ∼35 AU to the Kuiper belt without excessive inclination excitation. In particular, we show that during the phase when Neptune's eccentricity is large, the region interior to its 1:2 mean motion resonance becomes unstable and disk particles can diffuse into this area. In addition, we perform numerical simulations where the planets are forced to evolve using fictitious analytic forces, in a way consistent with the direct N-body simulations of the Nice model. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ∼27 AU and an eccentricity of ∼0.3, and that subsequently Neptune's eccentricity damped in ∼1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level. In particular, our results explain, at least qualitatively: (1) the co-existence of resonant and non-resonant populations, (2) the eccentricity-inclination distribution of the Plutinos, (3) the peculiar semi-major axis—eccentricity distribution in the classical belt, (4) the outer edge at the 1:2 mean motion resonance with Neptune, (5) the bi-modal inclination distribution of the classical population, (6) the correlations between inclination and physical properties in the classical Kuiper belt, and (7) the existence of the so-called extended scattered disk. Nevertheless, we observe in the simulations a deficit of nearly-circular objects in the classical Kuiper belt.     

The edge of the Kuiper belt: the Planet X scenario

Our goal is to determine whether or not the observed sudden termination of the Edgeworth-Kuiper belt can be the result of perturbations from a hypothetical planet. We investigate the effects that such an object would produce on the primordial orbital distribution if the trans-neptunian objects, for a range of masses and orbital parameters of the hypothetical planet. In this numerical investigation, the motion of the hypothetical planet was influenced by the existing planets but not by its interaction with the disk. We find that no set of parameters produce results that match the observed data. Dynamical interaction with the disk is likely to be important so that the orbit of the hypothetical planet changes significantly during the integration interval. This is also discussed. The overall conclusion is that none of the models for the hypothetical planet that were investigated can reproduce the observed features of the Edgeworth-Kuiper belt starting from any probable primordial distribution.     

Origin of scattered disk resonant TNOs: Evidence for an ancient excited Kuiper belt of 50 AU radius

Resonance occupation of trans-neptunian objects (TNOs) in the scattered disk (>48 AU) was investigated by integrating the orbits of 85 observed members for 4 Gyr. Twenty seven TNOs were locked in the 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4 resonances. We then explored mechanisms for the origin of the resonant structure in the scattered disk, in particular the long-term 9:4, 5:2, and 8:3 resonant TNOs (median 4 Gyr), by performing large scale simulations involving Neptune scattering and planetary migration over an initially excited planetesimals disk (wide range of eccentricities and inclinations). To explain the formation of Gyr-resident populations in such distant resonances, our results suggest the existence of a primordial planetesimal disk of at least 45-50 AU radius that suffered a dynamical perturbation leading to 0.1-0.3 or greater eccentricities and a range of inclinations up to ∼20° during early stages of the Solar System history, before planetary migration.     

Considerations on the magnitude distributions of the Kuiper belt and of the Jupiter Trojans

By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the ‘Nice model’. This leads us to predict that the HD of the hot population should have the same slope of the HD of the cold population for 6.5<H<9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H∼9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. [Bernstein, G.M., and 5 colleagues, 2004. Astron. J. 128, 1364-1390] and Fuentes and Holman [Fuentes, C.I., Holman, M.J., 2008. Astron. J. 136, 83-97]. Finally, we show that the constraint on the total mass of the primordial trans-neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU.     

Evolution of the obliquities of the giant planets in encounters during migration

Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] have proposed that the current orbital architecture of the outer Solar System could have been established if it was initially compact and Jupiter and Saturn crossed the 2:1 orbital resonance by divergent migration. The crossing led to close encounters among the giant planets, but the orbital eccentricities and inclinations were damped to their current values by interactions with planetesimals. Brunini [Brunini, A., 2006. Nature 440, 1163-1165] has presented widely publicized numerical results showing that the close encounters led to the current obliquities of the giant planets. We present a simple analytic argument which shows that the change in the spin direction of a planet relative to an inertial frame during an encounter between the planets is very small and that the change in the obliquity (which is measured from the orbit normal) is due to the change in the orbital inclination. Since the inclinations are damped by planetesimal interactions on timescales much shorter than the timescales on which the spins precess due to the torques from the Sun, especially for Uranus and Neptune, the obliquities should return to small values if they are small before the encounters. We have performed simulations using the symplectic integrator SyMBA, modified to include spin evolution due to the torques from the Sun and mutual planetary interactions. Our numerical results are consistent with the analytic argument for no significant remnant obliquities.     

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

A hybrid scenario for gas giant planet formation in rings

The core-accretion mechanism for gas giant formation may be too slow to create all observed gas giant planets during reasonable gas disk lifetimes, but it has yet to be firmly established that the disk instability model can produce permanent bound gaseous protoplanets under realistic conditions. Based on our recent simulations of gravitational instabilities in disks around young stars, we suggest that, even if instabilities due to disk self-gravity do not produce gaseous protoplanets directly, they may create persistent dense rings that are conducive to accelerated growth of gas giants through core accretion. The rings occur at and near the boundary between stable and unstable regions of the disk and appear to be produced by resonances with discrete spiral modes on the unstable side.     

On the collisional disruption of porous icy targets simulating Kuiper belt objects

We present results from 27 impact experiments using porous (porosity ranging from 0.39 to 0.54) ice targets and solid ice projectiles at impact speeds ranging from 90 to 155 m/s. These targets were designed to simulate Kuiper Belt Objects (KBOs) in structure. We measured a specific energy for shattering, , of 2.1×¹⁰⁵ erg/g for those snowball targets hit by intact ice projectiles; this is of the same order as that measured for solid ice targets. The fragment mass distribution follows a power law, although the exponent is not simply related to the largest fragment size as assumed by fragmentation models. We provide the first measurement of the three-dimensional mass-velocity distribution for disrupted ice targets and find that fragment speeds range from ∼2 to ∼20 m/s. The fraction of collisional kinetic energy that is partitioned into ejecta speeds is between 1 and 15% (although it should be noted that the lower limit is more reliable than the upper).     

The formation of the Oort cloud in open cluster environments

We study the influence of an open cluster environment on the formation and current structure of the Oort cloud. To do this, we have run 19 different simulations of the formation of the Oort cloud for 4.5 Gyrs. In each simulation, the Solar System spends its first 100 Myrs in a different open cluster environment before transitioning to its current field environment. We find that, compared to forming in the field environment, the inner Oort cloud is preferentially loaded with comets while the Sun resides in the open cluster and that most of this material remains locked in the interior of the cloud for the next 4.4 Gyrs. In addition, the outer Oort cloud trapping efficiencies we observe in our simulations are lower than previous formation models by about a factor of 2, possibly implying an even more massive early planetesimal disk. Furthermore, some of our simulations reproduce the orbits of observed extended scattered disk objects, which may serve as an observational constraint on the Sun's early environment. Depending on the particular open cluster environment, the properties of the inner Oort cloud and extended scattered disk can vary widely. On the other hand, the outer portions of the Oort cloud in each of our simulations are all similar.     

The plane of the Edgeworth-Kuiper belt

We examine possible locations for the primordial disk of the Edgeworth-Kuiper Belt (EKB), using several subsets of the known objects as markers of the total mass distribution. Using a secular perturbation theory, we find that the primordial plane of the EKB could have remained thin enough to escape detection only if it is clustered very closely about the invariable plane of the Solar System.     

Towards initial mass functions for asteroids and Kuiper Belt Objects

Our goal is to understand primary accretion of the first planetesimals. Some examples are seen today in the asteroid belt, providing the parent bodies for the primitive meteorites. The primitive meteorite record suggests that sizeable planetesimals formed over a period longer than a million years, each of which being composed entirely of an unusual, but homogeneous, mixture of millimeter-size particles. We sketch a scenario that might help explain how this occurred, in which primary accretion of 10-100 km size planetesimals proceeds directly, if sporadically, from aerodynamically-sorted millimeter-size particles (generically “chondrules”). These planetesimal sizes are in general agreement with the currently observed asteroid mass peak near 100 km diameter, which has been identified as a “fossil” property of the pre-erosion, pre-depletion population. We extend our primary accretion theory to make predictions for outer Solar System planetesimals, which may also have a preferred size in the 100 km diameter range. We estimate formation rates of planetesimals and explore parameter space to assess the conditions needed to match estimates of both asteroid and Kuiper Belt Object (KBO) formation rates. For parameters that satisfy observed mass accretion rates of Myr-old protoplanetary nebulae, the scenario is roughly consistent with not only the “fossil” sizes of the asteroids, and their estimated production rates, but also with the observed spread in formation ages of chondrules in a given chondrite, and with a tolerably small radial diffusive mixing during this time between formation and accretion. As previously noted, the model naturally helps explain the peculiar size distribution of chondrules within such objects. The optimum range of parameters, however, represents a higher gas density and fractional abundance of solids, and a smaller difference between Keplerian and pressure-supported orbital velocities, than “canonical” models of the solar nebula. We discuss several potential explanations for these differences. The scenario also produces 10-100 km diameter primary KBOs, and also requires an enhanced abundance of solids to match the mass production rate estimates for KBOs (and presumably the planetesimal precursors of the ice giants themselves). We discuss the advantages and plausibility of the scenario, outstanding issues, and future directions of research.     

The evidence of an early stellar encounter in Edgeworth-Kuiper belt

We investigate the influence of a stellar fly-by encounter on the Edgeworth-Kuiper belt objects through numerical orbital calculations, in order to explain both mass depletion and high orbital inclinations of the classical Edgeworth-Kuiper belt (CEKB) objects, which have semimajor axis of 42-48 AU and perihelia beyond 35 AU. The observationally inferred total mass of the CEKB is ∼1/10 Earth masses, which is only ∼0.02 of that extrapolated from the minimum-mass solar nebula model. The CEKB consists of bimodal population: “hot population” with inclinations i?0.2-0.6 radians and “cold population” with i?0.1. The observationally suggested difference in size and color of objects between the two populations may imply different origins of the two populations. We find that both the depletion of solid materials in the CEKB and the formation of the hot population are accounted for by a single close stellar encounter with pericenter distance of 80-100 AU and inclination relative to the initial protoplanetary disk ?50°-70°. Such a stellar encounter highly pumps up eccentricities of most objects in the CEKB and then their perihelia migrate within 35 AU. These objects would be removed by Neptune's perturbations after Neptune is formed at or migrates to the current position (30 AU). Less than 10% of the original objects remain in stable orbits with small eccentricities and perihelion distances larger than 35 AU, in the CEKB, which is consistent with the observation. We find that i of the remaining objects are as large as that of the observed hot population. The only problem is how to stop Neptune's migration at ∼30 AU, which is addressed in a separate paper. The depletion by the stellar encounter extends deeply into ∼30-35 AU, which provides the basis of the formation model for the cold population through Neptune's outward migration by Levison and Morbidelli (2003, Nature, 426, 419-421). The combination of our model with Levison and Morbidelli's model could consistently explain the mass depletion, truncation at 50 AU, bimodal distribution in i, and differences in size and color between the hot and the cold populations in the CEKB.     

The origin of the Kuiper Belt high-inclination population

I simulate the orbital evolution of the four major planets and a massive primordial planetesimal disk composed of 10⁴ objects, which perturb the planets but not themselves. As Neptune migrates by energy and angular momentum exchange with the planetesimals, a large number of primordial Neptune-scattered objects are formed. These objects may experience secular, Kozai, and mean motion resonances that induce temporary decrease of their eccentricities. Because planets are migrating, some planetesimals can escape those resonances while in a low-eccentricity incursion, thus avoiding the return path to Neptune close encounter dynamics. In the end, this mechanism produces stable orbits with high inclination and moderate eccentricities. The population so formed together with the objects coming from the classical resonance sweeping process, originates a bimodal distribution for the Kuiper Belt orbits. The inclinations obtained by the simulations can attain values above 30° and their distribution resembles a debiased distribution for the high-inclination population coming from the real classical Kuiper Belt.     

Oligarchic growth of giant planets

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as “oligarchic growth.” Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk (∼10 × minimum-mass) is required to produce giant planet core-sized bodies (∼10 M_⊕) within the lifetime of the nebular gas (?10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of noninteracting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System’s two outer giant planets at their current locations on any timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.     

Planetesimal formation by turbulent concentration

The formation of 1-1000 km diameter planetesimals from dust grains in a protoplanetary disk is a key step in planet formation. Conventional models for planetesimal formation involve pairwise sticking of dust grains, or the sedimentation of dust grains to a thin layer at the disk midplane followed by gravitational instability. Each of these mechanisms is likely to be frustrated if the disk is turbulent. Particles with stopping times comparable to the turnover time of the smallest eddies in a turbulent disk can become concentrated into dense clumps that may be the precursors of planetesimals. Such particles are roughly millimeter-sized for a typical protoplanetary disk. To survive to become planetesimals, clumps need to form in regions of low vorticity to avoid rotational breakup. In addition, clumps must have sufficient self gravity to avoid break up due to the ram pressure of the surrounding gas. Given these constraints, the rate of planetesimal formation can be estimated using a cascade model for the distribution of particle concentration and vorticity within eddies of various sizes in a turbulent disk. We estimate planetesimal formation rates and planetesimal diameters as a function of distance from a star for a range of protoplanetary disk parameters. For material with a solar composition, the dust-to-gas ratio is too low to allow efficient planetesimal formation, and most solid material will remain in small particles. Enhancement of the dust-to-gas ratio by 1-2 orders of magnitude, either vertically or radially, allows most solid material to be converted into planetesimals within the typical lifetime of a disk. Such dust-to-gas ratios may occur near the disk midplane as a result of vertical settling of short-lived clumps prior to clump breakup. Planetesimal formation rates are sensitive to the assumed size and rotational speed of the largest eddies in the disk, and formation rates increase substantially if the largest eddies rotate more slowly than the disk itself. Planetesimal formation becomes more efficient with increasing distance from the star unless the disk surface density profile has a slope of −1.5 or steeper as a function of distance. Planetesimal formation rates typically increase by an order-of-magnitude or more moving outward across the snow line for a solid surface density increase of a factor of 2. In all cases considered, the modal planetesimal size increases with roughly the square root of distance from the star. Typical modal diameters are 100 km and 400 km in the regions corresponding to the asteroid belt and Kuiper belt in the Solar System, respectively.     

Planetary growth with collisional fragmentation and gas drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 10⁷ years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.     

The color of the Kuiper belt Core

Recent dynamical analyses of the Kuiper belt have introduced a rigorous classification scheme, determined the mean orbital plane, and identified “Core” and “Halo” populations as a function of inclination with respect to this plane (Elliot, J.L., Kern, S.D., Clancy, K.B., Gulbis, A.A.S., Millis, R.L., Buie, M.W., Wasserman, L.H., Chiang, E.I., Jordan, A.B., Trilling, D.E., Meech, K.J., 2005. Astron. J. 129, 1117-1162). Here, we use new observations and existing data to investigate the colors of Kuiper belt objects (KBOs) within this framework. With respect to the bulk KBO color distribution (all objects for which we have B-V and V-R colors; median B-R=1.56), we find that the population of objects classified following (Elliot, J.L., Kern, S.D., Clancy, K.B., Gulbis, A.A.S., Millis, R.L., Buie, M.W., Wasserman, L.H., Chiang, E.I., Jordan, A.B., Trilling, D.E., Meech, K.J., 2005. Astron. J. 129, 1117-1162) as Classical tends to be red (B-R>1.56) while the Scattered Near population is mostly neutral (B-R<1.56). Colors of Scattered Extended and Resonant objects are consistent with the bulk distribution. Separating objects into specific resonances demonstrates that the color of the Resonant sample is dominated by KBOs in the 3:2 resonance, which is consistent with previous findings. Unlike the objects in the 3:2 resonance, however, the majority of objects in the 5:2 resonance are neutral and all but one of the objects in the 4:3, 5:3, 7:4, 2:1, and 7:3 resonances are red. In particular, the objects in the 7:4 resonance are remarkably red. We find that the colors of KBOs in the Core (low-inclination) and Halo (high-inclination) are statistically different, with Core objects being primarily red and Halo objects having a slight tendency to be neutral. Notably, virtually all of the non-Resonant Core objects are red. This combination of low inclination, unperturbed orbits and red colors in the Core may be indicative of a relic grouping of objects.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

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

Impact angle influence in high velocity dust collisions during planetesimal formation

We have examined the influence of impact angle in collisions between small dust aggregates and larger dust targets through laboratory experiments. Targets consisted of μm-sized quartz dust and had a porosity of about 67%; the projectiles, between 1 and 5 mm in diameter, were slightly more compact (64% porosity). The collision velocity was centered at 20 m/s and impact angles range from 0° to 45°. At a given impact angle, the target gained mass for projectiles smaller than a threshold size, which decreases with increasing angle from about 3 mm to 1 mm. The fact that growth is possible up to the largest angles studied supports the idea of planetesimal formation by sweep-up of small dust aggregates.