N-BODY SIMULATION OF PLANETESIMAL FORMATION THROUGH GRAVITATIONAL INSTABILITY AND COAGULATION. II. ACCRETION MODEL

Journal Article · Thu Oct 01 04:00:00 EDT 2009 · Astrophysical Journal

DOI:https://doi.org/10.1088/0004-637X/703/2/1363· OSTI ID:21371902

Michikoshi, Shugo; Kokubo, Eiichiro ^[1]; Inutsuka, Shu-ichiro ^[2]

---

1. Center for Computational Astrophysics, National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588 (Japan)
2. Department of Physics, Kyoto University, Kyoto 606-8502 (Japan)

The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local N-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of nonaxisymmetric structures; the creation of planetesimal seeds; and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to L {sup 3/2} {sub y}, where L{sub y} is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of L{sub x} , where L{sub x} is the size of the computational domain in the radial direction if L{sub x} is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.

• OSTI ID: 21371902
• Journal Information: Astrophysical Journal, Journal Name: Astrophysical Journal Journal Issue: 2 Vol. 703; ISSN ASJOAB; ISSN 0004-637X
• Country of Publication: United States
• Language: English

Related Subjects

79 ASTRONOMY AND ASTROPHYSICS
COLLISIONS
DUSTS
GRAVITATION
GRAVITATIONAL INSTABILITY
INSTABILITY
MASS
MOTION
PLANETS
PLASMA INSTABILITY
PROTOPLANETS
ROTATION
SATELLITES
SIMULATION