Abstract

High fidelity charged particle beam dynamics, especially if collisional, require specialized numerical methods. The challenges include the ability to deal with very large particle numbers, long-range electromagnetic forces, and vast spatial and timescales. We developed a novel algorithm to address these challenges. The main characteristics of the algorithm are Strang splitting to separate near and far forces, the fast multipole method to lower the computational cost of the far region, and the Sim\`o integrator to capture all close encounters efficiently in the near region. The algorithm is fully adaptive both in space and time, while maintaining symplecticity to machine precision. Its performance is illustrated with two challenging examples from nonlinear multiparticle beam dynamics, including the first electron cooling simulations based on first principles.

Highlights

  • Charged particle beams are ensembles of close-by electrons or ions in directed motion

  • The vast majority of these codes are based upon variants of what generically may be called mean-field, collisionless codes. While this is sufficient for many applications, several recent beam physics applications emphasized the need for new, collisional algorithms and codes based on them

  • We briefly summarized two such applications here: ultracold electron beam generation based on arrays of sharp nanotips and electron cooling of ion beams

Read more

Summary

INTRODUCTION

Charged particle beams are ensembles of close-by electrons or ions in directed motion. Based on a theorem of Simò [22] and aided by differential algebraic methods [23], we previously developed a collisional numerical integrator, the Simò integrator, with some provably optimal behavior [24] This integrator is variable order and adaptive with automatic selection of particle-by-particle optimal order and time step to minimize computational cost given an a priori user-selected tolerance for error. The resulting algorithm is an accurate and efficient collisional simulation of charged particle beams in external electromagnetic fields that is symplectic to machine precision and is provably the best in reaching an a priori set error level with minimized computational cost.

EQUATIONS OF MOTION AND SPLITTING METHOD
FAST MULTIPOLE METHOD AND THE LONG-RANGE INTERACTIONS
THE SIMÒ INTEGRATOR AND THE SHORT-RANGE INTERACTIONS
IMPLEMENTATION AND PARALLELIZATION OF THE ALGORITHM
Splitting reconciliation
Parallelization of the Simò integrator
Parallelization of PHAD
Dynamics of beamlets
Comparison of PHAD and the standalone Simò integrator simulations
PHAD simulations of beamlets dynamics
Electron cooling of ion beams
Density modulations for the coherent electron cooling
Conventional electron cooling of ion beams
Findings
SUMMARY AND CONCLUSIONS
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.