We present the results of both global cylindrical disc simulations and local shearing box simulations of protoplanets interacting with a disc undergoing MHD turbulence with zero net flux magnetic fields. We investigate the nature of the disc response and conditions for gap formation. This issue is an important one for determining the type and nature of the migration of the protoplanet, with the presence of a deep gap being believed to enable slower migration. For both types of simulation we find a common pattern of behaviour for which the main parameter determining the nature of the response is $M_p R^3/(M_* H^3)$, with $M_p$, $M_*$, $R$, and $H$ being the protoplanet mass, the central mass, the orbital radius and the disc semi-thickness respectively. We find that as this parameter is increased towards 0.1, the presence of the protoplanet is first indicated by the appearance of the well known trailing wake which, although it may appear erratic on account of the turbulence, appears to be well defined. Once the above parameter exceeds a number around unity a gap starts to develop inside which the magnetic energy density tends to be concentrated in the high density wakes. This gap formation condition can be understood from simple dimensional considerations of the conditions for nonlinearity, and the balance of angular momentum transport due to Maxwell and Reynolds' stresses with that due to tidal torques. An important result is that the basic flow morphology in the vicinity of the protoplanet is very similar in both the local and global simulations. This indicates that local shearing box simulations, which are computationally less demanding, capture much of the physics of disc-planet interaction. Thus they may provide a useful tool for studying the local interaction between forming protoplanets and turbulent, protostellar discs.

We consider two protoplanets gravitationally interacting with each other and a protoplanetary disc. The two planets orbit interior to a tidally maintained disc cavity while the disc interaction indices inward migration. When the migration is slow enough, the more rapidly migrating outer protoplanet approaches and becomes locked in a 2:1 commensurability with the inner one. This is maintained in subsequent evolution. We study this evolution using a simple anaytic model, full hydrodynamic 2D simulations of the disc planet system and longer time N body integrations incorporating simple prescriptions for the effect of the disc on the planet orbits. The eccentricity of the protoplanets are found to be determined by the migration rate induced in the outer planet orbit by the external disc. We apply our results to the recently discovered resonant planets around GJ876. Simulation shows that a disc with parameters expected for protoplanetary discs causes trapping in the 2:1 commensurability when the planets orbit in an inner cavity and that eccentricities in the observed range may be obtained.

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