Abstract: A translation of the penelope physics subroutines to C++, designed as an extension of the Geant4 toolkit, is presented. The Fortran code system penelope performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. Penelope implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of Geant4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original penelope programs. The combined code, named PenG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of Geant4.

Abstract: Purpose: In skin high-dose-rate (HDR) brachytherapy, sources are located outside, in contact with, or implanted at some depth below the skin surface. Most treatment planning systems use the TG-43 formalism, which is based on single-source dose superposition within an infinite water medium without accounting for the true geometry in which conditions for scattered radiation are altered by the presence of air. The purpose of this study is to evaluate the dosimetric limitations of the TG-43 formalism in HDR skin brachytherapy and the potential clinical impact. Methods: Dose rate distributions of typical configurations used in skin brachytherapy were obtained: a 5 cm x 5 cm superficial mould; a source inside a catheter located at the skin surface with and without backscatter bolus; and a typical interstitial implant consisting of an HDR source in a catheter located at a depth of 0.5 cm. Commercially available HDR Co-60 and Ir-192 sources and a hypothetical Yb-169 source were considered. The Geant4Monte Carlo radiation transport code was used to estimate dose rate distributions for the configurations considered. These results were then compared to those obtained with the TG-43 dose calculation formalism. In particular, the influence of adding bolus material over the implant was studied. Results: For a 5 cm x 5 cm Ir-192 superficial mould and 0.5 cm prescription depth, dose differences in comparison to the TG-43 method were about -3%. When the source was positioned at the skin surface, dose differences were smaller than -1% for Co-60 and Ir-192, yet -3% for Yb-169. For the interstitial implant, dose differences at the skin surface were -7% for Co-60, -0.6% for Ir-192, and -2.5% for Yb-169. Conclusions: This study indicates the following: (i) for the superficial mould, no bolus is needed; (ii) when the source is in contact with the skin surface, no bolus is needed for either Co-60 and Ir-192. For lower energy radionuclides like Yb-169, bolus may be needed; and (iii) for the interstitial case, at least a 0.1 cm bolus is advised for Co-60 to avoid underdosing superficial target layers. For Ir-192 and Yb-169, no bolus is needed. For those cases where no bolus is needed, its use might be detrimental as the lack of radiation scatter may be beneficial to the patient, although the 2% tolerance for dose calculation accuracy recommended in the AAPM TG-56 report is not fulfilled.

Abstract: A two-body fully relativistic three-dimensional scattering framework has been utilised to develop an alternative Compton scattering computational model to those adapted from Ribberfors' work for Monte Carlo modelling of Compton scattering. Using a theoretical foundation that ensures the conservation of energy and momentum in the relativistic impulse approximation, this new model, the Monash University Compton scattering model, develops energy and directional algorithms for both the scattered photon and ejected Compton electron from first principles. The Monash University Compton scattering model was developed to address the limitation of the Compton electron directionality algorithms of other computational models adapted from Ribberfors' work. Here the development of the Monash University Compton scattering model, including its implementation in a Geant4 low energy electromagnetic physics class, G4LowEPComptonModel, is outlined. Assessment of the performance of G4LowEPComptonModel was undertaken in two steps: (1) comparison with respect to the two standard Compton scattering classes of Geant4 version 9.5, G4LivermoreComptonModel and G4PenelopeComptonModel, and (2) experimental comparison with respect to Compton electron kinetic energy spectra obtained from the Compton scattering of 662 key photons off the K-shell of gold. Both studies illustrate that the Monash University Compton scattering model, and in turn G4LowEPComptonModel, is a viable replacement for the majority of computational models that have been adapted from Ribberfors' work. It was also shown that the Monash University Compton scattering model is able to reproduce the Compton scattering triply differential cross-section Compton electron kinetic energy spectra of 662 keV photons K-shell scattering off of gold to within experimental uncertainty.