View on GitHub

GEARS | 齿轮组

Geant4 Example Application with Rich features and Small footprints

Download this project as a .zip file Download this project as a tar.gz file

Stopping power Bragg curve Absorption range Physics list

Alpha interactions

Alpha particles, or ions in general, are heavy (compared to elemental particles) electrically charged particles. According to Knoll’s Radiation Detection and Measurement, they interact with matter primarily through coulomb forces between their positive charge and the negative charge of the orbital electrons within the absorber atoms. Although interactions of the particle with nuclei (as in Rutherford scattering or alpha-particle-induced reactions) are also possible, such encounters occur only rarely and they are not normally significant in the response of radiation detectors.

Common interactions that can be assigned to an alpha in a physics list include:

Stopping power

The stopping power S = -dE/dx is simply the energy loss in a unit track length. The GEARS output provides two variables to roughly calculate this: de and dl. de is the energy deposited in a step point. dl is the corresponding step length.

// draw dE/dx VS kinetic energy of alpha particles
root[] t->Draw("de/dl:k", "pdg==22")

An alternative simulation tool dedicated to ion interaction in matter is SRIM. It may produce better results than Geant4 in some cases.

The NIST stopping-power and range tables for helium ions is a great resource to verify your Geant4 simulation of alpha particles.

Bragg curve

The energy of an alpha particle decreases as it moves deeper and deeper into certain material. As it slows down, the surrounding electrons have more time to interact with it. Consequently, its dE/dx goes up. However, when it becomes really slow, it can pick up electrons and get neutralized. Consequently, its dE/dx quickly drops to zero. Such behavior can be shown clearly by the Bragg curve:

Bragg curve simulated Bragg curve

Such a curve can be produced by

$ root air.root
root[] t->Draw("de/dl:l","trk==1 && de>15","colz")

air.root can be created by

$ gears BraggCurve.mac

In BraggCurve.mac, 5.5 MeV alphas are shot to air. The simulated world is defined in air.tg

Range

The spatial distributions of MeV alphas in air can be draw with the following commands:

$ root air.root
// draw tracks of the first 100 simulated alpha particles
root[] t->Draw("y:x","trk==1","l",100,0)
// display end points of alpha tracks
root[] t->Draw("y:z","trk==1")

alpha in air end points of alpha tracks

One can see that there are often a group of hits very close to each other around the end points of alpha tracks. This is due to the large dE/dx at low energies.

Impact of physics list

In range.mac, 5.5 MeV alphas are shot to a CsI scintillation crystal instead of air. The simulated world is defined in CsI.tg. Without specifying a reference physics list (that is, to use the default one specified by GEARS, QGSP_BERT), all alpha particles deposit its full energy in just one step. The step lengths are all about 28 um. To do a detailed simulation, one needs to specify the maximal step length of alphas in CsI. This can be done manually using the macro command /run/setCut 1 um or to select a more suitable physics list:

# LBE: low background experiment (maximal step length: 1 um)
# LIV: Livermore data based EM models
# comment out the following line and run again to see the difference
/physics_lists/select LBE_LIV

Run gears range.mac with and without this line and draw the alpha hits distributions to see the difference.