Publication Date


Document Type


First Advisor

Erdelyi, Bela

Degree Name

Ph.D. (Doctor of Philosophy)

Legacy Department

Department of Physics


Accurate and efficient computational methods are essential to study the dynamics of charged particle beams. Most of the available numerical methods in beam physics are based on the collisionless model which is sufficient for applications that does not depend on Coulomb collisions of the simultaneously interacting particles. However, the development of collisional numerical methods has become greatly important in recent years for applications where collisions play an important role such as the electron cooling of high energy hadron beams. While collisional methods can provide very accurate simulations, their algorithms are very complex, and they face efficiency challenges. In this work, we present our development of new collisional numerical methods which are the first methods that can provide an accurate microscopic description of beam dynamics with high computational efficiency. The first method is the Simo` integrator which solves the N-body problem of beams by direct integration of the particles' equations of motion in the presence of external electromagnetic fields. Its development included very unique techniques to obtain accuracy while resolving all the efficiency challenges known to N-body integrators. Consequently, the Simo` integrator is the first large-scale collisional numerical method in beam physics that is accurate up to machine precision with a relatively high efficiency. Then, we incorporate the Simo` integrator to model collisions into our other collisional method referred to as the Particles' High-Order Adaptive Dynamics (PHAD). PHAD employs an advanced version of the fast multipole method (FMM) along with Strang splitting method, and the addition of the Simo` integrator makes PHAD the first most efficient, numerically symplectic, collisional method in beam physics. For an enhanced performance, the algorithms of both the Simo` integrator and PHAD were fully parallelized on a large-scale high-performance hybrid cluster. We present simulations performed by our codes of three complicated beam dynamics problems. One application is for the electron cooling of ion beams to which our simulations demonstrate and give the first insight of the microscopic description of electron cooling with accurate prediction of cooling time. The other application illustrates density modulations of electron beams due to ions from a collisional picture of the dynamics and provide conditions to obtain a strong modulation signal necessary for variants of coherent electron cooling systems. The last application considers microscopic simulations of the relaxation of certain beam perturbations which illustrates finite N effects in contrast to the kinetic limit of the collisionless methods, and the resulted relaxation times are important for applications like the beam echo.


190 pages




Northern Illinois University

Rights Statement

In Copyright

Rights Statement 2

NIU theses are protected by copyright. They may be viewed from Huskie Commons for any purpose, but reproduction or distribution in any format is prohibited without the written permission of the authors.

Media Type


Included in

Physics Commons