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.

Analysing the pattern of energy depositions may help elucidate differences in the severity of radiation-induced DNA strand breakage for different radiation qualities. It is often claimed that energy deposition (ED) sites from photon radiation form a uniform random pattern, but there is indication of differences in RBE values among different photon sources used in brachytherapy. The aim of this work is to analyse the spatial patterns of EDs from 103Pd, 125I, 192Ir, 137Cs sources commonly used in brachytherapy and a 60Co source as a reference radiation. The results suggest that there is both a non-uniform and a uniform random component to the frequency distribution of distances to the nearest neighbour ED. The closest neighbouring EDs show high spatial correlation for all investigated radiation qualities, whilst the uniform random component dominates for neighbours with longer distances for the three higher mean photon energy sources (192Ir, 137Cs, and 60Co). The two lower energy photon emitters (103Pd and 125I) present a very small uniform random component. The ratio of frequencies of clusters with respect to 60Co differs up to 15% for the lower energy sources and less than 2% for the higher energy sources when the maximum distance between each pair of EDs is 2 nm. At distances relevant to DNA damage, cluster patterns can be differentiated between the lower and higher energy sources. This may be part of the explanation to the reported difference in RBE values with initial DSB yields as an endpoint for these brachytherapy sources.

BackgroundThe globally dominant treatment with curative intent for locally advanced esophageal squamous cell carcinoma (ESCC) is neoadjuvant chemoradiotherapy (nCRT) with subsequent esophagectomy. This multimodal treatment leads to around 60% overall 5-year survival, yet with impaired post-surgical quality of life. Observational studies indicate that curatively intended chemoradiotherapy, so-called definitive chemoradiotherapy (dCRT) followed by surveillance of the primary tumor site and regional lymph node stations and surgery only when needed to ensure local tumor control, may lead to similar survival as nCRT with surgery, but with considerably less impairment of quality of life. This trial aims to demonstrate that dCRT, with selectively performed salvage esophagectomy only when needed to achieve locoregional tumor control, is non-inferior regarding overall survival, and superior regarding health-related quality of life (HRQOL), compared to nCRT followed by mandatory surgery, in patients with operable, locally advanced ESCC.MethodsThis is a pragmatic open-label, randomized controlled phase III, multicenter trial with non-inferiority design with regard to the primary endpoint overall survival and a superiority hypothesis for the experimental intervention dCRT with regard to the main secondary endpoint global HRQOL one year after randomization. The control intervention is nCRT followed by preplanned surgery and the experimental intervention is dCRT followed by surveillance and salvage esophagectomy only when needed to secure local tumor control. A target sample size of 1200 randomized patients is planned in order to reach 462 events (deaths) during follow-up.Clinical Trial Registrationwww.ClinicalTrials.gov, identifier: NCT04460352.

The linear-quadratic (LQ) parameterization of survival fraction [SF( D)] inherently assumes that all cells in a population receive the same dose ( D), albeit the distribution of specific energy z over the individual cells f( z, D) can be very wide. From these microdosimetric distributions, which are target size dependent, we estimate the size of the cellular sensitive volume by analyzing its influence on the LQ parameterization of cell survival. A Monte Carlo track structure code was used to simulate detailed tracks from a Co source as well as proton and carbon ions of various energies. From these tracks, f( z, D) distributions were calculated for spherical targets with diameters ranging from 10 nm to 12 μm. A cell survival function based on f( z, D) was fitted to experimental LQ α values, revealing an intrinsic limitation that target size imposes on the usage of f( z, D) to describe the linear term of the LQ parameterization. The results indicate that such threshold volume arises naturally from the relationship between the particle's probability of no-hit and the probability of cell survival. Further analysis led to the proposal of a radiobiological property [Formula: see text], defined as the mean lineal energy corresponding to the target size that allows equivalence between the mean inactivation dose (MID) and the mean specific energy [Formula: see text]. The fact that [Formula: see text] is an increasing continuous function of target size within the range of biological targets of interest in radiobiology, ensures the uniqueness of [Formula: see text] for any radiation quality, thus, its potential usefulness in modeling. In conclusion, an accurate estimation of such threshold volumes may be useful for improving modeling of cell survival curves.