A 3D Monte Carlo code for the modeling of plasma dynamics and beam formation mechanism in electron cyclotron resonance ion sources

The code here presented is the first part of a Monte Carlo (MC) self-consistent 3D plasma simulator. It is yet able to solve the equation of motion for thousands of independent charged particles. The procedure allows to understand the consequences of each phenomenon introduced in the evolution steps of the code. MC random selection of starting parameters is used for each particles; the environmental conditions enclosed in the simulation are ECRIS magnetic field, resonant electromagnetic wave, initial plasma density distribution and MC calculation of Spitzer collision. The results of the first simulations explain some typical effects as the hollow beam formation and the main plasma deconfinement mechanism.