Recent work has suggested that the net gravitational force acting on a massive and luminous perturber travelling through a gaseous and opaque medium can have same direction as the perturber's motion (an effect sometimes called negative dynamical friction). Analytic results were obtained using a linear analysis and were later confirmed by means of non-linear numerical simulations which did not resolve the flow within the Bondi sphere of the perturber, hence effectively restricted to weakly perturbed regions of the flow (paper~I). Here we present high resolution simulations, using either 3D Cartesian or 2D cylindrical meshes that resolve the flow within the Bondi sphere. We perform a systematic study of the force as a function of the perturber's mass and luminosity, in the subsonic regime. We find that perturbers with mass $M$ smaller than a few $M_c\sim \chi c_s/G$ are subjected to a thermal force with a magnitude in good agreement with linear theory ($\chi$ being the thermal diffusivity of the medium, $c_s$ the adiabatic sound speed and $G$ the gravitational constant), while for larger masses, the thermal forces are only a fraction of the linear estimate that decays as $M^{-1}$. Our analysis confirms the possibility of negative friction (hence a propulsion) on sufficiently luminous, low-mass embryos embedded in protoplanetary discs. Finally, we give an approximate expression of the total force at low Mach number, valid both for sub-critical ($M<M_c$) and super-critical ($M>M_c$) perturbers.

As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.

Using linear perturbation theory, we evaluate the time-dependent force exerted on an eccentric and inclined low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity $\chi$. We assume the eccentricity and inclination to be small compared to the size of the thermal lobes $\lambda\sim(\chi/\Omega)^{1/2}$, itself generally much smaller than the scalelength of pressure $H$. When the planet is non-luminous, we find that its eccentricity and inclination are vigorously damped by the disc, over a timescale shorter by a factor $H/\lambda$ than the damping timescale in adiabatic discs. On the contrary, when the luminosity-to-mass ratio of the planet exceeds a threshold that depends on the disc's properties, its eccentricity and inclination undergo an exponential growth. In the limit of a large luminosity, the growth rate of the eccentricity is 2.5~times larger than that of the inclination, in agreement with previous numerical work. Depending on their luminosity, planetary embryos therefore exhibit much more diverse behaviours than the mild damping of eccentricity and inclination considered hitherto.

We investigate by means of numerical simulations with nested Cartesian meshes the force exerted on a massive and luminous perturber moving at constant speed through a homogeneous and opaque gas, taking into account thermal diffusion in the gas and the radiative feedback from the perturber. The force arising from the release of energy into the ambient medium by the perturber, or heating force, is directed along the direction of motion and induces an acceleration of the perturber. Its value is compared to analytic estimates in the low and high Mach number limits, and found to match those accurately. In addition, the drag exerted on a non-luminous perturber significantly departs from the adiabatic expression when thermal diffusion is taken into account. In the limit of a vanishing velocity, this drag tends to a finite value which we determine using linear perturbation theory and corroborate with numerical simulations. The drag on a non-luminous perturber in a non-adiabatic gas therefore behaves like dry or solid friction. We work out the luminosity threshold to get a net acceleration of the perturber and find it to be generally much smaller than the luminosity of accreting low-mass planetary embryos embedded in a gaseous protoplanetary disc at a few astronomical units. We also present in some detail our implementation of nested meshes, which runs in parallel over several \emph{Graphics Processing Units} (GPUs).

We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour found using the FARGO3D code with a polar mesh. We make use of the time dependence of the torque exerted by the disc on the planet as a mean to quantify the numerical viscosity of the code. We find that the numerical viscosity of the Keplerian flow can be as low as a few $10^{-8}r^2\Omega$, $r$ and $\Omega$ being respectively the local orbital radius and frequency, for fifth order schemes and resolution of $\sim 10^{-2}r$. Although for a single disc problem a solution of low numerical viscosity can be obtained at lower computational cost with FARGO3D (which is nearly an order of magnitude faster than a fifth order method), discontinuous Galerkin methods appear promising to obtain solutions of low numerical viscosity in more complex situations where the flow cannot be captured on a polar or spherical mesh concentric with the disc.

Using linear perturbation theory, we investigate the torque exerted on a low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity. When the planet does not release energy into the ambient disc, the main effect of thermal diffusion is the softening of the enthalpy peak near the planet, which results in the appearance of two cold and dense lobes on either side of the orbit, of size smaller than the thickness of the disc. The lobes exert torques of opposite sign on the planet, each comparable in magnitude to the one-sided Lindblad torque. When the planet is offset from corotation, the lobes are asymmetric and the planet experiences a net torque, the `cold' thermal torque, which has a magnitude that depends on the relative value of the distance to corotation to the size of the lobes $\sim\sqrt{\chi/\Omega_p}$, $\chi$ being the thermal diffusivity and $\Omega_p$ the orbital frequency. We believe that this effect corresponds to the phenomenon named `cold finger' recently reported in numerical simulations, and we argue that it constitutes the dominant mode of migration of sub-Earth-mass objects. When the planet is luminous, the heat released into the ambient disc results in an additional disturbance that takes the form of hot, low-density lobes. They give a torque, named heating torque in previous work, that has an expression similar, but of opposite sign, to the cold thermal torque.

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited for non-magnetic problems. Benchmarked results of (hydrodynamical) planet-disk interaction and the of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficiently in straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community.

The migration of planets on nearly circular, non-inclined orbits in protoplanetary discs is entirely described by the disc's torque. This torque is a complex function of the disc parameters, and essentially amounts to the sum of two components: the Lindblad torque and the corotation torque. Known torque formulae do not reproduce accurately the torque actually experienced in numerical simulations by low- and intermediate- mass planets in radiative discs. One of the main reasons for this inaccuracy is that these formulae have been worked out in two-dimensional analyses. Here we revisit the torque formula and update many of its dimensionless coefficients by means of tailored, three- dimensional numerical simulations. In particular, we derive the dependence of the Lindblad torque on the temperature gradient, the dependence of the corotation torque on the radial entropy gradient (and work out a suitable expression of this gradient in a three-dimensional disc). We also work out the dependence of the corotation torque on the radial temperature gradient, overlooked so far. Corotation torques are known to scale very steeply with the width of the horseshoe region. We extend the expression of this width to the domain of intermediate mass planets, so that our updated torque formula remains valid for planets up to typically several tens of Earth masses, provided these relatively massive planets do not significantly deplete their coorbital region. Our torque expression can be applied to low- and intermediate-mass planets in optically thick protoplanetary discs, as well as protomoons embedded in circumplanetary discs.

We study the evolution of the eccentricity and inclination of protoplanetary embryos and low-mass protoplanets (from a fraction of an Earth mass to a few Earth masses) embedded in a protoplanetary disc, by means of three dimensional hydrodynamics calculations with radiative transfer in the diffusion limit. When the protoplanets radiate in the surrounding disc the energy released by the accretion of solids, their eccentricity and inclination experience a growth toward values which depend on the luminosity to mass ratio of the planet, which are comparable to the disc's aspect ratio and which are reached over timescales of a few thousand years. This growth is triggered by the appearance of a hot, under-dense region in the vicinity of the planet. The growth rate of the eccentricity is typically three times larger than that of the inclination. In long term calculations, we find that the excitation of eccentricity and the excitation of inclination are not independent. In the particular case in which a planet has initially a very small eccentricity and inclination, the eccentricity largely overruns the inclination. When the eccentricity reaches its asymptotic value, the growth of inclination is quenched, yielding an eccentric orbit with a very low inclination. As a side result, we find that the eccentricity and inclination of non-luminous planets are damped more vigorously in radiative discs than in isothermal discs.

We study gap formation in gaseous protoplanetary discs by a Jupiter mass planet. The planet's orbit is circular and inclined relative to the midplane of the disc. We use the impulse approximation to estimate the gravitational tidal torque between the planet and the disc, and infer the gap profile. For low-mass discs, we provide a criterion for gap opening when the orbital inclination is $\leq 30^{\circ}$. Using the FARGO3D code, we simulate the disc response to an inclined massive planet. The dependence of the depth and width of the gap obtained in the simulations on the inclination of the planet is broadly consistent with the scaling laws derived in the impulse approximation. Although we mainly focus on planets kept on fixed orbits, the formalism permits to infer the temporal evolution of the gap profile in cases where the inclination of the planet changes with time. This study may be useful to understand the migration of massive planets on inclined orbit, because the strength of the interaction with the disc depends on whether a gap is opened or not.

Planet migration is inherently a three-dimensional (3D) problem, because Earth-size planetary cores are deeply embedded in protoplanetary disks. Simulations of these 3D disks remain challenging due to the steep requirement in resolution. Using two different hydrodynamics code, FARGO3D and PEnGUIn, we simulate disk-planet interaction for a 1 to 5 Earth-mass planet embedded in an isentropic disk. We measure the torque on the planet and ensure that the measurements are converged both in resolution and between the two codes. We find that the torque is independent of the smoothing length of the planet's potential ($r_{\rm s}$), and that it has a weak dependence on the adiabatic index of the gaseous disk ($\gamma$). The torque values correspond to an inward migration rate qualitatively similar to previous linear calculations. We perform additional simulations with explicit radiative transfer using FARGOCA, and again find agreement between 3D simulations and existing torque formulae. We also present the flow pattern around the planets, and show that active flow is present within the planet's Hill sphere, and meridional vortices are shed downstream. The vertical flow speed near the planet is faster for a smaller $r_{\rm s}$ or $\gamma$, up to supersonic speeds for the smallest $r_{\rm s}$ and $\gamma$ in our study.

We consider the gravitational force exerted on a point-like perturber of mass $M$ travelling within a uniform gaseous, opaque medium at constant velocity $V$. The perturber irradiates the surrounding gas with luminosity $L$. The diffusion of the heat released is modelled with a uniform thermal diffusivity $\chi$. Using linear perturbation theory, we show that the force exerted by the perturbed gas on the perturber differs from the force without radiation (or standard dynamical friction). Hot, underdense gas trails the mass, which gives rise to a new force component, the heating force, with direction $+V$, thus opposed to the standard dynamical friction. In the limit of low Mach numbers, the heating force has expression $F_\mathrm{heat}=\gamma(\gamma-1)GML/(2\chi c_s^2)$, $c_s$ being the sound speed and $\gamma$ the ratio of specific heats. In the limit of large Mach numbers, $F_\mathrm{heat}=(\gamma-1)GML/(\chi V^2)f(r_\mathrm{min}V/4\chi)$, where $f$ is a function that diverges logarithmically as $r_\mathrm{min}$ tends to zero. Remarkably, the force in the low Mach number limit does not depend on the velocity. The equilibrium speed, when it exists, is set by the cancellation of the standard dynamical friction and heating force. In the low Mach number limit, it scales with the luminosity to mass ratio of the perturber. Using the above results suggests that Mars- to Earth-sized planetary embryos heated by accretion in a gaseous protoplanetary disc should have eccentricities and inclinations that amount to a sizeable fraction of the disc's aspect ratio, for conditions thought to prevail at a few astronomical units.

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.

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

We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.

We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on both "Graphical Processing Units" (GPUs) or "Central Processing unit" (CPUs), achieving large speed up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for the CPU, and have them translated automatically for the GPU.

We study the horseshoe dynamics of a low-mass planet in a three-dimensional, globally isothermal, inviscid disk. We find, as reported in previous work, that the boundaries of the horseshoe region (separatrix sheets) have cylindrical symmetry about the disk's rotation axis. We interpret this feature as arising from the fact that the whole separatrix sheets have a unique value of Bernoulli's constant, and that this constant does not depend on altitude, but only on the cylindrical radius, in barotropic disks. We next derive an expression for the torque exerted by the horseshoe region onto the planet, or horseshoe drag. Potential vorticity is not materially conserved as in two-dimensional flows, but it obeys a slightly more general conservation law (Ertel's theorem) which allows to obtain an expression for the horseshoe drag identical to the expression in a two-dimensional disk. Our results are illustrated and validated by three-dimensional numerical simulations. The horseshoe region is found to be slightly more narrow than previously extrapolated from two-dimensional analyses with a suitable softening length of the potential. We discuss the implications of our results for the saturation of the corotation torque, and the possible connection to the flow at the Bondi scale, which the present analysis does not resolve.

Planetary systems are born in the disks of gas, dust and rocky fragments that surround newly formed stars. Solid content assembles into ever-larger rocky fragments that eventually become planetary embryos. These then continue their growth by accreting leftover material in the disc. Concurrently, tidal effects in the disc cause a radial drift in the embryo orbits, a process known as migration. Fast inward migration is predicted by theory for embryos smaller than three to five Earth masses. With only inward migration, these embryos can only rarely become giant planets located at Earth's distance from the Sun and beyond, in contrast with observations. Here we report that asymmetries in the temperature rise associated with accreting infalling material produce a force (which gives rise to an effect that we call "heating torque") that counteracts inward migration. This provides a channel for the formation of giant planets and also explains the strong planet-metallicity correlation found between the incidence of giant planets and the heavy-element abundance of the host stars.

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analyzed. We also review its potential for public outreach in science museums.

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

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

Embedded planets disturb the density structure of the ambient disk and gravitational back-reaction will induce possibly a change in the planet's orbital elements. The accurate determination of the forces acting on the planet requires careful numerical analysis. Recently, the validity of the often used fast orbital advection algorithm (FARGO) has been put into question, and special numerical resolution and stability requirements have been suggested. In this paper we study the process of planet-disk interaction for small mass planets of a few Earth masses, and reanalyze the numerical requirements to obtain converged and stable results. One focus lies on the applicability of the FARGO-algorithm. Additionally, we study the difference of two and three-dimensional simulations, compare global with local setups, as well as isothermal and adiabatic conditions. We study the influence of the planet on the disk through two- and three-dimensional hydrodynamical simulations. To strengthen our conclusions we perform a detailed numerical comparison where several upwind and Riemann-solver based codes are used with and without the FARGO-algorithm. With respect to the wake structure and the torque density acting on the planet we demonstrate that the FARGO-algorithm yields correct results, and that at a fraction of the regular cpu-time. We find that the resolution requirements for achieving convergent results in unshocked regions are rather modest and depend on the pressure scale height of the disk. By comparing the torque densities of 2D and 3D simulations we show that a suitable vertical averaging procedure for the force gives an excellent agreement between the two. We show that isothermal and adiabatic runs can differ considerably, even for adiabatic indices very close to unity.

We give an expression for the Lindblad torque acting on a low-mass planet embedded in a protoplanetary disk that is valid even at locations where the surface density or temperature profile cannot be approximated by a power law, such as an opacity transition. At such locations, the Lindblad torque is known to suffer strong deviation from its standard value, with potentially important implications for type I migration, but the full treatment of the tidal interaction is cumbersome and not well suited to models of planetary population synthesis. The expression that we propose retains the simplicity of the standard Lindblad torque formula and gives results that accurately reproduce those of numerical simulations, even at locations where the disk temperature undergoes abrupt changes. Our study is conducted by means of customized numerical simulations in the low-mass regime, in locally isothermal disks, and compared to linear torque estimates obtained by summing fully analytic torque estimates at each Lindblad resonance. The functional dependence of our modified Lindblad torque expression is suggested by an estimate of the shift of the Lindblad resonances that mostly contribute to the torque, in a disk with sharp gradients of temperature or surface density, while the numerical coefficients of the new terms are adjusted to seek agreement with numerics. As side results, we find that the vortensity related corotation torque undergoes a boost at an opacity transition that can counteract migration, and we find evidence from numerical simulations that the linear corotation torque has a non-negligible dependency upon the temperature gradient, in a locally isothermal disk.

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).

Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be highly non spherical due to a hydrodynamical instability of the shocked accretion flow. Here we report the first laboratory experiment of a shallow water analogue, based on the physics of hydraulic jumps. Both the experiment and its shallow water modeling demonstrate a robust linear instability and nonlinear properties of symmetry breaking, in a system which is one million times smaller and about hundred times slower than its astrophysical analogue.

Based on the viscous disk theory, a number of recent studies have suggested there is large scale meridional circulation in protoplanetary disks. Such a flow could account for the presence of crystalline silicates, including calcium- and aluminum-rich inclusions (CAIs), at large distances from the sun. This paper aims at examining whether such large-scale flows exist in turbulent protoplanetary disks. High-resolution global hydrodynamical and magnetohydrodynamical (MHD) numerical simulations of turbulent protoplanetary disks were used to infer the properties of the flow in such disks. By performing hydrodynamic simulations using explicit viscosity, we demonstrate that our numerical setup does not suffer from any numerical artifact. The aforementioned meridional circulation is easily recovered in viscous and laminar disks and is quickly established. In MHD simulations, the magnetorotational instability drives turbulence in the disks. Averaging out the turbulent fluctuations on a long timescale, the results fail to show any large-scale meridional circulation. A detailed analysis of the simulations show that this lack of meridional circulation is due to the turbulent stress tensor having a vertical profile different from the viscous stress tensor. A simple model is provided that successfully accounts for the structure of the flow in the bulk of the disk. In addition to those results, possible deviations from standard vertically averaged alpha disk models are suggested by the simulations and should be the focus of future work. Global MHD numerical simulations of fully ionized and turbulent protoplanetary disks are not consistent with the existence of a large-scale meridional flow. As a consequence, the presence of crystalline silicates at large distance for the central star cannot be accounted for by that process as suggested by recent models based on viscous disk theory.

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.

The lower limit to the distribution of orbital periods P for the current population of close-in exoplanets shows a distinctive discontinuity located at approximately one Jovian mass. Most smaller planets have orbital periods longer than P~2.5 days, while higher masses are found down to P~1 day. We analyze whether this observed mass-period distribution could be explained in terms of the combined effects of stellar tides and the interactions of planets with an inner cavity in the gaseous disk. We performed a series of hydrodynamical simulations of the evolution of single-planet systems in a gaseous disk with an inner cavity mimicking the inner boundary of the disk. The subsequent tidal evolution is analyzed assuming that orbital eccentricities are small and stellar tides are dominant. We find that most of the close-in exoplanet population is consistent with an inner edge of the protoplanetary disk being located at approximately P>2 days for solar-type stars, in addition to orbital decay having been caused by stellar tides with a specific tidal parameter on the order of Q'*=10^7. The data is broadly consistent with planets more massive than one Jupiter mass undergoing type II migration, crossing the gap, and finally halting at the interior 2/1 mean-motion resonance with the disk edge. Smaller planets do not open a gap in the disk and remain trapped in the cavity edge. CoRoT-7b appears detached from the remaining exoplanet population, apparently requiring additional evolutionary effects to explain its current mass and semimajor axis.

We provide torque formulae for low mass planets undergoing type I migration in gaseous disks. These torque formulae put special emphasis on the horseshoe drag, which is prone to saturation: the asymptotic value reached by the horseshoe drag depends on a balance between coorbital dynamics (which tends to cancel out or saturate the torque) and diffusive processes (which tend to restore the unperturbed disk profiles, thereby desaturating the torque). We entertain here the question of this asymptotic value, and we derive torque formulae which give the total torque as a function of the disk's viscosity and thermal diffusivity. The horseshoe drag features two components: one which scales with the vortensity gradient, and one which scales with the entropy gradient, and which constitutes the most promising candidate for halting inward type I migration. Our analysis, which is complemented by numerical simulations, recovers characteristics already noted by numericists, namely that the viscous timescale across the horseshoe region must be shorter than the libration time in order to avoid saturation, and that, provided this condition is satisfied, the entropy related part of the horseshoe drag remains large if the thermal timescale is shorter than the libration time. Side results include a study of the Lindblad torque as a function of thermal diffusivity, and a contribution to the corotation torque arising from vortensity viscously created at the contact discontinuities that appear at the horseshoe separatrices. For the convenience of the reader mostly interested in the torque formulae, section 8 is self-contained.

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

We evaluate the horseshoe drag exerted on a low-mass planet embedded in a gaseous disk, assuming the disk's flow in the coorbital region to be adiabatic. We restrict this analysis to the case of a planet on a circular orbit, and we assume a steady flow in the corotating frame. We also assume that the corotational flow upstream of the U-turns is unperturbed, so that we discard saturation effects. In addition to the classical expression for the horseshoe drag in barotropic disks, which features the vortensity gradient across corotation, we find an additional term which scales with the entropy gradient, and whose amplitude depends on the perturbed pressure at the stagnation point of the horseshoe separatrices. This additional torque is exerted by evanescent waves launched at the horseshoe separatrices, as a consequence of an asymmetry of the horseshoe region. It has a steep dependence on the potential's softening length, suggesting that the effect can be extremely strong in the three dimensional case. We describe the main properties of the coorbital region (the production of vortensity during the U-turns, the appearance of vorticity sheets at the downstream separatrices, and the pressure response), and we give torque expressions suitable to this regime of migration. Side results include a weak, negative feed back on migration, due to the dependence of the location of the stagnation point on the migration rate, and a mild enhancement of the vortensity related torque at large entropy gradient.

We investigate the unsaturated horseshoe drag exerted on a low-mass planet by an isothermal gaseous disk. In the globally isothermal case, we use a formal- ism, based on the use of a Bernoulli invariant, that takes into account pressure effects, and that extends the torque estimate to a region wider than the horse- shoe region. We find a result that is strictly identical to the standard horseshoe drag. This shows that the horseshoe drag accounts for the torque of the whole corotation region, and not only of the horseshoe region, thereby deserving to be called corotation torque. We find that evanescent waves launched downstream of the horseshoe U-turns by the perturbations of vortensity exert a feed-back on the upstream region, that render the horseshoe region asymmetric. This asymmetry scales with the vortensity gradient and with the disk's aspect ratio. It does not depend on the planetary mass, and it does not have any impact on the horseshoe drag. Since the horseshoe drag has a steep dependence on the width of the horseshoe region, we provide an adequate definition of the width that needs to be used in horseshoe drag estimates. We then consider the case of locally isothermal disks, in which the tempera- ture is constant in time but depends on the distance to the star. The horseshoe drag appears to be different from the case of a globally isothermal disk. The difference, which is due to the driving of vortensity in the vicinity of the planet, is intimately linked to the topology of the flow. We provide a descriptive inter- pretation of these effects, as well as a crude estimate of the dependency of the excess on the temperature gradient.

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.

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

We 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.

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

The increase of computational resources has recently allowed high resolution, three dimensional calculations of planets embedded in gaseous protoplanetary disks. They provide estimates of the planet migration timescale that can be compared to analytical predictions. While these predictions can result in extremely short migration timescales for cores of a few Earth masses, recent numerical calculations have given an unexpected outcome: the torque acting on planets with masses between 5 M_Earth and 20 M_Earth is considerably smaller than the analytic, linear estimate. These findings motivated the present work, which investigates existence and origin of this discrepancy or ``offset'', as we shall call it, by means of two and three dimensional numerical calculations. We show that the offset is indeed physical and arises from the coorbital corotation torque, since (i) it scales with the disk vortensity gradient, (ii) its asymptotic value depends on the disk viscosity, (iii) it is associated to an excess of the horseshoe zone width. We show that the offset corresponds to the onset of non-linearities of the flow around the planet, which alter the streamline topology as the planet mass increases: at low mass the flow non-linearities are confined to the planet's Bondi sphere whereas at larger mass the streamlines display a classical picture reminiscent of the restricted three body problem, with a prograde circumplanetary disk inside a ``Roche lobe''. This behavior is of particular importance for the sub-critical solid cores (M <~ 15 M_Earth) in thin (H/r <~0.06) protoplanetary disks. Their migration could be significantly slowed down, or reversed, in disks with shallow surface density profiles.

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.

The discovery of close orbiting extrasolar giant planets led to extensive studies of disk planet interactions and the forms of migration that can result as a means of accounting for their location. Early work established the type I and type II migration regimes for low mass embedded planets and high mass gap forming planets respectively. While providing an attractive means of accounting for close orbiting planets intially formed at several AU, inward migration times for objects in the earth mass range were found to be disturbingly short, making the survival of giant planet cores an issue. Recent progress in this area has come from the application of modern numerical techniques which make use of up to date supercomputer resources. These have enabled higher resolution studies of the regions close to the planet and the initiation of studies of planets interacting with disks undergoing MHD turbulence. This work has led to indications of how the inward migration of low to intermediate mass planets could be slowed down or reversed. In addition, the possibility of a new very fast type III migration regime, that can be directed inwards or outwards, that is relevant to partial gap forming planets in massive disks has been investigated.

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

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

Among those objects formed out of collisional debris during galaxy mergers, the prominent gaseous accumulations observed near the tip of some long tidal tails are the most likely to survive long enough to form genuine recycled galaxies. Using simple numerical models, Bournaud, Duc & Masset (2003) claimed that tidal objects as massive as 10^9 Msun could only form, in these simulations, within extended dark matter (DM) haloes. We present here a new set of simulations of galaxy collisions to further investigate the structure of tidal tails. First of all, we checked that massive objects are still produced in full N-body codes that include feedback and a large number of particles. Using a simpler N-body code with rigid haloes, we noticed that dissipation and self-gravity in the tails, although important, are not the key factors. Exploiting toy models, we found that, for truncated DM haloes, material is stretched along the tail, while, within extended haloes, the tidal field can efficiently carry away from the disk a large fraction of the gas, while maintaining its surface density to a high value. This creates a density enhancement near the tip of the tail. Only later-on, self-gravity takes over; the gas clouds collapse and start forming stars. Thus, such objects were fundamentally formed following a kinematical process, according to a top-down scenario, contrary to the less massive Super Star Clusters that are also present around mergers. This conclusion leads us to introduce a restrictive definition for Tidal Dwarf Galaxies (TDGs) and their progenitors, considering only the most massive ones, initially mostly made of gas, that were able to pile up in the tidal tails. More simulations will be necessary to precisely determine the fate of these proto--TDGs and estimate their number.

[Abridged] The torque exerted by an external potential on a two-dimensional gaseous disk at non-co-orbital corotation resonances is studied by means of numerical simulations. The degree of saturation of these resonances is important in determining whether an eccentric giant planet embedded in a protoplanetary disk experiences an eccentricity excitation or damping. We perform calculations restricted to one or two resonances to investigate the properties of two neighboring corotation resonances, as well as the properties of a corotation resonance that overlaps a Lindblad resonance. We find that these properties hardly differ from the case of an isolated corotation resonance. In particular, although the torque of two neighboring corotation resonances may differ from the sum of the torques of the corresponding resonances considered as isolated, it never exceeds the sum of the fully unsaturated isolated corotation resonances, and it saturates in a fashion similar to an isolated resonance. Similarly, the presence of an underlying Lindblad resonance hardly affects the corotation torque, even if that resonance implies a torque strong enough to significantly redistribute the azimuthally averaged surface density profile, in which case the corotation torque scales with the resulting vortensity gradient. This set of numerical experiments thus essentially validates previous analytic studies. As a side result, we show that corotation libration islands misrepresented by a mesh of too low resolution can lead to a strongly overestimated corotation torque. This may be an explanation why the eccentricity of embedded Jupiter-sized planets was never observed to undergo an excitation in the numerical simulations performed so far.

Most Tidal Dwarf Galaxies (TDGs) so-far discussed in the literature may be considered as young ones or even newborns, as they are still physically linked to their parent galaxies by an umbilical cord: the tidal tail at the tip of which they are usually observed. Old Tidal Dwarf Galaxies, completely detached from their progenitors, are still to be found. Using N--body numerical simulations, we have shown that tidal objects as massive as 10^9 solar masses may be formed in interacting systems and survive for more than one Gyr. Old TDGs should hence exist in the Universe. They may be identified looking at a peculiarity of their "genetic identity card": a relatively high abundance in heavy elements, inherited from their parent galaxies. Finally, using this technique, we revisit the dwarf galaxies in the local Universe trying to find arguments pro and con a tidal origin.

In several interacting systems, gas accumulations as massive as 10^9 solar masses are observed near the tip of tidal tails, and are thought to be possible progenitors of Tidal Dwarf Galaxies. Using N-body simulations of galaxy interactions, we show that the existence of such features requires that dark matter haloes around spiral galaxies extend at least ten times further than the stellar disks. The massive gas clouds formed in our simulations have a kinematical origin and gravitationally collapse into dwarf galaxies that often survive for a few billion years.

[Abridged] We evaluate the coorbital corotation torque on a migrating protoplanet. The coorbital torque is assumed to come from orbit crossing fluid elements which exchange angular momentum with the planet when they execute a U-turn at the end of horseshoe streamlines. When the planet migrates inward, the fluid elements of the inner disk undergo one such exchange as they pass to the outer disk. The angular momentum they gain is removed from the planet, and this corresponds to a negative contribution to the corotation torque, which scales with the drift rate. In addition, the material trapped in the coorbital region drifts radially with the planet giving a positive contribution to the corotation torque, which also scales with the drift rate. These two contributions do not cancel out if the coorbital region is depleted, in which case there is a net corotation torque which scales with the drift rate and the mass deficit in the coorbital region, and which has same sign as the drift rate. This leads to a positive feedback on the migrating planet. In particular, if the coorbital mass deficit is larger than the planet mass, the migration rate undergoes a runaway which can vary the protoplanet semi-major axis by 50% over a few tens of orbits. This can happen only if the planet mass is sufficient to create a dip or gap in its surrounding region, and if the surrounding disk mass is larger than the planet mass. This typically corresponds to planet masses in the sub-Saturnian to Jovian mass range embedded in massive protoplanetary disks. Runaway migration is a good candidate to account for the orbital characteristics of close orbiting giant planets, most of which have sub-Jovian masses. Further, we show that in the runaway regime, migration can be directed outwards.

NGC 2915 is a blue compact dwarf galaxy embedded in an extended, low surface brightness HI disk exhibiting a two-armed spiral structure and a central bar-like component. Commonly accepted mechanisms are unable to explain the existence of these patterns and Bureau et al. proposed disk dark matter (scaling with the HI distribution) or a rotating triaxial dark halo as alternative solutions. In an attempt to explore these mechanisms, hydrodynamical simulations were run for each case and compared to observations using customized column density and kinematic constraints. The spiral structure can be accounted for both by an unseen bar or triaxial halo, the former fitting the observations slightly better. However, the large bar mass or halo pattern frequency required make it unlikely that the spiral wave is driven by an external perturber. In particular, the spin parameter is much higher than predicted by current cold dark matter (CDM) structure formation scenarios. The massive disk models show that when the observed gas surface density is scaled up by a factor about 10, the disk develops a spiral structure resembling closely the observed one, in perturbed density as well as perturbed velocity. This is consistent with more limited studies in other galaxies and suggests that the disk of NGC 2915 contains much more mass than is visible, tightly linked to the neutral hydrogen. A classic (quasi-)spherical halo is nevertheless still required, as increasing the disk mass further to fit the circular velocity curve would make the disk violently unstable. Scaling the observed surface density profile by an order of magnitude brings the disk and halo masses to comparable values within the disk radius.

We evaluate the coorbital corotation torque on a planet on a fixed circular orbit embedded in a viscous protoplanetary disk, for the case of a steady flow in the planet frame. This torque can be evaluated just from the flow properties at the separatrix between the librating (horseshoe) and circulating streamlines. A stationary solution is searched for the flow in the librating region. When used to evaluate the torque exerted by the circulating material of the outer and inner disk on the trapped material of the librating region, this solution leads to an expression of the coorbital corotation torque in agreement with previous estimates. An analytical expression is given for the corotation torque as a function of viscosity. Lastly, we show that additional terms in the torque expression can play a crucial role. In particular, they introduce a coupling with the disk density profile perturbation (the `dip' which surrounds the planet) and add to the corotation torque a small, positive fraction of the one-sided Lindblad torque. As a consequence, the migration could well be directed outwards in very thin disks (aspect ratio smaller than a few percent). This 2D analysis is especially relevant for mildly embedded protoplanets (sub-Saturn sized objects).

The torque felt by a non-accreting protoplanet on a circular orbit embedded in a uniform surface density protoplanetary disk is analyzed by means of time-dependent numerical simulations. Varying the viscosity enables one to disentangle the Lindblad torque (which is independent of viscosity) from the corotation torque, which saturates at low viscosity and is unsaturated at high viscosity. The dependence of the corotation torque upon the viscosity and upon the width of the librating zone is compared with previous analytical expressions, and shown to be in agreement with those. The effect of the potential smoothing respectively on the Lindblad torque and on the corotation torque is investigated, and the question of whether 3D effects and their impact on the total torque sign and magnitude can be modeled by an adequate smoothing prescription in a 2D simulation is addressed. As a side result, this study shows that the total torque acting on a Neptune-sized protoplanet is positive in a sufficiently thin, viscous disk (H/r < 4%, \alpha > 10^-2), but the inward migration time of smaller bodies is still very short, making it unlikely that they reach the torque reversal mass before having migrated all the way to the central object.

We present new results related to the coupled evolution of a two giant planet system embedded in a protoplanetary disk, in which a Saturn mass protoplanet is trapped in an outer mean motion resonance with a Jupiter mass protoplanet. The gaps opened in the disk by the two planets overlap, therefore the two planet system exchanges angular momentum with the disk at the inner planet's Inner Lindblad Resonances (ILRs) and at the outer planet Outer Lindblad Resonances (OLRs). Since the torques are proportional to the square of the planet masses, and since the inner planet is about three times more massive than the outer one, the ILR torques are favored by a factor ~10 with respect to the one planet case. In the case presented here, this leads to a positive differential Lindblad torque and consequently an outwards migration. We briefly discuss the long-term behavior of the system, which could account for the high eccentricities of the extra-solar planets with semi-major axis a>0.2 AU.

We present a mechanism related to the migration of giant protoplanets embedded in a protoplanetary disc whereby a giant protoplanet is caught up, before having migrated all the way to the central star, by a lighter outer giant protoplanet. This outer protoplanet may get captured into the 2:3 resonance with the more massive one, in which case the gaps that the two planets open in the disc overlap. Two effects arise, namely a squared mass weighted torque imbalance and an increased mass flow through the overlapping gaps from the outer disc to the inner disc, which both play in favour of an outwards migration. Indeed under the conditions presented here, which describe the evolution of a pair of protoplanets respectively Jupiter and Saturn sized, the migration is reversed, while the planets semi-major axis ratio is constant and the eccentricities are confined to small values by the disc material. The long-term behaviour of the system is briefly discussed, and could account for the high eccentricities observed for the extrasolar planets with semi-major axis a>0.2 AU.

We present an efficient and simple modification of the standard transport algorithm used in explicit eulerian fixed polar grid codes, aimed at getting rid of the average azimuthal velocity when applying the Courant condition. This results in a much larger timestep then the usual procedure, and it is particularly well-suited to the description of a Keplerian disk where one is traditionally limited by the very demanding Courant condition on the fast orbital motion at the inner boundary. In this modified algorithm, the timestep is limited by the perturbed velocity and by the shear arising from the differential rotation. FARGO stands for ``Fast Advection in Rotating Gaseous Objects''. The speed-up resulting from the use of the FARGO algorithm is problem dependent. In the example presented here, which shows the evolution of a Jupiter sized protoplanet embedded in a minimum mass protoplanetary nebula, the FARGO algorithm is about an order of magnitude faster than a traditional transport scheme, with a much smaller numerical diffusivity.

We investigate the gravitational interaction of a Jovian mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions have been investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Nonlinear hydrodynamic simulations are employed using three independent numerical codes. A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near circular orbit after a characteristic viscous time scale of approximately 10,000 initial orbital periods. This was found to be independent of whether the protoplanet was allowed to accrete mass or not. Inward migration is helped through the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion onto the central star.Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is approximately 5 Jupiter masses, because of the reduced accretion rates obtained for planets of increasing mass. Assuming some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets (Marcy, Cochran, & Mayor 1999; Marcy & Butler 1998) as well as their inward orbital migration can be accounted for by consideration of disc--protoplanet interactions during the late stages of giant planet formation. Maximally accreting protoplanets reached about four Jovian masses on reaching the neighbourhood of the central star.

We propose a new dynamical mechanism to account for the bar off-centering observed in many barred galaxies. It is based on the presence of a m=1 density wave, non-linearly excited by the strong m=2 component due to the linearly unstable bar and by a m=3 mode. N-body simulations suggest that this mechanism is at work in the centers of galaxies, and provide a natural explanation by spontaneously producing bar off-centering which amounts to several tens of parsecs.

We present the results of two-dimensional numerical simulations of stellar galactic disks, aimed at studying the non-linear coupling between bar and spiral waves and modes, in disks with realistically peaked rotation profiles. The power spectrum analysis of the perturbed density in the disk, for azimuthal numbers ranging from m=0 to m=4, shows an unambiguous signature of non-linear coupling between the bar and spiral waves, or between spiral waves only, with a very sharp selection of the frequencies which optimize the coupling efficiency. It turns out that non-linear coupling can be quite efficient, and even more relevant than the Swing mechanism to account for the dynamics of the galaxy beyond the corotation of the bar. Non-linear coupling is also responsible for a number of other behaviors observed in our runs, such as harmonic or sub-harmonic excitation, and the excitation of m=1 spiral waves.

We present a new mechanism to explain the frequently observed and thus certainly permanent warping of spiral galaxies. We consider the possibility of non-linear coupling between the spiral wave of the galaxy and two warp waves, such that the former, which is linearly unstable and extracts energy and angular momentum from the inner regions of the galactic disk, can continuously feed the latter. We derive an expression for the coupling coefficient in the WKB approximation. We show that the coupling is too weak in the stellar disk, except at the Outer Lindblad Resonance where the spiral slows down and is efficiently coupled to warp waves. There, the spiral can be almost completely converted into ``transmitted'' warps, which we can observe in HI, and a ``reflected'' one, which we can observe as a corrugation. Our mechanism reproduces the observed amplitudes of the warp and of the corrugation, and might explain related phenomena such as the behavior of the line of nodes of the warp. Furthermore we show that the energy and momentum fluxes of observed spirals and warps are of the same order of magnitude, adding a strong point in favor of this model.

We show that the propagation of warps in gaseous disks can be strongly affected by compressional effects, when the thickness of the disk is taken into account. The physical reason is that, in realistic self-gravitating disks, the sound time through the disk is comparable with the rotation time; thus the vertical hydrostatic equilibrium cannot be maintained adiabatically as the wave propagates (an implicit hypothesis in the thin-disk approximation) and the disk cannot move up and down solidly to follow the warp perturbation. There results, together with the main vertical motion, a strong horizontal one which significantly modifies the dispersion relation. We then turn to the case of a disk composed of two fluids with different temperatures: this can correspond either to the combined motions of the gas and the stars in a galactic disk, or to their coupling with a flattened massive halo (assuming that, as for spiral waves, stars are conveniently represented by a fluid if one stays away from Lindblad resonances). We find, in addition to the usual warps, a short-wavelength wave which might explain the ``corrugations'' observed in many galactic disks. Finally, as a side result of this analysis, we discuss a possible weak amplification of m>1 warps.