Authors: Harold F. Levison, Martin J. Duncan, & Edward Thommes

Status: To appear in Astronomical Journal

Abstract: We presented the first particle based, Lagrangian code that can follow the collisional/accretional/dynamical evolution of a large number of km-sized planetesimals through the entire growth process to become planets. We refer to it as the Lagrangian Integrator for Planetary Accretion and Dynamics or LIPAD. LIPAD is built on top of SyMBA, which is a symplectic N-body integrator (Duncan, Levison & Lee 1998). In order to handle the very large number of planetesimals required by planet formation simulations, we introduce the concept of a tracer particle. Each tracer is intended to represent a large number of disk particles on roughly the same orbit and size as one another, and is characterized by three numbers: the physical radius, the bulk density, and the total mass of the disk particles represented by the tracer. We developed statistical algorithms that follow the dynamical and collisional evolution of the tracers due to the presence of one another. The tracers mainly dynamically interact with the larger objects (planetary embryos)} in the normal N-body way. LIPAD's greatest strength is that it can accurately model the wholesale redistribution of planetesimals due to gravitational interaction with the embryos, which has recently been shown to significantly affect the growth rate of planetary embryos. We verify the code via a comprehensive set of tests which compare our results with those of Eulerian and/or direct N-body codes.

Click here for pdf file this paper.
Click here for pdf with embedded animations. It will only work with acrobat reader.

Click here for the animation assoaciated with Figure 6.
Click here for the animation assoaciated with Figure 8.
Click here for the animation assoaciated with Figure 12.
Click here for the animation assoaciated with Figure 13.
Click here for the animation assoaciated with Figure 14.