Franck – Hertz experiment: Monte Carlo simulation

In this section you can find a simulation of the scattering of electrons in Neon gas. Electrons are accelerated by an electric potential (increasing in time from 0 V to 80 V) in this region they perform some scattering events. Then they reach a repulsive region and only the most energetic particles can overcome it. In the end the simultaion results are multiplied by V3/2 to take into account the space charge effect (Child-Langmuir law).

Show the code

Here a comparison of the experiment with the result of the simulation:

Some info about the code:

  • Source is planar, negative value of Z are forbidden.
  • Scattering angle is isotropic – means cosθ = 2*R-1  (R: random number from 0 to 1).
  • Random freepath d = -1/Σ*log(1-R)  (Σ: macroscopic cross section).
  • I have considered a polinomial interpolation of the following cross sections: elastic, excitation(16.86 eV – 18.38 eV), ionization (lxcat). with some corrections to take into account the other excitations.
  • The instrumantation have an additional potential called “contact potential”, in the code is “-U_cont”.

If you need some kind of introduction to monte carlo procedure you can take a look to these other pages: thermal neutron diffusion, neutron slowing down, iterated source method, Irradiation damage.