Tiny meteoroids entering the Earth's atmosphere and inducing meteor showers have long been thought to originate partly from cometary dust. Together with other dust particles, they form a huge cloud around the Sun, the zodiacal cloud. From our previous studies of the zodiacal light, as well as other independent methods (dynamical studies, infrared observations, data related to Earth's environment), it is now established that a significant fraction of dust particles entering the Earth's atmosphere comes from Jupiter-family comets (JFCs). This paper relies on our understanding of key properties of the zodiacal cloud and of comet 67P/Churyumov-Gerasimenko, extensively studied by the Rosetta mission to a JFC. The interpretation, through numerical and experimental simulations of zodiacal light local polarimetric phase curves, has recently allowed us to establish that interplanetary dust is rich in absorbing organics and consists of fluffy particles. The ground-truth provided by Rosetta presently establishes that the cometary dust particles are rich in organic compounds and consist of quite fluffy and irregular aggregates. Our aims are as follows: (1) to make links, back in time, between peculiar micrometeorites, tiny meteoroids, interplanetary dust particles, cometary dust particles, and the early evolution of the Solar System, and (2) to show how detailed studies of such meteoroids and of cometary dust particles can improve the interpretation of observations of dust in protoplanetary and debris disks. Future modeling of dust in such disks should favor irregular porous particles instead of more conventional compact spherical particles.

A key challenge for protoplanetary disks and planet formation models is to be able to make a reliable connection between observed structures in the disks emission, like bright and dark rings or asymmetries, and the supposed existence of planets triggering these structures. The observation of N dark rings of emission is often interpreted as evidence for the presence of N planets which clear dust gaps around their orbit and form dust-trapping pressure maxima in the disk. The vast majority of the models that studied the impact of planets on the dynamics of dust and gas in a protoplanetary disk assumed planets on fixed orbits. Here we go a different route and examine how the large-scale inward migration of a single planet structures the dust content of a massive disk. In many circumstances, the migration of a partial gap-opening planet with a mass comparable to Saturn is found to run away intermittently. By means of 2D gas and dust hydrodynamical simulations, we show that intermittent runaway migration can form multiple dust rings and gaps across the disk. Each time migration slows down, a pressure maximum forms beyond the planet gap that traps the large dust. Post-processing of our simulations results with 3D dust radiative transfer calculations confirms that intermittent runaway migration can lead to the formation of multiple sets of bright and dark rings of continuum emission in the (sub)millimeter beyond the planet location.

Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets. Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters. Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context. Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates.

The dissipation of the kinetic energy of wave-like tidal flows within the convective envelope of low-mass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of short-period exoplanetary systems. Although low-mass stars are magnetically active objects, the question of how the star's magnetic field impacts large-scale tidal flows and the excitation, propagation and dissipation of tidal waves still remains open. Our goal is to investigate the impact of stellar magnetism on the forcing of tidal waves, and their propagation and dissipation in the convective envelope of low-mass stars as they evolve. We have estimated the amplitude of the magnetic contribution to the forcing and dissipation of tidally induced magneto-inertial waves throughout the structural and rotational evolution of low-mass stars (from M to F-type). For this purpose, we have used detailed grids of rotating stellar models computed with the stellar evolution code STAREVOL. The amplitude of dynamo-generated magnetic fields is estimated via physical scaling laws at the base and the top of the convective envelope. We find that the large-scale magnetic field of the star has little influence on the excitation of tidal waves in the case of nearly-circular orbits and coplanar hot-Jupiter planetary systems, but that it has a major impact on the way waves are dissipated. Our results therefore indicate that a full magneto-hydrodynamical treatment of the propagation and dissipation of tidal waves is needed to properly assess the impact of star-planet tidal interactions throughout the evolutionary history of low-mass stars hosting short-period massive planets.

While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features can be well explained with the two planets in formation, providing strong evidence that previously proposed theories of planet-disk interaction are in action, including resonant migration, particle trapping, size segregation, and filtration. Our simulations suggest that the two planets are likely in 2:1 mean motion resonance and can remain dynamically stable over million-year timescales. The growth of the planets at $10^{-8}-10^{-7}~M_{\rm Jup}~{\rm yr}^{-1}$, rates comparable to the estimates from H$\alpha$ observations, does not destabilize the resonant configuration. Large grains are filtered at the gap edge and only small, (sub-)$\mu$m grains can flow to the circumplanetary disks and the inner circumstellar disk. With the sub-millimeter continuum ring observed outward of the two directly imaged planets, PDS 70 provides the first observational evidence of particle filtration by gap-opening planets. The observed sub-millimeter continuum emission at the vicinity of the planets can be reproduced when (sub-)$\mu$m grains survive over multiple circumplanetary disk gas viscous timescales and accumulate therein. One such possibility is if (sub-)$\mu$m grains grow in size and remain trapped in pressure bumps, similar to what we find happening in circumstellar disks. We discuss potential implications to planet formation in the solar system and mature extrasolar planetary systems.

We report the analysis, conducted as part of the MaTYSSE programme, of a spectropolarimetric monitoring of the ~0.8 Myr, ~1.4 MSun disc-less weak-line T Tauri star V410 Tau with the ESPaDOnS instrument at the Canada-France-Hawaii Telescope and NARVAL at the T\'elescope Bernard Lyot, between 2008 and 2016. With Zeeman-Doppler Imaging, we reconstruct the surface brightness and magnetic field of V410 Tau, and show that the star is heavily spotted and possesses a ~550 G relatively toroidal magnetic field. We find that V410 Tau features a weak level of surface differential rotation between the equator and pole ~5 times weaker than the solar differential rotation. The spectropolarimetric data exhibit intrinsic variability, beyond differential rotation, which points towards a dynamo-generated field rather than a fossil field. Long-term variations in the photometric data suggest that spots appear at increasing latitudes over the span of our dataset, implying that, if V410 Tau has a magnetic cycle, it would have a period of more than 8 years. Having derived raw radial velocities (RVs) from our spectra, we filter out the stellar activity jitter, modeled either from our Doppler maps or using Gaussian Process Regression. Thus filtered, our RVs exclude the presence of a hot Jupiter-mass companion below ~0.1 au, which is suggestive that hot Jupiter formation may be inhibited by the early depletion of the circumstellar disc, which for V410 Tau may have been caused by the close (few tens of au) M dwarf stellar companion.

Chondrules are silicate spheroids found in meteorites, serving as important fossil records of the early solar system. In order to form chondrules, chondrule precursors must be heated to temperatures much higher than the typical conditions in the current asteroid belt. One proposed mechanism for chondrule heating is the passage through bow shocks of highly eccentric planetesimals in the protoplanetary disk in the early solar system. However, it is difficult for planetesimals to gain and maintain such high eccentricities. In this paper, we present a new scenario in which planetesimals in the asteroid belt region are excited to high eccentricities by the Jovian sweeping secular resonance in a depleting disk, leading to efficient formation of chondrules. We study the orbital evolution of planetesimals in the disk using semi-analytic models and numerical simulations. We investigate the dependence of eccentricity excitation on the planetesimal's size as well as the physical environment, and calculate the probability for chondrule formation. We find that 50 - 2000 km planetesimals can obtain eccentricities larger than 0.6 and cause effective chondrule heating. Most chondrules form in high velocity shocks, in low density gas, and in the inner disk. The fraction of chondrule precursors which become chondrules is about 4 - 9 % between 1.5 - 3 AU. Our model implies that the disk depletion timescale is $\tau_\mathrm{dep}\approx 1~\mathrm{Myr}$, comparable to the age spread of chondrules; and that Jupiter formed before chondrules, no more than 0.7 Myr after the formation of the CAIs.

The Solar System gas giants are each surrounded by many moons, with at least 50 prograde satellites thought to have formed from circumplanetary material. Just like the Sun is not the only star surrounded by planets, extrasolar gas giants are likely surrounded by satellite systems. Here, we report on ALMA observations of four <40 Myr old stars with directly-imaged companions: PZ Tel, AB Pic, 51 Eri, and $\kappa$ And. Continuum emission at 1.3 mm is undetected for any of the systems. Since these are directly-imaged companions, there is knowledge of their temperatures, masses and locations. These allow for upper limits on the amount of circumplanetary dust to be derived from detailed radiative transfer models. These protolunar disc models consider two disc sizes: 0.4 and 0.04 times the exoplanet's Hill radius. The former is representative of hydrodynamic simulations of circumplanetary discs while the latter a case with significant radial drift of solids. The more compact case is also motivated by the semi-major axis of Callisto, enclosing Jupiter's Galilean satellites. All upper limits fall below the expected amount of dust required to explain regular satellite systems ($10^{-4}$ times the mass of their central planet). Upper limits are compared with viscous evolution and debris disc models. Our analysis suggests that the non detections can be interpreted as evidence of dust growth beyond metre sizes to form moonetesimals in timescales <10 Myr. This sample increases by 50% the number of ALMA non-detections of young companions available in the literature.

AB Aur is a Herbig Ae star hosting a well-known transitional disk. Because of its proximity and low inclination angle, it is an excellent object to study planet formation. Our goal is to investigate the chemistry and dynamics of the molecular gas component in the AB Aur disk, and its relation with the prominent horseshoe shape observed in continuum mm emission. We used the NOEMA interferometer to map with high angular resolution the J = 3-2 lines of HCO+ and HCN. By combining both, we can gain insight into the AB Aur disk structure. Chemical segregation is observed in the AB Aur disk: HCO+ shows intense emission toward the star position, at least one bright molecular bridge within the dust cavity, and ring-like emission at larger radii, while HCN is only detected in an annular ring that is coincident with the dust ring and presents an intense peak close to the dust trap. We use HCO+ to investigate the gas dynamics inside the cavity. The observed bright HCO+ bridge connects the compact central source with the outer dusty ring. This bridge can be interpreted as an accretion flow from the outer ring to the inner disk/jet system proving gas accretion through the cavity.

Rings and radial gaps are ubiquitous in protoplanetary disks, yet their possible connection to planet formation is currently subject to intense debates. In principle, giant planet formation leads to wide gaps which separate the gas and dust mass reservoir in the outer disk, while lower mass planets lead to shallow gaps which are manifested mainly on the dust component. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the star HD169142, host to a prominent disk with deep wide gaps that sever the disk into inner and outer regions. The new ALMA high resolution images allow for the outer ring to be resolved as three narrow rings. The HD169142 disk thus hosts both the wide gaps trait of transition disks and a narrow ring system similar to those observed in the TW Hya and HL Tau systems. The mass reservoir beyond a deep gap can thus host ring systems. The observed rings are narrow in radial extent (width/radius of 1.5/57.3, 1.8/64.2 and 3.4/76.0, in au) and have asymmetric mutual separations: the first and middle ring are separated by 7 au while the middle and outermost ring are distanced by ~12 au. Using hydrodynamical modeling we found that a simple explanation, involving a single migrating low mass planet (10 M$_\oplus$), entirely accounts for such an apparently complex phenomenon. Inward migration of the planet naturally explains the ring's asymmetric mutual separation. The isolation of HD169142's outer rings thus allows a proof of concept to interpret the detailed architecture of the outer region of protoplanetary disks with low mass planet formation of mini-Neptune's size, i.e. as in the protosolar nebula.

Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario.

The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques.

The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains.

The energy dissipation of wave-like tidal flows in the convective envelope of low-mass stars is one of the key physical mechanisms that shape the orbital and rotational dynamics of short-period planetary systems. Tidal flows, and the excitation, propagation, and dissipation of tidally-induced inertial waves can be modified by stellar magnetic fields (e.g., Wei 2016, 2018, Lin and Ogilvie 2018). It is thus important to assess for which stars, at which location of their internal structure, and at which phase of their evolution, one needs to take into account the effects of magnetic fields on tidal waves. Using scaling laws that provide the amplitude of dynamo-generated magnetic fields along the rotational evolution of these stars (e.g., Christensen et al. 2009, Brun et al. 2015), combined with detailed grids of stellar rotation models (e.g., Amard et al. 2016), we examine the influence of a magnetic field on tidal forcing and dissipation near the tachocline of solar-like stars. We show that full consideration of magnetic fields is required to compute tidal dissipation, but not necessarily for tidal forcing.

When a planet becomes massive enough, it gradually carves a partial gap around its orbit in the protoplanetary disk. A pressure maximum can be formed outside the gap where solids that are loosely coupled to the gas, typically in the pebble size range, can be trapped. The minimum planet mass for building such a trap, which is called the pebble isolation mass (PIM), is important for two reasons: it marks the end of planetary growth by pebble accretion, and the trapped dust forms a ring that may be observed with millimetre observations. We study the effect of disk turbulence on the pebble isolation mass and find its dependence on the gas turbulent viscosity, aspect ratio, and particles Stokes number. By means of 2D gas hydrodynamical simulations, we found the minimum planet mass to form a radial pressure maximum beyond the orbit of the planet, which is the necessary condition to trap pebbles. We then carried out 2D gas plus dust hydrodynamical simulations to examine how dust turbulent diffusion impacts particles trapping at the pressure maximum. We finally provide a semi-analytical calculation of the PIM based on comparing the radial drift velocity of solids and the root mean square turbulent velocity fluctuations around the pressure maximum. From our results of gas simulations, we provide an expression for the PIM versus disk aspect ratio and turbulent viscosity. Our gas plus dust simulations show that the effective PIM can be nearly an order of magnitude larger in high-viscosity disks because turbulence diffuse particles out of the pressure maximum. This is quantified by our semi-analytical calculation, which gives an explicit dependence of the PIM with Stokes number of particles. We conclude than disk turbulence can significantly alter the PIM, depending on the level of turbulence in regions of planet formation.

Star-planet tidal interactions play a significant role in the dynamical evolution of close-in planetary systems. We investigate the propagation and dissipation of tidal inertial waves in a stellar/planetary convective region. We take into account a latitudinal differential rotation for the background flow, similar to what is observed in the envelope of low-mass stars like the Sun. Previous works have shown that differential rotation significantly alters the propagation and dissipation properties of inertial waves. In particular, when the Doppler-shifted tidal frequency vanishes in the fluid, a critical layer forms where tidal dissipation can be greatly enhanced. Our present work develops a local analytic model to better understand the propagation and dissipation properties of tidally forced inertial waves at critical layers.

One serious challenge for planet formation is the rapid inward drift of pebble-sized dust particles in protoplanetary disks. Dust trapping at local maxima in the disk gas pressure has received much theoretical attention but still lacks observational support. The cold dust emission in the AB Aur disk forms an asymmetric ring at a radius of about 120 au, which is suggestive of dust trapping in a gas vortex. We present high spatial resolution (0".58x0".78 ~ 80x110 au) NOEMA observations of the 1.12 mm and 2.22 mm dust continuum emission from the AB Aur disk. Significant azimuthal variations of the flux ratio at both wavelengths indicate a size segregation of the large dust particles along the ring. Our continuum images also show that the intensity variations along the ring are smaller at 2.22 mm than at 1.12 mm, contrary to what dust trapping models with a gas vortex have predicted. Our two-fluid (gas+dust) hydrodynamical simulations demonstrate that this feature is well explained if the gas vortex has started to decay due to turbulent diffusion, and dust particles are thus losing the azimuthal trapping on different timescales depending on their size. The comparison between our observations and simulations allows us to constrain the size distribution and the total mass of solid particles in the ring, which we find to be of the order of 30 Earth masses, enough to form future rocky planets.

We report the results of an extended spectropolarimetric and photometric monitoring of the weak-line T Tauri star TAP 26, carried out within the MaTYSSE programme with the ESPaDOnS spectropolarimeter at the 3.6 m Canada-France-Hawaii Telescope. Applying Zeeman-Doppler Imaging to our observations, concentrating in 2015 November and 2016 January and spanning 72 d in total, 16 d in 2015 November and 13 d in 2016 January, we reconstruct surface brightness and magnetic field maps for both epochs and demonstrate that both distributions exhibit temporal evolution not explained by differential rotation alone. We report the detection of a hot Jupiter (hJ) around TAP 26 using three different methods, two using Zeeman-Doppler Imaging (ZDI) and one Gaussian-Process Regression (GPR), with a false-alarm probability smaller than 6.10^-4. However, as a result of the aliasing related to the observing window, the orbital period cannot be uniquely determined; the orbital period with highest likelihood is 10.79 +/- 0.14 d followed by 8.99 +/- 0.09 d. Assuming the most likely period, and that the planet orbits in the stellar equatorial plane, we obtain that the planet has a minimum mass M.sin(i) of 1.66 +/- 0.31 M_Jup and orbits at 0.0968 +/- 0.0032 au from its host star. This new detection suggests that disc type II migration is efficient at generating newborn hJs, and that hJs may be more frequent around young T Tauri stars than around mature stars (or that the MaTYSSE sample is biased towards hJ-hosting stars).

Tidal interactions in close star-planet or binary star systems may excite inertial waves (their restoring force is the Coriolis force) in the convective region of the stars. The dissipation of these waves plays a prominent role in the long-term orbital and rotational evolution of the bodies involved. If the primary star rotates as a solid body, inertial waves have a Doppler-shifted frequency restricted to the range $[-2\Omega, 2\Omega]$ ($\Omega$ being the angular velocity of the star), and they can propagate in the entire convective region. However, turbulent convection can sustain differential rotation with an equatorial acceleration (as in the Sun) or deceleration that modifies the frequency range and propagation domain of inertial waves and allows corotation resonances for non-axisymmetric oscillations. In this work, we perform numerical simulations of tidally excited inertial waves in a differentially rotating convective envelope with a conical (or latitudinal) rotation profile. The tidal forcing that we adopt contains spherical harmonics that correspond to the case of a circular and coplanar orbit. We study the viscous dissipation of the waves as a function of tidal frequency for various stellar masses and differential rotation parameters, as well as its dependence on the turbulent viscosity coefficient. We compare our results with previous studies assuming solid-body rotation and point out the potential key role of corotation resonances in the dynamical evolution of close-in star-planet or binary systems.

Context: Quantifying the gas content inside the dust gaps of transition disks is important to establish their origin. Aims: We seek to constrain the surface density of warm gas in the disk of HD 139614, a Herbig Ae star with a transition disk exhibiting a dust gap from 2.3 to 6 AU. Methods: We have obtained ESO/VLT CRIRES high-resolution spectra of CO ro-vibrational emission. We derived constraints on the disk's structure by modeling the line-profiles, the spectroastrometric signal, and the rotational diagrams using flat Keplerian disk models. Results: We detected v=1-0 12CO, 2-1 12CO, 1-0 13CO, 1-0 C18O, and 1-0 C17O ro-vibrational lines. 12CO v=1-0 lines have an average width of 14 km/s, Tgas of 450 K and an emitting region from 1 to 15 AU. 13CO and C18O lines are on average 70 and 100 K colder, 1 and 4 km/s narrower, and are dominated by emission at R>6 AU. The 12CO v=1-0 line-profile indicates that if there is a gap in the gas it must be narrower than 2 AU. We find that a drop in the gas surface density (delta_gas) at R<5-6 AU is required to simultaneously reproduce the line-profiles and rotational diagrams of the three CO isotopologs. Delta_gas can range from 10^-2 to 10^-4 depending on the gas-to-dust ratio of the outer disk. We find that at 1<R<6 AU the gas surface density profile is flat or increases with radius. We derive a gas column density at 1<R<6 AU of NH=3x10^19 - 10^21 cm^-2. We find a 5sigma upper limit on NCO at R<1 AU of 5x10^15 cm^-2 (NH<5x10^19 cm^-2). Conclusions: The dust gap in the disk of HD 139614 has gas. The gas surface density in the disk at R<6 AU is significantly lower than the surface density expected from HD 139614's accretion rate assuming a viscous alpha-disk model. The gas density drop, the non-negative density gradient of the gas inside 6 AU, and the absence of a wide (>2 AU) gas gap suggest the presence of an embedded <2 MJ planet at around 4 AU.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Hot Jupiters are giant Jupiter-like exoplanets that orbit 100x closer to their host stars than Jupiter does to the Sun. These planets presumably form in the outer part of the primordial disc from which both the central star and surrounding planets are born, then migrate inwards and yet avoid falling into their host star. It is however unclear whether this occurs early in the lives of hot Jupiters, when still embedded within protoplanetary discs, or later, once multiple planets are formed and interact. Although numerous hot Jupiters were detected around mature Sun-like stars, their existence has not yet been firmly demonstrated for young stars, whose magnetic activity is so intense that it overshadows the radial velocity signal that close-in giant planets can induce. Here we show that hot Jupiters around young stars can be revealed from extended sets of high-resolution spectra. Once filtered-out from the activity, radial velocities of V830 Tau derived from new data collected in late 2015 exhibit a sine wave of period 4.93 d and semi-amplitude 75 m/ s, detected with a false alarm probability <0.03%. We find that this signal is fully unrelated to the 2.741-d rotation period of V830 Tau and we attribute it to the presence of a 0.77 Jupiter mass planet orbiting at a distance of 0.057 au from the host star. Our result demonstrates that hot Jupiters can migrate inwards in <2 Myr, most likely as a result of planet-disc interactions, and thus yields strong support to the theory of giant planet migration in gaseous protoplanetary discs.

The gravito-inertial waves propagating over a shellular baroclinic flow inside a rotating spherical shell are analysed using the Boussinesq approximation. The wave properties are examined by computing paths of characteristics in the non-dissipative limit, and by solving the full dissipative eigenvalue problem using a high-resolution spectral method. Gravito-inertial waves are found to obey a mixed-type second-order operator and to be often focused around short-period attractors of characteristics or trapped in a wedge formed by turning surfaces and boundaries. We also find eigenmodes that show a weak dependence with respect to viscosity and heat diffusion just like truly regular modes. Some axisymmetric modes are found unstable and likely destabilized by baroclinic instabilities. Similarly, some non-axisymmetric modes that meet a critical layer (or corotation resonance) can turn unstable at sufficiently low diffusivities. In all cases, the instability is driven by the differential rotation. For many modes of the spectrum, neat power laws are found for the dependence of the damping rates with diffusion coefficients, but the theoretical explanation for the exponent values remains elusive in general. The eigenvalue spectrum turns out to be very rich and complex, which lets us suppose an even richer and more complex spectrum for rotating stars or planets that own a differential rotation driven by baroclinicity.

The growing body of observational data on extrasolar planets and protoplanetary disks has stimulated intense research on planet formation and evolution in the past few years. The extremely diverse, sometimes unexpected physical and orbital characteristics of exoplanets lead to frequent updates on the mainstream scenarios for planet formation and evolution, but also to the exploration of alternative avenues. The aim of this review is to bring together classical pictures and new ideas on the formation, orbital and internal evolutions of planets, highlighting the key role of the protoplanetary disk in the various parts of the theory. We begin by briefly reviewing the conventional mechanism of core accretion by the growth of planetesimals, and discuss a relatively recent model of core growth through the accretion of pebbles. We review the basic physics of planet-disk interactions, recent progress in this area, and discuss their role in observed planetary systems. We address the most important effects of planets internal evolution, like cooling and contraction, the mass-luminosity relation, and the bulk composition expressed in the mass-radius and mass-mean density relations.

Aims. We investigate the feasibility of planetesimal growth in circumbinary protoplanetary disks around the observed systems Kepler- 16 and Kepler-34 under the gravitational influence of a precessing eccentric gas disk. Methods. We embed the results of our previous hydrodynamical simulations of protoplanetary disks around binaries into an N-body code to perform 3D, high-resolution, inter-particle gravity-enabled simulations of planetesimal growth and dynamics that include the gravitational force imparted by the gas. Results. Including the full, precessing asymmetric gas disk generates high eccentricity orbits for planetesimals orbiting at the edge of the circumbinary cavity, where the gas surface density and eccentricity have their largest values. The gas disk is able to efficiently align planetesimal pericenters in some regions leading to phased, non-interacting orbits. Outside of these areas eccentric planetesimal orbits become misaligned and overlap leading to crossing orbits and high relative velocities during planetesimal collisions. This can lead to an increase in the number of erosive collisions that far outweighs the number of collisions that result in growth. Gravitational focusing from the static axisymmetric gas disk is weak and does not significantly alter collision outcomes from the gas free case. Conclusions. Due to asymmetries in the gas disk, planetesimals are strongly perturbed onto highly eccentric orbits. Where planetesimals orbits are not well aligned, orbit crossings lead to an increase in the number of erosive collisions. This makes it difficult for sustained planetesimal accretion to occur at the location of Kepler-16b and Kepler-34b and we therefore rule out in-situ growth. This adds further support to our initial suggestions that most circumbinary planets should form further out in the disk and migrate inwards.

Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as the result of an enhanced gas density vortex where the dust is trapped potentially promoting the rapid growth to the planetesimal scale. AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry. We used the NOEMA interferometer to image the 1mm continuum dust emission and the 13CO J=2->1, C18O J=2->1, SO J=56->45 and H2CO J=303->202 rotational lines. Line integrated intensity ratio images are built to investigate the chemical changes within the disk. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO and SO in the dust trap. Our model shows that for densities >10^7 cm^-3, the SO molecules are depleted (directly frozen or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment. Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk physical structure and initial conditions are important drawbacks.

Motivated by lopsided structures observed in some massive transition discs, we have carried out 2D numerical simulations to study vortex structure in massive discs, including the effects of disc self-gravity and the indirect force which is due to the displacement of the central star from the barycenter of the system by the lopsided structure. When only the indirect force is included, we confirm the finding by Mittal & Chiang (2015) that the vortex becomes stronger and can be more than two pressure scale heights wide, as long as the disc-to-star mass ratio is >1%. Such wide vortices can excite strong density waves in the disc and therefore migrate inwards rapidly. However, when disc self-gravity is also considered in simulations, self-gravity plays a more prominent role on the vortex structure. We confirm that when the disc Toomre Q parameter is smaller than pi/(2h), where h is the disc's aspect ratio, the vortices are significantly weakened and their inward migration slows down dramatically. Most importantly, when the disc is massive enough (e.g. Q~3), we find that the lopsided gas structure orbits around the star at a speed significantly slower than the local Keplerian speed. This sub-Keplerian pattern speed can lead to the concentration of dust particles at a radius beyond the lopsided gas structure (as shown in Paper II). Overall, disc self-gravity regulates the vortex structure in massive discs and the radial shift between the gas and dust distributions in vortices within massive discs may be probed by future observations.

Star-planet tidal interactions may result in the excitation of inertial waves in the convective region of stars. In low-mass stars, their dissipation plays a prominent role in the long-term orbital evolution of short-period planets. Turbulent convection can sustain differential rotation in their envelope, with an equatorial acceleration (as in the Sun) or deceleration, which can modify the waves' propagation properties. We explore in this first paper the general propagation properties of free linear inertial waves in a differentially rotating homogeneous fluid inside a spherical shell. We assume that the angular velocity background flow depends on the latitudinal coordinate only, close to what is expected in the external convective envelope of low-mass stars. We use i) an analytical approach in the inviscid case to get the dispersion relation, from which we compute the characteristic trajectories along which energy propagates. This allows us to study the existence of attractor cycles and infer the different families of inertial modes; ii) high-resolution numerical calculations based on a spectral method for the viscous problem. We find that modes that propagate in the whole shell (D modes) behave the same way as with solid-body rotation. However, another family of inertial modes exists (DT modes), which can propagate only in a restricted part of the convective zone. Our study shows that they are less common than D modes and that the characteristic rays and shear layers often focus towards a wedge - or point-like attractor. More importantly, we find that for non-axisymmetric oscillation modes, shear layers may cross a corotation resonance with a local accumulation of kinetic energy. Their damping rate scales very differently from what we obtain for standard D modes and we show an example where it is independent of viscosity (Ekman number) in the astrophysical regime in which it is small.

We investigate the dynamics of large dust grains in massive lopsided transition discs via 2D hydrodynamical simulations including both gas and dust. Our simulations adopt a ring-like gas density profile that becomes unstable against the Rossby-wave instability and forms a large crescent-shaped vortex. When gas self-gravity is discarded, but the indirect force from the displacement of the star by the vortex is included, we confirm that dust grains with stopping times of order the orbital time, which should be typically a few centimetres in size, are trapped ahead of the vortex in the azimuthal direction, while the smallest and largest grains concentrate towards the vortex centre. We obtain maximum shift angles of about 25 degrees. Gas self-gravity accentuates the concentration differences between small and large grains. At low to moderate disc masses, the larger the grains, the farther they are trapped ahead of the vortex. Shift angles up to 90 degrees are reached for 10 cm-sized grains, and we show that such large offsets can produce a double-peaked continuum emission observable at mm/cm wavelengths. This behaviour comes about because the large grains undergo horseshoe U-turns relative to the vortex due to the vortex's gravity. At large disc masses, since the vortex's pattern frequency becomes increasingly slower than Keplerian, small grains concentrate slightly beyond the vortex and large grains form generally non-axisymmetric ring-like structures around the vortex's radial location. Gas self-gravity therefore imparts distinct trapping locations for small and large dust grains which may be probed by current and future observations, and which suggest that the formation of planetesimals in vortices might be more difficult than previously thought.

Tidally-excited inertial waves in stellar convective regions are a key mechanism for tidal dissipation in stars and therefore the evolution of close-in binary or planetary systems. As a first step, we explore here the impact of latitudinal differential rotation on the properties of free inertial modes and identify the different families of modes. We show that they differ from the case of solid-body rotation. Using an analytical approach as well as numerical calculations, we conclude that critical layers (where the Doppler-shifted frequency vanishes) could play a very important role for tidal dissipation.

The Kepler mission's discovery of a number of circumbinary planets orbiting close (a_p < 1.1 au) to the stellar binary raises questions as to how these planets could have formed given the intense gravitational perturbations the dual stars impart on the disk. The gas component of circumbinary protoplanetary disks is perturbed in a similar manner to the solid, planetesimal dominated counterpart, although the mechanism by which disk eccentricity originates differs. This is the first work of a series that aims to investigate the conditions for planet formation in circumbinary protoplanetary disks. We present a number of hydrodynamical simulations that explore the response of gas disks around two observed binary systems: Kepler-16 and Kepler-34. We probe the importance of disk viscosity, aspect-ratio, inner boundary condition, initial surface density gradient, and self-gravity on the dynamical evolution of the disk, as well as its quasi steady-state profile. We find there is a strong influence of binary type on the mean disk eccentricity, e_d, leading to e_d = 0.02 - 0.08 for Kepler-16 and e_d = 0.10 - 0.15 in Kepler-34. The value of alpha-viscosity has little influence on the disk, but we find a strong increase in mean disk eccentricity with increasing aspect-ratio due to wave propagation effects. The choice of inner boundary condition only has a small effect on the surface density and eccentricity of the disk. Our primary finding is that including disk self-gravity has little impact on the evolution or final state of the disk for disks with masses less than 12.5 times that of the minimum-mass solar nebula. This finding contrasts with the results of self-gravity relevance in circumprimary disks, where its inclusion is found to be an important factor in describing the disk evolution.

Planets are thought to form in the gas and dust disks around young stars. In particular, it has been proposed that giant planets can form through the gravitational instability of massive extended disks around intermediate-mass stars. However, we still lack direct observations to constrain this mechanism. We have spatially resolved the 8.6 and 11.2 $\mu$m emission of a massive protoplanetary disk seen edge on around an A star, Gomez's Hamburger (GoHam), using VISIR at the Very Large Telescope. A compact region situated at a projected distance of $350\pm50$ AU south of the central star is found to have a reduced emission.This asymmetry is fully consistent with the presence of a cold density structure, or clump, identified in earlier CO observations, and we derive physical characteristics consistent with those observations: a mass of 0.8-11.4 Jupiter masses (for a dust-to-gas mass ratio of 0.01), a radius of about 10$^2$ astronomical units, and a local density of about $10^{7}$ cm$^{-3}$. Based on this evidence, we argue that this clump, which we call GoHam b, is a promising candidate for a young protoplanet formed by gravitational instability that might be representative of the precursors of massive planets observed around A stars, such as HR 8799 or Beta pictoris. More detailed studies at high angular resolution are needed to better constrain the physical properties of this object to confirm this proposal.

While recent observational progress is converging on the detection of compact regions of thermal emission due to embedded protoplanets, further theoretical predictions are needed to understand the response of a protoplanetary disk to the planet formation radiative feedback. This is particularly important to make predictions for the observability of circumplanetary regions. In this work we use 2D hydrodynamical simulations to examine the evolution of a viscous protoplanetary disk in which a luminous Jupiter-mass planet is embedded. We use an energy equation which includes the radiative heating of the planet as an additional mechanism for planet formation feedback. Several models are computed for planet luminosities ranging from $10^{-5}$ to $10^{-3}$ Solar luminosities. We find that the planet radiative feedback enhances the disk's accretion rate at the planet's orbital radius, producing a hotter and more luminous environement around the planet, independently of the prescription used to model the disk's turbulent viscosity. We also estimate the thermal signature of the planet feedback for our range of planet luminosities, finding that the emitted spectrum of a purely active disk, without passive heating, is appreciably modified in the infrared. We simulate the protoplanetary disk around HD 100546 where a planet companion is located at about 68 au from the star. Assuming the planet mass is 5 Jupiter masses and its luminosity is $\sim 2.5 \times 10^{-4} \, L_\odot$, we find that the radiative feedback of the planet increases the luminosity of its $\sim 5$ au circumplanetary disk from $10^{-5} \, \rm L_\odot$ (without feedback) to $10^{-3} \, \rm L_\odot$, corresponding to an emission of $\sim 1 \, \rm mJy$ in $L^\prime$ band after radiative transfer calculations, a value that is in good agreement with HD 100546b observations.

We propose a mechanism by which dust rings in protoplanetary disks can form and be long-lasting compared to gas rings. This involves the existence of a pressure maximum which traps dust either in between two gap-opening planets or at the outermost gap edge of a single or multiple planet system, combined with the decoupling of large dust particles from the gas. We perform 2D gas hydrodynamical simulations of disks with one and two giant planets which may open deep or partial gaps. A gas ring forms in between two planets such that the surface mass density is higher than on either side of it. This ring is a region of pressure maximum where we expect large grains, which are marginally coupled to the gas and would otherwise be subject to radial drift, to collect. Such a pressure maximum also occurs at the outermost gap edge in a disk with one or more planets. We infer the dust evolution in these regions as the gas disk evolves, to understand the longer term behavior of the resulting dust rings. Over time the gas surface density in the ring(s) decreases, which may cause the larger trapped particles to decouple. Consequently, these particles are expected to stay in ring structure(s) longer than the gas. For a Minimum Mass Solar Nebula model, we expect that millimeter and centimeter-sized grains in the outer O(10) au would most likely undergo this trapping and decoupling process.

While many intermediate- and high-mass main sequence stars are rapidly and differentially rotating, the effects of rotation on oscillation modes are poorly known. In this communication we present a first study of axisymmetric gravito-inertial modes in the radiative zone of a differentially rotating star. We consider a simplified model where the radiative zone of the star is a linearly stratified rotating fluid within a spherical shell, with differential rotation due to baroclinic effects. We solve the eigenvalue problem with high-resolution spectral computations and determine the propagation domain of the waves through the theory of characteristics. We explore the propagation properties of two kinds of modes: those that can propagate in the entire shell and those that are restricted to a subdomain. Some of the modes that we find concentrate kinetic energy around short-period shear layers known as attractors. We describe various geometries for the propagation domains, conditioning the surface visibility of the corresponding modes.

Oscillations have been detected in a variety of stars, including intermediate- and high-mass main sequence stars. While many of these stars are rapidly and differentially rotating, the effects of rotation on oscillation modes are poorly known. In this communication we present a first study on axisymmetric gravito-inertial modes in the radiative zone of a differentially rotating star. These modes probe the deep layers of the star around its convective core. We consider a simplified model where the radiative zone of a star is a linearly stratified rotating fluid within a spherical shell, with differential rotation due to baroclinic effects. We solve the eigenvalue problem with high-resolution spectral simulations and determine the propagation domain of the waves through the theory of characteristics. We explore the propagation properties of two kinds of modes: those that can propagate in the entire shell and those that are restricted to a subdomain. Some of the modes that we find concentrate kinetic energy around short-period shear layers known as attractors. We characterise these attractors by the dependence of their Lyapunov exponent with the \BV frequency of the background and the oscillation frequency of the mode. Finally, we note that, as modes associated with short-period attractors form dissipative structures, they could play an important role for tidal interactions but should be dismissed in the interpretation of observed oscillation frequencies.

Star-planet tidal interactions may result in the excitation of inertial waves in the convective region of stars. Their dissipation plays a prominent role in the long-term orbital evolution of short-period planets. If the star is assumed to be rotating as a solid-body, the waves' Doppler-shifted frequency is restricted to $[-2 \Omega, 2 \Omega]$ ($\Omega$ being the angular velocity of the star) and they can propagate in the entire convective region. However, turbulent convection can sustain differential rotation with an equatorial acceleration (as in the Sun) or deceleration that may modify waves propagation. We thus explore the properties of inertial modes of oscillation in a conically differentially rotating background flow whose angular velocity depends on the latitudinal coordinate only, close to what is expected in the external convective envelope of low-mass stars. We find that their frequency range is broadened by differential rotation, and that they may propagate only in a restricted part of the envelope. In some cases, inertial waves form shear layers around short-period attractor cycles. In others, they exhibit a remarkable behavior when a turning surface or a corotation layer exists in the star. We discuss how all these cases can impact tidal dissipation in stars.

Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors towards high energy, potentially destructive collisions that inhibit the growth of aspiring planets. In this work, we conduct high resolution, inter-particle gravity enabled N-body simulations to investigate the feasibility of planetesimal growth in the Kepler-34 system. We improve upon previous work by including planetesimal disk self-gravity and an extensive collision model to accurately handle inter-planetesimal interactions. We find that super-catastrophic erosion events are the dominant mechanism up to and including the orbital radius of Kepler-34(AB)b, making in-situ growth unlikely. It is more plausible that Kepler-34(AB)b migrated from a region beyond 1.5 AU. Based on the conclusions that we have made for Kepler-34 it seems likely that all of the currently known circumbinary planets have also migrated significantly from their formation location with the possible exception of Kepler-47(AB)c.

We perform one of the first studies into the nonlinear evolution of tidally excited inertial waves in a uniformly rotating fluid body, exploring a simplified model of the fluid envelope of a planet (or the convective envelope of a solar-type star) subject to the gravitational tidal perturbations of an orbiting companion. Our model contains a perfectly rigid spherical core, which is surrounded by an envelope of incompressible uniform density fluid. The corresponding linear problem was studied in previous papers which this work extends into the nonlinear regime, at moderate Ekman numbers (the ratio of viscous to Coriolis accelerations). By performing high-resolution numerical simulations, using a combination of pseudo-spectral and spectral element methods, we investigate the effects of nonlinearities, which lead to time-dependence of the flow and the corresponding dissipation rate. Angular momentum is deposited non-uniformly, leading to the generation of significant differential rotation in the initially uniformly rotating fluid, i.e. the body does not evolve towards synchronism as a simple solid body rotator. This differential rotation modifies the properties of tidally excited inertial waves, changes the dissipative properties of the flow, and eventually becomes unstable to a secondary shear instability provided that the Ekman number is sufficiently small. Our main result is that the inclusion of nonlinearities eventually modifies the flow and the resulting dissipation from what linear calculations would predict, which has important implications for tidal dissipation in fluid bodies. We finally discuss some limitations of our simplified model, and propose avenues for future research to better understand the tidal evolution of rotating planets and stars.

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.

The Kepler mission is dramatically increasing the number of planets known in multi-planetary systems. Many adjacent planets have orbital period ratios near resonant values, with a tendency to be larger than required for exact first-order mean-motion resonances. This intriguing feature has been shown to be a natural outcome of orbital circularization of resonant planetary pairs due to star-planet tidal interactions. However, this feature holds in multi-planetary systems with periods longer than ten days, for which tidal circularization is unlikely to provide efficient divergent evolution of the planets orbits. Gravitational interactions between planets and their parent protoplanetary disk may instead provide efficient divergent evolution. For a planet pair embedded in a disk, we show that interactions between a planet and the wake of its companion can reverse convergent migration, and significantly increase the period ratio from a near-resonant value. Divergent evolution due to wake-planet interactions is particularly efficient when at least one of the planets opens a partial gap around its orbit. This mechanism could help account for the diversity of period ratios in Kepler's multiple systems comprising super-Earth to sub-jovian planets with periods greater than about ten days. Diversity is also expected for planet pairs massive enough to merge their gap. The efficiency of wake-planet interactions is then much reduced, but convergent migration may stall with a variety of period ratios depending on the density structure in the common gap. This is illustrated for the Kepler-46 system, for which we reproduce the period ratio of Kepler-46b and c.

Recent observations from NASA's Kepler mission detected the first planets in circumbinary orbits. The question we try to answer is where these planets formed in the circumbinary disk and how far inside they migrated to reach their present location. We investigate the first and more delicate phase of planet formation when planetesimals accumulate to form planetary embryos. We use the hydrodynamical code FARGO to study the evolution of the disk and of a test population of planetesimals embedded in it. With this hybrid hydrodynamical--N--body code we can properly account for the gas drag force on the planetesimals and for the gravitational force of the disk on them. The numerical simulations show that the gravity of the eccentric disk on the planetesimal swarm excites their eccentricities to values much larger than those induced by the binary perturbations only within 10 AU from the stars. Moreover, the disk gravity prevents a full alignment of the planetesimal pericenters. Both these effects lead to large impact velocities, beyond the critical value for erosion. Planetesimals accumulation in circumbinary disks appears to be prevented close to the stellar pair by the gravitational perturbations of the circumbinary disk. The observed planets possibly formed in the outer regions of the disk and then migrated inside by tidal interaction with the disk.

[Abridged] The migration of low mass planets has been studied in hydrodynamical disc models for more than three decades, but the impact of a magnetic field in the protoplanetary disc is less known. When the disc's magnetic field is strong enough to prevent horseshoe motion, the corotation torque is replaced by a torque arising from magnetic resonances. For weak enough magnetic fields, horseshoe motion and a corotation torque exist, and recent turbulent MHD simulations have reported the existence of a new component of the corotation torque in the presence of a mean toroidal field. The aim of this paper is to investigate the physical origin and the properties of this new corotation torque. We performed MHD simulations of a low mass planet embedded in a 2D laminar disc threaded by a weak toroidal magnetic field, with the effects of turbulence modelled by a viscosity and a resistivity. We confirm that the interaction between the magnetic field and the horseshoe motion results in an additional corotation torque on the planet, which we dub the MHD torque excess. It is caused by the accumulation of the magnetic field along the downstream separatrices of the planet's horseshoe region, which gives rise to an azimuthally asymmetric underdense region at that location. The properties of the MHD torque excess are characterised by varying the slope of the density, temperature and magnetic field profiles, as well as the diffusion coefficients and the strength of the magnetic field. The sign of the MHD torque excess depends on the density and temperature gradients only, and is positive for profiles expected in protoplanetary discs. Its magnitude is in turn mainly determined by the strength of the magnetic field and the turbulent resistivity. The MHD torque excess can be strong enough to reverse migration, even when the magnetic pressure is less than one percent of the thermal pressure.

We investigate the properties of small-amplitude inertial waves propagating in a differentially rotating incompressible fluid contained in a spherical shell. For cylindrical and shellular rotation profiles and in the inviscid limit, inertial waves obey a second-order partial differential equation of mixed type. Two kinds of inertial modes therefore exist, depending on whether the hyperbolic domain where characteristics propagate covers the whole shell or not. The occurrence of these two kinds of inertial modes is examined, and we show that the range of frequencies at which inertial waves may propagate is broader than with solid-body rotation. Using high-resolution calculations based on a spectral method, we show that, as with solid-body rotation, singular modes with thin shear layers following short-period attractors still exist with differential rotation. They exist even in the case of a full sphere. In the limit of vanishing viscosities, the width of the shear layers seems to weakly depend on the global background shear, showing a scaling in E^{1/3} with the Ekman number E, as in the solid-body rotation case. There also exist modes with thin detached layers of width scaling with E^{1/2} as Ekman boundary layers. The behavior of inertial waves with a corotation resonance within the shell is also considered. For cylindrical rotation, waves get dramatically absorbed at corotation. In contrast, for shellular rotation, waves may cross a critical layer without visible absorption, and such modes can be unstable for small enough Ekman numbers.

Low-mass objects embedded in isothermal protoplanetary discs are known to suffer rapid inward Type I migration. In non-isothermal discs, recent work has shown that a decreasing radial profile of the disc entropy can lead to a strong positive corotation torque which can significantly slow down or reverse Type I migration in laminar viscous disc models. It is not clear however how this picture changes in turbulent disc models. The aim of this study is to examine the impact of turbulence on the torque experienced by a low-mass planet embedded in a non-isothermal protoplanetary disc. We particularly focus on the role of turbulence on the corotation torque whose amplitude depends on the efficiency of diffusion processes in the planet's horseshoe region. We performed 2D numerical simulations using a grid-based hydrodynamical code in which turbulence is modelled as stochastic forcing. In order to provide estimations for the viscous and thermal diffusion coefficients as a function of the amplitude of turbulence, we first set up non-isothermal disc models for different values of the amplitude of the turbulent forcing. We then include a low-mass planet and determine the evolution of its running time-averaged torque. We show that in non-isothermal discs, the entropy-related corotation torque can indeed remain unsaturated in the presence of turbulence. For turbulence amplitudes that do not strongly affect the disc temperature profile, we find that the running time-averaged torque experienced by an embedded protoplanet is in fairly good agreement with laminar disc models with appropriate values for the thermal and viscous diffusion coefficients. In discs with turbulence driven by stochastic forcing, the corotation torque therefore behaves similarly as in laminar viscous discs and can be responsible for significantly slowing down or reversing Type I migration.

We study planetesimal evolution in circumbinary disks, focusing on the three systems Kepler 16, 34 and 35 where planets have been discovered recently. We show that for circumbinary planetesimals, in addition to secular forcing, eccentricities evolve on a dynamical timescale, which leads to orbital crossings even in the presence of gas drag. This makes the current locations of the circumbinary Kepler planets hostile to planetesimal accretion. We then present results from simulations including planetesimal formation and dust accretion, and show that even in the most favourable case of 100% efficient dust accretion, in situ growth starting from planetesimals smaller than ~10 km is difficult for Kepler 16b, Kepler 34b and Kepler 35b. These planets were likely assembled further out in the disk, and migrated inward to their current location.

We present two-dimensional hydrodynamical simulations of the tidal interaction between a supermassive black hole binary with moderate mass ratio, and the fossil gas disc where it is embedded. Our study extends previous one-dimensional height-integrated disc models, which predicted that the density of the gas disc between the primary and the secondary black holes should rise significantly during the ultimate stages of the binary's hardening driven by the gravitational radiation torque. This snow-plough mechanism, as we call it, would lead to an increase in the bolometric luminosity of the system prior to the binary merger, which could be detected in conjunction with the gravitational wave signal. We argue here that the snow-plough mechanism is unlikely to occur. In two-dimensions, when the binary's hardening timescale driven by gravitational radiation becomes shorter than the disc's viscous drift timescale, fluid elements in the inner disc get funneled to the outer disc through horseshoe trajectories with respect to the secondary. Mass leakage across the secondary's gap is thus found to be effective and, as a result, the predicted accretion disc luminosity will remain at roughly the same level prior to merger.

Planetary migration is the process by which a forming planet undergoes a drift of its semi-major axis caused by the tidal interaction with its parent protoplanetary disc. One of the key quantities to assess the migration of embedded planets is the tidal torque between the disc and planet, which has two components: the Lindblad torque and the corotation torque. We review the latest results on both torque components for planets on circular orbits, with a special emphasis on the various processes that give rise to additional, large components of the corotation torque, and those contributing to the saturation of this torque. These additional components of the corotation torque could help address the shortcomings that have recently been exposed by models of planet population syntheses. We also review recent results concerning the migration of giant planets that carve gaps in the disc (type II migration) and the migration of sub-giant planets that open partial gaps in massive discs (type III migration).

Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the disc's radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the disc's longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($\pi$). The disc's radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.

The migration of low-mass planets is driven by the differential Lindblad torque and the corotation torque in non-magnetic viscous models of protoplanetary discs. The corotation torque has recently received detailed attention as it may slow down, stall, or reverse migration. In laminar viscous disc models, the long-term evolution of the corotation torque is intimately related to viscous and thermal diffusion processes in the planet's horseshoe region. This paper examines the properties of the corotation torque in discs where MHD turbulence develops as a result of the magnetorotational instability, considering a weak initial toroidal magnetic field. We present results of 3D MHD simulations carried out with two different codes. Non-ideal MHD effects and the disc's vertical stratification are neglected, and locally isothermal disc models are considered. The running time-averaged torque exerted by the disc on a fixed planet is evaluated in three disc models. We first present results with an inner disc cavity (planet trap). As in viscous disc models, the planet is found to experience a positive running time-averaged torque over several hundred orbits, which highlights the existence of an unsaturated corotation torque maintained in the long term in MHD turbulent discs. Two disc models with initial power-law density and temperature profiles are also adopted, in which the time-averaged torque is found to be in decent agreement with its counterpart in laminar viscous disc models with similar viscosity at the planet location. Detailed analysis of the averaged torque density distributions indicates that the differential Lindblad torque takes very similar values in MHD turbulent and laminar viscous discs, and there exists an unsaturated corotation torque in MHD turbulent discs. This analysis also reveals the existence of an additional corotation torque in weakly magnetized discs.

We study the numerical convergence of hydrodynamical simulations of self-gravitating accretion discs, in which a simple cooling law is balanced by shock heating. It is well-known that there exists a critical cooling time scale for which shock heating can no longer compensate for the energy losses, at which point the disc fragments. The numerical convergence of previous results of this critical cooling time scale was questioned recently using Smoothed Particle Hydrodynamics (SPH). We employ a two-dimensional grid-based code to study this problem, and find that for smooth initial conditions, fragmentation is possible for slower cooling as the resolution is increased, in agreement with recent SPH results. We show that this non-convergence is at least partly due to the creation of a special location in the disc, the boundary between the turbulent and the laminar region, when cooling towards a gravito-turbulent state. Converged results appear to be obtained in setups where no such sharp edges appear, and we then find a critical cooling time scale of ~ 4 $\Omega^{-1}$, where $\Omega$ is the local angular velocity.

The observation of massive exoplanets at large separation from their host star, like in the HR 8799 system, challenges theories of planet formation. A possible formation mechanism involves the fragmentation of massive self-gravitating discs into clumps. While the conditions for fragmentation have been extensively studied, little is known of the subsequent evolution of these giant planet embryos, in particular their expected orbital migration. Assuming a single planet has formed by fragmentation, we investigate its interaction with the gravitoturbulent disc it is embedded in. Two-dimensional hydrodynamical simulations are used with a simple prescription for the disc cooling. A steady gravitoturbulent disc is first set up, after which simulations are restarted including a planet with a range of masses approximately equal to the clump's initial mass expected in fragmenting discs. Planets rapidly migrate inwards, despite the stochastic kicks due to the turbulent density fluctuations. We show that the migration timescale is essentially that of type I migration, with the planets having no time to open a gap. In discs with aspect ratio ~ 0.1 at their forming location, planets with a mass comparable to, or larger than Jupiter's can migrate in as short as 10000 years, that is, about 10 orbits at 100 AU. Massive planets formed at large separation from their star by gravitational instability are thus unlikely to stay in place, and should rapidly migrate towards the inner parts of protoplanetary discs, regardless of the planet mass.

We investigate the evolution of a system of two super-Earths with masses < 4 Earth masses embedded in a turbulent protoplanetary disk. The aim is to examine whether or not resonant trapping can occur and be maintained in presence of turbulence and how this depends on the amplitude of the stochastic density fluctuations in the disk. We have performed 2D numerical simulations using a grid-based hydrodynamical code in which turbulence is modelled as stochastic forcing. We assume that the outermost planet is initially located just outside the position of the 3:2 mean motion resonance (MMR) with the inner one and we study the dependance of the resonance stability with the amplitude of the stochastic forcing. For systems of two equal-mass planets we find that in disk models with an effective viscous stress parameter {\alpha} 10^{-3}, damping effects due to type I migration can counteract the effects of diffusion of the resonant angles, in such a way that the 3:2 MMR can possibly remain stable over the disk lifetime. For systems of super-Earths with mass ratio q=m_i/m_o < 1/2, where m_i (m_o) is the mass of the innermost (outermost) planet, the 3:2 MMR is broken in turbulent disks with effective viscous stresses 2x10^{-4}< {\alpha}< 10^{-3} but the planets become locked in stronger p+1:p resonances, with p increasing as the value for {\alpha} increases. For {\alpha}> 2x10^{-3}, the evolution can eventually involve temporary capture in a 8:7 commensurability but no stable MMR is formed. Our results suggest that for values of the viscous stress parameter typical to those generated by MHD turbulence, MMRs between two super-Earths are likely to be disrupted by stochastic density fluctuations. For lower levels of turbulence however, as is the case in presence of a dead-zone, resonant trapping can be maintained in systems with moderate values of the planet mass ratio.

The massive stars in the Galactic center inner arcsecond share analogous properties with the so-called Hot Jupiters. Most of these young stars have highly eccentric orbits, and were probably not formed in-situ. It has been proposed that these stars acquired their current orbits from the tidal disruption of compact massive binaries scattered toward the proximity of the central supermassive black hole. Assuming a binary star formed in a thin gaseous disk beyond 0.1 pc from the central object, we investigate the relevance of disk-satellite interactions to harden the binding energy of the binary, and to drive its inward migration. A massive, equal-mass binary star is found to become more tightly wound as it migrates inwards toward the central black hole. The migration timescale is very similar to that of a single-star satellite of the same mass. The binary's hardening is caused by the formation of spiral tails lagging the stars inside the binary's Hill radius. We show that the hardening timescale is mostly determined by the mass of gas inside the binary's Hill radius, and that it is much shorter than the migration timescale. We discuss some implications of the binary's hardening process. When the more massive (primary) components of close binaries eject most their mass through supernova explosion, their secondary stars may attain a range of eccentricities and inclinations. Such processes may provide an alternative unified scenario for the origin of the kinematic properties of the central cluster and S-stars in the Galactic center as well as the high velocity stars in the Galactic halo.

We study the effects of diffusion on the non-linear corotation torque, or horseshoe drag, in the two-dimensional limit, focusing on low-mass planets for which the width of the horseshoe region is much smaller than the scale height of the disc. In the absence of diffusion, the non-linear corotation torque saturates, leaving only the Lindblad torque. Diffusion of heat and momentum can act to sustain the corotation torque. In the limit of very strong diffusion, the linear corotation torque is recovered. For the case of thermal diffusion, this limit corresponds to having a locally isothermal equation of state. We present some simple models that are able to capture the dependence of the torque on diffusive processes to within 20% of the numerical simulations.

Hydrodynamical simulations of two giant planets embedded in a gaseous disk have shown that in case of a smooth convergent migration they end up trapped into a mean motion resonance. These findings have led to the conviction that the onset of dynamical instability causing close encounters between the planets can occur only after the dissipation of the gas when the eccentricity damping is over. We show that a system of three giant planets may undergo planet-planet scattering when the gaseous disk, with density values comparable to that of the Minimum Mass Solar Nebula, is still interacting with the planets. The hydrodynamical code FARGO--2D--1D is used to model the evolution ofthe disk and planets, modified to properly handle close encounters between the massive bodies. Our simulations predict a variety of different outcomes of the scattering phase which includes orbital exchange, planet merging and scattering of a planet in a hyperbolic orbit. This implies thatthe final fate of a multiplanet system under the action of the disk torques is not necessarily a packed resonant configuration.

In order to reproduce the statistical properties of the observed exoplanets, population synthesis models have shown that the migration of protoplanets should be significantly slowed down, and that processes stalling migration should be at work. Much current theoretical efforts have thus been dedicated to find physical effects that slow down, halt or even reverse migration. Many of these studies rely on the horseshoe drag, whose long term-evolution (saturated or not) is intimately related to the disk viscosity in laminar disk models. We investigate how the horseshoe drag exerted on a low-mass planet is altered by a more realistic treatment of the turbulence in protoplanetary disks. Two-dimensional hydrodynamic simulations are performed with a turbulence model that reproduces the main turbulence properties of three-dimensional magnetohydrodynamic calculations. We find that the horseshoe drag can remain unsaturated on the long term, depending on the turbulence strength. We show that the desaturation of the horseshoe drag by turbulence can be modeled by vortensity diffusion across the time-averaged planet's horseshoe region. At low-turbulence, the running-time averaged torque is in good agreement with the total torque obtained for an equivalent laminar model, with a similar vortensity's diffusion coefficient. At high-turbulence, differences arise due to the time-evolution of the averaged density profile with turbulence.

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.

We study the evolution of circumstellar massive disks around the primary star of a binary system focusing on the computation of disk eccentricity. In particular, we concentrate on its dependence on the binary eccentricity. Self-gravity is included in our numerical simulations. Our standard model assumes a semimajor axis for the binary of 30 AU, the most probable value according to the present binary statistics.

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.

We investigate the tidal interaction between a low-mass planet and a self-gravitating protoplanetary disk, by means of two-dimensional hydrodynamic simulations. We first show that considering a planet freely migrating in a disk without self-gravity leads to a significant overestimate of the migration rate. The overestimate can reach a factor of two for a disk having three times the surface density of the minimum mass solar nebula. Unbiased drift rates may be obtained only by considering a planet and a disk orbiting within the same gravitational potential. In a second part, the disk self-gravity is taken into account. We confirm that the disk gravity enhances the differential Lindblad torque with respect to the situation where neither the planet nor the disk feels the disk gravity. This enhancement only depends on the Toomre parameter at the planet location. It is typically one order of magnitude smaller than the spurious one induced by assuming a planet migrating in a disk without self-gravity. We confirm that the torque enhancement due to the disk gravity can be entirely accounted for by a shift of Lindblad resonances, and can be reproduced by the use of an anisotropic pressure tensor. We do not find any significant impact of the disk gravity on the corotation torque.

We consider the angular momentum exchange at the corotation resonance between a two-dimensional gaseous disk and a uniformly rotating external potential, assuming that the disk flow is adiabatic. We first consider the linear case for an isolated resonance, for which we give an expression of the corotation torque that involves the pressure perturbation, and which reduces to the usual dependence on the vortensity gradient in the limit of a cold disk. Although this expression requires the solution of the hydrodynamic equations, it provides some insight into the dynamics of the corotation region. In the general case, we find an additional dependence on the entropy gradient at corotation. This dependence is associated to the advection of entropy perturbations. These are not associated to pressure perturbations. They remain confined to the corotation region, where they yield a singular contribution to the corotation torque. In a second part, we check our torque expression by means of customized two-dimensional hydrodynamical simulations. In a third part, we contemplate the case of a planet embedded in a Keplerian disk, assumed to be adiabatic. We find an excess of corotation torque that scales with the entropy gradient, and we check that the contribution of the entropy perturbation to the torque is in agreement with the expression obtained from the linear analysis. We finally discuss some implications of the corotation torque expression for the migration of low mass planets in the regions of protoplanetary disks where the flow is radiatively inefficient on the timescale of the horseshoe U-turns.