• B. Erdelyi
  Northern Illinois University, DeKalb, Illinois

When faced with the challenge of the design optimization of a charged particle beam system involving beam-material interactions, a framework is needed that seamlessly integrate the following tasks: 1) high order accurate and efficient beam optics, 2) a suite of codes that model the atomic and nuclear interactions between the beam and matter, and 3) the option to run many different optimization strategies at the code language level with a variety of user-defined objectives. To this end, we developed a framework in COSY Infinity with these characteristics and which can be run in two modes: map mode and a hybrid map-Monte Carlo mode. The code, its applications to the FRIB, and plans involving large-scale computing will be presented.

• T. Flisgen, G. Pöplau, U. van Rienen
  Rostock University, Faculty of Computer Science and Electrical Engineering, Rostock

Demanding applications such as heavy ion fusion, high energy colliders and free electron lasers require the study of beam phenomena like space-charge induced instabilities, emittance growth and halo formation. Numerical simulations for instance with GPT (General Particle Tracer, Pulsar Physics) calculate the mutual Coulomb interactions of the tracked particles *. The direct summation of the forces is rather costly and scales with O(N²). In this paper we investigate a new approach for the efficient calculation of particle-particle interactions: the fast summation by Nonequispaced Fast Fourier Transform (NFFT) **, whereas the NFFT is a generalization of the well known Fast Fourier Transformation (FFT). We describe the algorithm and discuss the performance and accuracy of this method for several particle distributions.

• L.L. Bandura
  ANL, Argonne
• B. Erdelyi
  Northern Illinois University, DeKalb, Illinois

While COSY INFINITY provides powerful DA methods for the simulation of fragment separator beam dynamics, the master version of COSY does not currently take into account beam-material interactions. These interactions are key for accurately simulating the dynamics from heavy ion fragmentation and fission. In order to model the interaction with materials such as the target or absorber, much code development was needed. There were four auxiliary codes implemented in COSY for the simulation of beam-material interactions. These include EPAX for returning the cross sections of isotopes produced by fragmentation and MCNPX for the cross sections of isotopes produced by the fission and fragmentation of a 238U beam. ATIMA is implemented to calculate energy loss and energy and angular straggling. GLOBAL returns the charge state. The extended version can be run in map mode or hybrid map-Monte Carlo mode, providing an integrated beam dynamics-nuclear processes design optimization and simulation framework that is efficient and accurate. The code, its applications, and plans for large-scale computational runs for optimization of separation purity of rare isotopes at FRIB will be presented.

• H. Okuno
  RIKEN Nishina Center, Wako
• A. Adelmann
  PSI, Villigen
• J.J. Yang
  CIAE, Beijing

Since 1997 RIKEN Nishina center has been constructing a next-generation exotic beam facility, RI beam factory (RIBF), based on a powerful heavy ion driver accelerator . Its accelerator complex was successfully commissioned at the end of 2006 and started supplying heavy ion beams in 2007. The four ring cyclotrons (RRC, fRC, IRC and SRC) connected in series accelerate the energy of the heavy ion beams up to 400 MeV/u for the lighter ions such as argon and 345 MeV/u for heavier ions such as uranium. Intensity upgrade plans are under way, including the construction of a new 28 GHz superconducting ECR ion source. The new ECR will take all the succeeding accelerators and beam transport lines to a space charge dominant regime, which should be carefully reconsidered to avoid emittance growth due to space charge forces. Beam dynamics in the low energy cyclotron, RRC was studied by OPAL-cycl a flavor of the OPAL. The simulation results clearly show vortex motions in the isochronous field, resulting in round beam formation in the first 10 turns after the injection point. The possible increase of beam loss at beam extraction will be also discussed in this paper.