Dynamical Models of the Solar System Formation and Evolution

Mark Lewis has extended his them is work by completing a series of N-body simulations of a narrow ring that: is in the location of Saturn's F-ring and is perturbed by a single satellite comparable to Prometheus, but on a circular orbit. We had previously shown how the satellite perturbations can cause a broadly distributed sparse population of ring particles to become concentrated into narrow rings that can be maintained outside of any resonance location. For low optical depths, the collisions between ring particles are highly localized in the peaks of the satellite wakes. The inelastic collisions therefore occur at a particular orbital phase angle so as to damp the azimuthal component of the relative velocities. Since particle positions are not changed by collisions, the semimajor axes of the particles are shifted toward the actual particle positions where the collisions occur. Thus, negative radial diffusion can occur while conserving orbital angular momentum so long as the forced eccentricity is continually re-excited by the satellite. We speculated that the separation between the final ringlets was largely determined by the magnitude of the forced eccentricities induced by the satellite at closest approach. We carried out a series of simulations with a variety of different satellite masses in order to vary the magnitude of the forced eccentricity. We found that indeed the final spacing of the ringlets does increase with the magnitude of the forced eccentricity (Lewis and Stewart 2002). This occurs because neighboring eccentric ringlets drift out of phase with one another due to Keplerian shear and eventually collide with one another, leading to a smaller number of more widely spaced ringlets, The time scale required to form narrow ringlets in these simulations is much shorter than one would expect from standard theories based upon the orbit-averaged torque produced by multiple passes by the satellite. We find that the initial ringlets form in less than two synodic periods and the final state is typically reached in 10 to 20 synodic periods. These studies move us closer to understanding the significantly more complex system of Saturn's F ring, where the perturbation magnitude varies over short temporal and spatial time scales due to the orbital eccentricities of the perturbing satellite. We are currently extending the simulation to allow for an eccentric orbit of the satellite.