Author: Berz, M.
Paper Title Page
TUSCI1
The Fast Multipole Method in the Differential Algebra Framework for the Computation of 3D Space Charge Fields and Transfer Maps  
 
  • H. Zhang, M. Berz
    MSU, East Lansing, Michigan, USA
 
  A new algorithm which combines the classical multiple level fast multipole algorithm with differential algebraic tools for three dimensional electrostatic field calculation of N-body systems is presented. The whole charged domain is hierarchically decomposed into boxes of varying sizes according to the local charge density. Each box has a near region, where the interaction is calculated directly, and its a region, where the interaction is represented by expansions in multipoles at infinity. The total field is the summation of the near region field and the far region field. This algorithm can treat any arbitrary charge distribution with an efficiency of O(N). Differential algebra simplifies the math especially for non-point charge macro particles, and allows the calculation of high order field derivatives, which is essential to generate the transfer map. Examples of applying the methods in beam dynamics simulations are given.  
slides icon Slides TUSCI1 [3.302 MB]  
 
TUADI1 Storage Ring EDM Simulation: Methods and Results 99
 
  • Y. Senichev, A. Lehrach, R. Maier, D. Zyuzin
    FZJ, Jülich, Germany
  • S.N. Andrianov, A.N. Ivanov
    St. Petersburg State University, St. Petersburg, Russia
  • M. Berz, K. Makino
    MSU, East Lansing, Michigan, USA
 
  The idea of Electric Dipole Moment search using the electrostatic storage ring with polarized beam is based on accumulation of additional tiny spin rotation, about one-billionth of radians per second, occurred only in the presence of EDM. This method can be realized under condition of the long-time spin coherency ~1000 seconds. During this time each particle performs about 109 turns in ring moving on different trajectories. At such conditions the spin-rotation aberrations associated with various types of space and time dependent nonlinearities start playing a crucial role. To design such a ring the computer simulation is necessary taking into account all factors affecting the spin. We used COSY-Infinity and integrating program with symplectic Runge-Kutta methods in composition with analytic methods. We developed a new lattice based on the alternating spin rotating. As a result, we can achieve the SCT of ~5000 seconds. The difficulties of these studies are still in the fact that the aberrations growth is observed in the scale of 109 turns and few million particles. For this simulation we use a supercomputer with parallel computing process.  
slides icon Slides TUADI1 [0.951 MB]  
 
THSCI1
Rigorous Fixed Point Enclosures and their Application to High-Order Transfer Maps  
 
  • A. Wittig, M. Berz
    MSU, East Lansing, Michigan, USA
 
  Funding: This work was supported in part by the US Department of Energy and the Studienstiftung des Deutschen Volkes.
A common task in the design and analysis of an accelerator is the study of the transfer map of the system. Of particular interest is the estimation of the region of stability of a given system. Typically, this is done using symplectic particle tracking and visual analysis of the resulting Poincare maps for signatures of chaoticity and island structures near high-period fixed points. We describe a method to compute rigorous enclosures of all periodic points of a given order in a given map based on Taylor Model methods. We then apply this algorithm to a real world transfer map of the Tevatron accelerator to rigorously identify islands and resonances in its transfer map. This mathematically rigorous method to locate resonances in the transfer map automatically yields all regions where resonances up to a certain order appear. The island structure exhibited by the map in those regions is then studied further by computing the invariant manifolds associated with the hyperbolic periodic points of the map. This manifold structure can provide further insight into the dynamics of the map, including the emergence of chaotic motion at the appearance of crossings of the manifolds.
 
slides icon Slides THSCI1 [3.230 MB]  
 
THADC2
Precision Nonlinear Dynamics in Electrostatic Rings for EDM Studies  
 
  • M. Berz, R. Jagasia, K. Makino
    MSU, East Lansing, Michigan, USA
 
  Funding: US Department of Energy
Recent interest electric dipole moments of hadrons and light nuclei necessitates the design of storage ring designs relying on electrostatic elements for deflection and focusing. Since the effects of interest are exceedingly small and difficult to extract from the dynamics, unusually high precision in the understanding of the linear and nonlinear optics is required. This leads to the need of precise treatment of a various effects that are often of subdued importance, including the detailed shape of fringe fields and their influence on steering, spin flip devices, and rf cavities out of synch with the frequency of the ring. In all studies, it is important to perform symplectic tracking for the full spin-orbit dynamics. At the same time it is necessary to preserve the particle’s energy during transversal of electric potentials, errors in which are similarly detrimental as errors in symplecticity. Unfortunately it is known that schemes to preserve symplecticity do not simultaneously preserve energy and vice versa. We present some hybrid approaches based on high-order maps that minimize the violations of either invariant while minimally impacting spin-orbit motion.
 
 
THADC3
Nonlinear, Nonscaling CW FFAG Design and Modeling Using Map Methods  
 
  • K. Makino, M. Berz
    MSU, East Lansing, Michigan, USA
  • C. Johnstone
    Fermilab, Batavia, USA
 
  Funding: US Department of Energy, FNAL, Particle Accelerator Corporation
FFAGs are being considered for various new classes of accelerators, including proton drivers for muon colliders and neutrino factories, for accelerator driven subcritical reactors, and for medical applications. Their key advantages are compactness, CW operation and large acceptance. However, due to the complicated field arrangements, beam dynamics simulations are challenging for conventional simulation codes. Among various levels of sophistication of dealing with beam elements, one of the most advanced modes of the code COSY INFINITY computes the 3D field distributions along the trajectory of the particles, not only at the point of interest but also in the neighborhood with functional dependencies, using DA PDE solvers based on Differential Algebraic techniques. The scheme allows the code developer and user to describe rather complicated field arrangements with rather limited effort. Of particular advantage is the seamless integration with map-based symplectic tracking schemes. The methods are illustrated in design studies of FFAGs for various different machine types, including both conventional strong focusing and continuously varying combined function setups.
 
slides icon Slides THADC3 [3.843 MB]