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 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.     

Embedded star clusters and the formation of the Oort cloud: III. Evolution of the inner cloud during the Galactic phase

In a previous publication [Brasser, R., Duncan, M.J., Levison, H.F., 2006. Icarus 184, 59-82], models of the inner Oort cloud were built which included the effect of an embedded star cluster on cometary orbits about the Sun. The main conclusions of that paper were that the formation efficiency is about 10% and the median distance of the cloud to the Sun only depends on the mean density of gas and stars the Sun encountered. Here we report on the results of simulations which followed the ensuing dynamical evolution of these comet clouds in the current Galactic environment once the Sun left the embedded star cluster. The goal is to determine whether or not the dynamical influence of passing Galactic field stars and the Galactic tidal field is sufficient to replenish the current outer cloud (semi-major axis a>20,000 AU) with enough material from the inner cloud (a<20,000 AU). Since visible new comets come directly from the outer cloud, a mass estimate only exists for the latter, with a lower limit of 1 M_⊕ [Francis, P.J., 2005. Astrophys. J. 635, 1348-1361]. Knowing the amount of expansion of the inner cloud may therefore yield an estimate of the mass of said (unseen) inner cloud. Our results indicate that typically only 10% of the comets from the inner cloud land in the outer cloud and are bound after 4.5 Gyr. If one assumes that in the extreme case all or the majority of the current population of the outer cloud has come from the inner cloud, then a typical value of the mass of the inner cloud is about 10 M_⊕. The results of [Brasser, R., Duncan, M.J., Levison, H.F., 2006. Icarus 184, 59-82] showed that ∼10% of comets from the Jupiter-Saturn region were implanted in the inner Oort cloud, which implies an uncomfortably large value of about 100 M_⊕ for the mass of solids in the primordial Jupiter-Saturn region. This extreme case might be remedied in two says: either the effect of Giant Molecular Cloud complexes on the inner Oort cloud must be much more severe than originally thought, or there was a two-stage formation process for the Oort cloud, in which the outer cloud was largely populated by comets scattered once the Sun had left its primordial birth cluster.     

Embedded star clusters and the formation of the Oort Cloud

Observations suggest most stars originate in clusters embedded in giant molecular clouds [Lada, C.J., Lada, E.A., 2003. Annu. Rev. Astron. Astrophys. 41, 57-115]. Our Solar System likely spent 1-5 Myrs in such regions just after it formed. Thus the Oort Cloud (OC) possibly retains evidence of the Sun's early dynamical history and of the stellar and tidal influence of the cluster. Indeed, the newly found objects (90377) Sedna and 2000 CR₁₀₅ may have been put on their present orbits by such processes [Morbidelli, A., Levison, H.F., 2004. Astron. J. 128, 2564-2576]. Results are presented here of numerical simulations of the orbital evolution of comets subject to the influence of the Sun, Jupiter and Saturn (with their current masses on orbits appropriate to the period before the Late Heavy Bombardment (LHB) [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]), passing stars and tidal force associated with the gas and stars of an embedded star cluster. The cluster was taken to be a Plummer model with 200-400 stars, with a range of initial central densities. The Sun's orbit was integrated in the cluster potential together with Jupiter and Saturn and the test particles. Stellar encounters were incorporated by directly integrating the effects of stars passing within a sphere centred on the Sun of radius equal to the Plummer radius for low-density clusters and half a Plummer radius for high-density clusters. The gravitational influence of the gas was modeled using the tidal force of the cluster potential. For a given solar orbit, the mean density, 〈ρ〉, was computed by orbit-averaging the density of material encountered. This parameter proved to be a good measure for predicting the properties of the OC. On average 2-18% of our initial sample of comets end up in the OC after 1-3 Myr. A comet is defined to be part of the OC if it is bound and has q>35 AU. Our models show that the median distance of an object in the OC scales approximately as 〈ρ〉^−1/2 when . Our models easily produce objects on orbits like that of (90377) Sedna [Brown, M.E., Trujillo, C., Rabinowitz, D., 2004. Astrophys. J. 617, 645-649] within ∼1 Myr in cases where the mean density is or higher; one needs mean densities of order to create objects like 2000 CR₁₀₅ by this mechanism, which are reasonable (see, e.g., Guthermuth, R.A., Megeath, S.T., Pipher, J.L., Williams, J.P., Allen, L.E., Myers, P.C., Raines, S.N., 2005. Astrophys. J. 632, 397-420). Thus the latter object may also be part of the OC. Close stellar passages can stir the primordial Kuiper Belt to sufficiently high eccentricities (e?0.05; Kenyon, S.J., Bromley, B.C., 2002. Astron. J. 123, 1757-1775) that collisions become destructive. From the simulations performed it is determined that there is a 50% or better chance to stir the primordial Kuiper Belt to eccentricities e?0.05 at 50 AU when . The orbit of the new object 2003 UB₃₁₃ [Brown, M.E., Trujillo, C.A., Rabinowitz, D.L., 2005. Astrophys. J. 635, L97-L100] is only reproduced for mean cluster densities of the order of , but in the simulations it could not come to be on its current orbit by this mechanism without disrupting the formation of bodies in the primordial Kuiper Belt down to 20 AU. It is therefore improbable that the latter object is created by this mechanism.     

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 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.     

Orbit computation for transneptunian objects

Using statistical orbital ranging, we systematically study the orbit computation problem for transneptunian objects (TNOs). We have automated orbit computation for large numbers of objects, and, more importantly, we are able to obtain orbits even for the most sparsely observed objects (observational arcs of a few days). For such objects, the resulting orbit distributions include a large number of high-eccentricity orbits, in which TNOs can be perturbed by close encounters with Neptune. The stability of bodies on the computed orbits has therefore been ascertained by performing a study of close encounters with the major planets. We classify TNO orbit distributions statistically, and we study the evolution of their ephemeris uncertainties. We find that the orbital element distributions for the most numerous single-apparition TNOs do not support the existence of a postulated sharp edge to the belt beyond 50 AU. The technique of statistical ranging provides ephemeris predictions more generally than previously possible also for poorly observed TNOs.     

The fossilized size distribution of the main asteroid belt

Planet formation models suggest the primordial main belt experienced a short but intense period of collisional evolution shortly after the formation of planetary embryos. This period is believed to have lasted until Jupiter reached its full size, when dynamical processes (e.g., sweeping resonances, excitation via planetary embryos) ejected most planetesimals from the main belt zone. The few planetesimals left behind continued to undergo comminution at a reduced rate until the present day. We investigated how this scenario affects the main belt size distribution over Solar System history using a collisional evolution model (CoEM) that accounts for these events. CoEM does not explicitly include results from dynamical models, but instead treats the unknown size of the primordial main belt and the nature/timing of its dynamical depletion using innovative but approximate methods. Model constraints were provided by the observed size frequency distribution of the asteroid belt, the observed population of asteroid families, the cratered surface of differentiated Asteroid (4) Vesta, and the relatively constant crater production rate of the Earth and Moon over the last 3 Gyr. Using CoEM, we solved for both the shape of the initial main belt size distribution after accretion and the asteroid disruption scaling law . In contrast to previous efforts, we find our derived function is very similar to results produced by numerical hydrocode simulations of asteroid impacts. Our best fit results suggest the asteroid belt experienced as much comminution over its early history as it has since it reached its low-mass state approximately 3.9-4.5 Ga. These results suggest the main belt's wavy-shaped size-frequency distribution is a “fossil” from this violent early epoch. We find that most diameter D?120 km asteroids are primordial, with their physical properties likely determined during the accretion epoch. Conversely, most smaller asteroids are byproducts of fragmentation events. The observed changes in the asteroid spin rate and lightcurve distributions near D∼100-120 km are likely to be a byproduct of this difference. Estimates based on our results imply the primordial main belt population (in the form of D<1000 km bodies) was 150-250 times larger than it is today, in agreement with recent dynamical simulations.     

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.     

The dynamics of Jupiter and Saturn in the gaseous protoplanetary disk

We study the possibility that the mutual interactions between Jupiter and Saturn prevented Type II migration from driving these planets much closer to the Sun. Our work extends previous results by Masset and Snellgrove [Masset, F., Snellgrove, M., 2001. Mon. Not. R. Astron. Soc. 320, L55-L59], by exploring a wider set of initial conditions and disk parameters, and by using a new hydrodynamical code that properly describes for the global viscous evolution of the disk. Initially both planets migrate towards the Sun, and Saturn's migration tends to be faster. As a consequence, they eventually end up locked in a mean motion resonance. If this happens in the 2:3 resonance, the resonant motion is particularly stable, and the gaps opened by the planets in the disk may overlap. This causes a drastic change in the torque balance for the two planets, which substantially slows down the planets' inward migration. If the gap overlap is substantial, planet migration may even be stopped or reversed. As the widths of the gaps depend on disk viscosity and scale height, this mechanism is particularly efficient in low viscosity, cool disks. The initial locking of the planets in the 2:3 resonance is a likely outcome if Saturn formed at the edge of Jupiter's gap, but also if Saturn initially migrated rapidly from further away. We also explore the possibility of trapping in other resonances, and the subsequent evolutions. We discuss the compatibility of our results with the initial conditions adopted in Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] and Gomes et al. [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469] to explain the current orbital architecture of the giant planets and the origin of the Late Heavy Bombardment of the Moon.     

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.     

The evolution of the water distribution in a viscous protoplanetary disk

Astronomical observations have shown that protoplanetary disks are dynamic objects through which mass is transported and accreted by the central star. This transport causes the disks to decrease in mass and cool over time, and such evolution is expected to have occurred in our own solar nebula. Age dating of meteorite constituents shows that their creation, evolution, and accumulation occupied several Myr, and over this time disk properties would evolve significantly. Moreover, on this timescale, solid particles decouple from the gas in the disk and their evolution follows a different path. It is in this context that we must understand how our own solar nebula evolved and what effects this evolution had on the primitive materials contained within it. Here we present a model which tracks how the distribution of water changes in an evolving disk as the water-bearing species experience condensation, accretion, transport, collisional destruction, and vaporization. Because solids are transported in a disk at different rates depending on their sizes, the motions will lead to water being concentrated in some regions of a disk and depleted in others. These enhancements and depletions are consistent with the conditions needed to explain some aspects of the chemistry of chondritic meteorites and formation of giant planets. The levels of concentration and depletion, as well as their locations, depend strongly on the combined effects of the gaseous disk evolution, the formation of rapidly migrating rubble, and the growth of immobile planetesimals. Understanding how these processes operate simultaneously is critical to developing our models for meteorite parent body formation in the Solar System and giant planet formation throughout the galaxy. We present examples of evolution under a range of plausible assumptions and demonstrate how the chemical evolution of the inner region of a protoplanetary disk is intimately connected to the physical processes which occur in the outer regions.     

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.     

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.     

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.     

Rotation rate and velocity dispersion of planetary ring particles with size distribution II. Numerical simulation for gravitating particles

We examine rotation rates of gravitating particles in low optical depth rings, on the basis of the evolution equation of particle rotational energy derived by Ohtsuki [Ohtsuki, K., 2006. Rotation rate and velocity dispersion of planetary ring particles with size distribution. I. Formulation and analytic calculation. Icarus 183, 373-383]. We obtain the rates of evolution of particle rotation rate and velocity dispersion, using three-body orbital integration that takes into account distribution of random velocities and rotation rates. The obtained stirring and friction rates are used to calculate the evolution of velocity dispersion and rotation rate for particles in one- and two-size component rings as well as those with a narrow size distribution, and agreement with N-body simulation is confirmed. Then, we perform calculations to examine equilibrium rotation rates and velocity dispersion of gravitating ring particles with a broad size distribution, from 1 cm up to 10 m. We find that small particles spin rapidly with 〈ω²〉^1/2/Ω?¹⁰²-¹⁰³, where ω and Ω are the particle rotation rate and its orbital angular frequency, respectively, while the largest particles spin slowly, with 〈ω²〉^1/2/Ω?1. The vertical scale height of rapidly rotating small particles is much larger than that of slowly rotating large particles. Thus, rotational states of ring particles have vertical heterogeneity, which should be taken into account in modeling thermal infrared emission from Saturn's rings.     

Rotation rate and velocity dispersion of planetary ring particles with size distribution: I. Formulation and analytic calculation

We derive an equation for the evolution of rotational energy of Keplerian particles in a dilute disk due to mutual collisions. Three-dimensional Keplerian motion of particles is taken into account precisely, on the basis of Hill's approximation. The Rayleigh distribution of particles' orbital eccentricities and inclinations, and the Gaussian distribution of their rotation rates are also taken into account. Performing appropriate variable transformation, we show that the equation can be expressed with two terms. The first term, which we call collisional stirring term, represents energy exchange between rotation and random motion via collisions. The second term, which we call rotational friction term, tends to equalize the mean rotational energy of particles with different sizes. The equation can describe the evolution of rotational energy of Keplerian particles with an arbitrary size distribution. We analytically evaluate the rates of stirring and friction for the random kinetic energy and rotational energy due to inelastic collisions, for non-gravitating particles in a dilute disk. Using these results, we discuss equilibrium states in a disk of spinning, non-gravitating Keplerian particles.     

Evolution of a Keplerian disk of colliding and fragmenting particles: a kinetic model with application to the Edgeworth-Kuiper belt

We present a kinetic model of a disk of solid particles, orbiting a primary and experiencing inelastic collisions. In distinction to other collisional models that use a 2D (mass-semimajor axis) binning and perform a separate analysis of the velocity (eccentricity, inclination) evolution, we choose mass and orbital elements as independent variables of a phase space. The distribution function in this space contains full information on the combined mass, spatial, and velocity distributions of particles. A general kinetic equation for the distribution function is derived, valid for any set of orbital elements and for any collisional outcome, specified by a single kernel function. The first implementation of the model utilizes a 3D phase space (mass-semimajor axis-eccentricity) and involves averages over the inclination and all angular elements. We assume collisions to be destructive, simulate them with available material- and size-dependent scaling laws, and include collisional damping. A closed set of kinetic equations for a mass-semimajor axis-eccentricity distribution is written and transformation rules to usual mass and spatial distributions of the disk material are obtained. The kinetic “core” of our approach is generic. It is possible to add inclination as an additional phase space variable, to include cratering collisions and agglomeration, dynamical friction and viscous stirring, gravity of large perturbers, drag forces, and other effects into the model. As a specific application, we address the collisional evolution of the classical population in the Edgeworth-Kuiper belt (EKB). We run the model for different initial disk's masses and radial profiles and different impact strengths of objects. Our results for the size distribution, collisional timescales, and mass loss are in agreement with previous studies. In particular, collisional evolution is found to be most substantial in the inner part of the EKB, where the separation size between the survivors over EKB's age and fragments of earlier collisions lies between a few and several tens of km. The size distribution in the EKB is not a single Dohnanyi-type power law, reflecting the size dependence of the critical specific energy in both strength and gravity regimes. The net mass loss rate of an evolved disk is nearly constant and is dominated by disruption of larger objects. Finally, assuming an initially uniform distribution of orbital eccentricities, we show that an evolved disk contains more objects in orbits with intermediate eccentricities than in nearly circular or more eccentric orbits. This property holds for objects of any size and is explained in terms of collisional probabilities. The effect should modulate the eccentricity distribution shaped by dynamical mechanisms, such as resonances and truncation of perihelia by Neptune.     

The primordial excitation and clearing of the asteroid belt—Revisited

We have performed new simulations of two different scenarios for the excitation and depletion of the primordial asteroid belt, assuming Jupiter and Saturn on initially circular orbits as predicted by the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. First, we study the effects of sweeping secular resonances driven by the depletion of the solar nebula. We find that these sweeping secular resonances are incapable of giving sufficient dynamical excitation to the asteroids for nebula depletion timescales consistent with estimates for solar-type stars, and in addition cannot cause significant mass depletion in the asteroid belt or produce the observed radial mixing of different asteroid taxonomic types. Second, we study the effects of planetary embryos embedded in the primordial asteroid belt. These embedded planetary embryos, combined with the action of jovian and saturnian resonances, can lead to dynamical excitation and radial mixing comparable to the current asteroid belt. The mass depletion driven by embedded planetary embryos alone, even in the case of an eccentric Jupiter and Saturn, is roughly 10-20× less than necessary to explain the current mass of the main belt, and thus a secondary depletion event, such as that which occurs naturally in the Nice Model, is required. We discuss the implications of our new simulations for the dynamical and collisional evolution of the main belt.