“…The results of test simulations are found to converge when the number of macroparticles is more than 40000 and step size is less than 2 mm. We have performed the simulation for the following choice of parameters: total beam current I = 10 mA, mesh size N x = N y = 128, number of macroparticles for each species N = 77000, step size ∆s = 1 mm in the axial direction, percentage species fractions [p, H + 2 , H + 3 ] as [80, 15,5] and initial rms beam size 2016 JINST 11 P04020 X rms = 1.25 mm at s = 0 for all the species. The charge on the macroparticles of each species is determined from the species fraction and the total number of macroparticles.…”