No Arabic abstract
In the classical core-accretion planet formation scenario, rapid inward migration and accretion timescales of kilometer size planetesimals may not favor the formation of massive cores of giant planets before the dissipation of protoplanetary disks. On the other hand, the existence of pressure maxima in the disk could act as migration traps and locations for solid material accumulation, favoring the formation of massive cores. We aim to study the radial drift of pebbles and planetesimals and planet migration at pressure maxima in a protoplanetary disk and their implications for the formation of massive cores as triggering a gaseous runaway accretion phase. The time evolution of a viscosity driven accretion disk is solved numerically introducing a a dead zone as a low-viscosity region in the protoplanetary disk. A population of pebbles and planetesimals evolving by radial drift and accretion by the planets is also considered. Finally, the embryos embedded in the disk grow by the simultaneous accretion of pebbles, planetesimals and the surrounding gas. Our simulations show that the pressure maxima generated at the edges of the low-viscosity region of the disk act as planet migration traps, and that the pebble and planetesimal surface densities are significantly increased due to the radial drift towards pressure maxima locations. However, our simulations also show that migration trap locations and solid material accumulation locations are not exactly at the same positions. Thus, a planets semi-major axis oscillations around zero torque locations, predicted by MHD and HD simulations, are needed for the planet to accrete all the available material accumulated at the pressure maxima. Pressure maxima generated at the edges of a low-viscosity region of a protoplanetary disk seem to be preferential locations for the formation and trap of massive cores.
Recent observations of protoplanetary disks have revealed ring-like structures that can be associated to pressure maxima. Pressure maxima are known to be dust collectors and planet migration traps. Most of planet formation works are based either on the pebble accretion model or on the planetesimal accretion model. However, recent studies proposed the possible formation of Jupiter by the hybrid accretion of pebbles and planetesimals. We aim to study the full process of planet formation consisting of dust evolution, planetesimal formation and planet growth at a pressure maximum in a protoplanetary disk. We compute, through numerical simulations, the gas and dust evolution, including dust growth, fragmentation, radial drift and particle accumulation at a pressure bump. We also consider the formation of planetesimals by streaming instability and the formation of a moon-size embryo that grows into a giant planet by the hybrid accretion of pebbles and planetesimals. We find that pressure maxima in protoplanetary disks are efficient collectors of dust drifting inwards. The condition of planetesimal formation by streaming instability is fulfilled due to the large amount of dust accumulated at the pressure bump. Then, a massive core is quickly formed (in $sim 10^4$ yr) by the accretion of pebbles. After the pebble isolation mass is reached, the growth of the core slowly continues by the accretion of planetesimals. The energy released by planetesimal accretion delays the onset of runaway gas accretion, allowing a gas giant to form after $sim$1 Myr of disk evolution. The pressure maximum also acts as a migration trap. Pressure maxima in protoplanetary disks are preferential locations for dust traps, planetesimal formation by streaming instability and planet migration traps. All these conditions allow the fast formation of a giant planet by the hybrid accretion of pebbles and planetesimals.
Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.
The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems.
In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-type stars have been observed. Here, we quantitatively estimate the total mass of ejected materials during the giant impact stages. We found that $sim$0.4 times the Earths mass is ejected in total throughout the giant impact stage. Ejected materials are ground down by collisional cascade until micron-sized grains are blown out by radiation pressure. The depletion timescale of these ejected materials is determined primarily by the mass of the largest body among them. We conducted high-resolution simulations of giant impacts to accurately obtain the mass of the largest ejected body. We then calculated the evolution of the debris disks produced by a series of giant impacts and depleted by collisional cascades to obtain the infrared excess evolution of the debris disks. We found that the infrared excess is almost always higher than the stellar infrared flux throughout the giant impact stage ($sim$100 Myr) and is sometimes $sim$10 times higher immediately after a giant impact. Therefore, giant impact stages would explain the infrared excess from most observed warm debris disks. The observed fraction of stars with warm debris disks indicates that the formation probability of our solar system-like terrestrial planets is approximately 10%.
In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $sim 1.0 {rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($sim 4.0$ mm) being a factor of $sim10$ larger than other vortex regions. Both large ($sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.