No Arabic abstract
In this Thesis I studied the formation of the four giant planets of the Solar System in the framework of the nucleated instability hypothesis. The model considers that solids and gas accretion are coupled in an interactive fashion, taking into account detailed constitutive physics for the envelope. The accretion rate of the core corresponds to the oligarchic growth regime. I also considered that accreted planetesimals follow a size distribution. One of the main results of this Thesis is that I was able to compute the formation of Jupiter, Saturn, Uranus and Neptune in less than 10 million years, which is considered to be the protoplanetary disk mean lifetime.
The terrestrial planets are believed to have formed by violent collisions of tens of lunar- to Mars-size protoplanets at time t<200 Myr after the protoplanetary gas disk dispersal (t_0). The solar system giant planets rapidly formed during the protoplanetary disk stage and, after t_0, radially migrated by interacting with outer disk planetesimals. An early (t<100 Myr) dynamical instability is thought to have occurred with Jupiter having gravitational encounters with a planetary-size body, jumping inward by ~0.2-0.5 au, and landing on its current, mildly eccentric orbit. Here we investigate how the giant planet instability affected formation of the terrestrial planets. We study several instability cases that were previously shown to match many solar system constraints. We find that resonances with the giant planets help to remove solids available for accretion near ~1.5 au, thus stalling the growth of Mars. It does not matter, however, whether the giant planets are placed on their current orbits at t_0 or whether they realistically evolve in one of our instability models; the results are practically the same. The tight orbital spacing of Venus and Earth is difficult to reproduce in our simulations, including cases where bodies grow from a narrow annulus at 0.7-1 au, because protoplanets tend to radially spread during accretion. The best results are obtained in the narrow-annulus model when protoplanets emerging from the dispersing gas nebula are assumed to have (at least) the Mars mass. This suggests efficient accretion of the terrestrial protoplanets during the first ~10 Myr of the solar system.
The solar systems dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the Early Instability scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar systems terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system.
Forming gas giant planets by the accretion of 100 km diameter planetesimals, a typical size that results from self-gravity assisted planetesimal formation, is often thought to be inefficient. Many models therefore use small km-sized planetesimals, or invoke the accretion of pebbles. Furthermore, models based on planetesimal accretion often use the ad hoc assumption of planetesimals distributed radially in a minimum mass solar nebula fashion. We wish to investigate the impact of various initial radial density distributions in planetesimals with a dynamical model for the formation of planetesimals on the resulting population of planets. In doing so, we highlight the directive role of the early stages of dust evolution into pebbles and planetesimals in the circumstellar disk on the following planetary formation. We have implemented a two population model for solid evolution and a pebble flux regulated model for planetesimal formation into our global model for planet population synthesis. This framework is used to study the global effect of planetesimal formation on planet formation. As reference, we compare our dynamically formed planetesimal surface densities with ad-hoc set distributions of different radial density slopes of planetesimals. Even though required, it is not solely the total planetesimal disk mass, but the planetesimal surface density slope and subsequently the formation mechanism of planetesimals, that enables planetary growth via planetesimal accretion. Highly condensed regions of only 100 km sized planetesimals in the inner regions of circumstellar disks can lead to gas giant growth. Pebble flux regulated planetesimal formation strongly boosts planet formation, because it is a highly effective mechanism to create a steep planetesimal density profile. We find this to lead to the formation of giant planets inside 1 au by 100 km already by pure planetesimal accretion.
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 the protoplanetary disc surrounding AB Aurigae have revealed the possible presence of two giant planets in the process of forming. The young measured age of $1-4$Myr for this system allows us to place strict time constraints on the formation histories of the observed planets. Hence we may be able to make a crucial distinction between formation through core accretion (CA) or the gravitational instability (GI), as CA formation timescales are typically Myrs whilst formation through GI will occur within the first $approx10^4-10^5$yrs of disc evolution. We focus our analysis on the $4-13$M$_{rm Jup}$ planet observed at $Rapprox30$AU. We find CA formation timescales for such a massive planet typically exceed the systems age. The planets high mass and wide orbit may instead be indicative of formation through GI. We use smoothed particle hydrodynamic simulations to determine the systems critical disc mass for fragmentation, finding $M_{rm d,crit}=0.3$M$_{odot}$. Viscous evolution models of the discs mass history indicate that it was likely massive enough to exceed $M_{rm d,crit}$ in the recent past, thus it is possible that a young AB Aurigae disc may have fragmented to form multiple giant gaseous protoplanets. Calculations of the Jeans mass in an AB Aurigae-like disc find that fragments may initially form with masses $1.6-13.3$M$_{rm Jup}$, consistent with the planets which have been observed. We therefore propose that the inferred planets in the disc surrounding AB Aurigae may be evidence of planet formation through GI.