MONTE CARLO SIMULATION OF ELECTRICAL DISCHARGES IN GASES

New Monte Carlo methods of calculation for tracing the motion of electrons in gases are described. Two problems have been treated by these methods, and the techniques applied to electron swarms in neon. Firstly, the probability of back-scatter of electrons to the cathode has been estimated for values of E/p0 between 10 and 200 V cm-1 torr-1 (E is the electric field and p0 the pressure reduced to 0°C) and for emission energies up to the lowest excitation potential. Secondly, Townsend primary ionization coefficients α/p0 and drift velocities have been evaluated for values of E/p0 from 20 to 400 V cm-1 torr-1. It is shown that the use of the ergodic hypothesis is not valid, yielding values of α/p0 which are up to 40% too high in neon, and that entire avalanches must be simulated before meaningful results can be obtained. The present results agree well with those obtained recently at Swansea using the Boltzmann equation, indicating that the Lorentz approximation employed in the latter method leads to only small errors in estimates of the mean properties of the electron swarms. The equilibrium of the swarm with the field is considered in detail. It is shown that, for uniform fields, deviations from exponential growth of current with increasing electrode separation occur if the interelectrode voltage is not at least three times the mean energy of electrons impinging on the anode.