Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Monte Carlo from Scratch

From the lab to the computer

Monte Carlo simulation of radiation transport is a scientific instrument by its own right. Cultivating the skill and the art is no less than gaining mastery in instrumentation. Monte Carlo is never a plug-and-play software tool. Sadly, Monte Carlo is often abused as a black box. Throughout the years I have also been ‘offered’ co-authorship of to-be papers which were to become publishable if I ran Monte Carlo for them. No, Monte Carlo cannot make an otherwise unpublishable work publishable.

A piece of study begins with a study or experimental design, considering the choice of samples and the control of variables. Letting all variables lose is not going to lead us to anything conclusive.

Embarking on a Monte Carlo project, I would think of the experimental setup the similar way as I would do when setting up a physical experiment in the radiation laboratory. The difference is that now:

we need not worry about ALARA or ALARP, radiation exposure or radiation protection; we can do cleaner experiments, we shall be able to control controlled variables much better; we can do experiments which would otherwise be


Experiments


Irradiation condition


Emission and transmission


Shielding


A source can be very simple


A source can be very complex


Photon cross-sections


A simple estimation


FLUKA input decks


Geometry and materials


The input file


A simple simulation


100% escaped


Carbon instead of vacuum


80% transmission


More histories


70% transmission


Lower cutoffs


Carbon beam instead of carbon shell


Beam-start position


Pion+


Alpha


Helium-3


Deuterium


Tritium


Gamma


Electron


Positron


Proton


Proton


Proton


Counts of different particles


So many neutron


Pion-


Carbon-12


Nine protons


The 9 protons


2nd-generation protons


Two gammas


Pion+


The pion+


The other pion+


Particles no longer tracked


A history out of many


FLUKA


Fortran coding is optional


Housebrands