共查询到20条相似文献，搜索用时 343 毫秒
1.
We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (106-107 years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 105 years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested. 相似文献
2.
John Chambers 《Icarus》2006,180(2):496-513
A new semi-analytic model for the oligarchic growth phase of planetary accretion is developed. The model explicitly calculates damping and excitation of planetesimal eccentricities e and inclinations i due to gas drag and perturbations from embryos. The effects of planetesimal fragmentation, enhanced embryo capture cross sections due to atmospheres, inward planetesimal drift, and embryo-embryo collisions are also incorporated. In the early stages of oligarchic growth, embryos grow rapidly as e and i fall below their equilibrium values. The formation of planetesimal collision fragments also speeds up embryo growth as fragments have low-e, low-i orbits, thereby optimizing gravitational focussing. At later times, the presence of thick atmospheres captured from the nebula aids embryo growth by increasing their capture cross sections. Planetesimal drift due to gas drag can lead to substantial inward transport of solid material. However, inward drift is greatly reduced when embryo atmospheres are present, as the drift timescale is no longer short compared to the accretion timescale. Embryo-embryo collisions increase embryo growth rates by 50% compared to the case where growth is solely due to accretion of planetesimals. Formation of 0.1-Earth-mass protoplanets at 1 AU and 10-Earth-mass cores at 5 AU requires roughly 0.1 and 1 million years respectively, in a nebula where the local solid surface density is 7 g cm−2 at each of these locations. 相似文献
3.
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 107 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. 相似文献
4.
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 ∼1027 g (∼ Mars mass) are found to form in ∼105 years at Jupiter's position (5.2 AU), assuming a large enough value of the surface density of solid material (25 g/cm2) 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×106 years. The planetary embryo is massive enough to start a rapid gas accretion and forms a gas giant planet. 相似文献
5.
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 R0. 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. 相似文献
6.
Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M⊙, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities. 相似文献
7.
The validity of the two-body approximation in calculating encounters between planetesimals has been evaluated as a function of the ratio of unperturbed planetesimal velocity (with respect to a circular orbit) to mutual escape velocity when their surfaces are in contact (V/Ve). Impact rates as a function of V/Ve are calculated to within ～20% by numerical integration of the equations of motion. It is found that when V/Ve > 0.4, the two-body approximation is a good one. At low velocities (V/Ve < 0.1) two-body “collision-course” trajectories fail to lead to impacts. On the other hand, at these low velocities many impacts result from encounter trajectories with unperturbed separation distances far beyond the two-body gravitational radius. These two effects tend to cancel, and the resulting impact rates remain within a factor of ～3 of the two-body value in spite of these major differences in the nature of the impact trajectories. Therefore, on the average, the two-body approximation is useful well below the value of V/Ve for which it fails to describe individual encounters, and the required corrections are not large. As a consequence of this “anomalous gravitational focusing” planetesimals will continue to interact even when their orbits are noncrossing. This reduces the difficulty with premature isolation of planetesimal embryos during accumulation. Quantitatively, when 0.06 ? V/Ve ? 0.2, the impact rate varies approximately with the fifth power of the radius of the larger body, and is about a factor of 3 above that predicted using the conventional two-body gravitational cross-section formula. At lower values of V/Ve , the impact rate increases less rapidly. Finally, at the lowest values of V/Ve (<.02), the impact rate increases only in proportion to the geometric cross section, as a consequence of the swarm being essentially two dimensional for large unperturbed encounter distances. The gravitational enhancement in effective cross section is thereby limited to a value of about 3000. This leads to an optimal size for growth of planetesimals from a swarm of given eccentricity, and places a limit on the extent of runaway accretion. 相似文献
8.
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 105 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. 相似文献
9.
This paper investigates the surface density evolution of a planetesimal disk due to the effect of type-I migration by carrying out N-body simulation and through analytical method, focusing on terrestrial planet formation. The coagulation and the growth of the planetesimals take place in the abundant gas disk except for a final stage. A protoplanet excites density waves in the gas disk, which causes the torque on the protoplanet. The torque imbalance makes the protoplanet suffer radial migration, which is known as type-I migration. Type-I migration time scale derived by the linear theory may be too short for the terrestrial planets to survive, which is one of the major problems in the planet formation scenario. Although the linear theory assumes a protoplanet being in a gas disk alone, Kominami et al. [Kominami, J., Tanaka, H., Ida, S., 2005. Icarus 167, 231-243] showed that the effect of the interaction with the planetesimal disk and the neighboring protoplanets on type-I migration is negligible. The migration becomes pronounced before the planet's mass reaches the isolation mass, and decreases the solid component in the disk. Runaway protoplanets form again in the planetesimal disk with decreased surface density. In this paper, we present the analytical formulas that describe the evolution of the solid surface density of the disk as a function of gas-to-dust ratio, gas depletion time scale and semimajor axis, which agree well with our results of N-body simulations. In general, significant depletion of solid material is likely to take place in inner regions of disks. This might be responsible for the fact that there is no planet inside Mercury's orbit in our Solar System. Our most important result is that the final surface density of solid components (Σd) and mass of surviving planets depend on gas surface density (Σg) and its depletion time scale (τdep) but not on initial Σd; they decrease with increase in Σg and τdep. For a fixed gas-to-dust ratio and τdep, larger initial Σd results in smaller final Σd and smaller surviving planets, because of larger Σg. To retain a specific amount of Σd, the efficient disk condition is not an initially large Σd but the initial Σd as small as the specified final one and a smaller gas-to-dust ratio. To retain Σd comparable to that of the minimum mass solar nebula (MMSN), a disk must have the same Σd and a gas-to-dust ratio that is smaller than that of MMSN by a factor of 1.3×(τdep/1 Myr) at ∼1 AU. (Equivalently, type-I migration speed is slower than that predicted by the linear theory by the same factor.) The surviving planets are Mars-sized ones in this case; in order to form Earth-sized planets, their eccentricities must be pumped up to start orbit crossing and coagulation among them. At ∼5 AU, Σd of MMSN is retained under the same condition, but to form a core massive enough to start runaway gas accretion, a gas-to-dust ratio must be smaller than that of MMSN by a factor of 3×τdep/1 Myr. 相似文献
10.
The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System. 相似文献
11.
P. Thebault 《Celestial Mechanics and Dynamical Astronomy》2011,111(1-2):29-49
HD 196885 Ab is the most ??extreme?? planet-in-a-binary discovered to date, whose orbit places it at the limit for orbital stability. The presence of a planet in such a highly perturbed region poses a clear challenge to planet-formation scenarios. We investigate this issue by focusing on the planet-formation stage that is arguably the most sensitive to binary perturbations: the mutual accretion of kilometre-sized planetesimals. To this effect we numerically estimate the impact velocities dv amongst a population of circumprimary planetesimals. We find that most of the circumprimary disc is strongly hostile to planetesimal accretion, especially the region around 2.6 AU (the planet??s location) where binary perturbations induce planetesimal-shattering dv of more than 1 kms?1. Possible solutions to the paradox of having a planet in such accretion-hostile regions are (1) that initial planetesimals were very big, at least 250 km (2) that the binary had an initial orbit at least twice the present one, and was later compacted due to early stellar encounters (3) that planetesimals did not grow by mutual impacts but by sweeping of dust (the ??snowball?? growth mode identified by Xie et al., in Astrophys J 724:1153, 2010b), or (4) that HD 196885 Ab was formed not by core-accretion but by the concurrent disc instability mechanism. All of these 4 scenarios remain however highly conjectural. 相似文献
12.
We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×1024 g), and vary in mass between 0.23M⊕ and 3.85M⊕. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M⊕ protoplanets which can serve as the accretion seeds for 2+M⊕ planets with large water contents. 相似文献
13.
14.
Edward R.D. Scott 《Icarus》2006,185(1):72-82
Thermal models and radiometric ages for meteorites show that the peak temperatures inside their parent bodies were closely linked to their accretion times. Most iron meteorites come from bodies that accreted <0.5 Myr after CAIs formed and were melted by 26Al and 60Fe, probably inside 2 AU. Rare carbon-rich differentiated meteorites like ureilites probably also come from bodies that formed <1 Myr after CAIs, but in the outer part of the asteroid belt. Chondrite groups accreted intermittently from diverse batches of chondrules and other materials over a 4 Myr period starting 1 Myr after CAI formation when planetary embryos may already have formed at ∼1 AU. Meteorite evidence precludes accretion of late-forming chondrites on the surface of early-formed bodies; instead chondritic and non-chondritic meteorites probably formed in separate planetesimals. Maximum metamorphic temperatures in chondrite groups are correlated with mean chondrule age, as expected if 26Al and 60Fe were the predominant heat sources. Because late-forming bodies could not accrete close to large, early-formed bodies, planetesimal formation may have spread across the nebula from regions where the differentiated bodies formed. Dynamical models suggest that the asteroids could not have accreted in the main belt if Jupiter formed before the asteroids. Therefore Jupiter probably reached its current mass >3-5 Myr after CAIs formed. This precludes formation of Jupiter via a gravitational instability <1 Myr after the solar nebula formed, and strongly favors core accretion. Jupiter probably formed too late to make chondrules by generating shocks directly, or indirectly by scattering Ceres-sized bodies across the belt. Nevertheless, shocks formed by gravitational instabilities or Ceres-sized bodies scattered by planetary embryos may have produced some chondrules. The minimum lifetime for the solar nebula of 3-5 Myr inferred from the total spread of CAI and chondrule ages may exceed the median lifetime of 3 Myr for protoplanetary disks, but is well within the 1-10 Myr observed range. Shorter formation times for extrasolar planets may help to explain their unusual orbits compared to those of solar giant planets. 相似文献
15.
Richard Greenberg 《Earth, Moon, and Planets》1980,22(1):63-66
Safronov's (1972) demonstration that relative velocities of planetesimals would be comparable to the dominant size bodies' escape velocities, combined with a plausible size distribution that has most mass in the largest bodies, yielded his evolution model with limited growth of the largest planetesimal with respect to its next largest neighbors. A numerical simulation of planetesimal accretion (Greenberget al., 1978) suggests that at least over one stage of collisional accretion, velocities were much lower than the escape velocity of the largest bodies, because the bulk of the mass still resided in km-scale bodies. The low velocities at this early stage may conceivably have permitted early runaway growth, which, in turn, would have kept the velocities low and permitted continued runaway growth of the largest bodies.Paper presented at the European Workshop on Planetary Sciences, organised by the Laboratorio di Astrofisica Spaziale di Frascati, and held between April 23–27, 1979, at the Accademia Nazionale del Lincei in Rome, Italy. 相似文献
16.
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. 相似文献
17.
S.J. Weidenschilling 《Icarus》2011,214(2):671-684
The present size frequency distribution (SFD) of bodies in the asteroid belt appears to have preserved some record of the primordial population, with an excess of bodies of diameter D ∼ 100 km relative to a simple power law. The survival of Vesta’s basaltic crust also implies that the early SFD had a shallow slope in the range ∼10-100 km. (Morbidelli, A., Bottke, W.F., Nesvorny, D., Levison, H.F. [2009]. Icarus 204, 558-573) were unable to produce these features by accretion from an initial population of km-sized planetesimals. They concluded that bodies with sizes in the range ∼100-1000 km and a SFD similar to the current population were produced directly from solid particles of sub-meter scale, without experiencing accretion through intermediate sizes. We present results of new accretion simulations in the primordial asteroid region. The requisite SFD can be produced from an initial population of planetesimals of sizes ?0.1 km, smaller than the usual assumption of km-sized bodies. The bump at D ∼ 100 km is produced by a transition from dispersion-dominated runaway growth to a regime dominated by Keplerian shear, before the formation of large protoplanetary embryos. Thus, accretion of the asteroids from an initial population of small (sub-km) planetesimals cannot be ruled out. 相似文献
18.
The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation. 相似文献
19.
The effects of metallicity and grain growth and settling on the early evolution of gaseous protoplanets 总被引：1，自引：0，他引：1
Giant protoplanets formed by gravitational instability in the outer regions of circumstellar disks go through an early phase of quasi-static contraction during which radii are large (∼1 AU) and internal temperatures are low (<2000 K). The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of isolated, non-accreting planets of 3, 5, and 7 MJ. First, we pick three different overall metallicities for the planet and simply scale the opacity accordingly. We show that higher metallicity results in slower contraction as a result of higher opacity. It is found that the pre-collapse time scale is proportional to the metallicity. In this scenario, survival of giant planets formed by gravitational instability is predicted to be more likely around low-metallicity stars, since they evolve to the point of collapse to small size on shorter time scales. But metal-rich planets, as a result of longer contraction times, have the best opportunity to capture planetesimals and form heavy-element cores. Second, we investigate the effects of opacity reduction as a result of grain growth and settling, for the same three planetary masses and for three different values of overall metallicity. When these processes are included, the pre-collapse time scale is found to be of order 1000 years for the three masses, significantly shorter than the time scale calculated without these effects. In this case the time scale is found to be relatively insensitive to planetary mass and composition. However, the effects of planetary rotation and accretion of gas and dust, which could increase the timescale, are not included in the calculation. The short time scale we find would preclude metal enrichment by planetesimal capture, as well as heavy-element core formation, over a large range of planetary masses and metallicities. 相似文献
20.
The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their ‘pre-collapse’ stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be and 5.67×103 years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions. 相似文献