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.

We investigate the response of an accretion disk to the presence of a perturbing protoplanet embedded in the disk through time dependent hydrodynamical simulations. The disk is treated as a two-dimensional viscous fluid and the planet is kept on a fixed orbit. We run a set of simulations varying the planet mass, and the viscosity and temperature of the disk. All runs are followed until they reach a quasi-equilibrium state. We find that for planetary masses above a certain minimum mass, already 3 M_Jup for a viscosity of nu = 10^{-5}, the disk makes a transition from a nearly circular state into an eccentric state. Increasing the planetary mass leads to a saturation of disk eccentricity with a maximum value of around 0.25. The transition to the eccentric state is driven by the excitation of an m=2 spiral wave at the outer 1:3 Lindblad resonance. The effect occurs only if the planetary masses are large enough to clear a sufficiently wide and deep gap to reduce the damping effect of the outer 1:2 Lindblad resonance. An increase in viscosity and temperature in the disk, which both tend to close the gap, have an adverse influence on the disk eccentricity. In the eccentric state the mass accretion rate onto the planet is greatly enhanced, an effect that may ease the formation of massive planets beyond about 5 M_Jup that are otherwise difficult to reach.