• R. D. Ryne, E. W. Bethel, I. V. Pogorelov, J. Qiang, J. M. Shalf, C. Siegerist, M. Venturini
  LBNL, Berkeley, California
• D. T. Abell
  Tech-X, Boulder, Colorado
• A. Adelmann
  PSI, Villigen
• J. F. Amundson, P. Spentzouris
  Fermilab, Batavia, Illinois
• A. Dragt
  University of Maryland, College Park, Maryland
• C. Mottershead, N. Neri, P. L. Walstrom
  LANL, Los Alamos, New Mexico
• V. Samulyak
  BNL, Upton, Long Island, New York

MaryLie/IMPACT (ML/I) is a 3D parallel Particle-In-Cell code that combines the nonlinear optics capabilities of MaryLie 5.0 with the parallel particle-in-cell space-charge capability of IMPACT. In addition to combining the capabilities of these codes, ML/I has a number of powerful features, including a choice of Poisson solvers, a fifth-order rf cavity model, multiple reference particles for rf cavities, a library of soft-edge magnet models, representation of magnet systems in terms of coil stacks with possibly overlapping fields, and wakefield effects. The code allows for map production, map analysis, particle tracking, and 3D envelope tracking, all within a single, coherent user environment. ML/I has a front end that can read both MaryLie input and MAD lattice descriptions. The code can model beams with or without acceleration, and with or without space charge. Developed under a US DOE Scientific Discovery through Advanced Computing (SciDAC) project, ML/I is well suited to large-scale modeling, simulations having been performed with up to 100M macroparticles. ML/I uses the H5Part* library for parallel I/O. The code inherits the powerful fitting/optimizing capabilities of MaryLie, augmented for the new features of ML/I. The combination of soft-edge magnet models, high-order capability, and fitting/optimization, makes it possible to simultaneously remove third-order aberrations while minimizing fifth-order, in systems with overlapping, realistic magnetic fields. Several applications will be presented, including aberration correction in a magnetic lens for radiography, linac and beamline simulations of an e-cooling system for RHIC, design of a matching section across the transition of a superconducting linac, and space-charge tracking in the damping rings of the International Linear Collider.

*ICAP 2006 paper ID 1222, A. Adelmann et al., "H5Part: A Portable High Performance Parallel Data Interface for Electromagnetics Simulations"