Collision probabilities in the presence of nebular gas drag

We are developing a model to determine what fraction of the planetesimals would have hit a protoplanet on their sunward journey as opposed to having a close approach and passing into an inferior orbit. The model involves direct numerical integration of restricted-three-body orbits using a predictor-corrector integrator. A simple gas drag law with a v(exp 2) dependence was also included in the equations of motion. Runs of 100 to 500 particles were already performed, while some future runs may require several times this number in order to get good impact statistics. All planetesimals start in superior orbits with semi-major axes 5 to 10 R(sub H) from the protoplanets, where R(sub H) is the protoplanet's Hill Sphere radius. The orbit is followed until the planetesimal passed into an inferior orbit at least 10 R(sub H) from the protoplanet. This process typically requires 10(exp 4) to 10(exp 5) orbits.