We investigate the stability of those low-mass Trojan planets that form in a protoplanetary disc and subsequently accrete gas to become gas giants. We calculate their evolution before, during, and after gas disc dispersal. A two-dimensional hydrodynamics code combined with an N-body solver is used to evolve the system of disc and planets. Gas disc dispersal is simulated in a simple manner by assuming global exponential decay of the disc mass, leading to the stalling of migration after semi-major axes have approximately halved from their initial values. We consider Trojan pairs with different initial masses and gas accretion rates and gas disc models with different masses and viscosities. An N-body code, adapted to model disc forces, is used to examine large-scale migration and the formation of very short period Trojan planets. For each combination of planetary pair and disc model that we consider in our hydrodynamic simulations, each Trojan system remains stable before, during, and after disc dispersal. The long-term stability of these systems in the absence of gas is tested using N-body simulations, and all systems remain stable for those evolution times equal to 10^9 years. Eccentricities remain low (e<0.02) in all cases. Increases in the amplitude of libration about the L4/L5 Lagrange points accompany the inward migration, and during very large-scale migration Trojan systems may be disrupted prior to the onset of disc dispersal. The stability of Trojan pairs during rapid type I migration, during the transition to type II migration with the accompanying gap formation in the gas disc, and during gas loss when the disc disperses, indicates that isolated Trojan planet systems are very stable. If a common mechanism exists for their formation, we suggest they may be readily observed in nature.

Protoplanet eccentricities of e >~ H/r can slow or reverse migration, but previous 2D studies have shown that gravitational scattering cannot maintain significant planet eccentricities against disc-induced damping. We simulate the evolution of low-mass protoplanetary swarms in three dimensions. The aim is to examine both protoplanet survival rates and the dynamical structure of the resulting planetary systems, and to compare them with 2D simulations. We present results from a 3D hydrodynamic simulation of eight protoplanets embedded in a protoplanetary disc. We also present a suite of simulations performed using an N-body code, modified to include prescriptions for planetary migration and for eccentricity and inclination damping. These prescriptions were obtained by fitting analytic formulae to hydrodynamic simulations of planets embedded in discs with initially eccentric and/or inclined orbits. As was found in two dimensions, differential migration produces groups of protoplanets in stable, multiple mean-motion resonances that migrate in lockstep, preventing prolonged periods of gravitational scattering. In almost all simulations, this leads to large-scale migration of the protoplanet swarm into the central star in the absence of a viable stopping mechanism. The evolution involves mutual collisions, occasional instances of large-scale scattering, and the frequent formation of the long-lived, co-orbital planet systems that arise in > 30% of all runs. Disc-induced damping overwhelms eccentricity and inclination growth due to planet-planet interactions. Co-orbital planets are a natural outcome of dynamical relaxation in a strongly dissipative environment, and if observed in nature would imply that such a period of evolution commonly arises during planetary formation.

Young planets embedded in their protoplanetary disk interact gravitationally with it leading to energy and angular momentum exchange. This interaction determines the evolution of the planet through changes to the orbital parameters. We investigate changes in the orbital elements of a 20 Earth--mass planet due to the torques from the disk. We focus on the non-linear evolution of initially non-vanishing eccentricity $e$ and/or inclination $i$. We treat the disk as a two- or three-dimensional viscous fluid and perform hydrodynamical simulations with an embedded planet. We find rapid exponential decay of the planet orbital eccentricity and inclination for small initial values of $e$ and $i$, in agreement with linear theory. For larger values of $e > 0.1$ the decay time increases and the decay rate scales as $\dot{e} \propto e^{-2}$, consistent with existing theoretical models. For large inclinations ($i$ > 6 deg) the inclination decay rate shows an identical scaling $di/dt \propto i^{-2}$. We find an interesting dependence of the migration on the eccentricity. In a disk with aspect ratio $H/r=0.05$ the migration rate is enhanced for small non-zero eccentricities ($e < 0.1$), while for larger values we see a significant reduction by a factor of $\sim 4$. We find no indication for a reversal of the migration for large $e$, although the torque experienced by the planet becomes positive when $e \simeq 0.3$. This inward migration is caused by the persisting energy loss of the planet. For non gap forming planets, eccentricity and inclination damping occurs on a time scale that is very much shorter than the migration time scale. The results of non linear hydrodynamic simulations are in very good agreement with linear theory for small $e$ and $i$.

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.