No Arabic abstract
We compare the planet-to-star mass-ratio distribution measured by gravitational microlensing to core accretion theory predictions from population synthesis models. The core accretion theorys runaway gas accretion process predicts a dearth of intermediate-mass giant planets that is not seen in the microlensing results. In particular, the models predict $sim10,times$ fewer planets at mass ratios of $10^{-4} leq q leq 4 times 10^{-4}$ than inferred from microlensing observations. This tension implies that gas giant formation may involve processes that have hitherto been overlooked by existing core accretion models or that the planet-forming environment varies considerably as a function of host-star mass. Variation from the usual assumptions for the protoplanetary disk viscosity and thickness could reduce this discrepancy, but such changes might conflict with microlensing results at larger or smaller mass ratios, or with other observations. The resolution of this discrepancy may have important implications for planetary habitability because it has been suggested that the runaway gas accretion process may have triggered the delivery of water to our inner solar system. So, an understanding of giant planet formation may help us to determine the occurrence rate of habitable planets.
Thanks to recent high resolution ALMA observations, there is an accumulating evidence for presence of giant planets with masses from $sim 0.01$ Jupiter mass to a few Jupiter mass with separations up to $ 100$~AU in the annular structures observed in young protoplanetary discs. We point out that these observations set unique live constraints on the process of gas accretion onto sub-Jovian planets that were not previously available. Accordingly, we use a population synthesis approach in a new way: we build time-resolved models and compare the properties of the synthetic planets with the ALMA data at the same age. Applying the widely used gas accretion formulae leads to a deficit of sub-Jovian planets and an over-abundance of a few Jupiter mass planets compared to observations. We find that gas accretion rate onto planets needs to be suppressed by about an order of magnitude to match the observed planet mass function. This slower gas giant growth predicts that the planet mass should correlate positively with the age of the protoplanetary disc, albeit with a large scatter. This effect is not clearly present in the ALMA data but may be confirmed in the near future with more observations.
(abridged) 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. 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 R_0. 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 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. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius. 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.
We conduct the first microlensing simulation in the context of planet formation model. The planet population is taken from the Ida & Lin core accretion model for $0.3M_odot$ stars. With $6690$ microlensing events, we find for a simplified Korea Microlensing Telescopes Network (KMTNet) the fraction of planetary events is $2.9%$ , out of which $5.5%$ show multiple-planet signatures. The number of super-Earths, super-Neptunes and super-Jupiters detected are expected to be almost equal. Our simulation shows that high-magnification events and massive planets are favored by planet detections, which is consistent with previous expectation. However, we notice that extremely high-magnification events are less sensitive to planets, which is possibly because the 10 min sampling of KMTNet is not intensive enough to capture the subtle anomalies that occur near the peak. This suggests that while KMTNet observations can be systematically analyzed without reference to any follow-up data, follow-up observations will be essential in extracting the full science potential of very high-magnification events. The uniformly high-cadence observations expected for KMTNet also result in $sim 55%$ of all detected planets being non-caustic-crossing, and more low-mass planets even down to Mars-mass being detected via planetary caustics. We also find that the distributions of orbital inclinations and planet mass ratios in multiple-planet events agree with the intrinsic distributions.
After protoplanets have acquired sufficient mass to open partial gaps in their natal protostellar disks, residual gas continues to diffuse onto horseshoe streamlines under effect of viscous dissipation, and meander in and out of the planets Hill sphere. Within the Hill sphere, the horseshoe streamlines intercept gas flow in circumplanetary disks. The host stars tidal perturbation induces a barrier across the converging streamlines interface. Viscous transfer of angular momentum across this tidal barrier determines the rate of mass diffusion from the horseshoe streamlines onto the circumplanetary disks, and eventually the accretion rate onto the protoplanets. We carry out a series of numerical simulations to test the influence of this tidal barrier on super thermal planets. In weakly viscous disks, protoplanets accretion rate steeply decreases with their masses above the thermal limit. As their growth timescale exceeds the gas depletion time scale, their masses reach asymptotic values comparable to that of Jupiter. In relatively thick and strongly viscous disks, protoplanets asymptotic masses exceed several times that of Jupiter. Two dimensional numerical simulations show that such massive protoplanets strongly excite the eccentricity of nearby horseshoe streamlines, destabilize orderly flow, substantially enhance the diffusion rate across the tidal barrier, and elevate their growth rate until their natal disk is severely depleted. In contrast, eccentric streamlines remain stable in three dimensional simulations. Based on the upper falloff in the observe mass distribution of known exoplanets, we suggest their natal disks had relatively low viscosity alpha sim 0.001, modest thickness H/R sim 0.03 to 0.05, and limited masses comparable to that of minimum mass solar nebula model.
Context: Planet formation by pebble accretion is an alternative to planetesimal-driven core accretion. In this scenario, planets grow by accreting cm-to-m-sized pebbles instead of km-sized planetesimals. One of the main differences with planetesimal-driven core accretion is the increased thermal ablation experienced by pebbles. This provides early enrichment to the planets envelope, which changes the process of core growth. Aims: We aim to predict core masses and envelope compositions of planets that form by pebble accretion and compare mass deposition of pebbles to planetesimals. Methods: We model the early growth of a proto-planet by calculating the structure of its envelope, taking into account the fate of impacting pebbles or planetesimals. The region where high-Z material can exist in vapor form is determined by the vapor pressure. We include enrichment effects by locally modifying the mean molecular weight. Results: In the pebble case, three phases of core growth can be identified. In the first phase, pebbles impact the core without significant ablation. During the second phase, ablation becomes increasingly severe. A layer of high-Z vapor starts to form around the core that absorbs a small fraction of the ablated mass. The rest either rains out to the core or mixes outwards instead, slowing core growth. In the third phase, the vapor inner region expands outwards, absorbing an increasing fraction of the ablated material as vapor. Rainout ends before the core mass reaches 0.6 M_Earth, terminating direct core growth. In the case of icy H2O pebbles, this happens before 0.1 M_Earth. Conclusions: Our results indicate that pebble accretion can directly form rocky cores up to only 0.6 M_Earth, and is unable to form similarly sized icy cores. Subsequent core growth can proceed indirectly when the planet cools, provided it is able to retain its high-Z material.