The stochastic nature of ionizing radiation interactions causes a microdosimetric spread in energy depositions for cell or cell nucleus-sized volumes. The magnitude of the spread may be a confounding factor in dose response analysis. The aim of this work is to give values for the microdosimetric spread for a range of doses imparted by (125)I and (192)Ir brachytherapy radionuclides, and for a (60)Co source. An upgraded version of the Monte Carlo code PENELOPE was used to obtain frequency distributions of specific energy for each of these radiation qualities and for four different cell nucleus-sized volumes. The results demonstrate that the magnitude of the microdosimetric spread increases when the target size decreases or when the energy of the radiation quality is reduced. Frequency distributions calculated according to the formalism of Kellerer and Chmelevsky using full convolution of the Monte Carlo calculated single track frequency distributions confirm that at doses exceeding 0.08 Gy for (125)I, 0.1 Gy for (192)Ir, and 0.2 Gy for (60)Co, the resulting distribution can be accurately approximated with a normal distribution. A parameterization of the width of the distribution as a function of dose and target volume of interest is presented as a convenient form for the use in response modelling or similar contexts.

The MC tool can predict SSB and DSB yields for light ions of various LET and estimate RBEDSB (direct). In addition, it can calculate the frequencies of different DNA lesion sizes, which is of interest in the context of biologically relevant absolute dosimetry of particle beams.

The f cluster characterization of ionizing radiation at a nanometer scale can effectively be used to calculate particle and energy dependent α and β values to predict RBE values with potential applications to, e.g., treatment planning systems in radiotherapy.

A discrepancy between the Monte Carlo derived relative standard deviation (microdosimetric spread) and experimental data was reported by Villegas et al (2013 Phys. Med. Biol. 58 6149–62) suggesting wall effects as a plausible explanation. The comment by Lindborg et al (2015 Phys. Med. Biol. 60 8621–4) concludes that this is not a likely explanation. A thorough investigation of the Monte Carlo (MC) transport code used for track simulation revealed a critical bug. The corrected MC version yielded values that are now within experimental uncertainty. Other microdosimetric findings are hereby communicated.

The magnitude of the spread in specific energy deposition per cell may be a confounding factor in dose-response analysis motivating derivation of explicit data for the most common brachytherapy isotopes (125)I and (192)Ir, and for (60)Co radiation frequently used as reference in RBE studies. The aim of this study is to analyse the microdosimetric spread as given by the frequency distribution of specific energy for a range of doses imparted by (125)I, (192)Ir and (60)Co sources. An upgraded version of the Monte Carlo code PENELOPE was used for scoring energy deposition distributions in liquid water for each of the radiation qualities. Frequency distributions of specific energy were calculated according to the formalism of Kellerer and Chmelevsky. Results indicate that the magnitude of the microdosimetric spread increases with decreasing target size and decreasing energy of the radiation quality. Within the clinical relevant dose range (1 to 100 Gy), the spread does not exceed 4 % for (60)Co, 5 % for (192)Ir and 6 % for (125)I. The frequency distributions can be accurately approximated with symmetrical normal distributions at doses down to 0.2 Gy for (60)Co, 0.1 Gy for (192)Ir and 0.08 Gy for (125)I.