Soft-edged magnet models for higher-order beam-optics map codes

Continuously varying surface and volume source-density distributions are used to model magnetic fields inside of cylindrical volumes. From these distributions, a package of subroutines computes on-axis generalized gradients and their derivatives at arbitrary points on the magnet axis for input to the numerical map-generating subroutines of the Lie-algebraic map code Marylie. In the present version of the package, the magnet menu includes: (1) cylindrical current-sheet or radially thick current distributions with either open boundaries or with a surrounding cylindrical boundary with normal field lines (which models high-permeability iron), (2) Halbach-type permanent multipole magnets, either as sheet magnets or as radially thick magnets, (3) modeling of arbitrary fields inside a cylinder by use of a fictitious current sheet. The subroutines provide on-axis gradients and their z derivatives to essentially arbitrary order, although in the present third- and fifth-order Marylie only the zeroth through sixth derivatives are needed. The formalism is especially useful in beam-optics applications, such as magnetic lenses, where realistic treatment of fringe-field effects is needed.