The relative biological effectiveness for carbon, nitrogen, and oxygen ion beams using passive and scanning techniques evaluated with fully 3D silicon microdosimeters

The rationale for a carbon ion radiation therapy facility in Australia

Australia has taken a collaborative nationally networked approach to achieve particle therapy capability. This supports the under-construction proton therapy facility in Adelaide, other potential proton centres and an under-evaluation proposal for a hybrid carbon ion and proton centre in western Sydney. A wide-ranging overview is presented of the rationale for carbon ion radiation therapy, applying observations to the case for an Australian facility and to the clinical and research potential from such a national centre.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

LETd Optimization Verification With an SOI Microdosimeter

PURPOSE

A first of its kind experimental verification of dose-averaged linear energy transfer (LETd) optimized treatment plans for proton therapy has been carried out using a silicon-on-insulator microdosimeter at the Massachusetts General Hospital (MGH), Boston, USA.

METHODS AND MATERIALS

Three clinical treatment plans of a typical ependymoma structure set were designed using the standard clinical approach, the proposed protocol approach, and a one-field approach. The plans were then reoptimized to reduce the LETd-weighted dose in the brain stem. All six plans were delivered in a solid water phantom and the experimental yD‾ measured.

RESULTS

After LETd optimization, a reduction in yD‾ was found within the brain stem by an average of 12%, 19%, and 4% for the clinical, protocol, and one-field plans, respectively, while maintaining adequate coverage of the tumor structure. The experimental LETd-weighted doses were in agreement with the treatment planning system calculations and Monte Carlo simulations and reinforced the improvement of the optimization.

CONCLUSIONS

This work demonstrates the first experimental verification of the clinical implementation of LETd optimization for patient treatment with proton therapy.

Evaluation of a therapeutic carbon beam using a G2000 glass scintillator

A G2000 glass scintillator (G2000-SC) was used to determine the carbon profile and range of a 290-MeV/n carbon beam used in heavy-ion therapy because it was sensitive enough to detect single-ion hits at hundreds of mega electron Volts. An electron-multiplying charge-coupled device camera was used to detect the ion luminescence generated during the irradiation of G2000-SC with the beam. The resulting image showed that the position of the Bragg peak can be determined. The beam passes through the 112-mm-thick water phantom and stops 5.73 ± 0.03 mm from the incident side to the G2000-SC. Additionally, the location of the Bragg peak was simulated when irradiating G2000-SC with the beam using the Monte Carlo code particle and heavy ion transport system (PHITS). Simulation results show that the incident beam stops at 5.60 mm after entering G2000-SC. The beam stop location obtained from images and the PHITS code is defined at 80% distal fall-off from the Bragg peak position. Consequently, G2000-SC provided effective profile measurements of therapeutic carbon beams.

Development of a more accurate Geant4 quantum molecular dynamics model for hadron therapy

Objective. Although in heavy-ion therapy, the quantum molecular dynamics (QMD) model is one of the most fundamental physics models providing an accurate daughter-ion production yield in the final state, there are still non-negligible differences with the experimental results. The aim of this study is to improve fragment production in water phantoms by developing a more accurate QMD model in Geant4. Approach. A QMD model was developed by implementing modern Skyrme interaction parameter sets, as well as by incorporating with an ad hoc α-cluster model in the initial nuclear state. Two adjusting parameters were selected that can significantly affect the fragment productions in the QMD model: the radius to discriminate a cluster to which nucleons belong after the nucleus–nucleus reaction, denoted by R, and the squared standard deviation of the Gaussian packet, denoted by L. Squared Mahalanobis’s distance of fragment yields and angular distributions with 1, 2, and the higher atomic number for the produced fragments were employed as objective functions, and multi-objective optimization (MOO), which make it possible to compare quantitatively the simulated production yields with the reference experimental data, was performed. Main results. The MOO analysis showed that the QMD model with modern Skyrme parameters coupled with the proposed α-cluster model, denoted as SkM* α, can drastically improve light fragments yields in water. In addition, the proposed model reproduced the kinetic energy distribution of the fragments accurately. The optimized L in SkM* α was confirmed to be realistic by the charge radii analysis in the ground state formation. Significance. The proposed framework using MOO was demonstrated to be very useful in judging the superiority of the proposed nuclear model. The optimized QMD model is expected to improve the accuracy of heavy-ion therapy dosimetry.

Microdosimetry performance of the first multi-arrays of 3D-cylindrical microdetectors

The present work reports on the microdosimetry measurements performed with the two first multi-arrays of microdosimeters with the highest radiation sensitive surface covered so far. The sensors are based on new silicon-based radiation detectors with a novel 3D cylindrical architecture. Each system consists of arrays of independent microdetectors covering 2 mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×2 mm and 0.4 mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×12 cm radiation sensitive areas, the sensor distributions are arranged in layouts of 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×11 microdetectors and 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×3 multi-arrays, respectively. We have performed proton irradiations at several energies to compare the microdosimetry performance of the two systems, which have different spatial resolution and detection surface. The unitcell of both arrays is a 3D cylindrical diode with a 25 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm diameter and a 20 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm depth that results in a welldefined and isolated radiation sensitive micro-volume etched inside a silicon wafer. Measurements were carried out at the Accélérateur Linéaire et Tandem à Orsay (ALTO) facility by irradiating the two detection systems with monoenergetic proton beams from 6 to 20 MeV at clinical-equivalent fluence rates. The microdosimetry quantities were obtained with a spatial resolution of 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm and 600 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm for the 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×11 system and for the 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×3 multi-array system, respectively. Experimental results were compared with Monte Carlo simulations and an overall good agreement was found. The good performance of both microdetector arrays demonstrates that this architecture and both configurations can be used clinically as microdosimeters for measuring the lineal energy distributions and, thus, for RBE optimization of hadron therapy treatments. Likewise, the results have shown that the devices can be also employed as a multipurpose device for beam monitoring in particle accelerators.

Predicting the Biological Effects of Human Salivary Gland Tumour Cells for Scanned 4He-, 12C-, 16O-, and 20Ne-Ion Beams Using an SOI Microdosimeter

Experimental microdosimetry along with the microdosimetric kinetic (MK) model can be utilized to predict the biological effects of ions. To predict the relative biological effectiveness (RBE) of ions and the survival fraction (SF) of human salivary gland tumour (HSGc-C5) cells, microdosimetric quantities measured by a silicon-on-insulator (SOI) MicroPlus-mushroom microdosimeter along the spread-out Bragg peak (SOBP) delivered by pencil beam scanning of 4He, 12C, 16O, and 20Ne ions were used. The MK model parameters of HSGc-C5 cells were obtained from the best fit of the calculated SF for the different linear energy transfer (LET) of these ions and the formerly reported in vitro SF for the same LET and ions used for calculations. For a cube-shaped target of 10 × 10 × 6 cm3, treatment plans for 4He, 12C, 16O, and 20Ne ions were produced with proprietary treatment planning software (TPS) aiming for 10% SF of HSGc-C5 cells over the target volume and were delivered to a polymethyl methacrylate (PMMA) phantom. Afterwards, the saturation-corrected dose-mean lineal energy derived based on the measured microdosimetry spectra, along with the physical dose at various depths in PMMA phantoms, was used for the estimation of the SF, RBE, and RBE-weighted dose using the MK model. The predicted SF, RBE, and the RBE-weighted dose agreed with what was planned by the TPS within 3% at most depths for these ions.

Silicon 3D Microdosimeters for Advanced Quality Assurance in Particle Therapy

The Centre for Medical Radiation Physics introduced the concept of Silicon On Insulator (SOI) microdosimeters with 3-Dimensional (3D) cylindrical sensitive volumes (SVs) mimicking the dimensions of cells in an array. Several designs of high-definition 3D SVs fabricated using 3D MEMS technology were implemented. 3D SVs were fabricated in different sizes and configurations with diameters between 18 and 30 µm, thicknesses of 2–50 µm and at a pitch of 50 µm in matrices with volumes of 20 × 20 and 50 × 50. SVs were segmented into sub-arrays to reduce capacitance and avoid pile up in high-dose rate pencil beam scanning applications. Detailed TCAD simulations and charge collection studies in individual SVs have been performed. The microdosimetry probe (MicroPlus) is composed of the silicon microdosimeter and low-noise front–end readout electronics housed in a PMMA waterproof sheath that allows measurements of lineal energies as low as 0.4 keV/µm in water or PMMA. Microdosimetric quantities measured with SOI microdosimeters and the MicroPlus probe were used to evaluate the relative biological effectiveness (RBE) of heavy ions and protons delivered by pencil-beam scanning and passive scattering systems in different particle therapy centres. The 3D detectors and MicroPlus probe developed for microdosimetry have the potential to provide confidence in the delivery of RBE optimized particle therapy when introduced into routine clinical practice.

Response of SOI microdosimeter in fast neutron beams: experiment and Monte Carlo simulations

In this study, Monte Carlo codes, Geant4 and MCNP6, were used to characterize the fast neutron therapeutic beam produced at iThemba LABS in South Africa. Experimental and simulation results were compared using the latest generation of Silicon on Insulator (SOI) microdosimeters from the Centre for Medical Radiation Physics (CMRP). Geant4 and MCNP6 were able to successfully model the neutron gantry and simulate the expected neutron energy spectrum produced from the reaction by protons bombarding a ⁹Be target. The neutron beam was simulated in a water phantom and its characteristics recorded by the silicon microdosimeters; bare and covered by a ¹⁰B enriched boron carbide converter, at different positions. The microdosimetric quantities calculated using Geant4 and MCNP6 are in agreement with experimental measurements. The thermal neutron sensitivity and production of ¹⁰B capture products in the p+ boron-implanted dopant regions of the Bridge microdosimeter is investigated. The obtained results are useful for the future development of dedicated SOI microdosimeters for Boron Neutron Capture Therapy (BNCT). This paper provides a benchmark comparison of Geant4 and MCNP6 capabilities in the context of further applications of these codes for neutron microdosimetry.

Estimating the biological effects of helium, carbon, oxygen, and neon ion beams using 3D silicon microdosimeters

In this study, the survival fraction (SF) and relative biological effectiveness (RBE) of pancreatic cancer cells exposed to spread-out Bragg peak helium, carbon, oxygen, and neon ion beams are estimated from the measured microdosimetric spectra using a microdosimeter and the application of the microdosimetric kinetic (MK) model. To measure the microdosimetric spectra, a 3D mushroom silicon-on-insulator microdosimeter connected to low noise readout electronics (MicroPlus probe) was used. The parameters of the MK model were determined for pancreatic cancer cells such that the calculated SFs reproduced previously reported in vitro SF data. For a cuboid target of 10 × 10 × 6 cm³, treatment plans of helium, carbon, oxygen, and neon ion beams were designed using in-house treatment planning software (TPS) to achieve a 10% SF of pancreatic cancer cells throughout the target. The physical doses and microdosimetric spectra of the planned fields were measured at different depths in polymethyl methacrylate phantoms. The biological effects, such as SF, RBE, and RBE-weighted dose at different depths along the fields were predicted from the measurements. The predicted SFs at the target region were generally in good agreement with the planned SF from the TPS in most cases.

Reconstruction of dose distributions for fine carbon-ion beams using iterative approximation toward carbon-knife

For the practical application of carbon-knife with fine carbon-ion beams, the quantification of the dose distribution is essential and requires a high spatial resolution. We propose a novel method to quantify dose distributions with a spatial resolution smaller than the dosimeter size. The proposed method innovates the iterative reconstruction technique. Using a diode dosimeter with a sensitive area of 1 mm², two-dimensional dose-area-product (DAP) distributions were measured at a 0.1 mm step at the surface and near the Bragg peak depths for fine carbon-ion beams of ∼1 mm size at the full width at half maximum (FWHM). Then, the dose distributions were reconstructed with a spatial resolution of 0.1 × 0.1 mm² from the measured DAP distributions. However, an unnaturally high noise was observed in the reconstructed dose distributions, which were considered to originate from the measurement reproducibility errors of the DAP distributions estimated to be 2.5%-3%. Therefore, a low-pass filtering process was implemented to reduce the errors on the reconstructed dose distributions. The optimum cut-off frequencies of the low-pass filter were estimated depending on the amplitude of the induced noise. Using the filtering process with the obtained optimum cut-off frequency, the dose distribution was quantified with an average error of approximately 3% or less with respect to the peak value, when the actual measurement had an error of 3%. In the reconstructed dose rate distributions, a steep penumbra P_80-20 ∼ 0.2 mm was observed at the surface, and a dose rate at the center axis of ∼90 Gy s^-1 and a beam size of ∼1.1 mm at FWHM near the Bragg peak were obtained. The proposed method is expected to be useful for the measurement-based determination of microbeam models for commissioning and dose distribution calculations toward carbon-knife applications.

Neutron activation of gadolinium for ion therapy: a Monte Carlo study of charged particle beams

This study investigates the photon production from thermal neutron capture in a gadolinium (Gd) infused tumor as a result of secondary neutrons from particle therapy. Gadolinium contrast agents used in MRI are distributed within the tumor volume and can act as neutron capture agents. As a result of particle therapy, secondary neutrons are produced and absorbed by Gd in the tumor providing potential enhanced localized dose in addition to a signature photon spectrum that can be used to produce an image of the Gd enriched tumor. To investigate this imaging application, Monte Carlo (MC) simulations were performed for 10 different particles using a 5-10 cm spread out-Bragg peak (SOBP) centered on an 8 cm3, 3 mg/g Gd infused tumor. For a proton beam, 1.9 × 106 neutron captures per RBE weighted Gray Equivalent dose (GyE) occurred within the Gd tumor region. Antiprotons ([Formula: see text]), negative pions (- π), and helium (He) ion beams resulted in 10, 17 and 1.3 times larger Gd neutron captures per GyE than protons, respectively. Therefore, the characteristic photon based spectroscopic imaging and secondary Gd dose enhancement could be viable and likely beneficial for these three particles.

Radioresistant tumours: From identification to targeting

From surviving fraction to tumour curability, definitions of tumour radioresistance may vary depending on the view angle. Yet, mechanisms of radioresistance have been identified and involve tumour-specific oncogenic signalling pathways, tumour metabolism and proliferation, tumour microenvironment/hypoxia, genomics. Correlations between tumour biology (histology) and imaging allow theragnostic approaches that use non-invasive biological imaging using tracer functionalization of tumour pathway biomarkers, imaging of hypoxia, etc. Modelling dose prescription function based on their tumour radio-resistant factor enhancement ratio, related to metabolism, proliferation, hypoxia is an area of investigation. Yet, the delivery of dose painting by numbers/voxel-based radiotherapy with low lineal energy transfer particles may be limited by the degree of modulation complexity needed to achieve the doses needed to counteract radioresistance. Higher lineal energy transfer particles or combinations of different particles, or combinations with drugs and devices such as done with radioenhancing nanoparticles may be promising.

The impact of sensitive volume thickness for silicon on insulator microdosimeters in hadron therapy

Compact silicon on insulator (SOI) microdosimeters have been used to characterise the radiation field of many different hadron therapy beams. SOI devices are particularly attractive in hadron therapy fields due to their spatial resolution being well suited to the sharp dose gradients at the end of the primary beam's range. Due to the small size of SOI's sensitive volumes (SVs), which are usually  ∼1-10 [Formula: see text]m thick, the fabrication of these devices can present challenges which are not as common for more conventional thickness silicon devices such as silicon spectroscopy detectors. Microdosimetry is the study of the energy deposition in micrometre sized volumes representing biological sites and is a powerful approach to estimate the biological effect of radiation on the micron-scale level, in a cell. However, cell sizes vary extensively translating in different energy deposition spectra. This work studies SV thicknesses between 1 and 100 [Formula: see text]m using Geant4 and examines the impact of SV dimensions on microdosimetric quantities. The quantities studied were the frequency mean lineal energy, [Formula: see text], and the dose mean lineal energy, [Formula: see text]. Additionally the relative biological effectiveness (RBE), estimated by the microdosimetric kinetic model (MKM), is also investigated. To study the impact of the SV thickness, SOI microdosimeters were irradiated with proton, [Formula: see text] and [Formula: see text] ion beams with ranges of  ∼160 mm, with the microdosimeter being set at various positions along the Bragg curve. It was found that [Formula: see text] was influenced the least in proton beams and increased for heavier ion beams. Conversely, [Formula: see text] was impacted by the SV thickness the most in proton beams and [Formula: see text] was the least. Similar to [Formula: see text], protons were impacted the most by the SV thickness when estimating the RBE using the MKM. The cause of these differences was largely due to the different densities of the delta electron track structure for the case of [Formula: see text] and the energy transferred to the medium from the primary beam for [Formula: see text].

SOI microdosimetry and modified MKM for evaluation of relative biological effectiveness for a passive proton therapy radiation field

With more patients receiving external beam radiation therapy with protons, it becomes increasingly important to refine the clinical understanding of the relative biological effectiveness (RBE) for dose delivered during treatment. Treatment planning systems used in clinics typically implement a constant RBE of 1.1 for proton fields irrespective of their highly heterogeneous linear energy transfer (LET). Quality assurance tools that can measure beam characteristics and quantify or be indicative of biological outcomes become necessary in the transition towards more sophisticated RBE weighted treatment planning and for verification of the Monte Carlo and analytical based models they use. In this study the RBE for the CHO-K1 cell line in a passively delivered clinical proton spread out Bragg peak (SOBP) is determined both in vitro and using a silicon-on-insulator (SOI) microdosimetry method paired with the modified microdosimetric kinetic model. The RBE along the central axis of a SOBP with 2 Gy delivered at the middle of the treatment field was found to vary between 1.11-1.98 and the RBE for 10% cell survival between 1.07-1.58 with a 250 kVp x-ray reference radiation and between 1.19-2.34 and 0.95-1.41, respectively, for a Co60 reference. Good agreement was found between RBE values calculated from the SOI-microdosimetry-MKM approach and in vitro. A strong correlation between proton lineal energy and RBE was observed particularly in the distal end and falloff of the SOBP.

Optimisation of the design of SOI microdosimeters for hadron therapy quality assurance

Silicon-on-insulator (SOI) microdosimeters offer a promising method for routine quality assurance (QA) for hadron therapy due to their ease of operation and high spatial resolution. However, one complication which has been shown previously is that the traditional use of the mean chord length, [Formula: see text], calculated using Cauchy's formula, for SOI devices in clinical carbon ion fields is not appropriate due to the strong directionality of the radiation field. In a previous study, we demonstrated that the mean path length, [Formula: see text], which is the mean path of charged particles in the sensitive volume (SV), is a more appropriate method to obtain microdosimetric quantities and biological relevant values, namely the relative biological effectiveness (RBE) by means of the microdosimetric kinetic model. The previous work, which was limited to mono-energetic [Formula: see text] ion beams typical of heavy ion therapy (HIT), is extended here to investigate the [Formula: see text] in a pristine proton beam as well as for spread out Bragg peaks (SOBP) for both proton and carbon ion clinical beams. In addition, the angular dependence of the SOI device for a number of different SV designs is also investigated to quantify the effects which the alignment has on the [Formula: see text]. It is demonstrated that the [Formula: see text] can be accurately estimated along the depth of a pristine or SOBP using the energy deposition spectra for both proton and [Formula: see text] ion beams. This observation allows a quick and accurate estimation of the [Formula: see text] for experimental use.