Jeremy Emmett edited Analytical Estimates.tex  over 9 years ago

Commit id: f72485a716ab0406764c21a7242e93b0b22b05ef

deletions | additions      

       

Migration occurs because of an interchange of angular momentum between disk particles and planets, such that the total angular momentum of the system remains constant. Depending on the relative angular momenta of these components, a planet can either undergo an outward or inward migration. Ida et al. (2000b) and Gomes (2004) offer an analytical explanation of this process, based on Rutherford Scattering. That is, during an encounter between two particles, the velocity vectors of the two particles instantaneously change. If the planet is placed on a circular orbit, then particles with a z-component of angular momentum larger than that of the planet will drive an outward migration. This is because, when crossing near the planet, their tangential velocities are larger than the circular velocity of the planet, thus accelerating it into a higher orbit. Particles with a z-component of angular momentum smaller than that of the planet, however, will drive an inward migration of the planet. This is because, when crossing near the planet, their tangential velocities are smaller than the circular velocity of the planet, thus negatively accelerating it into a lower orbit.  Though the disk has some 'thickness' due to the fact that particles were initially given non-zero z-components in their positions, and therefore begin on trajectories that are inclined, we will assume a simple model in which all particles follow perfectly circular orbits in the plane of the system ($e,i = 0$)  In this case, all angular momentum is in the z-component, with a magnitude of $a\times m\times v$ for each particle.  In this particular simulation, the mass of Neptune was 4.5e-5, and the disk particle mass was approximately $3.33\times 10^{-6}$ (about 7.4\% the mass of Neptune). All bodies were initially placed on circular orbits ranging between a =0 and a = 1, with Neptune at a radius of a=0.5. Therefore, particles which were initially in the vicinity of proximity to  Neptune had a angular momentum of approximately $L=0.5\times 3.33\times 10^{-6}*v_n$ whereas Neptune had an angular momentum of approximately $L=0.5\times 4.5\times 10^{-5}*v_n$ These differ only by a factor of the difference in their masses, so all particles which were initially near Neptune had an angular momentum smaller than that of Neptune. Therefore, we would expect to have observed an initial inward migration of Neptune. In order for Neptune to reverse this trend, and eventually undergo an outward migration, the angular momentum (and thus the orbital velocity) of the disk particles in the vicinity of Neptune would have to increase substantially. In other words, For particles near Neptune, we would have to see, through some process, a depletion in the number of those with angular momenta less than Neptune's and an increase in the number of those with angular momenta greater than Neptune's.