Planetesimal-driven planet migration in the presence of a gas disk

We report here on an extension of a previous study by Kirsh et al. (Kirsh, D.R., Duncan, M., Brasser, R., Levison, H.F. [2009]. Icarus 199, 197–209) of planetesimal-driven migration using our N-body code SyMBA (Duncan, M.J., Levison, H.F., Lee, M.H. [1998]. Astron. J. 116, 2067–2077). The previous work focused on the case of a single planet of mass Mem, immersed in a planetesimal disk with a power-law surface density distribution and Rayleigh distributed eccentricities and inclinations. Typically 104–105 equal-mass planetesimals were used, where the gravitational force (and the back-reaction) on each planetesimal by the Sun and planet were included, while planetesimal–planetesimal interactions were neglected. The runs reported on here incorporate the dynamical effects of a gas disk, where the Adachi et al. (Adachi, I., Hayashi, C., Nakazawa, K. [1976]. Prog. Theor. Phys. 56, 1756–1771) prescription of aerodynamic gas drag is implemented for all bodies. In some cases the Papaloizou and Larwood (Papaloizou, J.C.B., Larwood, J.D. [2000]. Mon. Not. R. Astron. Soc. 315, 823–833) prescription of Type-I migration for the planet are implemented, as well as a mass distribution.