Eccentricity and inclination of massive planets inside low-density cavities: results of 3D simulations

We study the evolution of eccentricity and inclination of massive planets in low-density cavities of protoplanetary discs using three-dimensional (3D) simulations. When the planet's orbit is aligned with the equatorial plane of the disc, the eccentricity increases to high values of 0.7-0.9 due to the resonant interaction with the inner parts of the disc. For planets on inclined orbits, the eccentricity increases due to the Kozai-Lidov mechanism, where the disc acts as an external massive body, which perturbs the planet's orbit. At small inclination angles, ${\lesssim}30^\circ$, the resonant interaction with the inner disc strongly contributes to the eccentricity growth, while at larger angles, eccentricity growth is mainly due to the Kozai-Lidov mechanism. We conclude that planets inside low-density cavities tend to acquire high eccentricity if favourable conditions give sufficient time for growth. The final value of the planet's eccentricity after the disc dispersal depends on the planet's mass and the properties of the cavity and protoplanetary disc.