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.     

Sintering of highly porous silica-particle samples: analogues of early Solar-System aggregates

I describe a new method to make particle layers which consist of SiO₂ spheres with 0.78 μm radius. The layers were produced by sedimentation of aggregates which had grown in ballistic particle collisions, and the layers had a porosity of 0.95. They were used for experiments on sintering, i.e., the samples were heated in an oven at varying temperatures and heating durations, and the samples were analyzed by scanning electron microscopy. Based on the change of particle diameter, surface diffusion sintering and viscous flow are identified as important transformation mechanisms. The first effect dominated at the start of restructuring and the latter at higher temperatures. The neck growth of adjacent particles was fitted to a surface diffusion sintering model and predicts neck radii as a heating temperature and duration function. Between the temperature range of neck formation and of melting, further restructuring occurred which lead to dissolution of particulate structure and to densification and which resulted in a porous object consisting of straight elongated substructures which connected kinks of higher material density. The thermal transformation is important for the change of strength, collisional behavior, light-scattering properties, and thermal conductivity with relevance to dust aggregates, planetesimals, comets, interplanetary dust particles, and regolith-covered celestial bodies.     

Accretion rates of planetesimals by protoplanets embedded in nebular gas

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals’ eccentricity, inclination, planet’s radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals’ kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with 1-10 Earth masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.     

Radiative cooling within illuminated layers of dust on (pre)-planetary surfaces and its effect on dust ejection

Temperature gradients in dust beds embedded in a low pressure gaseous environment induce a lift of particles under certain conditions. This effect can erode planetesimals and enables entrainment of dust into the martian atmosphere. Here, we consider a numerical model to calculate the temperature profile in a dust bed which is subject to illumination. We consider the situation when the illumination is switched on and heats the dust bed’s surface and when it is switched off again after a certain time. The calculations focus on the heat transfer by infrared radiation within the dust layer. We find that radiative transfer within the dust bed modifies the absolute temperatures and temperature gradients significantly. This is important for effects which are sensitive to absolute temperatures, i.e. ice sublimation or melting of solids. For low thermal conductivity dust beds of 0.001 W m⁻¹ K⁻¹ it determines the temperature structure of the dust. For higher thermal conductivities the modifications are moderate with respect to dust eruptions as the order of magnitude of temperature gradients stays the same.     

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.     

Accretion among preplanetary bodies: The many faces of runaway growth

When preplanetary bodies reach proportions of ∼1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. The numerical model is novel as it attempts to preserve the individual particle nature of the bodies (like N-body codes); yet it is statistical in nature since it must incorporate the very large number of planetesimals. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R₀. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We associate the initial runaway growth phase with increasing GF-factors for the largest body. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high-mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two-component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. We identify the peak of GF with the transition between the runaway growth and oligarchy accretion stages. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius, especially at the outer disk. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion, which could lead to a very fast growth.     

Embedded star clusters and the formation of the Oort cloud: II. The effect of the primordial solar nebula

This paper deals with Oort cloud formation while the Sun was in an embedded cluster and surrounded by its primordial nebula. This work is a continuation of Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82], building on the model presented therein, and adding the aerodynamic drag and gravitational potential of the primordial solar nebula. Results are presented of numerical simulations of comets subject to the gravitational influence of the Sun, Jupiter, Saturn, star cluster and primordial solar nebula; some of the simulations included the gravitational influence of Uranus and Neptune as well. The primordial solar nebula was approximated by the minimum-mass Hayashi model [Hayashi, C., Nakozawa, K., Nakagawa, Y., 1985. In: Black, D.C., Matthews, M.S. (Eds.). Protostars and Planets II. Univ. of Arizona Press, Tucson, AZ] whose inner and outer radii have been truncated at various distances from the Sun. A comet size of 1.7 km was used for most of our simulations. In all of our simulations, the density of the primordial solar nebula decayed exponentially with an e-folding time of 2 Myr. It turns out that when the primordial solar nebula extends much beyond Saturn or Neptune, virtually no material will end up in the Oort cloud (OC) during this phase. Instead, the majority of the material will be on circular orbits inside of Jupiter if the inner edge of the disk is well inside Jupiter's orbit. If the disk's inner edge is beyond Jupiter's orbit, most comets end up on orbits in exterior mean-motion resonances with Saturn when Uranus and Neptune are not present. In those cases where the outer edge of the disk is close to Saturn or Neptune, the fraction of material that ends up in the subsequently formed OC is much less than that found in Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82] for the same cluster densities. This implies that for comets of roughly 2 km in size, the presence of the primordial solar nebula hinders OC formation. A byproduct of some of our simulations are endresults with a substantial fraction of the comets in the Uranus-Neptune scattered disk. A subsequent followup of this material is planned for the near future. In order to determine the effect of the size of the comets on OC formation efficiency, a set of runs with the same initial conditions but different cometary radii have been performed as well, from which it is determined that the threshold comet size to begin producing significant Oort clouds is roughly 20 km. This implies that the presence of the primordial solar nebula acts as a size-sorting mechanism, with large bodies unaffected by the gas drag and ending up in the OC while small bodies remain trapped in the planetary region, in the models studied.     

Radial diffusion rate of planetesimals due to gravitational encounters

We study the rate of radial diffusion of planetesimals due to mutual gravitational encounters under Hill’s approximations in the three-body problem. Planetesimals orbiting a central star radially migrate inward and outward as a result of mutual gravitational encounters and transfer angular momentum. We calculate the viscosity in a disk of equal-sized planetesimals due to their mutual gravitational encounters using three-body orbital integrations, and obtain a semianalytic expression that reproduces the numerical results. We find that the viscosity is independent of the velocity dispersion of planetesimals when the velocity dispersion is so small that Kepler shear dominates planetesimals’ relative velocities. On the other hand, in high-velocity cases where random velocities dominate the relative velocities, the viscosity is a decreasing function of the velocity dispersion, and is found to agree with previous estimates under the two-body approximation neglecting the solar gravity. We also calculate the rate of radial diffusion of planetesimals due to gravitational scattering by a massive protoplanet. Using these results, we discuss a condition for formation of nonuniform radial surface density distribution of planetesimals by gravitational perturbation of an embedded protoplanet.     

Formation of gas giant planets: core accretion models with fragmentation and planetary envelope

We have calculated formation of gas giant planets based on the standard core accretion model including effects of fragmentation and planetary envelope. The accretion process is found to proceed as follows. As a result of runaway growth of planetesimals with initial radii of ∼10 km, planetary embryos with a mass of ∼10²⁷ g (∼ Mars mass) are found to form in ∼10⁵ years at Jupiter's position (5.2 AU), assuming a large enough value of the surface density of solid material (25 g/cm²) in the accretion disk at that distance. Strong gravitational perturbations between the runaway planetary embryos and the remaining planetesimals cause the random velocities of the planetesimals to become large enough for collisions between small planetesimals to lead to their catastrophic disruption. This produces a large number of fragments. At the same time, the planetary embryos have envelopes, that reduce energies of fragments by gas drag and capture them. The large radius of the envelope increases the collision rate between them, resulting in rapid growth of the planetary embryos. By the combined effects of fragmentation and planetary envelope, the largest planetary embryo with 21M_⊕ forms at 5.2 AU in 3.8×10⁶ years. The planetary embryo is massive enough to start a rapid gas accretion and forms a gas giant planet.     

The distributions and ages of refractory objects in the solar nebula

Refractory objects such as Calcium, Aluminum-rich Inclusions, Amoeboid Olivine Aggregates, and crystalline silicates, are found in primitive bodies throughout our Solar System. It is believed that these objects formed in the hot, inner solar nebula and were redistributed during the mass and angular momentum transport that took place during its early evolution. The ages of these objects thus offer possible clues about the timing and duration of this transport. Here we study how the dynamics of these refractory objects in the evolving solar nebula affected the age distribution of the grains that were available to be incorporated into planetesimals throughout the Solar System. It is found that while the high temperatures and conditions needed to form these refractory objects may have persisted for millions of years, it is those objects that formed in the first 10⁵ years that dominate (make up over 90%) those that survive throughout most of the nebula. This is due to two effects: (1) the largest numbers of refractory grains are formed at this time period, as the disk is rapidly drained of mass during subsequent evolution and (2) the initially rapid spreading of the disk due to angular momentum transport helps preserve this early generation of grains as opposed to later generations. This implies that most refractory objects found in meteorites and comets formed in the first 10⁵ years after the nebula formed. As these objects contained live ²⁶Al, this constrains the time when short-lived radionuclides were introduced to the Solar System to no later than 10⁵ years after the nebula formed. Further, this implies that the t=0 as defined by meteoritic materials represents at most, the instant when the solar nebula finished accreting significant amounts of materials from its parent molecular cloud.     

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.     

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.     

N-Body simulations of growth from 1 km planetesimals at 0.4 AU

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M_⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10⁵ particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.     

Oligarchic growth with migration and fragmentation

In the core-accretion model, giant-planet cores form by oligarchic growth from a population of planetesimals prior to the dispersal of the disk gas. Once a core reaches a critical mass of roughly 10 Earth masses, it begins to accrete a gaseous envelope, forming a giant planet. Collisions between planetesimals cause fragmentation. Planetesimal fragments are more easily captured by cores, speeding up growth, but fragments are also lost by radial drift, reducing the total solid mass in the disk. Interaction with the gas causes cores to undergo inward type-I migration. Migration allows a core to accrete planetesimals from a larger region, but migrating cores may be lost if they reach the star. Thus, migration and fragmentation have both a positive and a negative impact on core formation. Here we describe results of new simulations of oligarchic growth that include fragmentation and/or migration. In the absence of migration, cores grow until they reach their isolation mass, which increases with distance from the star, or until the disk gas disperses. Fragmentation increases the maximum core mass by increasing growth rates in the outer disk, allowing objects to reach their isolation mass during the disk lifetime. When migration is present, cores migrate inwards rapidly when they approach 1 Earth mass. Most migrating cores are lost. Migrating cores gain little extra mass since they are passing through regions that have been depleted by earlier generations of cores. For a disk viscosity parameter alpha=1e−3 and planetesimal radius = 10 km, the maximum core mass is roughly 4 and 0.5 Earth masses with/without fragmentation, respectively, with little dependence on the disk mass. Formation and survival of 10-Earth-mass cores, in the presence of migration, requires large alpha (1e−2) and a massive disk (0.1 solar masses). When alpha is large, type-I migration rates decrease rapidly with time, allowing large, late-forming cores to survive. The addition of a stochastic (random-walk) migration component makes little difference to the outcome, provided that stochastic migration affects only cores larger than 0.01 Earth masses. Stochastic migration becomes increasingly important if it also affects lower-mass objects.     

The influence of grain rotation on the structure of dust aggregates

We study the effect of rotation during the collision between dust aggregates, in order to address a mismatch between previous model calculations of Brownian motion driven aggregation and experiments. We show that rotation during the collision does influence the shape and internal structure of the aggregates formed. The effect is limited in the ballistic regime when aggregates can be considered to move on straight lines during a collision. However, if the stopping length of an aggregate becomes smaller than its physical size, extremely elongated aggregates can be produced. We show that this effect may have played a role in the inner regions of the solar nebula where densities were high.     

Making the Earth: Combining dynamics and chemistry in the Solar System

No terrestrial planet formation simulation completed to date has considered the detailed chemical composition of the planets produced. While many have considered possible water contents and late veneer compositions, none have examined the bulk elemental abundances of the planets produced as an important check of formation models. Here we report on the first study of this type. Bulk elemental abundances based on disk equilibrium studies have been determined for the simulated terrestrial planets of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58]. These abundances are in excellent agreement with observed planetary values, indicating that the models of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58] are successfully producing planets comparable to those of the Solar System in terms of both their dynamical and chemical properties. Significant amounts of water are accreted in the present simulations, implying that the terrestrial planets form “wet” and do not need significant water delivery from other sources. Under the assumption of equilibrium controlled chemistry, the biogenic species N and C still need to be delivered to the Earth as they are not accreted in significant proportions during the formation process. Negligible solar photospheric pollution is produced by the planetary formation process. Assuming similar levels of pollution in other planetary systems, this in turn implies that the high metallicity trend observed in extrasolar planetary systems is in fact primordial.     

Models of the collisional damping scenario for ice-giant planets and Kuiper belt formation

Chiang et al. [Chiang, E., Lithwick, Y., Murray-Clay, R., Buie, M., Grundy, W., Holman, M., 2007. In: Protostars and Planets V, pp. 895-911] have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ∼5 primordial ice-giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich et al. [Goldreich, P., Lithwick, Y., Sari, R., 2004. Astrophys. J. 614, 497-507; Annu. Rev. Astron. Astrophys. 42, 549-601]. Here we test this hypothesis with a series of numerical simulations using a new code designed to incorporate the dynamical effects of collisions. We find that we cannot reproduce the observed Solar System. In particular, Goldreich et al. [Goldreich, P., Lithwick, Y., Sari, R., 2004. Astrophys. J. 614, 497-507; Annu. Rev. Astron. Astrophys. 42, 549-601] and Chiang et al. [Chiang, E., Lithwick, Y., Murray-Clay, R., Buie, M., Grundy, W., Holman, M., 2007. In: Protostars and Planets V, pp. 895-911] argue that during the instability, all but two of the ice giants would be ejected from the Solar System by Jupiter and Saturn, leaving Uranus and Neptune behind. We find that ejections are actually rare and that instead the systems spread outward. This always leads to a configuration with too many planets that are too far from the Sun. Thus, we conclude that both Goldreich et al.'s scheme for the formation of Uranus and Neptune and Chiang et al.'s Kuiper belt formation scenario are not viable in their current forms.     

The contribution of small grains to the opacity of protoplanetary atmospheres

I compute the opacity of grains in a protoplanetary atmosphere. The grain size distribution at different levels in the atmosphere is calculated using a simple microphysical model of grain growth via collisions and destruction via vaporization at high temperatures. The Rosseland mean opacity of the resulting distribution is then computed. For most cases examined, the grain opacity is significantly lower than earlier estimates.     

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.