We study the formation of the protoplanetary disk by the collapse of a primordial molecular cloud, and how its evolution leads to the selection of specific types of planets. We use a hydrodynamical code that accounts for the dynamics, thermodynamics, geometry, and composition of the disk to numerically model its evolution as it is fed by the infalling cloud material. As the mass accretion rate of the disk onto the star determines its growth, we can calculate the stellar characteristics by interpolating its radius, luminosity, and temperature over the stellar mass from pre-calculated stellar evolution models. The density and midplane temperature of the disk then allow us to model the interactions between the disk and potential planets and determine their migration. At the end of the collapse phase, when the disk reaches its maximum mass, it pursues its viscous spreading, similarly to the evolution from a minimum mass solar nebula (MMSN). In addition, we establish a timeline equivalence between the MMSN and a "collapse-formed disk" that would be older by about 2 Myr. We can save various types of planets from a fatal type-I inward migration: in particular, planetary embryos can avoid falling on the star by becoming trapped at the heat transition barriers and at most sublimation lines (except the silicates one). One of the novelties concerns the possible trapping of putative giant planets around a few astronomical units from the star around the end of the infall. Moreover, trapped planets may still follow the traps outward during the collapse phase and inward after it. Finally, this protoplanetary disk formation model shows the early possibilities of trapping planetary embryos at disk stages that are anterior by a few million years to the initial state of the MMSN approximation.

Planet traps are necessary to prevent forming planets from falling onto their host star by type I migration. Surface mass density and temperature gradient irregularities favor the apparition of traps and deserts. Such features are found at the dust sublimation lines and heat transition barriers. We study how planets may remain trapped or escape as they grow and as the disk evolves. We model the temporal viscous evolution of a protoplanetary disk by coupling its dynamics, thermodynamics, geometry and composition. The resulting mid-plane density and temperature profiles allow the modeling of the interactions of such an evolving disk with potential planets, even before the steady state is reached. We follow the viscous evolution of a MMSN and compute the Lindblad and corotation torques that such a disk would exert on potential planets of various masses located within the planetary formation region. We determine the position of planet traps and deserts in relationship with the sublimation lines, shadowed regions and heat transition barriers. Planets that are a few tens of Earth masses can be trapped at the sublimation lines until they reach a certain mass while planets more massive than 100ME can only be trapped permanently at the heat transition barriers. Coupling a bimodal planetary migration model with a self-consistent evolved disk, we were able to distinguish several potential planet populations after 5 million years of evolution: two populations of giant planets that could stay trapped around 5.5 and 9 au and possibly open gaps, some super-Earths trapped around 5 and 7.5 au and a population of close-in super-Earths trapped inside 1 au. The traps corresponding to the last group could help validating the in-situ formation scenarios of the observed close-in super-Earths.

Aims. We track the time evolution of planet traps and snowlines in a viscously evolving protoplanetary disk using an opacity table that accounts for the composition of the dust material. Methods. We coupled a dynamical and thermodynamical disk model with a temperature-dependent opacity table (that accounts for the sublimation of the main dust components) to investigate the formation and evolution of snowlines and planet traps during the first million years of disk evolution. Results. Starting from a minimum mass solar nebula (MMSN), we find that the disk mid-plane temperature profile shows several plateaux (0.1-1 AU wide) at the different sublimation temperatures of the species that make up the dust. For water ice, the correspond- ing plateau can be larger than 1 AU, which means that this is a snow "region" rather than a snow "line". As a consequence, the surface density of solids in the snow region may increase gradually, not abruptly. Several planet traps and desert regions appear naturally as a result of abrupt local changes in the temperature and density profiles over the disk lifetime. These structures are mostly located at the edges of the temperature plateaux (surrounding the dust sublimation lines) and at the heat-transition barrier where the disk stellar heating and viscous heating are of the same magnitude (around 10 AU after 1 Myr). Conclusions. Several traps are identified: although a few appear to be transient, most of them slowly migrate along with the heat- transition barrier or the dust sublimation lines. These planet traps may temporarily favor the growth of planetary cores.

Observations of protoplanetary disks show that some characteristics seem recurrent, even in star formation regions that are physically distant such as surface mass density profiles varying as $r^{-1}$, or aspect ratios about 0.03 to 0.23. Accretion rates are also recurrently found around $10^{-8} - 10^{-6}~\mathrm{M_{\odot}~yr^{-1}}$ for disks already evolved (Isella et al., 2009, Andrews et al., 2009, 2010). Several models have been developed in order to recover these properties. However, most of them usually simplify the disk geometry if not its mid-plane temperature. This has major consequences for modeling the disk evolution over million years and consequently planet migration. In the present paper, we develop a viscous evolution hydrodynamical numerical code that determines simultaneously the disk photosphere geometry and the mid-plane temperature. We then compare our results of long-term simulations with similar simulations of disks with a constrained geometry along the Chiang & Goldreich (1997) prescription (dlnH/dlnr = 9/7). We find that the constrained geometry models provide a good approximation of the disk surface density evolution. However, they differ significantly regarding the temperature time evolution. In addition, we find that shadowed regions naturally appear at the transition between viscously dominated and radiation dominated regions that falls in the region of planetary formation. We show that $\chi$ (photosphere height to pressure scale height ratio) cannot be considered as a constant, consistently with Watanabe et al. (2008). Comparisons with observations show that all disk naturally evolve toward a shallow surface density disk ($\Sigma \propto r^{-1}$). The mass flux across the disk stabilizes in about 1 million year typically.