No Arabic abstract
Many massive objects have been found in the outer region of the Solar system. How they were formed and evolved has not been well understood, although there have been intensive studies on accretion process of terrestrial planets. One of the mysteries is the existence of binary planetesimals with near-equal mass components and highly eccentric orbits. These binary planetesimals are quite different from the satellites observed in the asteroid belt region. The ratio of the Hill radius to the physical radius of the planetesimals is much larger for the outer region of the disk, compared to the inner region of the disk. The Hill radius increases with the semi major axis. Therefore, planetesimals in the outer region can form close and eccentric binaries, while those in the inner region would simply collide. In this paper, we carried out $N$-body simulations in different regions of the disk and studied if binaries form in the outer region of the disk. We found that large planetesimals tend to form binaries. A significant fraction of large planetesimals are components of the binaries. Planetesimals that become the components of binaries eventually collide with a third body, through three-body encounters. Thus, the existence of binaries can enhance the growth rate of planetesimals in the Trans-Neptunian Object (TNO) region.
Most of planet formation models that incorporate planetesimal fragmentation consider a catastrophic impact energy threshold for basalts at a constant velocity of 3 km/s during all the process of the formation of the planets. However, as planets grow the relative velocities of the surrounding planetesimals increase from velocities of the order of m/s to a few km/s. In addition, beyond the ice line where giant planets are formed, planetesimals are expected to be composed roughly by 50 percentage of ices. We aim to study the role of planetesimal fragmentation on giant planet formation considering planetesimal catastrophic impact energy threshold as a function of the planetesimal relative velocities and compositions. We improve our model of planetesimal fragmentation incorporating a functional form of the catastrophic impact energy threshold with the planetesimal relative velocities and compositions. We also improve in our model the accretion of small fragments produced by the fragmentation of planetesimals during the collisional cascade considering specific pebble accretion rates. We find that a more accurate and realistic model for the calculation of the catastrophic impact energy threshold tends to slow down the formation of massive cores. Only for reduced grain opacity values at the envelope of the planet, the cross-over mass is achieved before the disk time-scale dissipation. While planetesimal fragmentation favors the quick formation of massive cores of 5-10 Earth masses the cross-over mass could be inhibited by planetesimal fragmentation. However, grain opacity reduction or pollution by the accreted planetesimals together with planetesimal fragmentation could explain the formation of giant planets with low-mass cores.
Planetesimals are compact astrophysical objects roughly 1-1000 km in size, massive enough to be held together by gravity. They can grow by accreting material to become full-size planets. Planetesimals themselves are thought to form by complex physical processes from small grains in protoplanetary disks. The streaming instability (SI) model states that mm/cm-size particles (pebbles) are aerodynamically collected into self-gravitating clouds which then directly collapse into planetesimals. Here we analyze ATHENA simulations of the SI to characterize the initial properties (e.g., rotation) of pebble clouds. Their gravitational collapse is followed with the PKDGRAV N-body code, which has been modified to realistically account for pebble collisions. We find that pebble clouds rapidly collapse into short-lived disk structures from which planetesimals form. The planetesimal properties depend on the clouds scaled angular momentum, l=L/(M R_H^2 Omega, where L and M are the angular momentum and mass, R_H is the Hill radius, and Omega is the orbital frequency. Low-l pebble clouds produce tight (or contact) binaries and single planetesimals. Compact high-l clouds give birth to binary planetesimals with attributes that closely resemble the equal-size binaries found in the Kuiper belt. Significantly, the SI-triggered gravitational collapse can explain the angular momentum distribution of known equal-size binaries -- a result pending verification from studies with improved resolution. About 10% of collapse simulations produce hierarchical systems with two or more large moons. These systems should be found in the Kuiper belt when observations reach the threshold sensitivity.
The equation of state calculated by Saumon and collaborators has been adopted in most core-accretion simulations of giant-planet formation performed to date. Since some minor errors have been found in their original paper, we present revised simulations of giant-planet formation that considers a corrected equation of state. We employ the same code as Fortier and collaborators in repeating our previous simulations of the formation of Jupiter. Although the general conclusions of Fortier and collaborators remain valid, we obtain significantly lower core masses and shorter formation times in all cases considered. The minor errors in the previously published equation of state have been shown to affect directly the adiabatic gradient and the specific heat, causing an overestimation of both the core masses and formation times.
We study the collisional evolution of km-sized planetesimals in tight binary star systems to investigate whether accretion towards protoplanets can proceed despite the strong gravitational perturbations from the secondary star. The orbits of planetesimals are numerically integrated in two dimensions under the influence of the two stars and gas drag. The masses and orbits of the planetesimals are allowed to evolve due to collisions with other planetesimals and accretion of collisional debris. In addition, the mass in debris can evolve due to planetesimal-planetesimal collisions and the creation of new planetesimals. We show that it is possible in principle for km-sized planetesimals to grow by two orders of magnitude in size if the efficiency of planetesimal formation is relatively low. We discuss the limitations of our two-dimensional approach.
(Abridged) We present local 2D and 3D hybrid numerical simulations of particles and gas in the midplane of protoplanetary disks (PPDs) using the Athena code. The particles are coupled to gas aerodynamically, with particle-to-gas feedback included. Magnetorotational turbulence is ignored as an approximation for the dead zone of PPDs, and we ignore particle self-gravity to study the precursor of planetesimal formation. Our simulations include a wide size distribution of particles, ranging from strongly coupled particles with dimensionless stopping time tau_s=Omega t_stop=1e-4 to marginally coupled ones with tau_s=1 (where Omega is the orbital frequency, t_stop is the particle friction time), and a wide range of solid abundances. Our main results are: 1. Particles with tau_s>=0.01 actively participate in the streaming instability, generate turbulence and maintain the height of the particle layer before Kelvin-Helmholtz instability is triggered. 2. Strong particle clumping as a consequence of the streaming instability occurs when a substantial fraction of the solids are large (tau_s>=0.01) and when height-integrated solid to gas mass ratio Z is super-solar. 3. The radial drift velocity is reduced relative to the conventional Nakagawa-Sekiya-Hayashi (NSH) model, especially at high Z. We derive a generalized NSH equilibrium solution for multiple particle species which fits our results very well. 4. Collision velocity between particles with tau_s>=0.01 is dominated by differential radial drift, and is strongly reduced at larger Z. 5. There exist two positive feedback loops with respect to the enrichment of local disk solid abundance and grain growth. All these effects promote planetesimal formation.