mass (solid curve) and the mean mass (dashed curve) of the system. thanthisrangearenotstatisticallyvalidsinceeachmassbinoften has only a few bodies. First, the distribution tends to relax to a తͷ༷ࢠ ฏۉ ࠷େͷఱମ ඍͷత ɹˠݪ͕࢝ੜ͢Δ 20 KOKUBO AND IDA FIG. 3. Snapshots of a planetesimal system on the a–e plane. The circles represent planetesimals and their radii are proportional to the radii of planetesi- mals. The system initially consists of 3000 equal-mass (1023 g) planetesimals. FIG. 4. Time evolution of the maximum mass (solid curve) and the mean mass (dashed curve) of the system. thanthisrangearenotstatisticallyvalidsinceeachmassbinoften has only a few bodies. First, the distribution tends to relax to a decreasing function of mass through dynamical friction among (energy equipartition of) bodies (t = 50,000, 100,000 years). Second, the distributions tend to flatten (t = 200,000 years). This is because as a runaway body grows, the system is mainly heated by the runaway body (Ida and Makino 1993). In this case, the eccentricity and inclination of planetesimals are scaled by the يಓܘ<"6> يಓ৺ ࣭ྔ<H> ࣌ؒ<> <,PLVCP*EB >
of a planetesimal system on the a–e plane. The cir- cles represent planetesimals and their radii are proportional to the radii of planetesimals. The system initially consists of 4000 planetesimals whose to- tal mass is 1.3 × 1027 g. The initial mass distribution is given by the power- FIG. 8. The number of bodies in linear mass bins is plotted for t = 100,000, 200,000, 300,000, 400,000, and 500,000 years. In Fig. 10, we plot the maximum mass and the mean mass of يಓ৺ ֤ॴͰඍ͕త ɹˠαΠζͷݪ͕࢝ฒͿ ՉతͱΑͿ ʹ ֤يಓͰͷݪ࢝ ࣭ྔ [kg] ܗ࣌ؒ [yr] ٿيಓ 1×1024 7×105 يಓ 3×1025 4×107 ఱԦيಓ 8×1025 2×109 يಓܘ<"6> <,PLVCP*EB >
0:6, which means that the typical result- ing system consists of two Earth-sized planets and a smaller planet. In this model, we obtain hna i ’ 1:8 Æ 0:7. In other words, one or two planets tend to form outside the initial distribution of protoplanets. In most runs, these planets are smaller scattered planets. Thus we obtain a high efficiency of h fa i ¼ 0:79 Æ 0:15. The accretion timescale is hTacc i ¼ 1:05 Æ 0:58 ð Þ ; 108 yr. These results are consistent with Agnor et al. (1999), whose initial con- clus sam ha1 i larg are 0:06 hM2 0:05 Fig. 2.—Snapshots of the system on the a-e (left) and a-i (right) planes at t ¼ 0, 106, 10 are proportional to the physical sizes of the planets. KOKUBO, KOMINAM 1134 ͍࣌ؒΛ͔͚ͯݪ࢝ಉ࢜ͷيಓ͕ཚΕΔ ɹˠޓ͍ʹিಥɾ߹ମͯ͠ΑΓେ͖ͳఱମʹ <,PLVCP*EB > (c) Hidenori Genda