We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.

A massive planet in a protoplanetary disc will open a gap in the disc material. A steep gap edge can be hydrodynamically unstable, which results in the formation of vortices that can act as tracers for the presence of planets in observational results. However, in a viscous disc, the potential formation of these vortices is dependent on the timescale over which the massive planet accretes mass and with a sufficiently long timescale it is possible for no vortices to form. Hence, there is a connection between the presence of vortices and the growth timescale of the planet and it may therefore be possible to exclude a planetary interpretation of observed structure from the absence of vortices. We have investigated the effect of the planet growth timescale on vortex formation for a range of planet masses and viscosities and have found an approximate relation between the planet mass, viscosity and planet growth timescale for which vortices are not formed within the disc. We then interpret these results in the light of recent observations. We have also found that planets do not need to be close to a Jupiter mass to form vortices in the disc if these discs have low viscosity, as these can be caused by planets as small as a few Neptune masses.

Convergent migration involving multiple planets embedded in a viscous protoplanetary disc is expected to produce a chain of planets in mean motion resonances, but the multiplanet systems observed by the Kepler spacecraft are generally not in resonance. We demonstrate that under equivalent conditions, where in a viscous disc convergent migration will form a long-term stable system of planets in a chain of mean motion resonances, migration in an inviscid disc often produces a system which is highly dynamically unstable. In particular, if planets are massive enough to significantly perturb the disc surface density and drive vortex formation, the smooth capture of planets into mean motion resonances is disrupted. As planets pile up in close orbits, not protected by resonances, close encounters increase the probability of planet-planet collisions, even while the gas disc is still present. While inviscid discs often produce unstable non-resonant systems, stable, closely packed, non-resonant systems can also be formed. Thus, when examining the expectation for planet migration to produce planetary systems in mean motion resonances, the effective turbulent viscosity of the protoplanetary disc is a key parameter.

Motivated by the stability of dust laden vortices, in this paper we study the terminal velocity approximation equations for a gas coupled to a pressureless dust fluid and present a numerical solver for the equations embedded in the FARGO3D hydrodynamics code. We show that for protoplanetary discs it is possible to use the baricenter velocity in the viscous stress tensor, making it trivial to simulate viscous dusty protoplanetary discs with this model. We also show that the terminal velocity model breaks down around shocks, becoming incompatible with the two fluid model it is derived from. Finally we produce a set of test cases for numerical schemes and demonstrate the performance of our code on these tests. Our implementation embedded in FARGO3D using an unconditionally stable explicit integrator is fast, and exhibits the desired second order spatial convergence for smooth problems.

We present the highest resolution study to date of super-Earths migrating in inviscid and low-viscosity discs, motivated by the connection to laminar, wind-driven models of protoplanetary discs. Our models unveil the critical role of vortices in determining the migration behaviour for partial gap-opening planets. Vortices form in pressure maxima at gap edges, and prevent the disc-feedback stopping of migration for intermediate planets in low-viscosity and inviscid discs, contrary to the concept of the `inertial limit' or `disc feedback' halting predicted from analytical models. Vortices may also form in the corotation region, and can dramatically modify migration behaviour through direct gravitational interaction with the planet. These features become apparent at high resolution, and for all but the highest viscosities there exist significant difficulties in obtaining numerically converged results. The migration of partial gap-opening planets, however, clearly becomes chaotic for sufficiently low viscosities. At moderate viscosity, a smooth disc-feedback regime is found in which migration can slow substantially, and the migration time-scale observed corresponds to migration being driven by diffusive relaxation of the gap edges. At high viscosity classical Type I migration is recovered. For Jupiter-analogue planets in inviscid discs, a wide, deep gap is formed. Transient Type II migration occurs over radial length-scales corresponding to the gap width, beyond which migration can stall. Finally, we examine the particle trapping driven by structures left in inviscid discs by a migrating planet, and find that particle traps in the form of multiple rings and vortices can persist long after the planet has passed. In this case, the observation of particle traps by submillimetre interferometers such as ALMA cannot be used to infer the current presence of an adjacent planet.

We study the hydrodynamical stability of the laminar flows associated with warped astrophysical discs using numerical simulations of warped shearing boxes. We recover linear growth rates reported previously due to a parametric resonance of inertial waves, and show that the nonlinear saturated state can significantly reduce the laminar flows, meaning that the warp would evolve on much longer time scales than would be concluded from the internal torques due to these laminar flows. Towards larger warp amplitudes, we find first of all a reversal of angular momentum flux, indicating that the mass distribution would evolve in an anti-diffusive manner, and second that the linear growth rates disappear, possibly because of the very strong shear in the laminar flows in this regime. For discs with small enough viscosity, a nonlinear state can still be found when linear growth rates are absent by introducing a large enough perturbation, either by starting from a nonlinear state obtained at smaller warp amplitude, or by starting from a state with no laminar flows.

A massive planet in a protoplanetary disc will open a gap in the disc material which acts as a transition between Type I and Type II planetary migration. Type II migration is slower than Type I migration, however it is still desirable to slow down Type II migration to allow gas giant planets with semi-major axis in the range 5 to 10AU to exist, similarly to our Solar system. We investigate a method of slowing down and reversing Type II migration by heating the outer gap edge due to incident radiation from the central star. Using an approximate vertically averaged heating method we find that Type II migration can be slowed or in extreme cases reversed if we assume near maximum allowed irradiation from the central star. Therefore, we believe this is a very interesting phenomenon that should be investigated in greater detail using three dimensional hydrodynamic and radiative transfer simulations.

Recent developments in non-ideal magnetohydrodynamic simulations of protoplanetary disks suggest that instead of being traditional turbulent (viscous) accretion disks, they have a largely laminar flow with accretion driven by large-scale wind torques. These disks are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields. We have examined the dynamics of the corotation region of a low mass planet embedded in such a disk and the evolution of the associated migration torque. These disks lack strong turbulence and associated turbulent diffusion, and the presence of a magnetic field and radial gas flow presents a situation outside the applicability of previous corotation torque theory. We summarize the analytical analysis of the corotation torque, give details on the numerical methods used, and in particular the relative merits of different numerical schemes for the inviscid problem.

We examine the migration of low mass planets in laminar protoplanetary discs, threaded by large scale magnetic fields in the dead zone that drive radial gas flows. As shown in Paper I, a dynamical corotation torque arises due to the flow-induced asymmetric distortion of the corotation region and the evolving vortensity contrast between the librating horseshoe material and background disc flow. Using simulations of laminar torqued discs containing migrating planets, we demonstrate the existence of the four distinct migration regimes predicted in Paper I. In two regimes, the migration is approximately locked to the inward or outward radial gas flow, and in the other regimes the planet undergoes outward runaway migration that eventually settles to fast steady migration. In addition, we demonstrate torque and migration reversals induced by midplane magnetic stresses, with a bifurcation dependent on the disc surface density. We develop a model for fast migration, and show why the outward runaway saturates to a steady speed, and examine phenomenologically its termination due to changing local disc conditions. We also develop an analytical model for the corotation torque at late times that includes viscosity, for application to discs that sustain modest turbulence. Finally, we use the simulation results to develop torque prescriptions for inclusion in population synthesis models of planet formation.

Motivated by models suggesting that the inner planet forming regions of protoplanetary discs are predominantly lacking in viscosity-inducing turbulence, and are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields, we examine the dynamics of the corotation region of a low-mass planet in such an environment. The corotation torque in an inviscid, isothermal, dead zone ought to saturate, with the libration region becoming both symmetrical and of a uniform vortensity, leading to fast inward migration driven by the Lindblad torques alone. However, in such a low viscosity situation, the material on librating streamlines essentially preserves its vortensity. If there is relative radial motion between the disc gas and the planet, the librating streamlines will no longer be symmetrical. Hence, if the gas is torqued by a large scale magnetic field so that it undergoes a net inflow or outflow past the planet, driving evolution of the vortensity and inducing asymmetry of the corotation region, the corotation torque can grow, leading to a positive torque. In this paper we treat this effect by applying a symmetry argument to the previously studied case of a migrating planet in an inviscid disc. Our results show that the corotation torque due to a laminar Hall-induced magnetic field in a dead zone behaves quite differently from that studied previously for a viscous disc. Furthermore, the magnetic field induced corotation torque and the dynamical corotation torque in a low viscosity disc can be regarded as one unified effect.

Young protoplanetary discs and the outer radii of active galactic nucleii may be subject to gravitational instability and, as a consequence, fall into a `gravitoturbulent' state. While in this state, appreciable angular momentum can be transported. Alternatively, the gas may collapse into bound clumps, the progenitors of planets or stars. In this paper, we numerically characterize the properties of 3D gravitoturbulence, focussing especially on its dependence on numerical parameters (resolution, domain size) and its excitation of small-scale dynamics. Via a survey of vertically stratified shearing box simulations with PLUTO and RODEO, we find (a) evidence that certain gravitoturbulent properties are independent of horizontal box size only when the box is larger than $\simeq 40 H_0$, where $H_0$ is the height scale, (b) at high resolution, small-scale isotropic turbulence appears off the midplane around $z\simeq 0.5 -1 H_0$, and (c) this small-scale dynamics results from a parametric instability, involving the coupling of inertial waves with a large-scale axisymmetric epicyclic mode. This mode oscillates at a frequency close to $\Omega$ and is naturally excited by gravito-turbulence, via a nonlinear process to be determined. The small-scale turbulence we uncover has potential implications for a wide range of disc physics, e.g. turbulent saturation levels, fragmentation, turbulent mixing, and dust settling.

It is known that an embedded massive planet will open a gap in a protoplanetary disc via angular momentum exchange with the disc material. The resulting surface density profile of the disc is investigated for one dimensional and two dimensional disc models and, in agreement with previous work, it is found that one dimensional gaps are significantly deeper than their two dimensional counterparts for the same initial conditions. We find, by applying one dimensional torque density distributions to two dimensional discs containing no planet, that the excitement of the Rossby wave instability and the formation of Rossby vortices play a critical role in setting the equilibrium depth of the gap. Being a two dimensional instability, this is absent from one dimensional simulations and does not limit the equilibrium gap depth there. We find similar gap depths between two dimensional gaps formed by torque density distributions, in which the Rossby wave instability is present, and two dimensional planet gaps, in which no Rossby wave instability is present. This can be understood if the planet gap is maintained at marginal stability, even when there is no obvious Rossby wave instability present. Further investigation shows the final equilibrium gap depth is very sensitive to the form of the applied torque density distribution, and using improved one dimensional approximations from three dimensional simulations can go even further to reducing the discrepancy between one and two dimensional models, especially for lower mass planets. This behaviour is found to be consistent across discs with varying parameters.

We present a new method for numerical hydrodynamics which uses a multidimensional generalisation of the Roe solver and operates on an unstructured triangular mesh. The main advantage over traditional methods based on Riemann solvers, which commonly use one-dimensional flux estimates as building blocks for a multidimensional integration, is its inherently multidimensional nature, and as a consequence its ability to recognise multidimensional stationary states that are not hydrostatic. A second novelty is the focus on Graphics Processing Units (GPUs). By tailoring the algorithms specifically to GPUs we are able to get speedups of 100-250 compared to a desktop machine. We compare the multidimensional upwind scheme to a traditional, dimensionally split implementation of the Roe solver on several test problems, and we find that the new method significantly outperforms the Roe solver in almost all cases. This comes with increased computational costs per time step, which makes the new method approximately a factor of 2 slower than a dimensionally split scheme acting on a structured grid.

We present a study of four different formation scenarios that may be able to explain the origin of the recently announced planet (`Proxima b') orbiting Proxima Centauri. The aim is to examine how the formation scenarios differ in their predictions for the multiplicity of the Proxima system, the water/volatile content of Proxima b and its orbital eccentricity, so that these can be tested by future observations. A scenario of in situ formation via giant impacts from a locally enhanced disc of planetary embryos/planetesimals predicts that Proxima b will be in a multiplanet system with a measurably finite eccentricity. Assuming that the local solid enhancement needed to form a Proxima b analogue arises because of the inwards drift of solids in the form of small planetesimals/boulders, this scenario also results in Proxima b analogues being only moderately endowed with water/volatiles. A scenario in which multiple embryos form, migrate and mutually collide within a gas disc results in Proxima b being a member of a multiple system, possibly displaying mean motion resonances, but where the constituent members are Ocean planets due to accretion occurring mainly outside of the snowline. A scenario in which a single accreting embryo forms outside the snowline, and migrates inwards while accreting planetesimals/pebbles results in Proxima b being an isolated Ocean planet on a circular orbit. A scenario in which Proxima b formed via pebble accretion interior to the snowline produces a dry planet on a circular orbit. Future observations that characterise the physical and orbital properties of Proxima b, and any additional planets in the system, will provide valuable insights into the formation history of this neighbouring planetary system.

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($\alpha$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

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.

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.

In this note we discuss the main results of a study of a massive binary with unequal mass ratio, q, embedded in an accretion disk, with its orbital rotation being opposed to that of the disk. When the mass ratio is sufficiently large, a gap opens in the disk, but the mechanism of gap formation is very different from the prograde case. Inward migration occurs on a timescale of t_ev ~ M_p/(dot M), where M_p is the mass of the less massive component (the perturber), and dot M is the accretion rate. When q<< 1, the accretion takes place mostly onto the more massive component, with the accretion rate onto the perturber being smaller than, or of order of, q^(1/3)M. However, this rate increases when supermassive binary black holes are considered and gravitational wave emission is important. We estimate a typical duration of time for which the accretion onto the perturber and gravitational waves could be detected.

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.

One of the striking discoveries of protoplanetary disc research in recent years are the spiral arms seen in several transitional discs in polarised scattered light. An interesting interpretation of the observed spiral features is that they are density waves launched by one or more embedded (proto-)planets in the disc. In this paper we investigate whether planets can be held responsible for the excitation mechanism of the observed spirals. We use locally isothermal hydrodynamic simulations as well as analytic formulae to model the spiral waves launched by planets. Then H-band scattered light images are calculated using a 3D continuum radiative transfer code to study the effect of surface density and pressure scale height perturbation on the detectability of the spirals. We find that a relative change of about 3.5 in the surface density is required for the spirals to be detected with current telescopes in the near-infrared for sources at the distance of typical star-forming regions (140pc). This value is a factor of eight higher than what is seen in hydrodynamic simulations. We also find that a relative change of only 0.2 in pressure scale height is sufficient to create detectable signatures under the same conditions. Therefore, we suggest that the spiral arms observed to date in protoplanetary discs are the results of changes in the vertical structure of the disc (e.g. pressure scale height perturbation) instead of surface density perturbations.

Supermassive black hole binaries may form as a consequence of galaxy mergers. Both prograde and retrograde orbits have been proposed. We study a binary of a small mass ratio, q, in a retrograde orbit immersed in and interacting with a gaseous accretion disk in order to estimate time scales for inward migration leading to coalescence and the accretion rate to the secondary component. We employ both semi-analytic methods and two dimensional numerical simulations, focusing on the case where the binary mass ratio is small but large enough to significantly perturb the disk. We develop the theory of type I migration for this case and determine conditions for gap formation finding that then inward migration occurs on a time scale equal to the time required for one half of the secondary mass to be accreted through the unperturbed disk, with accretion onto the secondary playing only a minor role. The semi-analytic and fully numerical approaches are in good agreement, the former being applicable over long time scales. Inward migration induced by interaction with the disk alleviates the final parsec problem. Accretion onto the secondary does not significantly affect the orbital evolution, but may have observational consequences for high accretion efficiency. The binary may then appear as two sources of radiation rotating around each other. This study should be extended to consider orbits with significant eccentricity and the effects of gravitational radiation at small length scales. Note too that torques acting between a circumbinary disk and a retrograde binary orbit may cause the mutual inclination to increase on a timescale that can be similar to, or smaller than that for orbital evolution, depending on detailed parameters. This is also an aspect for future study (abridged).

We study torques on migrating low-mass planets in locally isothermal discs. Previous work on low-mass planets generally kept the planet on a fixed orbit, after which the torque on the planet was measured. In addition to these static torques, when the planet is allowed to migrate it experiences dynamical torques, which are proportional to the migration rate and whose sign depends on the background vortensity gradient. We show that in discs a few times more massive than the Minimum Mass Solar Nebula, these dynamical torques can have a profound impact on planet migration. Inward migration can be slowed down significantly, and if static torques lead to outward migration, dynamical torques can take over, taking the planet beyond zero-torque lines set by saturation of the corotation torque in a runaway fashion. This means the region in non-isothermal discs where outward migration is possible can be larger than what would be concluded from static torques alone.

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.

Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.

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-36 system consists of two planets that are spaced unusually close together, near the 7:6 mean motion resonance. While it is known that mean motion resonances can easily form by convergent migration, Kepler-36 is an extreme case due to the close spacing and the relatively high planet masses of 4 and 8 times that of the Earth. In this paper, we investigate whether such a system can be obtained by interactions with the protoplanetary disc. These discs are thought to be turbulent and exhibit density fluctuations which might originate from the magneto-rotational instability. We adopt a realistic description for stochastic forces due to these density fluctuations and perform both long term hydrodynamical and N-body simulations. Our results show that planets in the Kepler-36 mass range can be naturally assembled into a closely spaced planetary system for a wide range of migration parameters in a turbulent disc similar to the minimum mass solar nebula. The final orbits of our formation scenarios tend to be Lagrange stable, even though large parts of the parameter space are chaotic and unstable.

Giant planets form in protoplanetary disks while these disks are still gas-rich, and can reveal their presence through the annular gaps they carve out. HD 100546 is a gas-rich disk with a wide gap between between a radius of ~1 and 13 AU, possibly cleared out by a planetary companion or planetary system. We want to identify the nature of the unseen companion near the far end of the disk gap. We use mid-infrared interferometry at multiple baselines to constrain the curvature of the disk wall at the far end of the gap. We use 2D hydrodynamical simulations of embedded planets and brown dwarfs to estimate viscosity of the disk and the mass of a companion close to the disk wall. We find that the disk wall at the far end of the gap is not vertical, but rounded-off by a gradient in the surface density. Such a gradient can be reproduced in hydrodynamical simulations with a single, heavy companion (>=30...80 M_Jup) while the disk has viscosity of at least \alpha{} >~ 5 10^-3. Taking into account the changes in the temperature structure after gap opening reduces the lower limit on the planet mass and disk viscosity to 20 M_Jup and \alpha{} = 2 10^-3. The object in the disk gap of HD 100546 that shapes the disk wall is most likely a 60(+20)(-40) M_Jup brown dwarf, while the disk viscosity is estimated to be at least \alpha{} = 2 10^-3. The disk viscosity is an important factor in estimating planetary masses from disk morphologies: more viscous disks need heavier planets to open an equally deep gap.

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 study numerical convergence in local two-dimensional hydrodynamical simulations of self-gravitating accretion discs with a simple cooling law. It is well-known that there exists a steady gravito-turbulent state, in which cooling is balanced by dissipation of weak shocks, with a net outward transport of angular momentum. Previous results indicated that if cooling is too fast (typical time scale 3/Omega, where Omega is the local angular velocity), this steady state can not be maintained and the disc will fragment into gravitationally bound clumps. We show that, in the two-dimensional local approximation, this result is in fact not converged with respect to numerical resolution and longer time integration. Irrespective of the cooling time scale, gravito-turbulence consists of density waves as well as transient clumps. These clumps will contract because of the imposed cooling, and collapse into bound objects if they can survive for long enough. Since heating by shocks is very local, the destruction of clumps is a stochastic process. High numerical resolution and long integration times are needed to capture this behaviour. We have observed fragmentation for cooling times up to 20/Omega, almost a factor 7 higher than in previous simulations. Fully three-dimensional simulations with a more realistic cooling prescription are necessary to determine the effects of the use of the two-dimensional approximation and a simple cooling law.

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 consider the radial migration of vortices in two-dimensional isothermal gaseous disks. We find that a vortex core, orbiting at the local gas velocity, induces velocity perturbations that propagate away from the vortex as density waves. The resulting spiral wave pattern is reminiscent of an embedded planet. There are two main causes for asymmetries in these wakes: geometrical effects tend to favor the outer wave, while a radial vortensity gradient leads to an asymmetric vortex core, which favors the wave at the side that has the lowest density. In the case of asymmetric waves, which we always find except for a disk of constant pressure, there is a net exchange of angular momentum between the vortex and the surrounding disk, which leads to orbital migration of the vortex. Numerical hydrodynamical simulations show that this migration can be very rapid, on a time scale of a few thousand orbits, for vortices with a size comparable to the scale height of the disk. We discuss the possible effects of vortex migration on planet formation scenarios.

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.

Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii with zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.

We study the collisional evolution of km-sized planetesimals in tight binary star systems to investigate whether accretion towards protoplanets can proceed despite the strong gravitational perturbations from the secondary star. The orbits of planetesimals are numerically integrated in two dimensions under the influence of the two stars and gas drag. The masses and orbits of the planetesimals are allowed to evolve due to collisions with other planetesimals and accretion of collisional debris. In addition, the mass in debris can evolve due to planetesimal-planetesimal collisions and the creation of new planetesimals. We show that it is possible in principle for km-sized planetesimals to grow by two orders of magnitude in size if the efficiency of planetesimal formation is relatively low. We discuss the limitations of our two-dimensional approach.

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.

One class of protoplanetary disc models, the X-wind model, predicts strongly subkeplerian orbital gas velocities, a configuration that can be sustained by magnetic tension. We investigate disc-planet interactions in these subkeplerian discs, focusing on orbital migration for low-mass planets and gap formation for high-mass planets. We use linear calculations and nonlinear hydrodynamical simulations to measure the torque and look at gap formation. In both cases, the subkeplerian nature of the disc is treated as a fixed external constraint. We show that, depending on the degree to which the disc is subkeplerian, the torque on low-mass planets varies between the usual Type I torque and the one-sided outer Lindblad torque, which is also negative but an order of magnitude stronger. In strongly subkeplerian discs, corotation effects can be ignored, making migration fast and inward. Gap formation near the planet's orbit is more difficult in such discs, since there are no resonances close to the planet accommodating angular momentum transport. In stead, the location of the gap is shifted inwards with respect to the planet, leaving the planet on the outside of a surface density depression. Depending on the degree to which a protoplanetary disc is subkeplerian, disc-planet interactions can be very different from the usual Keplerian picture, making these discs in general more hazardous for young planets.

We study the coorbital flow for embedded, low mass planets. We provide a simple semi-analytic model for the corotation region, which is subsequently compared to high resolution numerical simulations. The model is used to derive an expression for the half-width of the horseshoe region, x_s, which in the limit of zero softening is given by x_s/r_p = 1.68(q/h)^(1/2), where q is the planet to central star mass ratio, h is the disc aspect ratio and r_p the orbital radius. This is in very good agreement with the same quantity measured from simulations. This result is used to show that horseshoe drag is about an order of magnitude larger than the linear corotation torque in the zero softening limit. Thus the horseshoe drag, the sign of which depends on the gradient of specific vorticity, is important for estimates of the total torque acting on the planet. We further show that phenomena, such as the Lindblad wakes, with a radial separation from corotation of ~ a pressure scale height H can affect x_s, even though for low-mass planets x_s << H. The effect is to distort streamlines and to reduce x_s through the action of a back pressure. This effect is reduced for smaller gravitational softening parameters and planets of higher mass, for which x_s becomes comparable to H.

We study the torque on low mass protoplanets on fixed circular orbits, embedded in a protoplanetary disc in the isothermal limit. For low mass protoplanets and large viscosity the corotation torque behaves as expected from linear theory. However, when the viscosity becomes small enough to enable horseshoe turns to occur, the linear corotation torque exists only temporarily after insertion of a planet into the disc, being replaced by the horseshoe drag first discussed by Ward. This happens after a time that is equal to the horseshoe libration period reduced by a factor amounting to about twice the disc aspect ratio. This torque scales with the radial gradient of specific vorticity, as does the linear torque, but we find it to be many times larger. If the viscosity is large enough for viscous diffusion across the coorbital region to occur within a libration period, we find that the horseshoe drag may be sustained. If not, the corotation torque saturates leaving only the linear Lindblad torques. As the magnitude of the non linear coorbital torque (horseshoe drag) is always found to be larger than the linear torque, we find that the sign of the total torque may change even for for mildly positive surface density gradients. In combination with a kinematic viscosity large enough to keep the torque from saturating, strong sustained deviations from linear theory and outward or stalled migration may occur in such cases (abridged).

We study the disc planet interactions of low-mass protoplanets embedded in a circumstellar disc. We extend the standard theory of planet migration from the usual locally isothermal assumption to include non-barotropic effects, focusing on the validity of linear theory. We compared solutions of the linear equations with results from non-linear hydrodynamic simulations, where in both cases we adopted a background entropy gradient and solved the energy equation. We show that the migration behavior of embedded planets depends critically on the background radial entropy gradient in the disc. The presence of such a gradient not only changes the corotation torque on the planet, but also always guarantees a departure from linear behavior, which gives a singular density response at corotation, in the absence of thermal or viscous diffusion. A negative entropy gradient tends to give rise to positive, non-linear corotation torques apparently produced as material executes horseshoe turns at approximately constant entropy. These torques have no counterpart in linear theory, but can be strong enough to push the planet outwards until saturation starts to occur after a horseshoe libration period. Increased thermal diffusion acts to reduce these non-linear torques, but, at the same time, it can help to prevent their saturation. In combination with a small kinematic viscosity that is able to maintain a smooth density profile the positive torque could be sustained.

Observations of extrasolar planets reveal that planets can be found in close binary systems, where the semi-major axis of the binary orbit is less than 20 AU. The existence of these planets challenges planet formation theory, because the strong gravitational perturbations due to the companion increase encounter velocities between planetesimals and make it difficult for them to grow through accreting collisions. We study planetesimal encounter velocities in binary systems, where the planetesimals are embedded in a circumprimary gas disc that is allowed to evolve under influence of the gravitational perturbations of the companion star. We find that the encounter velocities between planetesimals of different size strongly depend on the gas disc eccentricity. In all cases studied, inclusion of the full gas dynamics increases the encounter velocity compared to the case of a static, circular gas disc. Full numerical parameter exploration is still impossible, but we derive analytical formulae to estimate encounter velocities between bodies of different sizes given the gas disc eccentricity. The gas dynamical evolution of a protoplanetary disc in a binary system tends to make planetesimal accretion even more difficult than in a static, axisymmetric gas disc.

We study the interaction of a low-mass planet with a protoplanetary disk with a realistic treatment of the energy balance by doing radiation-hydrodynamical simulations. We look at accretion and migration rates and compare them to isothermal studies. We used a three-dimensional version of the hydrodynamical method RODEO, together with radiative transport in the flux-limited diffusion approach. The accretion rate, as well as the torque on the planet, depend critically on the ability of the disk to cool efficiently. For densities appropriate to 5 AU in the solar nebula, the accretion rate drops by more than an order of magnitude compared to isothermal models, while at the same time the torque on the planet is positive, indicating outward migration. It is necessary to lower the density by a factor of 2 to recover inward migration and more than 2 orders of magnitude to recover the usual Type I migration. The torque appears to be proportional to the radial entropy gradient in the unperturbed disk. These findings are critical for the survival of protoplanets, and they should ultimately find their way into population synthesis models.

We study the accretion of dust particles of various sizes onto embedded massive gas giant planets, where we take into account the structure of the gas disk due to the presence of the planet. The accretion rate of solids is important for the structure of giant planets: it determines the growth rate of the solid core that may be present as well as their final enrichment in solids. We use the RODEO hydrodynamics solver to solve the flow equations for the gas, together with a particle approach for the dust. The solver for the particles' equations of motion is implicit with respect to the drag force, which allows us to treat the whole dust size spectrum. We find that dust accretion is limited to the smallest particle sizes. The largest particles get trapped in outer mean-motion resonances with the planet, while particles of intermediate size are pushed away from the orbit of the planet by the density structure in the gas disk. Only particles smaller than approximately s_max =10 micron may accrete on a planet with the mass of Jupiter. For a ten times less massive planet s_max=100 micron. The strongly reduced accretion of dust makes it very hard to enrich a newly formed giant planet in solids.

Aims: We investigate the effect of including a proper energy balance on the interaction of a low-mass planet with a protoplanetary disk. Methods: We use a three-dimensional version of the RODEO method to perform hydrodynamical simulations including the energy equation. Radiation is included in the flux-limited diffusion approach. Results: The sign of the torque depends sensitively on the ability of the disk to radiate away the energy generated in the immediate surroundings of the planet. In the case of high opacity, corresponding to the dense inner regions of protoplanetary disks, migration is directed \emph{outward}, instead of the usual inward migration that was found in locally isothermal disks. For low values of the opacity we recover inward migration, and we show that torques originating in the coorbital region are responsible for the change in migration direction.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

We study the dynamics of gas and dust in a protoplanetary disk in the presence of embedded planets. We investigate the conditions for dust-gap formation in terms of particle size and planetary mass. We also monitor the amount of dust that is accreted by the planet relative to the amount of gas, which is an important parameter in determining the enrichment of solids in giant planets compared to the solid content of the central star. We use a new two-fluid hydrodynamics code to solve the flow equations for both gas and dust. For the gas, we use a Godunov-type scheme with an approximate Riemann solver (the Roe solver). The dust is treated as a pressureless fluid by essentially the same numerical method as is used for the gas. We find that it only takes a planet of 0.05 Jupiter masses to open up a gap in a disk with a significant population of mm-sized particles. Dust particles larger than 150 micron participate in gap formation. We also find that the formation of the gap severely slows down dust accretion compared to that in the gas. Therefore, it is not possible to enrich a newly formed giant planet in solids, if these solids are contained in particles with sizes from 150 micron to approximately 10 cm.

In this paper we describe a new method for studying the hydrodynamical problem of a planet embedded in a gaseous disk. We use a finite volume method with an approximate Riemann solver (the Roe solver), together with a special way to integrate the source terms. This new source term integration scheme sheds new light on the Coriolis instability, and we show that our method does not suffer from this instability. The first results on flow structure and gap formation are presented, as well as accretion and migration rates. For Mpl < 0.1 M_J and Mpl > 1.0 M_J (M_J = Jupiter's mass) the accretion rates do not depend sensitively on numerical parameters, and we find that within the disk's lifetime a planet can grow to 3-4 M_J. In between these two limits numerics play a major role, leading to differences of more than 50 % for different numerical parameters. Migration rates are not affected by numerics at all as long as the mass inside the Roche lobe is not considered. We can reproduce the Type I and Type II migration for low-mass and high-mass planets, respectively, and the fastest moving planet of 0.1 M_J has a migration time of only 2.0 10^4 yr.

We investigate the interaction of gas and dust in a protoplanetary disk in the presence of a massive planet using a new two-fluid hydrodynamics code. In view of future observations of planet-forming disks we focus on the condition for gap formation in the dust fluid. While only planets more massive than 1 Jupiter mass (MJ) open up a gap in the gas disk, we find that a planet of 0.1 MJ already creates a gap in the dust disk. This makes it easier to find lower-mass planets orbiting in their protoplanetary disk if there is a significant population of mm-sized particles.