Eccentricity evolution during planet-disc interaction

During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.
18 pages, 13 figures, accepted for publication in MNRAS