Comparison between EGSnrc, Geant4, MCNP5 and Penelope for mono-energetic electron beams

A simple geometry is chosen to highlight similarities and differences of current electron transport algorithms implemented in four Monte Carlo codes commonly used in radiation physics. Energy deposited in a water-filled sphere by mono-energetic electron beams was calculated using EGSnrc, Geant4, MCNP5 and Penelope as the radius of the sphere varied from 0.25 cm to 4.5 cm for beam energies of 0.5 MeV, 1.0 MeV and 5.0 MeV. The calculations were performed in single-scattering mode (where applicable) and in condensed history mode. A good agreement is found for the single-scattering calculations except for the in-air case at 0.5 MeV where differences increase with decreasing radius up to 5% between EGSnrc and Penelope. Differences between results calculated with the default user settings when compared to their own single-scattering modes are under 5% for all codes when the sphere is surrounded by vacuum, however, large differences occur for Geant4, MCNP5 and Penelope when air is introduced around the sphere. Finally, the parameters associated with the multiple scattering algorithms were tuned reducing these differences below 10% for these codes at the expense of increased computation time.