No Arabic abstract
We study the stability of gaps opened by a giant planet in a self-gravitating protoplanetary disc. We find a linear instability associated with both the self-gravity of the disc and local vortensity maxima which coincide with gap edges. For our models, these edge modes develop and extend to twice the orbital radius of a Saturn mass planet in discs with disc-to-star mass ratio >0.06, corresponding to a Toomre Q < 1.5 at the outer disc boundary. Unlike the local vortex-forming instabilities associated with gap edges in weakly or non-self-gravitating low viscosity discs, the edge modes are global and exist only in sufficiently massive discs, but for the typical viscosity values adopted for protoplanetary discs. Analytic modelling and linear calculations show edge modes may be interpreted as a localised disturbance associated with a gap edge inducing activity in the extended disc, through the launching of density waves excited at Lindblad resonances. Nonlinear hydrodynamic simulations are performed to investigate the evolution of edge modes in disc-planet systems. The form and growth rates of unstable modes are consistent with linear theory. Their dependence on viscosity and gravitational softening is also explored. We also performed a first study of the effect of edge modes on planetary migration. We found that if edge modes develop, then the average disc-on-planet torque becomes more positive with increasing disc mass. In simulations where the planet was allowed to migrate, although a fast type III migration could be seen that was similar to that seen in non-self-gravitating discs, we found that it was possible for the planet to interact gravitationally with the spiral arms associated with an edge mode and that this could result in the planet being scattered outwards. Thus orbital migration is likely to be complex and non monotonic in massive discs of the type we consider.
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.
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 planets 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 carry out three-dimensional SPH simulations to study whether planets can survive in self-gravitating protoplanetary discs. The discs modelled here use a cooling prescription that mimics a real disc which is only gravitationally unstable in the outer regions. We do this by modelling the cooling using a simplified method such that the cooling time in the outer parts of the disc is shorter than in the inner regions, as expected in real discs. We find that both giant (> M_Sat) and low mass (< M_Nep) planets initially migrate inwards very rapidly, but are able to slow down in the inner gravitationally stable regions of the disc without needing to open up a gap. This is in contrast to previous studies where the cooling was modelled in a more simplified manner where regardless of mass, the planets were unable to slow down their inward migration. This shows the important effect the thermodynamics has on planet migration. In a broader context, these results show that planets that form in the early stages of the discs evolution, when they are still quite massive and self-gravitating, can survive.
The gravitational interaction between a protoplanetary disc and planetary sized bodies that form within it leads to the exchange of angular momentum, resulting in migration of the planets and possible gap formation in the disc for more massive planets. In this article, we review the basic theory of disc-planet interactions, and discuss the results of recent numerical simulations of planets embedded in protoplanetary discs. We consider the migration of low mass planets and recent developments in our understanding of so-called type I migration when a fuller treatment of the disc thermodynamics is included. We discuss the runaway migration of intermediate mass planets (so-called type III migration), and the migration of giant planets (type II migration) and the associated gap formation in the disc. The availability of high performance computing facilities has enabled global simulations of magnetised, turbulent discs to be computed, and we discuss recent results for both low and high mass planets embedded in such discs.
It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragments interaction with the disc, but with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment-fragment interactions while the disc is still present. We find that fragment-fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations - this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment-fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low.