A generalized Zel'dovich approximation to gravitational instability

The orbits of particles undergoing gravitational instability are parameterized by generalizing the Zel'dovich approximation to a series expansion of arbitrary accuracy in the nonlinear regime. The coefficients of this series are determined from an action principle, or, more generally, from moments of the equation of motion. It is shown that the series is more rapidly convergent than previous nonlinear approximations. The method therefore provides a practical means of determining particle orbits, even for highly nonlinear perturbations. As an alternative, we also outline how the nonlinear dynamics may be computed as a field theory in which the evolution of the density and the velocity is determined in fixed comoving Eulerian coordinates.