In order to efficiently perform particle simulations in systems with widely varying magnetization, we developed a drift-Lorentz mover, which interpolates between full particle dynamics and drift kinetics in such a way as to preserve a physically correct gyroradius and particle drifts for both large and small ratios of the timestep to the cyclotron period. In order to extend applicability of the mover to systems with plasma frequency exceeding the cyclotron frequency such as one may have with fully neutralized drift compression of a heavy-ion beam we have developed an implicit version of the mover. A first step in this direction, in which the polarization charge was added to the field solver, was described previously. Here we describe a fully implicit algorithm (which is analogous to the direct-implicit method for conventional particle-in-cell simulation), summarize a stability analysis of it, and describe several tests of the resultant code.

We present initial results for the self-consistent beam-cloud dynamics simulations for a sample LHC beam, using a newly developed set of modeling capability based on a merge of the three-dimensional parallel Particle-In-Cell (PIC) accelerator code WARP and the electron-cloud code POSINST. Although the storage ring model we use as a test bed to contain the beam is much simpler and shorter than the LHC, its lattice elements are realistically modeled, as is the beam and the electron cloud dynamics. The simulated mechanisms for generation and absorption of the electrons at the walls are based on previously validated models available in POSINST.

In order to efficiently perform particle simulations in systems with widely varying magnetization, we developed a drift-Lorentz mover, which interpolates between full particle dynamics and drift kinetics in such a way as to preserve a physically correct gyroradius and particle drifts for both large and small ratios of the timestep to the cyclotron period. In order to extend applicability of the mover to systems with plasma frequency exceeding the cyclotron frequency such as one may have with fully neutralized drift compression of a heavy-ion beam we have developed an implicit version of the mover. A first step in this direction, in which the polarization charge was added to the field solver, was described previously. Here we describe a fully implicit algorithm (which is analogous to the direct-implicit method for conventional
particle-in-cell simulation), summarize a stability analysis of it, and describe several tests of the resultant code.