Beam dynamics with the Hamilton-Jacobi equation

We describe a non-perturbative method to solve the Hamilton-Jacobi equation for invariant surfaces in phase space. The problem is formulated in action-angle variables with a general nonlinear perturbation. The solution of the Hamilton-Jacobi equation is regarded as the fixed point of a map on the Fourier coefficients of the generating function. Periodicity of the generator in the independent variable is enforced with a shooting method. We present two methods for finding the fixed point and hence the invariant surface. A solution by plain iteration is economical but has a restricted domain of convergence. The Newton iteration is costly but yields solutions up to the dynamic aperture. Examples of lattices with sextupoles for chromatic correction are discussed. 10 refs., 5 figs., 1 tab