SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

The harvest of exoplanet discoveries has opened the area of exoplanet characterisation. But this cannot be achieved without a careful analysis of the host star parameters. The system of HD219134 hosts two transiting exoplanets and at least two additional non-transiting exoplanets. We used the VEGA/CHARA interferometer to measure the angular diameter of HD219134, leading to a stellar radius of $R_{\star}=0.726\pm0.014 R_{\odot}$. We also derived the stellar density from the transits light curves ($\rho_{\star}=1.82\pm0.19 \rho_{\odot}$), which finally gives a direct estimate of the mass ($M_{\star}=0.696\pm0.078 M_{\odot}$) with a correlation of 0.46 between $R_{\star}$ and $M_{\star}$. This new mass is smaller than that derived from the C2kSMO stellar evolutionary model, which provides a mass range of 0.755$-$0.810 ($\pm 0.040$) $M_{\odot}$. This allows us to infer the mass, radius and density of the two transiting exoplanets of the system. We then use an inference model to obtain the internal parameters of these two transiting exoplanets. Moreover, we find that planet $b$ and $c$ have smaller radii than previously estimated ($1.500\pm0.057$ and $1.415\pm0.049 R_{\oplus}$, respectively); this clearly puts these planets out of the gap in the exoplanetary radii distribution and validates their super-Earth nature. Planet $b$ is more massive than planet $c$, but possibly less dense. We investigate whether this could be caused by partial melting of the mantle and find that tidal heating due to non-zero eccentricity of planet $b$ may be powerful enough. The system of HD219134 constitutes a very valuable benchmark for both stellar physics and exoplanetary science. The direct determination of the stellar density, radius and mass should be more extensively applied to provide accurate exoplanets properties and calibrate stellar models.

Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 100 Jupiter masses per Myr. Such high accretion rates would result in all planets with larger-than-10-Earth-mass cores forming Jupiter-like planets, in clear contrast to the ice giants in the Solar System and the observed exoplanet population. In this work, we use 3D hydrodynamical simulations, that include radiative transfer, to model the growth of the envelope on planets with different masses. We find that gas flows rapidly through the outer part of the envelope, but this flow does not drive accretion. Instead, gas accretion is the result of quasi-static contraction of the inner envelope, which can be orders of magnitude smaller than the mass flow through the outer atmosphere. For planets smaller than Saturn, we measure moderate gas accretion rates that are below 1 Jupiter mass per Myr. Higher mass planets, however, accrete up to 10 times faster and do not reveal a self-driven mechanism that can halt gas accretion. Therefore, the reason for the final masses of Saturn and Jupiter remains difficult to understand, unless their completion coincided with the dissipation of the Solar Nebula.

The characterization of exoplanets relies on that of their host star. However, stellar evolution models cannot always be used to derive the mass and radius of individual stars, because many stellar internal parameters are poorly constrained. Here, we use the probability density functions (PDFs) of directly measured parameters to derive the joint PDF of the stellar and planetary mass and radius. Because combining the density and radius of the star is our most reliable way of determining its mass, we find that the stellar (respectively planetary) mass and radius are strongly (respectively moderately) correlated. We then use a generalized Bayesian inference analysis to characterize the possible interiors of 55 Cnc e. We quantify how our ability to constrain the interior improves by accounting for correlation. The information content of the mass-radius correlation is also compared with refractory element abundance constraints. We provide posterior distributions for all interior parameters of interest. Given all available data, we find that the radius of the gaseous envelope is $0.08 \pm 0.05 R_p$. A stronger correlation between the planetary mass and radius (potentially provided by a better estimate of the transit depth) would significantly improve interior characterization and reduce drastically the uncertainty on the gas envelope properties.

In a recent paper (Crida et al. 2018, accepted on April 19), we presented a method to derive the mass and radius of a transiting exoplanet and their intrinsic correlation, that we applied to 55 Cnc e. We wrote: "More precise observations of the transit would be very useful in this particular case and would allow to increase significantly the gain on the planetary density precision." Three months later, Bourrier et al. (2018) published new observations of the system. Additionally, Gaia's DR2 was released on April 25, 2018. The purpose of this research note is solely to implement these precise new data in our pipeline to provide an up-to-date result of our model. We find $R_p = 1.947\pm 0.038\,R_\oplus$ and $M_p = 8.59\pm 0.43\,M_\oplus$. We also double the precision on the planetary density ($\rho_p = 1.164\pm 0.062\,\rho_\oplus$), which allows us to refine the estimates for all interior parameters of interest. In particular, we now find that the radius of the gaseous envelope is $0.03\pm 0.02\,R_p$.

Massive planets form within the lifetime of protoplanetary disks and undergo orbital migration due to planet-disk interactions. When the first planet reaches the inner edge of the disk its migration stops and the second planet is locked in resonance. We detail the resonant trapping comparing semi-analytical formulae and numerical simulations in the case of two equal-mass coplanar planets trapped in first order resonances. We describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using expanded and non-expanded Hamiltonians. We show that during convergent migration the planets evolve along these families of equilibria. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that the resonant configurations achieved by migration can become unstable after the removal of the disk. We probe the stability of the resonances as a function of planetary mass. To do this, we fictitiously increase the masses of the planets, adiabatically maintaining the low-amplitude libration regime until instability. We discuss two hypotheses for instability: low-order secondary resonance of the libration frequency with a fast synodic frequency, and minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low-amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter still play the decisive role in the destabilisation of resonant pairs. We show that, as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.

Context. Giant planets open gaps in their protoplanetary and subsequently suffer so-called type II migration. Schematically, planets are thought to be tightly locked within their surrounding disks, and forced to follow the viscous advection of gas onto the central star. This fundamental principle however has recently been questioned, as migrating planets were shown to decouple from the gas' radial drift. Aims. In this framework, we question whether the traditionally used linear scaling of migration rate of a giant planet with the disk's viscosity still holds. Additionally, we assess the role of orbit-crossing material as part of the decoupling mechanism. Methods. We have performed 2D (r, {\theta}) numerical simulations of point-mass planets embedded in locally isothermal {\alpha}-disks in steady-state accretion, with various values of {\alpha}. Arbitrary planetary accretion rates were used as a means to diminish or nullify orbit-crossing flows. Results. We confirm that the migration rate of a gap-opening planet is indeed proportional to the disk's viscosity, but is not equal to the gas drift speed in the unperturbed disk. We show that the role of gap-crossing flows is in fact negligible. Conclusions. From these observations, we propose a new paradigm for type II migration : a giant planet feels a torque from the disk that promotes its migration, while the gap profile relative to the planet is restored on a viscous timescale, thus limiting the planet migration rate to be proportional to the disk's viscosity. Hence, in disks with low viscosity in the planet region, type II migration should still be very slow. Key words. protoplanetary disks; planet-disk interactions; planets and satellites: formation

Rings are ubiquitous around giant planets in our Solar System. They evolve jointly with the nearby satellite system. They could form either during the giant planet formation process or much later, as a result of large scale dynamical instabilities either in the local satellite system, or at the planetary scale. We review here the main characteristics of rings in our solar system, and discuss their main evolution processes and possible origin. We also discuss the recent discovery of rings around small bodies.

The four massive Galilean satellites are believed to have formed within a circumplanetary disk during the last stages of Jupiter's formation. While the existence of a circum-jovian disk is supported by hydrodynamic simulations, no consensus exists regarding the origin and delivery mechanisms of the building blocks of the forming satellites. The opening of a gap in the circumsolar disk would have efficiently isolated Jupiter from the main sources of solid material. However, a reservoir of planetesimals should have existed at the outer edge of Jupiter's gap, where solids were trapped and accumulated over time. Here we show that the formation of Saturn's core within this reservoir, or its prompt inward migration, allows planetesimals to be redistributed from this reservoir towards Jupiter and the inner Solar System, thereby providing enough material to form the Galilean satellites and to populate the Main Belt with primitive asteroids. We find that the orbit of planetesimals captured within the circumjovian disk are circularized through friction with gas in a compact system comparable to the current radial extent of the Galilean satellites. The decisive role of Saturn in the delivery mechanism has strong implications for the occurrence of massive moons around extrasolar giant planets as they would preferentially form around planets within multiple planet systems.

The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be able to accrete gas from the protoplanetary disc and grow into a gas giant. Previous works have determined a scaling of the pebble isolation mass with cube of the disc aspect ratio. Here we expand on previous measurements and explore the dependency of the pebble isolation mass on all relevant parameters of the protoplanetary disc. We use 3D hydrodynamical simulations to measure the pebble isolation mass and derive a simple scaling law that captures the dependence on the local disc structure and the turbulent viscosity parameter $\alpha$. We find that small pebbles, coupled to the gas, with Stokes number $\tau_{\rm f}<0.005$ can drift through the partial gap at pebble isolation mass. However, as the planetary mass increases, particles must be decreasingly smaller to penetrate through the pressure bump. Turbulent diffusion of particles, however, can lead to an increase of the pebble isolation mass by a factor of two, depending on the strength of the background viscosity and on the pebble size. We finally explore the implications of the new scaling law of the pebble isolation mass on the formation of planetary systems by numerically integrating the growth and migration pathways of planets in evolving protoplanetary discs. Compared to models neglecting the dependence of the pebble isolation mass on the $\alpha$-viscosity, our models including this effect result in larger core masses for giant planets. These larger core masses are more akin to the core masses of the giant planets in the Solar System.

Magnetically-driven disk winds would alter the surface density slope of gas in the inner region of a protoplanetary disk $(r \lesssim 1 {\rm au})$. This in turn affects planet formation. Recently, the effect of disk wind torque has been considered with the suggestion that it would carve out the surface density of the disk from inside and would induce global gas flows (wind-driven accretion). We aim to investigate effects of global gas flows on type I migration and also examine planet formation. A simplified approach was taken to address this issue, and N-body simulations with isolation-mass planets were also performed. In previous studies, the effect of gas flow induced by turbulence-driven accretion has been taken into account for its desaturation effect of the corotation torque. If more rapid gas flows (e.g., wind-driven accretion) are considered, the desaturation effect can be modified. In MRI-inactive disks, in which the wind-driven accretion dominates the disk evolution, the gas flow at the midplane plays an important role. If this flow is fast, the corotation torque is efficiently desaturated. Then, the fact that the surface density slope can be positive in the inner region due to the wind torque can generate an outward migration region extended to super-Earth mass planets. In this case, we observe that no planets fall onto the central star in N-body simulations with migration forces imposed to reproduce such migration pattern. We also see that super-Earth mass planets can undergo outward migration. Relatively rapid gas flows affects type I migration and thus the formation of close-in planets.

This chapter of the book Planetary Ring Systems addresses the origin of planetary rings, one of the least understood processes related to planet formation and evolution. Whereas rings seem ubiquitous around giant planets, their great diversity of mass, structure and composition is a challenge for any formation scenario. Recent advances in our understanding of ring and satellite formation and destruction suggest that these processes are interconnected, so that rings and satellites may be two aspects of the same geological system. However, no single theory seems able to explain the origin of the different planetary rings known in our Solar System, and it now seems evident that rings may result from a variety of processes like giant collisions, tidal stripping of comets or satellites, as well as planet formation itself. In order to build any theory of ring formation it is important to specify physical processes that affect the long-term evolution of rings, as well as to describe the different observations that any ring formation model should explain. This is the topic of section 2. In section 3, we focus our attention on Saturn's rings and their main properties, and then discuss the pros and cons of a series of ring formation models. We also discuss the link between rings and satellites. In section 4, we extend the discussion to the other giant planets (Jupiter, Uranus, and Neptune). Section 5 is devoted to new types of rings -- the recent discovery of rings orbiting small outer Solar System bodies (Centaurs), and the possible rings around extrasolar planets. In section 6, we conclude and try to identify critical observations and theoretical advances needed to better understand the origin of rings and their significance in the global evolution of planets.

We aim to investigate the influence of the eccentricity and inclination damping due to planet-disc interactions on the final configurations of the systems, generalizing previous studies on the combined action of the gas disc and planet-planet scattering during the disc phase. Instead of the simplistic $K$-prescription, our n-body simulations adopt the damping formulae for eccentricity and inclination provided by the hydrodynamical simulations of our companion paper. We follow the evolution of $11000$ numerical experiments of three giant planets in the late stage of the gas disc, exploring different initial configurations, planetary mass ratios and disc masses. The dynamical evolutions of the planetary systems are studied along the simulations, with emphasis on the resonance captures and inclination-growth mechanisms. Most of the systems are found with small inclinations ($\le10^{\circ}$) at the dispersal of the disc. Even though many systems enter an inclination-type resonance during the migration, the disc usually damps the inclinations on a short timescale. Although the majority of the multiple systems in our results are quasi-coplanar, $\sim5\%$ of them end up with high mutual inclinations ($\ge10^{\circ}$). Half of these highly mutually inclined systems result from two- or three-body MMR captures, the other half being produced by orbital instability and/or planet-planet scattering. When considering the long-term evolution over $100$ Myr, destabilization of the resonant systems is common, and the percentage of highly mutually inclined systems still evolving in resonance drops to $30\%$. Finally, the parameters of the final system configurations are in very good agreement with the semi-major axis and eccentricity distributions in the observations, showing that planet-planet interactions during the disc phase could have played an important role in sculpting planetary systems.

Growing planets interact with their natal protoplanetary disc, which exerts a torque onto them allowing them to migrate in the disc. Small mass planets do not affect the gas profile and migrate in the fast type~I migration. Although type~I migration can be directed outwards for planets smaller than $20-30 M_\oplus$, planets above this mass should be lost into the central star long before the disc disperses. Massive planets open a gap and subsequently migrate in the slower, type~II migration, which could save them from migrating all the way to the star. Hence, growing giant planets can be saved if and only if they can reach the gap opening mass, because this extends their migration time-scale, allowing them to eventually survive at large orbits until the disc itself disperses. However, most of the previous studies only measured the torques on planets with fixed masses and orbits to determine the migration rate. Additionally, the transition between type~I and type~II migration itself is not well studied. Here we use isothermal 2D disc simulations with FARGO-2D1D to study the migration behaviour of gas accreting protoplanets in discs. We find that migrating giant planets always open gaps in the disc. We further show analytically and numerically that planets growing in the runaway gas accretion regime will reach gap opening masses before migrating all the way to the central star if the disc is not extremely viscous and/or thick. An accretion rate limited to the radial gas flow in the disc, in contrast, is not fast enough. When gas accretion by the planet is taken into account, the gap opening process is accelerated because the planet accretes material originating from its horseshoe region. This allows an accreting planet to transition to type~II migration before being lost even if gas fails to be provided for a rapid enough growth and the gap opening mass is not reached.

Aims: We investigate the evolution of protoplanetary discs (PPDs hereafter) with magnetically driven disc winds and viscous heating. Methods: We consider an initially massive disc with ~0.1 Msun to track the evolution from the early stage of PPDs. We solve the time evolution of surface density and temperature by taking into account viscous heating and the loss of the mass and the angular momentum by the disc winds within the framework of a standard alpha model for accretion discs. Our model parameters, turbulent viscosity, disc wind mass loss, and disc wind torque, which are adopted from local magnetohydrodynamical simulations and constrained by the global energetics of the gravitational accretion, largely depends on the physical condition of PPDs, particularly on the evolution of the vertical magnetic flux in weakly ionized PPDs. Results: Although there are still uncertainties concerning the evolution of the vertical magnetic flux remaining, surface densities show a large variety, depending on the combination of these three parameters, some of which are very different from the surface density expected from the standard accretion. When a PPD is in a "wind-driven accretion" state with the preserved vertical magnetic field, the radial dependence of the surface density can be positive in the inner region <1-10 au. The mass accretion rates are consistent with observations, even in the very low level of magnetohydrodynamical turbulence. Such a positive radial slope of the surface density gives a great impact on planet formation because (i)it inhibits the inward drift or even results in the outward drift of pebble/boulder-sized solid bodies, and (ii) it also makes the inward type-I migration of proto-planets slower or even reversed. Conclusions: The variety of our calculated PPDs should yield a wide variety of exoplanet systems.

Saturn formed beyond the snow line in the primordial solar nebula that made it possible for it to accrete a large mass. Disk instability and core accretion models have been proposed for Saturn's formation, but core accretion is favored on the basis of its volatile abundances, internal structure, hydrodynamic models, chemical characteristics of protoplanetary disk, etc. The observed frequency, properties and models of exoplanets provide additional supporting evidence for core accretion. The heavy elements with mass greater than 4He make up the core of Saturn, but are presently poorly constrained, except for carbon. The C/H ratio is super-solar, and twice that in Jupiter. The enrichment of carbon and other heavy elements in Saturn and Jupiter requires special delivery mechanisms for volatiles to these planets. In this chapter we will review our current understanding of the origin and evolution of Saturn and its atmosphere, using a multi-faceted approach that combines diverse sets of observations on volatile composition and abundances, relevant properties of the moons and rings, comparison with the other gas giant planet, Jupiter, analogies to the extrasolar giant planets, as well as pertinent theoretical models.

We present three-dimensional simulations with nested meshes of the dynamics of the gas around a Jupiter mass planet with the JUPITER and FARGOCA codes. We implemented a radiative transfer module into the JUPITER code to account for realistic heating and cooling of the gas. We focus on the circumplanetary gas flow, determining its characteristics at very high resolution ($80\%$ of Jupiter's diameter). In our nominal simulation where the temperature evolves freely by the radiative module and reaches 13000 K at the planet, a circumplanetary envelope was formed filling the entire Roche-lobe. Because of our equation of state is simplified and probably overestimates the temperature, we also performed simulations with limited maximal temperatures in the planet region (1000 K, 1500 K, and 2000 K). In these fixed temperature cases circumplanetary disks (CPDs) were formed. This suggests that the capability to form a circumplanetary disk is not simply linked to the mass of the planet and its ability to open a gap. Instead, the gas temperature at the planet's location, which depends on its accretion history, plays also fundamental role. The CPDs in the simulations are hot and cooling very slowly, they have very steep temperature and density profiles, and are strongly sub-Keplerian. Moreover, the CPDs are fed by a strong vertical influx, which shocks on the CPD surfaces creating a hot and luminous shock-front. In contrast, the pressure supported circumplanetary envelope is characterized by internal convection and almost stalled rotation.

Accurate stellar parameters are needed in numerous domains of astrophysics. The position of stars on the H-R diagram is an important indication of their structure and evolution, and it helps improve stellar models. Furthermore, the age and mass of stars hosting planets are required elements for studying exoplanetary systems. We aim at determining accurate parameters of a set of 18 bright exoplanet host and potential host stars from interferometric measurements, photometry, and stellar models. Using the VEGA/CHARA interferometer, we measured the angular diameters of 18 stars, ten of which host exoplanets. We combined them with their distances to estimate their radii. We used photometry to derive their bolometric flux and, then, their effective temperature and luminosity to place them on the H-R diagram. We then used the PARSEC models to derive their best fit ages and masses, with error bars derived from MC calculations. Our interferometric measurements lead to an average of 1.9% uncertainty on angular diameters and 3% on stellar radii. There is good agreement between measured and indirect estimations of angular diameters (from SED fitting or SB relations) for MS stars, but not as good for more evolved stars. For each star, we provide a likelihood map in the mass-age plane; typically, two distinct sets of solutions appear (an old and a young age). The errors on the ages and masses that we provide account for the metallicity uncertainties, which are often neglected by other works. From measurements of its radius and density, we also provide the mass of 55 Cnc independently of models. From the stellar masses, we provide new estimates of semi-major axes and minimum masses of exoplanets with reliable uncertainties. We also derive the radius, density, and mass of 55 Cnc e, a super-Earth that transits its stellar host. Our exoplanetary parameters reflect the known population of exoplanets.

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

Spatially resolving the inner dust cavity of the transitional disks is a key to understanding the connection between planetary formation and disk dispersal. The disk around the Herbig star HD 139614 is of particular interest since it presents a pretransitional nature with an au-sized gap, in the dust, that was spatially resolved by mid-IR interferometry. Using new NIR interferometric observations, we aim to characterize the 0.1-10~au region of the HD~139614 disk further and identify viable mechanisms for the inner disk clearing. We report the first multiwavelength radiative transfer modeling of the interferometric data acquired on HD~139614 with PIONIER, AMBER, and MIDI, complemented by Herschel/PACS photometries. We confirm a gap structure in the um-sized dust, extending from about 2.5 au to 6 au, and constrained the properties of the inner dust component: e.g., a radially increasing surface density profile, and a depletion of 10^3 relative to the outer disk. Since self-shadowing and photoevaporation appears unlikely to be responsible for the au-sized gap of HD~139614, we thus tested if dynamical clearing could be a viable mechanism using hydrodynamical simulations to predict the gaseous disk structure. Indeed, a narrow au-sized gap is expected when a single giant planet interacts with the disk. Assuming that small dust grains are well coupled to the gas, we found that a ~ 3~Mjup planet located at 4.5 au from the star could, in less than 1 Myr, reproduce most of the aspects of the dust surface density profile, while no significant depletion in gas occurred in the inner disk, in contrast to the dust. However, the dust-depleted inner disk could be explained by the expected dust filtration by the gap and the efficient dust growth/fragmentation in the inner disk regions. Our results support the hypothesis of a giant planet opening a gap and shaping the inner region of the HD~139614 disk.

For the very first time we present 3D simulations of planets embedded in stellar irradiated discs. It is well known that thermal effects could reverse the direction of planetary migration from inwards to outwards, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. When considering stellar irradiation in addition to viscous friction as a source of heating, the outer disc changes from a shadowed to a flared structure. Using a suited analytical formula it has been shown that in the flared part of the disc the migration is inwards; planets can migrate outwards only in shadowed regions of the disc, { because the radial gradient of entropy is stronger there}. In order to confirm this result numerically, we have computed the total torque acting on planets held on fixed orbits embedded in stellar irradiated 3D discs using the hydrodynamical code FARGOCA. We find qualitatively good agreement between the total torque obtained with numerical simulations and the one predicted by the analytical formula. For large masses (>20 Earth masses) we find quantitative agreement, and we obtain outwards migration regions for planets up to 60 Earth masses in the early stages of accretional discs. We find nevertheless that the agreement with the analytic formula is quite fortuitous because the formula underestimates the size of the horseshoe region; this error is compensated by imperfect estimates of other terms, most likely the cooling rate and saturation.

The migration strength and direction of embedded low-mass planets depends on the disc structure. In discs with an efficient radiative transport, the migration can be directed outwards for planets with more than 3-5 Earth masses. This is due to the entropy driven corotation torque, a process that extends the lifetimes of growing planetary embryos. We investigate the influence on the disc structure caused by a jump in the alpha parameter of the viscosity to model a dead-zone structure in the disc. We focus on M-dot discs, which have a constant net mass flux. Using the resulting disc structure, we investigate the consequences for the formation of planetesimals and determine the regions of outward migration for proto-planets. We performed numerical hydrosimulations of M-dot discs in the r-z-plane. We used the explicit/implicit hydrodynamical code FARGOCA that includes a full tensor viscosity and stellar irradiation as well as a two-temperature solver that includes radiation transport in the flux-limited diffusion approximation. Viscosity transitions inside the disc create transitions in density that stop inward migration for small planets through the so-called "planet trap" mechanism. This mechanism also works for planets down to M_P > 0.5M_E. Additionally, the viscosity transitions change the pressure gradient in the disc, which facilitates planetesimal formation via the streaming instability. However, a very steep transition in viscosity is needed to achieve in a pressure bump in the disc. The transition in viscosity facilitates planetesimal formation and can stop the migration of small-mass planets (M_P>0.5M_E), but still does not halt inward migration of smaller planets and planetesimals that are affected by gas drag. A very steep, probably unrealistic viscosity gradient is needed to trap planets of smaller masses and halt gas-drag-driven planetesimal migration at a pressure bump.

For 61 transiting hot Jupiters, the projection of the angle between the orbital plane and the stellar equator (called the spin-orbit angle) has been measured. For about half of them, a significant misalignment is detected, and retrograde planets have been observed. This challenges scenarios of the formation of hot Jupiters. In order to better constrain formation models, we relate the distribution of the real spin-orbit angle $\Psi$ to the projected one $\beta$. Then, a comparison with the observations is relevant. We analyse the geometry of the problem to link analytically the projected angle $\beta$ to the real spin-orbit angle $\Psi$. The distribution of $\Psi$ expected in various models is taken from the literature, or derived with a simplified model and Monte-Carlo simulations in the case of the disk-torquing mechanism. An easy formula to compute the probability density function (PDF) of $\beta$ knowing the PDF of $\Psi$ is provided. All models tested here look compatible with the observed distribution beyond 40 degrees, which is so far poorly constrained by only 18 observations. But only the disk-torquing mechanism can account for the excess of aligned hot Jupiters, provided that the torquing is not always efficient. This is the case if the exciting binaries have semi-major axes as large as 10000 AU. Based on comparison with the set of observations available today, scattering models and the Kozai cycle with tidal friction models can not be solely responsible for the production of all hot Jupiters. Conversely, the presently observed distribution of the spin-orbit angles is compatible with most hot Jupiters having been transported by smooth migration inside a proto-planetary disk, itself possibly torqued by a companion.

In this paper, we address the migration of small mass planets in 3D radiative disks. Indeed, migration of small planets is known to be too fast inwards in locally isothermal conditions. However, thermal effects could reverse its direction, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. This effect has been seen for masses larger than 5 Earth masses, but the minimum mass for this to happen has never been probed numerically, although it is of crucial importance for planet formation scenarios. We have extended the hydro-dynamical code FARGO to 3D, with thermal diffusion. With this code, we perform simulations of embedded planets down to 2 Earth masses. For a set of discs parameters for which outward migration has been shown in the range of $[5, 35]$ Earth masses, we find that the transition to inward migration occurs for masses in the range $[3, 5]$ Earth masses. The transition appears to be due to an unexpected phenomenon: the formation of an asymmetric cold and dense finger of gas driven by circulation and libration streamlines. We recover this phenomenon in 2D simulations where we control the cooling effects of the gas through a simple modeling of the energy equation.

The strength and direction of migration of embedded low mass planets depends on the disc's structure. It has been shown that, in discs with viscous heating and radiative transport, the migration can be directed outwards. In this paper we investigate the influence of a constant dM/dt-flux through the disc, as well as the influence of the disc's metallicity on the disc's thermodynamics. We focus on dM/dt discs, which have a net mass flux through them. Utilizing the resulting disc structure, we determine the regions of outward migration in the disc. We perform numerical hydrosimulations of dM/dt discs with viscous heating, radiative cooling and stellar irradiation in 2D in the r-z-plane. We use the explicit/implicit hydrodynamical code FARGOCA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration of embedded planets is studied by using torque formulae. For a disc of gas surface density and viscosity, we find that the discs thermal structure depends on the product Sigma_G \nu and the amount of heavy elements, while the migration of planets additionally to the mentioned quantities, depends on the amount of viscosity itself. As a result of this, the disc structure can not be approximated by simple power laws. During the lifetime of the disc, the structure of the disc changes significantly in a non-linear way in the inner parts. In the late stages of the discs evolution, outward migration is only possible if the metallicity of the disc is high. For low metallicity, planets would migrate inwards and could potentially be lost. The presented disc structures and migration maps have important consequences on the formation of planets, as they can give hints on the different formation mechanisms for different types of planets as a function of metallicity.

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.

In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in low-viscosity regime, or both. Here we explore the second way using global, three-dimensional isothermal hydrodynamical simulations with 8 levels of nested grids around the planet. In our simulations the vertical inflow from the circumstellar disk (CSD) to the CPD determines the shape of the CPD and its accretion rate. Even without prescribed viscosity Jupiter's mass-doubling time is $\sim 10^4$ years, assuming the planet at 5.2 AU and a Minimum Mass Solar Nebula. However, we show that this high accretion rate is due to resolution-dependent numerical viscosity. Furthermore, we consider the scenario of a layered CSD, viscous only in its surface layer, and an inviscid CPD. We identify two planet-accretion mechanisms that are independent of the viscosity in the CPD: (i) the polar inflow -- defined as a part of the vertical inflow with a centrifugal radius smaller than 2 Jupiter-radii and (ii) the torque exerted by the star on the CPD. In the limit of zero effective viscosity, these two mechanisms would produce an accretion rate 40 times smaller than in the simulation.

The great diversity of extrasolar planetary systems has challenged our understanding of how planets form, and how their orbits evolve as they form. Among the various processes that may account for this diversity, the gravitational interaction between planets and their parent protoplanetary disc plays a prominent role in shaping young planetary systems. Planet-disc forces are large, and the characteristic times for the evolution of planets orbital elements are much shorter than the lifetime of protoplanetary discs. The determination of such forces is challenging, because it involves many physical mechanisms and it requires a detailed knowledge of the disc structure. Yet, the intense research of the past few years, with the exploration of many new avenues, represents a very significant improvement on the state of the discipline. This chapter reviews current understanding of planet-disc interactions, and highlights their role in setting the properties and architecture of observed planetary systems.

In the Solar System, planets have a small inclination with respect to the equatorial plane of the Sun, but there is evidence that in extrasolar systems the inclination can be very high. This spin-orbit misalignment is unexpected, as planets form in a protoplanetary disc supposedly aligned with the stellar spin. Planet-planet interactions are supposed to lead to a mutual inclination, but the effects of the protoplanetary disc are still unknown. We investigate therefore planet-disc interactions for planets above 1M_Jup. We check the influence of the inclination i, eccentricity e, and mass M_p of the planet. We perform 3D numerical simulations of protoplanetary discs with embedded high-mass planets. We provide damping formulae for i and e as a function of i, e, and M_p that fit the numerical data. For highly inclined massive planets, the gap opening is reduced, and the damping of i occurs on time-scales of the order of 10^-4 deg/yr M_disc/(0.01 M_star) with the damping of e on a smaller time-scale. While the inclination of low planetary masses (<5M_Jup) is always damped, large planetary masses with large i can undergo a Kozai-cycle with the disc. These Kozai-cycles are damped in time. Eccentricity is generally damped, except for very massive planets (M_p = 5M_Jup) where eccentricity can increase for low inclinations. The dynamics tends to a final state: planets end up in midplane and can then, over time, increase their eccentricity as a result of interactions with the disc. The interactions with the disc lead to damping of i and e after a scattering event of high-mass planets. If i is sufficiently reduced, the eccentricity can be pumped up because of interactions with the disc. If the planet is scattered to high inclination, it can undergo a Kozai-cycle with the disc that makes it hard to predict the exact movement of the planet and its orbital parameters at the dispersal of the disc.

(Abridged) Circumstellar disks are believed to be the birthplace of planets and are expected to dissipate on a timescale of a few Myr. The processes responsible for the removal of the dust and gas will strongly modify the radial distribution of the dust and consequently the SED. In particular, a young planet will open a gap, resulting in an inner disk dominating the near-IR emission and an outer disk emitting mostly in the far-IR. We analyze a full set of data (including VLTI/Pionier, VLTI/Midi, and VLT/NaCo/Sam) to constrain the structure of the transition disk around TCha. We used the Mcfost radiative transfer code to simultaneously model the SED and the interferometric observations. We find that the dust responsible for the emission in excess in the near-IR must have a narrow temperature distribution with a maximum close to the silicate sublimation temperature. This translates into a narrow inner dusty disk (0.07-0.11 AU). We find that the outer disk starts at about 12 AU and is partially resolved by the Pionier, Sam, and Midi instruments. We show that the Sam closure phases, interpreted as the signature of a candidate companion, may actually trace the asymmetry generated by forward scattering by dust grains in the upper layers of the outer disk. These observations help constrain the inclination and position angle of the outer disk. The presence of matter inside the gap is difficult to assess with present-day observations. Our model suggests the outer disk contaminates the interferometric signature of any potential companion that could be responsible for the gap opening, and such a companion still has to be unambiguously detected. We stress the difficulty to observe point sources in bright massive disks, and the consequent need to account for disk asymmetries (e.g. anisotropic scattering) in model-dependent search for companions.

When a planetary tidal disk -like Saturn's rings- spreads beyond the Roche radius (inside which planetary tides prevent aggregation), satellites form and migrate away. Here, we show that most regular satellites in the solar system probably formed in this way. According to our analytical model, when the spreading is slow, a retinue of satellites appear with masses increasing with distance to the Roche radius, in excellent agreement with Saturn's, Uranus', and Neptune's satellite systems. This suggests that Uranus and Neptune used to have massive rings that disappeared to give birth to most of their regular satellites. When the spreading is fast, only one large satellite forms, as was the case for Pluto and Earth. This conceptually bridges the gap between terrestrial and giant planet systems.

We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micron HST/NICMOS data. While sub-millimeter studies suggested there is a dust-depleted cavity with r=0.35, we find scattered light as close as 0.1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks. We find two small-scaled spiral structures which asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h ~ 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5+3,-4 Mj, in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks, we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0.5. We reach 5-sigma contrasts limiting companions to planetary masses, 3-4 MJ at 1.0 and 2 MJ at 1.55 using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.

The strength and direction of migration of low mass planets depends on the disc's thermodynamics. In discs where the viscous heating is balanced by radiative transport, the migration can be directed outwards, a process which extends the lifetime of growing planetary embryos. We investigate the influence of opacity and stellar irradiation on the disc thermodynamics. Utilizing the resulting disc structure, we determine the regions of outward migration. We perform two-dimensional numerical simulations of equilibrium discs with viscous heating, radiative cooling and stellar irradiation. We use the hydrodynamical code NIRVANA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration is studied by using torque formulae. In the constant opacity case, we reproduce the analytical results of a black-body disc: the stellar irradiation dominates in the outer regions -- leading to flaring -- while the viscous heating dominates close to the star. We find that the inner edge of the disc should not be significantly puffed-up by the stellar irradiation. If the opacity depends on the local density and temperature, the structure of the disc is different, and several bumps in the aspect ratio H/r appear, due to transitions between different opacity regimes. The bumps in the disc can shield the outer disc from stellar irradiation. Stellar irradiation is an important factor for determining the disc structure and has dramatic consequences for the migration of embedded planets. Compared to discs with only viscous heating, a stellar irradiated disc features a much smaller region of outward migration for a range of planetary masses. This suggests that the region where the formation of giant planet cores takes place is smaller, which in turn might lead to a shorter growth phase.

With hundreds of exoplanets detected, it is necessary to revisit giant planets accretion models to explain their mass distribution. In particular, formation of sub-jovian planets remains unclear, given the short timescale for the runaway accretion of massive atmospheres. However, gas needs to pass through a circum-planetary disc. If the latter has a low viscosity (as expected if planets form in "dead zones"), it might act as a bottleneck for gas accretion. We investigate what the minimum accretion rate is for a planet under the limit assumption that the circum-planetary disc is totally inviscid, and the transport of angular momentum occurs solely because of the gravitational perturbations from the star. To estimate the accretion rate, we present a steady-state model of an inviscid circum-planetary disc, with vertical gas inflow and external torque from the star. Hydrodynamical simulations of a circum-planetary disc were conducted in 2D, in a planetocentric frame, with the star as an external perturber in order to measure the torque exerted by the star on the disc. The disc shows a two-armed spiral wave caused by stellar tides, propagating all the way in from the outer edge of the disc towards the planet. The stellar torque is small and corresponds to a doubling time for a Jupiter mass planet of the order of 5 Myrs. Given the limit assumptions, this is clearly a lower bound of the real accretion rate. This result shows that gas accretion onto a giant planet can be regulated by circum-planetary discs. This suggests that the diversity of masses of extra-solar planets may be the result of different viscosities in these discs.

We show that the existence of prograde equatorial satellites is consistent with a collisional tilting scenario for Uranus. In fact, if the planet was surrounded by a proto-satellite disk at the time of the tilting and a massive ring of material was temporarily placed inside the Roche radius of the planet by the collision, the proto-satellite disk would have started to precess incoherently around the equator of the planet, up to a distance greater than that of Oberon. Collisional damping would then have collapsed it into a thin equatorial disk, from which the satellites eventually formed. The fact that the orbits of the satellites are prograde requires Uranus to have had a non-negligible initial obliquity (comparable to that of Neptune) before it was finally tilted to 98 degrees.

With the uniquely high contrast within 0.1" (\Delta mag(L') = 5-6.5 magnitudes) available using Sparse Aperture Masking (SAM) with NACO at VLT, we detected asymmetry in the flux from the Herbig Fe star HD 142527 with a barycenter emission situated at a projected separation of 88+/-5 mas (12.8+/-1.5 AU at 145 pc) and flux ratios in H, K, and L' of 0.016+/-0.007, 0.012+/-0.008, 0.0086+/-0.0011 respectively (3-\sigma errors) relative to the primary star and disk. After extensive closure-phase modeling, we interpret this detection as a close-in, low-mass stellar companion with an estimated mass of ~0.1-0.4 M_Sun. HD 142527 has a complex disk structure, with an inner gap imaged in both the near and mid-IR as well as a spiral feature in the outer disk in the near-IR. This newly detected low-mass stellar companion may provide a critical explanation of the observed disk structure.

The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested in Canup (2010). However, the process driving mid-sized moons accretion from the icy debris disks has not been investigated in details. In particular, this process does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea's mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Qp<2000) would be required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites' cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process)while those still present in the rings are interpreted as today Saturn's rings' propellers and ring-moons (like Pan or Daphnis).

The discovery of multiple transiting planetary systems offers new possibilities for characterising exoplanets and understanding their formation. The Kepler-9 system contains two Saturn-mass planets, Kepler-9b and 9c. Using evolution models of gas giants that reproduce the sizes of known transiting planets and accounting for all sources of uncertainties, we show that Kepler-9b (respectively 9c) contains $45^{+17}_{-12}$\,\mearth\ (resp. $31^{+13}_{-10}$\,\mearth) of hydrogen and helium and $35^{+10}_{-15}$\,\mearth (resp. $24^{+10}_{-12}$\,\mearth) of heavy elements. More accurate constraints are obtained when comparing planets 9b and 9c: the ratio of the total mass fractions of heavy elements are $Z_{\rm b}/Z_{\rm c}=1.02\pm 0.14$, indicating that, although the masses of the planets differ, their global composition is very similar, an unexpected result for formation models. Using evolution models for super-Earths, we find that Kepler-9d must contain less than 0.1% of its mass in hydrogen and helium and predict a mostly rocky structure with a total mass between 4 and 16\,\mearth.

Hundred meter sized objects have been identified by the Cassini spacecraft in Saturn's A ring through the so-called "propeller" features they create in the ring. These moonlets should migrate, due to their gravitational interaction with the ring ; in fact, some orbital variation have been detected. The standard theory of type I migration of planets in protoplanetary disks can't be applied to the ring system, as it is pressureless. Thus, we compute the differential torque felt by a moonlet embedded in a two-dimensional disk of solid particles, with flat surface density profile, both analytically and numerically. We find that the corresponding migration rate is too small to explain the observed variations of the propeller's orbit in Saturn's A-ring. However, local density fluctuations (due to gravity wakes in the marginally gravitationally stable A-ring) may exert a stochastic torque on a moonlet. Our simulations show that this torque can be large enough to account for the observations, depending on the parameters of the rings. We find that on time scales of several years the migration of propellers is likely to be dominated by stochastic effects (while the former, non-stochastic migration dominates after ~ 10^{4-5} years). In that case, the migration rates provided by observations so far suggests that the surface density of the A ring should be of the order of 700 kg/m^2. The age of the propellers shouldn't exceed 1 to 100 million years, depending on the dominant migration regime.

We investigate the long-term and large-scale viscous evolution of dense planetary rings using a simple 1D numerical code. We use a physically realistic viscosity model derived from N-body simulations (Daisaka et al., 2001), and dependent on the disk's local properties (surface mass density, particle size, distance to the planet). Particularly, we include the effects of gravitational instabilities (wakes) that importantly enhance the disk's viscosity. We show that common estimates of the disk's spreading time-scales with constant viscosity significantly underestimate the rings' lifetime. With a realistic viscosity model, an initially narrow ring undergoes two successive evolutionary stages: (1) a transient rapid spreading when the disk is self-gravitating, with the formation of a density peak inward and an outer region marginally gravitationally stable, and with an emptying time-scale proportional to 1/M_0^2 (where M_0 is the disk's initial mass) (2) an asymptotic regime where the spreading rate continuously slows down as larger parts of the disk become not-self-gravitating due to the decrease of the surface density, until the disk becomes completely not-self-gravitating. At this point its evolution dramatically slows down, with an emptying time-scale proportional to 1/M_0, which significantly increases the disk's lifetime compared to the case with constant viscosity. We show also that the disk's width scales like t^{1/4} with the realistic viscosity model, while it scales like t^{1/2} in the case of constant viscosity, resulting in much larger evolutionary time-scales in our model. We find however that the present shape of Saturn's rings looks like a 100 million-years old disk in our simulations. Concerning Jupiter's, Uranus' and Neptune's rings that are faint today, it is not likely that they were much more massive in the past and lost most of their mass due to viscous spreading alone.

Recent observations of exoplanets by direct imaging, reveal that giant planets orbit at a few dozens to more than a hundred of AU from their central star. The question of the origin of these planets challenges the standard theories of planet formation. We propose a new way of obtaining such far planets, by outward migration of a pair of planets formed in the 10 AU region. Two giant planets in mean motion resonance in a common gap in the protoplanetary disk migrate outwards, if the inner one is significantly more massive than the outer one. Using hydrodynamical simulations, we show that their semi major axes can increase by almost one order of magnitude. In a flared disk, the pair of planets should reach an asymptotic radius. This mechanism could account for the presence of Fomalhaut b ; then, a second, more massive planet, should be orbiting Fomalhaut at about 75 AU.

We study the torque on low-mass planets embedded in protoplanetary discs in the two-dimensional approximation, incorporating non-isothermal effects. We couple linear estimates of the Lindblad (or wave) torque to a simple, but non-linear, model of adiabatic corotation torques (or horseshoe drag), resulting in a simple formula that governs Type I migration in non-isothermal discs. This formula should apply in optically thick regions of the disc, where viscous and thermal diffusion act to keep the horseshoe drag unsaturated. We check this formula against numerical hydrodynamical simulations, using three independent numerical methods, and find good agreement.

Numerical simulations of planets embedded in protoplanetary gaseous discs are a precious tool for studying the planetary migration ; however, some approximations have to be made. Most often, the selfgravity of the gas is neglected. In that case, it is not clear in the literature how the material inside the Roche lobe of the planet should be taken into account. Here, we want to address this issue by studying the influence of various methods so far used by different authors on the migration rate. We performed high-resolution numerical simulations of giant planets embedded in discs. We compared the migration rates with and without gas selfgravity, testing various ways of taking the circum-planetary disc (CPD) into account. Different methods lead to significantly different migration rates. Adding the mass of the CPD to the perturbing mass of the planet accelerates the migration. Excluding a part of the Hill sphere is a very touchy parameter that may lead to an artificial suppression of the type III, runaway migration. In fact, the CPD is smaller than the Hill sphere. We recommend excluding no more than a 0.6 Hill radius and using a smooth filter. Alternatively, the CPD can be given the acceleration felt by the planet from the rest of the protoplanetary disc. The gas inside the Roche lobe of the planet should be very carefully taken into account in numerical simulations without any selfgravity of the gas. The entire Hill sphere should not be excluded. The method used should be explicitly given. However, no method is equivalent to computing the full selfgravity of the gas.

The Minimum Mass Solar Nebula (MMSN) is a protoplanetary disk that contains the minimum amount of solids necessary to build the planets of the Solar System. Assuming that the giant planets formed in the compact configuration they have at the beginning of the "Nice model", Desch (2007) built a new MMSN. He finds a decretion disk, about ten times denser than the well-known Hayashi MMSN. The disk profile is almost stationary for about ten million years. However, a planet in a protoplanetary disk migrates. In a massive, long-lived disk, this question has to be addressed. With numerical simulations, we show that the four giant planets of the Solar System could not survive in this disk. In particular, Jupiter enters the type III, runaway regime, and falls into the Sun like a stone. Known planet-planet interaction mechanisms to prevent migration, fail in this nebula, in contrast to the Hayashi MMSN. Planetary migration constrains the construction of a MMSN. We show how this should be done self-consistently.

In this review, three major changes in our understanding of the early history of the Solar System are presented. 1) Early differentiation: A few recent results support the idea that protoplanet formation and differentiation occurred partly simultaneously than CAI formation. First, some iron meteorites, eucrites, and angrites older than the chondrules or even than the CAI have been found. Second, iron meteorites could be debris of early disrupted differentiated planetesimals, scattered from the terrestrial planet region to the Main Belt. Finally, chondrules contain fragments of planetesimal material. 2) Earth and Moon: An equilibration mechanism explains the identical Oxygen isotopic composition of the Earth and the Moon. In addition, it has been shown that the Earth and the Moon mantles have the same 182^W anomaly, in contrast to what was believed before. Consequently, the Moon forming impact should have occurred after the extinction of the 182Hf radioactivity, about 60 Myr after Solar System formation. This new datation is in agreement with new N-body numerical simulations of the last phase of terrestrial planets formation, in which giant impacts occur during about 100 Myr. 3) Giant planets and Nice model: The migration of the giant planets in the protoplanetary disc can be prevented if the planets are in resonance, close to each other. In the ``Nice model'', the 4 outer planets of the Solar System were in a compact configuration after the dissipation of gaseous disc. A few hundred million years later, a global instability drives the planets on their present orbits, producing the Late Heavy Bombardment. In this frame, a lot of characteristics of our Solar System can be explained.

The formation of resonant planets pairs in exoplanetary systems involves planetary migration inside the protoplanetary disc : an inwards migrating outer planet captures in Mean Motion Resonance an inner planet. During the migration of the resonant pair of planets, the eccentricities are expected to rise excessively, if no damping mechanism is applied on the inner planet. We express the required damping action to match the observations, and we show that the inner disk can play this role. This result applies for instance to the system GJ 876 : we reproduce the observed orbital elements through a fully hydrodynamical simulation of the evolution of the resonant planets.

In isothermal disks the migration of protoplanets is directed inward. For small planetary masses the standard type-I migration rates are so fast that this may result in an unrealistic loss of planets into the stars. We investigate the planet-disk interaction in non-isothermal disks and analyze the magnitude and direction of migration for an extended range of planet masses. We have performed detailed two-dimensional numerical simulations of embedded planets including heating/cooling effects as well as radiative diffusion for realistic opacities. In radiative disks, small planets with M_planet < 50 M_Earth do migrate outward with a rate comparable to absolute magnitude of standard type-I migration. For larger masses the migration is inward and approaches the isothermal, type-II migration rate. Our findings are particularly important for the first growth phase of planets and ease the problem of too rapid inward type-I migration.

The formation of resonant pairs of planets in exoplanetary systems involves planetary migration in the protoplanetary disc. After a resonant capture, the subsequent migration in this configuration leads to a large increase of planetary eccentricities if no damping mechanism is applied. This has led to the conclusion that the migration of resonant planetary systems cannot occur over large radial distances and has to be terminated sufficiently rapidly through disc dissipation. In this study, we investigate whether the presence of an inner disc might supply an eccentricity damping of the inner planet, and if this effect could explain the observed eccentricities in some systems. To investigate the influence of an inner disc, we first compute hydrodynamic simulations of giant planets orbiting with a given eccentricity around an inner gas disc, and measure the effect of the latter on the planetary orbital parameters. We then perform detailed long term calculations of the GJ 876 system. We also run N-body simulations with artificial forces on the planets mimicking the effects of the inner and outer discs. We find that the influence of the inner disc can not be neglected, and that it might be responsible for the observed eccentricities. In particular, we reproduce quite well the orbital parameters of a few systems engaged in 2:1 mean motion resonances: GJ 876, HD 73526, HD 82943 and HD 128311. Finally, we derive analytically the effect that the inner disc should have on the inner planet to reach a specific orbital configuration with a given damping effect of the outer disc on the outer planet.

A well-known bottleneck for the core-accretion model of giant-planet formation is the loss of the cores into the star by Type-I migration, due to the tidal interactions with the gas disk. It has been shown that a steep surface-density gradient in the disk, such as the one expected at the boundary between an active and a dead zone, acts as a planet trap and prevents isolated cores from migrating down to the central star. We study the relevance of the planet trap concept for the accretion and evolution of systems of multiple planetary embryos/cores. We performed hydrodynamical simulations of the evolution of systems of multiple massive objects in the vicinity of a planet trap. The planetary embryos evolve in 3 dimensions, whereas the disk is modeled with a 2D grid. Synthetic forces are applied onto the embryos to mimic the damping effect that the disk has on their inclinations. Systems with two embryos tend to acquire stable, separated and non-migrating orbits, with the more massive embryo placed at the planet trap and the lighter one farther out in the disk. Systems of multiple embryos are intrinsically unstable. Consequently, a long phase of mutual scattering can lead to accreting collisions among embryos; some embryos are injected into the inner part of the disk, where they can be evacuated into the star by Type I migration. The system can resume a stable, non-migrating configuration only when the number of surviving embryos decreases to a small value (~2-4). This can explain the limited number of giant planets in our solar system. These results should apply in general to any case in which the Type-I migration of the inner embryo is prevented by some mechanism, and not solely to the planet trap scenario.

We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonance (i.e. each planet is in resonance with its neighbor). In total we found 6 such configurations. For the gas disk parameters found in MC07, these configurations are characterized by a negligible migration rate. After the disappearance of the gas, and in absence of planetesimals, only two of these six configurations (the least compact ones) are stable for a time of hundreds of millions of years or more. The others become unstable on a timescale of a few My. Our preliminary simulations show that, when a planetesimal disk is added beyond the orbit of the outermost planet, the planets can evolve from the most stable of these configurations to their current orbits in a fashion qualitatively similar to that described in Tsiganis et al. (2005).

We study the possibility that the mutual interactions between Jupiter and Saturn prevented Type II migration from driving these planets much closer to the Sun. Our work extends previous results by Masset and Snellgrove (2001), by exploring a wider set of initial conditions and disk parameters, and by using a new hydrodynamical code that properly describes for the global viscous evolution of the disk. Initially both planets migrate towards the Sun, and Saturn's migration tends to be faster. As a consequence, they eventually end up locked in a mean motion resonance. If this happens in the 2:3 resonance, the resonant motion is particularly stable, and the gaps opened by the planets in the disk may overlap. This causes a drastic change in the torque balance for the two planets, which substantially slows down the planets' inward migration. If the gap overlap is substantial, planet migration may even be stopped or reversed. As the widths of the gaps depend on disk viscosity and scale height, this mechanism is particularly efficient in low viscosity, cool disks. We discuss the compatibility of our results with the initial conditions adopted in Tsiganis et al. (2005) and Gomes et al. (2005) to explain the current orbital architecture of the giant planets and the origin of the Late Heavy Bombardment of the Moon.

We study the effect of a Jovian planet on the gas distribution of a protoplanetary disc, using a new numerical scheme that allows us to take into consideration the global evolution of the disc, down to an arbitrarily small inner physical radius. We find that Jovian planets do not open cavities in the inner part of the disc (i.e. interior to their orbits) unless (a) the inner physical edge of the disc is close to the planet's location or (b) the planet is much more massive than the disc. In all other cases the planet simply opens a gap in the gas density distribution, whose global profile is essentially unchanged relative to the one that it would have if the planet were absent. We recognize, though, that the dust distribution can be significantly different from the gas distribution and that dust cavities might be opened in some situations, even if the gas is still present in the inner part of the disc. Concerning the migration of the planet, we find that classical Type-II migration (with speed proportional to the viscosity of the disc) occurs only if the gap opened by the planet is deep and clean. If there is still a significant amount of gas in the gap, the migration of the planet is generally slower than the theoretical Type-II migration rate. In some situations, migration can be stopped or even reversed. We develop a simple model that reproduces satisfactorily the migration rate observed in the simulations, for a wide range of disc viscosities and planet masses and locations relative to the inner disc edge. Our results are relevant for extra-solar planetary systems, as they explain (a) why some hot Jupiters did not migrate all the way down to their parent stars and (b) why the outermost of a pair of resonant planets is typically the most massive one.

Numerical simulations of planet-disk interactions are usually performed with hydro-codes that -- because they consider only an annulus of the disk, over a 2D grid -- can not take into account the global evolution of the disk. However, the latter governs planetary migration of type II, so that the accuracy of the planetary evolution can be questioned. To develop an algorithm that models the local planet-disk interactions together with the global viscous evolution of the disk, we surround the usual 2D grid with a 1D grid ranging over the real extension of the disk. The 1D and 2D grids are coupled at their common boundaries via ghost rings, paying particular attention to the fluxes at the interface, especially the flux of angular momentum carried by waves. The computation is done in the frame centered on the center of mass to ensure angular momentum conservation. The global evolution of the disk and the local planet-disk interactions are both well described and the feedback of one on the other can be studied with this algorithm, for a negligible additional computing cost with respect to usual algorithms.

[Abridged] The tidal torque exerted by a protoplanetary disk with power law surface density and temperature profiles onto an embedded protoplanetary embryo is generally a negative quantity that leads to the embryo inwards migration. Here we investigate how the tidal torque balance is affected at a disk surface density radial jump. The jump has two consequences : - it affects the differential Lindblad torque. In particular if the disk is merely empty on the inner side, the differential Lindblad torque almost amounts to the large negative outer Lindblad torque. - It affects the corotation torque, which is a quantity very sensitive to the local gradient of the disk surface density. In particular if the disk is depleted on the inside and if the jump occurs radially over a few pressure scale-heights, the corotation torque is a positive quantity that is much larger than in a power-law disk. We show by means of customized numerical simulations of low mass planets embedded in protoplanetary nebulae with a surface density jump that the second effect is dominant, that is that the corotation torque largely dominates the differential Lindblad torque on the edge of a central depletion. As a consequence the type I migration of low mass objects reaching the jump should be halted, and all these objects should be trapped there.

Although it is well known that a massive planet opens a gap in a proto-planetary gaseous disk, there is no analytic description of the surface density profile in and near the gap. The simplest approach, which is based upon the balance between the torques due to the viscosity and the gravity of the planet and assumes local damping, leads to gaps with overestimated width, especially at low viscosity. Here, we take into account the fraction of the gravity torque that is evacuated by pressure supported waves. With a novel approach, which consists of following the fluid elements along their trajectories, we show that the flux of angular momentum carried by the waves corresponds to a pressure torque. The equilibrium profile of the disk is then set by the balance between gravity, viscous and pressure torques. We check that this balance is satisfied in numerical simulations, with a planet on a fixed circular orbit. We then use a reference numerical simulation to get an ansatz for the pressure torque, that yields gap profiles for any value of the disk viscosity, pressure scale height and planet to primary mass ratio. Those are in good agreement with profiles obtained in numerical simulations over a wide range of parameters. Finally, we provide a gap opening criterion that simultaneously involves the planet mass, the disk viscosity and the aspect ratio.